Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources
PubMed2005

Crystal water dynamics of guanosine dihydrate: analysis of atomic displacement parameters, time profile of hydrogen-bonding probability, and translocation of water by MD simulation.

The dynamics of crystal water molecules of guanosine dihydrate are investigated in detail by molecular dynamics (MD) simulation. A 2 ns simulation is performed using a periodic boundary box composed of 4 x 5 x 8 crystallographic unit cells and using the particle-mesh Ewald method for calculation of electrostatic energy. The simulated average atomic positions and atomic displacement parameters are remarkably coincident with the experimental values determined by X-ray analysis, confirming the high accuracy of this simulation. The dynamics of crystal water are analyzed in terms of atomic displacement parameters, orientation vectors, order parameters, self-correlation functions of the orientation vectors, time profiles of hydrogen-bonding probability, and translocations. The simulation clarifies that the average structure is composed of various stable and transient structures of the molecules. The simulated guanosine crystal forms a layered structure, with four water sites per asymmetric unit, classified as either interlayer water or intralayer water. From a detailed analysis of the translocations of water molecules in the simulation, columns of intralayer water molecules along the c axis appear to represent a pathway for hydration and dehydration by a kind of molecular valve mechanism.

Computer Simulation↗
PubMed2006

Spectral signatures of heterogeneous protein ensembles revealed by MD Simulations of 2DIR spectra.

A model for the calculation of amide I FTIR and 2DIR spectra taking into account fluctuations in hydrogen bonding and structure from molecular dynamics (MD) simulations is tested on three systems. It is found that although the homogeneous lineshape approximation yields satisfactory FTIR spectra, 2DIR spectra are sensitive to the inhomogeneity naturally present in any solvated protein and the common approximations of a static structure and averaged-effect solvent are invalid. By building on the local amide Hamiltonian and incorporating site energy variation with electrostatic-based models and disorder from MD trajectories, good agreement is obtained between calculated and measured 2DIR spectra. The largest contribution to the observed inhomogeneity is found to be the fluctuating site energies, which in turn are most sensitive to the water solvent. With the ability to accurately predict 2DIR spectra from atomistic simulations, new opportunities to test force fields and mechanistic predictions from MD are revealed.

Computer Simulation↗
PubMed2005

MD simulations of the binding of alcohols and diols by a calixarene in water: connections between microscopic and macroscopic properties.

We report results of molecular dynamics (MD) simulations of the complexes of p-sulfonatocalix[4]arene with linear alcohols from ethanol to heptanol in water at 25 degrees C. We show that these complexes are of the inclusion type and are governed by van der Waals interactions between the calixarene cavity and the inserted alkyl chain of the alcohol. We establish a correlation between the experimental Delta(r)H degrees values and the number of atoms inserted into the calixarene cavity. We also focus on the desolvation of the host and guest to establish the importance, at the enthalpic level, of the formation of hydrogen bond bridges between the calixarene and the alcohol molecule. The fact that methanol is not complexed by p-sulfonatocalix[4]arene is explained by calculating the cost of the desolvation of the guest upon complexation. We complete this study by modeling the complexes formed with 1,4-butanediol and 1,5-pentanediol. To explain the difference between the thermodynamic properties for the binding of 1,4-butanediol and butanol, we examine the insertion rate and the solvation of each hydroxy group. We show a specific behavior of one of the two hydroxy groups at the structural and energetic levels.

Alcohols↗
PubMed2001

The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation.

