Solvation, molecular dynamics

Duan Y., Wang L. and Kouman P. A. The early stage of folding of villin headpiece sub-domain observed in a 200-nanosecond fully solvated molecular dynamics simulation. Proc. Natl. Acad. Sci., USA (1998) 95(17) 9897-9902. 

Stage of Folding of Villin Headpiece Subdomain Observed in a 200-Nanosecond Fnlly Solvated Molecular Dynamic Simulation. 

While numerous works of Monte Carlo simulations were carried out for the study of the clay-water system, especially in terms of clay swelling and interlayer cation solvation, molecular dynamics (MD) has only recently been used. While MC has an advantage in conformation sampling in the entire system in equilibrium by its stochastic method, MD is more suited to the dynamical or diffusion properties of the chemical system. Here, several examples of using MD techniques on clay minerals are reviewed for understanding the current activities in this field. 

Keywords liquid-liquid extraction interface ionophore calixarene counterion solvation molecular dynamics immiscible liquids phase separation recognition organization. 

Y. Duan, L. Wang, and P. A. Kollman, Proc. Natl. Acad. Sci. U.S.A., 95, 9897 (1998). The Early Stage of Folding of Villin Headpiece Subdomain Observed in a 200-Nanosecond Fully Solvated Molecular Dynamics Simulation. 

Specific solute-solvent interactions involving the first solvation shell only can be treated in detail by discrete solvent models. The various approaches like point charge models, siipennoleciilar calculations, quantum theories of reactions in solution, and their implementations in Monte Carlo methods and molecular dynamics simulations like the Car-Parrinello method are discussed elsewhere in this encyclopedia. Here only some points will be briefly mentioned that seem of relevance for later sections. 

Schreiber, H., Steinhauser, O. Cutoff size does strongly influence molecular dynamics results on solvated polypeptides. Biochem. 31 (1992) 5856-5860. 

Mavri, J., Berendsen, H.J.C., Van Gunsteren, W.F. Influence of solvent on intramolecular proton transfer in hydrogen malonate. Molecular dynamics study of tunneling by density matrix evolution and nonequilibrium solvation. J. Phys. Chem. 97 (1993) 13469-13476. 

Elamrani et al. 1996] Elamrani, S., Berry, M.B., Phillips Jr., G.N., McCammon, J.A. Study of Global Motions in Proteins by Weighted Masses Molecular Dynamics Adenylate Kinase as a Test Case. Proteins 25 (1996) 79-88 [Elcock et al. 1997] Elcock, A.H., Potter, M.J., McCammon, J.A. Application of Poisson-Boltzmann Solvation Forces to Macromolecular Simulations. In Computer Simulation of Biomoleeular Systems, Vol. 3, A.J. Wilkinson et al. eds., ESCOM Science Publishers B.V., Leiden... 

Abstract. Molecular dynamics (MD) simulations of proteins provide descriptions of atomic motions, which allow to relate observable properties of proteins to microscopic processes. Unfortunately, such MD simulations require an enormous amount of computer time and, therefore, are limited to time scales of nanoseconds. We describe first a fast multiple time step structure adapted multipole method (FA-MUSAMM) to speed up the evaluation of the computationally most demanding Coulomb interactions in solvated protein models, secondly an application of this method aiming at a microscopic understanding of single molecule atomic force microscopy experiments, and, thirdly, a new method to predict slow conformational motions at microsecond time scales. 

To enable an atomic interpretation of the AFM experiments, we have developed a molecular dynamics technique to simulate these experiments [49], Prom such force simulations rupture models at atomic resolution were derived and checked by comparisons of the computed rupture forces with the experimental ones. In order to facilitate such checks, the simulations have been set up to resemble the AFM experiment in as many details as possible (Fig. 4, bottom) the protein-ligand complex was simulated in atomic detail starting from the crystal structure, water solvent was included within the simulation system to account for solvation effects, the protein was held in place by keeping its center of mass fixed (so that internal motions were not hindered), the cantilever was simulated by use of a harmonic spring potential and, finally, the simulated cantilever was connected to the particular atom of the ligand, to which in the AFM experiment the linker molecule was connected. 

D. Beglov and B. Roux. Dominant solvations effects from the primary shell of hydration Approximation for molecular dynamics simulations. Biopolymers, 35 171-178, 1994. 

HyperChem uses th e ril 31 water m odel for solvation. You can place th e solute in a box of T1P3P water m oleeules an d impose periodic boun dary eon dition s. You may then turn off the boundary conditions for specific geometry optimi/.aiion or molecular dynamics calculations. However, th is produces undesirable edge effects at the solvent-vacuum interface. 

Often yon need to add solvent molecules to a solute before running a molecular dynamics simiilatmn (see also Solvation and Periodic Boundary Conditions" on page 62). In HyperChem, choose Periodic Box on the Setup m en ii to enclose a soln te in a periodic box filled appropriately with TIP3P models of water inole-cii les. 

Dang L X, J E Rice, J Caldwell and P A Kollman 1991. Ion Solvation in Polarisable Water Molecular Dynamics Simulations. Journal of the American Chemical Society 113 2481-2486. 

Cheatham T E III, J L Miller, T Fox, T A Darden and P A Kollman 1995. Molecular Dynamics Simulations on Solvated Biomolecular Systems The Particle Mesh Ewald Method Leads to Stable Trajectories of DNA, RNA and Proteins. Journal of the American Chemical Society 117 4193-4194. 

PLS (partial least-squares) algorithm used for 3D QSAR calculations PM3 (parameterization method three) a semiempirical method PMF (potential of mean force) a solvation method for molecular dynamics calculations... 

The idea of a finite simulation model subsequently transferred into bulk solvent can be applied to a macromolecule, as shown in Figure 5a. The alchemical transformation is introduced with a molecular dynamics or Monte Carlo simulation for the macromolecule, which is solvated by a limited number of explicit water molecules and otherwise surrounded by vacuum. Then the finite model is transferred into a bulk solvent continuum... 

For 25 years, molecular dynamics simulations of proteins have provided detailed insights into the role of dynamics in biological activity and function [1-3]. The earliest simulations of proteins probed fast vibrational dynamics on a picosecond time scale. Fifteen years later, it proved possible to simulate protein dynamics on a nanosecond time scale. At present it is possible to simulate the dynamics of a solvated protein on the microsecond time scale [4]. These gains have been made through a combination of improved computer processing (Moore s law) and clever computational algorithms [5]. 

RB Yelle, BW Beck, JB Koerner, CA Sacksteder, T Ichiye. Influence of the metal site on the structure and solvation of nibredoxm and its analogs A molecular dynamics study. Proteins accepted. 

Molecular dynamics simulations have also been used to interpret phase behavior of DNA as a function of temperature. From a series of simulations on a fully solvated DNA hex-amer duplex at temperatures ranging from 20 to 340 K, a glass transition was observed at 220-230 K in the dynamics of the DNA, as reflected in the RMS positional fluctuations of all the DNA atoms [88]. The effect was correlated with the number of hydrogen bonds between DNA and solvent, which had its maximum at the glass transition. Similar transitions have also been found in proteins. 

This chapter has given an overview of the structure and dynamics of lipid and water molecules in membrane systems, viewed with atomic resolution by molecular dynamics simulations of fully hydrated phospholipid bilayers. The calculations have permitted a detailed picture of the solvation of the lipid polar groups to be developed, and this picture has been used to elucidate the molecular origins of the dipole potential. The solvation structure has been discussed in terms of a somewhat arbitrary, but useful, definition of bound and bulk water molecules. 

---