Molecular dynamics simulation techniques for determination of molecular structures from nuclear magnetic resonance data.
Explore the source record for details and available documents.
Magnetic Resonance Spectroscopy↗
SEARCH · PubMed Health
Results for “Molecular Dynamic simulation”
Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Molecular dynamics simulations of "loop closing" in the enzyme triose phosphate isomerase.
We present molecular dynamics simulations on the active site region of dimeric triose phosphate isomerase (TIM) using the co-ordinates of native chicken muscle TIM as a starting point and performing simulations with no substrate, with dihydroxyacetone phosphate (DHAP), the natural substrate, and with dihydroxyacetone sulfate (DHAS), a substrate analog. Whereas most of the protein moves less than 1 A during the simulation, some residues in the active site loop move more than 8 A during the 10.5 picoseconds of dynamics for each of the simulations. Most interestingly, the nature of the loop motion depends on the substrate, with the largest motion found in the presence of DHAP, and only in the presence of DHAP does the loop move to "close off" the active site pocket. The final structure found for the DHAP-chicken TIM complex is qualitatively similar to that described by Alber et al. for DHAP-yeast TIM. Simulations on the monomeric protein gives insight into why the molecule is active only as a dimer.
Acetone↗
A sampling problem in molecular dynamics simulations of macromolecules.
Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.
Computer Simulation↗
Molecular dynamics simulations of cyclosporin A: the crystal structure and dynamic modelling of a structure in apolar solution based on NMR data.
The conformation of the immunosuppressive drug cyclosporin A (CPA), both in apolar solution and in crystalline state, has been studied by computer simulation techniques. Three molecular dynamics (MD) simulations have been performed: one modelling the crystal structure and two modelling the structure in apolar solution, using a restrained MD approach in which data from nuclear magnetic resonance (NMR) and infrared (IR) spectroscopy are taken into account. The simulation of the crystalline state (MDC) concerns a system of 4 unit cells containing 16 cyclosporin A molecules and 22 water molecules, which is simulated using crystalline periodic boundary conditions. The simulations modelling the apolar solvent conformation (MDS) concern one isolated cyclosporin A molecule. In these simulations an extra term in the interatomic potential function is used, which forces the molecule to satisfy a set of 57 atom-atom distance constraints originating from nuclear Overhauser effects (NOEs) obtained from NMR spectroscopy and one distance constraint deduced from IR spectroscopy. From a comparison of the results of the crystal simulation to those of the X-ray experiment in terms of structure, atomic fluctuations, hydrogen bond pattern, etc., it is concluded that the force field that is used yields an adequate representation of crystalline cyclosporin A. Secondly, it is shown that the dynamic modelling technique that is used to obtain a structure in a polar solution from NMR distance information works well. Starting from initial conformations which have a root mean square difference of 0.14 nm both distance restrained MD simulations converge to the same final solution structure. A comparison of the crystal structure of cyclosporin A and the one in apolar solution shows that there are significant differences. The overall difference in atomic positions is 0.09 nm for the C alpha atoms and 0.17 nm for all atoms. In apolar solution, the molecule is slightly more bent and the side chains of 1 MeBmt and 10 MeLeu adopt a different conformation.
Computer Simulation↗
Hydration of an alpha-helical peptide: comparison of theory and molecular dynamics simulation.
We present a statistical mechanical description of biomolecular hydration that accurately describes the hydrophobic and hydrophilic hydration of a model alpha-helical peptide. The local density of water molecules around a biomolecule is obtained by means of a potential-of-mean-force (PMF) expansion in terms of pair- and triplet-correlation functions of bulk water and dilute solutions of nonpolar atoms. The accuracy of the method is verified by comparing PMF results with the local density and site-site correlation functions obtained by molecular dynamics simulations of a model alpha-helix in solution. The PMF approach quantitatively reproduces all features of the peptide hydration determined from the molecular dynamics simulation. Regions of hydrophobic hydration near the C alpha and C beta atoms along the helix are well reproduced. The hydration of exposed polar groups at the N- and C-termini of the helix are also well described by the theory. A detailed comparison of the local hydration by means of site-site radial distribution functions evaluated with the PMF theory shows agreement with the molecular dynamics simulations. The formulation of this theory is general and can be applied to any biomolecular system. The accuracy, speed of computation, and local character of this theory make it especially suitable for studying large biomolecular systems.
Computer Simulation↗
Molecular dynamics simulations of oligonucleotides in solution: visualization of intrinsic curvature.
We have undertaken molecular dynamics simulations on the d(CGCAAAAAAGCG).d(CGCTTTTTTGCG) dodecamer in solution. In this study, we focus on aspects of conformation and dynamics, including the possibility of cross-strand hydrogen bonds. We compare our results with those from crystallography as well as infrared, Raman and NMR spectroscopy and cyclization kinetics. Our method of analysis allows us to visualise the curvature of the helix as a function of time during the simulation. We find that the major distortions of the helix axis path occur at the junctions between the (essentially straight) A-tract and the CG- and GC-tracts, although at one junction this is due to hyperflexibility (i.e., regions of high flexibility with no preferred direction of curvature), while at the other junction a static curvature is found (i.e., a preferred, sustained direction of curvature).
Base Sequence↗
Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.
Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.
Amino Acid Sequence↗
Mass-weighted molecular dynamics simulation of cyclic polypeptides.
A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.
Amino Acid Sequence↗
Structure and dynamics of the active site gorge of acetylcholinesterase: synergistic use of molecular dynamics simulation and X-ray crystallography.
The active site of acetylcholinesterase (AChE) from Torpedo californica is located 20 A from the enzyme surface at the bottom of a narrow gorge. To understand the role of this gorge in the function of AChE, we have studied simulations of its molecular dynamics. When simulations were conducted with pure water filling the gorge, residues in the vicinity of the active site deviated quickly and markedly from the crystal structure. Further study of the original crystallographic data suggests that a bis-quaternary decamethonium (DECA) ion, acquired during enzyme purification, residues in the gorge. There is additional electron density within the gorge that may represent small bound cations. When DECA and 2 cations are placed within the gorge, the simulation and the crystal structure are dramatically reconciled. The small cations, more so than DECA, appear to stabilize part of the gorge wall through electrostatic interactions. This part of the gorge wall is relatively thin and may regulate substrate, product, and water movement through the active site.
Acetylcholinesterase↗
Threonine6-bradykinin: molecular dynamics simulations in a biphasic membrane mimetic.
The natural peptide [Thr6]-bradykinin, Arg1-Pro2-Pro3-Gly4-Phe5-Thr6-Pro7-Phe8-Arg9, has been conformationally examined by molecular dynamics simulations using a two-phase box consisting of H2O and CCl4 to mimic the micellar environment utilized in the 1H-NMR studies. The different conformations generated from distance geometry calculations were refined with extensive molecular dynamics simulations. The resulting conformations provide additional structural insight into the differing biological activities of native bradykinin and [Thr6]-bradykinin, produced by the one conservative substitution Thr6 for Ser6. In addition, the simulations give some indication of the interaction of the peptide with the biphasic, hydrophilic/hydrophobic environment of the micelle. Such information is vital given the accumulating data indicating that the peptide first interacts with the membrane before the membrane-bound receptor. The structures of membrane-bound [Thr6]-bradykinin developed here provide experimental support for the interaction of residues 7 and 8 with the core of the membrane-bound receptor and the N-terminus and C-terminal arginine interacting with the extracellular portion of the receptor.
Bradykinin↗
Molecular dynamics simulations of a lipid bilayer and of hexadecane: an investigation of membrane fluidity.
Molecular dynamics simulations of a fluid-phase dipalmitoyl phosphatidylcholine lipid bilayer in water and of neat hexadecane are reported and compared with nuclear magnetic resonance spin-lattice relaxation and quasi-elastic neutron scattering data. On the 100-picosecond time scale of the present simulations, there is effectively no difference in the reorientational dynamics of the carbons in the membrane interior and in pure hexadecane. Given that the calculated fast reorientational correlation times and the "microscopic" lateral diffusion of the lipids show excellent agreement with the experimental results, it is concluded that the apparently high viscosity of the membrane is more closely related to molecular interactions on the surface rather than in the interior.
1,2-Dipalmitoylphosphatidylcholine↗
Molecular dynamics simulations of phospholipases A2.
An extensive molecular dynamics study of phospholipases A2 from pancreatic bovine and Crotalus atrox venom has shown that the well-conserved homologous core of the phospholipases A2, including the so called catalytic network, is very stable during the course of the calculations. The fluctuations which occur are located in segments which have significantly different three-dimensional conformations in the two phospholipases A2 studied, suggesting that a particularly stable core conformation gives rise to a large homologous family of similar three-dimensional structure. The calcium ion, which exhibits a crucial structural role in the monomeric phospholipases A2, appears not to be required to stabilize the C.atrox dimer. Moreover, the behaviour of the dimeric structure during the dynamics raises the question of a possible dissociation of the two subunits into functional monomers.
Animals↗
A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation.
Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.
Molecular Dynamics Simulation↗
Molecular dynamics simulations of the unfolding of barnase in water and 8 M aqueous urea.
Molecular dynamics simulations of barnase have been conducted both in water and in 8 M urea solution for 500 ps at 25 degrees C and for 2000 ps at 85 degrees C. The final structure of the aqueous simulation at room temperature matches closely the structure obtained by NMR and the experimentally observed protections from isotopic exchange. The comparison of the structures generated by the aqueous simulation at 85 degrees C reveals a trajectory composed of groups of geometrically related structures separated by narrow regions of rapid change in structure. The first of these regions displays changes in backbone rmsd to the crystal structure and solvent-accessible area suggestive of a transition state, while the properties observed during the final 300 ps of the simulation are consistent with a stable intermediate. These assignments were confirmed by calculation of the "progress along the reaction coordinate" phi-values using an empirical equation based on a linear response method. The pathway of unfolding defined in this fashion agrees well with the experimental results of site-directed mutagenesis in terms of secondary structure content of the transition state and the intermediate and reproduces the relative stability of the different elements of secondary structure. The results of the simulations in urea suggest a mechanism at the molecular level for its well-known enhancement of the denaturation of proteins. The analysis of radial distribution functions shows that the first solvation shell of the protein is enriched in urea relative to the bulk solvent. The displacement of water molecules allows greater exposure of hydrophobic side chains, as witnessed particularly in the analysis of solvent-accessible surface areas at the higher temperature. Almost all urea molecules in the first shell form at least one hydrogen bond with the protein. They provide a more favorable environment for accommodation of the remaining water molecules, and they facilitate the separation of secondary structure elements by acting as a bridge between groups previously forming intraprotein hydrogen bonds.
Bacterial Proteins↗
Molecular dynamics simulations of the N-linked oligosaccharide of the lectin from Erythrina corallodendron.
Molecular dynamics simulations have been used to model the flexibility of the seven-sugar oligosaccharide of the lectin from Erythrina corallodendron in three separate simulations: one of the isolated oligosaccharide in vacuo, one of the oligosaccharide in solution and one of the oligosaccharide linked to the protein in solution. Adiabatic conformational energy maps were prepared for each of the disaccharide linkages as a means of interpreting the observed dynamics and conformational averages in terms of intramolecular energy. The inclusion of aqueous solvent molecules appears to be necessary to reproduce the experimental conformational behavior, which also cannot be predicted well from conformational energy maps for the disaccharide linkages alone. The crystallographically determined conformation does not appear to be induced by the crystal dimerization, but is rather stable in solution. The build-up of fluctuations along the successive linkages of the oligosaccharide is significant and would be sufficient to prevent branch residues from being located in most crystal structure determinations. Good general agreement between the calculated solution structure and the average structure determined by NMR was found for most of the oligosaccharide linkages.
Carbohydrate Conformation↗
FOCUS: a program for analyzing molecular dynamics simulations, featuring digital signal-processing techniques.
FOCUS is a program for analyzing molecular dynamics simulations. It enables the researcher to monitor structural and energetic properties during the trajectory, and to calculate the corresponding statistical averages, correlation functions and Fourier transforms. In addition to these conventional methods, the program also utilizes novel methods based on digital signal-processing techniques to characterize the various motions. The characteristic frequencies in the system are revealed by the frequency distribution function g(v), which is calculated from the Fourier transform of the atomic coordinates. A filtering technique is employed to remove uninteresting motion (e.g., high-frequency bond stretching) while retaining and focusing on important motion (e.g., low-frequency conformational motion). The filtering technique enables fast display of slow events without getting a blurry or jittery picture due to the high-frequency motions. Another new way for analyzing the motion is by extracting "characteristic modes" and associated frequencies. This yields a pictorial description of the oscillatory motions in a manner analogous to normal mode analysis.
Acetamides↗
Molecular dynamics simulations of trp apo- and holorepressors: domain structure and ligand-protein interaction.
Molecular dynamics simulations of the apo- and holo-forms of the trp-repressor protein were performed under extensively solvated conditions in order to elucidate their dynamic structures and ligand-protein interactions. The root mean square fluctuations calculated from the trajectories agreed with those calculated from X-ray temperature factors. Distance, distance fluctuation, and dynamic cross-correlation maps were drawn to provide information on the dynamic structures and communications among the domains. A three-domain format has been proposed for the crystal structure (Zhang et al., Nature 327:591-597, 1987); namely, helices A-C and F of both subunits make up a central core, and D and E of each subunit forms a DNA binding head. The results of the simulations were mostly consistent with the three-domain format. However, helix F was more flexible and freer than other parts of the central core. The turn DE, the helix-turn-helix DNA binding motif, was free from interactions and correlations with other domains in both forms of the repressor. A comparison of the simulations of the aporepressor and holorepressor showed that tryptophan binding made the DNA-binding helix D more flexible but helix F less flexible. Several amino acid residues in contact with the bound tryptophan were identified as making concerted motions with it. Interaction energies between the corepressor and the amino acid residues of the protein were analyzed; the results were mostly consistent with the mutational experiments.
Amino Acid Sequence↗
Molecular dynamics simulation study of DNA dodecamer d(CGCGAATTCGCG) in solution: conformation and hydration.
A molecular dynamics simulation of the dodecamer duplex d(CGCGAATTCGCG) using the particle mesh Ewald sum assumed a B-conformation remarkably close to the observed X-ray structure. The Ewald summation method effectively eliminates the usual "cut-off" of long-range interactions and allowed us to evaluate the full effect of the electrostatic forces. This simulation showed remarkable agreement with the Dickerson X-ray structure in both average structure and B-factors; within the EcoRI site itself, the rms deviation between the average theoretical and observed structures was 1.1 A. The width of the minor groove fluctuated between a wide and narrow configuration with the latter corresponding closely to the X-ray structure. The simulation also suggested a strong sequence-dependent signature on the minor groove width in both wide and narrow conformers. Hydration shells in both the major and minor grooves were observed. The "spine of hydration" in the minor groove was clear. In the major groove the first hydration shell appears to be a ribbon-like structure that reproduces the principal features of observed X-ray structures; subtle variations of this hydration pattern suggest sequence dependencies. Sequence-dependent features were also examined for helical and other geometric parameters. The successful reproduction of many experimentally observed fine structural features shows that the Ewald summation significantly improves the fidelity of the calculations.
Crystallography, X-Ray↗