Continuous medium theory for nonequilibrium solvation: IV. Solvent reorganization energy of electron transfer based on conductor-like screening model

In this work the authors present some evidences of defects in the popular continuous medium theories for nonequilibrium solvation. Particular attention has been paid to the incorrect reversible work approach. After convincing reasoning, the nonequilibrium free energy has been formulated to an expression different from the traditional ones. In a series of recent works by the authors, new formulations and some analytical application models for ultrafast processes were developed. Here, the authors extend the new theory to the cases of discrete bound charge distributions and present the correct form of the nonequilibrium solvation energy in such cases. A numerical solution method is applied to the evaluation of solvent reorganization energy of electron transfer. The test calculation for biphenyl-cyclohexane-naphthalene anion system achieves excellent agreement with the experimental fitting. The central importance presented in this work is the very simple and a consistent form of nonequilibrium free energy for both continuous and discrete charge distributions, based on which the new models can be established.     

Thermodynamics for nonequilibrium solvation and numerical evaluation of solvent reorganization energy

This work presents a thermodynamic method for treating nonequilibrium solvation. By imposing an extra electric field onto the nonequilibrium solvation system, a virtual constrained equilibrium state is prepared. In this way, the free energy difference between the real nonequilibrium state and the con-strained equilibrium one is simply the potential energy of the nonequilibrium polarization in the extra electronic field, according to thermodynamics. Further, new expressions of nonequilibrium solvation energy and solvent reorganization energy have been formulated. Analysis shows that the present formulations will give a value of reorganization energy about one half of the traditional Marcus theory in polar solvents, thus the explanation on why the traditional theory tends to overestimate this quantity has been found out. For the purpose of numerical determination of solvent reorganization energy, we have modified Gamess program on the basis of dielectric polarizable continuum model. Applying the procedure to the well-investigated intramolecular electron transfer in biphenyl-androstane-naphthyl and biphenyl-androstane-phenanthryl systems, the numerical results of solvent reorganization energy have been found to be in good agreement with the experimental fittings.     

Ermod: Fast and versatile computation software for solvation free energy with approximate theory of solutions

ERmod is a software package to efficiently and approximately compute the solvation free energy using the method of energy representation. Molecular simulation is to be conducted at two condensed‐phase systems of the solution of interest and the reference solvent with test‐particle insertion of the solute. The subprogram ermod in ERmod then provides a set of energy distribution functions from the simulation trajectories, and another subprogram slvfe determines the solvation free energy from the distribution functions through an approximate functional. This article describes the design and implementation of ERmod, and illustrates its performance in solvent water for two organic solutes and two protein solutes. Actually, the free‐energy computation with ERmod is not restricted to the solvation in homogeneous medium such as fluid and polymer and can treat the binding into weakly ordered system with nano‐inhomogeneity such as micelle and lipid membrane. ERmod is available on web at http://sourceforge.net/projects/ermod . © 2014 Wiley Periodicals, Inc.     

Breaking the polar‐nonpolar division in solvation free energy prediction

Implicit solvent models divide solvation free energies into polar and nonpolar additive contributions, whereas polar and nonpolar interactions are inseparable and nonadditive. We present a feature functional theory (FFT) framework to break this ad hoc division. The essential ideas of FFT are as follows: (i) representability assumption: there exists a microscopic feature vector that can uniquely characterize and distinguish one molecule from another; (ii) feature‐function relationship assumption: the macroscopic features, including solvation free energy, of a molecule is a functional of microscopic feature vectors; and (iii) similarity assumption: molecules with similar microscopic features have similar macroscopic properties, such as solvation free energies. Based on these assumptions, solvation free energy prediction is carried out in the following protocol. First, we construct a molecular microscopic feature vector that is efficient in characterizing the solvation process using quantum mechanics and Poisson–Boltzmann theory. Microscopic feature vectors are combined with macroscopic features, that is, physical observable, to form extended feature vectors. Additionally, we partition a solvation dataset into queries according to molecular compositions. Moreover, for each target molecule, we adopt a machine learning algorithm for its nearest neighbor search, based on the selected microscopic feature vectors. Finally, from the extended feature vectors of obtained nearest neighbors, we construct a functional of solvation free energy, which is employed to predict the solvation free energy of the target molecule. The proposed FFT model has been extensively validated via a large dataset of 668 molecules. The leave‐one‐out test gives an optimal root‐mean‐square error (RMSE) of 1.05 kcal/mol. FFT predictions of SAMPL0, SAMPL1, SAMPL2, SAMPL3, and SAMPL4 challenge sets deliver the RMSEs of 0.61, 1.86, 1.64, 0.86, and 1.14 kcal/mol, respectively. Using a test set of 94 molecules and its associated training set, the present approach was carefully compared with a classic solvation model based on weighted solvent accessible surface area. © 2017 Wiley Periodicals, Inc.     

