Origin of hydrophobicity and enhanced water hydrogen bond strength near purely hydrophobic solutes

Does a pair of methane molecules aggregate in water?

Molecular dynamics (MD) simulations of methane-water mixtures were performed using ab initio force fields for the CH₄-H₂O, H₂O-H₂O, and CH₄-CH₄ interactions. Both methane and water molecules were polarizable. From these calculations, the potential of mean force (PMF) between two methane molecules was extracted. Our results are compared with PMFs from a density-functional-theory (DFT) based Born-Oppenheimer type MD (BOMD) simulation, from a Monte Carlo (MC) simulation with ab initio-based force fields, and from MD simulations with empirical force fields. Our PMF is qualitatively similar to that obtained from the simulations with empirical force fields but differs significantly from those resulting from the DFT-BOMD and MC simulations. The depth of the PMF global minimum obtained in the present work is in a much better agreement with the experimental estimate than the result of the DFT-BOMD simulation, possibly due to the inability of DFT to describe the dispersion interactions and the lack of extensive sampling in the BOMD simulations. Our work indicates that, for a pair of methane molecules, there are configurations where the solvent increases the attraction between the solutes, but there are also conformations in which the solvent causes a weak net repulsion. On average, the methane molecules are more likely to be in the configuration where they are separated by a water molecule than in the one in which they are in contact even though the minimum of the PMF at the latter configuration is deeper than that at the former. Finally, we found that the water structure around methane solutes does not show a greater tetrahedral ordering than in neat bulk water.

Hydrophobic chirality amplification in confined water cages

The manipulation of the transition states of a chemical process is essential to achieve the desired selectivity. In particular, transition states of chemical reactions can be significantly modified in a confined environment. We report a catalytic reaction with remarkable amplification of stereochemical information in a confined water cage. Surprisingly, this amplification is significantly dependent on droplet size. This water-induced chirality amplification stems from the hydrophobic hydration effects, which ensures high proximity of the catalyst and substrates presumably at the transition state, leading to higher enantioselectivity. Flow and batch reactors were evaluated to confirm the generality of this water-induced chirality amplification. Our observation on efficient chiral induction in confined water cages might lead to an understanding of the chirality amplification in the prebiotic era, which is a key feature for the chemical evolution of homochirality.

Confined environments may have a significant influence on the transition state and, consequently, on the stereochemical outcome of a chemical reaction. Here, the authors show a droplet size-dependent chirality amplification effect for two model reactions in confined organic droplets surrounded by water cages.

Factors Promoting the Formation of Clathrate-Like Ordering of Water in Biomolecular Structure at Ambient Temperature and Pressure

Clathrate hydrate forms when a hydrophobic molecule is entrapped inside a water cage or cavity. Although biomolecular structures also have hydrophobic patches, clathrate-like water is found in only a limited number of biomolecules. Also, while clathrate hydrates form at low temperature and moderately higher pressure, clathrate-like water is observed in biomolecular structure at ambient temperature and pressure. These indicate presence of other factors along with hydrophobic environment behind the formation of clathrate-like water in biomolecules. In the current study, we presented a systematic approach to explore the factors behind the formation of clathrate-like water in biomolecules by means of molecular dynamics simulation of a model protein, maxi, which is a naturally occurring nanopore and has clathrate-like water inside the pore. Removal of either confinement or hydrophobic environment results in the disappearance of clathrate-like water ordering, indicating a coupled role of these two factors. Apart from these two factors, clathrate-like water ordering also requires anchoring groups that can stabilize the clathrate-like water through hydrogen bonding. Our results uncover crucial factors for the stabilization of clathrate-like ordering in biomolecular structure which can be used for the development of new biomolecular structure promoting clathrate formation.

Understanding the effect of nanoconfinement on the structure of water hydrogen bond networks

Dynamic hydrogen bond trails in water confined between two phospholipid membranes traced by the information flow model.

Contrasting the hydration thermodynamics of methane and methanol

The hydration thermodynamics of methane and methanol depend on the cavity creation work and energy of van der Waals and H-bonding attractions.

Hydrophobic hydration affects growth of clathrate hydrate: insight from an NMR relaxometric and calorimetric study

Water tightly bound to the kinetic inhibitors of tetrahydrofuran hydrate is related to the hydrophobic hydration effect of the inhibitors.

