Solvation potentials for flexible chain molecules in solution: on the validity of a pairwise decomposition

The effects of a solvent on the conformation of a flexible n-site solute molecule can be described formally in terms of an n-body solvation potential. Given the practical difficulty in computing such multibody potentials, it is common to carry out a pairwise decomposition in which the n-body potential is approximated by a sum of two-body potentials. Here we investigate the validity of this two-site approximation for short interaction-site chain-in-solvent systems. Using exact expressions for the conformation of an isolated chain, we construct a mapping between the full chain-in-solvent system and its solvation potential representation. We present results for both hard-sphere and square-well systems with n=5 that show that the two-site approximation is sufficient to completely capture the effects of an explicit solvent on chain conformation for a wide range of conditions (which include varying the solvent diameter in the hard-sphere system and varying the chain-solvent coupling in the square-well system). In all cases, a set of two-site potentials (one for each distinct site-site pair) is required. We also show that these two-site solvation potentials can be used to accurately compute a multisite intramolecular correlation function.     

The impact of surface area,volume, curvature,and Lennard–Jones potential to solvation modeling

This article explores the impact of surface area, volume, curvature, and Lennard–Jones (LJ) potential on solvation free energy predictions. Rigidity surfaces are utilized to generate robust analytical expressions for maximum, minimum, mean, and Gaussian curvatures of solvent–solute interfaces, and define a generalized Poisson–Boltzmann (GPB) equation with a smooth dielectric profile. Extensive correlation analysis is performed to examine the linear dependence of surface area, surface enclosed volume, maximum curvature, minimum curvature, mean curvature, and Gaussian curvature for solvation modeling. It is found that surface area and surfaces enclosed volumes are highly correlated to each other's, and poorly correlated to various curvatures for six test sets of molecules. Different curvatures are weakly correlated to each other for six test sets of molecules, but are strongly correlated to each other within each test set of molecules. Based on correlation analysis, we construct twenty six nontrivial nonpolar solvation models. Our numerical results reveal that the LJ potential plays a vital role in nonpolar solvation modeling, especially for molecules involving strong van der Waals interactions. It is found that curvatures are at least as important as surface area or surface enclosed volume in nonpolar solvation modeling. In conjugation with the GPB model, various curvature‐based nonpolar solvation models are shown to offer some of the best solvation free energy predictions for a wide range of test sets. For example, root mean square errors from a model constituting surface area, volume, mean curvature, and LJ potential are less than 0.42 kcal/mol for all test sets. © 2016 Wiley Periodicals, Inc.     

Electrolytes in a nanometer slab-confinement: ion-specific structure and solvation forces

We study the liquid structure and solvation forces of dense monovalent electrolytes (LiCl, NaCl, CsCl, and NaI) in a nanometer slab-confinement by explicit-water molecular dynamics (MD) simulations, implicit-water Monte Carlo (MC) simulations, and modified Poisson-Boltzmann (PB) theories. In order to consistently coarse-grain and to account for specific hydration effects in the implicit methods, realistic ion-ion and ion-surface pair potentials have been derived from infinite-dilution MD simulations. The electrolyte structure calculated from MC simulations is in good agreement with the corresponding MD simulations, thereby validating the coarse-graining approach. The agreement improves if a realistic, MD-derived dielectric constant is employed, which partially corrects for (water-mediated) many-body effects. Further analysis of the ionic structure and solvation pressure demonstrates that nonlocal extensions to PB (NPB) perform well for a wide parameter range when compared to MC simulations, whereas all local extensions mostly fail. A Barker-Henderson mapping of the ions onto a charged, asymmetric, and nonadditive binary hard-sphere mixture shows that the strength of structural correlations is strongly related to the magnitude and sign of the salt-specific nonadditivity. Furthermore, a grand canonical NPB analysis shows that the Donnan effect is dominated by steric correlations, whereas solvation forces and overcharging effects are mainly governed by ion-surface interactions. However, steric corrections to solvation forces are strongly repulsive for high concentrations and low surface charges, while overcharging can also be triggered by steric interactions in strongly correlated systems. Generally, we find that ion-surface and ion-ion correlations are strongly coupled and that coarse-grained methods should include both, the latter nonlocally and nonadditive (as given by our specific ionic diameters), when studying electrolytes in highly inhomogeneous situations.     

具有多体效应的胶体聚团的特征

