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 1,387 records · Page 77Linked to original sources

PubMed2005

Molecular dynamics integration and molecular vibrational theory. III. The infrared spectrum of water.

The new symplectic molecular dynamics (MD) integrators presented in the first paper of this series were applied to perform MD simulations of water. The physical properties of a system of flexible TIP3P water molecules computed by the new integrators, such as diffusion coefficients, orientation correlation times, and infrared (IR) spectra, are in good agreement with results obtained by the standard method. The comparison between the new integrators' and the standard method's integration time step sizes indicates that the resulting algorithm allows a 3.0 fs long integration time step as opposed to the standard leap-frog Verlet method, a sixfold simulation speed-up. The accuracy of the method was confirmed, in particular, by computing the IR spectrum of water in which no blueshifting of the stretching normal mode frequencies is observed as occurs with the standard method.

Journal Article↗

PubMed2005

Stability and structure of polyelectrolyte multilayers deposited from salt free solutions.

Molecular-dynamics (MD) simulation results show that polyelectrolyte multilayers deposited from salt free solutions on charged planar surfaces are thermodynamically stable structures that form spontaneously regardless of the method of deposition. The simulation also shows that the polyelectrolyte multilayers are "fuzzy" in nature and molecules in one layer interpenetrate other layers. The influence of chain length, surface charge, and polymer charge is also investigated. Layer thickness was found to be independent of chain length. The ratio of surface to chain charge was found to influence the thickness of the first layer and the amount of polymer absorbed in the first few layers. The thickness of the subsequent layers was found to be independent of the charge ratio.

Journal Article↗

PubMed2005

Brownian dynamics of mixed surfactant micelles.

We investigate micelle formation in a system containing two or more different amphiphiles with different geometries using a stochastic molecular-dynamics (MD) simulation method. For a binary system containing two amphiphiles, we calculate the critical micelle concentration (CMC) and cluster distribution for the mixture at several mole fractions and compare the simulation results with those predicted by analytic theories in the dilute limit and with experiments. We find that the CMC obtained from molecular mean-field theory agrees well with our simulation results. Motivated by the industrial use of mixed surfactant systems, we then extend our studies to a system containing six different chain lengths drawn from a Poisson distribution. We find that unlike a binary mixture of amphiphiles, the different species cancel the effects of each other so that the cluster distribution for the mixture has a shape of a system consisted entirely of amphiphiles of length equal to the mean chain length of the Poisson distribution.

Journal Article↗

PubMed2006

Free-energy calculation of structure-H hydrates.

A molecular-dynamics (MD) simulation of structure-H hydrates was performed under constant pressure and temperature with 6120 TIP4P water molecules, 900 OPLS-UA methane molecules, and 180 large molecular guest substance (LMGS) molecules. The LMGS molecules were represented in the form of a one-site Lennard-Jones (LJ) model using the LJ parameters sigma and epsilon. In order to clarify the thermodynamic stability of structure-H hydrates, we calculated the free-energy difference, changing on the sigma and epsilon only of the LMGS molecules. In this simulation, stable crystals of structure-H hydrates and a minimum value of DeltaG were obtained at sigma approximately 6.2 A and large values of epsilon. All simulations were performed using the special-purpose computer hardware MDGRAPE-2.

Journal Article↗

PubMed2006

Calculating potentials of mean force and diffusion coefficients from nonequilibrium processes without Jarzynski's equality.

