Improved estimates for hydration free energy obtained by the reference interaction site model

Solvation entropy, enthalpy and free energy prediction using a multi-task deep learning functional in 1D-RISM

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.

Effect of Conformational Equilibrium on Solvation Properties of 1,2-DCE in Water: A Solvation Thermodynamics and 3D-RISM Study

The contributions of the enthalpy and entropy of solvation for the study of chemical and biological systems are important in the prediction, interpretation, and manipulation of these processes. The relation between solvation Gibbs energies, enthalpies, and entropies of solvation, and their rigorous relation with the conformational equilibrium, are derived for the first time and applied with a computational method, in accordance with the Solvation Thermodynamics previous results, to 1,2-dichloroethane solvation in water. The rigid conformer calculations in solution were performed by using PC+/3D-RISM approach, with the conformational averaged results for enthalpy and solvation Gibbs energy reproducing the experimental results quite successfully. A qualitative agreement in the entropy of solvation predictions was also observed.

Adsorption Energetics of Amino Acid Analogs on Polymer/Water Interfaces Studied by All-Atom Molecular Dynamics Simulation and a Theory of Solutions

Energetics of adsorption was addressed with all-atom molecular dynamics simulation on the interfaces of poly(2-methoxyethyl acrylate) (PMEA), poly(methyl methacrylate) (PMMA), and poly(butyl acrylate) (PBA) with water. A wide variety of adsorbate solutes were examined, and the free energy of adsorption was computed with the method of energy representation. It was found that the adsorption free energy was favorable (negative) for all the combinations of solute and polymer, and among PMEA, PMMA, and PBA, the strongest adsorption was observed on PMMA for the hydrophobic solutes and on PMEA for the hydrophilic ones. According to the decomposition of the adsorption free energy into the contributions from polymer and water, it was seen that the polymer contribution is larger in magnitude with the solute size. The total free energy of adsorption was correlated well with the solvation free energy in bulk water only for hydrophobic solutes. The roles of the intermolecular interaction components such as electrostatic, van der Waals, and excluded-volume were further studied, and the electrostatic component was influential only in determining the polymer dependences of the adsorption propensities of hydrophilic solutes. The extent of adsorption was shown to be ranked by the van der Waals component in the solute-polymer interaction separately over the hydrophilic and hydrophobic solutes, with the excluded-volume effect from water pointed out to also drive the adsorption.

Grid inhomogeneous solvation theory for cross-solvation in rigid solvents

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.

Solvation energetics of proteins and their aggregates analyzed by all-atom molecular dynamics simulations and the energy-representation theory of solvation

Solvation is a controlling factor for the structure and function of proteins. This article addresses the effects of solvation from an energetic perspective for the fluctuations and cosolvent-induced changes in protein structures and the equilibrium of aggregate formation for a peptide. A theoretical framework to analyze the solvation effects with an explicit solvent is introduced by adopting the energy-representation theory of solvation, and the connection of the solvation free energy to the protein structure and the aggregation tendency is quantitatively described in combination with all-atom molecular dynamics simulations. The interaction components that govern the solvation effects on the structural variations of proteins are further identified through correlation analysis, and a computational scheme to assess the shift of an aggregation equilibrium due to the addition of a cosolvent is provided.

Water Dissolved in a Variety of Polymers Studied by Molecular Dynamics Simulation and a Theory of Solutions

The performance of a polymer medium as a separation membrane is determined by the dissolution free energy ΔG and diffusion coefficient D of the permeant. In this work, ΔG and D of water are investigated with all-atom molecular dynamics simulation in a wide variety of polymer species in the amorphous state. The computed ΔG is shown to agree well with the experimental value for linear homopolymers, and the degrees of polymerization of the homopolymers do not affect ΔG when they are beyond ∼10. The copolymers of ethylene-vinylidene difluoride, ethylene-vinyl acetate, and ethylene-acrylamide are then examined by changing the repeating patterns of the constituent monomers in both the periodic and graft forms. It is found that ΔG is determined primarily by the overall compositions of the monomers and is not affected by the copolymerization topology (periodic or graft). The hydrophobicity of the copolymer is enhanced, furthermore, when the hydrophobicity and hydrophilicity of the ethylene and non-ethylene parts are well contrasted and those parts are fragmented along the polymer chain. According to the computed D, the diffusivity of water tends to be larger when the (co)polymer is more hydrophobic and ΔG is more positive. D is actually seen to vary by orders of magnitude with the polymer structures, while the effect of the polymer species on the water permeation is stronger for ΔG than for D.

