Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic 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 1,495 records · Page 83Linked to original sources
PubMed2000

Molecular dynamics simulations of the gramicidin A-dimyristoylphosphatidylcholine system with an ion in the channel pore region.

To investigate the process of ion permeation in an ion channel systematically, we performed molecular dynamics (MD) simulations on a gramicidin A (GA)-phospholipid model system with an ion in the channel pore region. Each of the three types of ions (Ca2+, Na+ Cl-) was placed at five different positions along the channel axis by replacing a water molecule. MD simulations were performed on each system at constant pressure and constant temperature. The MD trajectories showed that the Ca2+ and Na+ ions could stably fluctuate in the pore region, but the Cl- ion was pushed out because of the unfavorable interaction with the channel. This result is consistent with experimental data. It was also found that the conformation of the GA channel underwent a significant change due to the presence of the ion, and the two ends of the GA monomer were more flexible than its middle region. In particular, the dramatic change of local pore radius near the ion indicated this kind of deformation. The strong interaction between the ion and carbonyl oxygen atoms of GA was the major contributor to this change. Furthermore, it was found that the ethanolamine group of the GA molecule was the most flexible group in the GA channel and often observed to block the entrance of GA. These results imply that the deformation of channel structure plays a very important factor in ion permeation, and the ethanolamine group may play a key role in regulating ion entry into the pore. In conclusion, our results indicate that the ion has a dominant influence on the structure of the GA channel and that the flexibility of the ion channel is a crucial factor in the ion permeation process.

Calcium↗
PubMed1995

Crankshaft motions of the polypeptide backbone in molecular dynamics simulations of human type-alpha transforming growth factor.

Order parameters for the backbone N-H and C alpha-H bond vectors have been calculated from a 150 ps molecular dynamics (MD) simulation of human type-alpha transforming growth factor in H2O solvent. Two kinds of 'crankshaft motions' of the polypeptide backbone are observed in this MD trajectory. The first involves small-amplitude rocking of the rigid peptide bond due to correlated changes in the backbone dihedral angles psi i-1 and phi i. These high-frequency 'librational crankshaft' motions are correlated with systematically smaller values of motional order parameters for backbone N-H bond vectors compared to C alpha-H bond vectors. In addition, infrequent 'crankshaft flips' of the peptide bond from one local minimum to another are observed for several amino acid residues. These MD simulations demonstrate that comparisons of N-H and C alpha-H order parameters provide a useful approach for identifying crankshaft librational motions in proteins.

Amino Acid Sequence↗
PubMed2001

Unrestrained 5 ns molecular dynamics simulation of a cisplatin-DNA 1,2-GG adduct provides a rationale for the NMR features and reveals increased conformational flexibility at the platinum binding site.

A 5 ns unrestrained molecular dynamics (MD) simulation of the DNA duplex d(GCCG*G*ATCGC)-d(GCGATCCGGC), bearing a cis-Pt(NH(3))(2)(2+) unit crosslinking the two G* guanine bases, is reported. The MD trajectory was a posteriori correlated with NMR data determined for the same adduct, and it is shown that interproton distances and the characteristic chemical shifts are accounted for by the simulation. The simulation and its confrontation with the NMR data have confirmed the finding derived early from static models that the cytosine complementary to the 5' G*, C17, is mobile with respect to its adjacent bases. However, in contrast to our previous description of this mobility, which included rupture of the Watson-Crick hydrogen bonds and formation of non-Watson-Crick hydrogen bonds, the MD simulation indicated that the G*4-C17 pair moves continuously along a trajectory roughly perpendicular to the local helix axis, with retention of all three Watson-Crick hydrogen bonds. The simulation indicated the reversible formation of a hydrogen bond between the 5' oriented NH(3) ligand of platinum and the C3pG*4 phosphate group, in accord with our former prediction. Furthermore, the simulation has disclosed previously undetected BI <=> BII transitions at the G*5pA6 and A6pT7 steps, connected to formation/rupture of a hydrogen bond between the 3' oriented NH(3) ligand of platinum and the N7 atom of A6. All these conformational equilibria affect the form of the minor groove and increase the conformational flexibility at the platination site, and are thus likely to facilitate recognition by cellular proteins.