In general, the direct application of the Jarzynski equality (JE) to reconstruct potentials of mean force (PMFs) from a small number of nonequilibrium unidirectional steered molecular-dynamics (SMD) paths is hindered by the lack of sampling of extremely rare paths with negative dissipative work. Such trajectories that transiently violate the second law of thermodynamics are crucial for the validity of JE. As a solution to this daunting problem, we propose a simple and efficient method, referred to as the FR method, for calculating simultaneously both the PMF U(z) and the corresponding diffusion coefficient D(z) along a reaction coordinate z for a classical many-particle system by employing a small number of fast SMD pullings in both forward (F) and time reverse (R) directions, without invoking JE. By employing Crooks [Phys. Rev. E 61, 2361 (2000)] transient fluctuation theorem (that is more general than JE) and the stiff-spring approximation, we show that (i) the mean dissipative work W(d) in the F and R pullings is the same, (ii) both U(z) and W(d) can be expressed in terms of the easily calculable mean work of the F and R processes, and (iii) D(z) can be expressed in terms of the slope of W(d). To test its viability, the FR method is applied to determine U(z) and D(z) of single-file water molecules in single-walled carbon nanotubes (SWNTs). The obtained U(z) is found to be in very good agreement with the results from other PMF calculation methods, e.g., umbrella sampling. Finally, U(z) and D(z) are used as input in a stochastic model, based on the Fokker-Planck equation, for describing water transport through SWNTs on a mesoscopic time scale that in general is inaccessible to MD simulations.

Algorithms↗

PubMed2006

Electron hydration dynamics in water clusters: A direct ab initio molecular dynamics approach.

Electron attachment dynamics of excess electron in water cluster (H2O)n (n = 2 and 3) have been investigated by means of full-dimensional direct ab initio molecular dynamics (MD) method at the MP26-311++G(d,p) level. It was found that the hydrogen bond breaking due to the excess electron is an important process in the first stage of electron capture in water trimer. Time scale of electron localization and hydrogen bond breaking were determined by the direct ab initio MD simulation. The initial process of hydration in water cluster is clearly visualized in the present study. In n = 3, an excess electron is first trapped around the cyclic water trimer with a triangular form, where the excess electron is equivalently distributed on the three water molecules at time zero. After 50 fs, the excess electron is concentrated into two water molecules, while the potential energy of the system decreases by -1.5 kcal/mol from the vertical point. After 100 fs, the excess electron is localized in one of the water molecules and the potential energy decreases by -5.3 kcal/mol, but the triangular form still remained. After that, one of the hydrogen bonds in the triangular form is gradually broken by the excess electron, while the structure becomes linear at 100-300 fs after electron capture. The time scale of hydrogen bond breaking due to the excess electron is calculated to be about 300 fs. Finally, a dipole bound state is formed by the linear form of three water molecules. In the case of n = 2, the dipole bound anion is formed directly. The mechanism of electron hydration dynamics was discussed on the basis of theoretical results.

Journal Article↗

PubMed2005

Statistical coevolution analysis and molecular dynamics: identification of amino acid pairs essential for catalysis.

Molecular dynamics (MD) simulations of HhaI DNA methyltransferase and statistical coupling analysis (SCA) data on the DNA cytosine methyltransferase family were combined to identify residues that are coupled by coevolution and motion. The highest ranking correlated pairs from the data matrix product (SCA.MD) are colocalized and form stabilizing interactions; the anticorrelated pairs are separated on average by 30 A and form a clear focal point centered near the active site. We suggest that these distal anti-correlated pairs are involved in mediating active-site compressions that may be important for catalysis. Mutants that disrupt the implicated interactions support the validity of our combined SCA.MD approach.

Binding Sites↗

PubMed2007

Mechanism of methanol oxidation by quinoprotein methanol dehydrogenase.

