Molecular solvation in water-methanol and water-sorbitol mixtures: the roles of preferential hydration, hydrophobicity, and the equation of state

Molecular dynamics simulations of aqueous mixtures of methanol and sorbitol were performed over a wide range of binary composition, density (pressure), and temperature to study the equation of state and solvation of small apolar solutes. Experimentally, methanol is a canonical solubilizing agent for apolar solutes and a protein denaturant in mixed-aqueous solvents; sorbitol represents a canonical "salting-out" or protein-stabilizing cosolvent. The results reported here show increasing sorbitol concentration under isothermal, isobaric conditions results in monotonic increases in apolar solute excess chemical potential (mu2ex) over the range of experimentally relevant temperatures. For methanol at elevated temperatures, increasing cosolvent composition results in monotonically decreasing mu2ex. However, at lower temperatures mu2ex exhibits a maximum versus cosolvent concentration, as seen experimentally for Ar in ethanol-water solutions. Both density anomalies and hydrophobic effects--characterized by temperatures of density maxima and apolar solute solubility minima, respectively--are suppressed upon addition of either sorbitol or methanol at all temperatures and compositions simulated here. Thus, the contrasting effects of sorbitol and methanol on solute chemical potential cannot be explained by qualitative differences in their ability to enhance or suppress hydrophobic effects. Rather, we find mu2ex values across a broad range of temperatures and cosolvent composition can be quantitatively explained in terms of isobaric changes in solvent density--i.e., the equation of state--along with the corresponding packing fraction of the solvent. Analysis in terms of truncated preferential interaction parameters highlights that care must be taken in interpreting cosolvent effects on solvation in terms of local preferential hydration.     

A molecular dynamics study of the dielectric properties of aqueous solutions of alanine and alanine dipeptide

Molecular dynamics simulations were used to compute the frequency-dependent dielectric susceptibility of aqueous solutions of alanine and alanine dipeptide. We studied four alanine solutions, ranging in concentration from 0.13-0.55 mol/liter, and two solutions of alanine dipeptide (0.13 and 0.27 mol/liter). In accord with experiment we find a strong dielectric increment for both solutes, whose molecular origin is shown to be the zwitterionic nature of the solutes. The dynamic properties were analyzed based on a dielectric component analysis into solute, a first hydration shell, and all remaining (bulk) waters. The results of this three component decomposition were interpreted directly, as well as by uniting the solute and hydration shell component to a "suprasolute" component. In both approaches three contributions to the frequency-dependent dielectric properties can be discerned. The quantitatively largest and fastest component arises from bulk water [i.e., water not influenced by the solute(s)]. The interaction between waters surrounding the solute(s) (the hydration shell) and bulk water molecules leads to a relaxation process occurring on an intermediate time scale. The slowest relaxation process originates from the solute(s) and the interaction of the solute(s) with the first hydration shell and bulk water. The primary importance of the hydration shell is the exchange of shell and bulk waters; the self-contribution from bound water molecules is comparatively small. While in the alanine solutions the solute-water cross-terms are more important than the solute self-term, the solute contribution is larger in the dipeptide solutions. In the latter systems a much clearer separation of time scales between water and alanine dipeptide related properties is observed. The similarities and differences of the dielectric properties of the amino acid/peptide solutions studied in this work and of solutions of mono- and disaccharides and of the protein ubiquitin are discussed.     

Water structure about the dimer and hexamer repeat units of amylose from molecular dynamics computer simulations

We have analyzed a set of molecular dynamics (MD) trajectories of maltose in vacuum and water for solute imposed structuring on the solvent. To do this, we used a novel technique to calculate water probability densities to locate the areas in which the solvent is most populated in the maltose solution. We found that only the layer of water within the first maltose hydration shell has a probability density 50% and greater than that of bulk water. On investigating this water layer using Voronoi polyhedra (VP) analysis it was seen that only the waters adjacent to the hydrophobic (CH and CH₂) groups are more structured than bulk water. We found that in a maltose solution of approximately 1.0 g/cm³ the solute does not disrupt the structure of the surrounding water beyond the first hydration shell. Next we performed a 700‐ps MD simulation of a maltohexaose strand in a box of 4096 SPC/E waters. The water probability density calculations and the VP analysis of the maltohexaose solution show that the larger amylose repeat unit decreases the solvent configurational entropy of the water beyond the first hydration shell. Analysis of this trajectory reveals that the helical conformation of the maltohexaose strand is preserved via bridging intermolecular water hydrogen bonds, indicating that a single amylose helical turn in water is preserved by hydrophilic and not hydrophobic interactions. Using VP analysis we present a method to accurately determine the number of water molecules in the first hydration shell of dissolved solutes. In the case of maltose, there are 40 water molecules in this shell, while for maltohexaose the number is 98. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 445–456, 2001     

