1
|
Low K, Coote ML, Izgorodina EI. Explainable Solvation Free Energy Prediction Combining Graph Neural Networks with Chemical Intuition. J Chem Inf Model 2022; 62:5457-5470. [PMID: 36317829 DOI: 10.1021/acs.jcim.2c01013] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The prediction of a molecule's solvation Gibbs free (ΔGsolv) energy in a given solvent is an important task which has traditionally been carried out via quantum chemical continuum methods or force field-based molecular simulations. Machine learning (ML) and graph neural networks in particular have emerged as powerful techniques for elucidating structure-property relationships. This work presents a graph neural network (GNN) for the prediction of ΔGsolv which, in addition to encoding typical atom and bond-level features, incorporates chemically intuitive, solvation-relevant parameters into the featurization process: semiempirical partial atomic charges and solvent dielectric constant. Solute-solvent interactions are included via an interaction map layer which can be visualized to examine solubility-enhancing or -decreasing interactions learnt by the model. On a test set of small organic molecules, our GNN predicts ΔGsolv in water and cyclohexane with an accuracy comparable to polarizable and ab initio generated force field methods [mean absolute error (MAE) = 0.4 and 0.2 kcal mol-1, respectively], without the need for any molecular simulation. For the FreeSolv data set of hydration free energies, the test MAE is 0.7 kcal mol-1. Interpretability and applicability of the model is highlighted through several examples including rationalizing the increased solubility of modified diaminoanthraquinones in organic solvents. The clear explanations afforded by our GNN allow for easy understanding of the model's predictions, giving the experimental chemist confidence in employing ML models toward more optimized synthetic routes.
Collapse
Affiliation(s)
- Kaycee Low
- Monash Computational Chemistry Group, School of Chemistry, Monash University, Clayton, Victoria3800, Australia
| | - Michelle L Coote
- Institute for Nanoscale Science and Technology, College of Science and Engineering, Flinders University, Bedford Park, South Australia5042, Australia
| | - Ekaterina I Izgorodina
- Monash Computational Chemistry Group, School of Chemistry, Monash University, Clayton, Victoria3800, Australia
| |
Collapse
|
2
|
Ganyecz Á, Kállay M. Implementation and Optimization of the Embedded Cluster Reference Interaction Site Model with Atomic Charges. J Phys Chem A 2022; 126:2417-2429. [PMID: 35394778 PMCID: PMC9036516 DOI: 10.1021/acs.jpca.1c07904] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
In this work, we
implemented the embedded cluster reference interaction
site model (EC-RISM) originally developed by Kloss, Heil, and Kast
(J. Phys. Chem. B2008, 112, 4337–4343).
This method combines quantum mechanical calculations with the 3D reference
interaction site model (3D-RISM). Numerous options, such as buffer,
grid space, basis set, charge model, water model, closure relation,
and so forth, were investigated to find the best settings. Additionally,
the small point charges, which are derived from the solvent distribution
from the 3D-RISM solution to represent the solvent in the QM calculation,
were neglected to reduce the overhead without the loss of accuracy.
On the MNSOL[a], MNSOL, and FreeSolv databases, our implemented and
optimized method provides solvation free energies in water with 5.70,
6.32, and 6.44 kJ/mol root-mean-square deviations, respectively, but
with different settings, 5.22, 6.08, and 6.63 kJ/mol can also be achieved.
Only solvent models containing fitting parameters, like COSMO-RS and
EC-RISM with universal correction and directly used electrostatic
potential, perform better than our EC-RISM implementation with atomic
charges.
Collapse
Affiliation(s)
- Ádám Ganyecz
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest P.O. Box 91, H-1521 Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest P.O. Box 91, H-1521 Hungary
| |
Collapse
|
3
|
Chuev GN, Fedotova MV, Valiev M. Renormalized site density functional theory for models of ion hydration. J Chem Phys 2021; 155:064501. [PMID: 34391371 DOI: 10.1063/5.0060249] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile the co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. The renormalized formulation of SDFT (RSDFT), which we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system-a problem that historically has always been one of the most difficult cases for SDFT applications. Using a two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional three-dimensional reference interaction site model implementations and may prove useful in coarse grained simulations based on two-site solvent models.
Collapse
Affiliation(s)
- Gennady N Chuev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region 142290, Russia
| | - Marina V Fedotova
- G. A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Marat Valiev
- Molecular Sciences Software Group, Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
4
|
Hsu TY, Jeanmairet G. Assessing the correctness of pressure correction to solvation theories in the study of electron transfer reactions. J Chem Phys 2021; 154:131102. [PMID: 33832266 DOI: 10.1063/5.0048343] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Liquid state theories have emerged as a numerically efficient alternative to costly molecular dynamics simulations of electron transfer reactions in solution. In a recent paper [Jeanmairet et al., Chem. Sci. 10, 2130-2143 (2019)], we introduced the framework to compute the energy gap, free energy profile, and reorganization free energy using molecular density functional theory. However, this technique, as other molecular liquid state theories, overestimates the bulk pressure of the fluid. Because of the very high pressure, the predicted free energy is dramatically exaggerated. Several attempts were made to fix this issue, either based on simple a posteriori correction or by introducing bridge terms. By studying two model half reactions in water, Cl → Cl+ and Cl → Cl-, we assess the correctness of these two types of corrections to study electron transfer reactions. We found that a posteriori correction, because it violates the Variational principle, leads to an inconsistency in the definition of the reorganization free energy and should not be used to study electron transfer reactions. The bridge approach, because it is theoretically well grounded, is perfectly suitable for this type of systems.
Collapse
Affiliation(s)
- Tzu-Yao Hsu
- Sorbonne Université, CNRS, Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| | - Guillaume Jeanmairet
- Sorbonne Université, CNRS, Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| |
Collapse
|
5
|
Zheng D, Yuan Y, Wang F. Determining the hydration free energies of selected small molecules with MP2 and local MP2 through adaptive force matching. J Chem Phys 2021; 154:104113. [PMID: 33722038 DOI: 10.1063/5.0044712] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Force fields for seven small solute molecules, ethanol, 2-methyl-1-propanol, 2-butanol, cyclohexene, tetrahydropyran, 1,4-dioxane, and 1,4-butanediol, in dilute aqueous solutions were created with the adaptive force matching (AFM) method using MP2 or local MP2 as reference. The force fields provide a way to predict the hydration free energies (HFEs) of these molecules with only electronic structure calculations as reference. For six of the seven molecules, the predicted HFEs are in very good agreement with experiments. For 1,4-butanediol, the model created by force matching LMP2 provides a HFE that is too positive. Further investigation suggests that LMP2 may not be sufficiently accurate for computing HFEs for alcohols with AFM. Other properties, such as enthalpy of hydration, diffusion constants, and vibrational spectra, are also computed with the force field developed. The force fields developed by AFM provide a bridge for computing ensemble properties of the reference electronic structure method. With MP2 and LMP2 as reference methods, the computed properties of the small molecular solutes are found to be in good agreement with experiments.
Collapse
Affiliation(s)
- Dong Zheng
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Ying Yuan
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
6
|
Abstract
The electrostatic response underlying the 3D-RISM theory and its general relationship to models in which the solvent is represented in terms of a dielectric continuum are examined. It is found that the theory provides a coherent picture of solvation, although its behavior is not entirely consistent with the trends that are expected in the limit of a large solute. The electrostatic discrepancy is due to the nature of the isotropic pair additive site-site correlation function associated with the susceptibility response of the uniform fluid. The influence of the discrepancy in the magnitude of the solvation free energy is negligible for a solvent with a large dielectric constant.
Collapse
Affiliation(s)
- Benoît Roux
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637-1454, United States
| |
Collapse
|
7
|
Robert A, Luukkonen S, Levesque M. Pressure correction for solvation theories. J Chem Phys 2020; 152:191103. [DOI: 10.1063/5.0002029] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Anton Robert
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Sohvi Luukkonen
- Maison de la Simulation, CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Maximilien Levesque
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
- Aqemia, Paris, France
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
Roy D, Kovalenko A. Performance of 3D-RISM-KH in Predicting Hydration Free Energy: Effect of Solute Parameters. J Phys Chem A 2019; 123:4087-4093. [PMID: 30993994 DOI: 10.1021/acs.jpca.9b01623] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The three-dimensional reference interaction site model molecular solvation theory with the Kovalenko-Hirata closure relation has been shown to produce excellent solvation characteristics for a large class of (bio)chemical systems in solution. Correct calculation of hydration free energy is central to successful application of any solvation model. In order to find out the best possible force-field parameters to be used for hydration free energy calculation with the aforementioned theory, we have developed an extended database containing a large number of experimental solvation free energies available in the current literature and used a plethora of theoretical models for assessment. The general Amber force field was found to perform satisfactorily, whereas special care should be taken in solute charge assignment with the universal force field.
Collapse
Affiliation(s)
- Dipankar Roy
- Department of Mechanical Engineering , University of Alberta , 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Street NW , Edmonton , Alberta T6G 1H9 , Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering , University of Alberta , 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Street NW , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , 11421 Saskatchewan Drive , Edmonton , Alberta T6G 2M9 , Canada
| |
Collapse
|
11
|
Kido K. A noniterative mean‐field QM/MM‐type approach with a linear response approximation toward an efficient free‐energy evaluation. J Comput Chem 2019; 40:2072-2085. [DOI: 10.1002/jcc.25844] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Revised: 03/22/2019] [Accepted: 04/01/2019] [Indexed: 01/23/2023]
Affiliation(s)
- Kentaro Kido
- Nuclear Safety Research CenterJapan Atomic Energy Agency 2‐4 Shirane, Shirakata, Tokai‐mura Ibaraki 319‐1195 Naka‐gun Japan
| |
Collapse
|
12
|
Kılıç H. Experiment-based physicochemical aspects for the coulombic hydration kinetics and thermodynamics of a pyrimidine and thiopyrimidine. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.06.090] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
13
|
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
|
14
|
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
|
15
|
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
|
16
|
Riniker S. Molecular Dynamics Fingerprints (MDFP): Machine Learning from MD Data To Predict Free-Energy Differences. J Chem Inf Model 2017; 57:726-741. [DOI: 10.1021/acs.jcim.6b00778] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
17
|
Wang B, Zhao Z, Wei GW. Automatic parametrization of non-polar implicit solvent models for the blind prediction of solvation free energies. J Chem Phys 2016; 145:124110. [DOI: 10.1063/1.4963193] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Affiliation(s)
- Bao Wang
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA
| | - Zhixiong Zhao
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA
- School of Medicine, Foshan University, Foshan, Guangdong 528000, China
| | - Guo-Wei Wei
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, USA
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
18
|
Sumi T, Maruyama Y, Mitsutake A, Koga K. A reference-modified density functional theory: An application to solvation free-energy calculations for a Lennard-Jones solution. J Chem Phys 2016; 144:224104. [DOI: 10.1063/1.4953191] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
|
19
|
Zafar A, Reynisson J. Hydration Free Energy as a Molecular Descriptor in Drug Design: A Feasibility Study. Mol Inform 2016; 35:207-14. [PMID: 27492087 DOI: 10.1002/minf.201501035] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2015] [Accepted: 03/11/2016] [Indexed: 01/29/2023]
Abstract
In this work the idea was investigated whether calculated hydration energy (ΔGhyd ) can be used as a molecular descriptor in defining promising regions of chemical space for drug design. Calculating ΔGhyd using the Density Solvation Model (SMD) in conjunction with the density functional theory (DFT) gave an excellent correlation with experimental values. Furthermore, calculated ΔGhyd correlates reasonably well with experimental water solubility (r(2) =0.545) and also log P (r(2) =0.530). Three compound collections were used: Known drugs (n=150), drug-like compounds (n=100) and simple organic compounds (n=140). As an approximation only molecules, which do not de/protonate at physiological pH were considered. A relatively broad distribution was seen for the known drugs with an average at -15.3 kcal/mol and a standard deviation of 7.5 kcal/mol. Interestingly, much lower averages were found for the drug-like compounds (-7.5 kcal/mol) and the simple organic compounds (-3.1 kcal/mol) with tighter distributions; 4.3 and 3.2 kcal/mol, respectively. This trend was not observed for these collections when calculated log P and log S values were used. The considerable greater exothermic ΔGhyd average for the known drugs clearly indicates in order to develop a successful drug candidate value of ΔGhyd <-5 kcal/mol or less is preferable.
Collapse
Affiliation(s)
- Ayesha Zafar
- School of Chemical Sciences, University of Auckland, Private Bag 92019, Auckland 1142, New Zealand, Tel. +64-9-373-7599 ext. 83746, Fax. +64-9-373-7422
| | - Jóhannes Reynisson
- School of Chemical Sciences, University of Auckland, Private Bag 92019, Auckland 1142, New Zealand, Tel. +64-9-373-7599 ext. 83746, Fax. +64-9-373-7422.
| |
Collapse
|
20
|
Misin M, Fedorov MV, Palmer DS. Hydration Free Energies of Molecular Ions from Theory and Simulation. J Phys Chem B 2016; 120:975-83. [PMID: 26756333 DOI: 10.1021/acs.jpcb.5b10809] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a theoretical/computational framework for accurate calculation of hydration free energies of ionized molecular species. The method is based on a molecular theory, 3D-RISM, combined with a recently developed pressure correction (PC+). The 3D-RISM/PC+ model can provide ∼3 kcal/mol hydration free energy accuracy for a large variety of ionic compounds, provided that the Galvani potential of water is taken into account. The results are compared with direct atomistic simulations. Several methodological aspects of hydration free energy calculations for charged species are discussed.
Collapse
Affiliation(s)
| | | | - David S Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde , 295 Cathedral Street, Glasgow, G1 1XL, United Kingdom
| |
Collapse
|
21
|
Wang B, Wei GW. Parameter optimization in differential geometry based solvation models. J Chem Phys 2015; 143:134119. [PMID: 26450304 PMCID: PMC4602332 DOI: 10.1063/1.4932342] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2015] [Accepted: 09/22/2015] [Indexed: 01/01/2023] Open
Abstract
Differential geometry (DG) based solvation models are a new class of variational implicit solvent approaches that are able to avoid unphysical solvent-solute boundary definitions and associated geometric singularities, and dynamically couple polar and non-polar interactions in a self-consistent framework. Our earlier study indicates that DG based non-polar solvation model outperforms other methods in non-polar solvation energy predictions. However, the DG based full solvation model has not shown its superiority in solvation analysis, due to its difficulty in parametrization, which must ensure the stability of the solution of strongly coupled nonlinear Laplace-Beltrami and Poisson-Boltzmann equations. In this work, we introduce new parameter learning algorithms based on perturbation and convex optimization theories to stabilize the numerical solution and thus achieve an optimal parametrization of the DG based solvation models. An interesting feature of the present DG based solvation model is that it provides accurate solvation free energy predictions for both polar and non-polar molecules in a unified formulation. Extensive numerical experiment demonstrates that the present DG based solvation model delivers some of the most accurate predictions of the solvation free energies for a large number of molecules.
Collapse
Affiliation(s)
- Bao Wang
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA
| | - G W Wei
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
22
|
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
|
23
|
Phanich J, Rungrotmongkol T, Sindhikara D, Phongphanphanee S, Yoshida N, Hirata F, Kungwan N, Hannongbua S. A 3D-RISM/RISM study of the oseltamivir binding efficiency with the wild-type and resistance-associated mutant forms of the viral influenza B neuraminidase. Protein Sci 2015; 25:147-58. [PMID: 26044768 DOI: 10.1002/pro.2718] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2015] [Revised: 05/27/2015] [Accepted: 05/28/2015] [Indexed: 11/09/2022]
Abstract
The binding affinity of oseltamivir to the influenza B neuraminidase and to its variants with three single substitutions, E119G, R152K, and D198N, is investigated by the MM/3D-RISM method. The binding affinity or the binding free energy of ligand to receptor was found to be determined by a subtle balance of two major contributions that largely cancel out each other: the ligand-receptor interactions and the dehydration free energy. The theoretical results of the binding affinity of the drug to the mutants reproduced the observed trend in the resistivity, measured by IC50 ; the high-level resistance of E119G and R152K, and the low-level resistance of D198N. For E119G and R152K, reduction of the direct drug-target interaction, especially at the mutated residue, is the main source of high-level oseltamivir resistance. This phenomenon, however, is not found in the D198N strain, which is located in the framework of the active-site.
Collapse
Affiliation(s)
- Jiraphorn Phanich
- Department of Chemistry, Computational Chemistry Unit Cell, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Thanyada Rungrotmongkol
- Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.,Ph.D. Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Daniel Sindhikara
- Schrödinger, Inc, 120 West 45th Street, 17th Floor, New York, New York, 10036
| | - Saree Phongphanphanee
- Department of Materials Science, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Sciences, Kyushu University, Fukuoka, 812-8581, Japan
| | - Fumio Hirata
- College of Life Sciences, Ritsumeikan University, and Molecular Design Frontier Co. Ltd, Kusatsu, 525-8577, Japan
| | - Nawee Kungwan
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, 50200, Thailand
| | - Supot Hannongbua
- Department of Chemistry, Computational Chemistry Unit Cell, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| |
Collapse
|
24
|
Raevsky OA, Polianczyk DE, Grigorev VY, Raevskaja OE, Dearden JC. In silico Prediction of Aqueous Solubility: a Comparative Study of Local and Global Predictive Models. Mol Inform 2015; 34:417-30. [PMID: 27490387 DOI: 10.1002/minf.201400144] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2014] [Accepted: 03/05/2015] [Indexed: 11/07/2022]
Abstract
32 Quantitative Structure-Property Relationship (QSPR) models were constructed for prediction of aqueous intrinsic solubility of liquid and crystalline chemicals. Data sets contained 1022 liquid and 2615 crystalline compounds. Multiple Linear Regression (MLR), Support Vector Machine (SVM) and Random Forest (RF) methods were used to construct global models, and k-nearest neighbour (kNN), Arithmetic Mean Property (AMP) and Local Regression Property (LoReP) were used to construct local models. A set of the best QSPR models was obtained: for liquid chemicals with RMSE (root mean square error) of prediction in the range 0.50-0.60 log unit; for crystalline chemicals 0.80-0.90 log unit. In the case of global models the large number of descriptors makes mechanistic interpretation difficult. The local models use only one or two descriptors, so that a medicinal chemist working with sets of structurally-related chemicals can readily estimate their solubility. However, construction of stable local models requires the presence of closely related neighbours for each chemical considered. It is probable that a consensus of global and local QSPR models will be the optimal approach for construction of stable predictive QSPR models with mechanistic interpretation.
Collapse
Affiliation(s)
- Oleg A Raevsky
- Department of Computer-Aided Molecular Design, Institute of Physiologically Active Compounds, Russian Academy of Science, 142432, Russia, Chernogolovka, Severniy proezd 1 phone: +7 496 52 21867.
| | - Daniel E Polianczyk
- Department of Computer-Aided Molecular Design, Institute of Physiologically Active Compounds, Russian Academy of Science, 142432, Russia, Chernogolovka, Severniy proezd 1 phone: +7 496 52 21867
| | - Veniamin Yu Grigorev
- Department of Computer-Aided Molecular Design, Institute of Physiologically Active Compounds, Russian Academy of Science, 142432, Russia, Chernogolovka, Severniy proezd 1 phone: +7 496 52 21867
| | - Olga E Raevskaja
- Department of Computer-Aided Molecular Design, Institute of Physiologically Active Compounds, Russian Academy of Science, 142432, Russia, Chernogolovka, Severniy proezd 1 phone: +7 496 52 21867
| | - John C Dearden
- School of Pharmacy and Biomolecular Sciences, Liverpool John Moores University, Liverpool L3 3AF, UK
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
Frolov AI. Accurate Calculation of Solvation Free Energies in Supercritical Fluids by Fully Atomistic Simulations: Probing the Theory of Solutions in Energy Representation. J Chem Theory Comput 2015; 11:2245-56. [DOI: 10.1021/acs.jctc.5b00172] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Andrey I. Frolov
- Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya St. 1, 153045 Ivanovo, Russia
| |
Collapse
|
28
|
Zhang R, Chen J, Mi J. Evaluation of solvent effect on the amine-based CO2 absorbents by theory and experiment. Theor Chem Acc 2015. [DOI: 10.1007/s00214-015-1637-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
29
|
Chong SH, Ham S. Thermodynamic-Ensemble Independence of Solvation Free Energy. J Chem Theory Comput 2015; 11:378-80. [DOI: 10.1021/ct500876x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Song-Ho Chong
- Department of Chemistry, Sookmyung Women’s University, Cheongpa-ro 47-gil 100, Yongsan-Ku, Seoul 140-742, Korea
| | - Sihyun Ham
- Department of Chemistry, Sookmyung Women’s University, Cheongpa-ro 47-gil 100, Yongsan-Ku, Seoul 140-742, Korea
| |
Collapse
|
30
|
Wang X, Wang S, Xu Q, Mi J. Thermodynamics of Ice Nucleation in Liquid Water. J Phys Chem B 2015; 119:1660-8. [DOI: 10.1021/jp512280p] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xin Wang
- State Key Laboratory of Organic-Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| | - Shui Wang
- State Key Laboratory of Organic-Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| | - Qinzhi Xu
- Institute of Microelectronics of Chinese Academy of Sciences, Beijing, 100029, China
| | - Jianguo Mi
- State Key Laboratory of Organic-Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
31
|
Li B, Matveev AV, Krüger S, Rösch N. Uranyl solvation by a reference interaction site model. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2014.10.032] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
32
|
Zhao S, Liu Y, Chen X, Lu Y, Liu H, Hu Y. Unified Framework of Multiscale Density Functional Theories and Its Recent Applications. MESOSCALE MODELING IN CHEMICAL ENGINEERING PART II 2015. [DOI: 10.1016/bs.ache.2015.10.001] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
33
|
Frolov AI, Kiselev MG. Prediction of Cosolvent Effect on Solvation Free Energies and Solubilities of Organic Compounds in Supercritical Carbon Dioxide Based on Fully Atomistic Molecular Simulations. J Phys Chem B 2014; 118:11769-80. [DOI: 10.1021/jp505731z] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Andrey I. Frolov
- Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya
St. 1, 153045 Ivanovo, Russia
| | - Michael G. Kiselev
- Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya
St. 1, 153045 Ivanovo, Russia
| |
Collapse
|
34
|
Sergiievskyi VP, Jeanmairet G, Levesque M, Borgis D. Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections. J Phys Chem Lett 2014; 5:1935-1942. [PMID: 26273876 DOI: 10.1021/jz500428s] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Molecular density functional theory (MDFT) offers an efficient implicit-solvent method to estimate molecule solvation free-energies, whereas conserving a fully molecular representation of the solvent. Even within a second-order approximation for the free-energy functional, the so-called homogeneous reference fluid approximation, we show that the hydration free-energies computed for a data set of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by 2-3 orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canonical to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justification to empirical partial molar volume corrections that have been proposed recently.
Collapse
Affiliation(s)
- Volodymyr P Sergiievskyi
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Guillaume Jeanmairet
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Maximilien Levesque
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Daniel Borgis
- †École Normale Supérieure - PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
- ‡Maison de la Simulation, CNRS USR 3441, CEA Saclay, 91191 Gif-sur-Yvette, France
| |
Collapse
|
35
|
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
|
36
|
Anashkin IP, Klinov AV. Determining the parameters of the potential of intermolecular interaction by the Zeno line. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2013. [DOI: 10.1134/s0036024413110034] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
37
|
Liu Y, Wu J. Communication: Long-range angular correlations in liquid water. J Chem Phys 2013; 139:041103. [DOI: 10.1063/1.4817321] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
38
|
Joung IS, Luchko T, Case DA. Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions. J Chem Phys 2013; 138:044103. [PMID: 23387564 DOI: 10.1063/1.4775743] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Using the dielectrically consistent reference interaction site model (DRISM) of molecular solvation, we have calculated structural and thermodynamic information of alkali-halide salts in aqueous solution, as a function of salt concentration. The impact of varying the closure relation used with DRISM is investigated using the partial series expansion of order-n (PSE-n) family of closures, which includes the commonly used hypernetted-chain equation (HNC) and Kovalenko-Hirata closures. Results are compared to explicit molecular dynamics (MD) simulations, using the same force fields, and to experiment. The mean activity coefficients of ions predicted by DRISM agree well with experimental values at concentrations below 0.5 m, especially when using the HNC closure. As individual ion activities (and the corresponding solvation free energies) are not known from experiment, only DRISM and MD results are directly compared and found to have reasonably good agreement. The activity of water directly estimated from DRISM is nearly consistent with values derived from the DRISM ion activities and the Gibbs-Duhem equation, but the changes in the computed pressure as a function of salt concentration dominate these comparisons. Good agreement with experiment is obtained if these pressure changes are ignored. Radial distribution functions of NaCl solution at three concentrations were compared between DRISM and MD simulations. DRISM shows comparable water distribution around the cation, but water structures around the anion deviate from the MD results; this may also be related to the high pressure of the system. Despite some problems, DRISM-PSE-n is an effective tool for investigating thermodynamic properties of simple electrolytes.
Collapse
Affiliation(s)
- In Suk Joung
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey 08854, USA
| | | | | |
Collapse
|
39
|
Stukan MR, Asmadi A, Abdallah W. Bulk properties of SWM4-NDP water model at elevated temperature and pressure. J Mol Liq 2013. [DOI: 10.1016/j.molliq.2012.12.023] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
40
|
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
|
41
|
Fedotova MV, Kruchinin SE. The hydration of aniline and benzoic acid: Analysis of radial and spatial distribution functions. J Mol Liq 2013. [DOI: 10.1016/j.molliq.2012.11.031] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
42
|
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]
|
43
|
Karino Y, Matubayasi N. Interaction-component analysis of the urea effect on amino acid analogs. Phys Chem Chem Phys 2013; 15:4377-91. [DOI: 10.1039/c3cp43346c] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
44
|
Wang X, Chen J, Mi J. Structure-Solubility Correlation Model for Carbon Dioxide in Ionic Liquids. Ind Eng Chem Res 2012. [DOI: 10.1021/ie3028747] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Xin Wang
- State Key Laboratory
of Organic−Inorganic Composites, Beijing University of Chemical Technology, Beijing
100029, People’s Republic of China
| | - Jian Chen
- State Key Laboratory
of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084,
People’s Republic of China
| | - Jianguo Mi
- State Key Laboratory
of Organic−Inorganic Composites, Beijing University of Chemical Technology, Beijing
100029, People’s Republic of China
| |
Collapse
|
45
|
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]
|
46
|
Chen Z, Zhao S, Chun J, Thomas DG, Baker NA, Bates PW, Wei GW. Variational approach for nonpolar solvation analysis. J Chem Phys 2012; 137:084101. [PMID: 22938212 PMCID: PMC3436914 DOI: 10.1063/1.4745084] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2012] [Accepted: 07/27/2012] [Indexed: 12/23/2022] Open
Abstract
Solvation analysis is one of the most important tasks in chemical and biological modeling. Implicit solvent models are some of the most popular approaches. However, commonly used implicit solvent models rely on unphysical definitions of solvent-solute boundaries. Based on differential geometry, the present work defines the solvent-solute boundary via the variation of the nonpolar solvation free energy. The solvation free energy functional of the system is constructed based on a continuum description of the solvent and the discrete description of the solute, which are dynamically coupled by the solvent-solute boundaries via van der Waals interactions. The first variation of the energy functional gives rise to the governing Laplace-Beltrami equation. The present model predictions of the nonpolar solvation energies are in an excellent agreement with experimental data, which supports the validity of the proposed nonpolar solvation model.
Collapse
Affiliation(s)
- Zhan Chen
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, USA
| | | | | | | | | | | | | |
Collapse
|
47
|
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
|
48
|
Sergiievskyi VP, Frolov AI. A universal bridge functional for infinitely diluted solutions: A case study for Lennard-Jones spheres of different diameters. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2012. [DOI: 10.1134/s0036024412080122] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
49
|
Xu M, Zhang C, Du Z, Mi J. Structures and Surface Tensions of Fluids near Solid Surfaces: An Integral Equation Theory Study. J Phys Chem B 2012; 116:6514-21. [DOI: 10.1021/jp301991k] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Mengjin Xu
- The Key Laboratory of Carbon
Fiber and Functional Polymers, Ministry of Education, Beijing University of Chemical Technology, Beijing
100029, Peopleʼs Republic of China
| | - Chen Zhang
- The Key Laboratory of Carbon
Fiber and Functional Polymers, Ministry of Education, Beijing University of Chemical Technology, Beijing
100029, Peopleʼs Republic of China
| | - Zhongjie Du
- The Key Laboratory of Carbon
Fiber and Functional Polymers, Ministry of Education, Beijing University of Chemical Technology, Beijing
100029, Peopleʼs Republic of China
| | - Jianguo Mi
- State Key Laboratory of Organic−Inorganic
Composites, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
50
|
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
|