Continuous medium theory for nonequilibrium solvation: II. Interaction energy between solute charge and reaction field and single-sphere model for spectral shift

On the basis of continuous medium theory, a model for evaluation of spectral shifts in solution has been developed in this work. The interaction energy between solute dipole and reaction field and the self-energy of the reaction field have been formulated through derivations. Applying the interaction energy expression together with the point dipole approximation to the case of spherical cavity produces new formulations of spectral shifts. The same expression of electrostatic free energy of the nonequilibrium state is achieved by integrating the change of the electrostatic free energy for a charging process. Moreover, generalized formulations evaluating spectral shifts have been established in the charge-potential notation, and the reduction of them to the point dipole case consistently leads to the same formulations of spectral shifts as those by interaction energy approach. Mutual supports provide convincing evidences for the reliability of the present results. In this work, attentions are particularly paid to the conclusion of zero self-energy of the reaction field, which is different from the previous theory. Reasoning and arguments are given on this point. From the present derivations, it is concluded that the spectral shifts of light absorption and emission were theoretically exaggerated in the past, in particular, by a factor of 2 for the spectral shift sum.     

Application of the frozen atom approximation to the GB/SA continuum model for solvation free energy

The generalized Born/surface area (GB/SA) continuum model for solvation free energy is a fast and accurate alternative to using discrete water molecules in molecular simulations of solvated systems. However, computational studies of large solvated molecular systems such as enzyme-ligand complexes can still be computationally expensive even with continuum solvation methods simply because of the large number of atoms in the solute molecules. Because in such systems often only a relatively small portion of the system such as the ligand binding site is under study, it becomes less attractive to calculate energies and derivatives for all atoms in the system. To curtail computation while still maintaining high energetic accuracy, atoms distant from the site of interest are often frozen; that is, their coordinates are made invariant. Such frozen atoms do not require energetic and derivative updates during the course of a simulation. Herein we describe methodology and results for applying the frozen atom approach to both the generalized Born (GB) and the solvent accessible surface area (SASA) parts of the GB/SA continuum model for solvation free energy. For strictly pairwise energetic terms, such as the Coulombic and van-der-Waals energies, contributions from pairs of frozen atoms can be ignored. This leaves energetic differences unaffected for conformations that vary only in the positions of nonfrozen atoms. Due to the nonlocal nature of the GB analytical form, however, excluding such pairs from a GB calculation leads to unacceptable inaccuracies. To apply a frozen-atom scheme to GB calculations, a buffer region within the frozen-atom zone is generated based on a user-definable cutoff distance from the nonfrozen atoms. Certain pairwise interactions between frozen atoms in the buffer region are retained in the GB computation. This allows high accuracy in conformational GB comparisons to be maintained while achieving significant savings in computational time compared to the full (nonfrozen) calculation. A similar approach for using a buffer region of frozen atoms is taken for the SASA calculation. The SASA calculation is local in nature, and thus exact SASA energies are maintained. With a buffer region of 8 A for the frozen-atom cases, excellent agreement in differences in energies for three different conformations of cytochrome P450 with a bound camphor ligand are obtained with respect to the nonfrozen cases. For various minimization protocols, simulations run 2 to 10.5 times faster and memory usage is reduced by a factor of 1.5 to 5. Application of the frozen atom method for GB/SA calculations thus can render computationally tractable biologically and medically important simulations such as those used to study ligand-receptor binding conformations and energies in a solvated environment.     