At neutral pH, oxidation of CH(3)OH --> CH(2)O by an o-quinone requires general-base catalysis and the reaction is endothermic. The active-site -CO(2)(-) groups of Glu-171 and Asp-297 (Glu-171-CO(2)(-) and Asp-297-CO(2)(-)) have been considered as the required general base catalysts in the bacterial o-quinoprotein methanol dehydrogenase (MDH) reaction. Based on quantum mechanics/molecular mechanics (QM/MM) calculations, the free energy for MeOH reduction of o-PQQ when MeOH is hydrogen bonded to Glu-171-CO(2)(-) and the crystal water (Wat1) is hydrogen bonded to Asp-297-CO(2)(-) is DeltaG++ = 11.7 kcal/mol, which is comparable with the experimental value of 8.5 kcal/mol. The calculated DeltaG++ when MeOH is hydrogen bonded to Asp-297-CO(2)(-) is >50 kcal/mol. The Asp-297-CO(2)(-)...Wat1 complex is very stable. Molecular dynamics (MD) simulations on MDH.PQQ.Wat1 complex in TIP3P water for 5 ns does not result in interchange of Asp-297-CO(2)(-) bound Wat1 for a solvent water. Starting with Wat1 removed and MeOH hydrogen bonded to Asp-297-CO(2)(-), we find that MeOH returns to be hydrogen bonded to Glu-171-CO(2)(-) and Asp-297-CO(2)(-) coordinates to Ca(2+) during 3 ns simulation. The Asp-297-CO(2)(-)...Wat1 of reactant complex does play a crucial role in catalysis. By QM/MM calculation DeltaG++ = 1.1 kcal/mol for Asp-297-CO(2)(-) general-base catalysis of Wat1 hydration of the immediate CH(2)==O product --> CH(2)(OH)(2). By this means, the endothermic oxidation-reduction reaction is pulled such that the overall conversion of MeOH to CH(2)(OH)(2) is exothermic.

Alcohol Oxidoreductases↗

PubMed1992

Molecular dynamics studies on nucleoside 2',3'-cyclic phosphates.

2',3'-cyclic nucleotides are intermediates and substrates of Ribonuclease (RNase)-catalysed reactions. The characterization of the equilibrium conformation as well as the flexibility inherent in these molecules helps in understanding the enzymatic action of RNases. The present study explores parameters like phase angle, glycosydic torsion angle and hydrogen bond to find possible interrelationship between them through Molecular Dynamics (MD) simulations on 3'-GMP,3'-UMP, A greater than p, G greater than p, U greater than p, C greater than p, GpA greater than p and UpA greater than p. Interesting results of the effect of cyclisation and other constraints such as hydrogen bond between certain groups on the equilibrium ribose conformation have emerged from this study.

Catalysis↗

PubMed1997

Sequence-independent recombination triple helices: a molecular dynamics study.

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a triplex comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such triplex structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

Base Sequence↗

PubMed1998

The intrinsic curvature of a 51 bp K-DNA fragment of Leishmania tarentolae: a molecular model.

DNA intrinsic structure and curvature is a subject of debate because of the importance of these attributes in processes such as DNA packaging, transcription, and gene regulation. X-ray crystallography of DNA single crystals has provided a wealth of information about the local, short range conformational features of DNA. On the other hand, gel electrophoresis analysis of DNA has not only uncovered the macroscopic curvature of DNA but it also provides most of the available data on DNA intrinsic curvature. However, gel electrophoresis can not identify features of DNA structure at the nucleotide or atomic level. In order to address the problem of DNA intrinsic curvature in an attempt to bridge the gap between X-ray crystallography and gel electrophoresis, we use the computational method of molecular dynamics (MD). In this study, we report the results of 2.0 ns MD simulations on a 51 bp fragment of the K-DNA of Leishmania tarentolae containing several A-tracts. The K-DNA double helix is very stable and remains in an intermediate state between the canonical A and B forms of the duplex. The magnitude of global curvature (75 degrees) agrees well with the experimental estimate (72 degrees) available. Analysis of local (every base triplet) and sublocal (every helix turn) curvature shows that the 51 bp K-DNA fragment has curvature features also present in the Wedge, Junction and Calladine's models of DNA intrinsic curvature. We further characterize the flexibility of individual nucleotides in the molecule and find the sugar flexibility within the A-tracts to be strongly correlated with the pattern of A-tract cleavage by the hydroxyl radical. Differential curvature and flexibility at the 5' and 3'junctions between A-tracts and general-sequence DNA are found to modulate the global curvature of the K-DNA fragment.

Animals↗

PubMed1998

Ion-RNA interactions in the RNA pseudoknot of a ribosomal frameshifting site: molecular modeling studies.

