Molecular dynamics simulations field—parameterization

A force field for solid state modeling of fluoropolymers predicted a suitable helical conformation but required further improvement in describing intermole-cular effects. Though victory cannot yet be declared, the derived force fields improve substantially on those previously available. Preliminary molecular dynamics simulations with the interim force field indicate that modeling of PTFE chain behavior can now be done in an all-inclusive manner instead of the piecemeal focus on isolated motions and defects required previously. Further refinement of the force field with a backbone dihedral term capable of reproducing the complex torsional profile of perfluorocarbons has provided a parameterization that promises both qualitative and quantitative modeling of fluoropolymer behavior in the near future. 

The atomic radii may be further refined to improve the agreement between experimental and theoretical solvation free energies. Work on this direction has been done by Luque and Orozco (see [66] and references cited therein) while Barone et al. [67] defined a set of rules to estimate atomic radii. Further discussion on this point can be found in the review by Tomasi and co-workers [15], It must be noted that the parameterization of atomic radii on the basis of a good experiment-theory agreement of solvation energies is problematic because of the difficulty to separate electrostatic and non-electrostatic terms. The comparison of continuum calculations with statistical simulations provides another way to check the validity of cavity definition. A comparison between continuum and classical Monte Carlo simulations was reported by Costa-Cabral et al. [68] in the early 1980s and more recently, molecular dynamics simulations using combined quantum mechanics and molecular mechanics (QM/MM) force-fields have been carried out to analyze the case of water molecule in liquid water [69],... 

The semiempirical methods are dramatically (100- to 1,000-fold) faster and are suitable tor many applications. The results of AMI calculations often are used as the starting points tor parameterizations (e.g., atomic charges) of force fields in molecular dynamics simulations (27,28) and CoMFA quantitative structure-activity relationship (QSAR) modeling studies (29). 

Figure 15.3 Generalized flow diagram of physics-based models, (a] Sequence information for PDB ID 4JRC is built into a (b] linear, coarsegrained molecular structure, (c] The coarse-grained structure undergoes molecular dynamics simulations with a parameterized force field, (d] The optimal structure is chosen and then (e] converted to an all-atom structure, which is minimized to form (f] the final three-dimensional structure. Disclaimer this diagram is generalized and not comprehensive to all physics-based models.

At the molecular level, refined approaches could involve explicit correlations between sidechain fluctuations and proton mobility. More recently, molecular dynamics simulations, utilizing force-field parameterizations [60, 62-64,83,84,109] or ab initio calculations based on a quantum mechanical description, were employed in order to study such correlations and exam-... 

The difference from classical force field based simulations where the forces are calculated from pre-defined pair potentials is that the forces are derived from the global potential energy surface of an electronic structure theory. The vastly higher computational costs of an electronic structure calculation restrict the system size and the length of trajectories accessible by ab initio molecular dynamics simulations. However, it becomes clear that CPMD and AIMD are important steps towards general predictive methods, due to their independence from parameterizations. 

The utility of CP dynamics is that it allows one to follow the electronic structure throughout the course of the simulation. This is not the case when empirical force fields are used. One simple property that can be evaluated throughout the CP liquid water simulation is the average dipole moment of a water molecule. The value extracted from these simulations, 2.66 D, is in good agreement with their cited experimental value of 2.6 Within CP dynamics, polarization of the water molecules (their gas phase dipole moment is 1.85 D) is automatically achieved, bypassing the parameterization process that use of a polarizable, nonadditive force field would entail. The ab initio molecular dynamics simulation also allowed them to follow the nature the lowest unoccupied molecular orbital (LUMO) which plays a role in the conduction of an excess electron. 

The choice of the adjustable parameters used in conjunction with classical potentials can result to either effective potentials that implicitly include the nuclear quantization and can therefore be used in conjunction with classical simulations (albeit only for the conditions they were parameterized for) or transferable ones that attempt to best approximate the Born-Oppenheimer PES and should be used in nuclear quantum statistical simulations. Representative examples of effective force fields for water consist of TIP4P (Jorgensen et al. 1983), SPC/E (Berendsen et al. 1987) (pair-wise additive), and Dang-Chang (DC) (Dang and Chang 1997) (polarizable, many-body). The polarizable potentials contain - in addition to the pairwise additive term - a classical induction (polarization) term that explicitly (albeit approximately) accounts for many-body effects to infinite order. These effective potentials are fitted to reproduce bulk-phase experimental data (i.e., the enthalpy at T = 298 K and the radial distribution functions at ambient conditions) in classical molecular dynamics simulations of liquid water. Despite their simplicity, these models describe some experimental properties of liquid... 

Aqueous Interfaces Biomembranes Modeling CHAR-MM The Energy Function and Its Parameterization Hydrophobic Effect Molecular Dynamics Studies of Lipid Bilayers Molecular Dynamics Techniques and Applications to Proteins Permeation of Lipid Membranes Molecular Dynamics Simulations Protein Force Fields. 

In light of the differences between a Morse and a harmonic potential, why do force fields use the harmonic potential First, the harmonic potential is faster to compute and easier to parameterize than the Morse function. The two functions are similar at the potential minimum, so they provide similar values for equilibrium structures. As computer resources expand and as simulations of thermal motion (See Molecular Dynamics , page 69) become more popular, the Morse function may be used more often. 

---