Force fields conformational energy

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. 

Conformational Adjustments The conformations of protein and ligand in the free state may differ from those in the complex. The conformation in the complex may be different from the most stable conformation in solution, and/or a broader range of conformations may be sampled in solution than in the complex. In the former case, the required adjustment raises the energy, in the latter it lowers the entropy in either case this effect favors the dissociated state (although exceptional instances in which the flexibility increases as a result of complex formation seem possible). With current models based on two-body potentials (but not with force fields based on polarizable atoms, currently under development), separate intra-molecular energies of protein and ligand in the complex are, in fact, definable. However, it is impossible to assign separate entropies to the two parts of the complex. 

Variations of these releases are implemented in almost every commercial or academic software package, which cannot be fisted in this context. A comprehensive comparison of several force fields focusing the calculation of conformational energies of organic molecules has been published by Pettersson and liljefors [1]. 

I. Pettersson, T. Liljefors, Molecular mechanics calculated conformational energies of organic molecules a comparison of force fields, in Reviews in Computational Chemistry, Vbl. 9,... 

Fig. 4.S Variation in torsional energy (AMBER force field) with O-C-C-0 torsion angle (to) for OCH -CHjO fragment. The minimum energy conformations arise for to = 60° and 300°.

The most ambitious approaches to the protein folding problem attempt to solve it from firs principles (ab initio). As such, the problem is to explore the coirformational space of th molecule in order to identify the most appropriate structure. The total number of possibl conformations is invariably very large and so it is usual to try to find only the very lowes energy structure(s). Some form of empirical force field is usually used, often augmente with a solvation term (see Section 11.12). The global minimum in the energy function i assumed to correspond to the naturally occurring structure of the molecule. 

The consistent force field (CFF) was developed to yield consistent accuracy of results for conformations, vibrational spectra, strain energy, and vibrational enthalpy of proteins. There are several variations on this, such as the Ure-Bradley version (UBCFF), a valence version (CVFF), and Lynghy CFF. The quantum mechanically parameterized force field (QMFF) was parameterized from ah initio results. CFF93 is a rescaling of QMFF to reproduce experimental results. These force fields use five to six valence terms, one of which is an electrostatic term, and four to six cross terms. 

Empirical conformational energy program for peptides (ECEPP) is the name of both a computer program and the force field implemented in that program. This is one of the earlier peptide force fields that has seen less use with the introduction of improved methods. It uses three valence terms that are fixed, a van der Waals term, and an electrostatic term. 

Ah initio calculations of polymer properties are either simulations of oligomers or band-structure calculations. Properties often computed with ah initio methods are conformational energies, polarizability, hyperpolarizability, optical properties, dielectric properties, and charge distributions. Ah initio calculations are also used as a spot check to verify the accuracy of molecular mechanics methods for the polymer of interest. Such calculations are used to parameterize molecular mechanics force fields when existing methods are insulficient, which does not happen too often. 

The molecular mechanics force fields available include MM+, OPLS, BIO+, and AMBER. Parameters missing from the force field will be automatically estimated. The user has some control over cutoff distances for various terms in the energy expression. Solvent molecules can be included along with periodic boundary conditions. The molecular mechanics calculations tested ran without difficulties. Biomolecule computational abilities are aided by functions for superimposing molecules, conformation searching, and QSAR descriptor calculation. 

PC Model has some features that are not found in many other molecular mechanics programs. This is one of the few programs that outputs the energy given by the force field and the heat of formation and a strain energy. Atom types for describing transition structures in the MMX force field are included. There is a metal coordination option for setting up calculations with metal atoms. There are also molecular similarity and conformation search functions. 

---