Cosolvent preferential molecular interactions in aqueous solutions.

We extend the application of Kirkwood-Buff (KB) solution theory integrals to calculate cosolvent preferential interaction parameters from molecular simulations by deriving from a single simulation trajectory the excess chemical potential of the solute in addition to the solute-solvent molecular distribution functions. The solute excess chemical potential is derived from the potential distribution theorem (PDT) and used to define the local solvent domain around the solute, as distinguished from bulk solution. We show that this KB/PDT characterization of preferential molecular interactions resolves the problem of convergence of the preferential interaction parameter in the bulk solution limit, and as such, gives reliable estimates of preferential interaction parameters for methanol, ethanol, glycerol, and urea in aqueous cosolvent solutions with neopentane or tetramethyl ammonium ion as the solute. Preferential interaction parameters that are also calculated on the basis of cosolvent proximal distributions around the constituent methyl groups of the two solutes with the assumption of group additivity are in good agreement with those obtained by considering the molecular solute as a whole. The results suggest that this approach can be applied to estimate site-specific cosolvent preferential interaction parameters locally on the surface of complex, macromolecular solutes, such as proteins.

Non-van der Waals treatment of the hydrophobic solubilities of CF4.

A quasi-chemical theory implemented on the basis of molecular simulation is derived and tested for the hydrophobic hydration of CF4(aq). The theory formulated here subsumes a van der Waals treatment of solvation and identifies contributions to the hydration free energy of CF4(aq) that naturally arise from chemical contributions defined by quasi-chemical theory and fluctuation contributions analogous to Debye-Hückel or random phase approximations. The resulting Gaussian statistical thermodynamic model avoids consideration of hypothetical drying-then-rewetting problems and is physically reliable in these applications as judged by the size of the fluctuation contribution. The specific results here confirm that unfavorable tails of binding energy distributions reflect few-body close solute-solvent encounters. The solvent near-neighbors are pushed by the medium into unfavorable interactions with the solute, in contrast to the alternative view that a preformed interface is pulled by the solute-solvent attractive interactions into contact with the solute. The polyatomic model of CF4(aq) studied gives a satisfactory description of the experimental solubilities including the temperature dependence. The proximal distributions evaluated here for polyatomic solutes accurately reconstruct the observed distributions of water near these molecules which are nonspherical. These results suggest that drying is not an essential consideration for the hydrophobic solubilities of CF4, or of C(CH3)4 which is more soluble.

An analysis of molecular packing and chemical association in liquid water using quasichemical theory.

We calculate the hydration free energy of liquid TIP3P water at 298 K and 1 bar using a quasi-chemical theory framework in which interactions between a distinguished water molecule and the surrounding water molecules are partitioned into chemical associations with proximal (inner-shell) waters and classical electrostatic-dispersion interactions with the remaining (outer-shell) waters. The calculated free energy is found to be independent of this partitioning, as expected, and in excellent agreement with values derived from the literature. An analysis of the spatial distribution of inner-shell water molecules as a function of the inner-shell volume reveals that water molecules are preferentially excluded from the interior of large volumes as the occupancy number decreases. The driving force for water exclusion is formulated in terms of a free energy for rearranging inner-shell water molecules under the influence of the field exerted by outer-shell waters in order to accommodate one water molecule at the center. The results indicate a balance between chemical association and molecular packing in liquid water that becomes increasingly important as the inner-shell volume grows in size.

Effect of solute size and solute-water attractive interactions on hydration water structure around hydrophobic solutes.

Using Monte Carlo simulations, we investigated the influence of solute size and solute-water attractive interactions on hydration water structure around spherical clusters of 1, 13, 57, 135, and 305 hexagonally close-packed methanes and the single hard-sphere (HS) solute analogues of these clusters. We obtain quantitative results on the density of water molecules in contact with the HS solutes as a function of solute size for HS radii between 3.25 and 16.45 A. Analysis of these results based on scaled-particle theory yields a hydration free energy/surface area coefficient equal to 139 cal/(mol A2), independent of solute size, when this coefficient is defined with respect to the van der Waals surface of the solute. The same coefficient defined with respect to the solvent-accessible surface decreases with decreasing solute size for HS radii less than approximately 10 A. We also find that solute-water attractive interactions play an important role in the hydration of the methane clusters. Water densities in the first hydration shell of the three largest clusters are greater than bulk water density and are insensitive to the cluster size. In contrast, contact water densities for the HS analogues of these clusters decrease with solute size, falling below the bulk density of water for the two largest solutes. Thus, the large HS solutes dewet, while methane clusters of the same size do not.

Conformational equilibria of alkanes in aqueous solution: relationship to water structure near hydrophobic solutes.

Conformational free energies of butane, pentane, and hexane in water are calculated from molecular simulations with explicit waters and from a simple molecular theory in which the local hydration structure is estimated based on a proximity approximation. This proximity approximation uses only the two nearest carbon atoms on the alkane to predict the local water density at a given point in space. Conformational free energies of hydration are subsequently calculated using a free energy perturbation method. Quantitative agreement is found between the free energies obtained from simulations and theory. Moreover, free energy calculations using this proximity approximation are approximately four orders of magnitude faster than those based on explicit water simulations. Our results demonstrate the accuracy and utility of the proximity approximation for predicting water structure as the basis for a quantitative description of n-alkane conformational equilibria in water. In addition, the proximity approximation provides a molecular foundation for extending predictions of water structure and hydration thermodynamic properties of simple hydrophobic solutes to larger clusters or assemblies of hydrophobic solutes.

Hydration and conformational equilibria of simple hydrophobic and amphiphilic solutes.

We consider whether the continuum model of hydration optimized to reproduce vacuum-to-water transfer free energies simultaneously describes the hydration free energy contributions to conformational equilibria of the same solutes in water. To this end, transfer and conformational free energies of idealized hydrophobic and amphiphilic solutes in water are calculated from explicit water simulations and compared to continuum model predictions. As benchmark hydrophobic solutes, we examine the hydration of linear alkanes from methane through hexane. Amphiphilic solutes were created by adding a charge of +/-1e to a terminal methyl group of butane. We find that phenomenological continuum parameters fit to transfer free energies are significantly different from those fit to conformational free energies of our model solutes. This difference is attributed to continuum model parameters that depend on solute conformation in water, and leads to effective values for the free energy/surface area coefficient and Born radii that best describe conformational equilibrium. In light of these results, we believe that continuum models of hydration optimized to fit transfer free energies do not accurately capture the balance between hydrophobic and electrostatic contributions that determines the solute conformational state in aqueous solution.