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 109 records · Page 6Linked to original sources
PubMed2005

Fractal dimensions of silica gels generated using reactive molecular dynamics simulations.

We have used molecular dynamics simulations based on a three-body potential with charge transfer to generate nanoporous silica aerogels. Care was taken to reproduce the sol-gel condensation reaction that forms the gel backbone as realistically as possible and to thereby produce credible gel structures. The self-similarity of aerogel structures was investigated by evaluating their fractal dimension from geometric correlations. For comparison, we have also generated porous silica glasses by rupturing dense silica and computed their fractal dimension. The fractal dimension of the porous silica structures was found to be process dependent. Finally, we have determined that the effect of supercritical drying on the fractal nature of condensed silica gels is not appreciable.

Journal Article↗
PubMed2005

Time reversible and symplectic integrators for molecular dynamics simulations of rigid molecules.

Molecular dynamics integrators are presented for translational and rotational motion of rigid molecules in microcanonical, canonical, and isothermal-isobaric ensembles. The integrators are all time reversible and are also, in some approaches, symplectic for the microcanonical ensembles. They are developed utilizing the quaternion representation on the basis of the Trotter factorization scheme using a Hamiltonian formalism. The structure is similar to that of the velocity Verlet algorithm. Comparison is made with standard integrators in terms of stability and it is found that a larger time step is stable with the new integrators. The canonical and isothermal-isobaric molecular dynamics simulations are defined by using a chain thermostat approach according to generalized Nosé-Hoover and Andersen methods.

Journal Article↗
PubMed1995

Rigid-body motions of sub-units in DNA: a correlation analysis of a 200 ps molecular dynamics simulation.

A 200 ps molecular dynamics simulation of the B-form double stranded self-complementary octanucleotide d(CTGATCAG) is analyzed in terms of correlated motions using the canonical analysis approach. Each nucleotide is decomposed in three sub-units corresponding to the base, the sugar ring and the backbone respectively. The correlation between the full dynamics of two sub-units was found to decrease as their mutual distance increases. The interpretation of the full dynamics of sub-units as the superimposition of rigid-body motions (translation and orientation) and deformation shows that the main source of correlation is rigid-body motions. Correlation between sub-units deformation is weak and practically vanishes for sub-units belonging to non-adjacent nucleotides. It is also shown that the correlation is much more important for sub-units of the same strand than of opposite strands. We conclude that the internal dynamics of the octanucleotide may be well described by rigid-body motions, the sub-units deformation having only local influence whereas sub-units translation and rotation have repercussion to long distances. The results presented in this study suggest how the number of degrees of freedom may be reduced for simulating long-time dynamics of oligonucleotides.

Computer Simulation↗
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↗
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↗