Large-scale compensation of errors in pairwise-additive empirical force fields: comparison of AMBER intermolecular terms with rigorous DFT-SAPT calculations.

The intermolecular interaction energy components for several molecular complexes were calculated using force fields available in the AMBER suite of programs and compared with Density Functional Theory-Symmetry Adapted Perturbation Theory (DFT-SAPT) values. The extent to which such comparison is meaningful is discussed. The comparability is shown to depend strongly on the intermolecular distance, which means that comparisons made at one distance only are of limited value. At large distances the coulombic and van der Waals 1/r(6) empirical terms correspond fairly well with the DFT-SAPT electrostatics and dispersion terms, respectively. At the onset of electronic overlap the empirical values deviate from the reference values considerably. However, the errors in the force fields tend to cancel out in a systematic manner at equilibrium distances. Thus, the overall performance of the force fields displays errors an order of magnitude smaller than those of the individual interaction energy components. The repulsive 1/r(12) component of the van der Waals expression seems to be responsible for a significant part of the deviation of the force field results from the reference values. We suggest that further improvement of the force fields for intermolecular interactions would require replacement of the nonphysical 1/r(12) term by an exponential function. Dispersion anisotropy and its effects are discussed. Our analysis is intended to show that although comparing the empirical and non-empirical interaction energy components is in general problematic, it might bring insights useful for the construction of new force fields. Our results are relevant to often performed force-field-based interaction energy decompositions.

[1]  S. M. Cybulski,et al.  Ab initio calculations of dispersion coefficients for nucleic acid base pairs. , 2006, The Journal of chemical physics.

[2]  T. Halgren,et al.  Polarizable force fields. , 2001, Current opinion in structural biology.

[3]  D. Turner,et al.  Effects of Restrained Sampling Space and Nonplanar Amino Groups on Free-Energy Predictions for RNA with Imino and Sheared Tandem GA Base Pairs Flanked by GC, CG, iGiC or iCiG Base Pairs , 2009, Journal of chemical theory and computation.

[4]  Jirí Cerný,et al.  Density functional theory augmented with an empirical dispersion term. Interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations , 2007, J. Comput. Chem..

[5]  K. Berka,et al.  On the reliability of the AMBER force field and its empirical dispersion contribution for the description of noncovalent complexes. , 2010, Chemphyschem : a European journal of chemical physics and physical chemistry.

[6]  J. Koča,et al.  Docking-based development of purine-like inhibitors of cyclin-dependent kinase-2. , 2000, Journal of medicinal chemistry.

[7]  P. Kollman,et al.  Advancing beyond the atom‐centered model in additive and nonadditive molecular mechanics , 1997 .

[8]  Krzysztof Szalewicz,et al.  Dispersion energy from density-functional theory description of monomers. , 2003, Physical review letters.

[9]  Pavel Hobza,et al.  Highly accurate CCSD(T) and DFT-SAPT stabilization energies of H-bonded and stacked structures of the uracil dimer. , 2008, Chemphyschem : a European journal of chemical physics and physical chemistry.

[10]  R. Podeszwa,et al.  Three-body symmetry-adapted perturbation theory based on Kohn-Sham description of the monomers. , 2007, The Journal of chemical physics.

[11]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[12]  Pedro E. M. Lopes,et al.  Molecular modeling and dynamics studies with explicit inclusion of electronic polarizability: theory and applications , 2009, Theoretical chemistry accounts.

[13]  J. Šponer,et al.  Balance of Attraction and Repulsion in Nucleic-Acid Base Stacking: CCSD(T)/Complete-Basis-Set-Limit Calculations on Uracil Dimer and a Comparison with the Force-Field Description. , 2009, Journal of chemical theory and computation.

[14]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[15]  I. Adamovic,et al.  Dynamic polarizability, dispersion coefficient C6 and dispersion energy in the effective fragment potential method , 2005 .

[16]  J. Koča,et al.  Different Mechanisms of CDK5 and CDK2 Activation as Revealed by CDK5/p25 and CDK2/Cyclin A Dynamics* , 2006, Journal of Biological Chemistry.

[17]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[18]  Jirí Cerný,et al.  Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. , 2006, Physical chemistry chemical physics : PCCP.

