MMFF VI. MMFF94s option for energy minimization studies

This article describes the derivation of MMFF94s, which is the “s” (static) variant of MMFF94. MMFF94s incorporates altered out of plane bending parameters that yield planar (or nearly planar) energy‐minimized geometries at unstrained delocalized trigonal nitrogen centers. Some experimental and most theoretical structures show appreciable puckering at nitrogen in isolated structures. However, condensed‐phase effects or even strong intermolecular hydrogen bonding tends to reduce, but need not eliminate, such puckering; in contravention to results reported on the lower level Hartree–Fock surface, we show in the correlated LMP2/6‐31G** calculations for the Watson–Crick guanine–cytosine base pair that one of the hydrogen‐bonded NH2 groups remains appreciably puckered. The resultant MMFF94s geometries emulate the “time‐averaged” structures typically observed in crystallographic and most other experimental structure determinations. MMFF94s also employs modified torsion parameters for interactions that involve such centers, but is identical to MMFF94 for other systems. Isolated instances are found in which MMFF94s fails to locate a (probably shallow) local minimum found on the MMFF94 and reference MP2/6‐31G* surfaces. In general, however, MMFF94s describes conformation energies for delocalized trigonal nitrogen systems nearly as well as MMFF94 does. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 720–729, 1999

[1]  F. Winkler,et al.  Amide group deformation in medium‐ring lactams , 1975 .

[2]  A. Nonat,et al.  Determination of the equilibrium molecular structure of inverting molecules by microwave spectroscopy: Application to aniline , 1986 .

[3]  Martin Karplus,et al.  Ab initio studies of hydrogen bonding of N-methylacetamide: structure, cooperativity, and internal rotational barriers , 1992 .

[4]  Peter A. Kollman,et al.  Theoretical Investigation of the Hydrogen Bond Strengths in Guanine-Cytosine and Adenine-Thymine Base Pairs , 1994 .

[5]  J. Šponer,et al.  Nonplanar geometries of DNA bases. Ab initio second-order Moeller-Plesset study , 1994 .

[6]  Alexander D. MacKerell,et al.  An all-atom empirical energy function for the simulation of nucleic acids , 1995 .

[7]  Richard A. Friesner,et al.  Pseudospectral localized Mo/ller–Plesset methods: Theory and calculation of conformational energies , 1995 .

[8]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[9]  J. Šponer,et al.  Structures and Energies of Hydrogen-Bonded DNA Base Pairs. A Nonempirical Study with Inclusion of Electron Correlation , 1996 .

[10]  Thomas A. Halgren,et al.  Merck molecular force field. IV. conformational energies and geometries for MMFF94 , 1996 .

[11]  Thomas A. Halgren,et al.  Merck molecular force field. III. Molecular geometries and vibrational frequencies for MMFF94 , 1996, J. Comput. Chem..

[12]  T. Halgren.  Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular interactions , 1996 .

[13]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[14]  Richard A. Friesner,et al.  Accurate ab Initio Quantum Chemical Determination of the Relative Energetics of Peptide Conformations and Assessment of Empirical Force Fields , 1997 .

[15]  Jerzy Leszczynski,et al.  THE POTENTIAL ENERGY SURFACE OF GUANINE IS NOT FLAT : AN AB INITIO STUDY WITH LARGE BASIS SETS AND HIGHER ORDER ELECTRON CORRELATION CONTRIBUTIONS , 1998 .

[16]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[17]  T. Halgren.  MMFF VII. Characterization of MMFF94, MMFF94s, and other widely available force fields for conformational energies and for intermolecular‐interaction energies and geometries , 1999, Journal of computational chemistry.