Molecular Dynamics Simulation for All.
Ontology highlight
ABSTRACT:
SUBMITTER: Hollingsworth SA
PROVIDER: S-EPMC6209097 | biostudies-literature | 2018 Sep
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus-host interactions and evolution, and present future perspectives on this technique.
| S-EPMC3400276 | biostudies-literature
Project description:Allosteric networks allow enzymes to transmit information and regulate their catalytic activities over vast distances. In principle, molecular dynamics (MD) simulations can be used to reveal the mechanisms that underlie this phenomenon; in practice, it can be difficult to discern allosteric signals from MD trajectories. Here, we describe how MD simulations can be analyzed to reveal correlated motions and allosteric networks, and provide an example of their use on the coagulation enzyme thrombin. Methods are discussed for calculating residue-pair correlations from atomic fluctuations and mutual information, which can be combined with contact information to identify allosteric networks and to dynamically cluster a system into highly correlated communities. In the case of thrombin, these methods show that binding of the antagonist hirugen significantly alters the enzyme's correlation landscape through a series of pathways between Exosite I and the catalytic core. Results suggest that hirugen binding curtails dynamic diversity and enforces stricter venues of influence, thus reducing the accessibility of thrombin to other molecules.
| S-EPMC5492730 | biostudies-literature
Project description:In this study, we construct novel RNA nanoclusters, RNA nanotubes made of several nanorings up to the size of 20 nm, utilizing the molecular dynamics simulation, and study their structural properties [i.e., the root mean square deviation, the radius of gyration and the radial distribution function (RDF)] in physiological solutions that can be used for drug delivery into the human body. The patterns of energy and temperature variations of the systems are also discussed. Furthermore, we study the concentration of ions around the tube as a function of time at a particular temperature. We have found that when the temperature increases, the number of ions increases within a certain distance of the tube. We report that the number of ions within this distance around the tubes decreases in quenched runs. This indicates that some ions evaporate with decrease in temperature, as has been observed in the case of the nanoring. RDF plots also demonstrate a similar trend with temperature, as was found in the case of RNA nanorings.
| S-EPMC6345520 | biostudies-literature
Project description:Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.
| S-EPMC1304655 | biostudies-literature
Project description:The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.
| S-EPMC8340580 | 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 nonpolarizable force field for alkyl nitrates developed by Borodin et al. [J. Phys. Chem. B, 2008, 112, 734-742] has been employed to calculate selected properties of crystalline and liquid erythritol tetranitrate (ETN). The set of partial charges proposed by Borodin for pentaerythritol tetranitrate (PETN) was used except for a small correction to the H atom charges to ensure charge neutrality owing to the absence of the neopentyl carbon in ETN. The force field was used to compute the isothermal compression curve, lattice parameters, heat capacity, thermal expansivity, single crystal elastic constants, and Gruneisen parameters of crystalline ETN. The density- and temperature-dependent viscosities of liquid ETN are also reported. We anticipate that these data will be of some utility to the development of equations of state and thermomechanical models for ETN.
| S-EPMC11000242 | biostudies-literature
Project description:Neuroglobin (Ngb) is a monomeric protein that, despite the small sequence similarity with other globins, displays the typical globin fold. In the absence of exogenous ligands, the ferric and the ferrous forms of Ngb are both hexacoordinated to the distal and proximal histidines. In the ferrous form, oxygen, nitric oxide or carbon monoxide can displace the distal histidine, yielding a reversible adduct. Crystallographic data show that the binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities. Molecular dynamics (MD) simulations in solution show that the heme oscillates between two positions, much as the ones observed in the crystal structure, although the occupancy is different. The simulations also suggest that ligand binding in solution can affect the flexibility and conformation of residues connecting the C and D helices, referred to as the CD corner, which is coupled to the configuration adopted by the distal histidine. In this study, we report the results of 30 ns MD simulations of CO-bound Ngb in the crystal. Our goal was to compare the protein dynamical behavior in the crystal with the results supplied by the previous MD simulation of CO-bound Ngb in solution and the x-ray experimental data. The results show that the different environments (crystal or solution) affect the dynamics of the heme group and of the CD corner.
| S-EPMC2567928 | biostudies-literature
Project description:The generalized Newton-Euler inverse mass operator (GNEIMO) method is an advanced method for internal coordinates molecular dynamics (ICMD). GNEIMO includes several theoretical and algorithmic advancements that address longstanding challenges with ICMD simulations. In this article, we describe the GneimoSim ICMD software package that implements the GNEIMO method. We believe that GneimoSim is the first software package to include advanced features such as the equipartition principle derived for internal coordinates, and a method for including the Fixman potential to eliminate systematic statistical biases introduced by the use of hard constraints. Moreover, by design, GneimoSim is extensible and can be easily interfaced with third party force field packages for ICMD simulations. Currently, GneimoSim includes interfaces to LAMMPS, OpenMM, and Rosetta force field calculation packages. The availability of a comprehensive Python interface to the underlying C++ classes and their methods provides a powerful and versatile mechanism for users to develop simulation scripts to configure the simulation and control the simulation flow. GneimoSim has been used extensively for studying the dynamics of protein structures, refinement of protein homology models, and for simulating large scale protein conformational changes with enhanced sampling methods. GneimoSim is not limited to proteins and can also be used for the simulation of polymeric materials.
| S-EPMC4211970 | biostudies-literature
Project description:BACKGROUND:Kinase-inducible domain (KID) as transcriptional activator can stimulate target gene expression in signal transduction by associating with KID interacting domain (KIX). NMR spectra suggest that apo-KID is an unstructured protein. After post-translational modification by phosphorylation, KID undergoes a transition from disordered to well folded protein upon binding to KIX. However, the mechanism of folding coupled to binding is poorly understood. METHODOLOGY:To get an insight into the mechanism, we have performed ten trajectories of explicit-solvent molecular dynamics (MD) for both bound and apo phosphorylated KID (pKID). Ten MD simulations are sufficient to capture the average properties in the protein folding and unfolding. CONCLUSIONS:Room-temperature MD simulations suggest that pKID becomes more rigid and stable upon the KIX-binding. Kinetic analysis of high-temperature MD simulations shows that bound pKID and apo-pKID unfold via a three-state and a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound pKID folds in the order of KIX access, initiation of pKID tertiary folding, folding of helix alpha(B), folding of helix alpha(A), completion of pKID tertiary folding, and finalization of pKID-KIX binding. Our data show that the folding pathways of apo-pKID are different from the bound state: the foldings of helices alpha(A) and alpha(B) are swapped. Here we also show that Asn139, Asp140 and Leu141 with large Phi-values are key residues in the folding of bound pKID. Our results are in good agreement with NMR experimental observations and provide significant insight into the general mechanisms of binding induced protein folding and other conformational adjustment in post-translational modification.
| S-EPMC2717803 | biostudies-literature