Hydrophobic hydration in an orientational lattice model

To shed light on the microscopic mechanism of hydrophobic hydration, we study a simplified lattice model for water solutions in which the orientational nature of hydrogen bonding as well as the degeneracy related to proton distribution are taken into account. Miscibility properties of the model are looked at for both polar (hydrogen bonding) and nonpolar (non-hydrogen bonding) solutes. A quasichemical solution for the pure system is reviewed and extended to include the different kinds of solute. A Monte Carlo study of our model yields a novel feature for the local structure of the hydration layer: energy correlation relaxation times for solvation water are larger than the corresponding relaxation times for bulk water. This result suggests the presence of ordering of water particles in the first hydration shell. A nonassociating model solvent, represented by a lattice gas, presents opposite behavior, indicating that this effect is a result of the directionality of the interaction. In presence of polar solutes, we find an ordered mixed pseudophase at low temperatures, indicating the possibility of closed loops of immiscibility.     

纳米受限下溶质水化结构的分子模拟

利用分子动力学模拟研究了五种不同种类的溶质分子(K+, Mg2+, Cl-, K-和K0)在直径为0.60-1.28 nm的纳米碳管内的水化结构. 模拟结果揭示了单电荷溶质、双电荷溶质和中性溶质在受限条件下具有不同的水化行为. 单价溶质的配位数只有在直径不大于0.73 nm的纳米碳管内才会明显减少. 和带有电荷的溶质不同, 中性溶质的配位数对纳米碳管直径的改变非常敏感, 并且随着管径的减小而迅速减少. 模拟结果还表明带单价正电荷的溶质(K+)第一配位层水分子的取向结构会随着纳米碳管直径的改变发生变化, 而其他溶质配位层取向结构在本文所涉及的纳米碳管内都几乎和体相中一致. 在直径大于1.0 nm的纳米碳管中, K+的配位层取向结构有序度随着管径的减小而单调下降, 但是在直径小于1.0 nm的纳米碳管中, 随着碳管管径的减小而迅速上升. 在两个最窄的纳米碳管内, 其结构有度甚至高于体相. 双电荷溶质的水化结构在本文所研究的碳管直径范围内和体相完全一致, 即使在直径只有0.6 nm的碳管内也无任何改变.     

Orientational dynamics for an amphiphilic-solvent solution

In this work, we performed Monte Carlo simulations on a lattice model for spontaneous amphiphilic aggregation, in order to study the orientational and hydrogen-bonding dynamics of water on different regions inside the micellar solution. We employed an associating lattice gas model that mimics the aqueous solvent, which presents a rich phase diagram with first- and second-order transition lines. Even though this is a simplified model, it makes possible to investigate the orientational dynamics of water in an equilibrium solution of amphiphiles, as well as the influence of the different phases of the solvent in the interfacial and bulk water dynamics. By means of extensive simulations, we showed that, at high temperatures, the behavior of the orientational relaxation and hydrogen bonding of water molecules in the bulk, first, and second hydration shells are considerable different. We observe the appearance of a very slow component for water molecules in the first hydration shell of micelles when the system reaches a high-density phase, consistent with previous theoretical and experimental studies concerning biological water. Also, at high temperatures, we find that water molecules in the second hydration shell of micelles have an orientational decay similar to that of bulk water, but with a generally slower dynamics. Otherwise, at low temperatures, we have two components for the orientational relaxation of bulk water in the low density liquid phase, and only a single component in the high density liquid (HDL) phase, which reflect the symmetry properties of the different phases of the solvent model. In the very dense region of water molecules in the first hydration shell of micelles at low temperatures, we find two components for the orientational relaxation on both liquid phases, one of them much slower than that in the single component of bulk water in the HDL phase. This happens even though our model does not present any hindrance to the water rotational freedom caused by the presence of the amphiphiles.     

Fluctuating hydration structure around nanometer-size hydrophobic solutes. I. Caging and drying around C60 and C60H60 spheres

The hydration structure around nanometer-size hydrophobic solutes is studied with molecular dynamics simulation by taking aqueous solutions of C60 and C60H60 as examples. In the hydration shell around a single C60 or C60H60, dipoles of simulated water molecules tend to be aligned to form the vortexlike coherent pattern which lasts for 100 ps, while individual water molecules stay within the hydration shell for about 10 ps. This structural pattern organized by fluctuating and diffusively moving molecules should be called a "fluctuating cage". In the narrow region between a pair of C60 molecules or a pair of C60H60 molecules, water density strongly fluctuates and is correlated to the mean force between solutes. The fluctuating caging and drying between solutes affect the hydrophobic interaction and dynamical behaviors of solutes.     