Base Sequence↗
PubMed2006

Steered molecular dynamics simulations of Na+ permeation across the gramicidin A channel.

The potential of mean forces (PMF) governing Na+ permeation through gramicidin A (gA) channels with explicit water and membrane was characterized using steered molecular dynamics (SMD) simulations. Constant-force SMD with a steering force parallel to the channel axis revealed at least seven energy wells in each monomer of the channel dimer. Except at the channel dimer interface, each energy well is associated with at least three and at most four backbone carbonyl oxygens and two water oxygens in a pseudo-hexahedral or pseudo-octahedral coordination with the Na+ ion. Repeated constant-velocity SMD by dragging a Na+ ion from each energy well in opposite directions parallel to the channel axis allowed the computation of the PMF across the gA channel, revealing a global minimum corresponding to Na+ binding sites near the entrance of gA at +/-9.3 A from the geometric center of the channel. The effect of volatile anesthetics on the PMF was also analyzed in the presence of halothane molecules. Although the accuracy of the current PMF calculation from SMD simulations is not yet sufficient to quantify the PMF difference with and without anesthetics, the comparison of the overall PMF profiles nevertheless confirms that the anesthetics cause insignificant changes to the structural makeup of the free energy wells along the channel and the overall permeation barrier. On average, the PMF appears less rugged in the outer part of the channel in the presence of anesthetics, consistent with our earlier finding that halothane interaction with anchoring residues makes the gA channel more dynamic. A causal relationship was observed between the reorientation of the coordinating backbone carbonyl oxygen and Na+ transit from one energy well to another, suggesting the possibility that even minute changes in the conformation of pore-lining residues due to dynamic motion could be sufficient to trigger the ion permeation. Because some of the carbonyl oxygens contribute to Na+ coordination in two adjacent energy wells, our SMD results reveal that the atomic picture of ion "hopping" through a gA channel actually involves a Na+ ion being carried in a relay by the coordinating oxygens from one energy well to the next. Steered molecular dynamics complements other computational approaches as an attractive means for the atomistic interpretation of experimental permeation studies.

Gramicidin↗
PubMed2006

Aqueous ionic and complementary zwitterionic soluble surfactants: molecular dynamics simulations and sum frequency generation spectroscopy of the surfaces.

An aqueous ionic surfactant, 1-dodecyl-4-(dimethylamino)pyridinium (DMP) bromide, and the corresponding zwitterion 2-[4-(dimethylamino)pyridinio]dodecanoate (DPN) were explored by means of molecular dynamics (MD) simulations and, for the ionic system, by infrared-visible sum frequency generation (IR-vis SFG). The molecular structure of the interfacial layer was investigated for the ionic and zwitterionic systems as a function of surfactant concentration, both in water and in salt (KF or KBr) solutions, by MD simulations in a slab geometry. The buildup of the surface monolayer and a sublayer was monitored, and density and orientational profiles of the surfactants were evaluated. The difference between the ionic and zwitterionic systems and the effect of the added salt were analyzed at the molecular level. The results of MD simulations were compared to those of nonlinear optical spectroscopy measurements. IR-vis SFG was employed to study the DMP ionic surfactant in water and upon addition of simple salts. The influence of added salts on the different molecular moieties at the interface was quantified in detail experimentally.

Ions↗
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↗
PubMed1995

Molecular dynamics simulation of anesthetic-phospholipid bilayer interactions.

