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 793 records · Page 44Linked to original sources
PubMed1999

Molecular dynamics study of substance P peptides partitioned in a sodium dodecylsulfate micelle.

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) during the equilibration stage of the simulation. It was shown that when SP-Y8 was initially placed in an insertion (perpendicular) configuration, the peptide equilibrated to a surface-bound (parallel) configuration in approximately 450 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained from the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbones of all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions come from the hydrophobic interaction between the side chains of Lys-3, Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, from Lys-3, Phe-7, Leu-10, and Met-11 in SP-Y8, and the micellar interior. Significant interactions, electrostatic and hydrogen bonding, between the N-terminal residues, Arg-Pro-Lys, and the micellar headgroups were observed. These latter interactions served to affect both the structure and, especially, the flexibility, of the N-terminus. The results from simulation of the same peptides in a water/CCl4 biphasic cell were compared with the results of the present study, and the validity of using the biphasic system as an approximation for peptide-micelle or peptide-bilayer systems is discussed.

Amino Acid Sequence↗
PubMed2003

Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets.

Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg2+ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na+ ions with residency times of individual cations 1-2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.

Binding Sites↗
PubMed2006

Neutron diffraction and simulation studies of the exocyclic hydroxymethyl conformation of glucose.

The techniques of neutron diffraction with isotopic substitution (NDIS) and molecular dynamics (MD) simulations have been used to examine the rotational conformation of the exocyclic hydroxymethyl group of D-glucopyranose. First order H/D NDIS experiments were performed on the H6 position in 3m aqueous glucose solutions where the average coherent scattering length of the exchangeable hydrogen atoms was zero (i.e., all correlations between exchangeable hydrogen atoms and other atoms cancel and thus are not present in the scattering data). This H6 experimental result suggests that no single conformation for the C4-C5-C6-O6 dihedral reproduces the observed scattering data well, but that a mixture of the gg and gt conformations, which has been suggested by NMR experiments, gives a reasonable agreement between the MD and experimental data.

Computer Simulation↗
PubMed1996

Conformational studies of microcystin-LR using NMR spectroscopy and molecular dynamics calculations.

NMR spectroscopy in aqueous and dimethyl sulfoxide/water solutions is used to determine the three-dimensional structures of microcystin-LR, a cyclic cyanobacterial heptapeptide toxin which is a potent inhibitor of type 1 and type 2A protein phosphatases. The conformations of this toxic peptide are studied using a simulated annealing (SA) protocol followed by refined SA calculations in vacuo and free MD simulations in water. Only one conformational family in each solvent is found. The peptide ring has a saddle-shaped form, essentially the same in both solvents. The structural difference observed between the two solution structures is located to the part consisting of Mdha, Ala, and Leu. This peptide segment is not present in nodularin, a cyclic pentapeptide of similar toxicity. The Arg side chain is very flexible, while the side chain of Leu is well defined. The side chain of Adda, essential for toxicity, is constrained in the vicinity of the backbone ring but appears to be flexible in the more remote part.

Bacterial Toxins↗
PubMed2005

Solvation dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation.

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.

Journal Article↗
PubMed1991

Are time-averaged restraints necessary for nuclear magnetic resonance refinement? A model study for DNA.

A recently suggested method for refinement of structural data obtained from two-dimensional nuclear magnetic resonance experiments using molecular dynamics (MD) is explored. In this method, the time-averaged values of the appropriate internal co-ordinates of the molecule, calculated from the MD trajectory, are driven by restraints towards the experimental target values. This contrasts with most refinement procedures currently in use, where restraints are applied based on the instantaneous values of the appropriate co-ordinates. Both refinement methods are applied to the EcoRI restriction site DNA hexamer d(GAATTC)2, using target nuclear Overhauser enhancement distances derived from a one nanosecond unrestrained MD simulation of this structure. The resulting refined structures are compared to the results of the unrestrained MD trajectory, which serves as our "experimental" data. We show that although both methods can yield an average structure with the correct gross morphology, the new method allows both a much more realistic picture of inherent flexibility, and reproduces fine conformational detail better, such as sequence dependency. We also analyze the very long MD trajectory generated here (longer than any previously reported for a DNA oligomer), and find that significantly shorter simulations, typical of those frequently performed, may not yield acceptably reliable values for certain structural parameters.

DNA↗
PubMed2003

