Periodic boundaries molecular dynamics with

Trypsin in aqueous solution has been studied by a simulation with the conventional periodic boundary molecular dynamics method and an NVT ensemble.312 340 A total of 4785 water molecules were included to obtain a solvation shell four to five water molecules thick in the periodic box the analysis period was 20 ps after an equilibration period of 20 ps at 285 K. The diffusion coefficient for the water, averaged over all molecules, was 3.8 X 10-5 cm2/s. This value is essentially the same as that for pure water simulated with the same SPC model,341 3.6 X 10-5 cm2/s at 300 K. However, the solvent mobility was found to be strongly dependent on the distance from the protein. This is illustrated in Fig. 47, where the mean diffusion coefficient is plotted versus the distance of water molecules from the closest protein atom in the starting configuration the diffusion coefficient at the protein surface is less than half that of the bulk result. The earlier simulations of BPTI in a van der Waals solvent showed similar, though less dramatic behavior 193 i.e., the solvent molecules in the first and second solvation layers had diffusion coefficients equal to 74% and 90% of the bulk value. A corresponding reduction in solvent mobility is observed for water surrounding small biopolymers.163 Thus it... 

As electric fields and potential of molecules can be generated upon distributed p, the second order energies schemes of the SIBFA approach can be directly fueled by the density fitted coefficients. To conclude, an important asset of the GEM approach is the possibility of generating a general framework to perform Periodic Boundary Conditions (PBC) simulations. Indeed, such process can be used for second generation APMM such as SIBFA since PBC methodology has been shown to be a key issue in polarizable molecular dynamics with the efficient PBC implementation [60] of the multipole based AMOEBA force field [61]. 

B. MOLECULAR DYNAMICS WITH CONVENTIONAL PERIODIC BOUNDARY CONDITIONS... 

Molecular dynamics with periodic boundary conditions is presently the most widely used approach for studying the equilibrium and dynamic properties of pure bulk solvent,97 as well as solvated systems. However, periodic boundary conditions have their limitations. They introduce errors in the time development of equilibrium properties for times greater than that required for a sound wave to traverse the central cell. This is because the periodicity of information flow across the boundaries interferes with the time development of other processes. The velocity of sound through water at a density of 1 g/cm3 and 300 K is 15 A/ps for a cubic cell with a dimension of 45 A, the cycle time is only 3 ps and the time development of all properties beyond this time may be affected. Also, conventional periodic boundary methods are of less use for studies of chemical reactions involving enzyme and substrate molecules because there is no means for such a system to relax back to thermal equilibrium. This is not the case when alternative ensembles of the constant-temperature variety are employed. However, in these models it is not clear that the somewhat arbitrary coupling to a constant temperature heat bath does not influence the rate of reequilibration from a thermally perturbed... 

For an understanding of protein-solvent interactions it is necessary to explore the modifications of the dynamics and structure of the surrounding water induced by the presence of the biopolymer. The theoretical methods best suited for this purpose are conventional molecular dynamics with periodic boundary conditions and stochastic boundary molecular dynamics techniques, both of which treat the solvent explicitly (Chapt. IV.B and C). We focus on the results of simulations concerned with the dynamics and structure of water in the vicinity of a protein both on a global level (i.e., averages over all solvation sites) and on a local level (i.e., the solvent dynamics and structure in the neighborhood of specific protein atoms). The methods of analysis are analogous to those commonly employed in the determination of the structure and dynamics of water around small solute molecules.163 In particular, we make use of the conditional protein solute -water radial distribution function,... 

Vienna ab initio Simulation Package (VASP) (Kresse and FurethmuUer 2000), and Car-Parrinello MD Program (CAMP-Atami) (Ohnishi 1994) combine molecular dynamics with DFT under a periodic boundary condition with the orbitals expanded in the plane wave. They have been routinely used in industrial applications such as heterogeneous catalysis. 

A typical molecular dynamics simulation comprises an equflibration and a production phase. The former is necessary, as the name imphes, to ensure that the system is in equilibrium before data acquisition starts. It is useful to check the time evolution of several simulation parameters such as temperature (which is directly connected to the kinetic energy), potential energy, total energy, density (when periodic boundary conditions with constant pressure are apphed), and their root-mean-square deviations. Having these and other variables constant at the end of the equilibration phase is the prerequisite for the statistically meaningful sampling of data in the following production phase. 

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. 

The tests in the two previous paragraphs are often used because they are easy to perform. They are, however, limited due to their neglect of intermolecular interactions. Testing the effect of intennolecular interactions requires much more intensive simulations. These would be simulations of the bulk materials, which include many polymer strands and often periodic boundary conditions. Such a bulk system can then be simulated with molecular dynamics, Monte Carlo, or simulated annealing methods to examine the tendency to form crystalline phases. 

Force field calculations often truncate the non bonded potential energy of a molecular system at some finite distance. Truncation (nonbonded cutoff) saves computing resources. Also, periodic boxes and boundary conditions require it. However, this approximation is too crude for some calculations. For example, a molecular dynamic simulation with an abruptly truncated potential produces anomalous and nonphysical behavior. One symptom is that the solute (for example, a protein) cools and the solvent (water) heats rapidly. The temperatures of system components then slowly converge until the system appears to be in equilibrium, but it is not. 

In a normal molecular dynamics simulation with repeating boundary conditions (i.e., periodic boundary condition), the volume is held fixed, whereas at constant pressure the volume of the system must fluemate. In some simulation cases, such as simulations dealing with membranes, it is more advantageous to use the constant-pressure MD than the regular MD. Various schemes for prescribing the pressure of a molecular dynamics simulation have also been proposed and applied [23,24,28,29]. In all of these approaches it is inevitable that the system box must change its volume. 

The study of the enantioselective hydrosilylation reaction was performed with a series of combined quantum mechanics/molecular mechanics (QM/MM) calculations [26, 30] within the computational scheme of ab initio (AIMD) (Car-Parrinello) [62] molecular dynamics. The AIMD approach has been described in a number of excellent reviews [63-66], AIMD as well as hybrid QM/MM-AIMD calculations [26, 47] were performed with the ab initio molecular dynamics program CPMD [67] based on a pseudopotential framework, a plane wave basis set, and periodic boundary conditions. We have recently developed an interface to the CPMD package in which the coupling with a molecular mechanics force field has been implemented [26, 68],... 

---