Molecular origin of anticooperativity in hydrophobic association

The structuring of water molecules in the vicinity of nonpolar solutes is responsible for hydrophobic hydration and association thermodynamics in aqueous solutions. Here, we studied the potential of mean force (PMF) for the formation of a dimer and trimers of methane molecules in three specific configurations in explicit water to explain multibody effects in hydrophobic association on a molecular level. We analyzed the packing and orientation of water molecules in the vicinity of the solute to explain the effect of ordering of the water around nonpolar solutes on many-body interactions. Consistent with previous theoretical studies, we observed cooperativity, manifested as a reduction of the height of the desolvation barrier for the trimer in an isosceles triangle geometry, but for linear trimers, we observed only anticooperativity. A simple mechanistic picture of hydrophobic association is drawn. The free energy of hydrophobic association depends primarily on the difference in the number of water molecules in the first solvation shell of a cluster and that in the monomers of a cluster; this can be approximated by the molecular surface area. However, there are unfavorable electrostatic interactions between the water molecules from different parts of the solvation shell of a trimer because of their increased orientation induced by the nonpolar solute. These electrostatic interactions make an anticooperative contribution to the PMF, which is clearly manifested for the linear trimer where the multibody contribution due to changes in the molecular surface area is equal to zero. The information theory model of hydrophobic interactions of Hummer et al. also explains the anticooperativity of hydrophobic association of the linear trimers; however, it predicts anticooperativity with a qualitatively identical distance dependence for nonlinear trimers, which disagrees with the results of simulations.     

Hydration of aprotic donor solvents studied by means of FTIR spectroscopy

The paper attempts to explain the mutual influence of nonpolar and electron-donor groups on solute hydration, the problem of big importance for biological aqueous systems. Aprotic organic solvents have been used as model solutes, differing in electron-donating power. Hydration of acetonitrile, acetone, 2-butanone, and triethylamine has been studied by HDO and (partially) H2O spectra. The quantitative version of difference spectra method has been applied to determine solute-affected water spectra. Analysis of the data suggests that solvent-water interaction via the donor center of the solute is averaged between water-water interactions around the solute. Such behavior can be simply explained by the model of solute rotating in a cavity of water structure, which is formed by clathratelike hydrogen-bonded water network. On the basis of the band shape of solute-affected HDO spectra and the corresponding distribution of intermolecular distances, the criterion for hydrophobic type hydration has been proposed. From that point of view, all the studied solutes could be treated as hydrophobic ones. The limiting band position and the corresponding intermolecular distance of affected water, gained with increasing electron-donating power of solutes, has been inferred from the data obtained. These observations are important for interpretation of vibrational spectra of water as well as for volumetric measurements of solutions. The simple model of hydration, proposed to better justify the results, connects the values obtained from the methods providing microscopic and macroscopic characteristics of the system studied.     

Structural characteristics of a 0.23 mole fraction aqueous solution of tetrahydrofuran at 20 degrees C

Hydrogen/deuterium isotopic substitution neutron diffraction techniques were used to measure the structural correlation functions in a 0.23 mole fraction solution of tetrahydrofuran in water at room temperature. Empirical potential structure refinement (EPSR) was used to build a three-dimensional model of the liquid structure that is consistent with the experimental data. Detailed analysis shows a preference for nonpolar interactions between the cyclic ether molecules plus polar interactions between the ether and solvent water and hydrophobic hydration of the nonpolar regions of the solute. The increase in the number of hydrogen-bond-acceptor sites relative to the number of hydrogen-bond-donor sites in this system, compared to the balanced situation that would be found in pure water, has a marked compressive effect on the structure of the solvent. Despite the small size of the solvent water molecules, the 0.23 mole fraction aqueous solution is still found to contain small voids akin to those in pure liquid tetrahydrofuran. In contrast to the positive surface charge of the voids in the pure system, the average void in this aqueous solution is found to have a net negative charge. This is due to contributions from the water oxygen atoms that are negatively polarized by their intramolecular bonding.     

Orientational dynamics of water trapped between two nanoscopic hydrophobic solutes: A molecular dynamics simulation study