The structural and thermodynamic properties of a 6-residue beta-peptide that was designed to form a hairpin conformation have been studied by NMR spectroscopy and MD simulation in methanol solution. The predicted hairpin would be characterized by a 10-membered hydrogen-bonded turn involving residues 3 and 4, and two extended antiparallel strands. The interproton distances and backbone torsional dihedral angles derived from the NMR experiments at room temperature are in general terms compatible with the hairpin conformation. Two trajectories of system configurations from 100-ns molecular-dynamics simulations of the peptide in solution at 298 and 340 K have been analyzed. In both simulations reversible folding to the hairpin conformation is observed. Interestingly, there is a significant conformational overlap between the unfolded state of the peptide at each of the temperatures. As already observed in previous studies of peptide folding, the unfolded state is composed of a (relatively) small number of predominant conformers and in this case lacks any type of secondary-structure element. The trajectories provide an excellent ground for the interpretation of the NMR-derived data in terms of ensemble averages and distributions as opposed to single-conformation interpretations. From this perspective, a relative population of the hairpin conformation of 20% to 30% would suffice to explain the NMR-derived data. Surprisingly, however, the ensemble of structures from the simulation at 340 K reproduces more accurately the NMR-derived data than the ensemble from the simulation at 298 K, a question that needs further investigation.

Magnetic Resonance Spectroscopy↗
PubMed2004

Molecular dynamics (MD) simulations for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD) based on the MM-PBSA (molecular mechanics-Poisson-Boltzmann surface area) approach.

Molecular dynamics (MD) simulations were performed for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD). Binding free energies and various conformational properties were obtained using by the MM-PBSA (molecular mechanics Poisson-Boltzmann/surface area) approach. The calculated relative difference (DeltaDeltabinding) of binding free energy was in fine agreement with the experimentally determined value. The difference of rotameric distributions of guest N-acetylphenylalanine enantiomers complexed with the host, beta-CD, was observed after the conformational analyses, suggesting that the conformational changes of guest captured within host cavity would be a decisive factor for enantiodifferentiation at a molecular level.

Computer Simulation↗
PubMed2005

Cooperative binding of DNA and CBFbeta to the Runt domain of the CBFalpha studied via MD simulations.

The Runt domain (RD) is the DNA-binding region of the Runx genes. A related protein, known as core binding factor beta (CBFbeta) also binds to the RD to enhance RD-DNA interaction by 6- to 10-fold. Here, we report results from molecular dynamics (MD) simulations of RD alone, as a dimer in complexes with DNA and CBFbeta and in a ternary complex with DNA and CBFbeta. Consistent with the experimental findings, in the presence of CBFbeta the estimated free energy of binding of RD to the DNA is more favorable, which is shown to be due to more favorable intermolecular interactions and desolvation contributions. Also contributing to the enhanced binding are favorable intramolecular interactions between the 'wing' residues (RD residues 139-145) and the 'wing1' residues (RD residues 104-116). The simulation studies also indicate that the RD-CBFbeta binding is more favorable in the presence of DNA due to a more favorable RD-CBFbeta interaction energy. In addition, it is predicted that long-range interactions involving ionic residues contribute to binding cooperativity. Results from the MD calculations are used to interpret a variety of experimental mutagenesis data. A novel role for RD Glu116 to the RD-CBFbeta interaction is predicted.

Amino Acids↗
PubMed2006

Orientation and conformational preference of leucine-enkephalin at the surface of a hydrated dimyristoylphosphatidylcholine bilayer: NMR and MD simulation.

The morphogenic opiate pentapeptide leucine-enkephalin (lenk) in a hydrated dimyristoylphosphatidylcholine (DMPC) bilayer is studied using NMR spectroscopy and molecular dynamics simulation. Contrary to the frequent assumption that the peptide attains a single fixed conformation in the presence of membranes, we find that the lenk molecule is flexible, switching between specific bent conformations. The constraints to the orientation of the aromatic rings that are identified by the NMR experiment are found by the MD simulation to be related to the depth of the peptide in the bilayer. The motion of the N-H vectors of the peptide bonds with respect to the magnetic field direction as observed by MD largely explain the magnitude of the observed residual dipolar coupling (RDC), which are much reduced over the static (15)N-(1)H coupling. The measured RDCs are nevertheless significantly larger than the predicted ones, possibly due the absence of long-time motions in the simulations. The conformational behavior of lenk at the DMPC surface is compared to that in the aqueous solution, both in the neutral and in the zwitterionic forms.

Computer Simulation↗
PubMed2002

General framework for studying the dynamics of folded and nonfolded proteins by NMR relaxation spectroscopy and MD simulation.