Renormalized site density functional theory for models of ion hydration

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.

Accurate prediction of hydration free energies and solvation structures using molecular density functional theory with a simple bridge functional

This paper assesses the ability of molecular density functional theory to predict efficiently and accurately the hydration free energies of molecular solutes and the surrounding microscopic water structure. A wide range of solutes were investigated, including hydrophobes, water as a solute, and the FreeSolv database containing 642 drug-like molecules having a variety of shapes and sizes. The usual second-order approximation of the theory is corrected by a third-order, angular-independent bridge functional. The overall functional is parameter-free in the sense that the only inputs are bulk water properties, independent of the solutes considered. These inputs are the direct correlation function, compressibility, liquid-gas surface tension, and excess chemical potential of the solvent. Compared to molecular simulations with the same force field and the same fixed solute geometries, the present theory is shown to describe accurately the solvation free energy and structure of both hydrophobic and hydrophilic solutes. Overall, the method yields a precision of order 0.5 k_BT for the hydration free energies of the FreeSolv database, with a computer speedup of 3 orders of magnitude. The theory remains to be improved for a better description of the H-bonding structure and the hydration free energy of charged solutes.

Cavity Particle in Aqueous Solution with a Hydrophobic Solute: Structure, Energetics, and Functionals

Endpoints density functional theory (DFT) provides a framework for calculating the excess chemical potential of a solute in solution using solvent distribution functions obtained from both physical endpoints of a hypothetical charging process which transforms the solvent density from that of the pure liquid to the solution state. In this work, the endpoints DFT equations are formulated in terms of the indirect (solvent-mediated) contribution ω(x) to the solute-solvent potential of mean force, and their connections are established with the conventional DFT expressions which are based on the use of direct correlation functions. ω actually corresponds to the free-energy cost to move a cavity particle (a tagged solvent molecule which interacts with the other solvent molecules but not the solute) from the bulk to the configuration x of a solvent molecule relative to the solute and is a suitable variable to describe the solvent effects on the solute-solvent interactions. HNC and PY type approximations are then used to integrate the DFT charging integral involved in the exact expression for the excess chemical potential. With these approximations, molecular simulations are to be performed at the two endpoints of solute insertion: pure solvent without the solute and the solution system with the fully coupled solute-solvent interaction. An endpoints method thus utilizes the ensembles of intermolecular configurations of physical interest, which are often readily accessible with MD simulations given the present computational power. To illustrate properties of the formulation, we perform simulations of model systems consisting of a cavity particle in an aqueous solution containing a spherical hydrophobic solute of three different sizes from which ω(x) and the solute chemical potential can be calculated using endpoints DFT expressions. These are compared with corresponding results obtained using the approximations needed in order to evaluate the endpoints DFT charging integral when cavity particle simulation data is not available. We analyze a new approximation (two-points quadratic HNC) to the DFT charging integral which captures the correct behavior of the cavity distributions at both endpoints of the solute insertion. The behavior of the cavity particle in simple and complex liquids plays an important role in various theoretical treatments of the solute chemical potential. For pure Lennard-Jones fluids, the free energy to bring a cavity particle from the bulk to the center of a fluid particle is negative. However, for solutes of varying size, this is not generally true for Lennard-Jones fluids or the systems studied in this work. We carry out energetic and structural analyses of the cavity particle in aqueous solution with hydrophobic solutes of varying size and discuss the results in the context of the hydrophobic effect.

Multisolvent Models for Solvation Free Energy Predictions Using 3D-RISM Hydration Thermodynamic Descriptors

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.

Two dimensional fluid with one site-site associating point. Monte Carlo, integral equation and thermodynamic perturbation theory study

