Nonbond Interactions

The energy of interactions between nonbonded atoms is accounted for by the following  

---

Guenot, J., Kollman, P.A. Conformational and energetic effects of truncating nonbonded interactions in an aqueous protein dynamics simulation. J. Comput. Chem. 14 (1993) 295-311. 

In LN, the bonded interactions are treated by the approximate linearization, and the local nonbonded interactions, as well as the nonlocal interactions, are treated by constant extrapolation over longer intervals Atm and At, respectively). We define the integers fci,fc2 > 1 by their relation to the different timesteps as Atm — At and At = 2 Atm- This extrapolation as used in LN contrasts the modern impulse MTS methods which only add the contribution of the slow forces at the time of their evaluation. The impulse treatment makes the methods symplectic, but limits the outermost timestep due to resonance (see figures comparing LN to impulse-MTS behavior as the outer timestep is increased in [88]). In fact, the early versions of MTS methods for MD relied on extrapolation and were abandoned because of a notable energy drift. This drift is avoided by the phenomenological, stochastic terms in LN. 

Formally, we describe the LN method with the above force splitting below for the triplet protocol At, Atm, At. The fast, medium, and slow force components are distinguished by subscripts we take the medium forces as those nonbonded interactions within a 6 A region. 

In an atomic level simulation, the bond stretch vibrations are usually the fastest motions in the molecular dynamics of biomolecules, so the evolution of the stretch vibration is taken as the reference propagator with the smallest time step. The nonbonded interactions, including van der Waals and electrostatic forces, are the slowest varying interactions, and a much larger time-step may be used. The bending, torsion and hydrogen-bonding forces are treated as intermediate time-scale interactions. 

Ding H-Q, N Karasawa and W A Goddard III 1992a. Atomic Level Simulations on a Milhon Particles The Cell Multipole Method for Coulomb and London Nonbonding Interactions. Journal of Chemical Physics 97 4309-4315. 

I J, J C Cole, J P M Lommerse, R S Rowland, R Taylor and M L Verdonk 1997. Isostar A Libraij )f Information about Nonbonded Interactions. Journal of Computer-Aided Molecular Design 11 525-531. g G, W C Guida and W C Still 1989. An Internal Coordinate Monte Carlo Method for Searching lonformational Space. Journal of the American Chemical Scociety 111 4379-4386. leld C and A J Collins 1980. Introduction to Multivariate Analysis. London, Chapman Hall, ig C-W, R M Cooke, A E I Proudfoot and T N C Wells 1995. The Three-dimensional Structure of 1 ANTES. Biochemistry 34 9307-9314. 

MOMEC is a force field for describing transition metal coordination compounds. It was originally parameterized to use four valence terms, but not an electrostatic term. The metal-ligand interactions consist of a bond-stretch term only. The coordination sphere is maintained by nonbond interactions between ligands. MOMEC generally works reasonably well for octahedrally coordinated compounds. 

YETI is a force held designed for the accurate representation of nonbonded interactions. It is most often used for modeling interactions between biomolecules and small substrate molecules. It is not designed for molecular geometry optimization so researchers often optimize the molecular geometry with some other force held, such as AMBER, then use YETI to model the docking process. Recent additions to YETI are support for metals and solvent effects. 

Most of the methods proposed include a van der Waals term for describing nonbonded interactions between atoms in the two regions. This is usually represented by a Leonard-Jones 6-12 potential of the form... 

If the QM and MM regions are separate molecules, having nonbonded interactions only might be sufficient. If the two regions are parts of the same molecule, it is necessary to describe the bond connecting the two sections. In most... 

First determine what parameters will be used for describing bond lengths and angles. Then determine torsional, inversion, and nonbonded interaction parameters. 

Calculating nonbonded interactions only to a certain distance imparts an error in the calculation. If the cutoff radius is fairly large, this error will be very minimal due to the small amount of interaction at long distances. This is why many bulk-liquid simulations incorporate 1000 molecules or more. As the cutoff radius is decreased, the associated error increases. In some simulations, a long-range correction is included in order to compensate for this error. 

Nonbonded interactions are the forces be tween atoms that aren t bonded to one another they may be either attractive or repulsive It often happens that the shape of a molecule may cause two atoms to be close in space even though they are sep arated from each other by many bonds Induced dipole/induced dipole interactions make van der Waals forces in alkanes weakly attractive at most distances but when two atoms are closer to each other than the sum of their van der Waals radii nuclear-nuclear and electron-electron repulsive forces between them dominate the fvan derwaais term The resulting destabilization is called van der Waals strain... 

United atom force fields (see United versus All Atom Force Fields on page 28) are sometimes used for biomolecules to decrease the number of nonbonded interactions and the computation time. Another reason for using a simplified potential is to reduce the dimensionality of the potential energy surface. This, in turn, allows for more samples of the surface. 

HyperChem also provides a shifting potential for terminating nonbonded interactions (equation 15). 

HyperChem supplements the standard MM2 force field (see References on page 106) by providing additional parameters (force constants) using two alternative schemes (see the second part of this book. Theory and Methods). This extends the range of chemical compounds that MM-t can accommodate. MM-t also provides cutoffs for calculating nonbonded interactions and periodic boundary conditions. 

You can choose to calculate all nonbonded interactions or to truncate (cut off) the nonbonded interaction calculations using a switched or shifted function. Computing time for molecular mechanics calculations is largely a function of the number of nonbonded interactions, so truncating nonbonded interactions reduces computing time. You must also truncate nonbonded interactions for periodic boundary conditions to prevent interaction problems between nearest neighbor images. 

Eor small and medium-sized molecules, calculate all nonbonded interactions. 

Usually, atoms with a vicinal relationship or more are considered to be nonbonded. Sometimes, however, only atoms with a 1-5 relationship are considered to be fully nonbonded and the atoms with a relationship have scaled down nonbonded interactions or are deleted completely from the nonbonded computations, or different parameters are used. 

The above potential describes the monopole-monopole interactions of atomic charges Q and Qj a distance Ry apart. Normally these charge interactions are computed only for nonbonded atoms and once again the interactions might be treated differently from the more normal nonbonded interactions (1-5 relationship or more). The dielectric constant 8 used in the calculation is sometimes scaled or made distance-dependent, as described in the next section. 

---