Influence of glycerol on the cooling effect of pair hydrophobicity in water: relevance to proteins’ stabilization at low temperature

Glycerol reduces the cooling effect of pair hydrophobicity (reduction of hydrophobicity with decreasing temperature) in water.

High and low density patches in simulated liquid water

We present insights into the nature of structural heterogeneities in liquid water by characterizing the empty space within the hydrogen bond network. Using molecular dynamics simulations, we show that density fluctuations create regions of empty space characterized by a diverse morphology - from spherical to fractal-like voids. These voids allow for the identification of low and high density patches of the liquid, encompassing short (0.3-0.5 nm) as well as long (1-2 nm) length-scales. In addition, we show that the formation of these patches is coupled to collective fluctuations involving the topology of hydrogen-bonded rings of water molecules. In particular, water molecules in the high density patches tend to be slightly more tetrahedral - which is consistent with the predictions of the hydrophobic effect.

The Chaotropic Effect as an Assembly Motif in Chemistry

Following up on scattered reports on interactions of conventional chaotropic ions (for example, I- , SCN- , ClO4- ) with macrocyclic host molecules, biomolecules, and hydrophobic neutral surfaces in aqueous solution, the chaotropic effect has recently emerged as a generic driving force for supramolecular assembly, orthogonal to the hydrophobic effect. The chaotropic effect becomes most effective for very large ions that extend beyond the classical Hofmeister scale and that can be referred to as superchaotropic ions (for example, borate clusters and polyoxometalates). In this Minireview, we present a continuous scale of water-solute interactions that includes the solvation of kosmotropic, chaotropic, and hydrophobic solutes, as well as the creation of void space (cavitation). Recent examples for the association of chaotropic anions to hydrophobic synthetic and biological binding sites, lipid bilayers, and surfaces are discussed.

Methane Hydration-Shell Structure and Fragility

The influence of oily molecules on the structure of liquid water is a question of importance to biology and geology and many other fields. Previous experimental, theoretical, and simulation studies of methane in liquid water have reached widely conflicting conclusions regarding the structure of hydrophobic hydration-shells. Herein we address this question by performing Raman hydration-shell vibrational spectroscopic measurements of methane in liquid water from -10 °C to 300 °C (at 30 MPa, along a path that parallels the liquid-vapor coexistence curve). We show that, near ambient temperatures, methane's hydration-shell is slightly more tetrahedral than pure water. Moreover, the hydration-shell undergoes a crossover to a more disordered structure above ca. 85 °C. Comparisons with the crossover temperature of aqueous methanol (and other alcohols) reveal the stabilizing influence of an alcohol OH head-group on hydrophobic hydration-shell fragility.

Bridging solution properties to gas hydrate nucleation through guest dynamics

By investigating the aqueous solution properties of several hydrate guests with molecular simulations, we find that with increasing guest concentration, the guest's hydration shell becomes more ordered and the system entropy decreases. A common critical value of the self-diffusion coefficient of different guest molecules is identified, below which hydrates will nucleate very readily.

How do glycerol and dimethyl sulphoxide affect local tetrahedral structure of water around a nonpolar solute at low temperature? Importance of preferential interaction

Glycerol and dimethyl sulphoxide (DMSO) have vital roles in cryoprotection of living cells, tissues, etc. The above action has been directly linked with disruption of hydrogen (H-) bond structure and dynamics of water by these cosolvents at bulk region and around various complex units, such as peptide, amino acid, protein, and lipid membrane. However, the disruption of the local structure of the water solvent around a purely hydrophobic solute is still not studied extensively. The latter is also important in the context of stabilization of protein from cold denaturation. Through all-atom molecular dynamics simulation, we have investigated the comparative effect of glycerol and DMSO on the orientational order of water around a nonpolar solute at -5 °C. A steady reduction of the tetrahedral order of water is observed at bulk (>10 Å distance from the solute) and solute interface (<5.5 Å distance from the solute) with increasing the cosolvent concentration. Contrasting roles of glycerol and DMSO have been evidenced. While DMSO affects the H-bond structure of the interfacial water more than that of the bulk water, glycerol affects the water structure almost uniformly at all regions around the solute. Furthermore, while glycerol helps to retain water molecules at the interface, DMSO significantly reduces the water content in that region. We have put forward a plausible mechanism for these contrasting roles of these cosolvents. The solute-cosolvent hydrophobic-interaction-induced orientational alignment of an interfacial cosolvent molecule determines whether the involvement of the cosolvent molecules in H-bonding with solvent water in the interface is akin to the bulk region or not.

