Non-bonding

Wiliams D E 1965 Non-bonded potential parameters derived from crystalline aromatic hydrocarbons J. Chem. Phys. 45 3770... 

Atomistically detailed models account for all atoms. The force field contains additive contributions specified in tenns of bond lengtlis, bond angles, torsional angles and possible crosstenns. It also includes non-bonded contributions as tire sum of van der Waals interactions, often described by Lennard-Jones potentials, and Coulomb interactions. Atomistic simulations are successfully used to predict tire transport properties of small molecules in glassy polymers, to calculate elastic moduli and to study plastic defonnation and local motion in quasi-static simulations [fy7, ( ]. The atomistic models are also useful to interiDret scattering data [fyl] and NMR measurements [70] in tenns of local order. 

Figure C2.5.7. The native confonnation of fast folding sequence (N= 27) witli random bond potentials is shown. This stmcture has c = 22 non-bonded contacts, tlierefore it is not a maximally compact confonnation for which c = 28. The figure was created using RasMol 2.6 [8].

The very low bond dissociation enthalpy of fluorine is an important factor contributing to the greater reactivity of fluorine. (This low energy may be due to repulsion between non-bonding electrons on the two adjacent fluorine atoms.) The higher hydration and lattice enthalpies of the fluoride ion are due to the smaller size of this ion. 

Ligand-Protein Interactions The energy of formation of ligand-protein contacts can be computed as a sum of non-bonded (Lennard-Jones and electrostatic) terms similar to those used in a molecular dynamics simulation. 

In addition, the non-bonded forces can be divided into several regions according to pair distances. The near region is normally more important than the distant region because the non-bonded forces decay with distance. Since most of the CPU time in a MD simulation is spent in the calculation of these non-bonded interactions, the separation in pair distance results in valuable speedups. Using a 3-fold distance split, the non-bonded forces are separated in 3 regions near, medium, and fax distance zones. Thus, the Liouville operator can be express as a sum of five terms... 

To separate the non-bonded forces into near, medium, and far zones, pair distance separations are used for the van der Waals forces, and box separations are used for the electrostatic forces in the Fast Multipole Method,[24] since the box separation is a more convenient breakup in the Fast Multipole Method (FMM). Using these subdivisions of the force, the propagator can be factorized according to the different intrinsic time scales of the various components of the force. This approach can be used for other complex systems involving long range forces. 

C- cutoff radius method - all non-bonded forces between particles within angstroms of each other are computed explicitly... 

The representation of non-bonding orbitals on an atom again uses the concept of. T-systems, though they may have any kind of hybridization (p, sp etc.), In Figure 2-56 the three possibilities arc shown lone pairs, radicals, and orbitals without electrons can be accommodated by this eoneept. 

Chiral carbon atoms are common, but they are not the only possible centers of chirality. Other possible chiral tetravalent atoms are Si, Ge, Sn, N, S, and P, while potential trivalent chiral atoms, in which non-bonding electrons occupy the position of the fourth ligand, are N, P, As, Sb, S, Se, and Te. Furthermore, a center of chirality does not even have to be an atom, as shown in the structure represented in Figure 2-70b, where the center of chirality is at the center of the achiral skeleton of adamantane. 

Z-matriccs arc commonly used as input to quantum mechanical ab initio and serai-empirical) calculations as they properly describe the spatial arrangement of the atoms of a molecule. Note that there is no explicit information on the connectivity present in the Z-matrix, as there is, c.g., in a connection table, but quantum mechanics derives the bonding and non-bonding intramolecular interactions from the molecular electronic wavefunction, starting from atomic wavefiinctions and a crude 3D structure. In contrast to that, most of the molecular mechanics packages require the initial molecular geometry as 3D Cartesian coordinates plus the connection table, as they have to assign appropriate force constants and potentials to each atom and each bond in order to relax and optimi-/e the molecular structure. Furthermore, Cartesian coordinates are preferable to internal coordinates if the spatial situations of ensembles of different molecules have to be compared. Of course, both representations are interconvertible. 