We investigate thoroughly the effect of confinement and solute topology on the orientational dynamics of water molecule in the interplate region between two nanoscopic hydrophobic paraffinlike plates. Results are obtained from molecular dynamics simulations of aqueous solutions of paraffinlike plates in the isothermal-isobaric ensemble. An analysis of survival time auto correlation function shows that the residence time of the water molecule in the confined region between two model nanoscopic hydrophobic plates depends on solute surface topology (intermolecular distance within the paraffinlike plate). As expected, the extent of confinement also changes the residence time of water molecules considerably. Orientational dynamics was analyzed along three different directions, viz., dipole moment, HH, and perpendicular to molecular plane vectors. It has been demonstrated that the rotational dynamics of the confined water does not follow the Debye rotational diffusion model, and surface topology of the solute plate and the extent of confinement have considerable effect on the rotational dynamics of the confined water molecules.     

Different Views on the Stability of Protein Conformations and Hydrophobic Effects

Apparently contradictory statements about the thermodynamics of aqueous protein solutions and of hydrophobic effect are quoted and discussed. Some credibility is found in the divergent points of view and it is pointed out that they focus attention on different aspects of the complicated conditions in aqueous solutions, some of which are more important than others for the stability of protein conformations.The importance of characteristics of solvent water is emphasized, in particular (1) the strong mutual cohesion of water molecules, and (2) structural changes of water induced by (nonpolar) solute molecules. It is stressed that consideration of only one of these effects and an inexpedient choice of standard states are origins of confusion in the literature about aqueous systems. A simple approach to hydrophobic effects considering both of the above mentioned effects, is proposed.     

The Hydrophobic Effects: Our Current Understanding

Hydrophobic interactions are involved in and believed to be the fundamental driving force of many chemical and biological phenomena in aqueous environments. This review focuses on our current understanding on hydrophobic effects. As a solute is embedded into water, the interface appears between solute and water, which mainly affects the structure of interfacial water (the topmost water layer at the solute/water interface). From our recent structural studies on water and air-water interface, hydration free energy is derived and utilized to investigate the origin of hydrophobic interactions. It is found that hydration free energy depends on the size of solute. With increasing the solute size, it is reasonably divided into initial and hydrophobic solvation processes, and various dissolved behaviors of the solutes are expected in different solvation processes, such as dispersed and accumulated distributions in solutions. Regarding the origin of hydrophobic effects, it is ascribed to the structural competition between the hydrogen bondings of interfacial and bulk water. This can be applied to understand the characteristics of hydrophobic interactions, such as the dependence of hydrophobic interactions on solute size (or concentrations), the directional natures of hydrophobic interactions, and temperature effects on hydrophobic interactions.     

Dynamics of protein and peptide hydration

Biological processes often involve the surfaces of proteins, where the structural and dynamic properties of the aqueous solvent are modified. Information about the dynamics of protein hydration can be obtained by measuring the magnetic relaxation dispersion (MRD) of the water (2)H and (17)O nuclei or by recording the nuclear Overhauser effect (NOE) between water and protein protons. Here, we use the MRD method to study the hydration of the cyclic peptide oxytocin and the globular protein BPTI in deeply supercooled solutions. The results provide a detailed characterization of water dynamics in the hydration layer at the surface of these biomolecules. More than 95% of the water molecules in contact with the biomolecular surface are found to be no more than two-fold motionally retarded as compared to bulk water. In contrast to small nonpolar molecules, the retardation factor for BPTI showed little or no temperature dependence, suggesting that the exposed nonpolar residues do not induce clathrate-like hydrophobic hydration structures. New NOE data for oxytocin and published NOE data for BPTI were analyzed, and a mutually consistent interpretation of MRD and NOE results was achieved with the aid of a new theory of intermolecular dipolar relaxation that accounts explicitly for the dynamic perturbation at the biomolecular surface. The analysis indicates that water-protein NOEs are dominated by long-range dipolar couplings to bulk water, unless the monitored protein proton is near a partly or fully buried hydration site where the water molecule has a long residence time.     

Self-assembly does not account for the hydrophobic effect

Marmur has claimed that large values of activity coefficients for nonelectrolytes, particularly in the context of hydrophobic interactions between solutes in aqueous solution at ambient temperature and pressure, cannot be accounted for by thermodynamics, and has suggested that association (self-assembly) of solute molecules in solution solves this dilemma. We show that the analysis of Marmur is incorrect, specifically because the equilibrium in solution between monomeric solute molecules and associated solute molecules is entirely ignored. We show further that activity coefficients such as that for nitromethane solute in hexane solvent, 39.7, and that for solute hexane in solvent water, 4.48 x 10(5), can be calculated as 31.9 and 4.71 x 10(5), respectively, by methods based on well-known molecule-molecule interactions. No assumption of self-assembly is required.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