Hydrophobic Hydration and the Effect of NaCl Salt in the Adsorption of Hydrocarbons and Surfactants on Clathrate Hydrates

Adsorption of functional molecules on the surface of hydrates is key in the understanding of hydrate inhibitors. We investigate the adsorption of a hydrocarbon chain, nonionic and ionic surfactants, and ions at the hydrate-aqueous interface. Our results suggest a strong connection between the water ordering around solutes in bulk and the affinity for the hydrates surface. We distinguish two types of water ordering around solutes: (i) hydrophobic hydration where water molecules form a hydrogen bond network similar to clathrate hydrates, and (ii) ionic hydration where water molecules align according to the polarity of an ionic group. The nonionic surfactant and the hydrocarbon chain induce hydrophobic hydration and are favorably adsorbed on the hydrate surface. Adsorption of ions and the ionic headgroups on the hydrate surface is not favorable because ionic hydration and the hydrogen bond structure of hydrates are incompatible. The nonionic surfactant is adsorbed by the headgroup and tail while adsorption of the ionic surfactants is not favorable through the head. Water ordering is analyzed using the hydrogen bond and tetrahedral density profiles as a function of the distance to the chemical groups. The adsorption of solutes is studied through the free energy profiles as a function of the distance to the hydrate surface. Salt lowers the melting temperature of hydrates, disrupts hydrophobic hydration, reduces the solubility of solutes in the aqueous solution, and increases the propensity of solutes to be adsorbed on hydrate surfaces. Our studies are performed by the unbiased and steered molecular dynamics simulations. The results are in line with experiments on the effect of salt and alkanes in hydrate antiagglomeration.

The spatial range of protein hydration

Proteins interact with their aqueous surroundings, thereby modifying the physical properties of the solvent. The extent of this perturbation has been investigated by numerous methods in the past half-century, but a consensus has still not emerged regarding the spatial range of the perturbation. To a large extent, the disparate views found in the current literature can be traced to the lack of a rigorous definition of the perturbation range. Stating that a particular solvent property differs from its bulk value at a certain distance from the protein is not particularly helpful since such findings depend on the sensitivity and precision of the technique used to probe the system. What is needed is a well-defined decay length, an intrinsic property of the protein in a dilute aqueous solution, that specifies the length scale on which a given physical property approaches its bulk-water value. Based on molecular dynamics simulations of four small globular proteins, we present such an analysis of the structural and dynamic properties of the hydrogen-bonded solvent network. The results demonstrate unequivocally that the solvent perturbation is short-ranged, with all investigated properties having exponential decay lengths of less than one hydration shell. The short range of the perturbation is a consequence of the high energy density of bulk water, rendering this solvent highly resistant to structural perturbations. The electric field from the protein, which under certain conditions can be long-ranged, induces a weak alignment of water dipoles, which, however, is merely the linear dielectric response of bulk water and, therefore, should not be thought of as a structural perturbation. By decomposing the first hydration shell into polarity-based subsets, we find that the hydration structure of the nonpolar parts of the protein surface is similar to that of small nonpolar solutes. For all four examined proteins, the mean number of water-water hydrogen bonds in the nonpolar subset is within 1% of the value in bulk water, suggesting that the fragmentation and topography of the nonpolar protein-water interface has evolved to minimize the propensity for protein aggregation by reducing the unfavorable free energy of hydrophobic hydration.

Entropy connects water structure and dynamics in protein hydration layer

The enzyme Candida Antarctica lipase B (CALB) serves here as a model for understanding connections among hydration layer dynamics, solvation shell structure, and protein surface structure. The structure and dynamics of water molecules in the hydration layer were characterized for regions of the CALB surface, divided around each α-helix, β-sheet, and loop structure. Heterogeneous hydration dynamics were observed around the surface of the enzyme, in line with spectroscopic observations of other proteins. Regional differences in the structure of the biomolecular hydration layer were found to be concomitant with variations in dynamics. In particular, it was seen that regions of higher density exhibit faster water dynamics. This is analogous to the behavior of bulk water, where dynamics (diffusion coefficients) are connected to water structure (density and tetrahedrality) by excess (or pair) entropy, detailed in the Rosenfeld scaling relationship. Additionally, effects of protein surface topology and hydrophobicity on water structure and dynamics were evaluated using multiregression analysis, showing that topology has a somewhat larger effect on hydration layer structure-dynamics. Concave and hydrophobic protein surfaces favor a less dense and more tetrahedral solvation layer, akin to a more ice-like structure, with slower dynamics. Results show that pairwise entropies of local hydration layers, calculated from regional radial distribution functions, scale logarithmically with local hydration dynamics. Thus, the Rosenfeld relationship describes the heterogeneous structure-dynamics of the hydration layer around the enzyme CALB. These findings raise the question of whether this may be a general principle for understanding the structure-dynamics of biomolecular solvation.

