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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20Linked to original sources

PubMed2001

Molecular dynamics simulation of carboxy-myoglobin embedded in a trehalose-water matrix.

We report on a molecular dynamics (MD) simulation of carboxy-myoglobin (MbCO) embedded in a water-trehalose system. The mean square fluctuations of protein atoms, calculated at different temperatures in the 100-300 K range, are compared with those from a previous MD simulation on an H2O-solvated MbCO and with experimental data from Mössbauer spectroscopy and incoherent elastic neutron scattering on trehalose-coated MbCO. The results show that, for almost all the atomic classes, the amplitude of the nonharmonic motions stemming from the interconversion among the protein's conformational substates is reduced with respect to the H2O-solvated system, and their onset is shifted toward higher temperature. Moreover, our simulation shows that, at 300 K, the heme performs confined diffusive motions as a whole, leaving the underlying harmonic vibrations unaltered.

Animals↗

PubMed2004

Continuum description of ionic and dielectric shielding for molecular-dynamics simulations of proteins in solution.

We extend our continuum description of solvent dielectrics in molecular-dynamics (MD) simulations, which has provided an efficient and accurate solution of the Poisson equation, to ionic solvents as described by the linearized Poisson-Boltzmann (LPB) equation. We start with the formulation of a general theory for the electrostatics of an arbitrarily shaped molecular system, which consists of partially charged atoms and is embedded in a LPB continuum. This theory represents the reaction field induced by the continuum in terms of charge and dipole densities localized within the molecular system. Because these densities cannot be calculated analytically for systems of arbitrary shape, we introduce an atom-based discretization and a set of carefully designed approximations. This allows us to represent the densities by charges and dipoles located at the atoms. Coupled systems of linear equations determine these multipoles and can be rapidly solved by iteration during a MD simulation. The multipoles yield the reaction field forces and energies. Finally, we scrutinize the quality of our approach by comparisons with an analytical solution restricted to perfectly spherical systems and with results of a finite difference method.

Binding Sites↗

PubMed2006

The midpoint method for parallelization of particle simulations.

The evaluation of interactions between nearby particles constitutes the majority of the computational workload involved in classical molecular dynamics (MD) simulations. In this paper, we introduce a new method for the parallelization of range-limited particle interactions that proves particularly suitable to MD applications. Because it applies not only to pairwise interactions but also to interactions involving three or more particles, the method can be used for evaluation of both nonbonded and bonded forces in a MD simulation. It requires less interprocessor data transfer than traditional spatial decomposition methods at all but the lowest levels of parallelism. It gains an additional practical advantage in certain commonly used interprocessor communication networks by distributing the communication burden more evenly across network links and by decreasing the associated latency. When used to parallelize MD, it further reduces communication requirements by allowing the computations associated with short-range nonbonded interactions, long-range electrostatics, bonded interactions, and particle migration to use much of the same communicated data. We also introduce certain variants of this method that can significantly improve the balance of computational load across processors.

Algorithms↗

PubMed2004

Molecular dynamics simulations of peptides containing an unnatural amino acid: dimerization, folding, and protein binding.

We have performed molecular dynamics (MD) simulations to study the dimerization, folding, and binding to a protein of peptides containing an unnatural amino acid. NMR studies have shown that the substitution of one residue in a tripeptide beta-strand by the unnatural amino acid Hao (5-HO2CCONH-2-MeO-C6H3-CO-NHNH2) modifies the conformational flexibility of the beta-strand and the hydrogen-bonding properties of its two edges: The number of hydrogen-bond donors and acceptors increases at one edge, whereas at the other, they are sterically hindered. In simulations in chloroform, the Hao-containing peptide 9 (i-PrCO-Phe-Hao-Val-NHBu) forms a beta-sheet-like hydrogen-bonded dimer, in good agreement with the available experimental data. Addition of methanol to the solution induces instability of this beta-sheet, as confirmed by the experiments. MD simulations also reproduce the folding of the synthetic peptide 1a (i-PrCO-Hao-Ut-Phe-Ile-Leu-NHMe) into a beta-hairpin-like structure in chloroform. Finally, the Hao-containing peptide, Ac-Ala-Hao-Ala-NHMe, is shown to form a stable complex with the Ras analogue, Rap1 A, in water at room temperature. Together with the available experimental data, these simulation studies indicate that Hao-containing peptides may serve as inhibitors of beta-sheet interactions between proteins.

