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 847 records · Page 47Linked to original sources
PubMed2004

Unveiling functional protein motions with picosecond x-ray crystallography and molecular dynamics simulations.

A joint analysis of all-atom molecular dynamics (MD) calculations and picosecond time-resolved x-ray structures was performed to gain single-molecule insights into mechanisms of protein function. Ensemble-averaged MD simulations of the L29F mutant of myoglobin after ligand dissociation reproduce the direction, amplitude, and time scales of crystallographically determined structural changes. This close agreement with experiments at comparable resolution in space and time validates the individual MD trajectories. From 1,700 single-molecule trajectories, we identified and structurally characterized a conformational switch that directs dissociated ligands to one of two nearby protein cavities. Subsequent ligand migration proceeds through a network of transiently interconnected internal cavities, with passage between them involving correlated protein-ligand motions. The simulations also suggest how picosecond protein motions modulate the functional dissociation of oxygen and suppress the geminate recombination of toxic carbon monoxide.

Carbon Monoxide↗
PubMed2006

Changes in the level of poly(Phe) synthesis in Escherichia coli ribosomes containing mutants of L4 ribosomal protein from Thermus thermophilus can be explained by structural changes in the peptidyltransferase center: a molecular dynamics simulation analysis.

Data from polyphenylalanine [poly(Phe)] synthesis determination in the presence and in the absence of erythromycin have been used in conjunction with Molecular Dynamics Simulation analysis, in order to localize the functional sites affected by mutations of Thermus thermophilus ribosomal protein L4 incorporated in Escherichia coli ribosomes. We observed that alterations in ribosome capability to synthesize poly(Phe) in the absence of erythromycin were mainly correlated to shifts of A2062 and C2612 of 23S rRNA, while in the presence of erythromycin they were correlated to shifts of A2060 and U2584 of 23S rRNA. Our results suggest a means of understanding the role of the extended loop of L4 ribosomal protein in ribosomal peptidyltransferase center.

Amino Acid Sequence↗
PubMed2003

Molecular dynamics simulations of phospholipid bilayers with cholesterol.

To investigate the microscopic interactions between cholesterol and lipids in biological membranes, we have performed a series of molecular dynamics simulations of large membranes with different levels of cholesterol content. The simulations extend to 10 ns, and were performed with hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers. The bilayers contain 1024 lipids of which 0-40% were cholesterol and the rest DPPC. The effects of cholesterol on the structure and mesoscopic dynamics of the bilayer were monitored as a function of cholesterol concentration. The main effects observed are a significant ordering of the DPPC chains (as monitored by NMR type order parameters), a reduced fraction of gauche bonds, a reduced surface area per lipid, less undulations--corresponding to an increased bending modulus for the membrane, smaller area fluctuations, and a reduced lateral diffusion of DPPC-lipids as well as cholesterols.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2001

Key interactions in the immunoglobulin-like structure of apo-neocarzinostatin: evidence from nuclear magnetic resonance relaxation data and molecular dynamics simulations.

The three-dimensional structure of apo-neocarzinostatin (apo-NCS, MW: ca.11000, antitumoral chromophore carrier protein) is based on a seven-stranded antiparallel beta-sandwich, very similar to the immunoglobulin folding domain. We investigated the backbone dynamics of apo-NCS by (13)C-NMR relaxation measurements and molecular dynamics simulation. Model-free parameters determined from the experimental data are compared with a 1.5-nsec molecular simulation of apo-NCS in aqueous solution. This comparison provides an accurate description of both local and collective movements within the protein. This analysis enabled us to correlate dynamic processes with key interactions of this beta-protein. Local motions that could be relevant for the intermolecular association with the ligand are also described.

Antibiotics, Antineoplastic↗
PubMed1997

Solution structure of alpha-conotoxin MI determined by 1H-NMR spectroscopy and molecular dynamics simulation with the explicit solvent water.

The conformation of alpha-conotoxin MI, a potent antagonist of the nicotinic acetylcholine receptor, has been investigated in aqueous solution. Two-dimensional NMR experiments and simulated annealing calculations provide the overall topology of alpha-conotoxin MI; then molecular dynamics simulation with the explicit solvent water was followed in order to obtain a more reliable solution structure. The resulting conformation indicates the presence of a 3(10) helix and a type I beta-turn for residues Pro6-Cys8 and Gly9-Try12, respectively, and shows a significant structural similarity to that of alpha-conotoxin GI, which has biological activity similar to that of MI. The present study provides a molecular basis for the alpha-conotoxin-receptor interaction.

Computer Simulation↗
PubMed2004

Comparison of kinetic Monte Carlo and molecular dynamics simulations of diffusion in a model glass former.