Temperature-Dependent Hydrophobic Crossover Length Scale and Water Tetrahedral Order

Experimental Raman multivariate curve resolution and molecular dynamics simulations are performed to demonstrate that the vibrational frequency and tetrahedrality of water molecules in the hydration-shells of short-chain alcohols differ from those of pure water and undergo a crossover above 100 °C (at 30 MPa) to a structure that is less tetrahedral than pure water. Our results demonstrate that the associated crossover length scale decreases with increasing temperature, suggesting that there is a fundamental connection between the spectroscopically observed crossover and that predicted to take place around idealized purely repulsive solutes dissolved in water, although the water structure changes in the hydration-shells of alcohols are far smaller than those associated with an idealized "dewetting" transition.

How Water's Properties Are Encoded in Its Molecular Structure and Energies

How are water's material properties encoded within the structure of the water molecule? This is pertinent to understanding Earth's living systems, its materials, its geochemistry and geophysics, and a broad spectrum of its industrial chemistry. Water has distinctive liquid and solid properties: It is highly cohesive. It has volumetric anomalies-water's solid (ice) floats on its liquid; pressure can melt the solid rather than freezing the liquid; heating can shrink the liquid. It has more solid phases than other materials. Its supercooled liquid has divergent thermodynamic response functions. Its glassy state is neither fragile nor strong. Its component ions-hydroxide and protons-diffuse much faster than other ions. Aqueous solvation of ions or oils entails large entropies and heat capacities. We review how these properties are encoded within water's molecular structure and energies, as understood from theories, simulations, and experiments. Like simpler liquids, water molecules are nearly spherical and interact with each other through van der Waals forces. Unlike simpler liquids, water's orientation-dependent hydrogen bonding leads to open tetrahedral cage-like structuring that contributes to its remarkable volumetric and thermal properties.

Addressing a Common Misconception: Ammonium Acetate as Neutral pH "Buffer" for Native Electrospray Mass Spectrometry

Native ESI-MS involves the transfer of intact proteins and biomolecular complexes from solution into the gas phase. One potential pitfall is the occurrence of pH-induced changes that can affect the analyte while it is still surrounded by solvent. Most native ESI-MS studies employ neutral aqueous ammonium acetate solutions. It is a widely perpetuated misconception that ammonium acetate buffers the analyte solution at neutral pH. By definition, a buffer consists of a weak acid and its conjugate weak base. The buffering range covers the weak acid pK_a ± 1 pH unit. NH₄⁺ and CH₃-COO^- are not a conjugate acid/base pair, which means that they do not constitute a buffer at pH 7. Dissolution of ammonium acetate salt in water results in pH 7, but this pH is highly labile. Ammonium acetate does provide buffering around pH 4.75 (the pK_a of acetic acid) and around pH 9.25 (the pK_a of ammonium). This implies that neutral ammonium acetate solutions electrosprayed in positive ion mode will likely undergo acidification down to pH 4.75 ± 1 in the ESI plume. Ammonium acetate nonetheless remains a useful additive for native ESI-MS. It is a volatile electrolyte that can mimic the solvation properties experienced by proteins under physiological conditions. Also, a drop from pH 7 to around pH 4.75 is less dramatic than the acidification that would take place in pure water. It is hoped that the habit of referring to pH 7 solutions as ammonium acetate "buffer" will disappear from the literature. Ammonium acetate "solution" should be used instead. Graphical Abstract ᅟ.

Calcium-Mediated Control of S100 Proteins: Allosteric Communication via an Agitator/Signal Blocking Mechanism