In the next step, the suggested models are translated into 3D space by subsequently combining the templates. Again, each model is assessed and ranked according to various structural criteria, such as the geometric fit of the 3D templates and non-bonding interactions (steric clashes). If none of the solu-... 

Figure 2-99. Elimination of non-bonded interactions (close contacts).

The space filling model developed by Corey, Pauling, and Koltun is also known as the CPK model, or scale model [197], It shows the relative volume (size) of different elements or of different parts of a molecule (Figure 2-123d). The model is based on spheres that represent the "electron cloud . These atomic spheres can be determined from the van der Waals radii (see Section 2.10.1), which indicate the most stable distance between two atoms (non-bonded nuclei). Since the spheres are all drawn to the same scale, the relative size of the overlapping electron clouds of the atoms becomes evident. The connectivities between atoms, the bonds, are not visualized because they are located beneath the atom spheres and are not visible in a non-transparent display (see Section 2.10). In contrast to other models, the CPK model makes it possible to visualize a first impression of the extent of a molecule. 

Figure 7-8. Bonded (upper row) and non-bonded (lower row) contributions to a typioal molecular mechanics force field potential energy function. The latter two types of Interactions can also occur within the same molecule.

The PEF is a sum of many individual contributions, Tt can be divided into bonded (bonds, angles, and torsions) and non-bonded (electrostatic and van der Waals) contributions V, responsible for intramolecular and, in tlic case of more than one molecule, also intermoleculai interactions. Figure 7-8 shows schematically these types of interactions between atoms, which arc included in almost all force field implementations. 

Vn is often called the barrier of rotation. This is intuitive but misleading, because the exact energetic barrier of a particular rotation is the sum of all V components and other non-bonding interactions with the atoms under consideration. The multiplicity n gives the number of minima of the function during a 360° rotation of the dihedral angle o). The phase y defines the exact position of the minima. 

It is noteworthy that it is not obligatory to use a torsional potential within a PEF. Depending on the parameterization, it is also possible to represent the torsional barrier by non-bonding interactions between the atoms separated by three bonds. In fact, torsional potentials and non-bonding 1,4-interactions are in a close relationship. This is one reason why force fields like AMBER downscale the 1,4-non-bonded Coulomb and van der Waals interactions. 

N is the number of point charges within the molecule and Sq is the dielectric permittivity of the vacuum. This form is used especially in force fields like AMBER and CHARMM for proteins. As already mentioned, Coulombic 1,4-non-bonded interactions interfere with 1,4-torsional potentials and are therefore scaled (e.g., by 1 1.2 in AMBER). Please be aware that Coulombic interactions, unlike the bonded contributions to the PEF presented above, are not limited to a single molecule. If the system under consideration contains more than one molecule (like a peptide in a box of water), non-bonded interactions have to be calculated between the molecules, too. This principle also holds for the non-bonded van der Waals interactions, which are discussed in Section 7.2.3.6. 

The raie gas atoms reveal through their deviation from ideal gas behavior that electrostatics alone cannot account for all non-bonded interactions, because all multipole moments are zero. Therefore, no dipole-dipole or dipole-induced dipole interactions are possible. Van der Waals first described the forces that give rise to such deviations from the expected behavior. This type of interaction between two atoms can be formulated by a Lennaid-Jones [12-6] function Eq. (27)). 

A second idea to save computational time addresses the fact that hydrogen atoms, when involved in a chemical bond, show the fastest motions in a molecule. If they have to be reproduced by the simulation, the necessary integration time step At has to be at least 1 fs or even less. This is a problem especially for calculations including explicit solvent molecules, because in the case of water they do not only increase the number of non-bonded interactions, they also increase the number of fast-moving hydrogen atoms. This particular situation is taken into account... 

Figure 7-14. All-atom and united-atom representation of the amino acid isoleucine. In this example, 13 atoms, which are able to form explicit non-bonding interactions, are reduced to only four pseudo-atoms,...

---