SEARCH · PubMed Health

Results for “Molecular Dynamics Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 91 records · Page 5Linked to original sources

PubMed1989

Superoxide dismutase: fluctuations in the structure and solvation of the active site channel studied by molecular dynamics simulation.

The molecular dynamics (MD) simulation of superoxide dismutase (SOD) in water is carried out for a total of 23 ps. The simulation system is a 26 A sphere centered at the active site of SOD, including 1602 atoms from SOD and 1761 water molecules. There is no gross deviation from the x-ray structure for the average MD structure. The structure and potential fluctuations around the active site are examined. The results provide new insight to the interactions between SOD and its substrate superoxide.

Binding Sites↗

PubMed1999

Structural and dynamic differences of the estrogen receptor DNA-binding domain, binding as a dimer and as a monomer to DNA: molecular dynamics simulation studies.

Molecular dynamics (MD) simulations of the estrogen receptor DNA-binding domain (ERDBD) as a dimer in complex with its DNA response element (ERE) show a significant difference in both structure and dynamics, compared to a MD simulation of monomeric ERDBD bound to its half-site response element (EREH). The C-terminal zinc binding domain (ZnII), including a region (helix II) which is in a helical conformation in ERE-(ERDBD)2, is considerably more flexible in EREH-ERDBD than in the dimeric complex. In EREH-ERDBD, all helical hydrogen bonds in helix II are broken and the entire ZnII region is detached from a hydrogen bonding network that in ERE-(ERDBD)2 connects to other parts of the protein as well as to the DNA. The regions that become flexible in EREH-ERDBD are identical to the regions where the NMR solution structure of free ERDBD is poorly ordered. This strongly suggests that dimerisation of ERDBD is required for ordering of the ZnII region and that monomeric binding to DNA is not sufficient for the ordering. This contrasts to the glucocorticoid receptor DNA-binding domain (GRDBD) which has essentially the same mobility (uniform and limited), regardless of whether it is free as a monomer in solution, bound as a monomer to its half-site response element or in a dimeric complex with the full response element. The hydrogen bonding network that connects ZnII with other parts of the protein and to DNA is almost identical in ERDBD and GRDBD. However, in GRDBD there is also a serine (in the N-terminal zinc coordinating region) with a central role in this network, connecting to the ZnII region. This serine is replaced by a glycine in ERDBD and we suggest that this substitution is sufficient for destabilisation of the network, thus leading to a more flexible ZnII region, which becomes ordered first upon forming a complex with another ERDBD and DNA.

Amino Acid Sequence↗

PubMed2006

Local density inhomogeneities and dynamics in supercritical water: A molecular dynamics simulation approach.

Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.

Journal Article↗

PubMed2005

Odd-even variations in the wettability of n-alkanethiolate monolayers on gold by water and hexadecane: a molecular dynamics simulation study.

Molecular dynamics (MD) simulations were performed to investigate odd-even chain length dependencies in the wetting properties of self-assembled monolayers (SAMs) of n-alkanethiols [CH3(CH2)n-1SH] on gold by water and hexadecane. Experimentally, the contact angle of hexadecane on the SAMs depends on whether n is odd or even, while contact angles for water show no odd-even dependence. Our MD simulations of this system included a microscopic droplet of either 256 water molecules or 60 hexadecane molecules localized on an n-alkanethiolate SAM on gold with either an even or odd chain length. Contact angles calculated for these nanoscopic droplets were consistent with experimentally observed macroscopic trends in wettability, namely, that hexadecane is sensitive to structural differences between odd- and even-chained SAMs while water is not. Structural properties for the SAMs (including features such as chain tilt, chain twist, and terminal methyl group tilt) were calculated during the MD simulations and used to generate IR spectra of these films that compared favorably with experimental spectra. MD simulations of SAMs in contact with slabs of water and hexadecane revealed that the effects of these solvents on the structure of the SAM was restricted to the chain terminus and had no effect on the inner structure of the SAM. The density profiles for water and hexadecane on the SAMs were different in that water displayed a significant depletion in its density at the liquid/SAM interface from its bulk value, while no such depletion occurred for hexadecane. This difference in contact may explain the lack of an odd-even variation in the wetting characteristics of water on these surfaces, because the water molecules are positioned further away from the surface and, therefore, are not sensitive to the structural differences in the average orientations for the terminal methyl groups in odd- and even-chained SAMs. In contrast, the differences in the wetting properties of hexadecane on the odd- and even-chained SAMs may reflect the closer proximity of these molecules to the SAM surface and a resulting greater sensitivity to the differences in the terminal methyl group orientations in the SAMs. SAM-solvent interaction energies were calculated during the MD simulations, yielding interaction energies that differed on the even- and odd-chained surfaces by approximately 10% for hexadecane and negligibly for water, in accord with estimates using experimental wetting results.

