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 145 records · Page 8Linked to original sources
PubMed2002

Measuring kinetic coefficients by molecular dynamics simulation of zone melting.

Molecular dynamics simulations are performed to measure the kinetic coefficient at the solid-liquid interface in pure gold. Results are obtained for the (111), (100), and (110) orientations. Both Au(100) and Au(110) are in reasonable agreement with the law proposed for collision-limited growth. For Au(111), stacking fault domains form, as first reported by Burke, Broughton, and Gilmer [J. Chem. Phys. 89, 1030 (1988)]. The consequence on the kinetics of this interface is dramatic: the measured kinetic coefficient is three times smaller than that predicted by collision-limited growth. Finally, crystallization and melting are found to be always asymmetrical and here again the effect is much more pronounced for the (111) orientation.

Journal Article↗
PubMed2002

Complex formation in systems of oppositely charged polyelectrolytes: a molecular dynamics simulation study.

Results of molecular dynamics simulations for systems with two flexible, oppositely charged polymer chains are presented. The lengths N and interaction strength lambda of the chains are varied. We find that the chains remain separated for small values of lambda. For large interaction strengths, i.e., large Bjerrum lengths, we find glasslike structures and order on the length scale of a few monomer diameters. Between these two limits of the interaction strengths, the chains of various lengths collapse into compact complexes that exhibit self-similar structures. The scaling behavior of the radius of gyration is discussed as a function of chain length and interaction strength. In addition, the local structure of the collapsed systems is analyzed and the dependence of the density of the aggregate on the interaction strength is discussed.

Journal Article↗
PubMed2003

Nonequilibrium molecular dynamics simulations of a bubble.

Molecular dynamics simulations are performed to investigate the microscopic dynamics of a bubble when liquids are locally heated. We successfully observe that the heated atoms scatter the neighboring nonheated atoms and make a bubble, and then the bubble is cooled and compressed by the surrounding liquids. The bubble dynamics in this process agrees with the results of the Rayleigh-Plesset equation which describes the dynamics of a bubble in terms of macroscopic hydrodynamics. In this way, we clarify that the hydrodynamic description is reliable even for a microscopic bubble.

Journal Article↗
PubMed2007

Estimation of folding probabilities and phi values from molecular dynamics simulations of reversible Peptide folding.

Molecular dynamics simulations with an implicit model of the solvent have allowed to investigate the reversible folding of structured peptides. For a 20-residue antiparallel beta-sheet peptide, the simulation results have revealed multiple folding pathways. Moreover, the conformational heterogeneity of the denatured state has been shown to originate from high enthalpy, high entropy basins with fluctuating non-native secondary structure, as well as low enthalpy, low entropy traps. An efficient and simple approach to estimate folding probabilities from molecular dynamics simulations has allowed to isolate conformations in the transition state ensemble and to evaluate Phi values, i.e., the effects of mutations on the folding kinetics and thermodynamic stability. These molecular dynamics studies have provided evidence that, if interpreted by neglecting the non-native interactions, Phi values overestimate the amount of native-like structure in the transition state.

Computer Simulation↗
PubMed2003

Active site modeling in copper azurin molecular dynamics simulations.

Active site modeling in molecular dynamics simulations is investigated for the reduced state of copper azurin. Five simulation runs (5 ns each) were performed at room temperature to study the consequences of a mixed electrostatic/constrained modeling for the coordination between the metal and the polypeptide chain, using for the ligand residues a set of charges that is modified with respect to the apo form of the protein by the presence of the copper ion. The results show that the different charge values do not lead to relevant effects on the geometry of the active site of the protein, as long as bond distance constraints are used for all the five ligand atoms. The distance constraint on the O atom of Gly45 can be removed without altering the active site geometry. The coordination between Cu and the other axial ligand Met121 is outlined as being flexible. Differences are found between the bonds of the copper ion with the two apparently equivalent N(delta1) atoms of His46 and His117. The overall findings are discussed in connection with the issue of determining a model for the active site of azurin suitable to be used in molecular dynamics simulations under unfolding conditions.

Azurin↗
PubMed2004

Studies of translational diffusion in the smectic A phase of a Gay-Berne mesogen using molecular dynamics computer simulation.

