Force Fields MM

The MM+ force field uses special values for ihe bending force coii-sianis when Lhe aloins are in a Ihree- or fmir-niembered ring. 

Th c values of V I, V2, and V3, in keal/tnol, are listed in mmp-tor.txtfdbf). fhc MM+ force field uses special values for the torsional force constants when the atoms are in a four-member ring. 

For XH bonds, where X isany heavy atom, the hydrogen electron den sity is ri ot th ough t to be cen tered at th e position of th e hydrogen n ueleus but displaced alon g th e bon d sorn ewhat, towards X. The MM+ force field reduces the XH bond length by a factor of 0.9 I 5 strictly for th e purposes of calculatin g van der Waals in teraction s with hydrogen atoms. 

In computational chemistry it can be very useful to have a generic model that you can apply to any situation. Even if less accurate, such a computational tool is very useful for comparing results between molecules and certainly lowers the level of pain in using a model from one that almost always fails. The MM+ force field is meant to apply to general organic chemistry more than the other force fields of HyperChem, which really focus on proteins and nucleic acids. HyperChem includes a default scheme such that when MM+ fails to find a force constant (more generally, force field parameter), HyperChem substitutes a default value. This occurs universally with the periodic table so all conceivable molecules will allow computations. Whether or not the results of such a calculation are realistic can only be determined by close examination of the default parameters and the particular molecular situation. ... 

Finally, the parametrization of the van der Waals part of the QM-MM interaction must be considered. This applies to all QM-MM implementations irrespective of the quantum method being employed. From Eq. (9) it can be seen that each quantum atom needs to have two Lennard-Jones parameters associated with it in order to have a van der Walls interaction with classical atoms. Generally, there are two approaches to this problem. The first is to derive a set of parameters, e, and G, for each common atom type and then to use this standard set for any study that requires a QM-MM study. This is the most common aproach, and the derived Lennard-Jones parameters for the quantum atoms are simply the parameters found in the MM force field for the analogous atom types. For example, a study that employed a QM-MM method implemented in the program CHARMM [48] would use the appropriate Lennard-Jones parameters of the CHARMM force field [52] for the atoms in the quantum region. 

Some authors refer to improper dihedrals when discussing dihedral torsion where we would not normally expect full rotation, for example, any of the C-C-C-C linkages in the benzene ring of aspirin. Many MM force fields treat... 

Studies of chemical reactions in solution and in enzymes present an enormous challenge because of the enormous size and complexity of these systems. MM force fields have made a tremendous impact in certain areas, but they cannot... 

It is worth noting that much of the development work for the MM force fields has centred on low energy structures of molecules. Consequently, some of the force constants are less applicable to higher energy molecular structures that can occur in molecular dynamics simulations of liquid crystals. 

Despite advent of theoretical methods and techniques and faster computers, no single theoretical method seems to be capable of reliable computational studies of reactivities of biocatalysts. Ab initio quantum mechanical (QM) methods may be accurate but are still too expensive to apply to large systems like biocatalysts. Semi-empirical quantum methods are not as accurate but are faster, but may not be fast enough for long time simulation of large molecular systems. Molecular mechanics (MM) force field methods are not usually capable of dealing with bond-breaking and formation... 

In combined QM/MM potentials, the system is divided into a QM region and an MM region. The QM region typically includes atoms that are directly involved in the chemical step and they are treated explicitly by a quantum mechanical electronic structure method. The MM region consists of the rest of the system and is approximated by an MM force field. The QM/MM potential is given by ... 

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. 

The total electric field, E, is composed of the external electric field from the permanent charges E° and the contribution from other induced dipoles. This is the basis of most polarizable force fields currently being developed for biomolecular simulations. In the present chapter an overview of the formalisms most commonly used for MM force fields will be presented. It should be emphasized that this chapter is not meant to provide a broad overview of the field but rather focuses on the formalisms of the induced dipole, classical Drude oscillator and fluctuating charge models and their development in the context of providing a practical polarization model for molecular simulations of biological macromolecules [12-21], While references to works in which the different methods have been developed and applied are included throughout the text, the major discussion of the implementation of these models focuses... 

Fig. 1 Structures of the hairpins 3G (left) and 3GAGG (right) calculated using the MM force field in Hyperchem 5.01a. The stilbenedicarboxamide located at the top of the structure and the guanine bases are enlarged to show their location... 

In the pseudobond method of Yang and coworkers [47] a pseudobond is formed with one free-valence atom with an effective core potential (optimized to reproduce the length and strength of the real bond). This core potential can be applied in Hartree-Fock and density functional calculations and is designed to be independent of the choice of the MM force field. 

The total potential energy of the system in the combined QM/MM force field is calculated from the expectation value of the wave function, over Herr... 

Once the energy of the solute/solvent system has been determined, using some variant of a classical or a QM/MM force field, it can be used to find the free energy of solvation, which is the property of real interest. This is given, at constant temperature and pressure or at constant temperature and volume, respectively, by Eqs. (20) and (21) ... 

---