We present results from kinetic Monte Carlo (KMC) simulations of diffusion in a model glass former. We find that the diffusion constants obtained from KMC simulations have Arrhenius temperature dependence, while the correct behavior, obtained from molecular dynamics simulations, can be super-Arrhenius. We conclude that the discrepancy is due to undersampling of higher-lying local minima in the KMC runs. We suggest that the relevant connectivity of minima on the potential energy surface is proportional to the energy density of the local minima, which determines the "inherent structure entropy." The changing connectivity with potential energy may produce a correlation between dynamics and thermodynamics.

Journal Article↗
PubMed2000

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer: influence of the crystal environment.

Molecular dynamics (MD) simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure. Two crystal simulations were carried out; one consisted of one unit cell containing two duplexes, and the other of two unit cells containing four duplexes. Two solution simulations were also carried out, one starting from canonical B-DNA and the other starting from the crystal structure. For many helicoidal parameters, the results from the crystal and solution simulations were essentially identical. However, for other parameters, in particular, alpha, gamma, delta, (epsilon - zeta), phase, and helical twist, differences between crystal and solution simulations were apparent. Notably, during crystal simulations, values of helical twist remained comparable to those in the crystal structure, to include the sequence-dependent differences among base steps, in which values ranged from 20 degrees to 50 degrees per base step. However, in the solution simulations, not only did the average values of helical twist decrease to approximately 30 degrees per base step, but every base step was approximately 30 degrees, suggesting that the sequence-dependent information may be lost. This study reveals that MD simulations of the crystal environment complement solution simulations in validating the applicability of MD to the analysis of DNA structure.

Computer Simulation↗
PubMed2006

Car-parrinello molecular dynamics simulation of liquid formic acid.

First-principles molecular dynamics has been used to investigate the structural, vibrational, and energetic properties of formic acid, formic acid-formate anion dimers, and liquid formic acid in a periodically repeated box with 32 formic acid molecules. We found that in liquid formic acid the hydrogen-bonded clusters mainly consist of linear branching chains. From our simulation, we got good agreement with the available structural and dynamical data. We also studied the proton transfer in the cis-formic acid-formate anion dimer, and we showed that this proton transfer does not have any potential barrier. The hydrogen bonding statistics as well as the mean lifetime of the hydrogen bonds are analyzed.

Anions↗
PubMed1982

Jumping frequencies in membrane channels. Comparison between stochastic molecular dynamics simulation and rate theory.

Permeation of molecules through membrane channels involves local interactions with a limited number of ligand groups. A method for the molecular dynamics simulation of particle movement in small ligand systems is described. It is assumed that the ligand groups carry out thermal vibrations, whereas the rest of the channel molecule and the surroundings act as a heat bath which is coupled via random forces to the motions of the ligands. The simulation technique is applied to a simple system which contains some of the essential features influencing jumping rates in membrane channels, such as flexibility of ligand configuration or inertial effects in the motion of the ligands. Since the simulation is based on strictly microscopic parameters of the particle-ligand system, a rigorous test of the predictions of rate theory is possible. It is found that rate theory describes the general dependence of jumping frequency k' on temperature and on ligand binding strength rather well, although the values of k' obtained by computer simulation are 2-3 times smaller than those predicted by rate theory.

Biological Transport, Active↗
PubMed2000

Molecular dynamics simulation of the laser disintegration of aerosol particles.

The mechanisms of disintegration of submicrometer particles irradiated by short laser pulses are studied by a molecular dynamics simulation technique. Simulations at different laser fluences are performed for particles with homogeneous composition and particles with transparent inclusions. Spatially nonuniform deposition of laser energy is found to play a major role in defining the character and the extent of disintegration. The processes that contribute to the disintegration include overheating and explosive decomposition of the illuminated side of the particle, spallation of the backside of large particles, and disruption of the transparent inclusion caused by the relaxation of the laser-induced pressure. The observed mechanisms are related to the nature of the disintegration products and implications of the simulation results for aerosol time-of-flight mass spectrometry are discussed. Application of multiple laser pulses is predicted to be advantageous for efficient mass spectrometry sampling of aerosols with a large size to laser penetration depth ratio.

Aerosols↗
PubMed2006

Calculation of the surface tension of oxygen using molecular-dynamics simulations.

The surface tension of oxygen at the liquid-vapor interface is calculated for the temperature range of 60-90 K using molecular-dynamics simulations and is shown to be within 1.0% error of experimental values for most of the temperatures studied. The potential used here is the same as in our previous study on liquid oxygen alone [S. D. Bembenek and B. M. Rice, J. Chem. Phys. 113, 2354 (2000)] and was optimized with an innovative statistical-mechanical method. The potential does not use a long-range cutoff nor a tail correction, which are usually considered necessary to obtain accurate values for the surface tension. We reason that the accuracy in surface tension is directly related to our parametrization method for the potential.

Journal Article↗
PubMed2006

A novel algorithm to model the influence of host lattice flexibility in molecular dynamics simulations: loading dependence of self-diffusion in carbon nanotubes.