Molecular dynamics computer simulations are used to determine the self-diffusion coefficients for a Gay-Berne model mesogen GB (4.4,20,1,1) in the isotropic, nematic and smectic A phases along two isobars. The values of the parallel and perpendicular diffusion coefficients, D(parallel) and D(perpendicular), are calculated and compared in the different phases. For the phase sequence isotropic-smectic A, D(perpendicular)*> or =D(parallel)* over the whole smectic A range with the ratio D(parallel)*/D(perpendicular)* decreasing with decreasing temperature. At a higher pressure, a nematic phase is observed between these two phases and we find that D(parallel)*>D(perpendicular)* throughout the nematic region and the inequality D(parallel)*>D(perpendicular)* remains on entering the smectic A phase. However, the ratio D(parallel)*/D(perpendicular)* decreases with decreasing temperature within the smectic A range and eventually this ratio inverts such that D(perpendicular)*>D(parallel)* at low temperatures. The temperature dependence of the parallel diffusion coefficient in the smectic A phase for this model mesogen is compared to that predicted by a theoretical model for diffusion subject to a cosine potential.

Journal Article↗
PubMed2004

Dynamical behavior of one-dimensional water molecule chains in zeolites: nanosecond time-scale molecular dynamics simulations of bikitaite.

Nanosecond scale molecular dynamics simulations of the behavior of the one-dimensional water molecule chains adsorbed in the parallel nanochannels of bikitaite, a rare lithium containing zeolite, were performed at different temperatures and for the fully and partially hydrated material. New empirical potential functions have been developed for representing lithium-water interactions. The structure and the vibrational spectrum of bikitaite were in agreement both with experimental data and Car-Parrinello molecular dynamics results. Classical molecular dynamics simulations were extended to the nanosecond time scale in order to study the flip motion of water molecules around the hydrogen bonds connecting adjacent molecules in the chains, which has been observed by NMR experiments, and the dehydration mechanism at high temperature. Computed relaxation times of the flip motion follow the Arrhenius behavior found experimentally, but the activation energy of the simulated system is slightly underestimated. Based on the results of the simulations, it may be suggested that the dehydration proceeds by a defect-driven stepwise diffusion. The diffusive mechanism appears as a single-file motion: the molecules never pass one another, even at temperatures as high as about 1000 K, nor can they switch between different channels. However, the mean square displacement (MSD) of the molecules, computed with respect to the center of mass of the simulated system, shows an irregular trend from which the single-file diffusion cannot be clearly evidenced. If the MSDs are evaluated with respect to the center of mass of the molecules hosted in each channel, the expected dependence on the square root of time finally appears.

Journal Article↗
PubMed2005

Identifying and engineering ion pairs in adenylate kinases. Insights from molecular dynamics simulations of thermophilic and mesophilic homologues.

Molecular dynamics simulations were performed to study thermal stabilization of proteins via electrostatic interactions of ion pairs. Dynamic motions of four ion pairs previously proposed to be important in thermal stability of adenylate kinase from the thermophile Bacillus stearothermophilus were monitored during the simulation. One of the four ion pairs identified in the crystal structure, Lys180-Asp114, was not maintained in close contact suggesting that the ion pair does not contribute to thermal stability. Among the other three ion pairs, the ion pair Arg116-Glu198 was proposed to be the most important for stability. To predict behaviors of the ion pairs when engineered into a mesophilic homologue to increase stability, in silico mutants of adenylate kinase from the mesophile Bacillus subtilis were generated, and their molecular dynamics simulations were carried out. The ion pairs in the mutant simulations displayed similar behaviors to those in the simulation of the thermophilic protein. To validate the results of the simulations experimentally, the same mutants were produced in vitro and their thermal stabilities were measured using differential scanning calorimetry. In agreement with the simulations, the Lys180-Asp114 did not result in any increase in stability by itself or additive effect with other ion pairs, whereas a mutant with the Arg116-Glu198 exhibited the highest stability among the mutants having one of the four ion pairs. These results provide specific knowledge about stability in adenylate kinases and more generally suggest that molecular dynamics simulations can provide valuable information for identifying and engineering ion pairs.

Adenylate Kinase↗
PubMed2005

Hydration of Li+ -ion in atom-bond electronegativity equalization method-7P water: a molecular dynamics simulation study.