Amino Acids, Aromatic↗

PubMed2003

Comparing simulation and experiment of a 2D granular Couette shear device.

We present experiments along with molecular-dynamics (MD) simulations of a two-dimensional (2D) granular material in a Couette cell undergoing slow shearing. The grains are disks confined between an inner, rotating wheel and a fixed outer ring. The simulation results are compared to experimental studies and quantitative agreement is found. Tracking the positions and orientations of individual particles allows us to obtain density distributions, velocity and particle rotation rates for the system. The key issue of this paper is to show the extent to which quantitative agreement between an experiment and MD simulations is possible. Besides many differences in model details and the experiment, the qualitative features are nicely reproduced. We discuss the quantitative agreement/disagreement, give possible reasons, and outline further research perspectives.

Journal Article↗

PubMed2006

Adaptive resolution scheme for efficient hybrid atomistic-mesoscale molecular dynamics simulations of dense liquids.

The adaptive resolution scheme (AdResS) for efficient hybrid particle-based atomistic/mesoscale molecular dynamics (MD) simulations recently introduced by us, [J. Chem. Phys. 123, 224106 (2005)] is extended to high density molecular liquids with spherical boundaries between the atomistic and mesoscale regions. The key feature of this approach is that it allows for a dynamical change of the number of molecular degrees of freedom during the course of a MD simulation by an on-the-fly switching between the atomistic and mesoscopic levels of detail. Pressure and density variations occurring at the atomistic/mesoscale boundary in the original version are considerably reduced employing the improved methodology presented here.

Journal Article↗

PubMed2006

Molecular-dynamics simulations of pyronine 6G and rhodamine 6G dimers in aqueous solution.

We have carried out molecular-dynamics (MD) simulations on dimers of the positively charged laser dyes pyronine 6G (P6G) and rhodamine 6G (R6G) in aqueous solution, generating trajectories of 2.5 ns for various computational protocols. We discuss how the choice of atomic partial charges and the length of the trajectories affect the predicted structures of the dimers and compare our results to those of earlier MD-simulations, which were restricted to only 0.7 ns. Our results confirm that monomers of P6G easily undergo relative rotations within the dimer, but we found new conformations of the R6G dimer at longer simulation times. In addition, we analyzed in detail the energy change during the formation of dimers. With suitable corrections, the electrostatic energy from an Ewald treatment agrees with the results from an approach relying on a residue-based cutoff. For P6G, we show that the strong solvent-mediated electrostatic attraction between the monomers is counteracted by an almost equally large solvent-induced entropy contribution to yield a small driving force to dimer formation, in very good agreement with the free-energy change from a thermodynamic-integration procedure. Thus, earlier rationalizations of the dimer formation, based only on energy arguments, yield a qualitatively wrong picture.

Computer Simulation↗

PubMed2003

Molecular dynamics simulations of lignin peroxidase in solution.

The dynamical and structural properties of lignin peroxidase and its Trp171Ala mutant have been investigated in aqueous solution using molecular dynamics (MD) simulations. In both cases, the enzyme retained its overall backbone structure and all its noncovalent interactions in the course of the MD simulations. Very interestingly, the analysis of the MD trajectories showed the presence of large fluctuations in correspondence of the residues forming the heme access channel; these movements enlarge the opening and facilitate the access of substrates to the enzyme active site. Moreover, steered molecular dynamics docking simulations have shown that lignin peroxidase natural substrate (veratryl alcohol) can easily approach the heme edge through the access channel.

Benzyl Alcohols↗

PubMed2005

Structure elucidation and 3D solution conformation of the antibiotic enduracidin determined by NMR spectroscopy and molecular dynamics.

Enduracidin and ramoplanin belong to the large family of cyclodepsipeptide antibiotics, highly effective against Gram-positive bacteria. The primary and 3D solution structure of ramoplanin is already well known, and the primary structure of enduracidin has been determined by a combination of chemical and NMR spectroscopic methods. Both antibiotics share a similar peptide core of 17 amino acids and differ mainly in the length of the acyl chain and the presence of two D-mannose moieties in ramoplanin. Based on the high sequence homology with ramoplanin, the structure in solution of enduracidin is modeled as a cyclic peptide. The tertiary structure thus obtained was refined through molecular dynamics (MD) simulation, in which the interatomic NOE-derived distance restraints were imposed. MD simulations yielded a family of representative 3D structures (RMSD = 0.89), which highlighted a backbone geometry similar to that of ramoplanin in its beta-hairpin arrangement. In contrast, enduracidin displays a different arrangement of the side-chain and of the residues forming the hydrophobic core.

