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 901 records · Page 50Linked to original sources
PubMed2005

Molecular dynamics simulations of angiotensin II in aqueous and dimethyl sulfoxide environments.

Angiotensin II (Ang II) is an octapeptidic hormone, which plays an important role in the mechanisms of blood pressure control. In this work, extensive molecular dynamics (MD) simulations have been carried out on this peptide, both in aqueous and in dimethyl sulfoxide (DMSO) environments. Experimentally proposed models for the structure of angiotensin II in both environments are not consensual and the results obtained have provided some further insight about the structural properties of this hormone. In these simulations, the N-terminus of Ang II in the aqueous environment has been associated with a considerable larger flexibility than the correspondent C-terminus, but this was not found in the case of the DMSO environment. This is consistent with the assumption that the biological activity of Ang II is associated with its C-terminal residues embedded in a hydrophobic environment of its AT1 receptor. Other features detected in DMSO environment were an H(His6 imidazole)-O(Phe8 carboxylate) hydrogen bond and a salt-bridge structure involving the Asp1 and Arg2 side chains. An additional important conformational feature is the spatial proximity between Tyr4 and His6 in both water and DMSO environments. This molecular feature may trigger the interest for the synthetic chemists to apply rational design for the synthesis of novel AT1 antagonists.

Angiotensin II↗
PubMed2003

Triplex hydration: nanosecond molecular dynamics simulation of the solvated triplex formed by mixed sequences.

A theoretical model for the hydration pattern and motion of ions around the triple helical DNA with mixed sequences d(GACTGGTGAC)d(GTCACCAGTC)*d(GACTGGTGAC) in solution, during MD simulation, using the particle mesh Ewald sum method, is elaborated here. The AMBER 5.0 force field has been used during the simulation in solvent. The simulation studies support a dynamically stable atmosphere around the DNA triplex in solution over the entire length of the trajectory. The results have been compared with Hoogsteen triplexes and examined in the context of the observed behaviour of hydration in crystallographic data of duplexes. The dynamical organization of counterions and water molecules around the triplex formed by mixed sequences is described here. It has been observed that cations prefer to bind between two adjoining purines of the second and the third strands. The idea of localized complexes (mobile counterions in unspecific electronegative pockets around the DNA triplex with water molecules) may have important implications for understanding the specificity of the interactions of nucleic acids with proteins and other ligands.

Base Sequence↗
PubMed2004

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study.

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.

Computer Simulation↗
PubMed2005

A model for short alpha-neurotoxin bound to nicotinic acetylcholine receptor from Torpedo californica: comparison with long-chain alpha-neurotoxins and alpha-conotoxins.