Allosteric proteins possess dynamically coupled residues for the propagation of input signals to distant target binding sites. The input signals usually correspond to "effector is present" or "effector is not present". Many aspects of allosteric regulation remain incompletely understood. This work focused on S100A11, a dimeric EF-hand protein with two hydrophobic target binding sites. An annexin peptide (Ax) served as the target. Target binding is allosterically controlled by Ca²⁺ over a distance of ∼26 Å. Ca²⁺ promotes formation of a [Ca₄ S100 Ax₂] complex, where the Ax peptides are accommodated between helices III/IV and III'/IV'. Without Ca²⁺ these binding sites are closed, precluding interactions with Ax. The allosteric mechanism was probed by microsecond MD simulations in explicit water, complemented by hydrogen exchange mass spectrometry (HDX/MS). Consistent with experimental data, MD runs in the absence of Ca²⁺ and Ax culminated in target binding site closure. In simulations on [Ca₄ S100] the target binding sites remained open. These results capture the essence of allosteric control, revealing how Ca²⁺ prevents binding site closure. Both HDX/MS and MD data showed that the metalation sites become more dynamic after Ca²⁺ loss. However, these enhanced dynamics do not represent the primary trigger of the allosteric cascade. Instead, a labile salt bridge acts as an incessantly active "agitator" that destabilizes the packing of adjacent residues, causing a domino chain of events that culminates in target binding site closure. This agitator represents the starting point of the allosteric signal propagation pathway. Ca²⁺ binding rigidifies elements along this pathway, thereby blocking signal transmission. This blocking mechanism does not conform to the commonly held view that allosteric communication pathways generally originate at the sites where effectors interact with the protein.

Mapping Hydration Water around Alcohol Chains by THz Calorimetry

THz spectroscopy was used to probe changes that occur in the dynamics of the hydrogen bond network upon solvation of alcohol chains. The THz spectra can be decomposed into the spectrum of bulk water, tetrahedral hydration water, and more disordered (or interstitial) hydration water. The tetrahedrally ordered hydration water exhibits a band at 195 cm⁻¹ and is localized around the hydrophobic moiety of the alcohol. The interstitial component yields a band at 164 cm⁻¹ which is associated with hydration water in the first hydration shell. These temperature‐dependent changes in the low‐frequency spectrum of solvated alcohol chains can be correlated with changes of heat capacity, entropy, and free energy upon solvation. Surprisingly, not the tetrahedrally ordered component but the interstitial hydration water is found to be mainly responsible for the temperature‐dependent change in ΔC_p and ΔG. The solute‐specific offset in free energy is attributed to void formation and scales linearly with the chain length.

Influence of Order within Nonpolar Monolayers on Hydrophobic Interactions

We report an experimental investigation of the influence of molecular-level order (crystallinity) within nonpolar monolayers on hydrophobic interactions. The measurements were performed using gold film-supported monolayers formed from alkanethiols (CH₃(CH₂)_nSH, with n ranging from 3 to 17), which we confirmed by using polarization-modulation infrared reflection-adsorption spectroscopy to exhibit chain-length-dependent order (methylene peak moves from 2926 to 2919 cm^-1, corresponding to a transition from liquid- to quasi-crystalline-like order) in the absence of substantial changes in chain density (constant methyl peak intensity). By using monolayer-covered surfaces immersed in either aqueous triethanolamine (TEA, 10 mM, pH 7.0) or pure methanol, we quantified hydrophobic and van der Waals contributions to adhesive interactions between identical pairs of surfaces (measured using an atomic force microscope) as a function of the length and order of the aliphatic chains within the monolayers. In particular, we measured pull-off forces arising from hydrophobic adhesion to increase in a nonlinear manner with chain length (abrupt increase between n = 5 and 6 from 2.1 ± 0.3 to 14.1 ± 0.7 nN) and to correlate closely with a transition from a liquid-like to crystalline-like monolayer phase. In contrast, adhesion in methanol increased gradually with chain length from 0.3 ± 0.1 to 3.2 ± 0.3 nN for n = 3 to 7 and then did not change further with an increase in chain length. These results lead to the hypothesis that order within nonpolar monolayers influences hydrophobic interactions. Additional support for this hypothesis was obtained from measurements reported in this paper using long-chain alkanethiols (ordered) and alkenethiols (disordered). The results are placed into the context of recent spectroscopic studies of hydrogen bonding of water at nonpolar monolayers. Overall, our study provides new insight into factors that influence hydrophobic interactions at nonpolar monolayers.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.