Force fields interactions

Another area of rapid growth for particle separation has been that of Field-Flow Fractionation (FFF) originally developed by Giddings (12,13>1 1 ) (see also papers in this symposium series). Like HDC, the separation in field-flow fractionation (FFF) results from the combination of force field interactions and the convected motion of the particles, rather than a partitioning between phases. In FFF the force field is applied externally while in HDC it results from internal, interactions. 

Figure 3 Proper Torsions 5 H-C-C-H 5 H-C-C-C Force field interactions...

This book is concerned mainly with the study of the viscoelastic response of isotropic macromolecular systems to mechanical force fields. Owing to diverse influences on the viscoelastic behavior in multiphase systems (e.g., changes in morphology and interfaces by action of the force fields, interactions between phases, etc.), it is difficult to relate the measured rheological functions to the intrinsic physical properties of the systems and, as a result, the viscoelastic behavior of polymer blends and liquid crystals is not addressed in this book. 

In the past, force fields were parameterized based only on experimental data nowadays, most modem force fields include substantial quantum chemical information. According to the nature of the data used for parameterization, force fields can be classified as ab initio, semi-empirical, and empirical. Simple potentials, e.g., for argon, which require few parameters, can be fitted exclusively to macroscopic experimental data however, more complex force fields have numerous parameters and thus depend heavily on ab initio data. This chapter gives an introduction to the present state-of-the-art in this field. Attention is given to the way modeling and simulation on the scale of molecular force fields interact with other scales, which is mainly by parameter inheritance. Parameters are determined both bottom-up from quantum chemistry and top-down from experimental data. 

In many force fields, interactions between bonded atoms as well as 1 - 3 interactions are excluded, and the energy cannot be cast in the required form (1). This complication can be easily handled by considering the proper correction terms explicitly. 

The idea that unsymmetrical molecules will orient at an interface is now so well accepted that it hardly needs to be argued, but it is of interest to outline some of the history of the concept. Hardy [74] and Harkins [75] devoted a good deal of attention to the idea of force fields around molecules, more or less intense depending on the polarity and specific details of the structure. Orientation was treated in terms of a principle of least abrupt change in force fields, that is, that molecules should be oriented at an interface so as to provide the most gradual transition from one phase to the other. If we read interaction energy instead of force field, the principle could be reworded on the very reasonable basis that molecules will be oriented so that their mutual interaction energy will be a maximum. 

Kramers solution of the barrier crossing problem [45] is discussed at length in chapter A3.8 dealing with condensed-phase reaction dynamics. As the starting point to derive its simplest version one may use the Langevin equation, a stochastic differential equation for the time evolution of a slow variable, the reaction coordinate r, subject to a rapidly statistically fluctuating force F caused by microscopic solute-solvent interactions under the influence of an external force field generated by the PES F for the reaction... 

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. 

But the methods have not really changed. The Verlet algorithm to solve Newton s equations, introduced by Verlet in 1967 [7], and it s variants are still the most popular algorithms today, possibly because they are time-reversible and symplectic, but surely because they are simple. The force field description was then, and still is, a combination of Lennard-Jones and Coulombic terms, with (mostly) harmonic bonds and periodic dihedrals. Modern extensions have added many more parameters but only modestly more reliability. The now almost universal use of constraints for bonds (and sometimes bond angles) was already introduced in 1977 [8]. That polarisability would be necessary was realized then [9], but it is still not routinely implemented today. Long-range interactions are still troublesome, but the methods that now become popular date back to Ewald in 1921 [10] and Hockney and Eastwood in 1981 [11]. 

This situation, despite the fact that reliability is increasing, is very undesirable. A considerable effort will be needed to revise the shape of the potential functions such that transferability is greatly enhanced and the number of atom types can be reduced. After all, there is only one type of carbon it has mass 12 and charge 6 and that is all that matters. What is obviously most needed is to incorporate essential many-body interactions in a proper way. In all present non-polarisable force fields many-body interactions are incorporated in an average way into pair-additive terms. In general, errors in one term are compensated by parameter adjustments in other terms, and the resulting force field is only valid for a limited range of environments. 