We have carried out molecular dynamics simulations of a Li(+) ion in water over a wide range of temperature (from 248 to 368 K). The simulations make use of the atom-bond electronegativity equalization method-7P water model, a seven-site flexible model with fluctuating charges, which has accurately reproduced many bulk water properties. The recently constructed Li(+)-water interaction potential through fitting to the experimental and ab initio gas-phase binding energies and to the measured structures for Li(+)-water clusters is adopted in the simulations. ABEEM was proposed and developed in terms of partitioning the electron density into atom and bond regions and using the electronegativity equalization method (EEM) and the density functional theory (DFT). Based on a combination of the atom-bond electronegativity equalization method and molecular mechanics (ABEEM/MM), a new set of water-water and Li(+)-water potentials, successfully applied to ionic clusters Li(+)(H(2)O)(n)(n=1-6,8), are further investigated in an aqueous solution of Li(+) in the present paper. Two points must be emphasized in the simulations: first, the model allows for the charges on the interacting sites fluctuating as a function of time; second, the ABEEM-7P model has applied the parameter k(lp,H)(R(lp,H)) to explicitly describe the short-range interaction of hydrogen bond in the hydrogen bond interaction region, and has a new description for the hydrogen bond. The static, dynamic, and thermodynamic properties have been studied in detail. In addition, at different temperatures, the structural properties such as radial distribution functions, and the dynamical properties such as diffusion coefficients and residence times of the water molecules in the first hydration shell of Li(+), are also simulated well. These simulation results show that the ABEEM/MM-based water-water and Li(+)-water potentials appear to be robust giving the overall characteristic hydration properties in excellent agreement with experiments and other molecular dynamics simulations on similar system.

Journal Article↗
PubMed2004

Potential energy constrained molecular dynamics simulations.

A method for carrying out molecular dynamics simulations in which the potential energy U of the molecular system is constrained at its initial value is developed and thoroughly tested. The constraint is not introduced within the framework of the Lagrange multipliers technique, rather it is fulfilled in a natural way by carrying out the simulations in terms of suitable sets of delocalized coordinates. Such coordinates are defined by an appropriate tuning of the Baker, Kessi, and Delley internal delocalized nonredundant coordinates technique [J. Chem. Phys. 105, 192 (1996)]. The proposed method requires multiple evaluations of energy and gradients in each step of the molecular dynamics simulation, so that constant U simulations suffer some overhead compared to ordinary simulations. But the particular formulation of the delocalized coordinates and of the equations of motion greatly simplifies all the various steps required by the Baker's technique, thus allowing for the efficient implementation of the method itself. The technique is reliable and allows for very high accuracy in the potential energy conservation during the whole simulation. Moreover, it proved to be free of drift troubles which can occur when standard constraint methods are straightforwardly implemented without the application of appropriate correcting techniques.

Journal Article↗
PubMed2004

Refinement of homology-based protein structures by molecular dynamics simulation techniques.

The use of classical molecular dynamics simulations, performed in explicit water, for the refinement of structural models of proteins generated ab initio or based on homology has been investigated. The study involved a test set of 15 proteins that were previously used by Baker and coworkers to assess the efficiency of the ROSETTA method for ab initio protein structure prediction. For each protein, four models generated using the ROSETTA procedure were simulated for periods of between 5 and 400 nsec in explicit solvent, under identical conditions. In addition, the experimentally determined structure and the experimentally derived structure in which the side chains of all residues had been deleted and then regenerated using the WHATIF program were simulated and used as controls. A significant improvement in the deviation of the model structures from the experimentally determined structures was observed in several cases. In addition, it was found that in certain cases in which the experimental structure deviated rapidly from the initial structure in the simulations, indicating internal strain, the structures were more stable after regenerating the side-chain positions. Overall, the results indicate that molecular dynamics simulations on a tens to hundreds of nanoseconds time scale are useful for the refinement of homology or ab initio models of small to medium-size proteins.

Algorithms↗
PubMed1994

Molecular dynamics simulations of phospholipid bilayers.

Molecular dynamics (MD) simulations at 37 degrees C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a(t)/ah), correctly predict the aggregation preferences of the three lipids observed experimentally at 37 degrees C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.

Algorithms↗
PubMed2000

Molecular dynamics simulation of nucleic acids.