To probe the hypothesis of a lipid-mediated mechanism of general anesthetic action on a molecular level, and to help elucidate the nature of the interactions of bioactive compounds with membranes, the effects of trichloroethylene (TCE), an inhalational general anesthetic, on a dioleoylphosphatidylcholine (DOPC) lipid bilayer have been investigated by molecular dynamics (MD) simulations at 37 degrees C and 1 atm and the results compared with 31P and 2H NMR experimental studies (Ref 1). The model used included a single TCE molecule embedded in a lipid bilayer consisting of 24 DOPC molecules and an 8 A layer of explicit water of solvation in each polar head group region of the bilayer, together with constant-pressure periodic boundary conditions in three dimensions. A comparison of the bilayer properties calculated in the presence and absence of the anesthetic led to the detection of three major perturbations of the bilayer caused by the anesthetic at 1 atm: i) an increase in the ratio of the effective areas of hydrocarbon tails and the head group per lipid, predicting the tendency of lipids near the anesthetic site of action to form a hexagonal phase (HII); ii) a slight increase in the frequency of chain dihedral angles found in the gauche conformation; and iii) a significant increase in the lateral mean-square displacement of lipid molecules, an indication of increased lipid lateral diffusion and membrane fluidity. The pressure antagonism of these effects was also studied by MD simulations at pressures of 200 and 400 atm. The study of the pressure reversibility of these effects at 200 and 400 atm indicated that they were partially prevented at 200 atm and essentially blocked at 400 atm, suggesting their probable relevance to the pressure reversal effect seen with general anesthesia. These results may thus provide insights into the interaction between general anesthetics and similar small organic molecules with membranes.

Computer Graphics↗
PubMed2006

Local order in aqueous NaCl solutions and pure water: X-ray scattering and molecular dynamics simulations study.

The microstructures of pure water and aqueous NaCl solutions over a wide range of salt concentrations (0-4 m) under ambient conditions are characterized by X-ray scattering and molecular dynamics (MD) simulations. MD simulations are performed with the rigid SPC water model as a solvent, while the ions are treated as charged Lennard-Jones particles. Simulated data show that the first peaks in the O...O and O...H pair correlation functions clearly decrease in height with increasing salt concentration. Simultaneously, the location of the second O...O peak, the signature of the so-called tetrahedral structure of water, gradually disappears. Consequently, the degree of hydrogen bonding in liquid water decreases when compared to pure fluid. MD results also show that the hydration number around the cation decreases as the salt concentration increases, which is most likely because some water molecules in the first hydration shell are occasionally substituted by chlorine. In addition, the fraction of contact ion pairs increases and that of solvent-separated ion pairs decreases. Experimental data are analyzed to deduce the structure factors and the pair correlation functions of each system. X-ray results clearly show a perturbation of the association structure of the solvent and highlight the appearance of new interactions between ions and water. A model of intermolecular arrangement via MD results is then proposed to describe the local order in each system, as deduced from X-ray scattering data.

Computer Simulation↗
PubMed2005

Tests of the homogeneous nucleation theory with molecular-dynamics simulations. I. Lennard-Jones molecules.

Two kinds of the homogeneous nucleation theory exist at the present: the classical nucleation theory and the semiphenomenological model. To test them, we performed molecular-dynamics (MD) simulations of nucleation from vapor to liquid with 5000-20,000 Lennard-Jones-type molecules. Simulations were done for various values of supersaturation ratios (from 2 to 10) and temperatures (from 80 to 120 K). We compared the size distribution of clusters in MD simulations with those in the theoretical models because the number density of critical clusters governs the nucleation rate. We found that the semiphenomenological model achieves excellent agreements in size distributions of the clusters with all MD simulations we done. The classical theory underestimates the number density of the clusters in the temperature range of 80-100 K, but overestimates in 100-120 K. The semiphenomenological model also predicts well the nucleation rate in MD simulations, while the classical nucleation theory does not. Our results confirmed the validity of the semiphenomenological model for Lennard-Jones-type molecules.

Journal Article↗
PubMed2006

Structural and dynamical properties of a full-length HIV-1 integrase: molecular dynamics simulations.

The structural and dynamical properties of the complete full-length structure of HIV-1 integrase were investigated using Molecular Dynamics approach. Simulations were carried out for the three systems, core domain only (CORE), full-length structure without (FULL) and with a Mg2+ (FULL+ION) in its active site, aimed to investigate the difference in the molecular properties of the full-length models due to their different construction procedures as well as the effects of the two ends, C- and N-terminal, on those properties in the core domain. The full-length structure was prepared from the two experimental structures of two-domain fragment. The following properties were observed to differ significantly from the previous reports: (i) relative topology formed by an angle between the three domains; (ii) the cavity size defined by the catalytic triad, Asp64, Asp116, and Glu152; (iii) distances and solvation of the Mg2+; and (iv) conformation of the catalytic residues. In addition, the presence of the two terminal domains decreases the mobility of the central core domain significantly.

