1
|
Oliveira MP, Hünenberger PH. Force fields optimized against experimental data for large compound families using CombiFF: Validation considering non-target properties and polyfunctional compounds. J Mol Graph Model 2023; 118:108312. [PMID: 36252318 DOI: 10.1016/j.jmgm.2022.108312] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 08/09/2022] [Accepted: 08/16/2022] [Indexed: 10/14/2022]
Abstract
The CombiFF scheme is a workflow for the automated calibration of force-field parameters against condensed-phase experimental data considering simultaneously entire classes of organic molecules. The main steps of this scheme are: (i) selection of a molecule family; (ii) enumeration of all isomers; (iii) query for experimental data; (iv) automatic construction of the molecular topologies; (v) iterative refinement of the force-field parameters considering the entire family. In two recent articles, CombiFF was applied to the design of GROMOS-compatible united-atom force fields for the saturated acyclic haloalkanes and for saturated acyclic compounds involving eight common chemical functional groups of oxygen and nitrogen. This calibration and the subsequent initial validation involved two limitations: (i) the experimental data considered was restricted to values for the pure-liquid density ρliq and the vaporization enthalpy ΔHvap of the compounds; (ii) beyond monofunctional compounds, the training set only involved homo-polyhaloalkanes (possibly mixing halogen types) in the first study, and homo-polyfunctional compounds of the considered oxygen or nitrogen functional groups (no mixing of different group types) in the second one. The goal of this article is to further test the accuracy of CombiFF-generated force fields by extending the validation to: (i) nine additional properties that were not used as optimization targets (pure-liquid thermodynamic, dielectric and transport properties, as well as solvation properties); (ii) hetero-polyfunctional molecules that were not included in the calibration and initial validation sets. The results for the nine additional properties show good agreement with experiment, except for the shear viscosity and the dielectric permittivity. There, larger discrepancies are observed, likely due to the united-atom representation adopted for the aliphatic groups and to the implicit treatment of electronic polarization effects. The results for the hetero-polyfunctional molecules also show reasonable agreement with experiment in terms of the monitored properties.
Collapse
Affiliation(s)
- Marina P Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093, Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093, Zürich, Switzerland.
| |
Collapse
|
2
|
Kashefolgheta S, Wang S, Acree WE, Hünenberger PH. Evaluation of nine condensed-phase force fields of the GROMOS, CHARMM, OPLS, AMBER, and OpenFF families against experimental cross-solvation free energies. Phys Chem Chem Phys 2021; 23:13055-13074. [PMID: 34105547 PMCID: PMC8207520 DOI: 10.1039/d1cp00215e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 04/28/2021] [Indexed: 12/02/2022]
Abstract
Experimental solvation free energies are nowadays commonly included as target properties in the validation of condensed-phase force fields, sometimes even in their calibration. In a previous article [Kashefolgheta et al., J. Chem. Theory. Comput., 2020, 16, 7556-7580], we showed how the involved comparison between experimental and simulation results could be made more systematic by considering a full matrix of cross-solvation free energies . For a set of N molecules that are all in the liquid state under ambient conditions, such a matrix encompasses N×N entries for considering each of the N molecules either as solute (A) or as solvent (B). In the quoted study, a cross-solvation matrix of 25 × 25 experimental value was introduced, considering 25 small molecules representative for alkanes, chloroalkanes, ethers, ketones, esters, alcohols, amines, and amides. This experimental data was used to compare the relative accuracies of four popular condensed-phase force fields, namely GROMOS-2016H66, OPLS-AA, AMBER-GAFF, and CHARMM-CGenFF. In the present work, the comparison is extended to five additional force fields, namely GROMOS-54A7, GROMOS-ATB, OPLS-LBCC, AMBER-GAFF2, and OpenFF. Considering these nine force fields, the correlation coefficients between experimental values and simulation results range from 0.76 to 0.88, the root-mean-square errors (RMSEs) from 2.9 to 4.8 kJ mol-1, and average errors (AVEEs) from -1.5 to +1.0 kJ mol-1. In terms of RMSEs, GROMOS-2016H66 and OPLS-AA present the best accuracy (2.9 kJ mol-1), followed by OPLS-LBCC, AMBER-GAFF2, AMBER-GAFF, and OpenFF (3.3 to 3.6 kJ mol-1), and then by GROMOS-54A7, CHARM-CGenFF, and GROMOS-ATB (4.0 to 4.8 kJ mol-1). These differences are statistically significant but not very pronounced, and are distributed rather heterogeneously over the set of compounds within the different force fields.
Collapse
Affiliation(s)
- Sadra Kashefolgheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| | - Shuzhe Wang
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| | - William E. Acree
- Department of Chemistry, University of North Texas1155 Union Circle Drive #305070DentonTexas 76203USA
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| |
Collapse
|
3
|
Kashefolgheta S, Oliveira MP, Rieder SR, Horta BAC, Acree WE, Hünenberger PH. Evaluating Classical Force Fields against Experimental Cross-Solvation Free Energies. J Chem Theory Comput 2020; 16:7556-7580. [DOI: 10.1021/acs.jctc.0c00688] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Sadra Kashefolgheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - William E. Acree
- Department of Chemistry, University of North Texas, 1155 Union Circle Drive #305070, Denton, Texas 76203, United States
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
4
|
König G, Riniker S. On the faithfulness of molecular mechanics representations of proteins towards quantum-mechanical energy surfaces. Interface Focus 2020; 10:20190121. [PMID: 33184586 DOI: 10.1098/rsfs.2019.0121] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/06/2020] [Indexed: 12/11/2022] Open
Abstract
Force fields based on molecular mechanics (MM) are the main computational tool to study the relationship between protein structure and function at the molecular level. To validate the quality of such force fields, high-level quantum-mechanical (QM) data are employed to test their capability to reproduce the features of all major conformational substates of a series of blocked amino acids. The phase-space overlap between MM and QM is quantified in terms of the average structural reorganization energies over all energy minima. Here, the structural reorganization energy is the MM potential-energy difference between the structure of the respective QM energy minimum and the structure of the closest MM energy minimum. Thus, it serves as a measure for the relative probability of visiting the QM minimum during an MM simulation. We evaluate variants of the AMBER, CHARMM, GROMOS and OPLS biomolecular force fields. In addition, the two blocked amino acids alanine and serine are used to demonstrate the dependence of the measured agreement on the QM method, the phase, and the conformational preferences. Blocked serine serves as an example to discuss possible improvements of the force fields, such as including polarization with Drude particles, or using tailored force fields. The results show that none of the evaluated force fields satisfactorily reproduces all energy minima. By decomposing the average structural reorganization energies in terms of individual energy terms, we can further assess the individual weaknesses of the parametrization strategies of each force field. The dominant problem for most force fields appears to be the van der Waals parameters, followed to a lesser degree by dihedral and bonded terms. Our results show that performing a simple QM energy optimization from an MM-optimized structure can be a first test of the validity of a force field for a particular target molecule.
Collapse
Affiliation(s)
- Gerhard König
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.,Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
5
|
Ramos MC, Horta VAC, Horta BAC. Molecular Dynamics Simulations of PAMAM and PPI Dendrimers Using the GROMOS-Compatible 2016H66 Force Field. J Chem Inf Model 2019; 59:1444-1457. [DOI: 10.1021/acs.jcim.8b00911] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Mayk C. Ramos
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| | - Vitor A. C. Horta
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Juiz de Fora, 36036-900, Brazil
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| |
Collapse
|
6
|
Biomolecular force fields: where have we been, where are we now, where do we need to go and how do we get there? J Comput Aided Mol Des 2018; 33:133-203. [DOI: 10.1007/s10822-018-0111-4] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2017] [Accepted: 03/09/2018] [Indexed: 10/27/2022]
|
7
|
Horta BAC, Merz PT, Fuchs PFJ, Dolenc J, Riniker S, Hünenberger PH. A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set. J Chem Theory Comput 2016; 12:3825-50. [DOI: 10.1021/acs.jctc.6b00187] [Citation(s) in RCA: 88] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bruno A. C. Horta
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Pascal T. Merz
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Patrick F. J. Fuchs
- Institut Jacques Monod, UMR 7592 CNRS, Université Paris-Diderot, Sorbonne Paris Cité, F-75205 Paris, France
| | - Jozica Dolenc
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Chemistry,
Biology and Pharmacy Information Center, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | | |
Collapse
|
8
|
Zhang J, Tuguldur B, van der Spoel D. Force Field Benchmark of Organic Liquids. 2. Gibbs Energy of Solvation. J Chem Inf Model 2015; 55:1192-201. [DOI: 10.1021/acs.jcim.5b00106] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jin Zhang
- Department of Chemistry, Zhejiang University, Hangzhou 310027, China
- Uppsala Center for Computational Chemistry, Science for Life Laboratory,
Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Badamkhatan Tuguldur
- Uppsala Center for Computational Chemistry, Science for Life Laboratory,
Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
- Department of Biology,
School of Arts and Sciences, National University of Mongolia, Ulaanbaatar 14200, Mongolia
| | - David van der Spoel
- Uppsala Center for Computational Chemistry, Science for Life Laboratory,
Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
9
|
Huang W, Lin Z, van Gunsteren WF. Use of Enveloping Distribution Sampling to Evaluate Important Characteristics of Biomolecular Force Fields. J Phys Chem B 2014; 118:6424-30. [DOI: 10.1021/jp411005x] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Wei Huang
- Laboratory of Physical Chemistry,
Swiss Federal Institute of Technology, ETH, Vladimir-Prelog-Weg 2/HCI, 8093 Zürich, Switzerland
| | - Zhixiong Lin
- Laboratory of Physical Chemistry,
Swiss Federal Institute of Technology, ETH, Vladimir-Prelog-Weg 2/HCI, 8093 Zürich, Switzerland
| | - Wilfred F. van Gunsteren
- Laboratory of Physical Chemistry,
Swiss Federal Institute of Technology, ETH, Vladimir-Prelog-Weg 2/HCI, 8093 Zürich, Switzerland
| |
Collapse
|
10
|
Lin Z, van Gunsteren WF. Enhanced conformational sampling using enveloping distribution sampling. J Chem Phys 2013; 139:144105. [DOI: 10.1063/1.4824391] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Lin Z, van Gunsteren WF. Refinement of the application of the GROMOS 54A7 force field to β-peptides. J Comput Chem 2013; 34:2796-805. [DOI: 10.1002/jcc.23459] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2013] [Revised: 09/12/2013] [Accepted: 09/13/2013] [Indexed: 11/08/2022]
Affiliation(s)
- Zhixiong Lin
- Laboratory of Physical Chemistry; Swiss Federal Institute of Technology; ETH 8093 Zürich Switzerland
| | - Wilfred F. van Gunsteren
- Laboratory of Physical Chemistry; Swiss Federal Institute of Technology; ETH 8093 Zürich Switzerland
| |
Collapse
|
12
|
Pechlaner M, Sigel RKO, van Gunsteren WF, Dolenc J. Structure and Conformational Dynamics of the Domain 5 RNA Hairpin of a Bacterial Group II Intron Revealed by Solution Nuclear Magnetic Resonance and Molecular Dynamics Simulations. Biochemistry 2013; 52:7099-113. [DOI: 10.1021/bi400784r] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Maria Pechlaner
- Institute
of Inorganic Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Roland K. O. Sigel
- Institute
of Inorganic Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Wilfred F. van Gunsteren
- Laboratory
of Physical Chemistry, Swiss Federal Institute of Technology, CH-8093 Zurich, Switzerland
| | - Jožica Dolenc
- Laboratory
of Physical Chemistry, Swiss Federal Institute of Technology, CH-8093 Zurich, Switzerland
| |
Collapse
|