We review molecular dynamics simulations of nucleic acids, including those completed from 1995 to 2000, with a focus on the applications and results rather than the methods. After the introduction, which discusses recent advances in the simulation of nucleic acids in solution, we describe force fields for nucleic acids and then provide a detailed summary of the published literature. We emphasize simulations of small nucleic acids ( approximately 6 to 24 mer) in explicit solvent with counterions, using reliable force fields and modern simulation protocols that properly represent the long-range electrostatic interactions. We also provide some limited discussion of simulation in the absence of explicit solvent. Absent from this discussion are results from simulations of protein-nucleic acid complexes and modified DNA analogs. Highlights from the molecular dynamics simulation are the spontaneous observation of A B transitions in duplex DNA in response to the environment, specific ion binding and hydration, and reliable representation of protein-nucleic acid interactions. We close by examining major issues and the future promise for these methods.

Base Sequence↗
PubMed2004

Molecular alignment in a liquid induced by a nonresonant laser field: Molecular dynamics simulation.

We carried out molecular dynamics (MD) simulations for a dilute aqueous solution of pyrimidine in order to investigate the mechanisms of field-induced molecular alignment in a liquid phase. An anisotopically polarizable molecule can be aligned in a liquid phase by the interaction with a nonresonant intense laser field. We derived the effective forces induced by a nonresonant field on the basis of the concept of the average of the total potential over one optical cycle. The results of MD simulations show that a pyrimidine molecule is aligned in an aqueous solution by a linearly polarized field of light intensity I approximately 10(13) W/cm2 and wavelength lambda = 800 nm. The temporal behavior of field-induced alignment is adequately reproduced by the solution of the Fokker-Planck equation for a model system in which environmental fluctuations are represented by Gaussian white noise. From this analysis, we have revealed that the time required for alignment in a liquid phase is in the order of the reciprocals of rotational diffusion coefficients of a solute molecule. The degree of alignment is determined by the anisotropy of the polarizability of a molecule, light intensity, and temperature. We also discuss differences between the mechanisms of optical alignment in a gas phase and a liquid phase.

Journal Article↗
PubMed1995

Structure, thermodynamics and cooperativity of the glucocorticoid receptor DNA-binding domain in complex with different response elements. Molecular dynamics simulation and free energy perturbation studies.

Molecular dynamics simulations and free energy perturbation calculations have been performed on the glucocorticoid receptor DNA-binding domain (GRDBD) in complex with three different half-site response elements (RE). These are the glucocorticoid receptor RE (GREH), the estrogen receptor RE (EREH) and an "intermediate" RE (GRE2H), where the base-pair 6A.T in GREH has been replaced by 6G.C. The molecular dynamics/free energy perturbation calculations reveal the same order of stability of the three complexes, i.e. GREH-GRDBD > GRE2H-GRDBD > EREH-GRDBD, as found experimentally, with reasonable quantitative agreement of the relative stabilities with the experiments. The more favourable entropy in associating GRDBD to GREH than to GRE2H could be explained by the additional non-polar surface that is removed from the solvent in the former complex formation. Part of this additional surface originates from the 6T methyl group in GREH. Sequence-dependent hydration patterns of the DNA major groove, which is interrupted by the GRDBD binding, could also contribute to the relatively more favourable entropy in forming GREH-GRDBD than GRE2H-GRDBD. The tenfold reduced cooperativity in the dimeric GRDBD binding to a partially palindromic ERE (palERE) compared to palGRE/GRE2 could be attributed to conformational changes in the residues of the "D-box" (Cys476 to Cys482) observed in the molecular dynamics simulations of EREH-GRDBD, but not in the other two complexes. Moreover, the conformations of the base-pairs at positions 5 and 6 are distorted (compared to free B-DNA) in GREH and GRE2H, due to the binding of GRDBD. Since the corresponding conformational distortions were not found in EREH-GRDBD, it is plausible that the distortions mediate cooperative binding through a propagation of the distortions to the other half site, thus facilitating binding of the second GRDBD.

Amino Acid Sequence↗
PubMed2003

Plausible interaction of an alpha-fetoprotein cyclopeptide with the G-protein-coupled receptor model GPR30: docking study by molecular dynamics simulated annealing.