考察了多体效应（用Ｓｔｕｔｔｏｎ－Ｃｈｅｎ势，ＳＣ）对胶体聚团的影响并与双体（ＬＪ）势下的结果作了比较。研究表明，ＳＣ和ＬＪ势下簇团的性质有其相似的一面，如：随着剪应力的增加，系统里颗粒的平均势能增加，而每个簇团的颗粒数减少；在较强的剪应力场里，簇团沿剪应力方向（Ｘ轴）被明显拉长且其主轴偏离Ｘ轴等。但它们间的差异也是明显的，在剪应力下ＳＣ系统内颗粒排列得更合理，从而使得平均位能比ＬＪ系统低约１－３     

Effective multipoles and Yukawa electrostatics in dressed molecule theory

In this paper we derive the multipolar expansion of the screened Coulomb potential in electrolyte solutions with molecular solvent. The solute and solvent molecules can have arbitrary sizes, shapes, and internal charge distributions. We use the exact statistical mechanical definition of renormalized charge distributions coming from "dressed molecule theory" to determine the effective multipoles of a molecule immersed in an electrolyte. The effects of many-body correlations are fully included in our formally exact theory. We restrict ourselves to sufficiently dilute solutions so the screened Coulomb potential decays for large distances like a Yukawa function, exp(-kappa r)/r, where r is the distance and 1/kappa is the decay length (it is normally different from the Debye length). The resulting "Yukawa electrostatics" differ in many respects from ordinary, unscreened electrostatics. The "Yukawa charge" of a molecule (the lowest order moment in the multipolar expansion) is in general not equal to its Coulombic charge and it is not the integral of the renormalized charge distribution of the molecule. Moreover, as shown in this paper, the multipolar expansion of the Yukawa potential does not correspond, contrary to the case of the Coulomb potential, to its asymptotic expansion for large r. As a consequence, the charge term in the multipolar expansion is not the leading term in the asymptotic expansion. Instead, for large r values, multipoles of all orders contribute to the leading asymptotic term. Thus, the electrostatic potential from, for example, an electroneutral solvent molecule in an electrolyte solution has generally the same range as that from an ion. The proper asymptotic expansion for electrostatic interactions in electrolytes is derived. It is briefly shown how the multipole expansion formalism can also be applied in the Poisson-Boltzmann approximation for primitive model electrolytes.     

Identifying low variance pathways for free energy calculations of molecular transformations in solution phase

Improving the efficiency of free energy calculations is important for many biological and materials design applications, such as protein-ligand binding affinities in drug design, partitioning between immiscible liquids, and determining molecular association in soft materials. We show that for any pair potential, moderately accurate estimation of the radial distribution function for a solute molecule is sufficient to accurately estimate the statistical variance of a sampling along a free energy pathway. This allows inexpensive analytical identification of low statistical error free energy pathways. We employ a variety of methods to estimate the radial distribution function (RDF) and find that the computationally cheap two-body "dilute gas" limit performs as well or better than 3D-RISM theory and other approximations for identifying low variance free energy pathways. With a RDF estimate in hand, we can search for pairwise interaction potentials that produce low variance. We give an example of a search minimizing statistical variance of solvation free energy over the entire parameter space of a generalized "soft core" potential. The free energy pathway arising from this optimization procedure has lower curvature in the variance and reduces the total variance by at least 50% compared to the traditional soft core solvation pathway. We also demonstrate that this optimized pathway allows free energies to be estimated with fewer intermediate states due to its low curvature. This free energy variance optimization technique is generalizable to solvation in any homogeneous fluid and for any type of pairwise potential and can be performed in minutes to hours, depending on the method used to estimate g(r).     

Solvation of hydrophobes in water and simple liquids

The solvation of nonpolar molecules in water and that in simple liquids are compared and contrasted. First, solvation thermodynamics is reviewed in a way that focuses on how the enthalpy and entropy of solvation depend on the choice of microscopic volume change v in the solvation process--including special choices v being zero (fixed-volume condition) and v being the partial molecular volume of a solute molecule (fixed-pressure condition)--and how the solvation quantities are related with temperature derivatives of the solvation free energy. Second, the solvation free energy and the solvation enthalpy of a Lennard-Jones (LJ) atom in model water are calculated in the parameter space representing the solute size and the strength of the solute-solvent interaction, and the results are compared with those for an LJ atom in the LJ solvent. The solvation diagrams showing domains of different types of solvation in the parameter space are obtained both for the constant-volume condition and for the constant-pressure condition. Similarities between water and the simple liquid are found when the constant-volume solvation is considered while a significant difference manifests itself in the fixed-pressure solvation. The domain of solvation of hydrophobic character in the parameter space is large in the constant-volume solvation both for water and for the simple liquid. When switched to the constant-pressure condition accompanying a microscopic volume change, the hydrophobic domain remains large in water but it becomes significantly small in the simple liquid. The contrasting results are due to the smallness of the thermal pressure coefficient of water at low temperatures.     

