Molecular dynamics simulation practice

How can we apply molecular dynamics simulations practically. This section gives a brief outline of a typical MD scenario. Imagine that you are interested in the response of a protein to changes in the amino add sequence, i.e., to point mutations. In this case, it is appropriate to divide the analysis into a static and a dynamic part. What we need first is a reference system, because it is advisable to base the interpretation of the calculated data on changes compared with other simulations. By taking this relative point of view, one hopes that possible errors introduced due to the assumptions and simplifications within the potential energy function may cancel out. All kinds of simulations, analyses, etc., should always be carried out for the reference and the model systems, applying the same simulation protocols. 

V. MOLECULAR DYNAMICS SIMULATION PRACTICE A. Assigning Initial Values... 

Predicting the solvent or density dependence of rate constants by equation (A3.6.29) or equation (A3.6.31) requires the same ingredients as the calculation of TST rate constants plus an estimate of and a suitable model for the friction coefficient y and its density dependence. While in the framework of molecular dynamics simulations it may be worthwhile to numerically calculate friction coefficients from the average of the relevant time correlation fiinctions, for practical purposes in the analysis of kinetic data it is much more convenient and instructive to use experimentally detemiined macroscopic solvent parameters. 

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. 

We have presented a simple protocol to run MD simulations for systems of interest. There are, however, some tricks to improve the efficiency and accuracy of molecular dynamics simulations. Some of these techniques, which are discussed later in the book, are today considered standard practice. These methods address diverse issues ranging from efficient force field evaluation to simplified solvent representations. 

Chapter 15 gives an extensive and detailed review of theoretical and practical aspects of macromolecular transport in nanostructured media. Chapter 16 examines the change in transport properties of electrolytes confmed in nanostructures, such as pores of membranes. The confinment effect is also analyzed by molecular dynamic simulation. 

This equation means that when there is a free energy difference of a few fcs T the probability P( ) is reduced considerably, that is, those conformations with large A( ) are sampled very rarely. This is a very important observation in terms of numerical efficiency. At the transition region for example, the free energy is maximum and typically very few sample points are obtained during the course of molecular dynamics simulation. In turn this results in very large statistical errors. Those errors can only be reduced by increasing the simulation time, sometimes beyond what is practically feasible. 

The above discussion was based on the results of molecular dynamics simulations on unsaturated or conjugated hydrocarbons. Although the general features can be extended to molecular structures of more general types, in practice it is appropriate to consider the specific form of the electron orbitals involved. For instance, d d transitions in transition metal ion complexes involve orbitals mainly localized on the metal ion that, in the crystal field... 

The credit load for die computational chemistry laboratory course requires that the average student should be able to complete almost all of the work required for the course within die time constraint of one four-hour laboratory period per week. This constraint limits the material covered in the course. Four principal computational methods have been identified as being of primary importance in the practice of chemistry and thus in the education of chemistry students (1) Monte Carlo Methods, (2) Molecular Mechanics Methods, (3) Molecular Dynamics Simulations, and (4) Quantum Chemical Calculations. Clearly, other important topics could be added when time permits. These four methods are developed as separate units, in each case beginning with die fundamental principles including simple programming and visualization, and building to the sophisticated application of the technique to a chemical problem. 

Sen, S. and Nilsson, L. 1999. Some Practical Aspects of Free Energy Calculations from Molecular Dynamics Simulation , J. Comput. Chem., 20, 877. 

Besides the obvious practical importance of solvation dynamics in water, from a fundamental standpoint, water offers a unique opportunity to evaluate theoretical models for solvation. Only for water have extensive molecular dynamics simulations been accomplished (see below). Also semi-empirical models for solvation dynamics, such as MSA, can be carefully examined for water because the necessary information on the dielectric dispersion 2(high level of accuracy. The results of the solvation... 

Figure 2.1 Density of a liquid versus the coordinate normal to its surface (a) is a schematic plot (b) results from molecular dynamics simulations of a n-tridecane (C13H28) at 27°C adapted from Ref. [11]. Tridecane is practically not volatile. For this reason the density in the vapor phase is negligible.

Note that the water molecules are not aligned into a rigid, ordered structure in practice all of the molecules are moving rapidly and randomly. Molecular dynamics simulations represent all of the intermolecular interactions with classical potentials, generating forces and acceleration via Newton s laws. Such simulations give very good descriptions of many of the properties of bulk solutions. 

Sumpter and Thompson [70] used DMNA as a prototypical nitramine in one of the earliest molecular dynamics simulations of gas-phase decompositions via competing pathways. The studies focused on practical aspects of simulating unimolecular reactions in large molecules (e.g., the influence of the details of the potential energy surface) and the fundamental dynamics (e.g., IVR) on the decomposition reactions. They carried out simulations using various models for the potential energy surfaces and for various initial energy distributions. 

One of the features that makes Equation (1) such a good starting point for our work is that it can be, in principal, exact. It is possible to show, without ever explicitly evaluating T and 37, that these crucial functions really do exist and are well defined (31). These formal definitions are rarely, if ever, useful in practical numerical calculations, but one can also work backwards from the exact dynamics x(t) (e.g., from a molecular dynamics simulation) to derive what the friction in particular must look like (32). The analysis tells us, moreover, that the exact T and are actually related to one another (28,31). The requirement that the relaxed system must be in equilibrium at some temperature T can be shown to set the magnitude and correlations of the fluctuating force ... 

---