[19]  Harry A. Stern,et al.  Fluctuating Charge, Polarizable Dipole, and Combined Models: Parameterization from ab Initio Quantum Chemistry , 1999 .

[20]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[21]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[22]  Georg Jansen,et al.  First-order intermolecular interaction energies from Kohn–Sham orbitals , 2002 .

[23]  Peter Politzer,et al.  σ-hole bonding between like atoms; a fallacy of atomic charges , 2008, Journal of molecular modeling.

[24]  M. Schütz,et al.  Density-functional theory-symmetry-adapted intermolecular perturbation theory with density fitting: a new efficient method to study intermolecular interaction energies. , 2005, The Journal of chemical physics.

[25]  Piotr Cieplak,et al.  Polarization effects in molecular mechanical force fields , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[26]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[27]  N. Špačková,et al.  Theoretical Study of Binding of Hydrated Zn(II) and Mg(II) Cations to 5‘-Guanosine Monophosphate. Toward Polarizable Molecular Mechanics for DNA and RNA , 2003 .

[28]  William L. Jorgensen,et al.  PERFORMANCE OF THE AMBER94, MMFF94, AND OPLS-AA FORCE FIELDS FOR MODELING ORGANIC LIQUIDS , 1996 .

[29]  R. Podeszwa,et al.  Density-Fitting Method in Symmetry-Adapted Perturbation Theory Based on Kohn-Sham Description of Monomers , 2006, 2006 HPCMP Users Group Conference (HPCMP-UGC'06).

[30]  Eric Westhof,et al.  Halogen bonds in biological molecules. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

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

[32]  Robert Moszynski,et al.  Perturbation Theory Approach to Intermolecular Potential Energy Surfaces of van der Waals Complexes , 1994 .

[33]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[34]  Stefan Grimme,et al.  Accurate description of van der Waals complexes by density functional theory including empirical corrections , 2004, J. Comput. Chem..

[35]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[36]  Krzysztof Szalewicz,et al.  Intermolecular potentials based on symmetry-adapted perturbation theory with dispersion energies from time-dependent density-functional calculations. , 2005, The Journal of chemical physics.

[37]  Georg Jansen,et al.  Intermolecular induction and exchange-induction energies from coupled-perturbed Kohn–Sham density functional theory , 2002 .

[38]  Krzysztof Szalewicz,et al.  Intermolecular forces from asymptotically corrected density functional description of monomers , 2002 .

[39]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

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

[41]  Dynamics and Binding Modes of Free cdk2 and its Two Complexes with Inhibitors Studied by Computer Simulations , 2002, Journal of biomolecular structure & dynamics.

[42]  Rafał R Toczyłowski,et al.  An analysis of the electrostatic interaction between nucleic acid bases. , 2005, The Journal of chemical physics.

[43]  Georg Jansen,et al.  Intermolecular dispersion energies from time-dependent density functional theory , 2003 .

[44]  C. Boned,et al.  Molecular dynamics study of the repulsive form influence of the interaction potential on structural, thermodynamic, interfacial, and transport properties. , 2008, The Journal of chemical physics.

[45]  K. Tang,et al.  Erratum: A simple theoretical model for the van der Waals potential at intermediate distances. I. Spherically symmetric potentials , 1977 .

[46]  Pavel Hobza,et al.  Performance of empirical potentials (AMBER, CFF95, CVFF, CHARMM, OPLS, POLTEV), semiempirical quantum chemical methods (AM1, MNDO/M, PM3), and ab initio Hartree–Fock method for interaction of DNA bases: Comparison with nonempirical beyond Hartree–Fock results , 1997 .

[47]  Hirotaka Ode,et al.  Force field parameters for rotation around χ torsion axis in nucleic acids , 2008, J. Comput. Chem..

[48]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[49]  Georg Jansen,et al.  Interaction energy contributions of H-bonded and stacked structures of the AT and GC DNA base pairs from the combined density functional theory and intermolecular perturbation theory approach. , 2006, Journal of the American Chemical Society.

[50]  J. R. Carl,et al.  Atom dipole interaction model for molecular polarizability. Application to polyatomic molecules and determination of atom polarizabilities , 1972 .

[51]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[52]  Steve Scheiner,et al.  Comparison of Morokuma and perturbation theory approaches to decomposition of interaction energy. (NH4)+…NH3 , 1990 .

[53]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..