We describe here a general Amber force field (GAFF) for organic molecules. GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids, and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens. It uses a simple functional form and a limited number of atom types, but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges. The performance of GAFF in test cases is encouraging. In test I, 74 crystallographic structures were compared to GAFF minimized structures, with a root-mean-square displacement of 0.26 A, which is comparable to that of the Tripos 5.2 force field (0.25 A) and better than those of MMFF 94 and CHARMm (0.47 and 0.44 A, respectively). In test II, gas phase minimizations were performed on 22 nucleic acid base pairs, and the minimized structures and intermolecular energies were compared to MP2/6-31G* results. The RMS of displacements and relative energies were 0.25 A and 1.2 kcal/mol, respectively. These data are comparable to results from Parm99/RESP (0.16 A and 1.18 kcal/mol, respectively), which were parameterized to these base pairs. Test III looked at the relative energies of 71 conformational pairs that were used in development of the Parm99 force field. The RMS error in relative energies (compared to experiment) is about 0.5 kcal/mol. GAFF can be applied to wide range of molecules in an automatic fashion, making it suitable for rational drug design and database searching.

/pdf/development-and-testing-of-a-general-amber-force-field-28v6k2tn5g.pdf

Development and testing of a general amber force field.

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent−solvent, solvent−solute, and solute−solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volume...

/pdf/all-atom-empirical-potential-for-molecular-modeling-and-2nds0dxw12.pdf

All-atom empirical potential for molecular modeling and dynamics studies of proteins.

A general all-atom force field for atomistic simulation of common organic molecules, inorganic small molecules, and polymers was developed using state-of-the-art ab initio and empirical parametrization techniques. The valence parameters and atomic partial charges were derived by fitting to ab initio data, and the van der Waals (vdW) parameters were derived by conducting MD simulations of molecular liquids and fitting the simulated cohesive energies and equilibrium densities to experimental data. The combined parametrization procedure significantly improves the quality of a general force field. Validation studies based on large number of isolated molecules, molecular liquids and molecular crystals, representing 28 molecular classes, show that the present force field enables accurate and simultaneous prediction of structural, conformational, vibrational, and thermophysical properties for a broad range of molecules in isolation and in condensed phases. Detailed results of the parametrization and validation f...

COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase ApplicationsOverview with Details on Alkane and Benzene Compounds

The widely used CHARMM additive all-atom force field includes parameters for proteins, nucleic acids, lipids, and carbohydrates. In the present article, an extension of the CHARMM force field to drug-like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present article in the context of the model systems, pyrrolidine, and 3-phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform "all-CHARMM" simulations on drug-target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems.

/pdf/charmm-general-force-field-a-force-field-for-drug-like-4mhn8ejv47.pdf

CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields.

This article introduces MMFF94, the initial published version of the Merck molecular force field (MMFF). It describes the objectives set for MMFF, the form it takes, and the range of systems to which it applies. This study also outlines the methodology employed in parameterizing MMFF94 and summarizes its performance in reproducing computational and experimental data. Though similar to MM3 in some respects, MMFF94 differs in ways intended to facilitate application to condensed-phase processes in molecular-dynamics simulations. Indeed, MMFF94 seeks to achieve MM3-like accuracy for small molecules in a combined “organic/protein” force field that is equally applicable to proteins and other systems of biological significance. A second distinguishing feature is that the core portion of MMFF94 has primarily been derived from high-quality computational data—ca. 500 molecular structures optimized at the HF/6-31G* level, 475 structures optimized at the MP2/6-31G* level, 380 MP2/6-31G* structures evaluated at a defined approximation to the MP4SDQ/TZP level, and 1450 structures partly derived from MP2/6-31G* geometries and evaluated at the MP2/TZP level. A third distinguishing feature is that MMFF94 has been parameterized for a wide variety of chemical systems of interest to organic and medicial chemists, including many that feature frequently occurring combinations of functional groups for which little, if any, useful experimental data are available. The methodology used in parameterizing MMFF94 represents a fourth distinguishing feature. Rather than using the common “functional group” approach, nearly all MMFF parameters have been determined in a mutually consistent fashion from the full set of available computational data. MMFF94 reproduces the computational data used in its parameterization very well. In addition, MMFF94 reproduces experimental bond lengths (0.014 A root mean square [rms]), bond angles (1.2° rms), vibrational frequencies (61 cm−1 rms), conformational energies (0.38 kcal/mol/rms), and rotational barriers (0.39 kcal/mol rms) very nearly as well as does MM3 for comparable systems. MMFF94 also describes intermolecular interactions in hydrogen-bonded systems in a way that closely parallels that given by the highly regarded OPLS force field. © 1996 John Wiley & Sons, Inc.

Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94

