1
|
Fowles DJ, Palmer DS. Solvation entropy, enthalpy and free energy prediction using a multi-task deep learning functional in 1D-RISM. Phys Chem Chem Phys 2023; 25:6944-6954. [PMID: 36806875 DOI: 10.1039/d3cp00199g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
Simultaneous calculation of entropies, enthalpies and free energies has been a long-standing challenge in computational chemistry, partly because of the difficulty in obtaining estimates of all three properties from a single consistent simulation methodology. This has been particularly true for methods from the Integral Equation Theory of Molecular Liquids such as the Reference Interaction Site Model which have traditionally given large errors in solvation thermodynamics. Recently, we presented pyRISM-CNN, a combination of the 1 Dimensional Reference Interaction Site Model (1D-RISM) solver, pyRISM, with a deep learning based free energy functional, as a method of predicting solvation free energy (SFE). With this approach, a 40-fold improvement in prediction accuracy was delivered for a multi-solvent, multi-temperature dataset when compared to the standard 1D-RISM theory [Fowles et al., Digital Discovery, 2023, 2, 177-188]. Here, we report three further developments to the pyRISM-CNN methodology. Firstly, solvation free energies have been introduced for organic molecular ions in methanol or water solvent systems at 298 K, with errors below 4 kcal mol-1 obtained without the need for corrections or additional descriptors. Secondly, the number of solvents in the training data has been expanded from carbon tetrachloride, water and chloroform to now also include methanol. For neutral solutes, prediction errors nearing or below 1 kcal mol-1 are obtained for each organic solvent system at 298 K and water solvent systems at 273-373 K. Lastly, pyRISM-CNN was successfully applied to the simultaneous prediction of solvation enthalpy, entropy and free energy through a multi-task learning approach, with errors of 1.04, 0.98 and 0.47 kcal mol-1, respectively, for water solvent systems at 298 K.
Collapse
Affiliation(s)
- Daniel J Fowles
- Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, UK.
| | - David S Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, Scotland G1 1XL, UK.
| |
Collapse
|
2
|
Fan K, Zhang Y, Qiu Y, Zhang H. Impacts of targeting different hydration free energy references on the development of ion potentials. Phys Chem Chem Phys 2022; 24:16244-16262. [PMID: 35758314 DOI: 10.1039/d2cp01237e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H+). The choice of ΔG°(H+) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K+, Mg2+, Al3+, and Cl- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H+) for the cation and an increase in ΔG°(H+) for the anion. The increase in ΔG°(H+) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H+) value of -259.8 kcal mol-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H+) for Mg2+ model is better for a reasonable description of the metal binding network.
Collapse
Affiliation(s)
- Kun Fan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| |
Collapse
|
3
|
Waibl F, Kraml J, Hoerschinger VJ, Hofer F, Kamenik AS, Fernández-Quintero ML, Liedl KR. Grid inhomogeneous solvation theory for cross-solvation in rigid solvents. J Chem Phys 2022; 156:204101. [PMID: 35649837 DOI: 10.1063/5.0087549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Grid Inhomogeneous Solvation Theory (GIST) has proven useful to calculate localized thermodynamic properties of water around a solute. Numerous studies have leveraged this information to enhance structure-based binding predictions. We have recently extended GIST toward chloroform as a solvent to allow the prediction of passive membrane permeability. Here, we further generalize the GIST algorithm toward all solvents that can be modeled as rigid molecules. This restriction is inherent to the method and is already present in the inhomogeneous solvation theory. Here, we show that our approach can be applied to various solvent molecules by comparing the results of GIST simulations with thermodynamic integration (TI) calculations and experimental results. Additionally, we analyze and compare a matrix consisting of 100 entries of ten different solvent molecules solvated within each other. We find that the GIST results are highly correlated with TI calculations as well as experiments. For some solvents, we find Pearson correlations of up to 0.99 to the true entropy, while others are affected by the first-order approximation more strongly. The enthalpy-entropy splitting provided by GIST allows us to extend a recently published approach, which estimates higher order entropies by a linear scaling of the first-order entropy, to solvents other than water. Furthermore, we investigate the convergence of GIST in different solvents. We conclude that our extension to GIST reliably calculates localized thermodynamic properties for different solvents and thereby significantly extends the applicability of this widely used method.
Collapse
Affiliation(s)
- Franz Waibl
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| | - Johannes Kraml
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| | - Valentin J Hoerschinger
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| | - Florian Hofer
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| | - Anna S Kamenik
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| | - Monica L Fernández-Quintero
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| | - Klaus R Liedl
- Center for Molecular Biosciences Innsbruck, Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, Innsbruck, Austria
| |
Collapse
|
4
|
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
|
5
|
Qiu Y, Jiang Y, Zhang Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential. J Phys Chem B 2021; 125:13502-13518. [PMID: 34860517 DOI: 10.1021/acs.jpcb.1c09103] [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/28/2022]
Abstract
Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C4 term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C4 term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C4 term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.
Collapse
Affiliation(s)
- Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
6
|
Salassi S, Caselli L, Cardellini J, Lavagna E, Montis C, Berti D, Rossi G. A Martini Coarse Grained Model of Citrate-Capped Gold Nanoparticles Interacting with Lipid Bilayers. J Chem Theory Comput 2021; 17:6597-6609. [PMID: 34491056 PMCID: PMC8515808 DOI: 10.1021/acs.jctc.1c00627] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Indexed: 12/29/2022]
Abstract
Citrate capping is one of the most common strategies to achieve the colloidal stability of Au nanoparticles (NPs) with diameters ranging from a few to hundreds of nanometers. Citrate-capped Au nanoparticles (CNPs) represent a step of the synthesis of Au NPs with specific functionalities, as CNPs can be further functionalized via ligand-exchange reactions, leading to the replacement of citrate with other organic ligands. In vitro, CNPs are also used to address the fundamental aspects of NP-membrane interactions, as they can directly interact with cells or model cell membranes. Their affinity for the bilayer is again mediated by the exchange of citrate with lipid molecules. Here, we propose a new computational model of CNPs compatible with the coarse grained Martini force field. The model, which we develop and validate through an extensive comparison with new all-atom molecular dynamics (MD) simulations and UV-vis and Fourier transform infrared spectroscopy data, is aimed at the MD simulation of the interaction between citrate-capped NPs and model phosphatidylcholine lipid membranes. As a test application we show that, during the interaction between a single CNP and a flat planar 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer, the citrate coating is spontaneously replaced by lipids on the surface of Au NPs, while the NP size and shape determine the final structural configuration of the NP-bilayer complex.
Collapse
Affiliation(s)
- Sebastian Salassi
- Department
of Physics, University of Genoa, Via Dodecaneso 33, Genoa 16146, Italy
| | - Lucrezia Caselli
- Department
of Chemistry “Ugo Schiff”, University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
- CSGI,
Consorzio Sistemi a Grande Interfase and Department of Chemistry “Ugo
Schiff” University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
| | - Jacopo Cardellini
- Department
of Chemistry “Ugo Schiff”, University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
- CSGI,
Consorzio Sistemi a Grande Interfase and Department of Chemistry “Ugo
Schiff” University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
| | - Enrico Lavagna
- Department
of Physics, University of Genoa, Via Dodecaneso 33, Genoa 16146, Italy
| | - Costanza Montis
- Department
of Chemistry “Ugo Schiff”, University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
- CSGI,
Consorzio Sistemi a Grande Interfase and Department of Chemistry “Ugo
Schiff” University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
| | - Debora Berti
- Department
of Chemistry “Ugo Schiff”, University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
- CSGI,
Consorzio Sistemi a Grande Interfase and Department of Chemistry “Ugo
Schiff” University of Florence, Via della Lastruccia 3, Sesto Fiorentino, Florence 50019, Italy
| | - Giulia Rossi
- Department
of Physics, University of Genoa, Via Dodecaneso 33, Genoa 16146, Italy
| |
Collapse
|
7
|
Zhang Y, Jiang Y, Qiu Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Highly Charged Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4613-4629. [PMID: 34467756 DOI: 10.1021/acs.jcim.1c00723] [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/28/2022]
Abstract
Here, we developed nonbonded point charge models using a simple Lennard-Jones (LJ) 12-6 potential for highly charged metal cations (18 trivalent and 6 tetravalent ions) for use with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models simultaneously reproduce the hydration free energy (HFE) and ion-oxygen distance (IOD) in the first hydration shell with an error within 1 kcal/mol and 0.01 Å on average, respectively, and yield reasonable coordination numbers for most cations. Such performance is equivalent to the previously reported point charge models using a more complex 12-6-4 LJ-type potential, while the LJ R parameters of our models are much close to Shannon's revised effective ion radii than that of the 12-6-4 models. Our designed models overestimate the diffusion constants of several trivalent ions by 5-68%. The performance in predicting osmotic coefficients of trivalent chlorides in aqueous solution depends on the salt type. A calibration of cation-anion interacting LJ parameters reproduces the experimental osmotic coefficients of an AlCl3 solution at 0.2-3.0 mol/L. The effectiveness of our new models is further demonstrated by simulating a metalloprotein system with four force field/water combinations. This work facilitates accurate modeling of metal-containing systems by a variety of force fields and water models in aqueous solution.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
8
|
Zhang Y, Jiang Y, Peng J, Zhang H. Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4031-4044. [PMID: 34313132 DOI: 10.1021/acs.jcim.1c00580] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C4 term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M2+) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/Eb, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion-oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M2+ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion-water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park 16802, Pennsylvania, United States
| | - Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
9
|
Fowles DJ, Palmer DS, Guo R, Price SL, Mitchell JBO. Toward Physics-Based Solubility Computation for Pharmaceuticals to Rival Informatics. J Chem Theory Comput 2021; 17:3700-3709. [PMID: 33988381 PMCID: PMC8190954 DOI: 10.1021/acs.jctc.1c00130] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
![]()
We demonstrate that
physics-based calculations of intrinsic aqueous
solubility can rival cheminformatics-based machine learning predictions.
A proof-of-concept was developed for a physics-based approach via
a sublimation thermodynamic cycle, building upon previous work that
relied upon several thermodynamic approximations, notably the 2RT approximation, and limited conformational sampling. Here,
we apply improvements to our sublimation free-energy model with the
use of crystal phonon mode calculations to capture the contributions
of the vibrational modes of the crystal. Including these improvements
with lattice energies computed using the model-potential-based Ψmol method leads to accurate estimates of sublimation free
energy. Combining these with hydration free energies obtained from
either molecular dynamics free-energy perturbation simulations or
density functional theory calculations, solubilities comparable to
both experiment and informatics predictions are obtained. The application
to coronene, succinic acid, and the pharmaceutical desloratadine shows
how the methods must be adapted for the adoption of different conformations
in different phases. The approach has the flexibility to extend to
applications that cannot be covered by informatics methods.
Collapse
Affiliation(s)
- Daniel J Fowles
- 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
| | - Rui Guo
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Sarah L Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - John B O Mitchell
- EaStCHEM School of Chemistry and Biomedical Sciences Research Complex, University of St Andrews, St Andrews, Scotland KY16 9ST, U.K
| |
Collapse
|
10
|
Peng J, Zhang Y, Jiang Y, Zhang H. Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References. J Chem Inf Model 2021; 61:2981-2997. [PMID: 34080414 DOI: 10.1021/acs.jcim.1c00281] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg2+ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl2 solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg2+ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl2 solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.
Collapse
Affiliation(s)
- Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
11
|
Fossat M, Zeng X, Pappu RV. Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Side Chains of Amino Acids. J Phys Chem B 2021; 125:4148-4161. [PMID: 33877835 PMCID: PMC8154595 DOI: 10.1021/acs.jpcb.1c01073] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Revised: 04/07/2021] [Indexed: 02/07/2023]
Abstract
Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged side chains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle-based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the side chains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD-derived trends in free energies of hydration. To this end, we performed temperature-dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of ca. -259 kcal/mol. Analysis of temperature-dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-.
Collapse
Affiliation(s)
| | | | - Rohit V. Pappu
- Department of Biomedical Engineering
and Center for Science & Engineering of Living Systems, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| |
Collapse
|
12
|
Tielker N, Eberlein L, Hessler G, Schmidt KF, Güssregen S, Kast SM. Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges? J Comput Aided Mol Des 2021; 35:453-472. [PMID: 33079358 PMCID: PMC8018924 DOI: 10.1007/s10822-020-00347-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Accepted: 09/26/2020] [Indexed: 01/26/2023]
Abstract
Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Gerhard Hessler
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - K Friedemann Schmidt
- R&D Preclinical Safety, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - Stefan Güssregen
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany.
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
13
|
Yadav AK, Bandyopadhyay P, Urbic T, Dill KA. Analytical 2-Dimensional Model of Nonpolar and Ionic Solvation in Water. J Phys Chem B 2021; 125:1861-1873. [PMID: 33539097 PMCID: PMC7958497 DOI: 10.1021/acs.jpcb.0c10329] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A goal in computational chemistry is computing hydration free energies of nonpolar and charged solutes accurately, but with much greater computational speeds than in today's explicit-water simulations. Here, we take one step in that direction: a simple model of solvating waters that is analytical and thus essentially instantaneous to compute. Each water molecule is a 2-dimensional dipolar hydrogen-bonding disk that interacts around small circular solutes with different nonpolar and charge interactions. The model gives good qualitative agreement with experiments. As a function of the solute radius, it gives the solvation free energy, enthalpy and entropy as a function of temperature for the inert gas series Ne, Ar, Kr, and Xe. For anions and cations, it captures relatively well the trends versus ion radius. This approach should be readily generalizable to three dimensions.
Collapse
Affiliation(s)
- Ajeet Kumar Yadav
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, SI-1000, Ljubljana, Slovenia
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, New York, New York 11794, United States
- Department of Physics and Astronomy, Stony Brook University, New York, New York 11794, United States
- Department of Chemistry, Stony Brook University, New York, New York 11794, United States
| |
Collapse
|
14
|
Han Y, Zhu L, Yao Y, Shi X, Zhang Y, Xiao H, Chen X. Strong bases behave as weak bases in nanoscale chemical environments: implication in humidity-swing CO 2 air capture. Phys Chem Chem Phys 2021; 23:14811-14817. [PMID: 34212971 DOI: 10.1039/d1cp01121a] [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/21/2022]
Abstract
Hydration of ions/molecules in nanometer-sized clusters or nanoscopic pores is ubiquitous and plays a key role in many chemical and physical systems. In this work, guanidine-H2O reactions with n = 1-8 water molecules were systematically studied by ab initio methods. The result suggests that the reduced availability of water molecules greatly inhibits the strong base guanidine from producing OH-. That is, guanidine exhibits the behavior of a weak bases in low-humidity nanoscale environments. Intriguingly, this effect is not limited to guanidine but could be applied to other strong bases. Furthermore, we demonstrate that the direction of guanidine-CO2 reactions can be controlled by changing the number of water molecules present, which in turn responds to the humidity change in air. These findings not only shed some light on unconventional chemical reactions of strong bases in atmospheric clusters and on solid porous surfaces, but also provide insights into the development of guanidine-based CO2 air-capture sorbents.
Collapse
Affiliation(s)
- Yingying Han
- Shaanxi Institute of Energy and Chemical Engineering, School of Chemical Engineering, Northwest University, Xi'an, 710069, China.
| | - Liangliang Zhu
- Shaanxi Institute of Energy and Chemical Engineering, School of Chemical Engineering, Northwest University, Xi'an, 710069, China.
| | - Yutong Yao
- Earth Engineering Center, Center for Advanced Materials for Energy and Environment, Department of Earth and Environmental Engineering, Columbia University, New York, NY10027, USA
| | - Xiaoyang Shi
- Earth Engineering Center, Center for Advanced Materials for Energy and Environment, Department of Earth and Environmental Engineering, Columbia University, New York, NY10027, USA
| | - Yayun Zhang
- State Key Laboratory of Chemical Engineering, East China University of Science and Technology, Shanghai, 200237, China
| | - Hang Xiao
- Shaanxi Institute of Energy and Chemical Engineering, School of Chemical Engineering, Northwest University, Xi'an, 710069, China.
| | - Xi Chen
- Earth Engineering Center, Center for Advanced Materials for Energy and Environment, Department of Earth and Environmental Engineering, Columbia University, New York, NY10027, USA
| |
Collapse
|
15
|
Miyata T. Sigma enlarging bridge function for heteronuclear Lennard-Jones diatomic solute solvated in a Lennard-Jones monatomic solvent in terms of the parameter transferability. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137777] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
16
|
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
|
17
|
Electric field assisted desalination of water using B- and N-doped-graphene sheets: A non-equilibrium molecular dynamics study. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.112574] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
18
|
Luukkonen S, Levesque M, Belloni L, Borgis D. Hydration free energies and solvation structures with molecular density functional theory in the hypernetted chain approximation. J Chem Phys 2020; 152:064110. [DOI: 10.1063/1.5142651] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Sohvi Luukkonen
- 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, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Luc Belloni
- LIONS, NIMBE, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, 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, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
19
|
Tielker N, Tomazic D, Eberlein L, Güssregen S, Kast SM. The SAMPL6 challenge on predicting octanol-water partition coefficients from EC-RISM theory. J Comput Aided Mol Des 2020; 34:453-461. [PMID: 31981015 PMCID: PMC7125249 DOI: 10.1007/s10822-020-00283-4] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Accepted: 01/08/2020] [Indexed: 12/14/2022]
Abstract
Results are reported for octanol–water partition coefficients (log P) of the neutral states of drug-like molecules provided during the SAMPL6 (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenge from applying the “embedded cluster reference interaction site model” (EC-RISM) as a solvation model for quantum-chemical calculations. Following the strategy outlined during earlier SAMPL challenges we first train 1- and 2-parameter water-free (“dry”) and water-saturated (“wet”) models for n-octanol solvation Gibbs energies with respect to experimental values from the “Minnesota Solvation Database” (MNSOL), yielding a root mean square error (RMSE) of 1.5 kcal mol−1 for the best-performing 2-parameter wet model, while the optimal water model developed for the pKa part of the SAMPL6 challenge is kept unchanged (RMSE 1.6 kcal mol−1 for neutral compounds from a model trained on both neutral and ionic species). Applying these models to the blind prediction set yields a log P RMSE of less than 0.5 for our best model (2-parameters, wet). Further analysis of our results reveals that a single compound is responsible for most of the error, SM15, without which the RMSE drops to 0.2. Since this is the only compound in the challenge dataset with a hydroxyl group we investigate other alcohols for which Gibbs energy of solvation data for both water and n-octanol are available in the MNSOL database to demonstrate a systematic cause of error and to discuss strategies for improvement.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Daniel Tomazic
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, R&D Integrated Drug Discovery, 65926, Frankfurt am Main, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
20
|
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
|
21
|
Lim H, Jung Y. Delfos: deep learning model for prediction of solvation free energies in generic organic solvents. Chem Sci 2019; 10:8306-8315. [PMID: 32110289 PMCID: PMC7017869 DOI: 10.1039/c9sc02452b] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Accepted: 08/19/2019] [Indexed: 11/21/2022] Open
Abstract
Prediction of aqueous solubilities or hydration free energies is an extensively studied area in machine learning applications in chemistry since water is the sole solvent in the living system. However, for non-aqueous solutions, few machine learning studies have been undertaken so far despite the fact that the solvation mechanism plays an important role in various chemical reactions. Here, we introduce Delfos (deep learning model for solvation free energies in generic organic solvents), which is a novel, machine-learning-based QSPR method which predicts solvation free energies for various organic solute and solvent systems. A novelty of Delfos involves two separate solvent and solute encoder networks that can quantify structural features of given compounds via word embedding and recurrent layers, augmented with the attention mechanism which extracts important substructures from outputs of recurrent neural networks. As a result, the predictor network calculates the solvation free energy of a given solvent-solute pair using features from encoders. With the results obtained from extensive calculations using 2495 solute-solvent pairs, we demonstrate that Delfos not only has great potential in showing accuracy comparable to that of the state-of-the-art computational chemistry methods, but also offers information about which substructures play a dominant role in the solvation process.
Collapse
Affiliation(s)
- Hyuntae Lim
- Department of Chemistry , Seoul National University , Seoul 08826 , Korea . ;
| | - YounJoon Jung
- Department of Chemistry , Seoul National University , Seoul 08826 , Korea . ;
| |
Collapse
|
22
|
An assessment of the sigma enlarging bridge function for a Lennard-Jones solution using a solvent-solvent correlation function obtained from molecular dynamics simulation. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111167] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
23
|
Tanimoto S, Yoshida N, Yamaguchi T, Ten-no SL, Nakano H. Effect of Molecular Orientational Correlations on Solvation Free Energy Computed by Reference Interaction Site Model Theory. J Chem Inf Model 2019; 59:3770-3781. [DOI: 10.1021/acs.jcim.9b00330] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Shoichi Tanimoto
- Department of Chemistry, Graduate School of Science, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Tsuyoshi Yamaguchi
- Graduate School of Engineering, Nagoya University, Chikusa-ku, Nagoya 464-8603, Japan
| | - Seiichiro L. Ten-no
- Graduate School of Science, Technology, and Innovation, Kobe University, Nada-ku, Kobe 657-8501, Japan
| | - Haruyuki Nakano
- Department of Chemistry, Graduate School of Science, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan
| |
Collapse
|
24
|
Pliego JR, Riveros JM. Hybrid discrete‐continuum solvation methods. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1440] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Josefredo R. Pliego
- Departamento de Ciências Naturais Universidade Federal de São João del‐Rei São João del‐Rei Brazil
| | - Jose M. Riveros
- Departamento de Química Fundamental Instituto de Química, Universidade de São Paulo São Paulo Brazil
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
Ganyecz Á, Kállay M, Csontos J. Thermochemistry of Uracil, Thymine, Cytosine, and Adenine. J Phys Chem A 2019; 123:4057-4067. [DOI: 10.1021/acs.jpca.9b02061] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ádám Ganyecz
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, PO Box 91, H-1521, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, PO Box 91, H-1521, Hungary
| | - József Csontos
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, PO Box 91, H-1521, Hungary
| |
Collapse
|
28
|
Hutchinson ST, Kobayashi R. Solvent-Specific Featurization for Predicting Free Energies of Solvation through Machine Learning. J Chem Inf Model 2019; 59:1338-1346. [DOI: 10.1021/acs.jcim.8b00901] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Samuel T. Hutchinson
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Rika Kobayashi
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
- ANU Supercomputer Facility, Leonard Huxley Bldg 56, Mills Rd, Canberra, ACT 2601, Australia
| |
Collapse
|
29
|
Sáenz-Tavera IDC, Rosas-García VM. Ab initio calculations and reduced density gradient analyses of the structure and energetics of hydrated calcium fluoride and calcium carbonate. Phys Chem Chem Phys 2019; 21:5744-5758. [DOI: 10.1039/c8cp06353b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
We studied microhydrated calcium fluoride, calcium carbonate and their ions at the MP2/6-311++G** level of theory and found that water–water non-covalent interactions destabilize the solvation shell, and are compensated by cooperative hydrogen bonds.
Collapse
Affiliation(s)
| | - Victor M. Rosas-García
- Universidad Autónoma de Nuevo León
- UANL
- Facultad de Ciencias Químicas
- Ave. Universidad S/N
- Cd. Universitaria
| |
Collapse
|
30
|
Milne AW, Jorge M. Polarization Corrections and the Hydration Free Energy of Water. J Chem Theory Comput 2018; 15:1065-1078. [DOI: 10.1021/acs.jctc.8b01115] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Andrew W. Milne
- Department of Chemical and Process Engineering, University of Strathclyde, Glasgow G1 1XQ, United Kingdom
| | - Miguel Jorge
- Department of Chemical and Process Engineering, University of Strathclyde, Glasgow G1 1XQ, United Kingdom
| |
Collapse
|
31
|
Ekimoto T, Yamane T, Ikeguchi M. Elimination of Finite-Size Effects on Binding Free Energies via the Warp-Drive Method. J Chem Theory Comput 2018; 14:6544-6559. [DOI: 10.1021/acs.jctc.8b00280] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Toru Ekimoto
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Tsutomu Yamane
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
- RIKEN Medical Sciences Innovation Hub Program, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| |
Collapse
|
32
|
Sosnin S, Misin M, Palmer DS, Fedorov MV. 3D matters! 3D-RISM and 3D convolutional neural network for accurate bioaccumulation prediction. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:32LT03. [PMID: 29964270 DOI: 10.1088/1361-648x/aad076] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
In this work, we present a new method for predicting complex physical-chemical properties of organic molecules. The approach utilizes 3D convolutional neural network (ActivNet4) that uses solvent spatial distributions around solutes as input. These spatial distributions are obtained by a molecular theory called three-dimensional reference interaction site model. We have shown that the method allows one to achieve a good accuracy of prediction of bioconcentration factor which is difficult to predict by direct application of methods of molecular theory or simulations. Our research demonstrates that combination of molecular theories with modern machine learning approaches can be effectively used for predicting properties that are otherwise inaccessible to purely theory-based models.
Collapse
Affiliation(s)
- Sergey Sosnin
- Center for Computational and Data-intensive Science and Engineering, Skolkovo Institute of Science and Technology, Nobelya Ulitsa 3 Moscow, 121205, Russia
| | | | | | | |
Collapse
|
33
|
Tielker N, Eberlein L, Güssregen S, Kast SM. The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory. J Comput Aided Mol Des 2018; 32:1151-1163. [DOI: 10.1007/s10822-018-0140-z] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2018] [Accepted: 07/28/2018] [Indexed: 12/15/2022]
|
34
|
Hofer TS, Hünenberger PH. Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration. J Chem Phys 2018; 148:222814. [DOI: 10.1063/1.5000799] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Centre for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | | |
Collapse
|
35
|
Zhang H, Yin C, Jiang Y, van der Spoel D. Force Field Benchmark of Amino Acids: I. Hydration and Diffusion in Different Water Models. J Chem Inf Model 2018; 58:1037-1052. [DOI: 10.1021/acs.jcim.8b00026] [Citation(s) in RCA: 61] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Chunhua Yin
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Yang Jiang
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, Beijing 100029, China
| | - David van der Spoel
- Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box
596, SE-75124 Uppsala, Sweden
| |
Collapse
|
36
|
Epsztein R, Shaulsky E, Dizge N, Warsinger DM, Elimelech M. Role of Ionic Charge Density in Donnan Exclusion of Monovalent Anions by Nanofiltration. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2018; 52:4108-4116. [PMID: 29510032 DOI: 10.1021/acs.est.7b06400] [Citation(s) in RCA: 109] [Impact Index Per Article: 18.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The main objective of this study is to examine how the charge densities of four monovalent anions-fluoride (F-), chloride (Cl-), bromide (Br-), and nitrate (NO3-)-influence their Donnan (charge) exclusion by a charged nanofiltration (NF) membrane. We systematically studied the rejection behavior of ternary ion solutions containing sodium cation (Na+) and two of the monovalent anions as a function of the pH with a polyamide NF membrane. In the solutions containing F- and Cl- or F- and Br-, F- rejection was higher than Cl- or Br- rejection only when the solution pH was higher than 5.5, suggesting that F- (which has a higher charge density) was repelled more strongly by the negatively charged membrane. The order of change in the activation energy for the transport of the four anions through the polyamide membrane as a response to the increase of the membrane negative charge was the following: F- > Cl- > NO3- > Br-. This order corroborates our main hypothesis that an anion with a smaller ionic radius, and hence a higher charge density, is more affected by the Donnan (charge)-exclusion mechanism in NF. We conclude with a proposed mechanism for the role of ionic charge density in the rejection of monovalent anions in NF.
Collapse
Affiliation(s)
- Razi Epsztein
- Department of Chemical and Environmental Engineering , Yale University , New Haven , Connecticut 06520-8286 , United States
| | - Evyatar Shaulsky
- Department of Chemical and Environmental Engineering , Yale University , New Haven , Connecticut 06520-8286 , United States
| | - Nadir Dizge
- Department of Chemical and Environmental Engineering , Yale University , New Haven , Connecticut 06520-8286 , United States
- Department of Environmental Engineering , Mersin University , Mersin 33343 , Turkey
| | - David M Warsinger
- Department of Chemical and Environmental Engineering , Yale University , New Haven , Connecticut 06520-8286 , United States
| | - Menachem Elimelech
- Department of Chemical and Environmental Engineering , Yale University , New Haven , Connecticut 06520-8286 , United States
| |
Collapse
|
37
|
Qiao B, Skanthakumar S, Soderholm L. Comparative CHARMM and AMOEBA Simulations of Lanthanide Hydration Energetics and Experimental Aqueous-Solution Structures. J Chem Theory Comput 2018; 14:1781-1790. [DOI: 10.1021/acs.jctc.7b01018] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Baofu Qiao
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - S. Skanthakumar
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - L. Soderholm
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
38
|
Young JM, Panagiotopoulos AZ. System-Size Dependence of Electrolyte Activity Coefficients in Molecular Simulations. J Phys Chem B 2018; 122:3330-3338. [DOI: 10.1021/acs.jpcb.7b09861] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jeffrey M. Young
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
39
|
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
|
40
|
Kuhn AB, Kube S, Karow-Zwick AR, Seeliger D, Garidel P, Blech M, Schäfer LV. Improved Solution-State Properties of Monoclonal Antibodies by Targeted Mutations. J Phys Chem B 2017; 121:10818-10827. [DOI: 10.1021/acs.jpcb.7b09126] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Alexander B. Kuhn
- Theoretical
Chemistry, Ruhr-University Bochum, D-44780 Bochum, Germany
| | | | | | | | | | | | - Lars V. Schäfer
- Theoretical
Chemistry, Ruhr-University Bochum, D-44780 Bochum, Germany
| |
Collapse
|
41
|
Zhang H, Jiang Y, Yan H, Cui Z, Yin C. Comparative Assessment of Computational Methods for Free Energy Calculations of Ionic Hydration. J Chem Inf Model 2017; 57:2763-2775. [DOI: 10.1021/acs.jcim.7b00485] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Haiyang Zhang
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Hai Yan
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Ziheng Cui
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Box 53, 100029 Beijing, China
| | - Chunhua Yin
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
42
|
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
|
43
|
Pliego JR. Cluster expansion of the solvation free energy difference: Systematic improvements in the solvation of single ions. J Chem Phys 2017; 147:034104. [DOI: 10.1063/1.4993770] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Josefredo R. Pliego
- Departamento de Ciências Naturais, Universidade Federal de São João del-Rei, 36301-160 São João del-Rei, MG, Brazil
| |
Collapse
|
44
|
Zhang H, Jiang Y, Yan H, Yin C, Tan T, van der Spoel D. Free-Energy Calculations of Ionic Hydration Consistent with the Experimental Hydration Free Energy of the Proton. J Phys Chem Lett 2017; 8:2705-2712. [PMID: 28561580 DOI: 10.1021/acs.jpclett.7b01125] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Computational free-energy correction strategies and the choice of experimental proton hydration free energy, ΔGs*(H+), are analyzed to investigate the apparent controversy in experimental thermodynamics of ionic hydration. Without corrections, the hydration free-energy (ΔGhyd) calculations match experiments with ΔGs*(H+) = -1064 kJ/mol as reference. Using the Galvani surface potential the resulting (real) ΔGhyd are consistent with ΔGs*(H+) = -1098 kJ/mol. When applying, in an ad hoc manner, the discrete solvent correction, ΔGhyd matching the "consensus" ΔGs*(H+) of -1112 kJ/mol are obtained. This analysis rationalizes reports on ΔGhyd calculations for ions using different experimental references. For neutral amino acid side chains ΔGhyd are independent of the water model, whereas there are large differences in ΔGhyd due to the water model for charged species, suggesting that long-range ordering of water around ions yields an important contribution to the ΔGhyd. These differences are reduced significantly when applying consistent corrections, but to obtain the most accurate results it is recommended to use the water model belonging to the force field.
Collapse
Affiliation(s)
- Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing , 100083 Beijing, China
| | - Yang Jiang
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology , Box 53, 100029 Beijing, China
| | - Hai Yan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing , 100083 Beijing, China
| | - Chunhua Yin
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing , 100083 Beijing, China
| | - Tianwei Tan
- Beijing Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology , Box 53, 100029 Beijing, China
| | - David van der Spoel
- Uppsala Center for Computational Chemistry, Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University , Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
45
|
Misin M, Vainikka PA, Fedorov MV, Palmer DS. Salting-out effects by pressure-corrected 3D-RISM. J Chem Phys 2016; 145:194501. [DOI: 10.1063/1.4966973] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Maksim Misin
- Department of Physics, SUPA, University of Strathclyde, 107 Rottenrow, Glasgow G4 0NG, United Kingdom
| | - Petteri A. Vainikka
- Department of Pure and Applied Chemistry, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, United Kingdom
| | - Maxim V. Fedorov
- Department of Physics, SUPA, University of Strathclyde, 107 Rottenrow, Glasgow G4 0NG, United Kingdom
- Skolkovo Institute of Science and Technology, 3 Nobel Street, Moscow 143026, Russian Federation
| | - David S. Palmer
- Department of Pure and Applied Chemistry, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, United Kingdom
| |
Collapse
|
46
|
Kobryn AE, Gusarov S, Kovalenko A. A closure relation to molecular theory of solvation for macromolecules. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:404003. [PMID: 27549008 DOI: 10.1088/0953-8984/28/40/404003] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We propose a closure to the integral equations of molecular theory of solvation, particularly suitable for polar and charged macromolecules in electrolyte solution. This includes such systems as oligomeric polyelectrolytes at a finite concentration in aqueous and various non-aqueous solutions, as well as drug-like compounds in solution. The new closure by Kobryn, Gusarov, and Kovalenko (KGK closure) imposes the mean spherical approximation (MSA) almost everywhere in the solvation shell but levels out the density distribution function to zero (with the continuity at joint boundaries) inside the repulsive core and in the spatial regions of strong density depletion emerging due to molecular associative interactions. Similarly to MSA, the KGK closure reduces the problem to a linear equation for the direct correlation function which is predefined analytically on most of the solvation shells and has to be determined numerically on a relatively small (three-dimensional) domain of strong depletion, typically within the repulsive core. The KGK closure leads to the solvation free energy in the form of the Gaussian fluctuation (GF) functional. We first test the performance of the KGK closure coupled to the reference interaction site model (RISM) integral equations on the examples of Lennard-Jones liquids, polar and nonpolar molecular solvents, including water, and aqueous solutions of simple ions. The solvation structure, solvation chemical potential, and compressibility obtained from RISM with the KGK closure favorably compare to the results of the hypernetted chain (HNC) and Kovalenko-Hirata (KH) closures, including their combination with the GF solvation free energy. We then use the KGK closure coupled to RISM to obtain the solvation structure and thermodynamics of oligomeric polyelectrolytes and drug-like compounds at a finite concentration in electrolyte solution, for which no convergence is obtained with other closures. For comparison, we calculate their solvation structure from molecular dynamics (MD) simulations. We further couple the 3D-RISM integral equation with the 3D-version of the KGK closure, and solve it for molecular mixtures as well as oligomeric polyelectrolytes and drug-like molecules in electrolyte solutions.
Collapse
Affiliation(s)
- Alexander E Kobryn
- National Institute for Nanotechnology, National Research Council Canada, 11421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada
| | | | | |
Collapse
|
47
|
Johnson J, Case DA, Yamazaki T, Gusarov S, Kovalenko A, Luchko T. Small molecule hydration energy and entropy from 3D-RISM. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:344002. [PMID: 27367817 PMCID: PMC5118872 DOI: 10.1088/0953-8984/28/34/344002] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Implicit solvent models offer an attractive way to estimate the effects of a solvent environment on the properties of small or large solutes without the complications of explicit simulations. One common test of accuracy is to compute the free energy of transfer from gas to liquid for a variety of small molecules, since many of these values have been measured. Studies of the temperature dependence of these values (i.e. solvation enthalpies and entropies) can provide additional insights into the performance of implicit solvent models. Here, we show how to compute temperature derivatives of hydration free energies for the 3D-RISM integral equation approach. We have computed hydration free energies of 1123 small drug-like molecules (both neutral and charged). Temperature derivatives were also used to calculate hydration energies and entropies of 74 of these molecules (both neutral and charged) for which experimental data is available. While direct results have rather poor agreement with experiment, we have found that several previously proposed linear hydration free energy correction schemes give good agreement with experiment. These corrections also provide good agreement for hydration energies and entropies though simple extensions are required in some cases.
Collapse
Affiliation(s)
- J Johnson
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854
| | - D A Case
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854
| | - T Yamazaki
- Vancouver Prostate Centre, 2660 Oak Street, Vancouver, BC, V6H 3Z6, Canada
| | - S Gusarov
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Dr., Edmonton, AB, T6G 2M9, Canada
| | - A Kovalenko
- National Institute for Nanotechnology, National Research Council of Canada, 11421 Saskatchewan Dr., Edmonton, AB, T6G 2M9, Canada
- Department of Mechanical Engineering, University of Alberta, 10-203 Donadeo Innovation Centre for Engineering, 9211-116 Str., Edmonton, AB, T6G 1H9, Canada
| | - T Luchko
- Department of Physics and Astronomy, California State University, Northridge, CA 91330
| |
Collapse
|
48
|
SAMPL5: 3D-RISM partition coefficient calculations with partial molar volume corrections and solute conformational sampling. J Comput Aided Mol Des 2016; 30:1115-1127. [PMID: 27585474 DOI: 10.1007/s10822-016-9947-7] [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: 07/06/2016] [Accepted: 08/18/2016] [Indexed: 12/15/2022]
Abstract
Implicit solvent methods for classical molecular modeling are frequently used to provide fast, physics-based hydration free energies of macromolecules. Less commonly considered is the transferability of these methods to other solvents. The Statistical Assessment of Modeling of Proteins and Ligands 5 (SAMPL5) distribution coefficient dataset and the accompanying explicit solvent partition coefficient reference calculations provide a direct test of solvent model transferability. Here we use the 3D reference interaction site model (3D-RISM) statistical-mechanical solvation theory, with a well tested water model and a new united atom cyclohexane model, to calculate partition coefficients for the SAMPL5 dataset. The cyclohexane model performed well in training and testing ([Formula: see text] for amino acid neutral side chain analogues) but only if a parameterized solvation free energy correction was used. In contrast, the same protocol, using single solute conformations, performed poorly on the SAMPL5 dataset, obtaining [Formula: see text] compared to the reference partition coefficients, likely due to the much larger solute sizes. Including solute conformational sampling through molecular dynamics coupled with 3D-RISM (MD/3D-RISM) improved agreement with the reference calculation to [Formula: see text]. Since our initial calculations only considered partition coefficients and not distribution coefficients, solute sampling provided little benefit comparing against experiment, where ionized and tautomer states are more important. Applying a simple [Formula: see text] correction improved agreement with experiment from [Formula: see text] to [Formula: see text], despite a small number of outliers. Better agreement is possible by accounting for tautomers and improving the ionization correction.
Collapse
|
49
|
Tielker N, Tomazic D, Heil J, Kloss T, Ehrhart S, Güssregen S, Schmidt KF, Kast SM. The SAMPL5 challenge for embedded-cluster integral equation theory: solvation free energies, aqueous pK a, and cyclohexane-water log D. J Comput Aided Mol Des 2016; 30:1035-1044. [PMID: 27554666 DOI: 10.1007/s10822-016-9939-7] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Accepted: 08/17/2016] [Indexed: 12/01/2022]
Abstract
We predict cyclohexane-water distribution coefficients (log D 7.4) for drug-like molecules taken from the SAMPL5 blind prediction challenge by the "embedded cluster reference interaction site model" (EC-RISM) integral equation theory. This task involves the coupled problem of predicting both partition coefficients (log P) of neutral species between the solvents and aqueous acidity constants (pK a) in order to account for a change of protonation states. The first issue is addressed by calibrating an EC-RISM-based model for solvation free energies derived from the "Minnesota Solvation Database" (MNSOL) for both water and cyclohexane utilizing a correction based on the partial molar volume, yielding a root mean square error (RMSE) of 2.4 kcal mol-1 for water and 0.8-0.9 kcal mol-1 for cyclohexane depending on the parametrization. The second one is treated by employing on one hand an empirical pK a model (MoKa) and, on the other hand, an EC-RISM-derived regression of published acidity constants (RMSE of 1.5 for a single model covering acids and bases). In total, at most 8 adjustable parameters are necessary (2-3 for each solvent and two for the pK a) for training solvation and acidity models. Applying the final models to the log D 7.4 dataset corresponds to evaluating an independent test set comprising other, composite observables, yielding, for different cyclohexane parametrizations, 2.0-2.1 for the RMSE with the first and 2.2-2.8 with the combined first and second SAMPL5 data set batches. Notably, a pure log P model (assuming neutral species only) performs statistically similarly for these particular compounds. The nature of the approximations and possible perspectives for future developments are discussed.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Daniel Tomazic
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Thomas Kloss
- IPhT, L'Orme des Merisiers, CEA-Saclay, 91191, Gif-sur-Yvette, France
| | | | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926, Frankfurt am Main, Germany
| | - K Friedemann Schmidt
- Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926, Frankfurt am Main, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
50
|
Horta BAC, Merz PT, Fuchs PFJ, Dolenc J, Riniker S, Hünenberger PH. A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set. J Chem Theory Comput 2016; 12:3825-50. [DOI: 10.1021/acs.jctc.6b00187] [Citation(s) in RCA: 88] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bruno A. C. Horta
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Pascal T. Merz
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Patrick F. J. Fuchs
- Institut Jacques Monod, UMR 7592 CNRS, Université Paris-Diderot, Sorbonne Paris Cité, F-75205 Paris, France
| | - Jozica Dolenc
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
- Chemistry,
Biology and Pharmacy Information Center, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory
of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | | |
Collapse
|