A number of issues need to be addressed before this method will become a routine tool applicable to problems as the conformational equilibrium of protein kinase. E.g. the accuracy of the force field, especially the combination of Poisson-Boltzmann forces and molecular mechanics force field, remains to be assessed. The energy surface for the opening of the two kinase domains in Pig. 2 indicates that intramolecular noncovalent energies are overestimated compared to the interaction with solvent. 

The problems that occur when one tries to estimate affinity in terms of component terms do not arise when perturbation methods are used with simulations in order to compute potentials of mean force or free energies for molecular transformations simulations use a simple physical force field and thereby implicitly include all component terms discussed earlier. We have used the molecular transformation approach to compute binding affinities from these first principles [14]. The basic approach had been introduced in early work, in which we studied the affinity of xenon for myoglobin [11]. The procedure was to gradually decrease the interactions between xenon atom and protein, and compute the free energy change by standard perturbation methods, cf. (10). An (issential component is to impose a restraint on the... 

MTS methods exploit the existence of different time scales arising from the many interactions present in the force field. For expository purposes assume... 

In order to represent 3D molecular models it is necessary to supply structure files with 3D information (e.g., pdb, xyz, df, mol, etc.. If structures from a structure editor are used directly, the files do not normally include 3D data. Indusion of such data can be achieved only via 3D structure generators, force-field calculations, etc. 3D structures can then be represented in various display modes, e.g., wire frame, balls and sticks, space-filling (see Section 2.11). Proteins are visualized by various representations of helices, / -strains, or tertiary structures. An additional feature is the ability to color the atoms according to subunits, temperature, or chain types. During all such operations the molecule can be interactively moved, rotated, or zoomed by the user. 

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. 

For each pair of interacting atoms (/r is their reduced mass), three parameters are needed D, (depth of the potential energy minimum, k (force constant of the par-tictilar bond), and l(, (reference bond length). The Morse ftinction will correctly allow the bond to dissociate, but has the disadvantage that it is computationally very expensive. Moreover, force fields arc normally not parameterized to handle bond dissociation. To circumvent these disadvantages, the Morse function is replaced by a simple harmonic potential, which describes bond stretching by Hooke s law (Eq. (20)). 

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. 

All of the contributions to the energy function presented above assume that pairwise interactions are sufficient to describe the situation within a molecule or molecular system. Whether or not multi-centered interactions are negligible is controversial. On the other hand, failure or success of a force field with its functional form and corresponding parameter set is not a matter of mathematics... 

The Universal Force Field, UFF, is one of the so-called whole periodic table force fields. It was developed by A. Rappe, W Goddard III, and others. It is a set of simple functional forms and parameters used to model the structure, movement, and interaction of molecules containing any combination of elements in the periodic table. The parameters are defined empirically or by combining atomic parameters based on certain rules. Force constants and geometry parameters depend on hybridization considerations rather than individual values for every combination of atoms in a bond, angle, or dihedral. The equilibrium bond lengths were derived from a combination of atomic radii. The parameters [22, 23], including metal ions [24], were published in several papers. 

Many problems in force field investigations arise from the calculation of Coulomb interactions with fixed charges, thereby neglecting possible mutual polarization. With that obvious drawback in mind, Ulrich Sternberg developed the COSMOS (Computer Simulation of Molecular Structures) force field [30], which extends a classical molecular mechanics force field by serai-empirical charge calculation based on bond polarization theory [31, 32]. This approach has the advantage that the atomic charges depend on the three-dimensional structure of the molecule. Parts of the functional form of COSMOS were taken from the PIMM force field of Lindner et al., which combines self-consistent field theory for r-orbitals ( nr-SCF) with molecular mechanics [33, 34]. 

---