A new force field has been developed for alkanes and cycloalkanes, excluding small rings, to improve the calculation of vibrational frequencies, rotational barriers, and numerous relatively small errors that were observed to result from the use of the MM3 force field. © 1996 by John Wiley & Sons, Inc.

An improved force field (MM4) for saturated hydrocarbons

The potential functions for simple amides, several peptides and a small protein have been worked out for the MM3 force field. Structures and energies were fit as previously with MM2, but additionally, we fit the vibrational spectra of the simple amides (average rms error over four compounds, 34 cm−1), and examined more carefully electrostatic interactions, including charge-charge and charge-dipole interactions. The parameters were obtained and tested by examining four simple amides, five electrostatic model complexes, two dipeptides, six crystalline cyclic peptides, and the protein Crambin. The average root-mean-square deviation from the X-ray structures for the six cyclic peptide crystals was only 0.10 A for the nonhydrogen atomic positions, and 0.011 A, 1.0°, and 4.9° for bond lengths, bond angles, and torsional angles, respectively. The parameter set was then further tested by minimizing the high resolution crystal structure of the hydrophobic protein Crambin. The resultant root-mean-square deviations for the non-hydrogen atomic data, in the presence of the crystal lattice, are 0.22 A, 0.023 A, 2.0°, and 6.4° for coordinates, bond lengths, bond angles, and torsional angles, respectively.

The MM3 force field for amides, polypeptides and proteins

An extensive quantum mechanical study of a water dimer suggests that the introduction of a diffuse function into the basis set, which significantly reduces the basis set superposition error (BSSE) in the hydrogen bonding energy calculation, is the key to better calculations of the potential energy surfaces of carbohydrates. This article examines the potential energy surfaces of selected d-aldo- and d-ketohexoses (a total of 82 conformers) by quantum mechanics (QM) and molecular mechanics (MM) methods. In contrast to the results with a smaller basis set (B3LYP/6-31G** 5d), we found at the higher level calculation (B3LYP/6-311++G(2d,2p)//B3LYP/6-31G** 5d) that, in most cases, the furanose forms are less stable than the pyranose forms. These discrepancies are mainly due to the fact that intramolecular hydrogen bonding energies are overestimated in the lower level calculations. The higher level QM calculations of the potential energy surfaces of d-aldo- and d-ketohexoses now are more comparable to the MM3 results. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1593–1603, 1999

Importance of selecting proper basis set in quantum mechanical studies of potential energy surfaces of carbohydrates

Ab initio calculations give, with an accuracy depending on the sophistication of the method, a bond length as an equilibrium value, re. The experimental bond lengths are always vibrationally averaged and may be expressed in different ways (rg, rz, ra, etc.). Since high-quality ab initio calculations now are capable of giving bond lengths that are approximately of experimental accuracy, it is important to be able to interconvert these values. We find that the bond lengths optimized at the TZ2P+f CCSD level may be considered as the converged re values and that the MM3 and MM4 force fields successfully convert rg to re values. We also evaluated the performance of quantum mechanics at the 6-31G* MP2 and the 6-31G* B3LYP levels and found that the bond lengths (re) at the 6-31G* B3LYP level are better than these at the 6-31G* MP2 level for molecules with only first-row atoms. However, the bond lengths for the bonds involving second-row atoms are too long at the 6-31G* B3LYP level, and for these, the 6-31G* MP2 ...

Systematic Comparison of Experimental, Quantum Mechanical, and Molecular Mechanical Bond Lengths for Organic Molecules

Extensive calculations on hydrogen bonded systems were carried out using the improved MM3 directional hydrogen bond potential. The resulting total function was reoptimized. Comparisons of the hydrogen bonding potential function from ab initio calculations (MP2/6-31G**); the original MM3(89); and the reoptimized MM3 force field MM3(96), for a variety of C, N, O, and Cl systems including the formamide dimer and formamide–water complex, are described herein. Hydrogen bonding is shown to be a far more complicated and ubiquitous phenomenon than is generally recognized. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1001–1016, 1998

Jenn-Huei Lii

Papers

An improved force field (MM4) for saturated hydrocarbons

The MM3 force field for amides, polypeptides and proteins

Importance of selecting proper basis set in quantum mechanical studies of potential energy surfaces of carbohydrates

Systematic Comparison of Experimental, Quantum Mechanical, and Molecular Mechanical Bond Lengths for Organic Molecules

Directional hydrogen bonding in the mm3 force field : ii