Amides↗

PubMed2001

Involvement of water in carbohydrate-protein binding.

The interactions of trimannosides 1 and 2 with Con A were studied to reveal the effects of displacement of well-ordered water molecules on the thermodynamic parameters of protein-ligand complexation. Trisaccharide 2 is a derivative of 1, in which the hydroxyl at C-2 of the central mannose unit is replaced by a hydroxyethyl moiety. Upon binding, this moiety displaces a conserved water molecule present in the Con A binding site. Structural studies by NMR spectroscopy and MD simulations showed that the two compounds have very similar solution conformational properties. MD simulations of the complexes of Con A with 1 and 2 demonstrated that the hydroxyethyl side chain of 2 can establish the same hydrogen bonds in a low energy conformation with the protein binding site as those mediated by the water molecule in the complex of 1 with Con A. Isothermal titration microcalorimetry (ITC) measurements showed that 2 has a more favorable entropy of binding compared to 1. This term, which was expected, arises from the return of the highly ordered water molecule to bulk solution. The favorable entropy term was, however, offset by a relatively large unfavorable enthalpy term. This observation was rationalized by comparing the extent of hydrogen bond and solvation changes during binding. It is proposed that an indirect interaction through a water molecule will provide a larger number of hydrogen bonds in the complex that have higher occupancies than in bulk solution, thereby stabilizing the complex.

Carbohydrate Conformation↗

PubMed2006

Computational assessment of the entropy of solvation of small-sized hydrophobic entities.

A high level polarizable force field is used to study the temperature dependence of hydrophobic hydration of small-sized molecules from computer simulations. Molecular dynamics (MD) simulations of liquid water at various temperatures form the basis of free energy perturbation calculations that consider the onset and growth of a repulsive sphere. This repulsive sphere acts as a model construct for the hydrophobic species. In the present study, an extension is pursued for seven independent target temperatures, ranging from close to the freezing point almost up to the boiling point of liquid water under standard conditions. Care is taken to maintain proper physico-chemical model descriptions by cross-checking with experimental water densities at the selected target temperatures. The polarizable force field description of molecular water turns out to be suitable throughout the entire temperature domain considered. Derivatives of the computed free energies of hydrophobic hydration with respect to the temperature give access to the changes in entropy. In practice the entropy differential is determined from the negative of the slope of tangential lines formed at a certain target temperature in the free energy profile. The obtained changes in entropy are negative for small-sized cavities, and hence reconfirm the basic ideas of the Lum-Chandler-Weeks theory on hydrophobic hydration of small-sized solutes.

Computer Simulation↗

PubMed2004

Characterization of the microscopic surface structure of the octadecylsilica stationary phase using a molecular-dynamics simulation.

The influence of the mobile-phase composition and temperature on the surface structure of the octadecylsilica (ODS) stationary-phase was investigated by applying a molecular-dynamics (MD) simulation. The molecular model to which the MD simulation was applied consisted of three parts: an amorphous silica base, dimethyloctadecylsilyl ligands and mobile-phase solvents. More detailed information on the effect of the mobile-phase composition was obtained by constructing larger molecular models than those used in our previous study. The thickness of the hydrocarbon layer of the stationary phase could be estimated based on the distance between the carbon atom located at the terminal end of the ODS ligand and the silica gel surface. The structural information obtained by the calculation showed good consistency with the experimentally observed values. The gauche fraction in the ODS ligand conformation could be also estimated to obtain a more detailed ligand conformation for each molecular model. It was found that as the temperature increased, the ligand conformation collapsed more. This trend was the same as the experimentally observed trends obtained by NMR, FT-IR, and Raman spectroscopic techniques.

Computer Simulation↗

PubMedDate not supplied

Improving the performance of molecular dynamics simulations on parallel clusters.

In this article a procedure is derived to obtain a performance gain for molecular dynamics (MD) simulations on existing parallel clusters. Parallel clusters use a wide array of interconnection technologies to connect multiple processors together, often at different speeds, such as multiple processor computers and networking. It is demonstrated how to configure existing programs for MD simulations to efficiently handle collective communication on parallel clusters with processor interconnections of different speeds.