Finite-temperature full configuration interaction

The exact basis-set values of various thermodynamic potentials of a molecule are evaluated by the finite-temperature full configuration-interaction (FCI) method using ab initio molecular integrals over Gaussian-type orbitals. The thermodynamic potentials considered are the grand partition function, grand potential, internal energy, entropy, and chemical potential in the grand canonical ensemble as well as the partition function, Helmholtz energy, internal energy, and entropy in canonical ensemble. Approximations to FCI that are accurate at low and high temperatures are proposed, implemented, and tested. The results of finite-temperature FCI and its approximations are compared with one another as well as with the results of finite-temperature zeroth-order many-body perturbation theory, in which the Fermi–Dirac statistics is exact. Analytical asymptotic properties in the low- or high-temperature limits of some of these thermodynamic potentials are also given.     

Unique continuation for the magnetic Schrödinger equation

The unique-continuation property from sets of positive measure is here proven for the many-body magnetic Schrödinger equation. This property guarantees that if a solution of the Schrödinger equation vanishes on a set of positive measure, then it is identically zero. We explicitly consider potentials written as sums of either one-body or two-body functions, typical for Hamiltonians in many-body quantum mechanics. As a special case, we are able to treat atomic and molecular Hamiltonians. The unique-continuation property plays an important role in density-functional theories, which underpins its relevance in quantum chemistry.     

Theoretical determination of the standard reduction potentials of pheophytin-a in N,N-dimethyl formamide and membrane

Quantum mechanical/molecular mechanics (QM/MM) calculations were performed on the neutral, anionic, and dianionic forms of Pheophytin-a (Pheo-a) in N,N-dimethyl formamide (DMF) in order to calculate the absolute free energy of reduction of Pheo-a in solution. The geometry of the solvated species was optimized by restricted open-shell density functional treatment (ROB3LYP) using the 6-31G(d) basis set for the molecular species while the primary solvent shell consisting of 45 DMF molecules was treated by the MM method using the universal force field (UFF). Electronic energies of the neutral, anionic, and dianionic species were obtained by carrying out single point density functional theory (DFT) calculations using the 6-311+G(2d,2p) basis set on the respective ONIOM optimized geometries. The CHARMM27 force field was used to account for the dynamical nature of the primary solvation shell of 45 DMF molecules. In the calculations using solvent shells, the required atomic charges for each solvent molecule were obtained from ROB3LYP/6-31G(d) calculation on a single isolated DMF molecule. Randomly sampled configurations generated by the Monte Carlo (MC) technique were used to determine the contribution of the primary shell to the free energy of solvation of the three species. The dynamical nature of the primary shell significantly corrects the free energy of solvation. Frequency calculations at the ROB3LYP/6-31G(d) level were carried out on the optimized geometries of truncated 47-atom models of the neutral and ionic species in vacuum so as to determine the differences in thermal energy and molecular entropy. The Born energy of ion-dielectric interaction, the Onsager energy of dipole-dielectric interaction, and the Debye-Hückel energy of ion-ionic cloud interaction for the pheophytin-solvent aggregate were added as perturbative corrections. The Born interaction also makes a large contribution to the absolute free energy of reduction. An implicit solvation model (DPCM) was also employed for the calculation of standard reduction potentials in DMF. Both the models were successful in reproducing the standard reduction potentials. An explicit solvent treatment(QM/MM/MC + Born + Onsager + Debye corrections) yielded the one electron reduction potential of Pheo-a as -0.92 +/- 0.27 V and the two electron reduction potential as -1.34 +/- 0.25 V at 298.15 K, while the implicit solvent treatment yielded the corresponding values as -1.03 +/- 0.17 and -1.30 +/- 0.17 V, respectively. The calculated values more or less agree with the experimental midpoint potentials of -0.90 and -1.25 V, respectively. Moreover, a numerical finite difference Poisson-Boltzmann solver (FDPB) along with the DPCM methodology was employed to obtain the reduction potential of pheophytin in the thylakoid membrane. The calculated reduction potential value of -0.58 V is in excellent agreement with the reported value -0.61 V.     