The three-dimensional (3-D) structure of a RNA pseudoknot that causes the efficient ribosomal frameshifting in the gag-pro region of mouse mammary tumor virus (MMTV) has been determined recently by nuclear magnetic resonance (NMR) studies. But since the structure refinement in the studies did not use metal ions and waters, it is not clear how metal ions participate in the stabilization of the pseudoknot, and what kind of ion-RNA interactions dominate in the tertiary contacts for the RNA pseudoknotting. Based on the reported structure data of the pseudoknot VPK of MMTV, we gradually refined the structure by restrained molecular dynamics (MD) using NMR distance restraints. Restrained MD simulation of the RNA pseudoknot was performed with sodium ions and water molecules. Our results are in good agreement with known NMR data and delineate the importance of the metal ion coordination in the stability of the pseudoknot. In the non-coaxially stacking pseudoknot, stem 1 (S1), stem 2 (S2), and the intervening A14 involves unconventional stacking of base pairs coordinated by Na+ and/or bridging water molecules. A6 and G7 of loop L1 make a perfect base stacking in the major groove and are further stabilized by coordinated Na+ ions and water molecules. The first 4-nucleotide (nt) ACUC of loop L2 form a sharp turn and the following 4-nt AAAA cross the minor groove of S1 and are steadied by interactions with the nucleotides of S , bridging water molecules and coordinated Na+ ions. Our studies suggest that the metal ion plays a crucial role in the RNA pseudoknotting of VPK. In the stacking interior of S1 and S2, the Na+ ion is positioned in the major groove and interacts directly with the carbonyl group O6 of G28 and carbonyl group O4 of U13 in the wobble base pair U13:G28. The ion-RNA interactions in MMTV VPK not only stabilize the RNA pseudoknot but also modify the electrostatic properties of the nucleotides at the critical parts of the pseudoknot VPK.

Animals↗

PubMed1999

The effect of structural water molecules on the normal mode spectrum of dTn . dAn x dTn DNA.

In this paper we present a theoretical treatment of triplex B type DNA hydration using normal mode calculation techniques. Discrete solvent is added as spines of hydration in the Watson-Crick and Crick-Hoogsteen grooves as well as water bridges between the Phosphate groups. The effect of binding the discrete structural waters on the normal mode of vibration of the system was studied by introducing a parameter, Xw, that is proportional to the degree of water binding and inversely proportional to the relative humidity (RH) of the system. We examined the variation of the dipole moments of characteristic modes with Xw. The results show that there is a direct relationship between the degree of binding of the water molecules to the atoms in the triple helix, the relative humidity of the system and the conformation and stability of the triple helix. At high RH and Xw = 0:0 the triple helix has mostly B type conformation characteristics, with C'2 -endo sugars. The emergence of normal modes of vibration characteristic to the A type conformation (C'3 - endo sugars) at Xw = 0:4 and 60% RH indicates a conformational shift towards A-type for some of the sugars between Xw = 0.2 (80% RH) and Xw = 0.4 (60% RH). These results are in agreement with the "economy of hydration hypothesis" of Saenger (Saenger et al., 1986) which maintains that the main difference in the hydration of A- and B- forms of DNA is the presence of water bridges between adjacent Phosphate groups in the low-hydration A-form but not in the B- form. Free energy calculations for the triplex DNA with structural waters show that there is a minimum of the free energy at Xw = 0.2 and the free energy increases with Xw and becomes larger than the free energy of the B conformation without structural waters for Xw equal to and larger than 0.4. This result indicates that the B conformation is more stable with bound structural water molecules (for degrees of water binding that are not over 20% higher than the degree of binding between bulk water molecules). The structural water molecules are bound much tighter in the A conformation than in the B conformation. The model predicts that the B to A transition occurs at higher relative humidities in D2O than in H2O. Part of these results (Dadarlat, 1997) have been subsequently confirmed by the experimental work and MD simulations of Ouali (Ouali et al., 1997). The experimental results showed that the N-type sugars corresponding to the A conformation are clearly detected below 75% RH.

Carbohydrate Metabolism↗

PubMed2000

