Applications of Molecular Dynamics Simulation in Protein Study.
Ontology highlight
ABSTRACT:
SUBMITTER: Sinha S
PROVIDER: S-EPMC9505860 | biostudies-literature | 2022 Aug
REPOSITORIES: biostudies-literature
Publications
Applications of Molecular Dynamics Simulation in Protein Study.
Sinha Siddharth S Tam Benjamin B Wang San Ming SM
Membranes 20220829 9
Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status ...[more]
Similar Datasets
Project description:Cholesteryl ester transfer protein (CETP) mediates the net transfer of cholesteryl esters (CEs) from atheroprotective high-density lipoproteins (HDLs) to atherogenic low-density lipoproteins (LDLs) or very-low-density lipoproteins (VLDLs). Inhibition of CETP raises HDL cholesterol (good cholesterol) levels and reduces LDL cholesterol (bad cholesterol) levels, making it a promising drug target for the prevention and treatment of coronary heart disease. Although the crystal structure of CETP has been determined, the molecular mechanism mediating CEs transfer is still unknown, even the structural features of CETP in a physiological environment remain elusive. We performed molecular dynamics simulations to explore the structural features of CETP in an aqueous solution. Results show that the distal portion flexibility of N-terminal β-barrel domain is considerably greater in solution than in crystal; conversely, the flexibility of helix X is slightly less. During the simulations the distal end of C-terminal β-barrel domain expanded while the hydrophilic surface increasing more than the hydrophobic surface. In addition, a new surface pore was generated in this domain. This surface pore and all cavities in CETP are stable. These results suggest that the formation of a continuous tunnel within CETP by connecting cavities is permitted in solution.
| S-EPMC3557553 | biostudies-literature
Project description:The structure and stability of membrane proteins can vary widely in different detergents and this variability has great practical consequences for working with membrane proteins. Nevertheless, the mechanisms that operate to alter the behavior of proteins in micelles are poorly understood and not predictable. Atomic simulations could provide considerable insight into these mechanisms. Building protein-micelle complexes for simulation is fraught with uncertainty, however, in part because it is often unknown how many detergent molecules are present in the complex. Here, we describe several convenient ways to employ Micelle Builder in CHARMM-GUI to rapidly construct protein-micelle complexes and performed simulations of the isolated voltage-sensor domain of voltage-dependent potassium-selective channel and an antimicrobial peptide papiliocin with varying numbers of detergents. We found that once the detergent number exceeds a threshold, protein-detergent interactions change very little and remain very consistent with experimental observations. Our results provide a platform for future studies of the interplays between protein structure and detergent properties at the atomic level. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
| S-EPMC4877248 | biostudies-literature
Project description:The estimation of protein model quality remains a challenging task and is important for protein structural model utilization. In the last decade, existing methods that rely on machine learning to deep learning have been developed and shown progressive improvement. Despite utilizing more sophisticated techniques and introducing new features, none of these methods employ explicit protein structure stability information. Hypothetically, protein model quality might be indicated by its structural stability in an in silico system disclosed by the structural difference from its initial structure. One of the possible methods to exploit such information is by implementing molecular dynamics simulations that have shown successful applications in many research fields. We present a novel approach by introducing explicit protein structure stability information using molecular dynamics simulation. Despite using only simple features, small data with no training process required, and a short molecular dynamics simulation time, our method shows comparable performance to the state-of-the-art deep learning-based method.
| S-EPMC9301944 | biostudies-literature
Project description:N-linked glycosylation is one of the most important, chemically complex, and ubiquitous post-translational modifications in all eukaryotes. The N-glycans that are covalently linked to proteins are involved in numerous biological processes. There is considerable interest in developments of general approaches to predict the structural consequences of site-specific glycosylation and to understand how these effects can be exploited in protein design with advantageous properties. In this study, the impacts of N-glycans on protein structure and dynamics are systematically investigated using an integrated computational approach of the Protein Data Bank structure analysis and atomistic molecular dynamics simulations of glycosylated and deglycosylated proteins. Our study reveals that N-glycosylation does not induce significant changes in protein structure, but decreases protein dynamics, likely leading to an increase in protein stability. Overall, these results suggest not only a common role of glycosylation in proteins, but also a need for certain proteins to be properly glycosylated to gain their intrinsic dynamic properties.
| S-EPMC4352867 | biostudies-literature
Project description:BackgroundTrp cage is a recently-constructed fast-folding miniprotein. It consists of a short helix, a 3,10 helix and a C-terminal poly-proline that packs against a Trp in the alpha helix. It is known to fold within 4 ns.ResultsHigh-temperature unfolding molecular dynamics simulations of the Trp cage miniprotein have been carried out in explicit water using the OPLS-AA force-field incorporated in the program GROMACS. The radius of gyration (Rg) and Root Mean Square Deviation (RMSD) have been used as order parameters to follow the unfolding process. Distributions of Rg were used to identify ensembles.ConclusionThree ensembles could be identified. While the native-state ensemble shows an Rg distribution that is slightly skewed, the second ensemble, which is presumably the Transition State Ensemble (TSE), shows an excellent fit. The denatured ensemble shows large fluctuations, but a Gaussian curve could be fitted. This means that the unfolding process is two-state. Representative structures from each of these ensembles are presented here.
| S-EPMC555740 | biostudies-literature
Project description:Technical strategies like amino acid substitution and residue modification have been widely used to characterize the importance of key amino acids and the role that each residue plays in the structural and functional properties of protein molecules. However, there is no systematic approach to assess the impact of the substituted/modified amino acids on the conformational dynamics of proteins. In this investigation to clarify the effects of residue modifications on the structural dynamics of human prion protein (PrP), a comparative molecular dynamics simulation study on the native and the amino acid-substituted analog at position 208 of PrP has been performed. It is believed that Arginine to Histidine mutation at position 208 is responsible for the structural transition of the native form of human prion protein to the pathogenic isoform causing Creutzfeldt-Jakob disease (CJD). So, three 10 ns molecular dynamics simulations on three model constructs have been performed. Simulation results indicated considerable differences of conformational fluctuations for Alanine substituted construct (PrPALA) and the analog form (PrPSB) comprising the neutralized state of the Arginine residue at position 208 of the human prion protein. According to our data, substitution of the Arginine residue by the uncharged state of this residue induces some reversible structural alterations in the intrinsically flexible loop area including residues 167-171 of PrP. Thus, deprotonation of Arg(208) is a weak perturbation to the structural fluctuations of the protein backbone and the resulting construct behaves almost identical as its native form. Otherwise, Alanine substitution at position 208 imposed an irreversible impact on the secondary and tertiary structure of the protein, which leads to conformational instabilities in the remote hot region comprising residues 190-195 of the C-terminal part of helix 2. Based on the results, it could be deduced that the observed conformational transitions upon Arg(208) to His point mutation, which is the main reason for CJD, may be mainly related to the structural instabilities due to the induced-conformational changes that caused alterations in local/spatial arrangements of the force distributions in the backbone of the human prion protein.
| S-EPMC4464497 | biostudies-literature
Project description:Auxin Binding Protein 1 (ABP1) is ubiquitous in green plants. It binds the phytohormone auxin with high specificity and affinity, but its role in auxin-induced processes is unknown. To understand the proposed receptor function of ABP1 we carried out a detailed molecular modeling study. Molecular dynamics simulations showed that ABP1 can adopt two conformations differing primarily in the position of the C-terminus and that one of them is stabilized by auxin binding. This is in agreement with experimental evidence that auxin induces changes at the ABP1 C-terminus. Simulations of ligand egress from ABP1 revealed three main routes by which an auxin molecule can enter or leave the ABP1 binding site. Assuming the previously proposed orientation of ABP1 to plant cell membranes, one of the routes leads to the membrane and the other two to ABP1's aqueous surroundings. A network of hydrogen-bonded water molecules leading from the bulk water to the zinc-coordinated ligands in the ABP1 binding site was formed in all simulations. Water entrance into the zinc coordination sphere occurred simultaneously with auxin egress. These results suggest that the hydrogen-bonded water molecules may assist in protonation and deprotonation of auxin molecules and their egress from the ABP1 binding site.
| S-EPMC2134879 | biostudies-literature
Project description:Biomolecules localize and function in microenvironments where their local concentration, spatial organization, and biochemical reactivity are regulated. To compartmentalize and control the local properties of the native microenvironment, cellular mimics and artificial bioreactors have been developed to study the properties of membraneless organelles or mimic the bio-environment for life origin. Here, we carried out molecular dynamics simulation with the Martini 3.0 model to reproduce the experimental salt concentration and pH dependency of different complex coacervates. We showed that coacervates inside vesicles are able to change their shape. In addition, we used these coacervate systems to explore the partitioning of the ubiquitous cytoskeletal protein actin and found that actin spontaneously partitions to all the coacervate peripheries. Therefore, we believe that our study can provide a better understanding of the versatile coacervate platform, where biomolecules partition and gather to fulfill their biological duties.
| S-EPMC9891350 | biostudies-literature
Project description:Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.
| S-EPMC1302733 | biostudies-literature
Project description:Molecular dynamics (MD) simulations of proteins provide important information to understand their functional mechanisms, which are, however, likely to be hidden behind their complicated motions with a wide range of spatial and temporal scales. A straightforward and intuitive analysis of protein dynamics observed in MD simulation trajectories is therefore of growing significance with the large increase in both the simulation time and system size. In this study, we propose a novel description of protein motions based on the hierarchical clustering of fluctuations in the inter-atomic distances calculated from an MD trajectory, which constructs a single tree diagram, named a "Motion Tree", to determine a set of rigid-domain pairs hierarchically along with associated inter-domain fluctuations. The method was first applied to the MD trajectory of substrate-free adenylate kinase to clarify the usefulness of the Motion Tree, which illustrated a clear-cut dynamics picture of the inter-domain motions involving the ATP/AMP lid and the core domain together with the associated amplitudes and correlations. The comparison of two Motion Trees calculated from MD simulations of ligand-free and -bound glutamine binding proteins clarified changes in inherent dynamics upon ligand binding appeared in both large domains and a small loop that stabilized ligand molecule. Another application to a huge protein, a multidrug ATP binding cassette (ABC) transporter, captured significant increases of fluctuations upon binding a drug molecule observed in both large scale inter-subunit motions and a motion localized at a transmembrane helix, which may be a trigger to the subsequent structural change from inward-open to outward-open states to transport the drug molecule. These applications demonstrated the capabilities of Motion Trees to provide an at-a-glance view of various sizes of functional motions inherent in the complicated MD trajectory.
| S-EPMC4492737 | biostudies-literature