In this paper we propose a model for the two dimensional fluid with one site-site associating point. We studied its structural and thermodynamic properties by the Monte Carlo computer simulations, the site-site integral equation theory (RISM), the Wertheim's thermodynamic perturbation theory (TPT) and the Wertheim's integral equation theory (WIET) for associative liquids. The model can have arbitrary position of the associating point from the center of particles. All particles have Lennard-Jones core while interactions between associating points are modeled as Gaussian like potential where the interaction depends only on the distance between sites. The methods were used to study the thermodynamic and structural properties as a function of the position of associating point, temperature and density. The accuracy of the analytic theories were checked by comparing the theoretical results with the corresponding Monte Carlo ones. The theories are quite accurate for cases when the associating point is on the surface and only dimers can be formed. In this case, the theories correctly predict the pair correlation functions of the model, internal energy, ratios of free and bonded particles and chemical potential. This is no longer true when associating point is away from the surface of particles and the higher clusters are formed.

Free-energy analysis of physisorption on solid-liquid interface with the solution theory in the energy representation

Physisorption of urea on its crystal in contact with water was subject to energetics analysis with all-atom molecular dynamics simulation. The transfer free energy of urea to an adsorption site was treated in the framework of the energy-representation theory of solutions, which allows a fast computation of the free energy in an inhomogeneous environment with solid-liquid interface. The preference of adsorption was then compared between the (001) and (110) faces, and it was found that the physisorption is more favorable on (001) than on (110) in correspondence to the hydrogen bonding between the adsorbed urea and the crystal urea. Among the terrace configurations of adsorption, the attractive interaction governs the preferable site with a minor role of the repulsive interaction. The effect of an edge was also treated by examining the terrace and step and was shown to be strongly operative on the (110) face when the CO group of the adsorbed urea points toward the edge. The present work demonstrates that the solution theory can be a framework for analyzing the energetics of physisorption and addressing the roles of the crystal and liquid at the interface through the systematic decomposition of free energy.

Structure and permeability of ionomers studied by atomistic molecular simulation combined with the theory of solutions in the energy representation

Ionomers play a key role in forming the catalyst layer of polymer electrolyte fuel cells. In the present work, we performed atomistic molecular dynamics simulations and free-energy calculations with the energy-representation method for sulfonated polyethersulfone (SPES) and its derivatives toward the rational design of ionomers for carbon alloy catalysts. It was observed that H₂O aggregates strongly in the branched SPES systems with fluorocarbons and is located homogeneously in the systems without fluorocarbons. The O₂ permeability was then examined within the framework of the solubility-diffusion mechanism. The permeability was seen to be large for the branched SPES with fluorocarbons, indicating that the performance of ionomers as a permeation medium for O₂ may be tuned by the flexibility and branching of the polymer chain.

Probabilistic analysis for identifying the driving force of protein folding

Toward identifying the driving force of protein folding, energetics was analyzed in water for Trp-cage (20 residues), protein G (56 residues), and ubiquitin (76 residues) at their native (folded) and heat-denatured (unfolded) states. All-atom molecular dynamics simulation was conducted, and the hydration effect was quantified by the solvation free energy. The free-energy calculation was done by employing the solution theory in the energy representation, and it was seen that the sum of the protein intramolecular (structural) energy and the solvation free energy is more favorable for a folded structure than for an unfolded one generated by heat. Probabilistic arguments were then developed to determine which of the electrostatic, van der Waals, and excluded-volume components of the interactions in the protein-water system governs the relative stabilities between the folded and unfolded structures. It was found that the electrostatic interaction does not correspond to the preference order of the two structures. The van der Waals and excluded-volume components were shown, on the other hand, to provide the right order of preference at probabilities of almost unity, and it is argued that a useful modeling of protein folding is possible on the basis of the excluded-volume effect.

Computing conformational free energy differences in explicit solvent: An efficient thermodynamic cycle using an auxiliary potential and a free energy functional constructed from the end points

