1
|
Ishizuka R. Martini 3D-OZ: A Theoretical Investigation of Solvation Shell Structures and Solvation Free Energies of Martini Coarse-Grained Proteins. J Phys Chem B 2024; 128:8522-8529. [PMID: 39180742 DOI: 10.1021/acs.jpcb.4c03897] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/26/2024]
Abstract
We investigate the properties of aqueous solutions using integral equation theories and molecular dynamics (MD) simulations within the framework of the MARTINI coarse-grained force field. The integral equation theory used in the present work is based on the Ornstein-Zernike equation coupled with the hypernetted chain (HNC) and Kovalenko-Hirata (KH) closures. Overall, the solvation shell structures and solvation thermodynamics in the HNC approximation are shown to be in better agreement with those from the MD simulation than the KH results. Especially, through the analysis of spatial distribution functions of water around a protein, we have demonstrated that the HNC approximation can provide the highly anisotropic structure of the solvation shell of the protein. On the other hand, the KH approximation works well for simple particle solutes, but the results for highly hydrated proteins deviate quite significantly from the MD results. We further explore in detail the reason underlying the deviation caused by the KH approximation. Lastly, a potential application of the integral equation theory with the MARTINI model is outlined.
Collapse
Affiliation(s)
- Ryosuke Ishizuka
- Zkanics F.P.S., Side-6 Senriyama-West, Suita, Osaka 565-0851, Japan
| |
Collapse
|
2
|
Subramanian V, Ratkova E, Palmer D, Engkvist O, Fedorov M, Llinas A. Multisolvent Models for Solvation Free Energy Predictions Using 3D-RISM Hydration Thermodynamic Descriptors. J Chem Inf Model 2020; 60:2977-2988. [PMID: 32311268 DOI: 10.1021/acs.jcim.0c00065] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The potential to predict solvation free energies (SFEs) in any solvent using a machine learning (ML) model based on thermodynamic output, extracted exclusively from 3D-RISM simulations in water is investigated. The models on multiple solvents take into account both the solute and solvent description and offer the possibility to predict SFEs of any solute in any solvent with root mean squared errors less than 1 kcal/mol. Validations that involve exclusion of fractions or clusters of the solutes or solvents exemplify the model's capability to predict SFEs of novel solutes and solvents with diverse chemical profiles. In addition to being predictive, our models can identify the solute and solvent features that influence SFE predictions. Furthermore, using 3D-RISM hydration thermodynamic output to predict SFEs in any organic solvent reduces the need to run 3D-RISM simulations in all these solvents. Altogether, our multisolvent models for SFE predictions that take advantage of the solvation effects are expected to have an impact in the property prediction space.
Collapse
Affiliation(s)
- Vigneshwari Subramanian
- Drug Metabolism and Pharmacokinetics, Research and Early Development-Respiratory, Inflammation and Autoimmune, Biopharmaceuticals R&D, AstraZeneca, Pepparedsleden 1, SE-431 83, Mölndal, Sweden.,Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, U.K
| | - Ekaterina Ratkova
- Medicinal Chemistry, Research and Early Development - Cardiovascular, Renal and Metabolism, Biopharmaceuticals R&D, AstraZeneca, Pepparedsleden 1, SE-431 83, Mölndal, Sweden
| | - David Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, U.K
| | - Ola Engkvist
- Hit Discovery, Discovery Sciences, R&D, AstraZeneca, Pepparedsleden 1, SE-431 83, Mölndal, Sweden
| | - Maxim Fedorov
- Center for Computational and Data-Intensive Science and Engineering, Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Moscow, 143026, Russia.,Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde, John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, U.K
| | - Antonio Llinas
- Drug Metabolism and Pharmacokinetics, Research and Early Development-Respiratory, Inflammation and Autoimmune, Biopharmaceuticals R&D, AstraZeneca, Pepparedsleden 1, SE-431 83, Mölndal, Sweden
| |
Collapse
|
3
|
Maruyama Y, Takano H, Mitsutake A. Analysis of molecular dynamics simulations of 10-residue peptide, chignolin, using statistical mechanics: Relaxation mode analysis and three-dimensional reference interaction site model theory. Biophys Physicobiol 2019; 16:407-429. [PMID: 31984194 PMCID: PMC6975981 DOI: 10.2142/biophysico.16.0_407] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Accepted: 08/29/2019] [Indexed: 01/03/2023] Open
Abstract
Molecular dynamics simulation is a fruitful tool for investigating the structural stability, dynamics, and functions of biopolymers at an atomic level. In recent years, simulations can be performed on time scales of the order of milliseconds using special purpose systems. Since the most stable structure, as well as meta-stable structures and intermediate structures, is included in trajectories in long simulations, it is necessary to develop analysis methods for extracting them from trajectories of simulations. For these structures, methods for evaluating the stabilities, including the solvent effect, are also needed. We have developed relaxation mode analysis to investigate dynamics and kinetics of simulations based on statistical mechanics. We have also applied the three-dimensional reference interaction site model theory to investigate stabilities with solvent effects. In this paper, we review the results for designing amino-acid substitution of the 10-residue peptide, chignolin, to stabilize the misfolded structure using these developed analysis methods.
Collapse
Affiliation(s)
- Yutaka Maruyama
- Architecture Development Team, FLAGSHIP 2020 Project, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Hiroshi Takano
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Kanagawa 223-8522, Japan
| | - Ayori Mitsutake
- Department of Physics, School of Science and Technology, Meiji University, Kawasaki, Kanagawa 214-8571, Japan
| |
Collapse
|
4
|
Nguyen C, Yamazaki T, Kovalenko A, Case DA, Gilson MK, Kurtzman T, Luchko T. A molecular reconstruction approach to site-based 3D-RISM and comparison to GIST hydration thermodynamic maps in an enzyme active site. PLoS One 2019; 14:e0219473. [PMID: 31291328 PMCID: PMC6619770 DOI: 10.1371/journal.pone.0219473] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2019] [Accepted: 06/24/2019] [Indexed: 11/25/2022] Open
Abstract
Computed, high-resolution, spatial distributions of solvation energy and entropy can provide detailed information about the role of water in molecular recognition. While grid inhomogeneous solvation theory (GIST) provides rigorous, detailed thermodynamic information from explicit solvent molecular dynamics simulations, recent developments in the 3D reference interaction site model (3D-RISM) theory allow many of the same quantities to be calculated in a fraction of the time. However, 3D-RISM produces atomic-site, rather than molecular, density distributions, which are difficult to extract physical meaning from. To overcome this difficulty, we introduce a method to reconstruct molecular density distributions from atomic-site density distributions. Furthermore, we assess the quality of the resulting solvation thermodynamics density distributions by analyzing the binding site of coagulation Factor Xa with both GIST and 3D-RISM. We find good qualitative agreement between the methods for oxygen and hydrogen densities as well as direct solute-solvent energetic interactions. However, 3D-RISM predicts lower energetic and entropic penalties for moving water from the bulk to the binding site.
Collapse
Affiliation(s)
- Crystal Nguyen
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California, United States of America
| | | | - Andriy Kovalenko
- National Institute for Nanotechnology, National Research Council of Canada, Edmonton, Alberta, Canada
- Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - David A. Case
- Department of Chemistry, Lehman College, The City University of New York, Bronx, New York, United States of America
| | - Michael K. Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California, United States of America
| | - Tom Kurtzman
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, United States of America
| | - Tyler Luchko
- Department of Physics and Astronomy, California State University, Northridge, California, United States of America
| |
Collapse
|
5
|
Cao S, Konovalov KA, Unarta IC, Huang X. Recent Developments in Integral Equation Theory for Solvation to Treat Density Inhomogeneity at Solute–Solvent Interface. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201900049] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Siqin Cao
- Department of Chemistrythe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
| | - Kirill A. Konovalov
- Department of Chemistrythe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
| | - Ilona Christy Unarta
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
- Bioengineering Graduate Programthe Hong Kong University of Science and TechnologyHong Kong of Chinese National EngineeringResearch Center for Tissue Restoration and Reconstructionthe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
| | - Xuhui Huang
- Department of Chemistrythe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
- Bioengineering Graduate Programthe Hong Kong University of Science and TechnologyHong Kong of Chinese National EngineeringResearch Center for Tissue Restoration and Reconstructionthe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- HKUST‐Shenzhen Research Institute Hi‐Tech Park, Nanshan Shenzhen 518057 China
| |
Collapse
|
6
|
Yildirim A, Wassenaar TA, van der Spoel D. Statistical efficiency of methods for computing free energy of hydration. J Chem Phys 2018; 149:144111. [DOI: 10.1063/1.5041835] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Ahmet Yildirim
- Department of Physics, Siirt University, Siirt 56100, Turkey
| | - Tsjerk A. Wassenaar
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands
| | - David van der Spoel
- Uppsala Centre for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala,
Sweden
| |
Collapse
|
7
|
Ansari SM, Palmer DS. Comparative Molecular Field Analysis Using Molecular Integral Equation Theory. J Chem Inf Model 2018; 58:1253-1265. [DOI: 10.1021/acs.jcim.7b00600] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Samiul M. Ansari
- Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, U.K
| | - David S. Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, U.K
| |
Collapse
|
8
|
Ansari SM, Sørensen J, Schiøtt B, Palmer DS. On the effect of mutations in bovine or camel chymosin on the thermodynamics of binding κ-caseins. Proteins 2018; 86:75-87. [DOI: 10.1002/prot.25410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2017] [Revised: 10/18/2017] [Accepted: 10/27/2017] [Indexed: 11/12/2022]
Affiliation(s)
- Samiul M. Ansari
- Department of Pure and Applied Chemistry; University of Strathclyde, Thomas Graham Building, 295 Cathedral Street; Glasgow G1 1XL Scotland
| | - Jesper Sørensen
- Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry; University of Aarhus, Langelandsgade 140; Aarhus DK 8000 Denmark
| | - Birgit Schiøtt
- Interdisciplinary Nanoscience Center (iNANO) and Department of Chemistry; University of Aarhus, Langelandsgade 140; Aarhus DK 8000 Denmark
| | - David S. Palmer
- Department of Pure and Applied Chemistry; University of Strathclyde, Thomas Graham Building, 295 Cathedral Street; Glasgow G1 1XL Scotland
| |
Collapse
|
9
|
Sumi T, Maruyama Y, Mitsutake A, Mochizuki K, Koga K. Application of reference‐modified density functional theory: Temperature and pressure dependences of solvation free energy. J Comput Chem 2017; 39:202-217. [DOI: 10.1002/jcc.25101] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Revised: 10/18/2017] [Accepted: 10/19/2017] [Indexed: 01/04/2023]
Affiliation(s)
- Tomonari Sumi
- Division of Superconducting and Functional MaterialsResearch Institute for Interdisciplinary Science, Okayama University, 3‐1‐1 Tsushima‐Naka, Kita‐kuOkayama700‐8530 Japan
- Department of Chemistry, Faculty of ScienceOkayama University, 3‐1‐1 Tsushima‐Naka, Kita‐kuOkayama700‐8530 Japan
| | - Yutaka Maruyama
- Co‐Design Team, FLAGSHIP 2020 Project, RIKEN Advanced Institute for Computational Science, 7‐1‐26, Minatojima‐minami‐machiKobe650‐0047 Japan
| | - Ayori Mitsutake
- Department of PhysicsKeio University, 3‐14‐1 Hiyoshi, Kohoku‐kuYokohama Kanagawa223–8522 Japan
| | - Kenji Mochizuki
- Division of Superconducting and Functional MaterialsResearch Institute for Interdisciplinary Science, Okayama University, 3‐1‐1 Tsushima‐Naka, Kita‐kuOkayama700‐8530 Japan
| | - Kenichiro Koga
- Division of Superconducting and Functional MaterialsResearch Institute for Interdisciplinary Science, Okayama University, 3‐1‐1 Tsushima‐Naka, Kita‐kuOkayama700‐8530 Japan
- Department of Chemistry, Faculty of ScienceOkayama University, 3‐1‐1 Tsushima‐Naka, Kita‐kuOkayama700‐8530 Japan
| |
Collapse
|
10
|
Ding L, Levesque M, Borgis D, Belloni L. Efficient molecular density functional theory using generalized spherical harmonics expansions. J Chem Phys 2017; 147:094107. [DOI: 10.1063/1.4994281] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Lu Ding
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Maximilien Levesque
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ. Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, École Normale Supérieure, CNRS, Processus d’Activation Sélective par Transfert d’Énergie Uni-Électronique ou Radiatif (PASTEUR), 75005 Paris, France
| | - Daniel Borgis
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ. Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, École Normale Supérieure, CNRS, Processus d’Activation Sélective par Transfert d’Énergie Uni-Électronique ou Radiatif (PASTEUR), 75005 Paris, France
| | - Luc Belloni
- LIONS, NIMBE, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
11
|
Effects of the crystal structure and thermodynamic stability on solubility of bioactive compounds: DFT study of isoniazid cocrystals. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.07.022] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
12
|
Sergiievskyi V, Jeanmairet G, Levesque M, Borgis D. Solvation free-energy pressure corrections in the three dimensional reference interaction site model. J Chem Phys 2015; 143:184116. [DOI: 10.1063/1.4935065] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
|
13
|
Palmer DS, Mišin M, Fedorov MV, Llinas A. Fast and General Method To Predict the Physicochemical Properties of Druglike Molecules Using the Integral Equation Theory of Molecular Liquids. Mol Pharm 2015. [PMID: 26212723 DOI: 10.1021/acs.molpharmaceut.5b00441] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report a method to predict physicochemical properties of druglike molecules using a classical statistical mechanics based solvent model combined with machine learning. The RISM-MOL-INF method introduced here provides an accurate technique to characterize solvation and desolvation processes based on solute-solvent correlation functions computed by the 1D reference interaction site model of the integral equation theory of molecular liquids. These functions can be obtained in a matter of minutes for most small organic and druglike molecules using existing software (RISM-MOL) (Sergiievskyi, V. P.; Hackbusch, W.; Fedorov, M. V. J. Comput. Chem. 2011, 32, 1982-1992). Predictions of caco-2 cell permeability and hydration free energy obtained using the RISM-MOL-INF method are shown to be more accurate than the state-of-the-art tools for benchmark data sets. Due to the importance of solvation and desolvation effects in biological systems, it is anticipated that the RISM-MOL-INF approach will find many applications in biophysical and biomedical property prediction.
Collapse
Affiliation(s)
- David S Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde , Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, U.K
| | - Maksim Mišin
- Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde , John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, U.K
| | - Maxim V Fedorov
- Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde , John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, U.K
| | - Antonio Llinas
- Respiratory, Inflammation and Autoimmune iMed, AstraZeneca R&D , Pepparedsleden 1, SE-431 83, Mölndal, Sweden
| |
Collapse
|
14
|
Ratkova EL, Palmer DS, Fedorov MV. Solvation thermodynamics of organic molecules by the molecular integral equation theory: approaching chemical accuracy. Chem Rev 2015; 115:6312-56. [PMID: 26073187 DOI: 10.1021/cr5000283] [Citation(s) in RCA: 135] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Ekaterina L Ratkova
- †G. A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, Akademicheskaya Street 1, Ivanovo 153045, Russia.,‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany
| | - David S Palmer
- ‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany.,§Department of Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, United Kingdom
| | - Maxim V Fedorov
- ‡The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig 04103, Germany.,∥Department of Physics, Scottish Universities Physics Alliance (SUPA), University of Strathclyde, John Anderson Building, 107 Rottenrow East, Glasgow G4 0NG, United Kingdom
| |
Collapse
|
15
|
Sumi T, Mitsutake A, Maruyama Y. A solvation-free-energy functional: A reference-modified density functional formulation. J Comput Chem 2015; 36:1359-69. [DOI: 10.1002/jcc.23942] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2015] [Revised: 02/26/2015] [Accepted: 03/29/2015] [Indexed: 01/21/2023]
Affiliation(s)
- Tomonari Sumi
- Department of Chemistry; Faculty of Science, Okayama University; 3-1-1 Tsushima-Naka Kita-ku Okayama 700-8530 Japan
| | - Ayori Mitsutake
- Department of Physics; Keio University; 3-14-1, Hiyoshi Kohoku-ku, Yokohama Kanagawa 223-8522 Japan
- JST, PREST; 3-14-1 Hiyoshi Kohoku-ku, Yokohama Kanagawa 223-8522 Japan
| | - Yutaka Maruyama
- Department of Physics; Keio University; 3-14-1, Hiyoshi Kohoku-ku, Yokohama Kanagawa 223-8522 Japan
| |
Collapse
|
16
|
|
17
|
Palmer DS, Mitchell JBO. Is Experimental Data Quality the Limiting Factor in Predicting the Aqueous Solubility of Druglike Molecules? Mol Pharm 2014; 11:2962-72. [DOI: 10.1021/mp500103r] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- David S. Palmer
- Department
of Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral
Street, Glasgow, Scotland G1 1XL, U.K
| | - John B. O. Mitchell
- Biomedical
Sciences Research Complex and EaStCHEM School of Chemistry, University of St. Andrews, Purdie Building, North Haugh, St. Andrews, Scotland KY16 9ST, U.K
| |
Collapse
|
18
|
Tanaka S, Nakano M. Classical density functional calculation of radial distribution functions of liquid water. Chem Phys 2014. [DOI: 10.1016/j.chemphys.2013.12.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
19
|
Narayana Moorthy NSH, Martins SA, Sousa SF, Ramos MJ, Fernandes PA. Classification study of solvation free energies of organic molecules using machine learning techniques. RSC Adv 2014. [DOI: 10.1039/c4ra07961b] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Classification models to predict the solvation free energies of organic molecules were developed using decision tree, random forest and support vector machine approaches and with MACCS fingerprints, MOE and PaDEL descriptors.
Collapse
Affiliation(s)
- N. S. Hari Narayana Moorthy
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
- 4169-007 Porto, Portugal
| | - Silvia A. Martins
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
- 4169-007 Porto, Portugal
| | - Sergio F. Sousa
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
- 4169-007 Porto, Portugal
| | - Maria J. Ramos
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
- 4169-007 Porto, Portugal
| | - Pedro A. Fernandes
- REQUIMTE
- Departamento de Química e Bioquímica
- Faculdade de Ciências
- Universidade do Porto
- 4169-007 Porto, Portugal
| |
Collapse
|
20
|
Palmer DS, Sørensen J, Schiøtt B, Fedorov MV. Solvent Binding Analysis and Computational Alanine Scanning of the Bovine Chymosin–Bovine κ-Casein Complex Using Molecular Integral Equation Theory. J Chem Theory Comput 2013; 9:5706-17. [DOI: 10.1021/ct400605x] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David S. Palmer
- Department
of Physics, University of Strathclyde, John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, United Kingdom
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| | - Jesper Sørensen
- Department
of Chemistry and Biochemistry, University of California, San Diego, Urey Hall, 9500 Gilman Drive, La Jolla, California 92093, United States
- The
Center for Insoluble Protein Structures (inSPIN) and the Interdisciplinary
Nanoscience Center (iNANO), Department of Chemistry, Aarhus University, Langelandsgade
140, DK-8000 Aarhus
C, Denmark
| | - Birgit Schiøtt
- The
Center for Insoluble Protein Structures (inSPIN) and the Interdisciplinary
Nanoscience Center (iNANO), Department of Chemistry, Aarhus University, Langelandsgade
140, DK-8000 Aarhus
C, Denmark
| | - Maxim V. Fedorov
- Department
of Physics, University of Strathclyde, John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, United Kingdom
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| |
Collapse
|
21
|
Liu Y, Zhao S, Wu J. A Site Density Functional Theory for Water: Application to Solvation of Amino Acid Side Chains. J Chem Theory Comput 2013; 9:1896-908. [DOI: 10.1021/ct3010936] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yu Liu
- Departments of Chemical and
Environmental Engineering and Mathematics, University of California,
Riverside, California 92521, United States
| | - Shuangliang Zhao
- State Key Laboratory of Chemical
Engineering, East China University of Science and Technology, Shanghai,
200238, P. R. China
| | - Jianzhong Wu
- Departments of Chemical and
Environmental Engineering and Mathematics, University of California,
Riverside, California 92521, United States
| |
Collapse
|
22
|
Kinoshita M. A new theoretical approach to biological self-assembly. Biophys Rev 2013; 5:283-293. [PMID: 28510109 DOI: 10.1007/s12551-013-0100-8] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2012] [Accepted: 01/02/2013] [Indexed: 11/25/2022] Open
Abstract
Upon biological self-assembly, the number of accessible translational configurations of water in the system increases considerably, leading to a large gain in water entropy. It is important to calculate the solvation entropy of a biomolecule with a prescribed structure by accounting for the change in water-water correlations caused by solute insertion. Modeling water as a dielectric continuum is not capable of capturing the physical essence of the water entropy effect. As a reliable tool, we propose a hybrid of the angle-dependent integral equation theory combined with a multipolar water model and a morphometric approach. Using our methods wherein the water entropy effect is treated as the key factor, we can elucidate a variety of processes such as protein folding, cold, pressure, and heat denaturating of a protein, molecular recognition, ordered association of proteins such as amyloid fibril formation, and functioning of ATP-driven proteins.
Collapse
Affiliation(s)
- Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto, 611-0011, Japan.
| |
Collapse
|
23
|
Sato H. A modern solvation theory: quantum chemistry and statistical chemistry. Phys Chem Chem Phys 2013; 15:7450-65. [DOI: 10.1039/c3cp50247c] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
24
|
Fedotova MV, Kruchinin SE. Hydration of methylamine and methylammonium ion: structural and thermodynamic properties from the data of the integral equation method in the RISM approximation. Russ Chem Bull 2012. [DOI: 10.1007/s11172-012-0034-5] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
25
|
Ratkova EL. A semiempirical model for estimating the hydration free energy of neutral nonpolar compounds. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2012. [DOI: 10.1134/s0036024412100196] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
26
|
Palmer DS, McDonagh JL, Mitchell JBO, van Mourik T, Fedorov MV. First-Principles Calculation of the Intrinsic Aqueous Solubility of Crystalline Druglike Molecules. J Chem Theory Comput 2012; 8:3322-37. [PMID: 26605739 DOI: 10.1021/ct300345m] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We demonstrate that the intrinsic aqueous solubility of crystalline druglike molecules can be estimated with reasonable accuracy from sublimation free energies calculated using crystal lattice simulations and hydration free energies calculated using the 3D Reference Interaction Site Model (3D-RISM) of the Integral Equation Theory of Molecular Liquids (IET). The solubilities of 25 crystalline druglike molecules taken from different chemical classes are predicted by the model with a correlation coefficient of R = 0.85 and a root mean square error (RMSE) equal to 1.45 log10S units, which is significantly more accurate than results obtained using implicit continuum solvent models. The method is not directly parametrized against experimental solubility data, and it offers a full computational characterization of the thermodynamics of transfer of the drug molecule from crystal phase to gas phase to dilute aqueous solution.
Collapse
Affiliation(s)
- David S Palmer
- Department of Physics, University of Strathclyde , John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, United Kingdom.,Max Planck Institute for Mathematics in the Sciences , Inselstrasse 22, DE-04103 Leipzig, Germany
| | - James L McDonagh
- Biomedical Sciences Research Complex and EaStCHEM School of Chemistry, University of St. Andrews , Purdie Building, North Haugh, St. Andrews, Scotland KY16 9ST, United Kingdom
| | - John B O Mitchell
- Biomedical Sciences Research Complex and EaStCHEM School of Chemistry, University of St. Andrews , Purdie Building, North Haugh, St. Andrews, Scotland KY16 9ST, United Kingdom
| | - Tanja van Mourik
- Biomedical Sciences Research Complex and EaStCHEM School of Chemistry, University of St. Andrews , Purdie Building, North Haugh, St. Andrews, Scotland KY16 9ST, United Kingdom
| | - Maxim V Fedorov
- Department of Physics, University of Strathclyde , John Anderson Building, 107 Rottenrow, Glasgow, Scotland G4 0NG, United Kingdom.,Max Planck Institute for Mathematics in the Sciences , Inselstrasse 22, DE-04103 Leipzig, Germany
| |
Collapse
|
27
|
Sergiievskyi VP, Fedorov MV. 3DRISM Multigrid Algorithm for Fast Solvation Free Energy Calculations. J Chem Theory Comput 2012; 8:2062-70. [DOI: 10.1021/ct200815v] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Maxim V. Fedorov
- Nanoscience
Division, Department
of Physics, Scottish Universities Physics Alliance (SUPA), Strathclyde University, Room JA 6.10, John Anderson
Building 107, Rottenrow East Glasgow, United Kingdom G4 0NG
| |
Collapse
|
28
|
Iida K, Sato H. Proton Transfer Step in the Carbon Dioxide Capture by Monoethanol Amine: A Theoretical Study at the Molecular Level. J Phys Chem B 2012; 116:2244-8. [DOI: 10.1021/jp212225r] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Kenji Iida
- Department of Molecular Engineering, Kyoto University, Kyoto, 615-8510, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Kyoto University, Kyoto, 615-8510, Japan
| |
Collapse
|
29
|
Fedotova MV, Kruchinin SE. Hydration of acetic acid and acetate ion in water studied by 1D-RISM theory. J Mol Liq 2011. [DOI: 10.1016/j.molliq.2011.09.015] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
30
|
Chen Z, Baker NA, Wei GW. Differential geometry based solvation model II: Lagrangian formulation. J Math Biol 2011; 63:1139-200. [PMID: 21279359 PMCID: PMC3113640 DOI: 10.1007/s00285-011-0402-z] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2010] [Revised: 12/24/2010] [Indexed: 10/18/2022]
Abstract
Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation models. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The optimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and PB equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of computation, thanks to the equivalence of the Laplace-Beltrami operator in the two representations. The coupled partial differential equations (PDEs) are solved with an iterative procedure to reach a steady state, which delivers desired solvent-solute interface and electrostatic potential for problems of interest. These quantities are utilized to evaluate the solvation free energies and protein-protein binding affinities. A number of computational methods and algorithms are described for the interconversion of Lagrangian and Eulerian representations, and for the solution of the coupled PDE system. The proposed approaches have been extensively validated. We also verify that the mean curvature flow indeed gives rise to the minimal molecular surface and the proposed variational procedure indeed offers minimal total free energy. Solvation analysis and applications are considered for a set of 17 small compounds and a set of 23 proteins. The salt effect on protein-protein binding affinity is investigated with two protein complexes by using the present model. Numerical results are compared to the experimental measurements and to those obtained by using other theoretical methods in the literature.
Collapse
Affiliation(s)
- Zhan Chen
- Department of Mathematics, Michigan State University, MI 48824, USA
| | - Nathan A. Baker
- Pacific Northwest National Laboratory,
902 Battelle Boulevard P.O. Box 999, MSIN K7-28, Richland, WA 99352 USA
| | - G. W. Wei
- Department of Mathematics, Michigan State University, MI 48824, USA
- Department of Electrical and Computer Engineering, Michigan State University, MI 48824, USA
| |
Collapse
|
31
|
Schneider G. From Hits to Leads: Challenges for the Next Phase of Machine Learning in Medicinal Chemistry. Mol Inform 2011; 30:759-63. [DOI: 10.1002/minf.201100070] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2011] [Accepted: 05/23/2011] [Indexed: 12/12/2022]
|
32
|
Palmer DS, Frolov AI, Ratkova EL, Fedorov MV. Toward a Universal Model To Calculate the Solvation Thermodynamics of Druglike Molecules: The Importance of New Experimental Databases. Mol Pharm 2011; 8:1423-9. [DOI: 10.1021/mp200119r] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David S. Palmer
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| | - Andrey I. Frolov
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| | - Ekaterina L. Ratkova
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| | - Maxim V. Fedorov
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| |
Collapse
|
33
|
Sergiievskyi VP. Model for calculating the free energy of hydration of bioactive compounds based on integral equation theory of liquids. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2011. [DOI: 10.1134/s1990793111020382] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
34
|
Ratkova EL, Fedorov MV. Combination of RISM and Cheminformatics for Efficient Predictions of Hydration Free Energy of Polyfragment Molecules: Application to a Set of Organic Pollutants. J Chem Theory Comput 2011; 7:1450-7. [DOI: 10.1021/ct100654h] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ekaterina L. Ratkova
- The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
| | - Maxim V. Fedorov
- The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
| |
Collapse
|
35
|
Frolov AI, Ratkova EL, Palmer DS, Fedorov MV. Hydration Thermodynamics Using the Reference Interaction Site Model: Speed or Accuracy? J Phys Chem B 2011; 115:6011-22. [DOI: 10.1021/jp111271c] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Andrey I. Frolov
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
| | - Ekaterina L. Ratkova
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
| | - David S. Palmer
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
| | - Maxim V. Fedorov
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
| |
Collapse
|
36
|
Sergiievskyi VP, Hackbusch W, Fedorov MV. Multigrid solver for the reference interaction site model of molecular liquids theory. J Comput Chem 2011; 32:1982-92. [DOI: 10.1002/jcc.21783] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2010] [Revised: 01/22/2011] [Accepted: 02/09/2011] [Indexed: 11/12/2022]
|
37
|
|
38
|
Yamazaki T, Kovalenko A. Spatial Decomposition of Solvation Free Energy Based on the 3D Integral Equation Theory of Molecular Liquid: Application to Miniproteins. J Phys Chem B 2010; 115:310-8. [DOI: 10.1021/jp1082938] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Takeshi Yamazaki
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, and Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada
| | - Andriy Kovalenko
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, and Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada
| |
Collapse
|
39
|
Palmer DS, Frolov AI, Ratkova EL, Fedorov MV. Towards a universal method for calculating hydration free energies: a 3D reference interaction site model with partial molar volume correction. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2010; 22:492101. [PMID: 21406779 DOI: 10.1088/0953-8984/22/49/492101] [Citation(s) in RCA: 89] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We report a simple universal method to systematically improve the accuracy of hydration free energies calculated using an integral equation theory of molecular liquids, the 3D reference interaction site model. A strong linear correlation is observed between the difference of the experimental and (uncorrected) calculated hydration free energies and the calculated partial molar volume for a data set of 185 neutral organic molecules from different chemical classes. By using the partial molar volume as a linear empirical correction to the calculated hydration free energy, we obtain predictions of hydration free energies in excellent agreement with experiment (R = 0.94, σ = 0.99 kcal mol (- 1) for a test set of 120 organic molecules).
Collapse
Affiliation(s)
- David S Palmer
- Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, DE-04103 Leipzig, Germany
| | | | | | | |
Collapse
|
40
|
Ratkova EL, Chuev GN, Sergiievskyi VP, Fedorov MV. An Accurate Prediction of Hydration Free Energies by Combination of Molecular Integral Equations Theory with Structural Descriptors. J Phys Chem B 2010; 114:12068-79. [DOI: 10.1021/jp103955r] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Ekaterina L. Ratkova
- The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany, and Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region, 142290, Russia
| | - Gennady N. Chuev
- The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany, and Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region, 142290, Russia
| | - Volodymyr P. Sergiievskyi
- The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany, and Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region, 142290, Russia
| | - Maxim V. Fedorov
- The Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany, and Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region, 142290, Russia
| |
Collapse
|