Journal Article↗

PubMed1997

Alamethicin channels - modelling via restrained molecular dynamics simulations.

Alamethicin channels have been modelled as approximately parallel bundles of transbilayer helices containing between N = 4 and 8 helices per bundle. Initial models were generated by in vacuo restrained molecular dynamics (MD) simulations, and were refined by 60 ps MD simulations with water molecules present within and at the mouths of the central pore. The helix bundles were stabilized by networks of H-bonds between intra-pore water molecules and Gln-7 side-chains. Channel conductances were predicted on the basis of pore radius profiles, and suggested that the N = 4 bundle formed an occluded pore, whereas pores with N > or = 5 helices per bundle were open. Continuum electrostatics calculations suggested that the N = 6 pore is cation-selective, whereas pores with N > or = 7 helices per bundle were predicted to be somewhat less ion-selective.

Alamethicin↗

PubMed2004

A molecular dynamics simulation study of polyamine- and sodium-DNA. Interplay between polyamine binding and DNA structure.

Four different molecular dynamics (MD) simulations have been performed for infinitely long ordered DNA molecules with different counterions, namely the two natural polyamines spermidine(3+) (Spd3+) and putrescine(2+) (Put2+), the synthetic polyamine diaminopropane(2+) (DAP2+), and the simple monovalent cation Na+. All systems comprised a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charge. The simulation setup mimics the DNA state in oriented DNA fibers, previously studied using NMR and other experimental methods. In this paper the interplay between polyamine binding and local DNA structure is analyzed by investigating how and if the minor groove width of DNA depends on the presence and dynamics of the counterions. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural [spermine(4+), Spd3+, Put2+] and the synthetic (DAP2+) polyamines.

Base Sequence↗

PubMed2005

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations.

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.

Computer Simulation↗

PubMed2006

Generalized correlation for biomolecular dynamics.

Correlated motions in biomolecules are often essential for their function, e.g., allosteric signal transduction or mechanical/thermodynamic energy transport. Because correlated motions in biomolecules remain difficult to access experimentally, molecular dynamics (MD) simulations are particular useful for their analysis. The established method to quantify correlations from MD simulations via calculation of the covariance matrix, however, is restricted to linear correlations and therefore misses part of the correlations in the atomic fluctuations. Herein, we propose a general statistical mechanics approach to detect and quantify any correlated motion from MD trajectories. This generalized correlation measure is contrasted with correlations obtained from covariance matrices for the B1 domain of protein G and T4 lysozyme. The new method successfully quantifies correlations and provides a valuable global overview over the functionally relevant collective motions of lysozyme. In particular, correlated motions of helix 1 together with the two main lobes of lysozyme are detected, which are not seen by the conventional covariance matrix. Overall, the established method misses more than 50% of the correlation. This failure is attributed to both, an interfering and unnecessary dependence on mutual orientations of the atomic fluctuations and, to a lesser extent, attributed to nonlinear correlations. Our generalized correlation measure overcomes these problems and, moreover, allows for an improved understanding of the conformational dynamics by separating linear and nonlinear contributions of the correlation.

Bacteriophage T4↗

PubMed2006

Molecular dynamics model of unliganded HIV-1 reverse transcriptase.

HIV-1 RT is one of the most important antiviral targets in the treatment of acquired immunodeficiency syndrome (AIDS). Several crystallographic structures are available for this enzyme, mostly with bound inhibitors. Despite their importance for structure based drug design towards new anti-HIV retrovirals, the X-ray structures of the unliganded enzyme could only be obtained incomplete, with a low resolution and until recently even the conformation of the p66 thumb was controversial. In this work we have aligned different X-ray RT structures, and built up a computational model of RT using homology modeling, which was afterwards refined and validated through MD simulations with explicit solvent. The model enzyme was structurally stable through the whole MD simulation, showing a RMSD of 2 Angstrom from the starting geometry. The Ramanchandram plot has improved along the simulation. Both intra-domain and interdomain movements were observed. The thumb kept its closed conformation through the whole simulation. A contact map, hydration sites study and a detailed analysis of the solvation of the nucleotide binding site are also presented.

Binding Sites↗