Ion pairs in aqueous electrolyte microdrops under conditions of a flat nanopore

The molecular mechanisms of aqueous solvent penetration into a flat nanopore with hydrophobic structureless walls containing a Na⁺Cl^? ion pair with nonfixed distance between ions is studied by computer simulations. A detailed many-body polycenter model of intermolecular interactions calibrated with respect to experimental data for the free energy of attachment of water vapor molecules and quantum-chemical calculations in clusters is used. The ion pair hydration results in its decomposition. Drawing the molecules into the gap between ions makes easier penetration of solvent and filling of the nanopore with electrolyte. The ion-pair dissociation is accompanied by dramatic changes in the chemical potential of molecules and electric properties of the whole system. The thermodynamic characteristics of decomposition are stable as regards variations in the pore width. The post-decomposition electric polarizability demonstrates strong anisotropy associated with the nanopore flatness.     

Crystallography Without a Lattice

Molecular structure determination in crystals depends on the presence of the crystal lattice. In liquids, there is no underlying lattice, so phase relationships between scattered X-rays or neutrons are not preserved. Hence, only pair-distance distributions can be obtained from simple diffraction experiments. Advances in the past 30 years in radiation sources, instrumentation, and computing have enabled us to go beyond this apparent limitation. We can now obtain detailed structural information on relatively complex liquids and thus see clearly for the first time how molecules actually interact in solution, and how the solvent is perturbed by the presence of the solute. Using as an example recent work on aqueous solutions of tertiary butanol, the ways in which we can obtain high quality structural information in the absence of a lattice are summarized. Not only can we see how intermolecular interactions are influenced by changes in temperature and concentration, we can also begin to see where the structural source of the entropic driving force for the hydrophobic interaction may be found.     

Thermodynamics of Interaction between Some Cellulose Ethers and SDS by Titration Microcalorimetry

The interaction between certain nonionic cellulose ethers (ethyl hydroxyethyl cellulose and hydroxypropyl methyl cellulose) and sodium dodecyl sulphate (SDS) has been investigated using isothermal titration microcalorimetry at temperatures between 25-50 degrees C. The observed heat flow curves have been interpreted in terms of a plausible mechanism of the interaction of the substituent groups with SDS monomers and clusters. The data have been related to changes occuring in the system at the macro- and microscopic levels with the addition of surfactants and with temperature. The process consists predominantly of polymer-surfactant interactions initially and surfactant-surfactant interactions at the later stages. A phenomenological model of the cooperative interaction (adsorption) process has been derived, and earlier published equilibrium binding data have been used to recover binding constants and Gibbs energy changes for this process. The adsorption enthalpies and entropies have been recovered along with the heat capacity change. The enthalpic cost of confining the nonpolar regions of the polymers in surfactant clusters is high, but the entropy gain from release of hydration shell water molecules as well as increased freedom of movement of these nonpolar regions in the clusters gives the process a strong entropic driving force. The process is entropy-driven initially and converts to being both enthalpy and entropy-driven at high SDS concentrations. An enthalpy-entropy compensation behavior is seen. Strongly negative heat capacity changes have been obtained resulting from the transfer of nonpolar groups from aqueous into nonpolar environments, as well as a reduction of conformational domains that the chains can populate. Changes in these two components cause the heat capacity change to become less negative at the higher binding levels. The system can be classified as exhibiting nonclassical hydrophobic binding at the later stages of binding. Copyright 1999 Academic Press.     

Spatial hydration structures and dynamics of phenol in sub- and supercritical water

The hydration structures and dynamics of phenol in aqueous solution at infinite dilution are investigated using molecular-dynamics simulation technique. The simulations are performed at several temperatures along the coexistence curve of water up to the critical point, and above the critical point with density fixed at 0.3 g/cm3. The hydration structures of phenol are characterized using the radial, cylindrical, and spatial distribution functions. In particular, full spatial maps of local atomic (solvent) density around a solute molecule are presented. It is demonstrated that in addition to normal H bonds with hydroxyl group of phenol, water forms pi-type complexes with the center of the benzene ring, in which H2O molecules act as H-bond donor. At ambient conditions phenol is solvated by 38 water molecules, which make up a large hydrophobic cavity, and forms on average 2.39 H bonds (1.55 of which are due to the hydroxyl group-water interactions and 0.84 are due to the pi complex) with its hydration shell. As temperature increases, the hydration structure of phenol undergoes significant changes. The disappearance of the pi-type H bonding is observed near the critical point. Self-diffusion coefficients of water and phenol are also calculated. Dramatic increase in the diffusivity of phenol in aqueous solution is observed near the critical point of simple point-charge-extended water and is related to the changes in water structure at these conditions.     

