Computational studies molecular dynamics simulations

P. A. J Hilbers and K. Esselink, Parallel computing and molecular dynamics simulations , Computer Simulations in Chemical Physics, Proc. of the NATO advanced study institute on new perspectives in computer simulations in chemical physics, 473-95, 1993. 

The protein folding problem - the ability to predict a protein fold from its sequence - is one of the major prizes in computational chemistry. Molecular dynamics simulations of solvated proteins is currently not a feasible approach to this problem. However, Duan and Kollman have shown that a 1 ps simulation on a small hydrated protein, here the 36 residue villin headpiece, is now possible using a massively parallel super computer.33 The native protein is estimated to fold in about 10-100 ps and so the simulation can only be used to study the early stages of protein folding. Nevertheless, starting from an extended structure the authors were able to observe hydrophobic collapse and secondary structure formation (helix 2 was well formed, helices 1 and 3 were partially formed and the loop connecting helices 1 and 2 was also partially... 

Figure 4.49 reprinted with permission from Pranata J and W L Jorgensen. Computational Studies on FK506 Conformational Search and Molecular Dynamics Simulations in Water. The Journal of the American Chemical Society 113 9483-9493. 1991 American Chemical Society. 

The rapid rise in computer speed over recent years has led to atom-based simulations of liquid crystals becoming an important new area of research. Molecular mechanics and Monte Carlo studies of isolated liquid crystal molecules are now routine. However, care must be taken to model properly the influence of a nematic mean field if information about molecular structure in a mesophase is required. The current state-of-the-art consists of studies of (in the order of) 100 molecules in the bulk, in contact with a surface, or in a bilayer in contact with a solvent. Current simulation times can extend to around 10 ns and are sufficient to observe the growth of mesophases from an isotropic liquid. The results from a number of studies look very promising, and a wealth of structural and dynamic data now exists for bulk phases, monolayers and bilayers. Continued development of force fields for liquid crystals will be particularly important in the next few years, and particular emphasis must be placed on the development of all-atom force fields that are able to reproduce liquid phase densities for small molecules. Without these it will be difficult to obtain accurate phase transition temperatures. It will also be necessary to extend atomistic models to several thousand molecules to remove major system size effects which are present in all current work. This will be greatly facilitated by modern parallel simulation methods that allow molecular dynamics simulations to be carried out in parallel on multi-processor systems [115]. 

Molecular dynamics simulation (MDS) is a powerful tool for the processing mechanism study of silicon surface fabrication. When a particle impacts with a solid surface, what will happen Depending on the interaction between cluster and surface, behaviors of the cluster fall into several categories including implantation [20,21], deposition [22,23], repulsion [24], and emission [25]. Owing to limitations of computer time, the cluster that can be simulated has a diameter of only a few nanometres with a small cohesive energy, which induces the cluster to fragment after collision. 

Tarek et al. [388] studied a system with some similarities to the work of Bocker et al. described earlier—a monolayer of n-tetradecyltrimethylammonium bromide. They also used explicit representations of the water molecules in a slab orientation, with the mono-layer on either side, in a molecular dynamics simulation. Their goal was to model more disordered, liquid states, so they chose two larger molecular areas, 0.45 and 0.67 nm molecule Density profiles normal to the interface were calculated and compared to neutron reflectivity data, with good agreement reported. The hydrocarbon chains were seen as highly disordered, and the diffusion was seen at both areas, with a factor of about 2.5 increase from the smaller molecular area to the larger area. They report no evidence of a tendency for the chains to aggregate into ordered islands, so perhaps this work can be seen as a realistic computer simulation depiction of a monolayer in an LE state. 

Equation (4-5) can be directly utilized in statistical mechanical Monte Carlo and molecular dynamics simulations by choosing an appropriate QM model, balancing computational efficiency and accuracy, and MM force fields for biomacromolecules and the solvent water. Our group has extensively explored various QM/MM methods using different quantum models, ranging from semiempirical methods to ab initio molecular orbital and valence bond theories to density functional theory, applied to a wide range of applications in chemistry and biology. Some of these studies have been discussed before and they are not emphasized in this article. We focus on developments that have not been often discussed. 

Sometimes the theoretical or computational approach to description of molecular structure, properties, and reactivity cannot be based on deterministic equations that can be solved by analytical or computational methods. The properties of a molecule or assembly of molecules may be known or describable only in a statistical sense. Molecules and assemblies of molecules exist in distributions of configuration, composition, momentum, and energy. Sometimes, this statistical character is best captured and studied by computer experiments molecular dynamics, Brownian dynamics, Stokesian dynamics, and Monte Carlo methods. Interaction potentials based on quantum mechanics, classical particle mechanics, continuum mechanics, or empiricism are specified and the evolution of the system is then followed in time by simulation of motions resulting from these direct... 

Molecular dynamics simulations are capable of addressing the self-assembly process at a rudimentary, but often impressive, level. These calculations can be used to study the secondary structure (and some tertiary structure) of large complex molecules. Present computers and codes can handle massive calculations but cannot eliminate concerns that boundary conditions may affect the result. Eventually, continued improvements in computer hardware will provide this added capacity in serial computers development of parallel computer codes is likely to accomplish the goal more quickly. In addition, the development of realistic, time-efficient potentials will accelerate the useful application of dynamic simulation to the self-assembly process. In addition, principles are needed to guide the selec-... 

To simulate the particle-particle collision, the hard-sphere model, which is based on the conservation law for linear momentum and angular momentum, is used. Two empirical parameters, a restitution coefficient of 0.9 and a friction coefficient of 0.3, are utilized in the simulation. In this study, collisions between spherical particles are assumed to be binary and quasi-instantaneous. The equations, which follow those of molecular dynamic simulation, are used to locate the minimum flight time of particles before any collision. Compared with the soft-sphere particle-particle collision model, the hard-sphere model accounts for the rotational particle motion in the collision dynamics calculation thus, only the translational motion equation is required to describe the fluid induced particle motion. In addition, the hard-sphere model also permits larger time steps in the calculation therefore, the simulation of a sequence of collisions can be more computationally effective. The details of this approach can be found in the literature (Hoomans et al., 1996 Crowe et al., 1998). 

The computational details for both the quantum mechanical and FEP studies using molecular dynamics simulations are described elsewhere.9 In brief, the gas-phase free energies (AGgas) were calculated using energies... 

One of the most convincing tests of the AG relationship appeared in the work of Scala et al.92 for the SPC/E model of water,57 which is known to reproduce many of water s distinctive properties in its super-cooled liquid state qualitatively. In this study, the dynamical quantity used to correlate with the configurational entropy was the self-diffusivity D. Scala et al. computed D via molecular dynamics simulations. The authors calculated the various contributions to the liquid entropy using the methods described above for a wide range of temperature and density [shown in Figure 12(a-c)]. 

None of the methods currently used to study molecular dynamics can span the whole time range of motions of interest, from picoseconds to seconds and minutes. However, the structural resolution of a method is of equal importance. A method has to not only provide information about the existence of motions with definite velocities but also to identify what structural element is moving and what is the mechanism of motion. Computer simulation of molecular dynamics has proved to be a very important tool for the development of theories concerning times and mechanisms of motions in proteins. In this approach, the initial coordinates and forces on each atom are input into the calculations, and classical equations of motions are solved by numerical means. The lengthy duration of the calculation procedure, even with powerful modem computers, does not permit the time interval investigated to be extended beyond hundreds of picoseconds. In addition, there are strong... 

---