A general framework is presented for the interpretation of NMR relaxation data of proteins. The method, termed isotropic reorientational eigenmode dynamics (iRED), relies on a principal component analysis of the isotropically averaged covariance matrix of the lattice functions of the spin interactions responsible for spin relaxation. The covariance matrix, which is evaluated using a molecular dynamics (MD) simulation, is diagonalized yielding reorientational eigenmodes and amplitudes that reveal detailed information about correlated protein dynamics. The eigenvalue distribution allows one to quantitatively assess whether overall and internal motions are statistically separable. To each eigenmode belongs a correlation time that can be adjusted to optimally reproduce experimental relaxation parameters. A key feature of the method is that it does not require separability of overall tumbling and internal motions, which makes it applicable to a wide range of systems, such as folded, partially folded, and unfolded biomolecular systems and other macromolecules in solution. The approach was applied to NMR relaxation data of ubiquitin collected at multiple magnetic fields in the native form and in the partially folded A-state using MD trajectories with lengths of 6 and 70 ns. The relaxation data of native ubiquitin are well reproduced after adjustment of the correlation times of the 10 largest eigenmodes. For this state, a high degree of separability between internal and overall motions is present as is reflected in large amplitude and collectivity gaps between internal and overall reorientational modes. In contrast, no such separability exists for the A-state. Residual overall tumbling motion involving the N-terminal beta-sheet and the central helix is observed for two of the largest modes only. By adjusting the correlation times of the 10 largest modes, a high degree of consistency between the experimental relaxation data and the iRED model is reached for this highly flexible biomolecule.

Computer Simulation↗
PubMed1999

Conformations of trypsin-bound amidine inhibitors of blood coagulant factor Xa by double REDOR NMR and MD simulations.

Double rotational-echo double resonance (double REDOR) has been used to investigate the bound conformations of (13)C,(15)N,(19)F-labeled factor Xa inhibitors to bovine trypsin. Carbon-fluorine dipolar couplings were measured by (13)C{(19)F} REDOR with natural-abundance background interferences removed by (13)C{(15)N} REDOR. The conformations of the bound inhibitors were characterized by molecular dynamics (MD) simulations of binding restrained by double REDOR-determined intramolecular C-F distances. A symmetrical bisamidine inhibitor and an asymmetrical monoamidine-monoamine inhibitor of the same general shape had distinctly different conformations in the bound state. According to the MD models, these differences arise from specific interactions of the amidine and amine groups with the active-site residues of trypsin and nearby water molecules.

Animals↗
PubMed2000

Conformational analysis of the Sd(a) determinant-containing tetrasaccharide and two mimics in aqueous solution by using 1H NMR ROESY spectroscopy in combination with MD simulations.

The conformational behaviour of the spacer-linked synthetic Sd(a) tetrasaccharide beta-D-GalpNAc-(1-->4)-[alpha-Neu5Ac-(2-->3)]-beta-D-Galp-(1-->4)- beta-D-GlcpNAc-(1-->O)(CH2)5NH2 (1) and the two mimics beta-D-Galp-(1-->4)-[alpha-Neu5Ac-(2-->3)]-beta-D-Galp-(1-->4)-bet a-D-GlcpNAc-(1-->O)(CH2)5NH2 (2) and beta-D-GlcpNAc-(1-->4)-[alpha-Neu5Ac-(2-->3)]-beta-D-Galp-(1-->4)- beta-D-GlcpNAc-(1-->O)(CH2)5NH2 (3) were investigated by 1H NMR spectroscopy in combination with molecular dynamics (MD) simulations in water. Experimental 2D 1H ROESY cross-peak intensities (ROEs) of the tetrasaccharides were compared with calculated ROEs derived from MD trajectories using the CROSREL program. Analysis of these data indicated that the oligosaccharidic skeletons of the compounds 1-3 are rather rigid, especially the beta-D-Hex(NAc)-(1-->4)-[alpha-Neu5Ac-(2-->3)]-beta-D-Galp fragments. The alpha-Neu5-Ac-(2-->3)-beta-D-Galp linkage occurred in two different energy minima in the three-dimensional structure of the compounds 1-3 in aqueous solution. Experimental data and dynamics simulations supported the finding that the higher energy rotamer (CHEAT forcefield) was abundant in compounds 1 and 3 due to the existence of a hydrogen bond between the carboxyl group of the sialic acid and the acetamido group of the terminal monosaccharide (GalNAc or GlcNAc) unit. The conformational similarity between 1 and 3 leads to the suggestion that also their activities will be alike.

