1
|
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
|
2
|
Sharapa DI, Genaev A, Cavallo L, Minenkov Y. A Robust and Cost‐Efficient Scheme for Accurate Conformational Energies of Organic Molecules. Chemphyschem 2018; 20:92-102. [DOI: 10.1002/cphc.201801063] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2018] [Indexed: 11/09/2022]
Affiliation(s)
- Dmitry I. Sharapa
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology (KIT) Hermann-von-Helmholtz Platz 1 Eggenstein-Leopoldshafen D-76344 Germany
| | - Alexander Genaev
- Vorozhtsov Novosibirsk Institute of Organic Chemistry Academician Lavrent'ev Ave., 9 Novosibirsk 630090 Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC)King Abdullah University of Science and Technology Thuwal- 23955-6900 Saudi Arabia
| | - Yury Minenkov
- Moscow Institute of Physics and Technology Institutskiy Pereulok 9, Dolgoprudny Moscow Region 141700 Russia
| |
Collapse
|
3
|
Oliver A, Canals V, Rosselló JL. A Bayesian Target Predictor Method based on Molecular Pairing Energies estimation. Sci Rep 2017; 7:43738. [PMID: 28263323 PMCID: PMC5338323 DOI: 10.1038/srep43738] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Accepted: 01/30/2017] [Indexed: 11/13/2022] Open
Abstract
Virtual screening (VS) is applied in the early drug discovery phases for the quick inspection of huge molecular databases to identify those compounds that most likely bind to a given drug target. In this context, there is the necessity of the use of compact molecular models for database screening and precise target prediction in reasonable times. In this work we present a new compact energy-based model that is tested for its application to Virtual Screening and target prediction. The model can be used to quickly identify active compounds in huge databases based on the estimation of the molecule’s pairing energies. The greatest molecular polar regions along with its geometrical distribution are considered by using a short set of smart energy vectors. The model is tested using similarity searches within the Directory of Useful Decoys (DUD) database. The results obtained are considerably better than previously published models. As a Target prediction methodology we propose the use of a Bayesian Classifier that uses a combination of different active compounds to build an energy-dependent probability distribution function for each target.
Collapse
Affiliation(s)
- Antoni Oliver
- Physics Department, Universitat de les Illes Balears, Palma de Mallorca, Spain
| | - Vincent Canals
- Physics Department, Universitat de les Illes Balears, Palma de Mallorca, Spain
| | - Josep L Rosselló
- Physics Department, Universitat de les Illes Balears, Palma de Mallorca, Spain
| |
Collapse
|
4
|
Ueno-Noto K, Takano K. Water molecules inside protein structure affect binding of monosaccharides with HIV-1 antibody 2G12. J Comput Chem 2016; 37:2341-8. [DOI: 10.1002/jcc.24447] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2016] [Revised: 06/03/2016] [Accepted: 06/16/2016] [Indexed: 11/07/2022]
Affiliation(s)
- Kaori Ueno-Noto
- Chemistry Section, Center for Natural Sciences, College of Liberal Arts and Sciences, Kitasato University; 1-15-1 Kitasato Minami-Ku, Sagamihara Kanagawa 252-0373 Japan
| | - Keiko Takano
- Department of Chemistry and Biochemistry; Graduate School of Humanities and Sciences, Ochanomizu University; 2-1-1 Otsuka, Bunkyo-Ku Tokyo 112-8610 Japan
| |
Collapse
|
5
|
Ahmed A, Sandler SI. Predictions of the physicochemical properties of amino acid side chain analogs using molecular simulation. Phys Chem Chem Phys 2016; 18:6559-68. [DOI: 10.1039/c5cp05393e] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
A candidate drug compound is released for clinical trails (in vivo activity) only if its physicochemical properties meet desirable bioavailability and partitioning criteria.
Collapse
Affiliation(s)
- Alauddin Ahmed
- Center for Molecular and Engineering Thermodynamics
- Department of Chemical and Biomolecular Engineering
- University of Delaware
- Newark
- USA
| | - Stanley I. Sandler
- Center for Molecular and Engineering Thermodynamics
- Department of Chemical and Biomolecular Engineering
- University of Delaware
- Newark
- USA
| |
Collapse
|
6
|
Caleman C, van Maaren PJ, Hong M, Hub JS, Costa LT, van der Spoel D. Force Field Benchmark of Organic Liquids: Density, Enthalpy of Vaporization, Heat Capacities, Surface Tension, Isothermal Compressibility, Volumetric Expansion Coefficient, and Dielectric Constant. J Chem Theory Comput 2012; 8:61-74. [PMID: 22241968 PMCID: PMC3254193 DOI: 10.1021/ct200731v] [Citation(s) in RCA: 463] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2011] [Indexed: 12/15/2022]
Abstract
The chemical composition of small organic molecules is often very similar to amino acid side chains or the bases in nucleic acids, and hence there is no a priori reason why a molecular mechanics force field could not describe both organic liquids and biomolecules with a single parameter set. Here, we devise a benchmark for force fields in order to test the ability of existing force fields to reproduce some key properties of organic liquids, namely, the density, enthalpy of vaporization, the surface tension, the heat capacity at constant volume and pressure, the isothermal compressibility, the volumetric expansion coefficient, and the static dielectric constant. Well over 1200 experimental measurements were used for comparison to the simulations of 146 organic liquids. Novel polynomial interpolations of the dielectric constant (32 molecules), heat capacity at constant pressure (three molecules), and the isothermal compressibility (53 molecules) as a function of the temperature have been made, based on experimental data, in order to be able to compare simulation results to them. To compute the heat capacities, we applied the two phase thermodynamics method (Lin et al. J. Chem. Phys.2003, 119, 11792), which allows one to compute thermodynamic properties on the basis of the density of states as derived from the velocity autocorrelation function. The method is implemented in a new utility within the GROMACS molecular simulation package, named g_dos, and a detailed exposé of the underlying equations is presented. The purpose of this work is to establish the state of the art of two popular force fields, OPLS/AA (all-atom optimized potential for liquid simulation) and GAFF (generalized Amber force field), to find common bottlenecks, i.e., particularly difficult molecules, and to serve as a reference point for future force field development. To make for a fair playing field, all molecules were evaluated with the same parameter settings, such as thermostats and barostats, treatment of electrostatic interactions, and system size (1000 molecules). The densities and enthalpy of vaporization from an independent data set based on simulations using the CHARMM General Force Field (CGenFF) presented by Vanommeslaeghe et al. (J. Comput. Chem.2010, 31, 671) are included for comparison. We find that, overall, the OPLS/AA force field performs somewhat better than GAFF, but there are significant issues with reproduction of the surface tension and dielectric constants for both force fields.
Collapse
Affiliation(s)
- Carl Caleman
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron Notkestraße 85, DE-22607 Hamburg, Germany
| | - Paul J. van Maaren
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Minyan Hong
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Jochen S. Hub
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Luciano T. Costa
- Departamento de Ciências Exatas, Federal University of Alfenas—MG Rua Gabriel Monteiro da Silva, 700 Alfenas—MG CEP:37130-000, Brazil
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
7
|
Cheong PHY, Legault CY, Um JM, Çelebi-Ölçüm N, Houk KN. Quantum mechanical investigations of organocatalysis: mechanisms, reactivities, and selectivities. Chem Rev 2011; 111:5042-137. [PMID: 21707120 PMCID: PMC3154597 DOI: 10.1021/cr100212h] [Citation(s) in RCA: 429] [Impact Index Per Article: 33.0] [Reference Citation Analysis] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Paul Ha-Yeon Cheong
- Oregon State University, Department of Chemistry, 153 Gilbert Hall, Corvallis, OR 97331-4003 USA
| | - Claude Y. Legault
- University of Sherbrooke, Department of Chemistry, 2500 boul. de l’Université, local D1-3029, Sherbrooke (Québec) J1K 2R1 CANADA
| | - Joann M. Um
- University of California Los Angeles, Department of Chemistry and Biochemistry, 607 Charles E. Young Drive East, Los Angeles, CA 90095-1569 USA
| | - Nihan Çelebi-Ölçüm
- University of California Los Angeles, Department of Chemistry and Biochemistry, 607 Charles E. Young Drive East, Los Angeles, CA 90095-1569 USA
| | - K. N. Houk
- University of California Los Angeles, Department of Chemistry and Biochemistry, 607 Charles E. Young Drive East, Los Angeles, CA 90095-1569 USA
| |
Collapse
|
8
|
Simón L, Goodman JM. How reliable are DFT transition structures? Comparison of GGA, hybrid-meta-GGA and meta-GGA functionals. Org Biomol Chem 2010; 9:689-700. [PMID: 20976314 DOI: 10.1039/c0ob00477d] [Citation(s) in RCA: 184] [Impact Index Per Article: 13.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
There have been many comparisons of computational methods applied to ground states, but studies of organic reactions usually require calculations on transition states, and these provide a different test of the methods. We present calculations of the geometries of nineteen covalent-bond forming transition states using HF and twelve different functionals, including GGA, hybrid-GGA and hybrid meta-GGA approaches. For the calculation of the TS geometries, the results suggest that B3LYP is only slightly less accurate than newer, computationally more expensive methods, and is less sensitive to choice of integration grid. We conclude that the use of B3LYP and related functionals is still appropriate for many studies of organic reaction mechanisms.
Collapse
Affiliation(s)
- Luis Simón
- Facultad de Ciencias Químicas, Universidad de Salamanca, Plaza de los Caídos 1-5, Salamanca, E37004, Spain.
| | | |
Collapse
|
9
|
Malo M, Brive L, Luthman K, Svensson P. Selective pharmacophore models of dopamine D(1) and D(2) full agonists based on extended pharmacophore features. ChemMedChem 2010; 5:232-46. [PMID: 20077461 DOI: 10.1002/cmdc.200900398] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
This study is focused on the identification of structural features that determine the selectivity of dopamine receptor agonists toward D(1) and D(2) receptors. Selective pharmacophore models were developed for both receptors. The models were built by using projected pharmacophoric features that represent the main agonist interaction sites in the receptor (the Ser residues in TM5 and the Asp in TM3), a directional aromatic feature in the ligand, a feature with large positional tolerance representing the positively charged nitrogen in the ligand, and sets of excluded volumes reflecting the shapes of the receptors. The sets of D(1) and D(2) ligands used for modeling were carefully selected from published sources and consist of structurally diverse, conformationally rigid full agonists as active ligands together with structurally related inactives. The robustness of the models in discriminating actives from inactives was tested against four ensembles of conformations generated by using different established methods and different force fields. The reasons for the selectivity can be attributed to both geometrical differences in the arrangement of the features, e.g., different tilt angels of the pi system, as well as shape differences covered by the different sets of excluded volumes. This work provides useful information for the design of new D(1) and D(2) agonists and also for comparative homology modeling of D(1) and D(2) receptors. The approach is general and could therefore be applied to other ligand-protein interactions for which no experimental protein structure is available.
Collapse
Affiliation(s)
- Marcus Malo
- Department of Chemistry, Medicinal Chemistry, University of Gothenburg, 41296 Göteborg, Sweden
| | | | | | | |
Collapse
|
10
|
James WH, Baquero EE, Shubert VA, Choi SH, Gellman SH, Zwier TS. Single-Conformation and Diastereomer Specific Ultraviolet and Infrared Spectroscopy of Model Synthetic Foldamers: α/β-Peptides. J Am Chem Soc 2009; 131:6574-90. [PMID: 19366210 DOI: 10.1021/ja901051v] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- William H. James
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - Esteban E. Baquero
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - V. Alvin Shubert
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - Soo Hyuk Choi
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - Samuel H. Gellman
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| | - Timothy S. Zwier
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907-2084, and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706
| |
Collapse
|