A mezoscopic model of nucleic acids. Part 2. An effective potential energy function for DNA.

In this study we present an effective Potential of Mean Force (PMF) designed for Lagrangian and Quaternion Molecular Dynamics (LQMD) of DNA. The DNA model is built from pseudoatoms as well as rigid and pseudo-elastic bodies described by a limited number of selected Cartesian and internal degrees of freedom. Phosphate groups, deoxyribose rings and nucleic acid bases are represented by pseudoparticles, some of them with internal degrees of freedom. PMF is defined as the sum of effective bonded and long-range potentials. The potentials were fitted to numerical free energy surfaces. Over 50 free energy surfaces, each depending on a conformational variable (pseudobond length, angle or dihedral angle) and the pseudorotation phase of the nearest neighbour deoxribose ring, were computed. The numerical free energy surfaces were obtained from probability distributions derived from a 1.5 ns conventional, microscopic MD simulation of the d(GpC)9 double helical DNA molecule. An umbrella sampling method was used to simulate transitions between the A and B DNA forms, and PMF reproduces these transitions.

Algorithms↗

PubMed2000

Computational studies of essential dynamics of Pseudomonas cepacia lipase.

In order to investigate the interfacial activation of a lipase from Pseudomonas cepacia (PcL), molecular dynamics (MD) simulations and essential dynamics (ED) analysis were performed in different solvent environments: vacuum and explicit water solvents. Starting from the active (open) structure of PcL, the essential dynamics analysis of the simulations revealed large correlated motions, which may be responsible for the activation of the enzyme. Fluctuations in the U1 (active-site lid) and U2 domains are found to be important in the activation of PcL. In contrast, the catalytic triad exhibits very little displacement. These results are consistent with the previous X-ray structural determination. A combined analysis of the trajectories showed some differences for the simulations in different solvent environments. It was found that the region around the helix alpha5 showed larger displacements in the water simulations. It can be concluded that the open structure of PcL becomes unstable in water solvents, leading to the closing of the so-called 'lid' region. The simulations and ED analysis on the closed structure of PgL provided additional information concerning the structural changes involved in the activation of the lipases. It was found that structural changes for PcL and PgL, which are responsible for the essential motions of the protein, showed contrasting behavior in the different solvent environments.

Binding Sites↗

PubMed2000

A 3D model for the human hepatic asialoglycoprotein receptor (ASGP-R).

The human hepatic Asialoglycoprotein Receptor (ASGP-R) consists of two different types of liver specific membrane glycoproteins that bind to terminal galactose and N-acetylgalactosamine residues of serum glycoproteins. The two different polypeptide chains are referred to as two receptor subunits, HH1 and HH2, which are both involved in the activity of the functional receptor. This receptor has served as a model for understanding receptor-mediated endocytosis and carbohydrate mediated recognition phenomena. Here models for the C-terminal extracellular region of both HH1 and HH2 subunit are presented. The standard homology building procedure was modified in order to make it suitable for the modeling problem at hand. The models for the extracellular regions of HH1 and HH2 were initially constructed by exploiting several fragments, belonging to proteins of known 3D structure, and showing high local sequence similarity with respect to the glycoproteins of interest. Putative binding sites were first hypothesized on the basis of the comparison with other complexes of lectins, the crystal structure of which was available in the Protein Data Bank. A model for the complex involving the HH2 subunit and the typical high affinity ligand N-acetylgalactosamine (NacGal) was refined as the first by a suitable combination of MD simulations and Energy Minimization calculations, since it seemed to quickly converge to a plausible structure. An intermediate model for HH1 was then rebuilt on the basis of the refined model for HH2. It was then submitted to a sequence of molecular dynamics simulations with templates which took into account the secondary structure prediction for a final refinement. The structures of small regions of the models, located around the binding sites, were compared with more recent crystallographic data regarding a complex involving the mutant of Mannose Binding Protein QPDWGH (1BCH entry in the Protein Data Bank) and NacGal. This mutant shows high local sequence similarity with HH1 and HH2 at the binding sites. On the basis of the above comparison, different locations of the binding sites were also considered. In addition to other expected interactions, two hydrophobic interactions were observed in the models with Trp residues (positions 243 in HH1 and 181 or 267 in HH2 respectively) and His residues (positions 256 in HHI and 184 in HH2.respectively). The quality of the models was evaluated by the Procheck program and they seemed plausible. This observation together with analogies found between binding sites of the models and IBCH supported the validity of the models. A further validation element arose by comparison between experimental binding data available in the literature about the homologous rat receptor subunits and theoretical interaction energies evaluated, by means of the DOCK 3.5 program, in models for the rat subunits obtained from the corresponding human ones. The new modeling procedure used here appears to be a well-suited method for structural analysis of small regions, located around the ligands, in proteins of unknown 3D structure.

