Molecular Mechanics Interaction

Arestraint (not to be confused with a Model Builder constraint) is a user-specified one-atom tether, two-atom stretch, three-atom bend, or four-atom torsional interaction to add to the list of molecular mechanics interactions computed for a molecule. These added interactions are treated no differently from any other stretch, bend, or torsion, except that they employ a quadratic functional form. They replace no interaction, only add to the computed interactions. 

Models used in the zeolite science can be divided into two categories (i) models that do not explicitly consider any electron in the system (molecular mechanics, interactions described with interatomic potential functions) and (ii) models that explicitly consider part of the electrons in the system (either at semiempirical level or at ab initio level). This text should serve as an introductory overview of quantum chemical approaches (excluding semiempirical methods) and models available for zeolite modeling. It is impossible to review the quantum chemical calculations in zeolite science on pages available here. Only qualitative description of methods will be given, avoiding mathematical equations. More details can be found, e. g., in Refs. [1-5]. 

All the above, and the fact that by definition the molecular mechanics interaction energy is proportional to the thermodynamic work of adhesion, hence kEtot = Wa, also brings the interesting consideration that... 

Molecular mechanics interaction energies of complexes of HIV protease inhibitors also showed good correlations between binding energy and inhibitory potency as well as good predictive ability. The Merck workers report that the interaction energies calculated with the Merck Molecular Force Field gave a correlation superior to that calculated via a commercial force field. 

The effective Hamiltonian can be divided into three terms (equation 2), which are generated from the interactions which occur within and between the components of the system (see Scheme 1). The contributions considered here include the completely quantum mechanical, Hqm (QM in Scheme 1), the purely molecular mechanical interactions, //mm (MM in Scheme 1), and the interactions between the QM and MM portions of the system, //qm/mm (indicated by the arrow in Scheme 1). 

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. 

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. 

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.

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]. 

Electrostatic terms other than the simple charge interactions above are commonly included in molecular mechanics calculations. particularly dipole-dipole interactions. More recently, second-order electrostatic interactions like those describing polarizability have been added to some force fields. 

IlyperChem avoids th e discon tin nily an d, in isotropy problem of th e implied cutoff by iin posing a sin oothed spherical cn toff within the implied cutoff. When a system is placed in a periodic box, a switched cnLoITis aiitoinatically added. The default outer radius, where the interaction is completely turned off, is the smallest of 1/2 R., 1/2 R.. and 1/2 R, so that the cutoff avoids discontinuities and is isotropic, fh is cutoff may be turned off or modified in the. Molecular Mechanics Options dialog box after solvation and before calcii lation. ... 

Halgren T A 1992. Representation of van der Waals (vdW) Interactions in Molecular Mechanics Force Fields Potential Form, Combination Rules, and vdW Parameters. Journal of the American Chemical Society 114 7827-7843. 

J G 1994. Extended Electron Distributions Applied to the Molecular Mechanics of Some termolecular Interactions. Journal of Computer-Aided Molecular Design 8 653-668. el A and M Karplus 1972. Calculation of Ground and Excited State Potential Surfaces of anjugated Molecules. 1. Formulation and Parameterisation. Journal of the American Chemical Society 1 5612-5622. 

---