Binding Sites↗
PubMed2004

2H-NMR study and molecular dynamics simulation of the location, alignment, and mobility of pyrene in POPC bilayers.

The alignment of pyrene in a 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer was investigated using two different approaches, namely solid-state (2)H-NMR spectroscopy and molecular dynamics (MD) simulations. Quadrupolar splittings from (2)H-NMR spectra of deuterated pyrene-d(10) in an oriented lipid bilayer give information about the orientation of C-D bonds with respect to the membrane normal. From MD simulations, geometric information is accessible via trajectories. By defining molecular and bond order parameters, the data from MD trajectories and NMR spectra can be compared straightforwardly. To ensure that the results from both methods are comparable, parameters of the experimental and the simulation setup were chosen to be as similar as possible. From simulations, we saw that pyrene prefers a position inside the lipid membrane near the headgroups and has no tendency to diffuse from one monolayer of the membrane to the other. The results from simulation and NMR show that the normal of the molecular plane is aligned nearly perpendicular to the bilayer normal. The long axis of pyrene lies preferentially parallel to the bilayer normal within a range of +/-30 degrees . The results from the two different methods are remarkably consistent. The good agreement can be explained by the fact that the different kind of motions of a pyrene molecule are already averaged within a few nanoseconds, which is the timescale covered by the MD simulation.

Computer Simulation↗
PubMed2005

The effect of the counterion on water mobility in reverse micelles studied by molecular dynamics simulations.