Effect of ionic solutes on the hydrogen bond network dynamics of water: power spectral analysis of aqueous NaCl solutions

To understand the modifications of the hydrogen bond network of water by ionic solutes, power spectra as well as static distributions of the potential energies of tagged solvent molecules and solute ions have been computed from molecular dynamics simulations of aqueous NaCl solutions. The key power spectral features of interest are the presence of high-frequency peaks due to localized vibrational modes, the existence of a multiple time scale or 1/falpha frequency regime characteristic of networked liquids, and the frequency of crossover from 1/falpha type behavior to white noise. Hydrophilic solutes, such as the sodium cation and the chloride anion, are shown to mirror the multiple time scale behavior of the hydrogen bond network fluctuations, unlike hydrophobic solutes which display essentially white noise spectra. While the power spectra associated with tagged H2O molecules are not very sensitive to concentration in the intermediate frequency 1/falpha regime, the crossover to white noise is shifted to lower frequencies on going from pure solvent to aqueous alkali halide solutions. This suggests that new and relatively slow time scales enter the picture, possibly associated with processes such as migration of water molecules from the hydration shell to the bulk or conversion of contact ion pairs into solvent-separated ion pairs which translate into variations in equilibrium transport properties of salt solutions with concentration. For anions, cations, and solvent molecules, the trends in the alpha exponents of the multiple time scale region and the self-diffusivities are found to be strongly correlated.     

Mode-coupling study on the dynamics of hydrophobic hydration II: Aqueous solutions of benzene and rare gases

The dynamic properties of both the solute and solvent of the aqueous solution of benzene, xenon and neon are calculated by the mode-coupling theory for molecular liquids based on the interaction-site model. The B-coefficients of the reorientational relaxation and the translational diffusion of the solvent are evaluated from their dependence on the concentration of the solute, and the reorientational relaxation time of water within the hydration shell is estimated based on the two-state model. The reorientational relaxation times of water in the bulk and within the hydration shell, that of solute, and the translational diffusion coefficients of solute and solvent, are calculated at 0-30 degrees C. The temperature dependence of these dynamic properties is in qualitative agreement with that of NMR experiment reported by Nakahara et al. (M. Nakahara, C. Wakai, Y. Yoshimoto and N. Matubayasi, J. Phys. Chem., 1996, 100, 1345-1349, ref. 36), although the agreement of the absolute values is not so good. The B-coefficients of the reorientational relaxation times for benzene, xenon and neon solution are correlated with the hydration number and the partial molar volume of the solute. The proportionality with the latter is better than that with the former. These results support the mechanism that the retardation of the mobility of water is caused by the cavity formation of the solute, as previously suggested by us (T. Yamaguchi, T. Matsuoka and S. Koda, J. Chem. Phys., 2004, 120, 7590-7601, ref. 34), rather than the conventional one that the rigid hydration structure formed around the hydrophobic solute reduces the mobility of water.     

Towards an understanding of the molecular basis of hydrophobicity

Summary The modern view is stressed that the structuring of water around nonpolar solutes, a process called hydrophobic hydration, actually favors the solubility of nonpolar solutes in water, its associated positive free energy of transfer arising from the enthalpic input required to create a cavity in water to accommodate the solute. The results of a series of molecular dynamics simulations of methane in SPC/E water at different temperatures are reported. These results show the existence of a larger fraction of broken hydrogen bonds in the hydration-shell water of the nonpolar solutes with respect to the bulk water, the difference increasing with a rise in temperature. This supports Muller's modified hydration-shell hydrogen-bond model predictions, where hydration-shell water molecules have lower free energies of hydrogen-bond breaking than those in the bulk.This paper is based on a presentation given at the 14th Molecular Graphics and Modelling Society Conference, held in Cairns, Australia, August 27 September 1, 1995.     

Ions in water: the microscopic structure of concentrated hydroxide solutions