Continuous medium theory for nonequilibrium solvation: III. Solvation shift by monopole approximation and multipole expansion in spherical cavity

According to the classical electrodynamics, a new and reasonable method about electrostatic energy decomposition of the solute-solvent system has been proposed in this work by introducing the concept of spring energy. This decomposition in equilibrium solvation gives the clear comprehension for different parts of total electrostatic free energy. Logically extending this cognition to nonequilibrium leads to the new formula of electrostatic free energy of nonequilibrium state. Furthermore, the general solvation shift for light absorption/emission has been reformulated and applied to the ideal sphere case with the monopole approximation and multipole expansion. Solvation shifts in vertical ionizations of atomic ions of some series of main group elements have been investigated with monopole approximation, and the variation tendency of the solvation shift versus atomic number has been discussed. Moreover, the solvation shift in photoionization of nitrate anion in glycol has been investigated by the multipole expansion method.     

Calculation of solvation free energy using RISM theory for peptide in salt solution

We developed a robust, highly efficient algorithm for solving the full reference interaction site model (RISM) equations for salt solutions near a solute molecule with many atomic sites. It was obtained as an extension of our previously reported algorithm for pure water near the solute molecule. The algorithm is a judicious hybrid of the Newton–Raphson and Picard methods. The most striking advantage is that the Jacobian matrix is just part of the input data and need not be recalculated at all. To illustrate the algorithm, we solved the full RISM equations for a dipeptide (NH₂(SINGLE BOND)CHCH₃(SINGLE BOND)CONH(SINGLE BOND)CHCH₃(SINGLE BOND)COOH) in a 1 M NaCl solution. The extended simple point charge (SPC/E) model was employed for water molecules. Two different conformations of the dipeptide were considered. It was assumed for each conformation that the dipeptide was present either as an un-ionized form or as a zwitterion. The structure of the salt solution near the dipeptide and salt effects on the solvation free energy were also discussed. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1724–1735, 1998     

New approach to free energy of solvation applying continuum models to molecular dynamics simulation

A new approach to the calculation of the free energy of solvation from trajectories obtained by molecular dynamics simulation is presented. The free energy of solvation is computed as the sum of three contributions originated at the cavitation of the solute by the solvent, the solute-solvent nonpolar (repulsion and dispersion) interactions, and the electrostatic solvation of the solute. The electrostatic term is calculated based on ideas developed for the broadly used continuum models, the cavitational contribution from the excluded volume by the Claverie-Pierotti model, and the Van der Waals term directly from the molecular dynamics simulation. The proposed model is tested for diluted aqueous solutions of simple molecules containing a variety of chemically important functions: methanol, methylamine, water, methanethiol, and dichloromethane. These solutions were treated by molecular dynamics simulations using SPC/E water and the OPLS force field for the organic molecules. Obtained free energies of solvation are in very good agreement with experimental data.     

Ab initio quantum mechanics-based free energy perturbation method for calculating relative solvation free energies

A free energy perturbation (FEP) method was developed that uses ab initio quantum mechanics (QM) for treating the solute molecules and molecular mechanics (MM) for treating the surroundings. Like our earlier results using AM1 semi empirical QMs, the ab initio QM/MM-based FEP method was shown to accurately calculate relative solvation free energies for a diverse set of small molecules that differ significantly in structure, aromaticity, hydrogen bonding potential, and electron density. Accuracy was similar to or better than conventional FEP methods. The QM/MM-based methods eliminate the need for time-consuming development of MM force field parameters, which are frequently required for drug-like molecules containing structural motifs not adequately described by MM. Future automation of the method and parallelization of the code for Linux 128/256/512 clusters is expected to enhance the speed and increase its use for drug design and lead optimization.     