Alkanes↗

PubMed2003

Test of mean-field equations for two types of hard-sphere systems by a Brownian-dynamics simulation and a molecular-dynamics simulation.

A mean-field nonlinear equation for the mean-square displacement, recently proposed by one of the present authors [M. Tokuyama, Phys. Rev. E 62, R5915 (2000); Physica A 289, 57 (2001)], for concentrated, equilibrium suspensions of hard spheres is extended to describe equilibrium atomic systems of hard spheres. The validity of two types of mean-field equations is investigated by two kinds of computer simulations; a Brownian-dynamics simulation on suspensions of hard spheres and a molecular-dynamics simulation on atomic systems of hard spheres. A good agreement between the mean-field equations and simulations is then shown for different volume fractions. The two types of model systems of hard spheres are thus shown to be identical to each other on the study of the liquid-solid transition. However, analyses suggest that a new interaction is indispensable to understand the mechanism for the liquid-glass transition in both systems.

Journal Article↗

PubMed2002

Calculating the bulk modulus for a lipid bilayer with nonequilibrium molecular dynamics simulation.

Nonequilibrium molecular dynamics (NEMD) computer simulations are used to calculated the bulk modulus for a dimyristoylphosphatidylcholine bilayer. A methodology is developed whereby NEMD can be effectively used to calculate material properties for complex systems that undergo long time-scale conformational changes. It is found that the bulk modulus upon expansion from a zero stress state agrees well with experimental estimates. However, it is also found that the modulus upon contraction from a zero stress state is larger. From a molecular perspective, it is possible to explain this phenomena by examining the molecular origins of the pressure response. The finding that the two moduli are not equal upon compression and expansion is in apparent contradiction to osmotic stress experiments where the area modulus was found to be the same upon expansion and contraction. This issue is addressed.

Algorithms↗

PubMed1996

H-bond stability in the tRNA(Asp) anticodon hairpin: 3 ns of multiple molecular dynamics simulations.

Multiple molecular dynamics trajectories of the solvated and neutralized 17-residue tRNA(Asp) anticodon hairpin were generated for a total of 3 ns. Explicit treatment of all long-ranged electrostatic interactions by the particle mesh Ewald algorithm, as implemented in the AMBER MD software package, effected a degree of structural stabilization not previously achieved by use of a long 16-A solvent interaction truncation scheme. The increased stability of this multiple molecular dynamics set was appropriate for an in-depth analysis of the six 500-ps-long trajectories and allowed the characterization of a number of key structural interactions. The dynamical behavior of the standard Watson-Crick base pairs, the noncanonical G30-U40 "wobble" base pair, and the psi 32-C38 pseudo-base pair is presented as well as that of two C--H... O hydrogen bonds found to contribute to the array of tertiary interactions that stabilize the seven-nucleotide native loop conformation. The least mobile residue in the loop is U33, which forms the U-turn motif and which participates in several hydrogen-bonding interactions, whereas the most mobile residue is the apical residue G34 at the wobble position, a factor undoubtedly important in its biological function. The set of multiple molecular dynamics trajectories obtained does not converge on a 500-ps time scale to a unique dynamical model but instead describes an ensemble of structural microstates accessible to the system under the present simulation protocol, which is the result of local structural heterogeneity rather than of global conformational changes.

Anticodon↗

PubMed2004

Calculating potentials of mean force from steered molecular dynamics simulations.

Steered molecular dynamics (SMD) permits efficient investigations of molecular processes by focusing on selected degrees of freedom. We explain how one can, in the framework of SMD, employ Jarzynski's equality (also known as the nonequilibrium work relation) to calculate potentials of mean force (PMF). We outline the theory that serves this purpose and connects nonequilibrium processes (such as SMD simulations) with equilibrium properties (such as the PMF). We review the derivation of Jarzynski's equality, generalize it to isobaric--isothermal processes, and discuss its implications in relation to the second law of thermodynamics and computer simulations. In the relevant regime of steering by means of stiff springs, we demonstrate that the work on the system is Gaussian-distributed regardless of the speed of the process simulated. In this case, the cumulant expansion of Jarzynski's equality can be safely terminated at second order. We illustrate the PMF calculation method for an exemplary simulation and demonstrate the Gaussian nature of the resulting work distribution.