Neutron-diffraction data on aqueous solutions of hydroxides, at solute concentrations ranging from 1 solute per 12 water molecules to 1 solute per 3 water molecules, are analyzed by means of a Monte Carlo simulation (empirical potential structure refinement), in order to determine the hydration shell of the OH- in the presence of the smaller alkali metal ions. It is demonstrated that the symmetry argument between H+ and OH- cannot be used, at least in the liquid phase at such high concentrations, for determining the hydroxide hydration shell. Water molecules in the hydration shell of K+ orient their dipole moment at about 45 degrees from the K+-water oxygen director, instead of radially as in the case of the Li+ and Na+ hydration shells. The K+-water oxygen radial distribution function shows a shallower first minimum compared to the other cation-water oxygen functions. The influence of the solutes on the water-water radial distribution functions is shown to have an effect on the water structure equivalent to an increase in the pressure of the water, depending on both ion concentration and ionic radius. The changes of the water structure in the presence of charged solutes and the differences among the hydration shells of the different cations are used to present a qualitative explanation of the observed cation mobility.     

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.     

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.     

How hydrophobic buckminsterfullerene affects surrounding water structure

The hydrophobic hydration of fullerenes in water is of significant interest as the most common Buckminsterfullerene (C60) is a mesoscale sphere; C60 also has potential in pharmaceutical and nanomaterial applications. We use an all-atom molecular dynamics simulation lasting hundreds of nanoseconds to determine the behavior of a single molecule of C60 in a periodic box of water, and compare this to methane. A C60 molecule does not induce drying at the surface; however, unlike a hard sphere methane, a hard sphere C60 solute does. This is due to a larger number of attractive Lennard-Jones interactions between the carbon atom centers in C60 and the surrounding waters. In these simulations, water is not uniformly arranged but rather adopts a range of orientations in the first hydration shell despite the spherical symmetry of both solutes. There is a clear effect of solute size on the orientation of the first hydration shell waters. There is a large increase in hydrogen-bonding contacts between waters in the C60 first hydration shell. There is also a disruption of hydrogen bonds between waters in the first and second hydration shells. Water molecules in the first hydration shell preferentially create triangular structures that minimize the net water dipole near the surface near both the methane and C60 surface, reducing the total energy of the system. Additionally, in the first and second hydration shells, the water dipoles are ordered to a distance of 8 A from the solute surface. We conclude that, with a diameter of approximately 1 nm, C60 behaves as a large hydrophobic solute.     

Communication: On the locality of hydrogen bond networks at hydrophobic interfaces

The formation of structured hydrogen bond networks in the solvation shells immediate to hydrophobic solutes is crucial for a large number of water mediated processes. A long lasting debate in this context regards the mutual influence of the hydrophobic solute into the bulk water and the role of the hydrogen bond network of the bulk in supporting the solvation structure around a hydrophobic molecule. In this context we present a molecular dynamics study of the solvation of various hydrophobic molecules where the effect of different regions around the solvent can be analyzed by employing an adaptive resolution method, which can systematically separate local and nonlocal factors in the structure of water around a hydrophobic molecule. A number of hydrophobic solutes of different sizes and two different model potential interactions between the water and the solute are investigated.     

Temperature and length scale dependence of tetraalkylammonium ion-amide interaction

We have studied the temperature and length scale dependence of the energetics of the pair interaction of well-established hydrophobic solutes tetraalkylammonium bromides with hydrophilic formamide (FA) and hydrophobic hexamethylphosphoric triamide (HMPT). Our results do indicate the anomalous length scale dependence of the tetraalkylammonium cation-amide interaction in water. As the cation size is increased, the unfavorable enthalpy of interaction is increased rather linearly until the maximum is reached, after which there appears to be a reversal of the trend. We believe that this phenomenon arises from the impossibility of water to maintain its H-bond network near large tetraalkylammonium cations that leads to the formation of a somewhat disordered solute hydration shell. The energetic cost for overlapping this shell with the amide hydration shell in water is noticeably smaller than that for tetraalkylammonium cations of a moderated size.     

A new angle on heat capacity changes in hydrophobic solvation

The differential solubility of polar and apolar groups in water is important for the self-assembly of globular proteins, lipid membranes, nucleic acids, and other specific biological structures through hydrophobic and hydrophilic effects. The increase in water's heat capacity upon hydration of apolar compounds is one signature of the hydrophobic effect and differentiates it from the hydration of polar compounds, which cause a decrease in heat capacity. Water structuring around apolar and polar groups is an important factor in their differential solubility and heat capacity effects. Here, it is shown that joint radial/angular distribution functions of water obtained from simulations reveal quite different hydration structures around polar and apolar groups: polar and apolar groups have a deficit or excess, respectively, of "low angle hydrogen bonds". Low angle hydrogen bonds have a larger energy fluctuation than high angle bonds, and analysis of these differences provides a physical reason for the opposite changes in heat capacity and new insight into water structure around solutes and the hydrophobic effect.