Short-chain alpha-neurotoxins from snakes are highly selective antagonists of the muscle-type nicotinic acetylcholine receptors (nAChR). Although their spatial structures are known and abundant information on topology of binding to nAChR is obtained by labeling and mutagenesis studies, the accurate structure of the complex is not yet known. Here, we present a model for a short alpha-neurotoxin, neurotoxin II from Naja oxiana (NTII), bound to Torpedo californica nAChR. It was built by comparative modeling, docking and molecular dynamics using 1H NMR structure of NTII, cross-linking and mutagenesis data, cryoelectron microscopy structure of Torpedo marmorata nAChR [Unwin, N., 2005. Refined structure of the nicotinic acetylcholine receptor at 4A resolution. J. Mol. Biol. 346, 967-989] and X-ray structures of acetylcholine-binding protein (AChBP) with agonists [Celie, P.H., van Rossum-Fikkert, S.E., van Dijk, W.J., Brejc, K., Smit, A.B., Sixma, T.K., 2004. Nicotine and carbamylcholine binding to nicotinic acetylcholine receptors as studied in AChBP crystal structures. Neuron 41 (6), 907-914] and antagonists: alpha-cobratoxin, a long-chain alpha-neurotoxin [Bourne, Y., Talley, T.T., Hansen, S.B., Taylor, P., Marchot, P., 2005. Crystal structure of Cbtx-AChBP complex reveals essential interactions between snake alpha-neurotoxins and nicotinic receptors. EMBO J. 24 (8), 1512-1522] and alpha-conotoxin [Celie, P.H., Kasheverov, I.E., Mordvintsev, D.Y., Hogg, R.C., van Nierop, P., van Elk, R., van Rossum-Fikkert, S.E., Zhmak, M.N., Bertrand, D., Tsetlin, V., Sixma, T.K., Smit, A.B., 2005. Crystal structure of nicotinic acetylcholine receptor homolog AChBP in complex with an alpha-conotoxin PnIA variant. Nat. Struct. Mol. Biol. 12 (7), 582-588]. In complex with the receptor, NTII was located at about 30 A from the membrane surface, the tip of its loop II plunges into the ligand-binding pocket between the alpha/gamma or alpha/delta nAChR subunits, while the loops I and III contact nAChR by their tips only in a 'surface-touch' manner. The toxin structure undergoes some changes during the final complex formation (for 1.45 rmsd in 15-25 ps according to AMBER'99 molecular dynamics simulation), which correlates with NMR data. The data on the mobility and accessibility of spin- and fluorescence labels in free and bound NTII were used in MD simulations. The binding process is dependent on spontaneous outward movement of the C-loop earlier found in the AChBP complexes with alpha-cobratoxin and alpha-conotoxin. Among common features in binding of short- and long alpha-neurotoxins is the rearrangement of aromatic residues in the binding pocket not observed for alpha-conotoxin binding. Being in general very similar, the binding modes of short- and long alpha-neurotoxins differ in the ways of loop II entry into nAChR.

Amino Acid Sequence↗
PubMed2001

Molecular dynamics simulations of the three dimensional model of plasmepsin II-peptidic inhibitor complexes.

Plasmodium falciparum is a major causative agent of malaria, a disease of worldwide importance. Inhibition of a hemoglobin degrading P. falciparum aspartic protease Plasmepsin II (Plm II) provides a viable strategy for antimalarial therapy. Linear peptidic inhibitors based on the 4(S)-amino-3(S)-hydroxy-5-phenylpentanoic acid at the P1-P1' positions are known which inhibit Plm II with improved selectivity over cathepsin D. A series of computations were performed in order to gain insight into the interactions of these inhibitors with Plm II. The docking and molecular dynamics simulations were performed on a model ligand/enzyme complex to optimize the variables involved in the generation of ligand/enzyme models. This protocol of docking and molecular dynamics (MD) simulation was then used to derive the ligand-enzyme complexes of the molecules used in the present study. Different modes of binding of pepstatin and the three linear inhibitors were studied. Molecular dynamics simulation was performed at 300K for 100ps with a time step of Ifs. The structural effects of ligand binding were analyzed on the basis of hydrogen bond interactions, interaction energies, hydrophobic contacts and RMS deviations in the resulting energy-minimized structures of the receptor-ligand complexes. The results indicate that hydrophobic and hydrogen bonding interactions are responsible for selective inhibition of Plm II and improved selectivity over cathepsin D. Hydrogen bonding interaction plays an important role for amino acid residues such as Asp-34, Asp-214, Thr-217, Ser-218, Val-78, Ser-79, Tyr-192 and Gly-216. The binding of the inhibitors to the enzyme, while producing no large distortions in the enzyme active site cleft, results in significant RMS deviations of the inhibitor, which represent the distortion of the inhibitor, effected by the proteinase. Thus, the information generated from this analysis should be useful for further work in the area of antimalarial research.

Aspartic Acid Endopeptidases↗
PubMed2006

Low-temperature protein dynamics: a simulation analysis of interprotein vibrations and the boson peak at 150 k.

An understanding of low-frequency, collective protein dynamics at low temperatures can furnish valuable information on functional protein energy landscapes, on the origins of the protein glass transition and on protein-protein interactions. Here, molecular dynamics (MD) simulations and normal-mode analyses are performed on various models of crystalline myoglobin in order to characterize intra- and interprotein vibrations at 150 K. Principal component analysis of the MD trajectories indicates that the Boson peak, a broad peak in the dynamic structure factor centered at about approximately 2-2.5 meV, originates from approximately 10(2) collective, harmonic vibrations. An accurate description of the environment is found to be essential in reproducing the experimental Boson peak form and position. At lower energies other strong peaks are found in the calculated dynamic structure factor. Characterization of these peaks shows that they arise from harmonic vibrations of proteins relative to each other. These vibrations are likely to furnish valuable information on the physical nature of protein-protein interactions.

Cold Temperature↗
PubMed2004

Insight into G[bond]T mismatch recognition using molecular dynamics with time-averaged restraints derived from NMR spectroscopy.

Molecular dynamics (MD) simulations were conducted for a G[bond]T mismatch-containing DNA decamer, d(CCATGCGTGG)(2), and its Watson-Crick parent sequence, d(CCACGCGTGG)(2). Dynamics in unrestrained MD trajectories were in poor agreement with prior (13)C NMR studies. However, the accuracy of the trajectories was improved by the use of time-averaged interatomic distance restraints derived from (1)H NMR. Postprocess smoothing of the trajectories further improved accuracy. Comparison of restrained and smoothed trajectories of the two DNA molecules revealed distinct differences in dynamics. The major groove width of the mismatched oligomer was more variable over the course of the simulation compared to its parent sequence. Greater variability in helical parameters stretch and opening for the mismatches indicated less kinetically stable base pairing. Interbase helical parameters rise, roll, and tilt were also more variable in certain base steps involving mismatched bases. These dynamic differences between normal and G[bond]T mismatched DNA reflect differences in local flexibility that may play a role in mismatch recognition by the MutS. A potential alternate G[bond]T mismatch binding mode for MutS is also proposed.

Base Pair Mismatch↗
PubMed2003

8-oxoguanine lesioned B-DNA molecule complexed with repair enzyme hOGG1: a molecular dynamics study.

The molecular dynamics (MD) simulation of DNA mutagenic oxidative lesion, 7,8-dihydro-8-oxoguanine (8-oxoG), complexed with the repair enzyme, human oxoguanine glycosylase 1 (hOGG1), was performed for 1 nanosecond (ns) in order to describe the dynamical process of DNA-enzyme complex formation. After 900 picoseconds of MD the lesioned DNA and enzyme formed a complex that lasted until the end of the simulation at 1 ns. The complex was mainly represented by the overlapping van der Waals surfaces of DNA and enzyme molecules. The amino group of arginine 324 was located close to the phosphodiester bond of the nucleotide with 8-oxoG enabling chemical reactions between amino acid and lesion. The broken hydrogen bonds resulting in locally collapsed B-DNA structure were observed at the lesion site. The phosphodiester bond at C5' of 8-oxoG was displaced to the position close to the amino group of arginine 324. The water-mediated hydrogen bond network was formed in each contact area between DNA and enzyme, further enhancing the stability of the complex. In the background simulation of the identical molecular system with the native DNA, neither the complex nor the water- mediated hydrogen bond network was observed.

Algorithms↗
PubMed1998

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P.

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.

Amino Acid Sequence↗
PubMed1996

Structural consequences of heme removal: molecular dynamics simulations of rat and bovine apocytochrome b5.

Molecular dynamics simulations of rat and bovine apocytochrome b5 were performed to investigate the structural and dynamical consequences of heme removal. A crystal structure is available for the bovine holoprotein, while experimental studies of apocytochrome b5 have focused on the rat protein. The rat and bovine proteins are 93% homologous by sequence, and the sequence differences (six residues) appear to have no effect on the structure of the native holoprotein, as seen by the correlation of a bovine simulation with rat holocytochrome b5 experimental data (Storch & Daggett, 1995). There was a marked effect, however, on the structure and dynamics of the apo form. The bovine apocytochrome b5 simulation displayed subtle inconsistencies when compared to the experimental results on the rat apoprotein. Therefore, the rat protein was constructed from the bovine crystal structure coordinates. The MD simulation of the rat apoprotein displayed greater deviations from the crystal structure, yet it was in much closer agreement to the experimental data for the apoprotein. Additionally, the six variant residues fall in the regions where the bovine protein deviated from experiment. The two hydrophobic cores of the rat protein behaved very differently. Core 2 was well maintained, retained native-like structure, and is in good agreement with NMR data (Moore & Lecomte, 1990). Conversely, core 1, which is normally constrained by the prosthetic heme group, exhibited conformational heterogeneity, increased mobility, and some loss of secondary structure. Thus, the model of rat apocytochrome b5 complements past studies by providing structural information about core 1 that has proved difficult to obtain by experiment. The bovine simulation serves as a prediction, since little to no experimental data exist for this form of the apoprotein.

Animals↗
PubMed2005

Nanometer-scale ion aggregates in aqueous electrolyte solutions: guanidinium sulfate and guanidinium thiocyanate.

Neutron diffraction experiments and molecular dynamics simulations are used to study the structure of aqueous solutions of two electrolytes: guanidinium sulfate (a mild protein conformation stabilizer) and guanidinium thiocyanate (a powerful denaturant). The MD simulations find the unexpected result that in the Gdm2SO4 solution the ions aggregated into mesoscopic (nanometer-scale) clusters, while no such aggregation is found in the GdmSCN solution. The neutron diffraction studies, the most direct experimental probe of solution structure, provide corroborating evidence that the predicted very strong ion pairing does occur in solutions of 1.5 m Gdm2SO4 but not in 3 m solutions of GdmSCN. A mechanism is proposed as to how this mesoscopic solution structure affects solution denaturant properties and suggests an explanation for the Hofmeister ordering of these solutions in terms of this ion pairing and the ability of sulfate to reverse the denaturant power of guanidinium.

Electrolytes↗
PubMed2003

Dynamics-driven reaction pathway in an intramolecular rearrangement.

A critical role is traditionally assigned to transition states (TSs) and minimum energy pathways, or intrinsic reaction coordinates (IRCs), in interpreting organic reactivity. Such an interpretation, however, ignores vibrational and kinetic energy effects of finite temperature. Recently it has been shown that reactions do not necessarily follow the intermediates along the IRC. We report here molecular dynamics (MD) simulations that show that dynamics effects may alter chemical reactions even more. In the heterolysis rearrangement of protonated pinacolyl alcohol Me3C-CHMe-OH2+ (Me, methyl), the MD pathway involves a stepwise route with C-O bond cleavage followed by methyl group migration, whereas the IRC pathway suggests a concerted mechanism. Dynamics effects may lead to new interpretations of organic reactivity.

Journal Article↗
PubMed2004

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Binding Sites↗
PubMed2005

Ab initio QM/MM simulation of Ag+ in 18.6% aqueous ammonia solution: structure and dynamics investigations.

Structure and dynamics investigations of Ag(+) in 18.6% aqueous ammonia solution have been carried out by means of the ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation method. The most important region, the first solvation shell, was treated by ab initio quantum mechanics at the Restricted Hartree-Fock (RHF) level using double-zeta plus polarization basis sets for ammonia and plus ECP for Ag(+). For the remaining region in the system, newly constructed three-body corrected potential functions were used. The average composition of the first solvation shell was found to be [Ag(NH(3))(2)(H(2)O)(2.8)](+). No ammonia exchange process was observed for the first solvation shell, whereas ligand exchange processes occurred with a very short mean residence time of 1.1 ps for the water ligands. No distinct second solvation shell was observed in this simulation.

Ammonia↗
PubMed1996

Hydration and DNA recognition by homeodomains.

A 2-nanosecond molecular dynamics (MD) simulation of an Antennapedia homeodomain-DNA complex in explicit solvent water at ambient temperature and pressure was performed to supplement experimental nuclear magnetic resonance (NMR) data on the structure and dynamics of this complex. In addition to direct protein-DNA contacts, the MD trajectory attributes an essential role for specific DNA recognition to hydration water molecules that mediate intermolecular contacts. The simulation provides a detailed description of the pathways of hydration water molecules exchanging in and out of the protein-DNA interface and indicates that the residence times of these "interior" waters are on the nanosecond time scale, near the lower end of the range determined by NMR.

Amino Acids↗
PubMed1999

Molecular dynamics study of kaliotoxin in water.

Kaliotoxin (KTX), a potassium channel blocker found in the venom of the scorpion Androctonous Mauretanicus is a 38 residue polypeptide with a well defined structure consisting of a alpha-helix and a three strand antiparallel beta-sheet interconnected by three disulfide bonds. Although the 3D structure has been determined by NMR, there is a number of features, mainly concerning the conformation and flexibility of the side chains, but also the long range order in the peptide and its fluctuations, that may have escaped the experimental study. These questions are now being addressed using molecular dynamics (MD) simulations. Accordingly, the present work reports the analysis of a 430 ps molecular dynamics trajectory of the polypeptide soaked with 4700 TIP3 water molecules inside a 56 A box. MD calculations were performed with periodic boundary conditions. Analysis of the conformational space sampled by each of the residues along the trajectory, suggests a special behavior of Pro17 and Lys19 both located on the helix. Furthermore, analysis of the relative movements of the secondary structure elements indicates that the alpha-helix and beta-sheets fluctuate in a correlated motion, preserving the tertiary structure of the polypeptide along the trajectory. Finally, analysis of the charge distribution was also examined. The direction of the dipole moment, computed from the center of masses appears to be an interesting feature of the structure probably related to the biological function of the molecule.

Computer Simulation↗
PubMed2000

A three-dimensional model of the delta-opioid pharmacophore: comparative molecular modeling of peptide and nonpeptide ligands.

A comparative molecular modeling study of delta-opioid ligands was performed under the assumption that potent peptide and nonpeptide agonists may have common three-dimensional (3D) arrangement of pharmacophore groups upon binding to the delta-receptor. Low-energy conformations of the agonists 7-spiroindanyloxymorphone (SIOM) and 2-methyl-4a-alpha-(3-hydroxyphenyl)-1,2,3,4,4a,5,12, 12a-alpha-octahydro-quinolino[2,3,3-g]isoquinoline (TAN-67), and a partial agonist oxomorphindole (OMI) were determined by high-temperature molecular dynamics (MD). A good spatial overlap was found for the pharmacophore groups of SIOM, TAN-67, and OMI, including the basic nitrogen, phenol hydroxyl, and two aromatic ring. Based on this overlap we proposed a 3D pharmacophore model for nonpeptide delta-opioid agonists with a distance of 7.0 +/- 1.3 A between the two aromatic rings and of 8.2 +/- 1.0 A between the nitrogen and phenyl ring. The potent and highly delta-opioid receptor selective agonist [(2S,3R)-TMT(1)]DPDPE, which shares global backbone constraints of the 14-membered disulfide cycle and a strong preference for the trans rotamer of the TMT(1) side chain, was chosen as a peptide template of the delta-opioid pharmacophore. Extensive MD simulations at 300 K with the AMBER force field were performed for [(2S,3R)-TMT(1)]DPDPE and the less potent [(2S, 3S)-TMT(1)]DPDPE analogue. Multiple MD trajectories were collected for each peptide starting from the x-ray structures of DPDPE and [L-Ala(3)]DPDPE and from models proposed in the literature. Low-energy MD conformations were filtered by the nonpeptide pharmacophore query and then directly superimposed with SIOM, OMI, and TAN-67. Two conformers of [(2S,3R)-TMT(1)]DPDPE that showed the best overlap with the nonpeptide pharmacophore (rms deviation </= 1. 0 A for N,O atoms and centroids of two aromatic rings) were selected as possible delta-receptor binding conformations. These conformations have similar backbone structures, and trans rotamers of the TMT(1) side-chain group. They are reasonably close to the crystal structure of [L-Ala(3)]DPDPE, and differ significantly from the crystal structure of DPDPE. The conformer with a gauche(-) rotamer of Phe(4) is most consistent with structure-activity relationships of delta-opioid peptides. The proposed 3D models were used for rational design of new nonpeptide delta-receptor ligands.

Analgesics↗
PubMed2005

The effect of salt concentration on DNA conformation transition: a molecular-dynamics study.

We performed three 3-ns molecular dynamics simulations of d(CGCGAATTCGCG)2 using the AMBER 8 package to determine the effect of salt concentration on DNA conformational transitions. All the simulations were started with A-DNA, with different salt concentrations, and converged with B-DNA with similar conformational parameters. However, the dynamic processes of the three MD simulations were very different. We found that the conformation transition was slow in the solution with higher salt concentration. To determine the cause of this retardation, we performed three additional 1.5-ns simulations starting with B-DNA and with the salt concentrations corresponding to the simulations mentioned above. However, astonishingly, there was no delayed conformation evolution found in any of the three simulations. Thus, our simulation conclusion is that higher salt concentrations slows the A --> B conformation transition, but have no effect on the final stable structure. [Figure: see text].

Carbohydrates↗