Nucleic Acid Force Fields

Data Correlation with Chemical Structure NMR Refinement Nucleic Acid Conformation and Flexibility Modeling Using Molecular Mechanics Nucleic Acid Force Fields Nucleic Acids Qualitative Modeling Object-oriented Programming. 

Weiner, Kolltnan, I ,.A, Nguyen, D. f. Case, n,A,. All all aloin force field for simulations of proteins and nucleic acids,/. Cornpui. Chem. 7 230 252, 1986. 

The AMBER (Assisted Model Building and Energy Refin emeni) is based on a force field developed for protein and nucleic acid computations by members of the Peter Kollman research group at the... 

OPTS (Optim i/.ed Potentials for Liquid Simulations) is based on a force field developed by the research group of Bill Jorgensen now at Yale University and previously at Purdue University. Like AMBER, the OPLS force field is designed for calculations on proteins an d nucleic acids. It in troduces non bonded in leraclion parameters that have been carefully developed from extensive Monte Carlo liquid sim u lation s of small molecules. These n on-bonded interactions have been added to the bonding interactions of AMBER to produce a new force field that is expected to be better than AMBER at describing simulations w here the solvent isexplic-... 

The OPLS force field is described in twtt papers, one discussing parameters for proteins W. L. Jorgensen and J. Tirado-Rives,/. Amer. (. hem. Soc., 110, 1557 (iy8K) and on e discii ssin g param eters for n iicleotide bases [J. Pranata, S. Wiersch ke, and W. L. Jorgen sen. , /.. Amer. Chem. Soc.. 117, 281(1 ( 1991)1. The force field uses the united atom concept ftir many, but not all. hydrttgens attached to carbons to allow faster calculation s on macromolecular systems. The amino and nucleic acid residue templates in HyperChein automatically switch to a united atom representation where appropriate when th e OPLS option is selected. 

Cornell W D, P Cieplak, CI Bayly, I R Gould, K M Merz Jr, D M Ferguson, D C Spellmeyer, T Fox, J W Caldwell and P A Kollman 1995. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids and Organic Molecules. Journal of the American Chemical Society 117 5179-5197. 

Assisted model building with energy refinement (AMBER) is the name of both a force field and a molecular mechanics program. It was parameterized specifically for proteins and nucleic acids. AMBER uses only five bonding and nonbonding terms along with a sophisticated electrostatic treatment. No cross terms are included. Results are very good for proteins and nucleic acids, but can be somewhat erratic for other systems. 

The AMBER and CHARMM force fields are best suited for protein and nucleic acid studies. 

ChemSketch has some special-purpose building functions. The peptide builder creates a line structure from the protein sequence defined with the typical three-letter abbreviations. The carbohydrate builder creates a structure from a text string description of the molecule. The nucleic acid builder creates a structure from the typical one-letter abbreviations. There is a function to clean up the shape of the structure (i.e., make bond lengths equivalent). There is also a three-dimensional optimization routine, which uses a proprietary modification of the CHARMM force field. It is possible to set the molecule line drawing mode to obey the conventions of several different publishers. 

Focuses on force field calculations for understanding the dynamic properties of proteins and nucleic acids. Provides a useful introduction to several computational techniques, including molecular mechanics minimization and molecular dynamics. Includes discussions of research involving structural changes and short time scale dynamics of these biomolecules, and the influence of solvent in these processes. 

Weiner, S.J. Kollman, P.A. Case, D.A. Singh, U.C., Ohio, C. Alagona, G. Profeta Jr., S. Weiner, P. Anew force field for molecular mechanical simulation of nucleic acids and proteins 7 Am. Chem. Soc. 106 765-784, 1984. 

In computational chemistry it can be very useful to have a generic model that you can apply to any situation. Even if less accurate, such a computational tool is very useful for comparing results between molecules and certainly lowers the level of pain in using a model from one that almost always fails. The MM+ force field is meant to apply to general organic chemistry more than the other force fields of HyperChem, which really focus on proteins and nucleic acids. HyperChem includes a default scheme such that when MM+ fails to find a force constant (more generally, force field parameter), HyperChem substitutes a default value. This occurs universally with the periodic table so all conceivable molecules will allow computations. Whether or not the results of such a calculation are realistic can only be determined by close examination of the default parameters and the particular molecular situation. ... 

To date, a number of simulation studies have been performed on nucleic acids and proteins using both AMBER and CHARMM. A direct comparison of crystal simulations of bovine pancreatic trypsin inliibitor show that the two force fields behave similarly, although differences in solvent-protein interactions are evident [24]. Side-by-side tests have also been performed on a DNA duplex, showing both force fields to be in reasonable agreement with experiment although significant, and different, problems were evident in both cases [25]. It should be noted that as of the writing of this chapter revised versions of both the AMBER and CHARMM nucleic acid force fields had become available. Several simulations of membranes have been performed with the CHARMM force field for both saturated [26] and unsaturated [27] lipids. The availability of both protein and nucleic acid parameters in AMBER and CHARMM allows for protein-nucleic acid complexes to be studied with both force fields (see Chapter 20), whereas protein-lipid (see Chapter 21) and DNA-lipid simulations can also be performed with CHARMM. 

Significant progress in the optimization of VDW parameters was associated with the development of the OPLS force field [53]. In those efforts the approach of using Monte Carlo calculations on pure solvents to compute heats of vaporization and molecular volumes and then using that information to refine the VDW parameters was first developed and applied. Subsequently, developers of other force fields have used this same approach for optimization of biomolecular force fields [20,21]. Van der Waals parameters may also be optimized based on calculated heats of sublimation of crystals [68], as has been done for the optimization of some of the VDW parameters in the nucleic acid bases [18]. Alternative approaches to optimizing VDW parameters have been based primarily on the use of QM data. Quantum mechanical data contains detailed information on the electron distribution around a molecule, which, in principle, should be useful for the optimization of VDW... 

WD Cornell, P Cieplak, Cl Bayly, IR Gould, KM Merz Jr, DM Fergusson, DC Spellmeyer, DC Fox, JW Caldwell, PA Kollman. A second generation force field for the simulation of proteins and nucleic acids. J Am Chem Soc 117 5179-5197, 1995. 

The methodological advances just presented have brought the field of nucleic acid force field calculations to a point where results from the calculations can be used with reasonable confidence to aid in the interpretation of experimental data as well as to be used for scientific investigations that are not accessible to experiment. Accordingly, a number of studies based on MD simulations, as well as other methods, have been undertaken to study a wide array of biologically relevant events associated with DNA. A brief overview of some of these efforts follows. 

Central to the quality of any computational smdy is the mathematical model used to relate the structure of a system to its energy. General details of the empirical force fields used in the study of biologically relevant molecules are covered in Chapter 2, and only particular information relevant to nucleic acids is discussed in this chapter. 

The second generation force fields for nucleic acids were designed to be used with an explicit solvent representation along with inclusion of the appropriate ions [28,29]. In addition, efforts were made to improve the representation of the conformational energetics of selected model compounds. Eor example, the availability of high level ab initio calculations on the conformational energetics of the model compound dimethylphosphate yielded... 

---