Molecular dynamics simulation of highly charged proteins: comparison of the particle-particle particle-mesh and reaction field methods for the calculation of electrostatic interactions.

Molecular dynamics (MD) simulations of the activation domain of porcine procarboxypeptidase B (ADBp) were performed to examine the effect of using the particle-particle particle-mesh (P3M) or the reaction field (RF) method for calculating electrostatic interactions in simulations of highly charged proteins. Several structural, thermodynamic, and dynamic observables were derived from the MD trajectories, including estimated entropies and solvation free energies and essential dynamics (ED). The P3M method leads to slightly higher atomic positional fluctuations and deviations from the crystallographic structure, along with somewhat lower values of the total energy and solvation free energy. However, the ED analysis of the system leads to nearly identical results for both simulations. Because of the strong similarity between the results, both methods appear well suited for the simulation of highly charged globular proteins in explicit solvent. However, the lower computational demand of the RF method in the present implementation represents a clear advantage over the P3M method.

Algorithms↗
PubMed2002

Molecular dynamics simulations of the dinuclear zinc-beta-lactamase from Bacteroides fragilis complexed with imipenem.

Herein, we present results from MD simulations of the Michaelis complex formed between the dizinc beta-lactamase from B. fragilis and imipenem. We considered two catalytically important configurations, which differ in the presence or absence of a hydroxide bridge connecting the two zinc ions in the active site. The structural and dynamical effects induced by substrate binding, the specific roles of the conserved residues and the zinc-bound water molecules, the near attack conformers of the Michaelis complex, and so forth, are discussed in detail. The relative stability of the two configurations was estimated from QM linear scaling calculations on the enzyme-substrate complex combined with Poisson-Boltzmann electrostatic calculations and normal mode calculations. Importantly, we find that the two configurations have similar energies, indicating that these two structures could readily be interchanged, thereby facilitating catalysis. The configuration with the hydroxide bound to the two zinc ions is predicted to be the resting form of the enzyme, while the configuration without the bridge is the reactive form that was found to place the hydroxide in position to attack the carbonyl of the beta-lactam ring. Thus, we propose that the enzyme initiates catalysis by converting from the hydroxide bridge form into the configuration that lacks the hydroxide bridge. This interconversion increases the nucleophilicity of the hydroxide ion and exposes it to the beta-lactam carbonyl, which ultimately facilitates nucleophilic attack. The implications of the observed modes of binding, the possible influence of mutating the Lys184 and Asn193 residues on substrate binding, and the reaction mechanism are also discussed in detail.

Algorithms↗
PubMed1996

Molecular dynamics simulation of conformational flexibility of alamethicin fragments in aqueous and membranous environment.

We present here results on molecular dynamics (MD) simulation on two fragments of channel forming antibiotic peptide Alamethicin, containing isoamino butyric acid (Aib). Simulations are carried out in aqueous and membranous environment in a bilayer of 39 molecules of Dimyristoyl phosphatidyl choline (DMPC). The peptides Boc-Pro-Aib-Ala-Aib-OBzl (Alam 1) and Boc-Leu-Aib-Pro-OBzl (Alam 2) were simulated from their crystallographic coordinates. The bilayers were built from two different conformations (A and B) of DMPC reported in crystal data. The P-N dipoles were arranged hexagonally with surface area per lipid molecule 66.5 A degrees 2 and P-P separation across the bilayer 34 A degrees. They were hydrated by 28.6 and 25.5 water molecules per DMPC molecule. Simulations are done using AMBER 4.0 package in constant number volume temperature (NVT) condition for 100 pico seconds (ps) in aqueous environment and 250 ps of equilibrated bilayer. Geometric parameters of lipids as: bilayer thickness, order parameter of the chains, transfraction of chain torsional angles were monitored. We also monitored geometric parameters of the peptides as backbone torsional angles, distances amongst C alpha atoms, angles between C alpha atoms, movement of center of gravity (CG) along and perpendicular to bilayer normal. We find that membrane bilayer is slightly disturbed due to the presence of peptides. In case of alam 2 in water angles phi 1 and phi 3 showed larger variation in water compared to same in the bilayer. The peptide conformation is more stable in DMPC bilayer. However the peptides showed movement along and perpendicular to bilayer normal. This we believe is due to hydrophobic nature of these peptides.

Alamethicin↗
PubMed2005