Many biomolecules undergo conformational changes associated with allostery or ligand binding. Observing these changes in computer simulations is difficult if their timescales are long. These calculations can be accelerated by observing the transition on an auxiliary free energy surface with a simpler Hamiltonian and connecting this free energy surface to the target free energy surface with free energy calculations. Here, we show that the free energy legs of the cycle can be replaced with energy representation (ER) density functional approximations. We compute: (1) The conformational free energy changes for alanine dipeptide transitioning from the right-handed free energy basin to the left-handed basin and (2) the free energy difference between the open and closed conformations of β-cyclodextrin, a "host" molecule that serves as a model for molecular recognition in host-guest binding. β-cyclodextrin contains 147 atoms compared to 22 atoms for alanine dipeptide, making β-cyclodextrin a large molecule for which to compute solvation free energies by free energy perturbation or integration methods and the largest system for which the ER method has been compared to exact free energy methods. The ER method replaced the 28 simulations to compute each coupling free energy with two endpoint simulations, reducing the computational time for the alanine dipeptide calculation by about 70% and for the β-cyclodextrin by > 95%. The method works even when the distribution of conformations on the auxiliary free energy surface differs substantially from that on the target free energy surface, although some degree of overlap between the two surfaces is required. © 2016 Wiley Periodicals, Inc.

Properties of the two-dimensional heterogeneous Lennard-Jones dimers: An integral equation study

Structural and thermodynamic properties of a planar heterogeneous soft dumbbell fluid are examined using Monte Carlo simulations and integral equation theory. Lennard-Jones particles of different sizes are the building blocks of the dimers. The site-site integral equation theory in two dimensions is used to calculate the site-site radial distribution functions and the thermodynamic properties. Obtained results are compared to Monte Carlo simulation data. The critical parameters for selected types of dimers were also estimated and the influence of the Lennard-Jones parameters was studied. We have also tested the correctness of the site-site integral equation theory using different closures.

Free-energy analysis of protein solvation with all-atom molecular dynamics simulation combined with a theory of solutions

The structure of a protein is strongly influenced by solvation. In the present review, the solvation effect is treated within the framework of statistical thermodynamics. The key quantity is the solvation free energy of protein and a fast method for its computation with explicit solvent is introduced. The applications of the method to the structural fluctuation of protein in water, structure determination of protein-protein complex, and urea effect on protein structure are then described, and the roles of solvent water and cosolvent are discussed from the standpoint of energetics.

Correlation analysis for heat denaturation of Trp-cage miniprotein with explicit solvent

Energetics was analyzed for Trp-cage miniprotein in water to elucidate the solvation effect in heat denaturation. The solvation free energy was computed for a set of protein structures at room and high temperatures with all-atom molecular dynamics simulation combined with the solution theory in the energy representation, and its correlations were investigated against the intramolecular (structural) energy of the protein and the average interaction energy of the protein with the solvent water. It was observed both at room and high temperatures that the solvation free energy is anticorrelated to the structural energy and varies in parallel to the electrostatic component of the protein-water interaction energy without correlations to the van der Waals and excluded-volume components. When the set of folded structures sampled at room temperature was compared with the set of unfolded ones at high temperature, it was found that the preference order of the two sets is in correspondence to the van der Waals and excluded-volume components in the sum of the protein intramolecular and protein-water intermolecular interactions and is not distinguished by the electrostatic component.

Fast and General Method To Predict the Physicochemical Properties of Druglike Molecules Using the Integral Equation Theory of Molecular Liquids

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.

Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections

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.

Properties of the Lennard-Jones dimeric fluid in two dimensions: an integral equation study

The thermodynamic and structural properties of the planar soft-sites dumbbell fluid are examined by Monte Carlo simulations and integral equation theory. The dimers are built of two Lennard-Jones segments. Site-site integral equation theory in two dimensions is used to calculate the site-site radial distribution functions for a range of elongations and densities and the results are compared with Monte Carlo simulations. The critical parameters for selected types of dimers were also estimated. We analyze the influence of the bond length on critical point as well as tested correctness of site-site integral equation theory with different closures. The integral equations can be used to predict the phase diagram of dimers whose molecular parameters are known.