Amino Acid Sequence↗

PubMed2001

Human aldehyde dehydrogenase catalytic activity and structural interactions with coenzyme analogs.

K(m) and V(max) values for 10 coenzyme analogs never previously studied with any aldehyde dehydrogenase and NADP(+) were compared with those for NAD(+) for three human aldehyde dehydrogenases (EC 1.2.1.3); the cytoplasmic E1 (the product of the aldh1 gene), the mitochondrial E2 (the product of the aldh2 gene) and the cytoplasmic E3 (the product of the aldh9 gene) isozymes. Structural information on changes in coenzyme-protein interactions were obtained via molecular dynamics (MD) studies with the E2 isozyme and quantum mechanical (QM) calculations were used to study changes in charge distribution of the pyridine ring and relative free energies of solvation of the purine ring in the analogs. E1 showed the broadest substrate specificity and was the only isozyme subject to substrate inhibition, both of which are suggested to be due to the two coenzyme conformations observed previously in the sheep crystal structure. NADP(+) selectivity is indicated to be influenced by Glu195 in E1 and E2. Substitutions in the purine ring affected K(m) but not V(max), with the changes in K(m) being dominated by the hydrophobicity of the purine ring as indicted by the QM calculations. Substitutions in the pyridine ring sometimes rendered the coenzymes inactive, with no consistent pattern observed for the three coenzymes. Structural analysis of the coenzyme analog-E2 MD simulations revealed different structural perturbations of the surrounding active site, though interactions with Asn169 and Glu399 were preserved in all cases.

Adult↗

PubMed2003

Molecular modeling of a Leishmania antigen eIF-4A: identification of a potential epitope implicated in the adjuvant effect.

In leishmaniasis, the development of an effective vaccine depends on its capacity to elicit Th1 immune responses. Despite many approaches, recent vaccines do not induce sufficient levels of protection and long-term memory. To overcome this problem, vaccines are used with adjuvants that drive immunity towards Th1 and enhance endogenous production of IL-12, a Th1-promoting cytokine. Although exogenous IL-12 may act as an effective adjuvant, it has an elevated cost and toxic effects. Therefore, it is important to develop cheap and safer adjuvants that trigger endogenous IL-12. Of particular interest is LmeIF a unique Leishmania protein that provides significant adjuvant effects by stimulating high IL-12 production. This investigation was designed to identify the structural factors responsible for the adjuvant effects of LmeIF by establishing the 3D models of LmeIF and MueIF (mouse) by homology modeling based on the X-ray structure of their homologs in yeast and comparing their stereo-electronic properties. The molecular electrostatic potential was computed around each model and used to localize the active site and the most different amino acids between LmeIF and MueIF. Sequence alignment of LmeIF with eIF-4 from other species showed three residues (Q186, A189, E214) in the active site which were peculiar to the Leishmania protein. Long MD simulation was carried out on LmeIF fragment 129-261 to compare its folding with the native protein. Despite a high degree of sequence similarity with different species, we have identified in LmeIF a set of residues unique to the protozoan parasite Leishmania which may be potentially responsible for its adjuvant property. Using LmeIF model, a plausible surface region for binding with its receptor was also identified.

Adjuvants, Immunologic↗