Density dependence of the entropy and the solvation shell structure in supercritical water via molecular dynamics simulation

We perform molecular dynamics simulations of supercritical water (SCW) with a wide range of densities along a near critical isotherm using the simple point charge extended (SPC/E) pair potential in order to study the entropy and the solvation shell structure around a central water molecule. It is shown that both the translational and orientational two-particle correlation entropy terms can serve as the metrics of the translational-orientational structural orders in water and it is revealed that the translational structural order is very sensitive to the density variation in the gas-like and liquid-like region, while the orientational structural order is much more dependent upon compression in the medium-density SCW region. The comparison of the magnitudes of the full thermodynamic excess entropy and two-particle correlation entropy confirms the recent findings that the many-body terms other than two-body ones also present significant and non-neglectable contributions to the full excess entropy for the highly anomalous fluids like water. The analysis of entropy terms as a function of intermolecular distance and the orientational distribution functions as well as the three-dimensional spatial distribution functions indicate that the structural order occurs only in a much more diffused first solvation shell due to the elongated hydrogen bonds under supercritical conditions. It is revealed that no obvious second or higher neighbor shells occur in SCW, in contrast with the feature of normal liquid water that the anomalous decrease of translational order upon compression occurs mainly in the second shell.     

A novel algorithm for creating coarse-grained, density dependent implicit solvent models

Implicit solvent simulations are those in which solvent molecules are not explicitly simulated, and the solute-solute interaction potential is modified to compensate for the implicit solvent effect. Implicit solvation is well known in Brownian dynamics of dilute solutions but offers promise to speed up many other types of molecular simulations as well, including studies of proteins and colloids where the local density can vary considerably. This work examines implicit solvent potentials within a more general coarse-graining framework. While a pairwise potential between solute sites is relatively simple and ubiquitous, an additional parametrization based on the local solute concentration has the possibility to increase the accuracy of the simulations with only a marginal increase in computational cost. We describe here a method in which the radial distribution function and excess chemical potential of solute insertion for a system of Lennard-Jones particles are first measured in a fully explicit, all-particle simulation, and then reproduced across a range of solute particle densities in an implicit solvent simulation.     

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.     

Molecular density functional theory of solvation: from polar solvents to water

A classical density functional theory approach to solvation in molecular solvent is presented. The solvation properties of an arbitrary solute in a given solvent, both described by a molecular force field, can be obtained by minimization of a position and orientation-dependent free-energy density functional. In the homogeneous reference fluid approximation, limited to two-body correlations, the unknown excess term of the functional approximated by the angular-dependent direct correlation function of the pure solvent. We show that this function can be extracted from a preliminary MD simulation of the pure solvent by computing the angular-dependent pair distribution function and solving subsequently the molecular Ornstein-Zernike equation using a discrete angular representation. The corresponding functional can then be minimized in the presence of an arbitrary solute on a three-dimensional cubic grid for positions and Gauss-Legendre angular grid for orientations to provide the solvation structure and free-energy. This two-step procedure is proved to be much more efficient than direct molecular dynamics simulations combined to thermodynamic integration schemes. The approach is shown to be relevant and accurate for prototype polar solvents such as the Stockmayer solvent or acetonitrile. For water, although correct for neutral or moderately charged solute, it tends to underestimate the tetrahedral solvation structure around H-bonded solutes, such as spherical ions. This can be corrected by introducing suitable three-body correlation terms that restore both an accurate hydration structure and a satisfactory energetics.     

The solvation of Li and Na in acetonitrile from ab initio-derived many-body ion–solvent potentials