We describe a novel algorithm that includes the effect of host lattice flexibility into molecular dynamics simulations that use rigid lattices. It uses a Lowe-Andersen thermostat for interface-fluid collisions to take the most important aspects of flexibility into account. The same diffusivities and other properties of the flexible framework system are reproduced at a small fraction of the computational cost of an explicit simulation. We study the influence of flexibility on the self-diffusion of simple gases inside single walled carbon nanotubes. Results are shown for different guest molecules (methane, helium, and sulfur hexafluoride), temperatures, and types of carbon nanotubes. We show, surprisingly, that at low loadings flexibility is always relevant. Notably, it has a crucial influence on the diffusive dynamics of the guest molecules.

Journal Article↗
PubMed1997

Effect of warmup protocol and sampling time on convergence of molecular dynamics simulations of a DNA dodecamer using AMBER 4.1 and particle-mesh Ewald method.

This report describes one 3000 ps and two 1500 ps molecular dynamic simulations on a TATA box containing dodecamer DNA duplex in a periodic box of TIP3P water molecules, using the AMBER 4.1 implementation of the particle-mesh Ewald method. We compare the effect of warmup protocol and simulation time length on the root-mean square deviation (RMSD) parameter. For the longer simulation, the RMSD computed for the 500-1000 ps time interval is representative of longer time intervals, including 500-3000 ps. The various warmup protocols do not appear to have a significant effect on the simulation results. Based on the present results, DNA sequence-dependent differences in RMSD, or related properties, should exceed two standard deviations before being attributed to non-simulation factors, such as warmup protocol and sampling time effects; we recommend a minimum criterion of at least a three standard deviation difference with a sampling period of at least 500-1000 ps. In addition, while end effects appear negligible there is a consistent dependence of RMSD on DNA helix length.

Adenoviridae↗
PubMed2005

Molecular dynamic simulation of disorder induced amorphization in pyrochlore.

The defect accumulation mechanism of amorphization has been studied for the La2Zr2O7 pyrochlore by means of classical molecular dynamic simulations. Present calculations show that the accumulation of cation Frenkel pairs is the main driving parameter for the amorphization process, while the oxygen atoms simply rearrange around cations. Under Frenkel pair accumulation, the structure follows the pyrochlore-fluorite-amorphous sequence. Present results consequently provide atomic-level interpretation to previous experimental irradiation observations of the two-step phase transition.

Journal Article↗
PubMed2000

The structure and dynamics of ACTH (1-10) on the surface of a sodium dodecylsulfate (SDS) micelle: a molecular dynamics simulation study.

ACTH (1-10), an adrenocorticotropin hormone fragment, was studied by molecular dynamics (MD) simulation in the NPT ensemble in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated during the equilibration stage of the simulation. The analyses of the trajectories from the subsequent unrestrained MD showed that ACTH (1-10) does not conform to a helical structure at the micelle-water interface; however, the structure is amphipathic. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide hydrogens of the peptide and the micelle head-groups. ACTH (1-10) was found to lie on the surface of the SDS micelle. Most of the hydrophobic interactions came from the side-chains of Met-4, Phe-7 and Trp-9. The peptide bonds were either hydrated or involved in intramolecular hydrogen bonding. Decreased hydration for the backbone of His-6 and Phe-7 was due to intermolecular hydrogen bonding with the SDS head-groups. The time correlation functions of the N-H bonds of the peptide in water and in the micelle showed that the motions of the peptide, except for the N- and C-termini, are significantly reduced when partitioned in the micelle.

Adrenocorticotropic Hormone↗
PubMed2003

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies.

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.

Carboxylic Ester Hydrolases↗
PubMed2002

Molecular dynamics simulations of a double unit cell in a protein crystal: volume relaxation at constant pressure and correlation of motions between the two unit cells.

Eight molecular dynamics simulations of a double crystal unit cell of ubiquitin were performed to investigate the effects of simulating at constant pressure and of simulating two unit cells compared to a single unit cell. To examine the influence of different simulation conditions, the constant-pressure and constant-volume simulations were each performed with and without counterions and using two different treatments of the long-range electrostatic interactions (lattice-sum and reaction-field methods). The constant-pressure simulations were analyzed in terms of unit cell deformation and accompanying protein deformations. Energetic and structural properties of the proteins in the simulations of the double unit cell were compared to the results of previously reported one-unit-cell simulations. Correlation between the two unit cells was also investigated based on relative translational and rotational movements of the proteins and on dipole fluctuations. The box in the constant-pressure simulations is found to deform slowly to reach convergence only after 5-10 ns. This deformation does not result from a distortion in the structure of the proteins but rather from changes in protein packing within the unit cell. The results of the double-unit-cell simulations are closely similar to the results of the single-unit-cell simulations, and little motional correlation is found between the two unit cells.

Chlorides↗
PubMed1996

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Muramidase↗