Dynamic simulation mechanisms

Progress in the theoretical description of reaction rates in solution of course correlates strongly with that in other theoretical disciplines, in particular those which have profited most from the enonnous advances in computing power such as quantum chemistry and equilibrium as well as non-equilibrium statistical mechanics of liquid solutions where Monte Carlo and molecular dynamics simulations in many cases have taken on the traditional role of experunents, as they allow the detailed investigation of the influence of intra- and intemiolecular potential parameters on the microscopic dynamics not accessible to measurements in the laboratory. No attempt, however, will be made here to address these areas in more than a cursory way, and the interested reader is referred to the corresponding chapters of the encyclopedia. 

At any geometry g.], the gradient vector having components d EjJd Q. provides the forces (F. = -d Ej l d 2.) along each of the coordinates Q-. These forces are used in molecular dynamics simulations which solve the Newton F = ma equations and in molecular mechanics studies which are aimed at locating those geometries where the F vector vanishes (i.e. tire stable isomers and transition states discussed above). 

Prenkel, D. Pree energy computation and first order phase transitions. In Molecular Dynamic Simulation of Statistical Mechanical Systems, Enrico Fermi Summer School, Varenna 1985, G. Ciccotti and W. Hoover, eds. North Holland, Amsterdam (1986) 43-65. 

Field, M.J., Bash, P.A., Karplus, M. A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations. J. Comput. Chem. 11 (1990) 700-733. 

Berendsen. H.J.C., Van Gunsteren, W.F. Practical algorithms for dynamic simulations, in Molecular Dynamics Simulations of Statistical Mechanical Systems, G. Ciccotti, ed., Soc. Italiana di Fisica, Bologna (1987) 43-65. 

Berendsen, H.J.C., Mavri, J. Quantum dynamics simulation of a small quantum system embedded in a classical environment. In Quantum mechanical simulation methods for studying biological systems, D. Bicout and M. Field, eds. Springer, Berlin (1996) 157-179. 

Zigmond, 1988). The ATP-hydrolysis that accompanies actin polymerization, ATP —> ADP + Pj, and the subsequent release of the cleaved phosphate (Pj) are believed to act as a clock (Pollard et ah, 1992 Allen et ah, 1996), altering in a time-dependent manner the mechanical properties of the filament and its propensity to depolymerize. Molecular dynamics simulations suggested a so-called back door mechanism for the hydrolysis reaction ATP ADP - - Pj in which ATP enters the actin from one side, ADP leaves from the same side, but Pj leaves from the opposite side, the back door (Wriggers and Schulten, 1997b). This hypothesis can explain the effect of the toxin phalloidin which blocks the exit of the putative back door pathway and, thereby, delays Pi release as observed experimentally (Dancker and Hess, 1990). 

The principal idea behind the CSP approach is to use input from Classical Molecular Dynamics simulations, carried out for the process of interest as a first preliminary step, in order to simplify a quantum mechanical calculation, implemented in a subsequent, second step. This takes advantage of the fact that classical dynamics offers a reasonable description of many properties of molecular systems, in particular of average quantities. More specifically, the method uses classical MD simulations in order to determine effective... 

Rappe A K, C J Casewit, K S Colwell, W A Goddard III and W M Skiff 1992. UFF, a Full Periodic Table Force Field for Molecular Mechanics and Molecular Dynamics Simulations. Journal of the American Chemical Society 114 10024-10035. 

The thermodynamic properties that we have considered so far, such as the internal energy, the pressure and the heat capacity are collectively known as the mechanical properties and can be routinely obtained from a Monte Carlo or molecular dynamics simulation. Other thermodynamic properties are difficult to determine accurately without resorting to special techniques. These are the so-called entropic or thermal properties the free energy, the chemical potential and the entropy itself. The difference between the mechanical emd thermal properties is that the mechanical properties are related to the derivative of the partition function whereas the thermal properties are directly related to the partition function itself. To illustrate the difference between these two classes of properties, let us consider the internal energy, U, and the Fielmholtz free energy, A. These are related to the partition function by ... 

Field M J, P A Bash and M Karplus 1990. A Combined Quantum Mechanical and Molecular Mechanical Potential for Molecular Dynamics Simulations. Journal of Computational Chemistry 11 700-733. 

Molecular dynamics is a simulation of the time-dependent behavior of a molecular system, such as vibrational motion or Brownian motion. It requires a way to compute the energy of the system, most often using a molecular mechanics calculation. This energy expression is used to compute the forces on the atoms for any given geometry. The steps in a molecular dynamics simulation of an equilibrium system are as follows ... 

Compute the forces on each atom from the energy expression. This is usually a molecular mechanics force held designed to be used in dynamical simulations. 

Before running a molecular dynamics simulation with solvent and a molecular mechanics method, choose the appropriate dielectric constant. You specify the type and value of the dielectric constant in the Force Field Options dialog box. The dielectric constant defines the screening effect of solvent molecules on nonbonded (electrostatic) interactions. 

You can use any ab initio SCF calculation and all the semi-empiri-cal methods, except Extended Hiickel, for molecular dynamics simulations. The procedures and considerations are similar for simulations using molecular mechanics methods (see Molecular Dynamics on page 69). 

Molecular Dynamics and Monte Carlo Simulations. At the heart of the method of molecular dynamics is a simulation model consisting of potential energy functions, or force fields. Molecular dynamics calculations represent a deterministic method, ie, one based on the assumption that atoms move according to laws of Newtonian mechanics. Molecular dynamics simulations can be performed for short time-periods, eg, 50—100 picoseconds, to examine localized very high frequency motions, such as bond length distortions, or, over much longer periods of time, eg, 500—2000 ps, in order to derive equiUbrium properties. It is worthwhile to summarize what properties researchers can expect to evaluate by performing molecular simulations ... 

Combined Quantum and Molecular Mechanical Simulations. A recentiy developed technique is one wherein a molecular dynamics simulation includes the treatment of some part of the system with a quantum mechanical technique. This approach, QM/MM, is similar to the coupled quantum and molecular mechanical methods introduced by Warshel and Karplus (45) and at the heart of the MMI, MMP2, and MM3 programs by AUinger (60). These latter programs use quantum mechanical methods to treat the TT-systems of the stmctures in question separately from the sigma framework. 

---