Carbohydrate Conformation↗
PubMed2006

Homology modeling of human serum carnosinase, a potential medicinal target, and MD simulations of its allosteric activation by citrate.

Recent biochemical and clinical evidence implicates human serum carnosinase in a variety of pathological conditions, such as neurological disorders and diabetic nephropathy, suggesting that this enzyme is of potential interest as a novel medicinal target. The present study was undertaken with a view to model the serum carnosinase and its catalytic site and to unravel the molecular mechanism by which citrate ions increase the catalytic efficiency of serum carnosinase. A homology model of the enzyme was obtained on the basis of beta-alanine synthetase, and its active center was found to bind known substrates carnosine, homocarnosine, and anserine in a binding mode conducive to catalysis. Citrate ions were shown to bind at only three well-defined sites involving both ion pairs and hydrogen bonds. Molecular dynamics simulations evidenced that citrate binding had a remarkable conformational influence on the 3D structure of carnosinase, increasing the binding affinity (i.e., binding score) of carnosine to the catalytic site. This is one of the first reports documenting the molecular mechanism of an allosteric enzyme activator using MD simulations.

Allosteric Regulation↗
PubMed2002

Conformational analysis of two xylose-containing N-glycans in aqueous solution by using 1H NMR ROESY and NOESY spectroscopy in combination with MD simulations.