Ultrafast Vibrational Dynamics and Local Interactions of Hydrated DNA

Biochemical processes occur mainly in aqueous environments, where interactions with water molecules play a key role for both the structure and function of biomolecules. Deoxyribonucleic acid (DNA), the basic carrier of genetic information, is characterized by an equilibrium double helix structure which is held together by intermolecular hydrogen bonds between base pairs and hydrated by an environment of water molecules with fluctuating hydrogen bonds. Basic vibrational motions of hydrated DNA and the fastest changes in the DNA–water interactions and hydration geometries occur in less than 1 ps. These processes can be accessed by mapping the vibrational dynamics of DNA and water in a time‐resolved way by nonlinear ultrafast vibrational spectroscopy. Recent studies provide a detailed understanding of DNA vibrations and their dynamics, and give insight into nonequilibrium properties and structures of hydrated DNA.     

Study on Volumetric,Compressibility and Viscometric Behavior of Cationic Surfactants (CTAB and DTAB) in Aqueous Glycyl Dipeptide: A Thermo-Acoustic Approach

This study aims to understand how glycyl dipeptide affected the compressibility, volumetric behavior and viscometric behavior of the cationic surfactants CTAB (Cetyltrimethylammonium bromide) and DTAB (dodecyltrimethylammonium bromide). Information on solute–solute, solute–solvent, and solvent–solvent interactions has been inferred using the quantification of density (ρ), speed of sound (u) and viscosity in aqueous media containing glycyl dipeptide in the temperature range 293.15–313.15 K at an interval of 5 K. The data from the aforementioned research have been used to enumerate numerous volumetric and compressibility metrics that aid in the collection of information about the interactional behavior of the system under consideration. The study suggests that CTAB interacts strongly compared to DTAB with dipeptide, and it also significantly dehydrates glycyl dipeptide. The difference in water–water interactions caused by the loss of hydrophobic hydration of the surfactant molecules upon the addition of cationic surfactants may be the cause of the variation in determined parameters with surfactant concentration. Consideration of the structural rearrangement of molecules that may occur in the system has been used to explain the results of viscosity and computed factors related to viscosity. The patterns of competitive intermolecular interactions in the ternary (dipeptide + water + surfactant) system have been used to analyze the trends of all the parameters. The study may be helpful to understand the stability and structural changes in protein–surfactant systems mediated through various interactions that may be present in the system.     

The effect of pressure on the hydration structure around hydrophobic solute: a molecular dynamics simulation study

Molecular dynamics simulations are performed to study the effects of pressure on the hydrophobic interactions between neopentane molecules immersed in water. Simulations are carried out for five different pressure values ranging from 1 atm to 8000 atm. From potential of mean force calculations, we find that with enhancement of pressure, there is decrease in the well depth of contact minimum (CM) and the relative stability of solvent separated minimum over CM increases. Lower clustering of neopentane at high pressure is also observed in association constant and cluster-structure analysis. Selected site-site radial distribution functions suggest efficient packing of water molecules around neopentane molecules at elevated pressure. The orientational profile calculations of water molecules show that the orientation of water molecules in the vicinity of solute molecule is anisotropic and this distribution becomes flatter as we move away from the solute. Increasing pressure slightly changes the water distribution. Our hydrogen bond properties and dynamics calculations reveal pressure-induced formation of more and more number of water molecules with five and four hydrogen bond at the expense of breaking of two and three hydrogen bonded water molecules. We also find lowering of water-water continuous hydrogen bond lifetime on application of pressure. Implication of these results for relative dispersion of hydrophobic molecules at high pressure are discussed.     

Clustering of Lennard-Jones particles in water: temperature and pressure effects

