1
|
Fortuna A, Costa PJ. Assessment of Halogen Off-Center Point-Charge Models Using Explicit Solvent Simulations. J Chem Inf Model 2023; 63:7464-7475. [PMID: 38010191 DOI: 10.1021/acs.jcim.3c01561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Compounds containing halogens can form halogen bonds (XBs) with biological targets such as proteins and membranes due to their anisotropic electrostatic potential. To accurately describe this anisotropy, off-center point-charge (EP) models are commonly used in force field methods, allowing the description of XBs at the molecular mechanics and molecular dynamics level. Various EP implementations have been documented in the literature, and despite being efficient in reproducing protein-ligand geometries and sampling of XBs, it is unclear how well these EP models predict experimental properties such as hydration free energies (ΔGhyd), which are often used to validate force field performance. In this work, we report the first assessment of three EP models using alchemical free energy calculations to predict ΔGhyd values. We show that describing the halogen anisotropy using some EP models can lead to a slight improvement in the prediction of the ΔGhyd when compared with the models without EP, especially for the chlorinated compounds; however, this improvement is not related to the establishment of XBs but is most likely due to the improvement of the sampling of hydrogen bonds. We also highlight the importance of the choice of the EP model, especially for the iodinated molecules, since a slight tendency to improve the prediction is observed for compounds with a larger σ-hole but significantly worse results were obtained for compounds that are weaker XB donors.
Collapse
Affiliation(s)
- Andreia Fortuna
- BioISI─Instituto de Biosistemas e Ciências Integrativas, Faculdade de Ciências, Universidade de Lisboa, Lisboa 1749-016, Portugal
- Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, Lisboa 1749-016, Portugal
- Research Institute for Medicines (iMed.ULisboa), Faculty of Pharmacy, University of Lisbon, Av. Professor Gama Pinto, Lisbon 1649-003, Portugal
| | - Paulo J Costa
- BioISI─Instituto de Biosistemas e Ciências Integrativas, Faculdade de Ciências, Universidade de Lisboa, Lisboa 1749-016, Portugal
- Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, Lisboa 1749-016, Portugal
| |
Collapse
|
2
|
Pal P, Chakraborty S, Jana B. Number of Hydrogen Bonds per Unit Solvent Accessible Surface Area: A Descriptor of Functional States of Proteins. J Phys Chem B 2022; 126:10822-10833. [PMID: 36524238 DOI: 10.1021/acs.jpcb.2c05367] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Proteins function close to native and near-native conformations. These states are evolutionarily selected to ensure the effect of mutations is minimized. The structural organization of a protein is hierarchical and modular, which reduces the dimensionality of the configurational space of the native states. Thus, finding appropriate descriptors that define the native state among all possible states of a protein is a problem of immense interest. The present study explores the correlation between solvent accessible surface areas (SASAs) and different intraprotein as well as protein-water hydrogen bonds of 55 single-chain globular proteins from four different structural classes (all α, all β, α+β, and α/β), 16 multichain proteins, and 4 macromolecular complexes. A systematic analysis of the solvent accessible surface area and intraprotein and protein-water hydrogen bonds suggests a linear relationship between SASAs and hydrogen bonds. The number of protein-water hydrogen bonds per unit SASA ranges from 3 to 4 for all the different structural protein classes. In contrast, the number of intramolecular hydrogen bonds per unit SASA, including the mainchain-mainchain, mainchain-sidechain, and sidechain-sidechain, varies between 0.75 to 2. The solvation free energy of a protein linearly decreases with SASA. Our study also shows that the solvation free energy/SASA varies from -75 to -105 kJ mol-1 nm-2 across all the native states studied here. The number conservancy of intraprotein hydrogen bonds per unit SASA possibly imparts structural stability to the native structure. On the other hand, 3-4 protein-water hydrogen bonds per unit SASA are possibly required to maintain a balance between the solubility and functionality of the native states. This study provides a basis for synthetic biologists to design new folds with improved functionality.
Collapse
Affiliation(s)
- Prasun Pal
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Sandipan Chakraborty
- Center for Innovation in Molecular and Pharmaceutical Sciences (CIMPS), Dr. Reddy's Institute of Life Sciences, University of Hyderabad Campus, Gachibowli, Hyderabad 500046, India
| | - Biman Jana
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
3
|
Spicher S, Plett C, Pracht P, Hansen A, Grimme S. Automated Molecular Cluster Growing for Explicit Solvation by Efficient Force Field and Tight Binding Methods. J Chem Theory Comput 2022; 18:3174-3189. [PMID: 35482317 DOI: 10.1021/acs.jctc.2c00239] [Citation(s) in RCA: 35] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
An automated and broadly applicable workflow for the description of solvation effects in an explicit manner is introduced. This method, termed quantum cluster growth (QCG), is based on the semiempirical GFN2-xTB/GFN-FF methods, enabling efficient geometry optimizations and MD simulations. Fast structure generation is provided using the intermolecular force field xTB-IFF. Additionally, the approach uses an efficient implicit solvation model for the electrostatic embedding of the growing clusters. The novel QCG procedure presents a robust cluster generation tool for subsequent application of higher-level (e.g., DFT) methods to study solvation effects on molecular geometries explicitly or to average spectroscopic properties over cluster ensembles. Furthermore, the computation of the solvation free energy with a supermolecular approach can be carried out with QCG. The underlying growing process is physically motivated by computing the leading-order solute-solvent interactions first and can account for conformational and chemical changes due to solvation for low-energy barrier processes. The conformational space is explored with the NCI-MTD algorithm as implemented in the CREST program, using a combination of metadynamics and MD simulations. QCG with GFN2-xTB yields realistic solution geometries and reasonable solvation free energies for various systems without introducing many empirical parameters. Computed IR spectra of some solutes with QCG show a better match to the experimental data compared to well-established implicit solvation models.
Collapse
Affiliation(s)
- Sebastian Spicher
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Christoph Plett
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Philipp Pracht
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
4
|
Caceres-Delpiano J, Wang LP, Essex JW. The automated optimisation of a coarse-grained force field using free energy data. Phys Chem Chem Phys 2021; 23:24842-24851. [PMID: 34723311 PMCID: PMC8579472 DOI: 10.1039/d0cp05041e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2020] [Accepted: 10/18/2021] [Indexed: 11/21/2022]
Abstract
Atomistic models provide a detailed representation of molecular systems, but are sometimes inadequate for simulations of large systems over long timescales. Coarse-grained models enable accelerated simulations by reducing the number of degrees of freedom, at the cost of reduced accuracy. New optimisation processes to parameterise these models could improve their quality and range of applicability. We present an automated approach for the optimisation of coarse-grained force fields, by reproducing free energy data derived from atomistic molecular simulations. To illustrate the approach, we implemented hydration free energy gradients as a new target for force field optimisation in ForceBalance and applied it successfully to optimise the un-charged side-chains and the protein backbone in the SIRAH protein coarse-grain force field. The optimised parameters closely reproduced hydration free energies of atomistic models and gave improved agreement with experiment.
Collapse
Affiliation(s)
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, California 95616, USA.
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Southapton, S017 1BJ, UK.
| |
Collapse
|
5
|
Fortuna A, Costa PJ. Optimized Halogen Atomic Radii for PBSA Calculations Using Off-Center Point Charges. J Chem Inf Model 2021; 61:3361-3375. [PMID: 34185532 DOI: 10.1021/acs.jcim.1c00177] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In force-field methods, the usage of off-center point charges, also called extra points (EPs), is a common strategy to tackle the anisotropy of the electrostatic potential of covalently bonded halogens (X), thus allowing the description of halogen bonds (XBs) at the molecular mechanics/molecular dynamics (MM/MD) level. Diverse EP implementations exist in the literature differing on the charge sets and/or the X-EP distances. Poisson-Boltzmann and surface area (PBSA) calculations can be used to obtain solvation free energies (ΔGsolv) of small molecules, often to compute binding free energies (ΔGbind) at the MM-PBSA level. This method depends, among other parameters, on the empirical assignment of atomic radii (PB radii). Given the multiplicity of off-center point-charge models and the lack of specific PB radii for halogens compatible with such implementations, in this work, we assessed the performance of PBSA calculations for the estimation of ΔGsolv values in water (ΔGhyd), also conducting an optimization of the halogen PB radii (Cl, Br, and I) for each EP model. We not only expand the usage of EP models in the scope of the general AMBER force field (GAFF) but also provide the first optimized halogen PB radii in the context of the CHARMM general force field (CGenFF), thus contributing to improving the description of halogenated compounds in PBSA calculations.
Collapse
Affiliation(s)
- Andreia Fortuna
- BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, 1749-016 Lisboa, Portugal.,Research Institute for Medicines (iMed.ULisboa), Faculty of Pharmacy, University of Lisbon, Av. Professor Gama Pinto, 1649-003 Lisbon, Portugal
| | - Paulo J Costa
- BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, 1749-016 Lisboa, Portugal
| |
Collapse
|
6
|
Sengupta A, Li Z, Song LF, Li P, Merz KM. Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models. J Chem Inf Model 2021; 61:869-880. [PMID: 33538599 DOI: 10.1021/acs.jcim.0c01390] [Citation(s) in RCA: 82] [Impact Index Per Article: 27.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Monovalent ions play significant roles in various biological and material systems. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB), with significantly improved descriptions of condensed phase water, have been developed. The pairwise interaction between the metal ion and water necessitates the development of ion parameters specifically for these water models. Herein, we parameterized the 12-6 and the 12-6-4 nonbonded models for 12 monovalent ions with the respective four new water models. These monovalent ions contain eight cations including alkali metal ions (Li+, Na+, K+, Rb+, Cs+), transition-metal ions (Cu+ and Ag+), and Tl+ from the boron family, along with four halide anions (F-, Cl-, Br-, I-). Our parameters were designed to reproduce the target hydration free energies (the 12-6 hydration free energy (HFE) set), the ion-oxygen distances (the 12-6 ion-oxygen distance (IOD) set), or both of them (the 12-6-4 set). The 12-6-4 parameter set provides highly accurate structural features overcoming the limitations of the routinely used 12-6 nonbonded model for ions. Specifically, we note that the 12-6-4 parameter set is able to reproduce experimental hydration free energies within 1 kcal/mol and experimental ion-oxygen distances within 0.01 Å simultaneously. We further reproduced the experimentally determined activity derivatives for salt solutions, validating the ion parameters for simulations of ion pairs. The improved performance of the present water models over our previous parameter sets for the TIP3P, TIP4P, and SPC/E water models (Li, P. et al J. Chem. Theory Comput. 2015 11 1645 1657) highlights the importance of the choice of water model in conjunction with the metal ion parameter set.
Collapse
Affiliation(s)
- Arkajyoti Sengupta
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Zhen Li
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Lin Frank Song
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Pengfei Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States.,Department of Chemistry and Biochemistry, Loyola University Chicago, Chicago, Illinois 60660, United States
| | - Kenneth M Merz
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States.,Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
7
|
Liu B, Zhang W, Guo S, Zuo Z. Discovery of novel modulators targeting human TRPC5: Docking-based virtual screening, molecular dynamics simulation and binding affinity predication. J Mol Graph Model 2020; 102:107795. [PMID: 33161371 DOI: 10.1016/j.jmgm.2020.107795] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Revised: 10/28/2020] [Accepted: 10/28/2020] [Indexed: 11/20/2022]
Abstract
Canonical transient receptor potential channel 5 (TRPC5) plays a key role in the regulation of central nervous system, cardiovascular system, kidney disease, cancer, and could be also involved in liver function, arthritis, diabetes-associated complications and so on. However, evidence of TRPC5 function on cellular or organismic levels is sparse. There is still a need for identifying novel and efficient TRPC5 channel modulators to study TRPC5 function. In this study, based on the hTRPC5 structure obtained by homology modeling and the predicted binding site, we have performed virtual screening of 212,736 compounds from the specs database(http://www.specs.net) to find potential hTRPC5 modulators. Lipinski and Veber rules, ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) and PAINS (Pan Assay Interference structures) filters were used to screen the large database. Further, multi-software combination docking, cluster analysis and interaction analysis were used to select 20 potential active candidates with novel skeleton. 4 Hits, bearing appreciable binding affinity with hTRPC5 were selected for 40ns all-atom molecular dynamics (MD) simulations under explicit water conditions. The MD simulation results suggested that the 4 Hits binding induces a slight structural change and stabilizes the hTRPC5 structure. In addition, decomposition free energy demonstrated that residues TRP434, LEU437, MET438, ALA441, ILE484, ILE487, LEU488, LEU491, LEU515, ILE517, LEU518, LEU521, PHE531, THR607, VAL610, ILE611, VAL615 played the critical role on system stability. 4 Hits, as potential modulators of hTRPC5, may be potential leads to develop effective therapeutics hTRPC5-associated diseases.
Collapse
Affiliation(s)
- Bin Liu
- State Key Laboratory of Phytochemistry and Plant Resources in West China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, PR China; University of Chinese Academy of Sciences, Beijing 100049, PR China
| | - Wei Zhang
- State Key Laboratory of Phytochemistry and Plant Resources in West China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, PR China
| | - Sheng Guo
- State Key Laboratory of Phytochemistry and Plant Resources in West China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, PR China
| | - Zhili Zuo
- State Key Laboratory of Phytochemistry and Plant Resources in West China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, PR China; University of Chinese Academy of Sciences, Beijing 100049, PR China.
| |
Collapse
|
8
|
Wang B, Li C, Xiangyu J, Zhu T, Zhang JZH. An Approach to Computing Solvent Reorganization Energy. J Chem Theory Comput 2020; 16:6513-6519. [PMID: 32941037 DOI: 10.1021/acs.jctc.0c00672] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A method for accurate calculation of reorganization free energy of an explicit solvent solvating a solute molecule is presented. The method relies on the knowledge of solvation free energy known either from experimental measurement or from accurate calculation. An important part of this approach lies in the calculation of entropy in solute-solvent interaction free energy using the interaction entropy method, combined with the calculation of enthalpy of solvent reorganization based on a finite number of solvent molecules near the solute molecule. This interaction entropy-solvent reorganization or IESR method enables one to compute both enthalpy and entropy components of the solvent reorganization energy from MD simulation of the solvated system. The calculated results are determined by the molecular force field only without any empirical parameter. The current method is applied to computing the reorganization energies of water solvent solvating a wide range of solute molecules including both hydrophilic and hydrophobic ones. The accuracy of the approach is indirectly verified by the excellent agreement of individual enthalpies and entropies of the solvation energies between theory and experiment.
Collapse
Affiliation(s)
- Bo Wang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Cuiyu Li
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Jia Xiangyu
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Tong Zhu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States.,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan 030006, Shanxi, China
| |
Collapse
|
9
|
He X, Man VH, Yang W, Lee TS, Wang J. A fast and high-quality charge model for the next generation general AMBER force field. J Chem Phys 2020; 153:114502. [PMID: 32962378 DOI: 10.1063/5.0019056] [Citation(s) in RCA: 222] [Impact Index Per Article: 55.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The General AMBER Force Field (GAFF) has been broadly used by researchers all over the world to perform in silico simulations and modelings on diverse scientific topics, especially in the field of computer-aided drug design whose primary task is to accurately predict the affinity and selectivity of receptor-ligand binding. The atomic partial charges in GAFF and the second generation of GAFF (GAFF2) were originally developed with the quantum mechanics derived restrained electrostatic potential charge, but in practice, users usually adopt an efficient charge method, Austin Model 1-bond charge corrections (AM1-BCC), based on which, without expensive ab initio calculations, the atomic charges could be efficiently and conveniently obtained with the ANTECHAMBER module implemented in the AMBER software package. In this work, we developed a new set of BCC parameters specifically for GAFF2 using 442 neutral organic solutes covering diverse functional groups in aqueous solution. Compared to the original BCC parameter set, the new parameter set significantly reduced the mean unsigned error (MUE) of hydration free energies from 1.03 kcal/mol to 0.37 kcal/mol. More excitingly, this new AM1-BCC model also showed excellent performance in the solvation free energy (SFE) calculation on diverse solutes in various organic solvents across a range of different dielectric constants. In this large-scale test with totally 895 neutral organic solvent-solute systems, the new parameter set led to accurate SFE predictions with the MUE and the root-mean-square-error of 0.51 kcal/mol and 0.65 kcal/mol, respectively. This newly developed charge model, ABCG2, paved a promising path for the next generation GAFF development.
Collapse
Affiliation(s)
- Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Viet H Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Wei Yang
- Department of Chemistry and Biochemistry and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306, USA
| | - Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| |
Collapse
|
10
|
Schindler CEM, Baumann H, Blum A, Böse D, Buchstaller HP, Burgdorf L, Cappel D, Chekler E, Czodrowski P, Dorsch D, Eguida MKI, Follows B, Fuchß T, Grädler U, Gunera J, Johnson T, Jorand Lebrun C, Karra S, Klein M, Knehans T, Koetzner L, Krier M, Leiendecker M, Leuthner B, Li L, Mochalkin I, Musil D, Neagu C, Rippmann F, Schiemann K, Schulz R, Steinbrecher T, Tanzer EM, Unzue Lopez A, Viacava Follis A, Wegener A, Kuhn D. Large-Scale Assessment of Binding Free Energy Calculations in Active Drug Discovery Projects. J Chem Inf Model 2020; 60:5457-5474. [PMID: 32813975 DOI: 10.1021/acs.jcim.0c00900] [Citation(s) in RCA: 136] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Accurate ranking of compounds with regards to their binding affinity to a protein using computational methods is of great interest to pharmaceutical research. Physics-based free energy calculations are regarded as the most rigorous way to estimate binding affinity. In recent years, many retrospective studies carried out both in academia and industry have demonstrated its potential. Here, we present the results of large-scale prospective application of the FEP+ method in active drug discovery projects in an industry setting at Merck KGaA, Darmstadt, Germany. We compare these prospective data to results obtained on a new diverse, public benchmark of eight pharmaceutically relevant targets. Our results offer insights into the challenges faced when using free energy calculations in real-life drug discovery projects and identify limitations that could be tackled by future method development. The new public data set we provide to the community can support further method development and comparative benchmarking of free energy calculations.
Collapse
Affiliation(s)
| | - Hannah Baumann
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Andreas Blum
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Dietrich Böse
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | | | - Lars Burgdorf
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | | | - Eugene Chekler
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Paul Czodrowski
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Dieter Dorsch
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | | | - Bruce Follows
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Thomas Fuchß
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Ulrich Grädler
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Jakub Gunera
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Theresa Johnson
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Catherine Jorand Lebrun
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Srinivasa Karra
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Markus Klein
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Tim Knehans
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Lisa Koetzner
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Mireille Krier
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | | | | | - Liwei Li
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Igor Mochalkin
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Djordje Musil
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Constantin Neagu
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | | | - Kai Schiemann
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Robert Schulz
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany.,Institute of Pharmacy, Freie Universität Berlin, Königin-Luise-Straße 2+4, 14195 Berlin, Germany
| | | | - Eva-Maria Tanzer
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | | | - Ariele Viacava Follis
- EMD Serono Research and Development Institute Inc., 45A Middlesex Turnpike, Billerica, Massachusetts 01821, United States
| | - Ansgar Wegener
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| | - Daniel Kuhn
- Merck KGaA, Frankfurter Strasse 250, 64293 Darmstadt, Germany
| |
Collapse
|
11
|
Scheen J, Wu W, Mey ASJS, Tosco P, Mackey M, Michel J. Hybrid Alchemical Free Energy/Machine-Learning Methodology for the Computation of Hydration Free Energies. J Chem Inf Model 2020; 60:5331-5339. [PMID: 32639733 DOI: 10.1021/acs.jcim.0c00600] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
A methodology that combines alchemical free energy calculations (FEP) with machine learning (ML) has been developed to compute accurate absolute hydration free energies. The hybrid FEP/ML methodology was trained on a subset of the FreeSolv database and retrospectively shown to outperform most submissions from the SAMPL4 competition. Compared to pure machine-learning approaches, FEP/ML yields more precise estimates of free energies of hydration and requires a fraction of the training set size to outperform standalone FEP calculations. The ML-derived correction terms are further shown to be transferable to a range of related FEP simulation protocols. The approach may be used to inexpensively improve the accuracy of FEP calculations and to flag molecules which will benefit the most from bespoke force field parametrization efforts.
Collapse
Affiliation(s)
- Jenke Scheen
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Wilson Wu
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Antonia S J S Mey
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| | - Paolo Tosco
- Cresset Group, New Cambridge House, Bassingbourn Road, Litlington, Cambridgeshire SG8 0SS, United Kingdom
| | - Mark Mackey
- Cresset Group, New Cambridge House, Bassingbourn Road, Litlington, Cambridgeshire SG8 0SS, United Kingdom
| | - Julien Michel
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| |
Collapse
|
12
|
Rasouli S, Moghbeli MR, Nikkhah SJ. Molecular dynamics simulation of polystyrene copolymer with octyl short-chain branches in toluene. J Mol Model 2020; 26:80. [PMID: 32180007 DOI: 10.1007/s00894-020-4339-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Accepted: 03/04/2020] [Indexed: 11/28/2022]
Abstract
In this study, dimensional, conformational and dynamic behaviors of a short-chain branched styrene/1-octene copolymer chain with different 1-octene percentages, i.e., 0, 2, 4 and 6%, in toluene are investigated at the temperature of 298.15 K via molecular dynamics simulation. The chain dimensions and flexibility in the solvent are evaluated by calculating the radius of gyration (Rg), end-to-end distance (<r0>), surface area (Ach), and volume (Vch) of the copolymer chain. The mean square displacement (MSD) and diffusivity coefficient for each copolymer chain are measured to determine its dynamic behavior and mobility in aromatic media. To consider the effect of increasing the 1-octene co-monomer percentage on the copolymer chain affinity to the solvent molecules, the interaction energy (Eint) and Flory-Huggins (FH) interaction parameter are calculated for each equilibrated solution model. The simulation results indicate that the co-monomer level increment in the copolymer structure reduces the chain Rg amount and its interaction with the solvent. The <r0> of the chain increases up to 4% co-monomer content, while further co-monomer content decreases the <r0> value. Additionally, the viscosity of the equilibrated dilute solutions is calculated via non-equilibrium molecular dynamics simulation (NEMD). Moreover, the steric hindrance of the copolymers and the solvent molecules capturing in the dilute solution is determined via radial distribution function (RDF) analysis. Helmholtz free energy and the system entropy changes are calculated to evaluate the tendency of the copolymer to the solvent molecules and its dilute solution irregularity, respectively. Graphical abstract The figure shows the variations trend of the poly(styrene-co-1-octene) chain dimensions in toluene aromatic solvent by increasing the 1-octene content (x), after the equilibration state. Red and blue colors represent the carbon atoms of the copolymer chain backbone and 1-octene side chains, respectively. The styrene rings and the hydrogen atoms of the chains were removed for better view.
Collapse
Affiliation(s)
- Sajad Rasouli
- Smart Polymers and Nanocomposites Research Group, School of Chemical Engineering, Iran University of Science and Technology, Tehran, 16846-13114, Iran
| | - Mohammad Reza Moghbeli
- Smart Polymers and Nanocomposites Research Group, School of Chemical Engineering, Iran University of Science and Technology, Tehran, 16846-13114, Iran.
| | - Sousa Javan Nikkhah
- Smart Polymers and Nanocomposites Research Group, School of Chemical Engineering, Iran University of Science and Technology, Tehran, 16846-13114, Iran
| |
Collapse
|
13
|
Abstract
We employed molecular dynamics simulations on the water solvation of conically shaped carbon nanoparticles. We explored the hydrophobic behaviour of the nanoparticles and investigated microscopically the cavitation of water in a conical confinement with different angles. We performed additional molecular dynamics simulations in which the carbon structures do not interact with water as if they were in vacuum. We detected a waving on the surface of the cones that resembles the shape agitations of artificial water channels and biological porins. The surface waves were induced by the pentagonal carbon rings (in an otherwise hexagonal network of carbon rings) concentrated near the apex of the cones. The waves were affected by the curvature gradients on the surface. They were almost undetected for the case of an armchair nanotube. Understanding such nanoscale phenomena is the key to better designed molecular models for membrane systems and nanodevices for energy applications and separation.
Collapse
|
14
|
Nunes R, Vila-Viçosa D, Costa PJ. Tackling Halogenated Species with PBSA: Effect of Emulating the σ-Hole. J Chem Theory Comput 2019; 15:4241-4251. [PMID: 31142112 DOI: 10.1021/acs.jctc.9b00106] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
To model halogen-bond phenomena using classical force fields, an extra point (EP) of charge is frequently introduced at a given distance from the halogen (X) to emulate the σ-hole. The resulting molecular dynamics (MD) trajectories can be used in subsequent molecular mechanics (MM) combined with Poisson-Boltzmann and surface area calculations (PBSA) to estimate protein-ligand binding free energies (Δ Gbind). While EP addition improves the MM/MD description of halogen-containing systems, its effect on the calculation of solvation free energies (Δ Gsolv) using the PBSA approach is yet to be assessed. As the PBSA calculations depend, among other parameters, on the empirical assignment of radii (PB radii), a problematic issue arises, since standard halogen radii are smaller than the typical X···EP distances, thus placing the EP within the solvent dielectric. Herein, we took a common literature EP parametrization scheme, which uses X···EP = Rmin and RESP charges in the context of GAFF, and performed a comprehensive study on the performance of PBSA (using three different setups) in the calculation of Δ Gsolv values for 142 halogenated compounds (bearing Cl, Br, or I) for which the experimental values are known. By conducting an optimization (minimizing the error against experimental values), we provide a new optimized set of halogen PB radii, for each PBSA setup, that should be used in the context of the aforementioned scenario. A simultaneous optimization of PB radii and X···EP distances shows that a wide range of distance/radius pairs can be used without significant loss of accuracy, therefore laying the basis for expanding this halogen radii optimization strategy to other force fields and EP implementations. As ligand Δ Gsolv estimation is an important term in the determination of protein-ligand Δ Gbind, this work is particularly relevant in the framework of structure-based virtual screening and related computer-aided drug design routines.
Collapse
|
15
|
Nguyen DD, Wei GW. DG-GL: Differential geometry-based geometric learning of molecular datasets. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3179. [PMID: 30693661 PMCID: PMC6598676 DOI: 10.1002/cnm.3179] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2018] [Revised: 11/21/2018] [Accepted: 12/06/2018] [Indexed: 05/11/2023]
Abstract
MOTIVATION Despite its great success in various physical modeling, differential geometry (DG) has rarely been devised as a versatile tool for analyzing large, diverse, and complex molecular and biomolecular datasets because of the limited understanding of its potential power in dimensionality reduction and its ability to encode essential chemical and biological information in differentiable manifolds. RESULTS We put forward a differential geometry-based geometric learning (DG-GL) hypothesis that the intrinsic physics of three-dimensional (3D) molecular structures lies on a family of low-dimensional manifolds embedded in a high-dimensional data space. We encode crucial chemical, physical, and biological information into 2D element interactive manifolds, extracted from a high-dimensional structural data space via a multiscale discrete-to-continuum mapping using differentiable density estimators. Differential geometry apparatuses are utilized to construct element interactive curvatures in analytical forms for certain analytically differentiable density estimators. These low-dimensional differential geometry representations are paired with a robust machine learning algorithm to showcase their descriptive and predictive powers for large, diverse, and complex molecular and biomolecular datasets. Extensive numerical experiments are carried out to demonstrate that the proposed DG-GL strategy outperforms other advanced methods in the predictions of drug discovery-related protein-ligand binding affinity, drug toxicity, and molecular solvation free energy. AVAILABILITY AND IMPLEMENTATION http://weilab.math.msu.edu/DG-GL/ Contact: wei@math.msu.edu.
Collapse
Affiliation(s)
- Duc Duy Nguyen
- Department of Mathematics, Michigan State University, East Lansing, 48824, Michigan
| | - Guo-Wei Wei
- Department of Mathematics, Michigan State University, East Lansing, 48824, Michigan
- Department of Electrical and Computer Engineering, Michigan State University, MI 48824, USA
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, Michigan
| |
Collapse
|
16
|
Jia X, Li P. Solvation Free Energy Calculation Using a Fixed-Charge Model: Implicit and Explicit Treatments of the Polarization Effect. J Phys Chem B 2019; 123:1139-1148. [PMID: 30628452 DOI: 10.1021/acs.jpcb.8b10479] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, IPolQ-Mod charges and the reference potential scheme are used to calculate the solvation free energies of a set of organic molecules. Both methods could capture the phase transfer of a solute with accompanying polarization cost utilizing a fixed-charge model. The IPolQ-Mod charges, which are the average of two charge sets fitted in a vacuum state and a condensed phase, take account of the polarization effect implicitly. For the reference potential method, the quantum mechanics polarization corrections are calculated explicitly by thermodynamic perturbation. The polarization effect captured by the IPolQ-Mod charges is an approximation to that of the reference potential method theoretically. In the present study, the reference potential method shows a slight improvement over the classical restrained electrostatic potential (RESP) charges, which perform pretty well in predicting the solvation free energy. However, IPolQ-Mod(MP2) shows a poor agreement with the experimental data. Compared with IPolQ-Mod(MP2), IPolQ-Mod(M06-2X) or IPolQ-Mod(ωB97X) is found to give more appropriate prediction of the molecule's dipole and the solvation free energies calculated by IPolQ-Mod(M06-2X) or IPolQ-Mod(ωB97X) are more compatible with those of the RESP charges. If the other force field parameters remain unchanged, M06-2X or ωB97X is recommended to derive the IPolQ-Mod charges.
Collapse
Affiliation(s)
- Xiangyu Jia
- NYU Shanghai , 1555 Century Avenue , Shanghai 200122 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , 3663 Zhongshan Road North , Shanghai 200127 , China
| | - Pengfei Li
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science , East China Normal University , Shanghai 200062 , China
| |
Collapse
|
17
|
Monte MJS, Almeida ARRP. Estimations of the thermodynamic properties of halogenated benzenes as they relate to their environment mobility. CHEMOSPHERE 2017; 189:590-598. [PMID: 28963976 DOI: 10.1016/j.chemosphere.2017.09.095] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2017] [Revised: 09/19/2017] [Accepted: 09/20/2017] [Indexed: 06/07/2023]
Abstract
In this work, several simple new equations for predicting important environmental mobility properties, at T = 298.15 K, were derived for halogenated benzenes: standard Gibbs energy of hydration, aqueous solubility, octanol-water partition coefficients, and Henry's law constants. A discussion on our previous estimates of other related properties (standard Gibbs energy and vapor pressure of sublimation and of vaporization) and their relation with entropy of fusion is also presented. As we aimed to estimate these properties for any of the ca. 1500 halogenated benzenes that may exist theoretically, an equation for estimating the temperature of fusion was also derived, since some of the proposed predictive equations (solubility of solids and Gibbs energy of sublimation) require its knowledge. For the other estimated properties just the number of each halogen that replaces hydrogen atoms in the halogenated benzene is needed. It was found that the coefficients that multiply the number of halogen atoms in the predictive equations vary linearly with the volume of the halogen atom.
Collapse
Affiliation(s)
- Manuel J S Monte
- Centro de Investigação em Química (CIQUP), Department of Chemistry and Biochemistry, Faculty of Science, University of Porto, Rua do Campo Alegre, 687, P-4169-007 Porto, Portugal.
| | - Ana R R P Almeida
- Centro de Investigação em Química (CIQUP), Department of Chemistry and Biochemistry, Faculty of Science, University of Porto, Rua do Campo Alegre, 687, P-4169-007 Porto, Portugal
| |
Collapse
|
18
|
Wang B, Wang C, Wu K, Wei G. Breaking the polar‐nonpolar division in solvation free energy prediction. J Comput Chem 2017; 39:217-233. [DOI: 10.1002/jcc.25107] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2017] [Revised: 09/02/2017] [Accepted: 10/22/2017] [Indexed: 11/12/2022]
Affiliation(s)
- Bao Wang
- Department of MathematicsMichigan State University Michigan48824
| | - Chengzhang Wang
- School of Statistics and MathematicsCentral University of Finance and EconomicsBeijing100081 China
| | - Kedi Wu
- Department of MathematicsMichigan State University Michigan48824
| | - Guo‐Wei Wei
- Department of MathematicsMichigan State University Michigan48824
- Department of Electrical and ComputerEngineering Michigan State University Michigan48824
- Department of Biochemistry and MolecularBiology Michigan State UniversityMichigan48824
| |
Collapse
|
19
|
Wang M, Li P, Jia X, Liu W, Shao Y, Hu W, Zheng J, Brooks BR, Mei Y. Efficient Strategy for the Calculation of Solvation Free Energies in Water and Chloroform at the Quantum Mechanical/Molecular Mechanical Level. J Chem Inf Model 2017; 57:2476-2489. [DOI: 10.1021/acs.jcim.7b00001] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Pengfei Li
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Xiangyu Jia
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Wei Liu
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman Oklahoma 73019, United States
| | - Wenxin Hu
- The Computer
Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer
Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Ye Mei
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
20
|
Mohamed NA, Bradshaw RT, Essex JW. Evaluation of solvation free energies for small molecules with the AMOEBA polarizable force field. J Comput Chem 2016; 37:2749-2758. [PMID: 27757978 PMCID: PMC5111595 DOI: 10.1002/jcc.24500] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Revised: 09/05/2016] [Accepted: 09/07/2016] [Indexed: 01/24/2023]
Abstract
The effects of electronic polarization in biomolecular interactions will differ depending on the local dielectric constant of the environment, such as in solvent, DNA, proteins, and membranes. Here the performance of the AMOEBA polarizable force field is evaluated under nonaqueous conditions by calculating the solvation free energies of small molecules in four common organic solvents. Results are compared with experimental data and equivalent simulations performed with the GAFF pairwise-additive force field. Although AMOEBA results give mean errors close to "chemical accuracy," GAFF performs surprisingly well, with statistically significantly more accurate results than AMOEBA in some solvents. However, for both models, free energies calculated in chloroform show worst agreement to experiment and individual solutes are consistently poor performers, suggesting non-potential-specific errors also contribute to inaccuracy. Scope for the improvement of both potentials remains limited by the lack of high quality experimental data across multiple solvents, particularly those of high dielectric constant. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Noor Asidah Mohamed
- Computational Systems Chemistry, School of ChemistryUniversity of SouthamptonHighfieldSouthamptonSO17 1BJUK
| | - Richard T. Bradshaw
- Computational Systems Chemistry, School of ChemistryUniversity of SouthamptonHighfieldSouthamptonSO17 1BJUK
| | - Jonathan W. Essex
- Computational Systems Chemistry, School of ChemistryUniversity of SouthamptonHighfieldSouthamptonSO17 1BJUK
| |
Collapse
|
21
|
Bradshaw RT, Essex JW. Evaluating Parametrization Protocols for Hydration Free Energy Calculations with the AMOEBA Polarizable Force Field. J Chem Theory Comput 2016; 12:3871-83. [PMID: 27341007 DOI: 10.1021/acs.jctc.6b00276] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hydration free energy (HFE) calculations are often used to assess the performance of biomolecular force fields and the quality of assigned parameters. The AMOEBA polarizable force field moves beyond traditional pairwise additive models of electrostatics and may be expected to improve upon predictions of thermodynamic quantities such as HFEs over and above fixed-point-charge models. The recent SAMPL4 challenge evaluated the AMOEBA polarizable force field in this regard but showed substantially worse results than those using the fixed-point-charge GAFF model. Starting with a set of automatically generated AMOEBA parameters for the SAMPL4 data set, we evaluate the cumulative effects of a series of incremental improvements in parametrization protocol, including both solute and solvent model changes. Ultimately, the optimized AMOEBA parameters give a set of results that are not statistically significantly different from those of GAFF in terms of signed and unsigned error metrics. This allows us to propose a number of guidelines for new molecule parameter derivation with AMOEBA, which we expect to have benefits for a range of biomolecular simulation applications such as protein-ligand binding studies.
Collapse
Affiliation(s)
- Richard T Bradshaw
- School of Chemistry, University of Southampton, Highfield Campus , Southampton SO17 1BJ, U.K
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Highfield Campus , Southampton SO17 1BJ, U.K
| |
Collapse
|
22
|
Mih N, Brunk E, Bordbar A, Palsson BO. A Multi-scale Computational Platform to Mechanistically Assess the Effect of Genetic Variation on Drug Responses in Human Erythrocyte Metabolism. PLoS Comput Biol 2016; 12:e1005039. [PMID: 27467583 PMCID: PMC4965186 DOI: 10.1371/journal.pcbi.1005039] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2016] [Accepted: 06/27/2016] [Indexed: 12/31/2022] Open
Abstract
Progress in systems medicine brings promise to addressing patient heterogeneity and individualized therapies. Recently, genome-scale models of metabolism have been shown to provide insight into the mechanistic link between drug therapies and systems-level off-target effects while being expanded to explicitly include the three-dimensional structure of proteins. The integration of these molecular-level details, such as the physical, structural, and dynamical properties of proteins, notably expands the computational description of biochemical network-level properties and the possibility of understanding and predicting whole cell phenotypes. In this study, we present a multi-scale modeling framework that describes biological processes which range in scale from atomistic details to an entire metabolic network. Using this approach, we can understand how genetic variation, which impacts the structure and reactivity of a protein, influences both native and drug-induced metabolic states. As a proof-of-concept, we study three enzymes (catechol-O-methyltransferase, glucose-6-phosphate dehydrogenase, and glyceraldehyde-3-phosphate dehydrogenase) and their respective genetic variants which have clinically relevant associations. Using all-atom molecular dynamic simulations enables the sampling of long timescale conformational dynamics of the proteins (and their mutant variants) in complex with their respective native metabolites or drug molecules. We find that changes in a protein's structure due to a mutation influences protein binding affinity to metabolites and/or drug molecules, and inflicts large-scale changes in metabolism.
Collapse
Affiliation(s)
- Nathan Mih
- Bioinformatics and Systems Biology Graduate Program, University of California, San Diego, La Jolla, California, United States of America
| | - Elizabeth Brunk
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
- * E-mail: (EB); (BOP)
| | - Aarash Bordbar
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
| | - Bernhard O. Palsson
- Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America
- Department of Pediatrics, University of California, San Diego, La Jolla, California, United States of America
- * E-mail: (EB); (BOP)
| |
Collapse
|
23
|
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
|
24
|
El Hage K, Bereau T, Jakobsen S, Meuwly M. Impact of Quadrupolar Electrostatics on Atoms Adjacent to the Sigma-Hole in Condensed-Phase Simulations. J Chem Theory Comput 2016; 12:3008-19. [PMID: 27158892 DOI: 10.1021/acs.jctc.6b00202] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Halogenation is one of the cases for which advanced molecular simulation methods are mandatory for quantitative and predictive studies. The present work provides a systematic investigation of the importance of higher-order multipoles on specific sites of halobenzenes, other than the halogen, for static and dynamic properties in condensed-phase simulations. For that purpose, solute-solvent interactions using point charge (PC), multipole (MTP), and hybrid point charge/multipole (HYB) electrostatic models are analyzed in regions of halogen bonding and extended to regions of π orbitals of phenyl carbons. Using molecular dynamics simulations and quantum chemical methods, it is found that the sigma-hole does not only affect the halogen and the carbon bound to it but its effect extends to the carbons adjacent to the CX bond. This effect increases with the magnitude of the positive potential of the sigma-hole. With the MTP and HYB3 models, all hydration free energies of the PhX compounds are reproduced within 0.1 kcal/mol. Analysis of pair distribution functions and hydration free energies of halogenated benzenes provides a microscopic explanation why "point charge"-based representations with off-site charges fail in reproducing thermodynamic properties of the sigma-hole. Application of the hybrid models to study protein-ligand binding demonstrates both their accuracy and computational efficiency.
Collapse
Affiliation(s)
- Krystel El Hage
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Tristan Bereau
- Max Planck Institute for Polymer Research , Ackermannweg 10, 55128 Mainz, Germany
| | - Sofie Jakobsen
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, 4056 Basel, Switzerland.,Department of Chemistry, Aarhus University , Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Markus Meuwly
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
25
|
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
|