The conformational behavior of the synthetic hexa- and heptasaccharide methyl beta-glycosides alpha-D-Manp-(1 --> 6)-[alpha-D-Manp-(1 --> 3)-][beta-D-Xylp-(1 --> 2)-]beta-D-Manp-(1 --> 4)-beta-D-GlcpNAc-(1 --> 4)-beta-D-GlcpNAc-(1 --> OMe and alpha-D-Manp-(1 --> 6)-[alpha-D-Manp-(1 --> 3)-][beta-D-Xylp-(1 --> 2)-]beta-D-Manp-(1 --> 4)-beta-D-GlcpNAc-(1 --> 4)-[alpha-L-Fucp-(1 --> 6)-]beta-D-GlcpNAc-(1 --> OMe, representing the xylosylated and the xylosylated alpha-(1 --> 6)-fucosylated core structures of N-glycans in alpha(D)-hemocyanin of the snail Helix pomatia, respectively, were investigated by 1H NMR spectroscopy in combination with molecular dynamics (MD) simulations in water. 1H and 13C chemical shifts of the oligosaccharides were assigned using 1H-(1)H COSY, TOCSY, and NOESY, and 1H-(13)C HMQC techniques. Experimental 2D 1H cross-peak intensities from one series of NOESY and one series of ROESY experiments of the two oligosaccharides were compared with calculated values derived from MD trajectories using the CROSREL program, yielding information about the conformation of each glycosidic linkage of the methyl glycosides. The flexibility of the linkages was described by generalized order parameters and internal rotation correlation times. Analysis of the data indicated that several conformations are likely to exist for the alpha-D-Man-(1 --> 6)-beta-D-Man, the alpha-L-Fuc-(1 --> 6)-beta-D-GlcNAc, and the alpha-D-Man-(1 --> 3)-beta-D-Man linkage, whereas the beta-D-Xyl-(1 --> 2)-beta-D-Man-(1 --> 4)-beta-D-GlcNAc-(1 --> 4)-beta-D-GlcNAc fragment occurs in one rigid conformation. No significant differences were found between the corresponding structural elements in both methyl glycosides. NOESY and ROESY experiments proved to be suitable for providing the experimental data required, however, due to more overlap within the ROESY spectra, reducing the accuracy of the analysis, NOESY spectral analysis is preferred.

Animals↗
PubMed2005

Solution conformation of Substance P antagonists-[D-Arg1, D-Trp7,9, Leu11]-SP, [D-Arg1, D-Pro2, D-Trp7,9, Leu11]-SP and [D-Pro2, D-Trp7,9]-SP by CD, NMR and MD simulations.

Substance P (SP) is an important neuropeptide involved in pain transmission and induction of inflammation. Its antagonists are being extensively investigated for their non-narcotic analgesic and anti-inflammatory activity. With a view towards better understanding the structural requirements of these analogs for efficient interaction with the SP receptor, the conformation of three SP antagonists [D-Arg1, D-Trp7,9, Leu11]-SP, [D-Arg1, D-Pro2, D-Trp7,9, Leu11]-SP and [D-Pro2, D-Trp7,9]-SP has been studied by CD, NMR and molecular dynamics (MD) simulations. All three peptides exhibit a high dependence of structure on the solvent. The molecules tend to adopt beta-turns in solvents like DMSO and H2O and form helices in a hydrophobic environment. A direct relation between the helix forming potential of these antagonists with their receptor binding potency has been observed.

Circular Dichroism↗
PubMed2005

A combination of docking, QM/MM methods, and MD simulation for binding affinity estimation of metalloprotein ligands.

To alleviate the problems in the receptor-based design of metalloprotein ligands due to inadequacies in the force-field description of coordination bonds, a four-tier approach was devised. Representative ligand-metalloprotein interaction energies are obtained by subsequent application of (1) docking with metal-binding-guided selection of modes, (2) optimization of the ligand-metalloprotein complex geometry by combined quantum mechanics and molecular mechanics (QM/MM) methods, (3) conformational sampling of the complex with constrained metal bonds by force-field-based molecular dynamics (MD), and (4) a single point QM/MM energy calculation for the time-averaged structures. The QM/MM interaction energies are, in a linear combination with the desolvation-characterizing changes in the solvent-accessible surface areas, correlated with experimental data. The approach was applied to structural correlation of published binding free energies of a diverse set of 28 hydroxamate inhibitors to zinc-dependent matrix metalloproteinase 9 (MMP-9). Inclusion of steps 3 and 4 significantly improved both correlation and prediction. The two descriptors explained 90% of variance in inhibition constants of all 28 inhibitors, ranging from 0.08 to 349 nM, with the average unassigned error of 0.318 log units. The structural and energetic information obtained from the time-averaged MD simulation results helped understand the differences in binding modes of related compounds.

Binding Sites↗
PubMed2005

MD simulation of the Na+-phenylalanine complex in water: competition between cation-pi interaction and aqueous solvation.

The competition between cation-pi interaction and aqueous solvation for the Na+ ion has been investigated by molecular dynamics simulations, using the phenylalanine amino acid as the test pi system. Starting from one of the best standard force fields, we have developed new parameters that significantly improve the agreement with experimental and high quality quantum mechanical results for the complexes of Na+ with phenylalanine, benzene, and water. The modified force field performs very well in forecasting energy and geometry of cation coordination for the complexes. Next, analysis of MD trajectories and steered MD simulations indicate that the Na+-phenylalanine complex survives for a significant time in aqueous solution and that the free energy barrier opposing dissociation of the complex is sizable. Finally, we analyze the role of different intermolecular interactions in determining the preference for cation-pi bonding with respect to aqueous solvation. We thus confirm that the Na+-phenylalanine stabilization energy may overcome the interactions with water.

Journal Article↗
PubMed2004

Conformational study of fragments of envelope proteins (gp120: 254-274 and gp41: 519-541) of HIV-1 by NMR and MD simulations.

The envelope proteins, gp 120 and gp41 of HIV-1, play a crucial role in receptor (CD4+ lymphocytes) binding and membrane fusion. The fragment 254-274 of gp120 is conserved in all strains of HIV and, as a part of the full gp120 protein, behaves as 'immunosilent', but as an individual fragment it is 'immunoreactive'. When this fragment binds to its receptor, it activates the fusion domain of gp41 allowing viral entry into the host CD4+ cells. The conformation of fragment 254-274 of the gp120 domain and fragment 519-541 of the gp41 domain was studied by NMR and MD simulations. The studies were carried out in three varied media--water, DMSO-d6 and hexafluoroacetone (HFA). The fusogenic nature of the gp41 domain peptide was investigated by 31P NMR experiments with model bilayers prepared from dimyristoyl-L-alpha-phosphatidylcholine (DMPC). The solvent was seen to exert a major effect on the structure of the two peptides. Fragment (254-274) of gp120 in DMSO-d6 had a type I beta-turn around the tetrad Val9-Ser10-Thr11-Gln12 while in HFA a helical structure spanning the region Ile5 to Gln12 was seen with the remaining part of the peptide in a random coil structure. It is possible that the beta-turn may constitute an initiation site for the formation of the helix. In water at pH 4.5, the peptide adopted a beta-sheet. The NMR results for fragment 519-541 of gp41 are conclusive of a beta-sheet structure in DMSO-d6, a conformation which may help in insertion into the membrane, a notion also put forward by others. The 31P NMR studies of DMPC vesicles with this fragment show its fusogenic nature, promoting fusion of unilamellar vesicles to larger agglomerates like multilamellar ones.

Amino Acid Sequence↗
PubMed2006

Comparative unfolding studies of psychrophilic and mesophilic uracil DNA glycosylase: MD simulations show reduced thermal stability of the cold-adapted enzyme.

Uracil DNA glycosylase (UDG) is a DNA repair enzyme involved in the base excision repair (BER) pathway, removing misincorporated uracil from the DNA strand. The native and mutant forms of Atlantic cod and human UDG have previously been characterized in terms of kinetic and thermodynamic properties as well as the determination of several crystal structures. This data shows that the cold-adapted enzyme is more catalytically efficient but at the same time less resistant to heat compared to its warm-active counterpart. In this study, the structure-function relationship is further explored by means of comparative molecular dynamics (MD) simulations at three different temperatures (375, 400 and 425K) to gain a deeper insight into the structural features responsible for the reduced thermostability of the cold-active enzyme. The simulations show that there are distinct structural differences in the unfolding pathway between the two homologues, particularly evident in the N- and C-terminals. Distortion of the mesophilic enzyme is initiated simultaneously in the N- and C-terminal, while the C-terminal part plays a key role for the stability of the psychrophilic enzyme. The simulations also show that at certain temperatures the cold-adapted enzyme unfolds faster than the warm-active homologues in accordance with the lower thermal stability found experimentally.

Acclimatization↗
PubMed2025

Computational identification of potential antifungal targets against Claviceps purpurea via MD simulation and MM/GBSA.

Ergot alkaloids produced by the fungus Claviceps purpurea pose significant risks to agriculture and human health. This study systematically investigates the pathogenicity of C. purpurea, analyzing five strains for their proteomic profiles, which revealed genetic variability in size and GC content. We identified proteins localized in various cellular compartments, contributing to our understanding of essential cellular processes. A focus on potential drug targets led to the identification of Alpha-N-acetylglucosaminidase, a hydrolase with significant mass and functional relevance, despite not matching a UniProt entry. The 3D structure prediction confirmed its integrity, making it a suitable target for further analysis. Molecular docking identified ligands CID:51,535,944 and CID:145,242,255 with strong binding affinities to Alpha-N-acetylglucosaminidase, highlighting interactions with key residues like TRP138, ARG651. Molecular docking interactions were validated and showed consistency through MD simulation analyses with greater RMSD, RMSF and PL contacts. This research enhances our understanding of C. purpurea, offering insights into its genetic diversity and cellular mechanisms while identifying promising therapeutic targets. The findings contribute to strategies for mitigating the economic and health impacts of C. purpurea infections, paving the way for innovative interventions in sustainable agriculture.

Molecular Docking Simulation↗