Computer Simulation↗

PubMed2003

Dynamic molecules: molecular dynamics for everyone. An internet-based access to molecular dynamic simulations: basic concepts.

Molecular dynamics is a rapidly developing field of science and has become an established tool for studying the dynamic behavior of biomolecules. Although several high quality programs for performing molecular dynamic simulations are freely available, only well-trained scientists are currently able to make use of the broad scientific potential that molecular dynamic simulations offer to gain insight into structural questions at an atomic level. The "Dynamic Molecules" approach is the first internet portal that provides an interactive access to set up, perform and analyze molecular dynamic simulations. It is completely based on standard web technologies and uses only publicly available software. The aim is to open molecular dynamics techniques to a broader range of users including undergraduate students, teachers and scientists outside the bioinformatics field. The time-limiting factors are the availability of free capacity on the computing server to run the simulations and the time required to transport the history file through the internet for the animation mode. The interactive access mode of the portal is acceptable for animations of molecules having up to about 500 atoms.

Computer Simulation↗

PubMed2004

Closing loop base pairs in RNA loop-loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations.

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.

Aptamers, Nucleotide↗

PubMed2002

Early steps of the intramolecular signal transduction in rhodopsin explored by molecular dynamics simulations.

We present molecular dynamics simulations of bovine rhodopsin in a membrane mimetic environment based on the recently refined X-ray structure of the pigment. The interactions between the protonated Schiff base and the protein moiety are explored both with the chromophore in the dark-adapted 11-cis and in the photoisomerized all-trans form. Comparison of simulations with Glu181 in different protonation states strongly suggests that this loop residue located close to the 11-cis bond bears a negative charge. Restrained molecular dynamics simulations also provide evidence that the protein tightly confines the absolute conformation of the retinal around the C12-C13 bond to a positive helicity. 11-cis to all-trans isomerization leads to an internally strained chromophore, which relaxes after a few nanoseconds by a switching of the ionone ring to an essentially planar all-trans conformation. This structural transition of the retinal induces in turn significant conformational changes of the protein backbone, especially in helix VI. Our results suggest a possible molecular mechanism for the early steps of intramolecular signal transduction in a prototypical G-protein-coupled receptor.

Animals↗

PubMed2003

Structure and dynamics of interfacial water in an Lalpha phase lipid bilayer from molecular dynamics simulations.

Based on molecular dynamics simulations, an analysis of structure and dynamics is performed on interfacial water at a liquid crystalline dipalmitoylphosphatidycholine/water system. Water properties relevant for understanding NMR relaxation are emphasized. The first and second rank orientational order parameters of the water O-H bonds were calculated, where the second rank order parameter is in agreement with experimental determined quadrupolar splittings. Also, two different interfacial water regions (bound water regions) are revealed with respect to different signs of the second rank order parameter. The water reorientation correlation function reveals a mixture of fast and slow decaying parts. The fast (ps) part of the correlation function is due to local anisotropic water reorientation whereas the much slower part is due to more complicated processes including lateral diffusion along the interface and chemical exchange between free and bound water molecules. The 100-ns-long molecular dynamics simulation at constant pressure (1 atm) and at a temperature of 50 degrees C of 64 lipid molecules and 64 x 23 water molecules lack a slow water reorientation correlation component in the ns time scale. The (2)H(2)O powder spectrum of the dipalmitoylphosphatidycholine/water system is narrow and consequently, the NMR relaxation time T(2) is too short compared to experimental results.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2003

The mechanism of vesicle fusion as revealed by molecular dynamics simulations.

We describe molecular dynamics simulations elucidating the molecular details of the process of fusion for small lipid vesicles. The simulations are based on a coarse grained (CG) lipid model that accurately represents the lamellar state of a variety of phospholipids and enables us to observe intermediate stages during fusion at near atomic detail. Simulations were conducted on a variety of systems containing common phospholipids such as phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysoPC, and mixtures of the above. The fusion intermediates found are in general agreement with the stalk-pore mechanism. Transient pores sometimes form adjacent to the stalk, however, resulting in the mixing of lipids from the outer and inner monolayers. The speed of stalk formation and the opening of the fusion pore can be modulated by altering the lipid composition in qualitative agreement with experimental observations.

Cytoplasmic Vesicles↗

PubMed2002