A new set of atomic radii for accurate estimation of solvation free energy by Poisson–Boltzmann solvent model

The Poisson–Boltzmann implicit solvent (PB) is widely used to estimate the solvation free energies of biomolecules in molecular simulations. An optimized set of atomic radii (PB radii) is an important parameter for PB calculations, which determines the distribution of dielectric constants around the solute. We here present new PB radii for the AMBER protein force field to accurately reproduce the solvation free energies obtained from explicit solvent simulations. The presented PB radii were optimized using results from explicit solvent simulations of the large systems. In addition, we discriminated PB radii for N‐ and C‐terminal residues from those for nonterminal residues. The performances using our PB radii showed high accuracy for the estimation of solvation free energies at the level of the molecular fragment. The obtained PB radii are effective for the detailed analysis of the solvation effects of biomolecules. © 2014 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Accurate estimation of solvation free energy using polynomial fitting techniques

This report details an approach to improve the accuracy of free energy difference estimates using thermodynamic integration data (slope of the free energy with respect to the switching variable λ) and its application to calculating solvation free energy. The central idea is to utilize polynomial fitting schemes to approximate the thermodynamic integration data to improve the accuracy of the free energy difference estimates. Previously, we introduced the use of polynomial regression technique to fit thermodynamic integration data (Shyu and Ytreberg, J Comput Chem, 2009, 30, 2297). In this report we introduce polynomial and spline interpolation techniques. Two systems with analytically solvable relative free energies are used to test the accuracy of the interpolation approach. We also use both interpolation and regression methods to determine a small molecule solvation free energy. Our simulations show that, using such polynomial techniques and nonequidistant λ values, the solvation free energy can be estimated with high accuracy without using soft‐core scaling and separate simulations for Lennard‐Jones and partial charges. The results from our study suggest that these polynomial techniques, especially with use of nonequidistant λ values, improve the accuracy for ΔF estimates without demanding additional simulations. We also provide general guidelines for use of polynomial fitting to estimate free energy. To allow researchers to immediately utilize these methods, free software and documentation is provided via http://www.phys.uidaho.edu/ytreberg/software . © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Implicit solvation in the self-consistent mean field theory method: sidechain modelling and prediction of folding free energies of protein mutants

The Atomic Solvation Parameter (ASP) model is one of the simplest models of solvation, in which the solvation free energy of a molecule is proportional to the solvent accessible surface area (SAS) of its atoms. However, until now this model had not been incorporated into the Self-Consistent Mean Field Theory (SCMFT) method for modelling sidechain conformations in proteins. The reason for this is that SAS is a many-body quantity and, thus, it is not obvious how to define it within the Mean Field (MF) framework, where multiple copies of each sidechain exist simultaneously. Here, we present a method for incorporating an SAS-based potential, such as the ASP model, into SCMFT. The theory on which the method is based is exact within the MF framework, that is, it does not depend on a pairwise or any other approximation of SAS. Therefore, SAS can be calculated to arbitrary accuracy. The method is computationally very efficient: only 7.6% slower on average than the method without solvation. We applied the method to the prediction of sidechain conformation, using as a test set high-quality solution structures of 11 proteins. Solvation was found to substantially improve the prediction accuracy of well-defined surface sidechains. We also investigated whether the methodology can be applied to prediction of folding free energies of protein mutants, using a set of barnase mutants. For apolar mutants, the modest correlation observed between calculated and observed folding free energies without solvation improved substantially when solvation was included, allowing the prediction of trends in the folding free energies of this type of mutants. For polar mutants, correlation was not significant even with solvation. Several other factors also responsible for the correlation were identified and analysed. From this analysis, future directions for applying and improving the present methodology are discussed.     

Localization and quantification of hydrophobicity: The molecular free energy density (MolFESD) concept and its application to sweetness recognition