Simulations of transient membrane behavior in cells subjected to a high-intensity ultrashort electric pulse.

A molecular dynamics (MD) scheme is combined with a distributed circuit model for a self-consistent analysis of the transient membrane response for cells subjected to an ultrashort (nanosecond) high-intensity (approximately 0.01-V/nm spatially averaged field) voltage pulse. The dynamical, stochastic, many-body aspects are treated at the molecular level by resorting to a course-grained representation of the membrane lipid molecules. Coupling the Smoluchowski equation to the distributed electrical model for current flow provides the time-dependent transmembrane fields for the MD simulations. A good match between the simulation results and available experimental data is obtained. Predictions include pore formation times of about 5-6 ns. It is also shown that the pore formation process would tend to begin from the anodic side of an electrically stressed membrane. Furthermore, the present simulations demonstrate that ions could facilitate pore formation. This could be of practical importance and have direct relevance to the recent observations of calcium release from the endoplasmic reticulum in cells subjected to such ultrashort, high-intensity pulses.

Animals↗
PubMed2001

Molecular dynamics simulation of the human U2B" protein complex with U2 snRNA hairpin IV in aqueous solution.

A 2200-ps molecular dynamics (MD) simulation of the U2 snRNA hairpin IV/U2B" complex was performed in aqueous solution using the particle mesh Ewald method to consider long-range electrostatic interactions. To investigate the interaction and recognition process between the RNA and protein, the free energy contributions resulting from individual amino acids of the protein component of the RNA/protein complex were calculated using the recently developed glycine-scanning method. The results revealed that the loop region of the U2 snRNA hairpin IV interacted mainly with three regions of the U2B" protein: 1) beta 1-helix A, 2) beta 2-beta 3, and 3) beta 4-helix C. U2 snRNA hairpin IV bound U2B" in a similar orientation as that previously described for U1 snRNA with the U1A' protein; however, the details of the interaction differed in several aspects. In particular, beta 1-helix A and beta 4-helix C in U2B" were not observed to interact with RNA in the U1A' protein complex. Most of the polar and charged residues in the interacting regions had larger mutant free energies than the nonpolar residues, indicating that electrostatic interactions were important for stabilizing the RNA/protein complex. The interaction was further stabilized by a network of hydrogen bonds and salt bridges formed between RNA and protein that was maintained throughout the MD trajectory. In addition to the direct interactions between RNA and the protein, solvent-mediated interactions also contributed significantly to complex stability. A detailed analysis of the ordered water molecules in the hydration of the RNA/protein complex revealed that bridged water molecules reside at the interface of RNA and protein as long as 2100 ps in the 2200-ps trajectory. At least 20 bridged water molecules, on average, contributed to the instantaneous stability of the RNA/protein complex. The stabilizing interaction energy due to bridging water molecules was obtained from ab initio Hartree-Fock and density functional theory calculations.

Autoantigens↗
PubMed2004

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide: can the alpha3 and beta2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide-MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the alpha1-alpha2 domain, the alpha3 and beta2m domains cannot be neglected.

Algorithms↗
PubMed2005

Ligand-binding regulation of LXR/RXR and LXR/PPAR heterodimerizations: SPR technology-based kinetic analysis correlated with molecular dynamics simulation.

Liver X receptor (LXR) and peroxisome proliferator-activated receptor (PPAR) are two members of nuclear receptors involved in the nutrient metabolisms of dietary fatty acid and cholesterol. They are found to be of cross-talk function in that LXR regulates fatty acid synthesis and PPAR controls fatty acid degradation. LXRs (LXRalpha and LXRbeta) function by forming obligate heterodimers with the retinoid X receptor (RXR), and subsequently binding to specific DNA response elements within the regulatory regions of their target genes. In this work, the kinetic features concerning LXR/RXR and LXR/PPAR interactions have been fully investigated using surface plasmon resonance (SPR) technology. It is found that LXRs could bind to all the three PPAR subtypes, PPARalpha, PPARgamma and PPARdelta with different binding affinities, and such receptor/receptor interactions could be regulated by ligand binding. Moreover, molecular dynamics (MD) simulations were performed on six typical complex models. The results revealed that ligands may increase the interaction energies between the receptor interfaces of the simulated receptor/receptor complexes. The MD results are in agreement with the SPR data. Further analyses on the MD results indicated that the ligand binding might increase the hydrogen bonds between the interfaces of the receptor/receptor complex.

