Force field methods, molecular modeling

With this information in hand, it seemed reasonable to attempt to use force field methods to model the transition states of more complex, chiral systems. To that end, transition state.s for the delivery of hydrogen atom from stannanes 69 71 derived from cholic acid to the 2.2,.3-trimethy 1-3-pentyl radical 72 (which was chosen as the prototypical prochiral alkyl radical) were modeled in a similar manner to that published for intramolecular free-radical addition reactions (Beckwith-Schicsscr model) and that for intramolecular homolytic substitution at selenium [32]. The array of reacting centers in each transition state 73 75 was fixed at the geometry of the transition state determined by ah initio (MP2/DZP) molecular orbital calculations for the attack of methyl radical at trimethyltin hydride (viz. rsn-n = 1 Si A rc-H = i -69 A 6 sn-H-C = 180°) [33]. The remainder of each structure 73-75 was optimized using molecular mechanics (MM2) in the usual way. In all, three transition state conformations were considered for each mode of attack (re or ) in structures 73-75 (Scheme 14). In general, the force field method described overestimates experimentally determined enantioseleclivities (Scheme 15), and the development of a flexible model is now being considered [33]. 

The range of systems that have been studied by force field methods is extremely varied. Some force fields liave been developed to study just one atomic or molecular sp>ecies under a wider range of conditions. For example, the chlorine model of Rodger, Stone and TUdesley [Rodger et al 1988] can be used to study the solid, liquid and gaseous phases. This is an anisotropic site model, in which the interaction between a pair of sites on two molecules dep>ends not only upon the separation between the sites (as in an isotropic model such as the Lennard-Jones model) but also upon the orientation of the site-site vector with resp>ect to the bond vectors of the two molecules. The model includes an electrostatic component which contciins dipwle-dipole, dipole-quadrupole and quadrupole-quadrupole terms, and the van der Waals contribution is modelled using a Buckingham-like function. 

Molecula.rMecha.nics. Molecular mechanics (MM), or empirical force field methods (EFF), ate so called because they are a model based on equations from Newtonian mechanics. This model assumes that atoms are hard spheres attached by networks of springs, with discrete force constants. 

Theoretical models include those based on classical (Newtonian) mechanical methods—force field methods known as molecular mechanical methods. These include MM2, MM3, Amber, Sybyl, UFF, and others described in the following paragraphs. These methods are based on Hook s law describing the parabolic potential for the stretching of a chemical bond, van der Waal s interactions, electrostatics, and other forces described more fully below. The combination assembled into the force field is parameterized based on fitting to experimental data. One can treat 1500-2500 atom systems by molecular mechanical methods. Only this method is treated in detail in this text. Other theoretical models are based on quantum mechanical methods. These include ... 

The various types of successful approaches can be classified into two groups empirical model calculations based on molecular force fields and quantum mechanical approximations. In the first class of methods experimental data are used to evaluate the parameters which appear in the model. The shape of the potential surfaces in turn is described by expressions which were found to be appropriate by semiclassicala> or quantum mechanical methods. Most calculations of this type are based upon the electrostatic model. Another more general approach, the "consistent force field method, was recently applied to the forces in hydrogen-bonded crystals 48> 49>. 

Modeling studies are most useful when experimentally detenriined structures are of modest quality and suitable force-field potentials and modeling software are available. Although statistical methods such as Monte Carlo and molecular dynamics would be preferred in solution or other disordered states, we feel that energy minimization criteria are valid for static, ordered structures such as crystals. 

Dynamic NMR gives information on the number and symmetries of conformations present in solution and on the energy barriers separating these conformations. This is particularly true for systems with barriers between about 25 and 90 kJ mol-1, a situation which often occurs in the medium ring. The interpretation of the NMR data can be carried out by the examination of molecular models, but this is a relatively crude and sometimes misleading method. Empirical force field (or molecular mechanics) calculations are much superior, even though the parametrization of heteroatoms may be open to question. Quantum mechanical calculations are not very suitable the semiempirical type, e.g. MINDO, do not reproduce conformational properties of even cyclohexane satisfactorily, and the ab initio... 

Molecular Mechanics. Molecular mechanics (MM), or empirical force field methods (EFF), are so called because they are a model based on equations from Newtonian mechanics. This model assumes that atoms are hard spheres attached by networks of springs, with discrete force constants. The force constants in the equations are adjusted empirically to repro duce experimental observations. The net result is a model which relates the "mechanical" forces within a structure to its properties. Force fields are made up of sets of equations each of which represents an element of the decomposition of the total energy of a system (not a quantum mechanical eneigy, but a classical mechanical one). The sum of the components is called the force field eneigy, or steric energy, which also routinely includes the electrostatic eneigy components. Typically, the steric energy is expressed as... 

Molecular mechanical methods that use force field parameters to model and assign secondary structures (Pittsyn and Finkelstein, 1983). 

Empirical force fields for the modelling of molecular shape are based on parameters measured in solid crystals, and not in isolated molecules. Although the method may therefore appear to simulate the shape of a single molecule, the reference structure is as observed in a crystal, surrounded by, and interacting with many neighbours. Force-field parameters are conditioned by the crystal environment and any attempt to derive them from first principles should take the influence of the environment into account. 

Ken pioneered the modeling of transition states with force field methods. Before modem tools existed for locating transition structures in all but the simplest reactions, his group used ab initio calculations to find the geometries of transition states and to determine force constants for distortions away from these preferred geometries. These force constants could then be used in standard molecular mechanics calculations, in order to predict how steric effects would affect the geometries and energies of the transition structures when substituent were present. 

---