While the hydrophobic effect is, for many systems, one of the most relevant interactions, it may be said that in the case of biological systems this effect becomes of determinant importance. Although the matter has been analyzed extensively, certain aspects are yet to be elucidated. Hence, the study on the behavior of the hydrophobic effect with temperature, and particularly with pressure deserves further investigation; model systems may help us in the task. We have analyzed the behavior of Lennard-Jones particles in water by means of molecular dynamics simulation under different conditions of size, concentration, temperature, and pressure. Following the formation of particle aggregates we can observe an increase of the hydrophobic effect with temperature and a strong weakening of the effect at high pressures. The results agree with the experimental evidence and show the ability of molecular dynamics simulation to account for the behavior of nonpolar substances under different conditions, provided that the intermolecular interactions used are adequate.     

Temperature-responsive chromatography for the separation of biomolecules

Temperature-responsive chromatography for the separation of biomolecules utilizing poly(N-isopropylacrylamide) (PNIPAAm) and its copolymer-modified stationary phase is performed with an aqueous mobile phase without using organic solvent. The surface properties and function of the stationary phase are controlled by external temperature changes without changing the mobile-phase composition. This analytical system is based on nonspecific adsorption by the reversible transition of a hydrophilic-hydrophobic PNIPAAm-grafted surface. The driving force for retention is hydrophobic interaction between the solute molecules and the hydrophobized polymer chains on the stationary phase surface. The separation of the biomolecules, such as nucleotides and proteins was achieved by a dual temperature- and pH-responsive chromatography system. The electrostatic and hydrophobic interactions could be modulated simultaneously with the temperature in an aqueous mobile phase, thus the separation system would have potential applications in the separation of biomolecules. Additionally, chromatographic matrices prepared by a surface-initiated atom transfer radical polymerization (ATRP) exhibit a strong interaction with analytes, because the polymerization procedure forms a densely packed polymer, called a polymer brush, on the surfaces. The copolymer brush grafted surfaces prepared by ATRP was an effective tool for separating basic biomolecules by modulating the electrostatic and hydrophobic interactions. Applications of thermally responsive columns for the separations of biomolecules are reviewed here.     

Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. 1. Approximate expression for the free energy of hydrophobic association based on a Gaussian-overlap model

A physics-based model is proposed to derive approximate analytical expressions for the cavity component of the free energy of hydrophobic association of spherical and spheroidal solutes in water. The model is based on the difference between the number and context of the water molecules in the hydration sphere of a hydrophobic dimer and of two isolated hydrophobic solutes. It is assumed that the water molecules touching the convex part of the molecular surface of the dimer and those in the hydration spheres of the monomers contribute equally to the free energy of solvation, and those touching the saddle part of the molecular surface of the dimer result in a more pronounced increase in free energy because of their more restricted mobility (entropy loss) and fewer favorable electrostatic interactions with other water molecules. The density of water in the hydration sphere around a single solute particle is approximated by the derivative of a Gaussian centered on the solute molecule with respect to its standard deviation. On the basis of this approximation, the number of water molecules in different parts of the hydration sphere of the dimer is expressed in terms of the first and the second mixed derivatives of the two Gaussians centered on the first and second solute molecules, respectively, with respect to the standard deviations of these Gaussians, and plausible analytical expressions for the cavity component of the hydrophobic-association energy of spherical and spheroidal solutes are introduced. As opposed to earlier hydration-shell models, our expressions reproduce the desolvation maxima in the potentials of mean force of pairs of nonpolar solutes in water, and their advantage over the models based on molecular-surface area is that they have continuous gradients in the coordinates of solute centers.     

Deuterium isotope effects on hydrophobic interactions: the importance of dispersion interactions in the hydrophobic phase

Hydrogen/deuterium isotope effects on hydrophobic binding were examined by means of reversed-phase chromatographic separation of protiated and deuterated isotopologue pairs for a set of 10 nonpolar and low-polarity compounds with 10 stationary phases having alkyl and aryl groups bonded to the silica surface. It was found that protiated compounds bind to nonpolar moieties attached to silica more strongly than deuterated ones, demonstrating that the CH/CD bonds of the solutes are weakened or have less restricted motions when bound in the stationary phase compared with the aqueous solvent (mobile phase). The interactions responsible for binding have been further characterized by studies of the effects of changes in mobile phase composition, temperature dependence of binding, and QSRR (quantitative structure-chromatographic retention relationship) analysis, demonstrating the importance of enthalpic effects in binding and differentiation between the isotopologues. To explain our results showing the active role of the hydrophobic (stationary) phase we propose a plausible model that includes specific contributions from aromatic edge-to-face attractive interactions and attractive interactions of aliphatic groups with the pi clouds of aromatic groups present as the solute or in the stationary phase.