DNA-Binding Proteins↗
PubMed2005

Efficient quantum monte carlo energies for molecular dynamics simulations.

A method is presented to treat electrons within the many-body quantum Monte Carlo (QMC) approach "on-the-fly" throughout a molecular dynamics (MD) simulation. Our approach leverages the large (10-100) ratio of the QMC electron to MD ion motion to couple the stochastic, imaginary-time electronic and real-time ionic trajectories. This continuous evolution of the QMC electrons results in highly accurate total energies for the full dynamical trajectory at a fraction of the cost of conventional, discrete sampling. We show that this can be achieved efficiently for both ground and excited states with only a modest overhead to an ab initio MD method. The accuracy of this dynamical QMC approach is demonstrated for a variety of systems, phases, and properties, including optical gaps of hot silicon quantum dots, dissociation energy of a single water molecule, and heat of vaporization of liquid water.

Journal Article↗
PubMedDate not supplied

A molecular dynamics simulation of the (dG)6 . (dC)6 minihelix including counterions and water.

The results of a 60 ps molecular dynamics (MD) simulation of (dG)6.(dC)6 including 10 Na+ counterions and 292 water molecules are presented. All backbone angles and helix parameters for the hexamer are reported in this paper along with trajectory plots of selected angles. Hydrogen bonding between the bases along the helical axis was observed to fluctuate with time, showing the dynamic nature of the base-pairing interaction. These fluctuations gave rise to unusual hydrogen-bonding patterns. Good intrastrand base stacking and no interstrand base stacking were also observed. The hexamer minihelix retains an essentially B-DNA conformation throughout the entire simulation even though some helix parameters and backbone angles do not have strict B-DNA values. The most striking feature obtained from the simulation was a high propeller twist, which resulted in a narrow minor groove for the minihelix. It is proposed that (dG)n.(dC)n sequences are resistant to DNAase I because of this narrow minor groove in dilute aqueous solution.

Base Sequence↗
PubMed2002

Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation.

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.

Computational Biology↗
PubMed2006

Folding mechanism of beta-hairpins studied by replica exchange molecular simulations.

The folding process of trpzip2 beta-hairpin is studied by the replica exchange molecular dynamics (REMD) and normal MD simulations, aiming to understand the folding mechanism of this unique small, stable, and fast folder, as well as to reveal the general principles in the folding of beta-hairpins. According to our simulations, the TS ensemble is mainly characterized by a largely formed turn and the interaction between the inner pair of hydrophobic core residues. The folding is a zipping up of hydrogen bonds. However, the nascent turn has to be stabilized by the partially formed hydrophobic core to cross the TS. Thus our folding picture is in essence a blend of hydrogen bond-centric and hydrophobic core-centric mechanism. Our simulations provide a direct evidence for a very recent experiment (Du et al., Proc Natl Acad Sci USA 2004;101:15915-15920), which suggests that the turn formation is the rate-limiting step for beta-hairpin folding and the unfolding is mainly determined by the hydrophobic interactions. Besides, the relationship between hydrogen bond stabilities and their relative importance in folding are investigated. It is found that the hydrogen bonds with higher stabilities need not play more important roles in the folding process, and vice versa.

Computer Simulation↗
PubMed2006

A combined molecular dynamics simulation and quantum chemical study on the mechanism for activation of the OxyR transcription factor by hydrogen peroxide.

Molecular dynamics (MD) simulations have been performed on the regulatory domain of the Escherichia coli OxyR transcription factor for the different chemical states along the mechanistic cycle for its activation by hydrogen peroxide. Conformational analysis indicates that His198 and Arg220 catalytic residues can be involved in the biochemical process of activation of OxyR. On the basis of the simulation data, a detailed mechanism for the oxidation process is suggested in which His198, in the presence of an arginine residue, functions as a unique acid-base catalyst in the successive oxidations of Cys199 and Cys208 by hydrogen peroxide. This mechanistic proposal has been tested by density functional theory (DFT-B3LYP) and ab initio (MP2) calculations on model systems. The two oxidations are both identified as nucleophilic substitution reactions of SN2 type with deprotonated cysteines functioning as nucleophiles. Both reactions have a calculated free energy of activation close to 15 kcal mol-1, which is consistent with the available experimental data on the kinetics of the activation process.

Amino Acids↗