Bringing the MMFF force field to the RDKit: implementation and validation

A general purpose force field such as MMFF94/MMFF94s, which can properly deal with a wide range of diverse structures, is very valuable in the context of a cheminformatics toolkit. Herein we present an open-source implementation of this force field within the RDKit. The new MMFF functionality can be accessed through a C++/C#/Python/Java application programming interface (API) developed along the lines of the one already available for UFF in the RDKit. Our implementation was fully validated against the official validation suite provided by the MMFF authors. All energies and gradients were correctly computed; moreover, atom type and force constants were correctly assigned for 3D molecules built from SMILES strings. To provide full flexibility, the available API provides direct access to include/exclude individual terms from the MMFF energy expression and to carry out constrained geometry optimizations. The availability of a MMFF-capable molecular mechanics engine coupled with the rest of the RDKit functionality and covered by the BSD license is appealing to researchers operating in both academia and industry.

[1]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[2]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

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

[4]  Klaus Gundertofte,et al.  A comparison of conformational energies calculated by several molecular mechanics methods , 1996, J. Comput. Chem..

[5]  Thomas A. Halgren,et al.  Merck molecular force field. V. Extension of MMFF94 using experimental data, additional computational data, and empirical rules , 1996, J. Comput. Chem..

[6]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[7]  Thomas A. Halgren,et al.  MMFF VII. Characterization of MMFF94, MMFF94s, and other widely available force fields for conformational energies and for intermolecular‐interaction energies and geometries , 1999, J. Comput. Chem..

[8]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

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

[10]  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.

[11]  Thomas A. Halgren,et al.  Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular. interactions , 1996, J. Comput. Chem..

[12]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[13]  Thomas A. Halgren.  MMFF VI. MMFF94s option for energy minimization studies , 1999, J. Comput. Chem..

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

[15]  T. Halgren.  MMFF VI. MMFF94s option for energy minimization studies , 1999, J. Comput. Chem..

[16]  Chris Morley,et al.  Open Babel: An open chemical toolbox , 2011, J. Cheminformatics.