Molecular origin of the hydrophobic effect: analysis using the angle-dependent integral equation theory

The molecular origin of the hydrophobic effect is investigated using the angle-dependent integral equation theory combined with the multipolar water model. The thermodynamic quantities of solvation (excess quantities) of a nonpolar solute are decomposed into the translational and orientational contributions. The translational contributions are substantially larger with the result that the temperature dependence of the solute solubility, for example, can well be reproduced by a model simple fluid where the particles interact through strongly attractive potential such as water and the particle size is as small as that of water. The thermodynamic quantities of solvation for carbon tetrachloride, whose molecular size is approximately 1.9 times larger than that of water, are roughly an order of magnitude smaller than those for water and extremely insensitive to the strength of solvent-solvent attractive interaction and the temperature. The orientational contributions to the solvation energy and entropy are further decomposed into the solute-water pair correlation terms and the solute-water-water triplet and higher-order correlation terms. It is argued that the formation of highly ordered structure arising from the enhanced hydrogen bonding does not occur in the vicinity of the solute. Our proposition is that the hydrophobic effect is ascribed to the interplay of the exceptionally small molecular size and the strongly attractive interaction of water, and not necessarily to its hydrogen-bonding properties.     

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.     

How ions affect the structure of water

We model ion solvation in water. We use the MB model of water, a simple two-dimensional statistical mechanical model in which waters are represented as Lennard-Jones disks having Gaussian hydrogen-bonding arms. We introduce a charge dipole into MB waters. We perform (NPT) Monte Carlo simulations to explore how water molecules are organized around ions and around nonpolar solutes in salt solutions. The model gives good qualitative agreement with experiments, including Jones-Dole viscosity B coefficients, Samoilov and Hirata ion hydration activation energies, ion solvation thermodynamics, and Setschenow coefficients for Hofmeister series ions, which describe the salt concentration dependence of the solubilities of hydrophobic solutes. The two main ideas captured here are (1) that charge densities govern the interactions of ions with water, and (2) that a balance of forces determines water structure: electrostatics (water's dipole interacting with ions) and hydrogen bonding (water interacting with neighboring waters). Small ions (kosmotropes) have high charge densities so they cause strong electrostatic ordering of nearby waters, breaking hydrogen bonds. In contrast, large ions (chaotropes) have low charge densities, and surrounding water molecules are largely hydrogen bonded.     

Evaluation of the contribution to hydration of nonelectrolytes from the hydrophobic effect

A new method was suggested for estimating the hydrophobic effect of contributions to the Gibbs energies and enthalpies of hydration of hydrocarbons, inorganic gases and rare gases. In accordance with this method the hydrophobic effect contribution to the Gibbs energy was evaluated from the difference between the hydration Gibbs energy of a solute and the non hydrophobic contribution. To estimate the latter value, the known dependence connecting the Gibbs energies of solvation of a solute in a number of aprotic solvents to the Hildebrand solubility parameter for these solvents was used. The non hydrophobic contribution to the Gibbs energy of hydration was calculated for various solutes from such dependences extended to water as solvent. The Hildebrand solubility parameter for water used in the calculation was corrected for the effect of association through hydrogen bonding. This correction was made by subtraction of the water self-association enthalpy from the enthalpy of vaporization of water. The evaluated Gibbs energies of the hydrophobic effect are positive for saturated hydrocarbons, inorganic gases and rare gases and linearly depend on the solute molecular refraction. The hydrophobic contribution to the hydration enthalpies of the solutes was calculated in the same manner as was made to calculate the hydrophobic contribution to Gibbs energies of hydration. Enthalpies of the hydrophobic effect for the solutes under study are negative.     

PM3-SM3: A general parameterization for including aqueous solvation effects in the PM3 molecular orbital model

Our recently proposed scheme for including aqueous solvation free energies in parameterized NDDO SCF models is extended to the Parameterized Model 3 semiempirical Hamiltonian. The solvation model takes accurate account of the hydrophobic effect for hydrocarbons, as well as electric polarization of the solvent, the free energy of cavitation, and dispersion interactions. Eight heteroatoms are included (along with H and C), and the new model is parameterized accurately for the water molecule itself, which allows meaningful treatments of specifically hydrogen bonded water molecules. The unphysical partial charges on nitrogen atoms predicted by the Parameterized Model 3 Hamiltonian limit the accuracy of the predicted solvation energies for some compounds containing nitrogen, but the model may be very useful for other systems, especially those for which PM3 is preferred over AM1 for the solute properties of the particular system under study. © 1992 by John Wiley & Sons, Inc.     

Microscopic probing of the size dependence in hydrophobic solvation

We report small angle x-ray scattering data demonstrating the direct experimental microscopic observation of the small-to-large crossover behavior of hydrophobic effects in hydrophobic solvation. By increasing the side chain length of amphiphilic tetraalkyl-ammonium (C(n)H(2n+1))(4)N(+) (R(4)N(+)) cations in aqueous solution we observe diffraction peaks indicating association between cations at a solute size between 4.4 and 5 A?, which show temperature dependence dominated by hydrophobic attraction. Using O K-edge x-ray absorption we show that small solutes affect hydrogen bonding in water similar to a temperature decrease, while large solutes affect water similar to a temperature increase. Molecular dynamics simulations support, and provide further insight into, the origin of the experimental observations.     

The effect of aqueous solutions of trimethylamine-N-oxide on pressure induced modifications of hydrophobic interactions

To understand the mechanism of protein protection by the osmolyte trimethylamine-N-oxide (TMAO) at high pressure, using molecular dynamics (MD) simulations, solvation of hydrophobic group is probed in aqueous solutions of TMAO over a wide range of pressures relevant to protein denaturation. The hydrophobic solute considered in this study is neopentane which is a considerably large molecule. The concentrations of TMAO range from 0 to 4 M and for each TMAO concentration, simulations are performed at five different pressures ranging from 1 atm to 8000 atm. Potentials of mean force are calculated and the relative stability of solvent-separated state over the associated state of hydrophobic solute are estimated. Results suggest that high pressure reduces association of hydrophobic solutes. From computations of site-site radial distribution function followed by analysis of coordination number, it is found that water molecules are tightly packed around the nonpolar particle at high pressure and the hydration number increases with increasing pressure. On the other hand, neopentane interacts preferentially with TMAO over water and although hydration of neopentane reduces in presence of this osmolyte, TMAO does not show any tendency to prevent the pressure-induced dispersion of neopentane moieties. It is also observed that TMAO molecules prefer a side-on orientation near the neopentane surface, allowing its oxygen atom to form favorable hydrogen bonds with water while maintaining some hydrophobic contacts with neopentane. Analysis of hydrogen-bond properties and solvation characteristics of TMAO reveals that TMAO can form hydrogen bonds with water and it reduces the identical nearest neighbor water molecules caused by high hydrostatic pressures. Moreover, TMAO enhances life-time of water-water hydrogen bonds and makes these hydrogen bonds more attractive. Implication of these results for counteracting effect of TMAO against protein denaturation at high pressures are discussed.     

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.     

Dependence of the number of hydrogen bonds per water molecule on its distance to a hydrophobic surface and a thereupon-based model for hydrophobic attraction

A water molecule in the vicinity of a hydrophobic surface forms fewer hydrogen bonds than a bulk molecule because the surface restricts the space available for other water molecules necessary for its hydrogen-bonding. In this vicinity, the number of hydrogen bonds per water molecule depends on its distance to the surface. Considering the number of hydrogen bonds per bulk water molecule (available experimentally) as the only reference quantity, we propose an improved probabilistic approach to water hydrogen-bonding that allows one to obtain an analytic expression for this dependence. (The original version of this approach [Y. S. Djikaev and E. Ruckenstein, J. Chem. Phys. 130, 124713 (2009)] provides the number of hydrogen bonds per water molecule in the vicinity of a hydrophobic surface as an average over all possible locations and orientations of the molecule.) This function (the number of hydrogen bonds per water molecule versus its distance to a hydrophobic surface) can be used to develop analytic models for the effect of hydrogen-bonding on the hydration of hydrophobic particles and their solvent-mediated interaction. Presenting a model for the latter, we also examine the temperature effect on the solvent-mediated interaction of two parallel hydrophobic plates.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Filled and empty states of carbon nanotubes in water: Dependence on nanotube diameter, wall thickness and dispersion interactions

We have carried out a series of molecular dynamics simulations of water containing a narrow carbon nanotube as a solute to investigate the filling and emptying of the nanotube and also the modifications of the density and hydrogen bond distributions of water inside and also in the vicinity of the outer surfaces of the nanotube. Our primary goal is to look at the effects of varying nanotube diameter, wall thickness and also solute-solvent interactions on the solvent structure in the confined region also near the outer surfaces of the solute. The thickness of the walls is varied by considering single and multi-walled nanotubes and the interaction potential is varied by tuning the attractive strength of the 12–6 pair interaction potential between a carbon atom of the nanotubes and a water molecule. The calculations are done for many different values of the tuning parameter ranging from fully Lennard-Jones to pure repulsive pair interactions. It is found that both the solvation characteristics and hydrogen bond distributions can depend rather strongly on the strength of the attractive part of the solute-water interaction potential. The thickness of the nanotube wall, however, is found to have only minor effects on the density profiles, hydrogen bond network and the wetting characteristics. This indicates that the long range electrostatic interactions between water molecules inside and on the outer side of the nanotube do not make any significant contribution to the overall solvation structure of these hydrophobic solutes. The solvation characteristics are primarily determined by the balance between the loss of energy due to hydrogen bond network disruption, cavity repulsion potential and offset of the same by attractive component of the solute-water interactions. Our studies with different system sizes show that the essential features of wetting and dewetting characteristics of narrow nanotubes for different diameter and interaction potentials are also present in relatively smaller systems consisting of about five hundred molecules. We dedicate this work to Professor Debashis Mukherjee on his 60th Birthday.     

Hydrogen bond strength and network structure effects on hydration of non-polar molecules

We measure the solvation free energy, Δμ*, for hard spheres and Lennard-Jones particles in a number of artificial liquids made from modified water models. These liquids have reduced hydrogen bond strengths or altered bond angles. By measuring Δμ* for a number of state points at P = 1 bar and different temperatures, we obtain solvation entropies and enthalpies, which are related to the temperature dependence of the solubilities. By resolving the solvation entropy into the sum of the direct solute-solvent interaction and a term depending on the solvent reorganisation enthalpy we show that, although the hydrophobic effect in water at 300 K arises mainly from the small molecular size, its temperature dependence is anomalously low because the reorganisation enthalpy of liquid water is unusually small. We attribute this to the strong tetrahedral network which results from both the molecular geometry and the hydrogen bond strength.     

Structural and dynamic properties of water within the solvation layer around various conformations of the glycine-based polypeptide

Several conformations of the solvated glycine-based polypeptides were investigated using molecular dynamics simulations. Some properties of water in the neighboring space around these molecules were investigated. It was found that water forms a well-defined layer-the first solvation shell-around the peptide molecule, and thickness of this layer is independent of the peptide structure and is equal to approximately 0.28 nm. Within this layer, water molecules show marked orientations relative to a peptide surface. Using the two-particle contribution to entropy as a measure of structural ordering of water, we found that the first solvation shell contributes 95% or more to the total water ordering around the peptide molecule. In investigating the dynamic properties of water, diffusion coefficients and lifetime of the hydrogen bond, clear differences between solvation layer and the bulk water were observed. It was found that the translational diffusion coefficient, D(T), decreases by 30% or more compared to bulk water; also, the lifetime of the water-water hydrogen bond clearly increases. The rotational diffusion coefficient, however, decreases only slightly, no more than approximately 10%. These differences correspond to the slightly higher energy of the hydrogen bond, and to its slightly distorted geometry. Analyzing the translational dynamics of water in the vicinity of the peptide molecule, it was deduced that the structure of the first solvation shell becomes more rigid than the structure of the bulk water. Investigation of a "pure hydrophobic" form of the polypeptide shows that the structure and the properties of water within the solvation shell are predominantly determined by the hydrophobic effect. The specific interactions between water molecules and various charge groups of the peptide molecule modifies this effect only slightly.     

Calculating the Gibbs energy of hydrogen bonding for proton acceptors with a solvent in methanol solutions

We propose a method for calculating the Gibbs energies of hydrogen bonding of solutes with associated solvents via the thermodynamic analysis of experimental values of solvation Gibbs energies. The method is applied to solutions of different proton acceptors in methanol. It is shown that the contribution of hydrogen bonding processes to the solvation Gibbs energy in methanol is in most cases very different in magnitude from the formation Gibbs energy of equimolar complexes of the solute and methanol. We demonstrate the need to include the contributions from solvophobic effects in investigating intermolecular interactions in associated solvents by means of thermodynamic data.     

Dependence of amide vibrations on hydrogen bonding

The effect of hydrogen bonding on the amide group vibrational spectra has traditionally been rationalized by invoking a resonance model where hydrogen bonding impacts the amide functional group by stabilizing its [(-)O-C=NH (+)] structure over the [O=C-NH] structure. However, Triggs and Valentini's UV-Raman study of solvation and hydrogen bonding effects on epsilon-caprolactum, N, N-dimethylacetamide (DMA), and N-methylacetamide (NMA) ( Triggs, N. E.; Valentini, J. J. J. Phys. Chem. 1992, 96, 6922-6931) casts doubt on the validity of this model by demonstrating that, contrary to the resonance model prediction, carbonyl hydrogen bonding does not impact the AmII' frequency of DMA. In this study, we utilize density functional theory (DFT) calculations to examine the impact of hydrogen bonding on the C=O and N-H functional groups of NMA, which is typically used as a simple model of the peptide bond. Our calculations indicate that, as expected, the hydrogen bonding frequency dependence of the AmI vibration predominantly derives from the C=O group, whereas the hydrogen bonding frequency dependence of the AmII vibration primarily derives from N-H hydrogen bonding. In contrast, the hydrogen bonding dependence of the conformation-sensitive AmIII band derives equally from both C=O and N-H groups and thus, is equally responsive to hydrogen bonding at the C=O or N-H site. Our work shows that a clear understanding of the normal mode composition of the amide vibrations is crucial for an accurate interpretation of the hydrogen bonding dependence of amide vibrational frequencies.     

Changes in water structure induced by the guanidinium cation and implications for protein denaturation

The effect of the guanidinium cation on the hydrogen bonding strength of water was analyzed using temperature-excursion Fourier transform infrared spectra of the OH stretching vibration in 5% H 2O/95% D 2O solutions containing a range of different guanidine-HCl and guanidine-HBr concentrations. Our findings indicate that the guanidinium cation causes the water H-bonds in solution to become more linear than those found in bulk water, and that it also inhibits the response of the H-bond network to increased temperature. Quantum chemical calculations also reveal that guanidinium affects both the charge distribution on water molecules directly H-bonded to it as well as the OH stretch frequency of H-bonds in which that water molecule is the donor. The implications of our findings to hydrophobic solvation and protein denaturation are discussed.     

Solute rotation and solvation dynamics in an alcohol-functionalized room temperature ionic liquid

Steady-state and time-resolved fluorescence behaviors of two dipolar solutes, coumarin 153 and 4-aminophthalimide, have been studied in an alcohol-functionalized room-temperature ionic liquid, 1-(hydroxyethyl)-3-methylimidazolium bis(trifluoromethanesulfonyl)imide. The steady-state fluorescence parameters have been exploited for the estimation of the polarity of this ionic liquid and to obtain information on the hydrogen bonding interaction between the ionic liquid and the probe molecules. The time-resolved measurements have been focused on the dynamics of solvation by studying the dynamic Stokes shift in the ps-ns time scale and solute rotation by measuring the time dependence of the fluorescence anisotropy. The time-resolved anisotropy studies reveal a significant slow down of the rotational motion of one of the probe molecules. The time-dependent fluorescence Stokes shift measurements suggest that the time-resolvable part of the dynamics is biphasic in nature, highly dependent on the probe molecule and the ultrafast component is comparatively less than that in other ionic liquids. The influence of the hydrogen bonding interaction between the probe molecules and the ionic liquids on the solute rotation and the various components of the solvation dynamics is carefully analyzed in an attempt to obtain further insight into the mechanism of solvation in these novel media.     

Role of solvation dynamics in excited state proton transfer of 1-naphthol in nanoscopic water clusters formed in a hydrophobic solvent

Excited state proton transfer (ESPT) in biologically relevant organic molecules in aqueous environments following photoexcitation is very crucial as the reorganization of polar solvents (solvation) in the locally excited (LE) state of the organic molecule plays an important role in the overall rate of the ESPT process. A clear evolution of the two photoinduced dynamics in a model ESPT probe 1-naphthol (NpOH) upon ultrafast photoexcitation is the motive of the present study. Herein, the detailed kinetics of the ESPT reaction of NpOH in water clusters formed in hydrophobic solvent are investigated. Distinct values of time constants associated with proton transfer and solvent relaxation have been achieved through picosecond-resolved fluorescence measurements. We have also used a model solvation probe Coumarin 500 (C500) to investigate the dynamics of solvation in the same environmental condition. The temperature dependent picosecond-resolved measurement of ESPT of NpOH and the dynamics of solvation from C500 identify the magnitude of intermolecular hydrogen bonding energy in the water cluster associated with the ultrafast ESPT process.