In this manuscript, the procedure of molecular dynamics simulated annealing is applied to locate a probable receptor and binding site of a cyclicpeptide that inhibits estrogen-stimulated proliferation of breast cancer. The hydrophilic cyclopeptide EMTOVNOGQ (O = 4-hydroxyproline), derived from alpha-fetoprotein, is an inhibitor of estrogen-stimulated proliferation of human breast cancer. This peptide has been shown to act through a mechanism different from that of estrogen; however, its receptor is unknown. We report computer experiments that suggest that this peptide may execute its actions by interacting with GPR30, a G-protein-coupled receptor. The subject of this work is the simulation, by molecular dynamics simulated annealing, of the interaction of cyclopeptide EMTOVNOGQ with receptor GPR30 protein. A conformational analysis of the cyclopeptide was undertaken and the final structure was docked on several sites of the GPR30 3D model. Our results show that the cyclopeptide interacts on the pocket located between TM6 and TM7 transmembrane helices of the G-protein, triggering a slight conformational change in the secondary structure of the receptor in the complex. Based on differences in accessible surface areas between GPR30 and its ligand, the residues in the interaction zone were identified. The cyclopeptide is stabilized in the active site by forming a network of hydrogen bonds between Glu, Thr, (1)Pro(OH) and GLn residues of the ligand and Arg-259, Cys-271, Asn-316, Asn-320 and Tyr-324 of the G-protein. Moreover, the study of the electrostatic surface potential on the GPR30 receptor shows that the active site is more positively charged than the other sites. Our modeling indicates a plausible interaction of the cyclopeptide with the seven transmembrane GPR30 protein. This may have profound implications for the treatment of breast cancer.

Amino Acid Sequence↗
PubMed2002

Contribution of energy values to the analysis of global searching molecular dynamics simulations of transmembrane helical bundles.

Molecular interactions between transmembrane alpha-helices can be explored using global searching molecular dynamics simulations (GSMDS), a method that produces a group of probable low energy structures. We have shown previously that the correct model in various homooligomers is always located at the bottom of one of various possible energy basins. Unfortunately, the correct model is not necessarily the one with the lowest energy according to the computational protocol, which has resulted in overlooking of this parameter in favor of experimental data. In an attempt to use energetic considerations in the aforementioned analysis, we used global searching molecular dynamics simulations on three homooligomers of different sizes, the structures of which are known. As expected, our results show that even when the conformational space searched includes the correct structure, taking together simulations using both left and right handedness, the correct model does not necessarily have the lowest energy. However, for the models derived from the simulation that uses the correct handedness, the lowest energy model is always at, or very close to, the correct orientation. We hypothesize that this should also be true when simulations are performed using homologous sequences, and consequently lowest energy models with the right handedness should produce a cluster around a certain orientation. In contrast, using the wrong handedness the lowest energy structures for each sequence should appear at many different orientations. The rationale behind this is that, although more than one energy basin may exist, basins that do not contain the correct model will shift or disappear because they will be destabilized by at least one conservative (i.e. silent) mutation, whereas the basin containing the correct model will remain. This not only allows one to point to the possible handedness of the bundle, but can be used to overcome ambiguities arising from the use of homologous sequences in the analysis of global searching molecular dynamics simulations. In addition, because clustering of lowest energy models arising from homologous sequences only happens when the estimation of the helix tilt is correct, it may provide a validation for the helix tilt estimate.

Amino Acid Sequence↗
PubMed2001

Direct comparison of experimental and calculated folding free energies for hydrophobic deletion mutants of chymotrypsin inhibitor 2: free energy perturbation calculations using transition and denatured states from molecular dynamics simulations of unfolding.

Previous molecular dynamics (MD) simulations of thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided transition-state models in good agreement with experiment. Unfortunately, however, the comparisons have been necessarily indirect. The simulations have provided detailed structural information but not energetics, while from experiment, structure is inferred from a ratio of free energy changes upon mutation (Phi values). Here, direct comparison with experimental free energies is obtained by performing free energy perturbation calculations of hydrophobic deletion mutants of CI2 using transition- and denatured-state structures from various denaturation MD simulations. The agreement between the calculated and experimental DeltaDeltaG and Phi values is quite good (R = 0.8-0.9). In addition, given the availability of realistic atomic models for the denatured protein, the common approach of using small peptides to represent the denatured state in stability calculations can now be evaluated. To this end, two different extended tripeptide models were used: one using the sequence from the protein with the residue to be mutated in the center and the other with this residue surrounded by Ala residues. The results for the two peptides agree neither with one another nor with the different full-length denatured-state models, which do provide results in good agreement with experiment. This finding is noteworthy because the denatured state of CI2 is very disrupted with little residual structure, such that the peptides might have been expected to serve as reasonable models. Overall the calculations presented here validate our previous MD-generated transition- and denatured-state models and therefore the simulated unfolding pathways and their relevance to refolding.

Alanine↗