Several Li⁺- and Na⁺-acetonitrile models were derived from ab initio calculations at the counterpoise-corrected MP2/TZV++(d,p) level for distorted ion-(MeCN)_n clusters with n=1, 4 and 6. Two different many-body ion-acetonitrile models were constructed: an effective three-body potential for use with the six-site effective pair model of Böhm et al., and an effective polarizable many-body model. The polarizable acetonitrile model used in the latter model is a new empirical model which was also derived in the present paper. Mainly for comparative purposes, two ion-acetonitrile pair potentials were also constructed from the ab initio cluster calculations: one pure pair potential and one effective pair potential. Using all these potential models, MD simulations in the NPT ensemble were performed for the pure acetonitrile liquid and for Li⁺(MeCN) and Na⁺(MeCN) solutions with 1 ion in 512 solvent molecules and with a simulation time of at least 120 ps per system. Thermodynamic properties, solvation-shell structure and the self-diffusion coefficient of the ions and of the solvent molecules were calculated and compared between the different models and with experimental data, where available. The Li⁺ ion is found to be four-coordinated when the new many-body potentials are used, in contrast to the six-coordinated structure obtained for the pure pair and effective pair potentials. The coordination number of Na⁺ is close to six for all the models derived here, although the coordination number becomes slightly smaller with the many-body potentials. For both ions, the solvent molecules in the first shell point their nitrogen ends towards the cation, while in the second shell the opposite orientation is the most common.     

Monte Carlo simulation of water solvent with biomolecules: Serine and the corresponding zwitterion

Monte Carlo simulation results are reported for clusters consisting of 250 water molecules surrounding serine, both in the neutral form and in two zwitterionic conformers (in order to gain some insight into conformational effects). Calculations were carried out at 300 K and using two-body potentials obtained by means of quantum-mechanical calculations. The spatial dependence of the average interaction energies was investigated. The solvation structure was investigated by means of radial distribution functions and probability density maps, which showed a few water molecules directly solvating the hydrophilic groups and, beyond them, a more or less rich and complex hydrogen-bonded network of solvent molecules.     

van der Waals forces between nanoclusters: importance of many-body effects

van der Waals interactions between nanoclusters have been calculated with a self-consistent, coupled dipole method. The method accounts for all many-body (MB) effects. Comparison is made between the exact potential energy, V, and the values obtained with two alternative methods: the sum of two-body interactions and the sum of two-body and three-body interactions. For all cases considered, the three-body term alone does not accurately represent the MB contributions to V. MB contributions are especially large for shape-anisotropic clusters.     

MLIMC: Machine Learning-Based Implicit-Solvent Monte Carlo

Monte Carlo (MC) methods are important computational tools for molecular structure optimizations and predictions. When solvent effects are explicitly considered, MC methods become very expensive due to the large degree of freedom associated with the water molecules and mobile ions. Alternatively implicit-solvent MC can largely reduce the computational cost by applying a mean field approximation to solvent effects and meanwhile maintains the atomic detail of the target molecule. The two most popular implicit-solvent models are the Poisson-Boltzmann (PB) model and the Generalized Born (GB) model in a way such that the GB model is an approximation to the PB model but is much faster in simulation time. In this work, we develop a machine learning-based implicit-solvent Monte Carlo (MLIMC) method by combining the advantages of both implicit solvent models in accuracy and efficiency. Specifically, the MLIMC method uses a fast and accurate PB-based machine learning (PBML) scheme to compute the electrostatic solvation free energy at each step. We validate our MLIMC method by using a benzene-water system and a protein-water system. We show that the proposed MLIMC method has great advantages in speed and accuracy for molecular structure optimization and prediction.     

Transient cavities and the excess chemical potentials of hard-spheroid solutes in dipolar hard-sphere solvents

Monte Carlo computer simulations are used to study transient cavities and the solvation of hard-spheroid solutes in dipolar hard-sphere solvents. The probability distribution of spheroidal cavities in the solvent is shown to be well described by a Gaussian function, and the variations of fit parameters with cavity elongation and solvent properties are analyzed. The excess chemical potentials of hard-spheroid solutes with aspect ratios x in the range of 15< or =x< or =5, and with volumes between 1 and 20 times that of a solvent molecule, are presented. It is shown that for a given molecular volume and solvent dipole moment (or temperature) a spherical solute has the lowest excess chemical potential and hence the highest solubility, while a prolate solute with aspect ratio x should be more soluble than an oblate solute with aspect ratio 1x. For a given solute molecule, the excess chemical potential increases with increasing temperature; this same trend can be observed in hydrophobic solvation. A scaled-particle theory based on the solvent equation of state and a fitted solute-solvent interfacial tension shows excellent agreement with the simulation results over the whole range of solute elongations and volumes considered. An information-theoretic model based on the solvent density and radial distribution function is less successful, being accurate only for small solute volumes and low solvent densities.