Mechanism of acetylcholinesterase inhibition by fasciculin: a 5-ns molecular dynamics simulation.

Our previous molecular dynamics simulation (10 ns) of mouse acetylcholinesterase (EC 3.1.1.7) revealed complex fluctuation of the enzyme active site gorge. Now we report a 5-ns simulation of acetylcholinesterase complexed with fasciculin 2. Fasciculin 2 binds to the gorge entrance of acetylcholinesterase with excellent complementarity and many polar and hydrophobic interactions. In this simulation of the protein-protein complex, where fasciculin 2 appears to sterically block access of ligands to the gorge, again we observe a two-peaked probability distribution of the gorge width. When fasciculin is present, the gorge width distribution is altered such that the gorge is more likely to be narrow. Moreover, there are large increases in the opening of alternative passages, namely, the side door (near Thr 75) and the back door (near Tyr 449). Finally, the catalytic triad arrangement in the acetylcholinesterase active site is disrupted with fasciculin bound. These data support that, in addition to the steric obstruction seen in the crystal structure, fasciculin may inhibit acetylcholinesterase by combined allosteric and dynamical means. Additional data from these simulations can be found at http://mccammon.ucsd.edu/.

Acetylcholinesterase↗

PubMed2003

Molecular dynamics simulations of metalloproteins.

Molecular dynamics simulations are now commonly applied to metalloproteins, despite the challenges introduced by the presence of metal ions. Force field parameters are nowadays available also for these 'exotic' atoms and several biological systems have been successfully studied. Some of the most relevant results and methodological advancements are reviewed.

Binding Sites↗

PubMed2002

Molecular dynamics simulations of biomolecules.

Molecular dynamics simulations are important tools for understanding the physical basis of the structure and function of biological macromolecules. The early view of proteins as relatively rigid structures has been replaced by a dynamic model in which the internal motions and resulting conformational changes play an essential role in their function. This review presents a brief description of the origin and early uses of biomolecular simulations. It then outlines some recent studies that illustrate the utility of such simulations and closes with a discussion of their ever-increasing potential for contributing to biology.

Models, Molecular↗

PubMed2001

Evidence of reduced flexibility in disulfide bridge-depleted azurin: a molecular dynamics simulation study.

Two molecular dynamics simulations have been performed for 2 ns, at room temperature, on fully hydrated wild type and Cys3Ala/Cys26Ala double-mutant azurin, to investigate the role of the unique disulfide bridge on the structure and dynamics of the protein. The results show that the removal of the [bond]SS[bond] bond does not affect the structural features of the protein, whereas alterations of the dynamical properties are observed. The root mean square fluctuations of the atomic positions are, on average, considerably reduced in the azurin mutant with respect to the wild type form. The number of intramolecular hydrogen bonds between protein backbone atoms that are lost during the simulation, with respect to the starting configuration, are reduced in the absence of the disulfide bond. The analysis of the dynamical cross-correlation map, characterising the protein co-ordinated internal motions, demonstrates in the mutated azurin a significant decrease in anti-correlated displacements between protein residues, with the only exception occurring in the region of the mutation sites. The overall findings show a relevant reduction in flexibility as a consequence of the disulfide bridge depletion in azurin, suggesting that the [bond]SS[bond] bond is a structural element which significantly contributes to the dynamic properties of the native protein.

Azurin↗

PubMed1999

Structure of dipalmitoylphosphatidylcholine/cholesterol bilayer at low and high cholesterol concentrations: molecular dynamics simulation.

By using molecular dynamics simulation technique we studied the changes occurring in membranes constructed of dipalmitoylphosphatidylcholine (DPPC) and cholesterol at 8:1 and 1:1 ratios. We tested two different initial arrangements of cholesterol molecules for a 1:1 ratio. The main difference between two initial structures is the average number of nearest-neighbor DPPC molecules around the cholesterol molecule. Our simulations were performed at constant temperature (T = 50 degrees C) and pressure (P = 0 atm). Durations of the runs were 2 ns. The structure of the DPPC/cholesterol membrane was characterized by calculating the order parameter profiles for the hydrocarbon chains, atom distributions, average number of gauche defects, and membrane dipole potentials. We found that adding cholesterol to membranes results in a condensing effect: the average area of membrane becomes smaller, hydrocarbon chains of DPPC have higher order, and the probability of gauche defects in DPPC tails is lower. Our results are in agreement with the data available from experiments.

1,2-Dipalmitoylphosphatidylcholine↗