1
|
D’Amore L, Hahn DF, Dotson DL, Horton JT, Anwar J, Craig I, Fox T, Gobbi A, Lakkaraju SK, Lucas X, Meier K, Mobley DL, Narayanan A, Schindler CE, Swope WC, in ’t Veld PJ, Wagner J, Xue B, Tresadern G. Collaborative Assessment of Molecular Geometries and Energies from the Open Force Field. J Chem Inf Model 2022; 62:6094-6104. [PMID: 36433835 PMCID: PMC9873353 DOI: 10.1021/acs.jcim.2c01185] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Force fields form the basis for classical molecular simulations, and their accuracy is crucial for the quality of, for instance, protein-ligand binding simulations in drug discovery. The huge diversity of small-molecule chemistry makes it a challenge to build and parameterize a suitable force field. The Open Force Field Initiative is a combined industry and academic consortium developing a state-of-the-art small-molecule force field. In this report, industry members of the consortium worked together to objectively evaluate the performance of the force fields (referred to here as OpenFF) produced by the initiative on a combined public and proprietary dataset of 19,653 relevant molecules selected from their internal research and compound collections. This evaluation was important because it was completely blind; at most partners, none of the molecules or data were used in force field development or testing prior to this work. We compare the Open Force Field "Sage" version 2.0.0 and "Parsley" version 1.3.0 with GAFF-2.11-AM1BCC, OPLS4, and SMIRNOFF99Frosst. We analyzed force-field-optimized geometries and conformer energies compared to reference quantum mechanical data. We show that OPLS4 performs best, and the latest Open Force Field release shows a clear improvement compared to its predecessors. The performance of established force fields such as GAFF-2.11 was generally worse. While OpenFF researchers were involved in building the benchmarking infrastructure used in this work, benchmarking was done entirely in-house within industrial organizations and the resulting assessment is reported here. This work assesses the force field performance using separate benchmarking steps, external datasets, and involving external research groups. This effort may also be unique in terms of the number of different industrial partners involved, with 10 different companies participating in the benchmark efforts.
Collapse
Affiliation(s)
- Lorenzo D’Amore
- Computational Chemistry, Janssen R&D, C/ Jarama 75A, 45007 Toledo, Spain
| | - David F. Hahn
- Computational Chemistry, Janssen R&D, Turnhoutseweg 30, Beerse B-2340, Belgium
| | - David L. Dotson
- The Open Force Field Initiative, Open Molecular Software Foundation, Davis, California 95616, USA
| | - Joshua T. Horton
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | - Jamshed Anwar
- Department of Chemistry, Lancaster University, Lancaster LA1 4YW, UK
| | - Ian Craig
- Molecular Modeling & Drug Discovery, BASF SE, 67056 Ludwigshafen, Germany
| | - Thomas Fox
- Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co KG, 88397 Biberach/Riss, Germany
| | - Alberto Gobbi
- Genentech, Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | | | - Xavier Lucas
- Roche Pharma Research and Early Development, Therapeutic Modalities, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070 Basel, Switzerland
| | - Katharina Meier
- Computational Life Science Technology Functions, Crop Science, R&D, Bayer AG, 40789 Monheim, Germany
| | - David L. Mobley
- Departments of Pharmaceutical Sciences and Chemistry, University of California 92617, Irvine, USA
| | - Arjun Narayanan
- Data and Computational Sciences, Vertex Pharmaceuticals, 50 Northern Ave, Boston, MA 02210, USA
| | | | - William C. Swope
- Genentech, Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | | | - Jeffrey Wagner
- The Open Force Field Initiative, Open Molecular Software Foundation, Davis, California, 95616, USA,Chemistry Department, The University of California at Irvine, Irvine, California, 92617, USA
| | - Bai Xue
- XtalPi Inc. Floor 3, International Biomedical Innovation Park II, No. 2 Hongliu Road, Fubao Community, Fubao Street, Futian District, Shenzhen, Guangdong, 518040 China
| | - Gary Tresadern
- Computational Chemistry, Janssen R&D, Turnhoutseweg 30, Beerse B-2340, Belgium
| |
Collapse
|
2
|
Lim VT, Hahn DF, Tresadern G, Bayly CI, Mobley DL. Benchmark assessment of molecular geometries and energies from small molecule force fields. F1000Res 2020; 9:Chem Inf Sci-1390. [PMID: 33604023 PMCID: PMC7863993 DOI: 10.12688/f1000research.27141.1] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 11/18/2020] [Indexed: 12/22/2022] Open
Abstract
Background: Force fields are used in a wide variety of contexts for classical molecular simulation, including studies on protein-ligand binding, membrane permeation, and thermophysical property prediction. The quality of these studies relies on the quality of the force fields used to represent the systems. Methods: Focusing on small molecules of fewer than 50 heavy atoms, our aim in this work is to compare nine force fields: GAFF, GAFF2, MMFF94, MMFF94S, OPLS3e, SMIRNOFF99Frosst, and the Open Force Field Parsley, versions 1.0, 1.1, and 1.2. On a dataset comprising 22,675 molecular structures of 3,271 molecules, we analyzed force field-optimized geometries and conformer energies compared to reference quantum mechanical (QM) data. Results: We show that while OPLS3e performs best, the latest Open Force Field Parsley release is approaching a comparable level of accuracy in reproducing QM geometries and energetics for this set of molecules. Meanwhile, the performance of established force fields such as MMFF94S and GAFF2 is generally somewhat worse. We also find that the series of recent Open Force Field versions provide significant increases in accuracy. Conclusions: This study provides an extensive test of the performance of different molecular mechanics force fields on a diverse molecule set, and highlights two (OPLS3e and OpenFF 1.2) that perform better than the others tested on the present comparison. Our molecule set and results are available for other researchers to use in testing.
Collapse
Affiliation(s)
- Victoria T. Lim
- Department of Chemistry, University of California, Irvine, CA, 92697, USA
| | - David F. Hahn
- Computational Chemistry, Janssen Research & Development, Beerse, B-2340, Belgium
| | - Gary Tresadern
- Computational Chemistry, Janssen Research & Development, Beerse, B-2340, Belgium
| | | | - David L. Mobley
- Department of Chemistry, University of California, Irvine, CA, 92697, USA
- Department of Pharmaceutical Sciences, University of California, Irvine, CA, 92697, USA
| |
Collapse
|
3
|
Valadez Huerta G, Raabe G. Genetic Parameterization of Interfacial Force Fields Based on Classical Bulk Force Fields and Ab Initio Data: Application to the Methanol-ZnO Interfaces. J Chem Inf Model 2020; 60:6033-6043. [DOI: 10.1021/acs.jcim.0c01093] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Gerardo Valadez Huerta
- Institut für Thermodynamik, Technische Universität Braunschweig, Hans-Sommer-Straße 5, D-38106 Braunschweig, Germany
| | - Gabriele Raabe
- Institut für Thermodynamik, Technische Universität Braunschweig, Hans-Sommer-Straße 5, D-38106 Braunschweig, Germany
| |
Collapse
|
4
|
Kulkarni A, Fingerhut R, Kohns M, Hasse H, Vrabec J. Correction to “Molecular Modeling and Simulation of Vapor–Liquid Equilibria of Ethylene Oxide, Ethylene Glycol, and Water as Well as their Binary Mixtures”. Ind Eng Chem Res 2020. [DOI: 10.1021/acs.iecr.0c04937] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Aditya Kulkarni
- Laboratory of Engineering Thermodynamics (LTD), Technische Universität Kaiserslautern (TUK), 67663 Kaiserslautern, Germany
| | - Robin Fingerhut
- Thermodynamics and Process Engineering, Technical University of Berlin, 10587 Berlin, Germany
| | - Maximilian Kohns
- Laboratory of Engineering Thermodynamics (LTD), Technische Universität Kaiserslautern (TUK), 67663 Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), Technische Universität Kaiserslautern (TUK), 67663 Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technical University of Berlin, 10587 Berlin, Germany
| |
Collapse
|
5
|
Tsimpanogiannis IN, Jamali SH, Economou IG, Vlugt TJH, Moultos OA. On the validity of the Stokes–Einstein relation for various water force fields. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1702729] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Chemical Process & Energy Resources Institute (CPERI), Centre for Research & Technology Hellas (CERTH) Thermi-Thessaloniki, Greece
| | - Seyed Hossein Jamali
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | | | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
6
|
Kaliannan NK, Henao Aristizabal A, Wiebeler H, Zysk F, Ohto T, Nagata Y, Kühne TD. Impact of intermolecular vibrational coupling effects on the sum-frequency generation spectra of the water/air interface. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1620358] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Naveen Kumar Kaliannan
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Andres Henao Aristizabal
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Hendrik Wiebeler
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Frederik Zysk
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
| | - Tatsuhiko Ohto
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, Japan
| | - Yuki Nagata
- Max-Planck Institute for Polymer Research, Mainz, Germany
| | - Thomas D. Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Chair of Theoretical Chemistry, University of Paderborn, Paderborn, Germany
- Paderborn Center for Parallel Computing and Institute for Lightweight Design, University of Paderborn, Paderborn, Germany
| |
Collapse
|
7
|
Zanette C, Bannan CC, Bayly CI, Fass J, Gilson MK, Shirts MR, Chodera JD, Mobley DL. Toward Learned Chemical Perception of Force Field Typing Rules. J Chem Theory Comput 2019; 15:402-423. [PMID: 30512951 PMCID: PMC6467725 DOI: 10.1021/acs.jctc.8b00821] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular mechanics force fields define how the energy and forces in a molecular system are computed from its atomic positions, thus enabling the study of such systems through computational methods like molecular dynamics and Monte Carlo simulations. Despite progress toward automated force field parametrization, considerable human expertise is required to develop or extend force fields. In particular, human input has long been required to define atom types, which encode chemically unique environments that determine which parameters will be assigned. However, relying on humans to establish atom types is suboptimal. Human-created atom types are often developed without statistical justification, leading to over- or under-fitting of data. Human-created types are also difficult to extend in a systematic and consistent manner when new chemistries must be modeled or new data becomes available. Finally, human effort is not scalable when force fields must be generated for new (bio)polymers, compound classes, or materials. To remedy these deficiencies, our long-term goal is to replace human specification of atom types with an automated approach, based on rigorous statistics and driven by experimental and/or quantum chemical reference data. In this work, we describe novel methods that automate the discovery of appropriate chemical perception: SMARTY allows for the creation of atom types, while SMIRKY goes further by automating the creation of fragment (nonbonded, bonds, angles, and torsions) types. These approaches enable the creation of move sets in atom or fragment type space, which are used within a Monte Carlo optimization approach. We demonstrate the power of these new methods by automating the rediscovery of human defined atom types (SMARTY) or fragment types (SMIRKY) in existing small molecule force fields. We assess these approaches using several molecular data sets, including one which covers a diverse subset of the DrugBank database.
Collapse
Affiliation(s)
- Camila Zanette
- Department of Pharmaceutical Sciences, University of California, Irvine
| | | | | | - Josh Fass
- Tri-Institutional Training Program in Computational Biology and Medicine, New York, NY 10065
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - Michael K. Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego
| | - Michael R. Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, CO 80309
| | - John D. Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - David L. Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine
- Department of Chemistry, University of California, Irvine
| |
Collapse
|
8
|
Tsimpanogiannis IN, Moultos OA, Franco LFM, Spera MBDM, Erdős M, Economou IG. Self-diffusion coefficient of bulk and confined water: a critical review of classical molecular simulation studies. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1511903] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Ioannis N. Tsimpanogiannis
- Environmental Research Laboratory, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Luís F. M. Franco
- School of Chemical Engineering, University of Campinas, Campinas, Brazil
| | | | - Máté Erdős
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Ioannis G. Economou
- Institute of Nanoscience and Nanotechnology, National Center for Scientific Research “Demokritos”, Aghia Paraskevi Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar, Doha, Qatar
| |
Collapse
|
9
|
Ojha D, Henao A, Kühne TD. Nuclear quantum effects on the vibrational dynamics of liquid water. J Chem Phys 2018; 148:102328. [PMID: 29544291 DOI: 10.1063/1.5005500] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Based on quantum-mechanical path-integral molecular dynamics simulations, the impact of nuclear quantum effects on the vibrational and hydrogen bond dynamics in liquid water is investigated. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time-series analysis, while the time scales of the vibrational spectral diffusion are determined from frequency-time correlation functions, joint probability distributions, and the slope of three-pulse photon echo. We find that the inclusion of nuclear quantum effects leads not only to a redshift of the vibrational frequency distribution by around 130 cm-1 but also to an acceleration of the vibrational dynamics by as much as 30%. In addition, quantum fluctuations also entail a significantly faster decay of correlation in the initial diffusive regime, which is in agreement with recent vibrational echo experiments.
Collapse
Affiliation(s)
- Deepak Ojha
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Andrés Henao
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Thomas D Kühne
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| |
Collapse
|
10
|
Mausbach P, Köster A, Rutkai G, Thol M, Vrabec J. Comparative study of the Grüneisen parameter for 28 pure fluids. J Chem Phys 2016; 144:244505. [DOI: 10.1063/1.4954282] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|