A method for the localization, the quantification, and the analysis of hydrophobicity of a molecule or a molecular fragment is presented. It is shown that the free energy of solvation for a molecule or the transfer free energy from one solvent to another can be represented by a surface integral of a scalar quantity, the molecular free energy surface density (MolFESD), over the solvent accessible surface of that molecule. This MolFESD concept is based on a model approach where the solvent molecules are considered to be small in comparison to the solute molecule, and the solvent can be represented by a continuous medium with a given dielectric constant. The transfer energy surface density for a 1-octanol/water system is empirically determined employing a set of atomic increment contributions and distance dependent membership functions measuring the contribution of the increments to the surface value of the MolFESD. The MolFESD concept can be well used for the quantification of the purely hydrophobic contribution to the binding constants of molecule-receptor complexes. This is demonstrated with the sweeteners sucrose and sucralose and various halogen derivatives. Therein the relative sweetness, which is assumed to be proportional to the binding constant, nicely correlates to the surface integral over the positive, hydrophobic part of the MolFESD, indicating that the sweetness receptor can be characterized by a highly flexible hydrophobic pocket instead of a localized binding site.     

Free energies of solvation in the context of protein folding: Implications for implicit and explicit solvent models

Implicit solvent models for biomolecular simulations have been developed to use in place of more expensive explicit models; however, these models make many assumptions and approximations that are likely to affect accuracy. Here, the changes in free energies of solvation upon folding of several fast folding proteins are calculated from previously run μs–ms simulations with a number of implicit solvent models and compared to the values needed to be consistent with the explicit solvent model used in the simulations. In the majority of cases, there is a significant and substantial difference between the values calculated from the two approaches that is robust to the details of the calculations. These differences could only be remedied by selecting values for the model parameters—the internal dielectric constant for the polar term and the surface tension coefficient for the nonpolar term—that were system‐specific or physically unrealistic. We discuss the potential implications of our findings for both implicit and explicit solvent simulations. © 2015 Wiley Periodicals, Inc.     

Correlation of the gibbs free energy of transfer of fluoride ion from water to other solvents with solvent actidity and polarity parameters

Linear free energy relationships between the free energy of transfer of fluoride ion from water to other solvents or solvent mixtures and each of several parameters describing the acidity of the solvent were found. The correlation is about of the same quality for the A parameter of Swain or the acceptor number, AN, and almost as good for the parameter of Taft. The correlation was somewhat poorer with the E _T (30) parameter.     

Application of reference‐modified density functional theory: Temperature and pressure dependences of solvation free energy

Recently, we proposed a reference‐modified density functional theory (RMDFT) to calculate solvation free energy (SFE), in which a hard‐sphere fluid was introduced as the reference system instead of an ideal molecular gas. Through the RMDFT, using an optimal diameter for the hard‐sphere reference system, the values of the SFE calculated at room temperature and normal pressure were in good agreement with those for more than 500 small organic molecules in water as determined by experiments. In this study, we present an application of the RMDFT for calculating the temperature and pressure dependences of the SFE for solute molecules in water. We demonstrate that the RMDFT has high predictive ability for the temperature and pressure dependences of the SFE for small solute molecules in water when the optimal reference hard‐sphere diameter determined for each thermodynamic condition is used. We also apply the RMDFT to investigate the temperature and pressure dependences of the thermodynamic stability of an artificial small protein, chignolin, and discuss the mechanism of high‐temperature and high‐pressure unfolding of the protein. © 2017 Wiley Periodicals, Inc.     

Electrostatic contribution to the ion solvation energy: cavity effects

Ion solvation process has been analysed for the spherically symmetrical system where an ion is located inside a cavity surrounded by an isotropic nonlocal dielectric medium. It has been proven that for any dielectric properties of the medium, the electric field outside the cavity as well as the ion solvation energy depend only on the total ion charge but not of the particular distribution of the ion charge density inside the cavity. These characteristics remain unchanged if the charge is displaced from the external boundary of the cavity into it. Analytical formulas for them have been derived for a particular model of the nonlocal dielectric function. Comparison of results for the solvation energy on the basis of this new theory and of the conventional approach (disregarding the existence of the cavity) shows a significant difference between their predictions if the ion charge is displaced inside the ion cavity.