In this study, mobility and structure of water molecules in Aerosol OT (bis(2-ethylhexyl) sulfosuccinate, AOT) reverse micelles with water content w0 = 5 and Na+, K+, Cs+ counterions have been explored with molecular dynamics (MD) simulations. Using the Faeder/Ladanyi model (J. Phys. Chem. B, 2000, 104, 1033) of the reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function, FS(Q,t), for water hydrogen atoms that could be measured in a quasielastic neutron scattering experiment. Separate intermediate scattering functions FRS(Q,t) and FCMS(Q,t) were determined for rotational and translational motion. We find that the decay of FCMS(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior arises from decreased water mobility for molecules close to the interface and from confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay, which is consistent with relatively rapid restricted rotation and slower rotational relaxation over the full angular range. Rotational relaxation is anisotropic, with the O-H bond short-time rotational mobility considerably higher than that of the molecular dipole. This behavior is related to the decreased density of water-water hydrogen bonds in the vicinity of the interface compared to core or bulk water. We find that the interfacial mobility of water molecules is quite different for the three counterion types, but that the core mobility exhibits weak counterion dependence. Differences in interfacial mobility are strongly correlated with structural features, especially ion-water coordination, and the extent of disruption by the counterions of the water hydrogen bond network.

Journal Article↗
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↗
PubMed2004

Effective force fields for condensed phase systems from ab initio molecular dynamics simulation: a new method for force-matching.

A novel least-squares fitting approach is presented to obtain classical force fields from trajectory and force databases produced by ab initio (e.g., Car-Parrinello) molecular dynamics (MD) simulations. The method was applied to derive effective nonpolarizable three-site force fields for liquid water at ambient conditions from Car-Parrinello MD simulations in the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure includes a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The various parameterizations of the water force field differ by an enforced smooth cut-off applied to the short-ranged interaction term. These were obtained by fitting to the trajectory and force data produced by Car-Parrinello MD simulations of systems of 32 and 64 H(2)O molecules. The new water force fields were developed assuming both flexible or rigid molecular geometry. The simulated structural and self-diffusion properties of liquid water using the fitted force fields are in close agreement with those observed in the underlying Car-Parrinello MD simulations. The resulting empirical models compare to experiment much better than many conventional simple point charge (SPC) models. The fitted potential is also shown to combine well with more sophisticated intramolecular potentials. Importantly, the computational cost of the new models is comparable to that for SPC-like potentials.

Journal Article↗
PubMed2005

Hydrogen bonding in supercritical tert-butanol assessed by vibrational spectroscopies and molecular-dynamics simulations.

We have investigated the state of aggregation in supercritical tert-butanol (T = 523 K,0.05 < rho < 0.4 g cm(-3)) by means of vibrational spectroscopies (infrared and Raman) and molecular-dynamics (MD) simulations. A quantitative band shape analysis of the spectra associated with the OH stretching mode of tert-butanol has been done using activities computed by ab initio calculations on small clusters. This allows us to determine the degree of hydrogen bonding and populations of oligomers. These latter quantities have been derived from MD simulations and very consistent results are found with experiments. These results show that hydrogen bond still exist in supercritical tert-butanol and that the fluid mainly consists of oligomers smaller than tetramers.

Journal Article↗
PubMed2006

Shear-induced parallel-to-perpendicular orientation transition in the amphiphilic lamellar phase: a nonequilibrium molecular-dynamics simulation study.

The present work is devoted to a study of the shear-induced parallel-to-perpendicular orientation transition in the lamellar system by the large-scale nonequilibrium molecular-dynamics (NEMD) simulation. An effective generic model-A2B2 tetramer for amphiphilies is used. The NEMD simulation produces unambiguous evidence that undulation instability along the vorticity direction sets in well above a critical shear rate and grows in magnitude as the shear rate is further increased. At a certain high shear rate, the coherent undulation instability grows so large that defects are nucleated and the global lamellar monodomain breaks into several aligned lamellar domains. Subsequently layers in these domains rotate into the perpendicular orientation with the rotation of chains towards the y direction, merge into a global perpendicular-aligned lamellar monodomain, and organize into a perfect well-aligned perpendicular lamellar phase by the migration and annihilation of edge dislocations and disclinations. The macroscopic observable viscosity as a function of time or shear rate is correlated with the structural response such as the mesoscopic domain morphology and the microscopic chain conformation. The onset of undulation instability concurs with the start-up of shear-thinning behavior. During the orientation transformation at the high shear rate, the complex time-dependent thixotropic behavior is observed. The smaller viscosity in the perpendicular lamellar phase gives an energetic reason for the shear-induced orientation transition.

Journal Article↗
PubMed2003

Improved prediction of HIV-1 protease-inhibitor binding energies by molecular dynamics simulations.

BACKGROUND: The accurate prediction of enzyme-substrate interaction energies is one of the major challenges in computational biology. This study describes the improvement of protein-ligand binding energy prediction by incorporating protein flexibility through the use of molecular dynamics (MD) simulations. RESULTS: Docking experiments were undertaken using the program AutoDock for twenty-five HIV-1 protease-inhibitor complexes determined by x-ray crystallography. Protein-rigid docking without any dynamics produced a low correlation of 0.38 between the experimental and calculated binding energies. Correlations improved significantly for all time scales of MD simulations of the receptor-ligand complex. The highest correlation coefficient of 0.87 between the experimental and calculated energies was obtained after 0.1 picoseconds of dynamics simulation. CONCLUSION: Our results indicate that relaxation of protein complexes by MD simulation is useful and necessary to obtain binding energies that are representative of the experimentally determined values.

Computer Simulation↗
PubMed1999

Molecular dynamics simulations of biomolecules: long-range electrostatic effects.

Current computer simulations of biomolecules typically make use of classical molecular dynamics methods, as a very large number (tens to hundreds of thousands) of atoms are involved over timescales of many nanoseconds. The methodology for treating short-range bonded and van der Waals interactions has matured. However, long-range electrostatic interactions still represent a bottleneck in simulations. In this article, we introduce the basic issues for an accurate representation of the relevant electrostatic interactions. In spite of the huge computational time demanded by most biomolecular systems, it is no longer necessary to resort to uncontrolled approximations such as the use of cutoffs. In particular, we discuss the Ewald summation methods, the fast particle mesh methods, and the fast multipole methods. We also review recent efforts to understand the role of boundary conditions in systems with long-range interactions, and conclude with a short perspective on future trends.

Computer Simulation↗