Approaches to hydration, old and new: Insights through Hofmeister effects

Hydration effects in colloidal interactions or problems involving electrolytes are usually taken care of by effective electrostatic potentials that subsume notions like hydrated ion size, Gurney potentials, soft and hard, chaotropic and cosmotropic ions, and inner and outer Helmholtz planes. Quantum fluctuation (dispersion) forces between ions and between ions and surfaces are missing from classical theories, at least not explicit in standard approaches to hydration. This paper outlines an evolving back-to-basics approach that allows these ion specific forces to be included in theories quantitatively. In this approach ab initio quantum mechanics is used to calculate dynamic polarisabilities of ions and to quantify bare ion radii. The ionic dispersion potentials between ions, and between ions and surfaces in water can then be given explicit analytic form from an extension of Lifshitz theory. They are included in the theory along with electrostatic potentials. In a first stage the primitive (continuum solvent) model provides a skeletal theory on which to build in hydration. Extension of the ab initio calculations to include “dressed” ions, i.e. water hydration shells for cosmotropic ions, quadrupolar and octupolar polarisability contributions and; for colloids, allowance for a surface hydration layer, permits quantification of Hofmeister effects and Gurney potentials. With these extensions, primary hydration forces (short range repulsion) arise due to an interplay between surface hydration layers and specific ion interactions. Apparent longer range “secondary hydration forces” are shown to be a consequence of ion-surface dispersion interactions and are not true “hydration forces”.     

An <Emphasis Type="Italic">ab initio</Emphasis>and AIM investigation into the hydration of 2-thioxanthine

Background  

Hydration is a universal phenomenon in nature. The interactions between biomolecules and water of hydration play a pivotal role in molecular biology. 2-Thioxanthine (2TX), a thio-modified nucleic acid base, is of significant interest as a DNA inhibitor yet its interactions with hydration water have not been investigated either computationally or experimentally. Here in, we reported an ab initio study of the hydration of 2TX, revealing water can form seven hydrated complexes.     

Water-induced interactions between carbon nanoparticles

Molecular dynamics simulations were carried out in order to study the hydration of C60 fullerenes, carbon nanotubes, and graphene sheets in aqueous solution and the nature of water-induced interactions between these carbon nanoparticles. The hydration of these nonpolar carbon nanoparticles does not exhibit classical hydrophobic character due to the high density of surface atoms (carbon) resulting in strong water-surface dispersion interactions. Water was found to wet the nanoparticle surfaces independent of nanoparticle surface curvature, with the decrease in the extent of water-water hydrogen bonding with decreasing surface curvature being offset by stronger water-surface interactions. While all carbon nanoparticles investigated are anticipated to aggregate in water due to strong direct nanoparticle-nanoparticle interactions, the water-induced interactions between nanoparticles were found to be repulsive and, in contrast to the wetting behavior, were observed to exhibit strong dependence on surface curvature. The strength of the water-induced interaction between carbon nanoparticles was found to correlate well with the number of hydration water molecules displaced upon particle aggregation, which, relative to the amount of direct nanoparticle-nanoparticle contact engendered upon aggregation, decreases with decreasing surface curvature.     

Role of attractive methane-water interactions in the potential of mean force between methane molecules in water

On the basis of a Gaussian quasichemical model of hydration, a model of non-van der Waals character, we explore the role of attractive methane-water interactions in the hydration of methane and in the potential of mean force between two methane molecules in water. We find that the hydration of methane is dominated by packing and a mean-field energetic contribution. Contributions beyond the mean-field term are unimportant in the hydration phenomena for a hydrophobic solute such as methane. Attractive solute-water interactions make a net repulsive contribution to these pair potentials of mean force. With no conditioning, the observed distributions of binding energies are super-Gaussian and can be effectively modeled by a Gumbel (extreme value) distribution. This further supports the view that the characteristic form of the unconditioned distribution in the high-epsilon tail is due to energetic interactions with a small number of molecules. Generalized extreme value distributions also effectively model the results with minimal conditioning, but in those cases the distributions are sufficiently narrow that the details of their shape are not significant.     

Specific interactions and single-lon activity coefficients in mixed electrolyte solutions

A reexamination of Guggenheim's theory of specific interactions for dilute aqueous electrolyte solutions shows that the parameters of the theory should include terms for interactions between like-charged ions. The equations developed are analogous to those in Pitzer's treatment of concentrated solutions, and the inclusion of ion-solvent interaction parameters allows comparison with the Stokes and Robinson hydration theory treatment. As with hydration theory, specific-interaction theory predicts differences in single-ion activity coefficients.     

尿素/氧化三甲胺混合溶剂影响单壁碳纳米管内部水合性质的分子动力学模拟

尿素是早已被人们认识的蛋白质变性剂,而氧化三甲胺则是最常用的蛋白质结构保护剂。虽然多年来被广泛应用在生物实验中,但是它们是如何在蛋白质结构形成中起作用,特别是氧化三甲胺是如何在高浓度尿素环境中起到抑制尿素蛋白变性作用的分子机制,至今仍然没有得到圆满解答。本文以单壁碳纳米管为模型疏水体系,采用分子动力学模拟研究尿素/氧化三甲胺混合溶液中纳米管内部水合性质,结果表明氧化三甲胺更易与水分子和尿素分子形成较强相互作用从而稳定了水溶液结构,这一结果亦表明了氧化三甲胺可以通过间接机制抵消尿素分子对于碳纳米管内部水合性质的影响。     

Mapping the hydration dynamics of ubiquitin

The nature of water's interaction with biomolecules such as proteins has been difficult to examine in detail at atomic resolution. Solution NMR spectroscopy is potentially a powerful method for characterizing both the structural and temporal aspects of protein hydration but has been plagued by artifacts. Encapsulation of the protein of interest within the aqueous core of a reverse micelle particle results in a general slowing of water dynamics, significant reduction in hydrogen exchange chemistry and elimination of contributions from bulk water thereby enabling the use of nuclear Overhauser effects to quantify interactions between the protein surface and hydration water. Here we extend this approach to allow use of dipolar interactions between hydration water and hydrogens bonded to protein carbon atoms. By manipulating the molecular reorientation time of the reverse micelle particle through use of low viscosity liquid propane, the T(1ρ) relaxation time constants of (1)H bonded to (13)C were sufficiently lengthened to allow high quality rotating frame nuclear Overhauser effects to be obtained. These data supplement previous results obtained from dipolar interactions between the protein and hydrogens bonded to nitrogen and in aggregate cover the majority of the molecular surface of the protein. A wide range of hydration dynamics is observed. Clustering of hydration dynamics on the molecular surface is also seen. Regions of long-lived hydration water correspond with regions of the protein that participate in molecular recognition of binding partners suggesting that the contribution of the solvent entropy to the entropy of binding has been maximized through evolution.     

Origin of short-range repulsion between hydrated phospholipid bilayers: a computer simulation study

The grand canonical Monte Carlo technique is used to simulate the pressure-distance dependence for supported dilauroylphosphatidylethanolamine (DLPE) membranes. The intra- and intermolecular interactions in the system are described with a combination of an AMBER-based force field for DLPE and a TIP4P model for water. To improve the balance between the pair interactions of like and unlike molecules, the water-lipid interaction potentials are scaled to reproduce the hydration level and intermembrane separation at full hydration. It is found that the short-range water-mediated repulsion originates from the hydration component of the intermembrane pressure, whereas the direct interaction between the membranes remains attractive throughout the pressure range studied (0-5 kbar).     

A method for calculating the Gibbs energies of hydrophobic effects and specific interactions of nonelectrolytes in aqueous solutions

The thermodynamic characteristics of hydrophobic hydration, the Gibbs energies of hydrophobic effect, were calculated. The method for calculations was based on the division of the Gibbs energy of hydration into contributions of nonspecific interactions, specific interactions between solutes and solvents (if they exist), and hydrophobic effect. In the absence of specific interactions between solutes and water, the Gibbs energy of hydrophobic effect depended linearly on the characteristic molecular volume of the solute for substances with different structures and properties. The universality of this dependence allows the suggestion to be made that it remains valid also in the presence of specific interactions. This allows the Gibbs energy of specific interactions in water to be determined for a wide range of compounds, in particular, for aliphatic alcohols.     

Charge separation in water molecule clusters under thermal fluctuations: 1. Intermolecular interactions

A detailed model of intermolecular interactions in water molecule clusters is developed that makes it possible to describe their disintegration to ions under conditions of finite temperatures by the stochastic simulation methods. In this model, the Hamiltonian in explicit form includes Coulomb, dispersion, exchange, and polarization interactions; many-particle covalent interactions and hydrogen bonds; the interaction of induced dipoles; charge transfers from ions to molecules; and the recombination of counterion charges, as well as the effect of an ion field on the unpaired interactions of molecules. The model is consistent with experimental data on the free energy and entropy of ion hydration in water vapors and the free energy of the hydration of a recombined ion pair.     

Formation and interaction of hydrated alkali metal ions at the graphite-water interface

Ion hydration at a solid surface ubiquitously exists in nature and plays important roles in many natural processes and technological applications. Aiming at obtaining a microscopic insight into the formation of such systems and interactions therein, we have investigated the hydration of alkali metal ions at a prototype surface-graphite (0001), using first-principles molecular dynamics simulations. At low water coverage, the alkali metal ions form two-dimensional hydration shells accommodating at most four (Li, Na) and three (K, Rb, Cs) waters in the first shell. These two-dimensional shells generally evolve into three-dimensional structures at higher water coverage, due to the competition between hydration and ion-surface interactions. Exceptionally K was found to reside at the graphite-water interface for water coverages up to bulk water limit, where it forms an "umbrellalike" surface hydration shell with an average water-ion-surface angle of 115 degrees . Interactions between the hydrated K and Na ions at the interface have also been studied. Water molecules seem to mediate an effective ion-ion interaction, which favors the aggregation of Na ions but prevents nucleation of K. These results agree with experimental observations in electron energy loss spectroscopy, desorption spectroscopy, and work function measurement. In addition, the sensitive dependence of charge transfer on dynamical structure evolution during the hydration process, implies the necessity to describe surface ion hydration from electronic structure calculations.     

Impacts of medium,substituents, and specific interactions with water on hydration of carbonyl compounds

Theoretical study on hydration of carbonyl compounds has been done at B3LYP/6-31++G** and MP2/6-31++G** levels. The variations in ΔG _hyd and hydration constants are explored in terms of medium effect, substituent effect, and hydrogen bonding abilities of carbonyl compounds and their hydrated products. The dielectric of medium decreases the ΔG _hyd values thereby favoring the process. The presence of electron-releasing substituents at the carbonyl carbon disfavors the hydration process, while that of electron-withdrawing substituents favor the process. Hydrogen bonding interactions stabilize the product to a larger extent than the carbonyl molecules, thereby favoring the hydration process. Linear correlation between the calculated log K _hyd values and the experimental values is seen in case of specific interactions with water (R = 0.976) than in the case without those interactions (R = 0.955).     

The hydration properties of carboxybetaine zwitterion brushes

Combined quantum mechanical calculations and classical molecular dynamics simulations were conducted to investigate the hydration properties of carboxybetaine zwitterion brushes with varying separation distances between the quaternary ammonium cation and carboxylic anion. The brushes consist of zwitterion trimers and are investigated to mimic interacting zwitterion chains grafted on a substrate as well as polymers with interacting zwitterion side chains. Our results show that the values of both positive and negative charges, their separation distances as well as chain interactions appear to play a critical role in the hydration properties of the zwitterions. The overall hydration property of these zwitterions is dictated by the competition between the strong hydration of the charged groups and the dehydration of the hydrocarbon chains. The strongest hydration occurs when the ? CH₂? unit in the hydrocarbon chain reaches 6–8 for these trimers. Further increase in the hydrocarbon chain length to 10–14 leads to significant and sudden dehydration of the trimers. The water structure and the water residence time surrounding the zwitterions also demonstrate substantial alteration at this length scale. This hydrophilic‐to‐hydrophobic transition is induced by the hydrophobic interactions of the trimer chains. Our hydration results could explain the observed trend of the superiority of the methylated carbohydrates and poly(ethylene glycol) as antifouling materials compared to corresponding hydroxyl‐terminated compounds. © 2015 Wiley Periodicals, Inc.     

On the molecular origins of volumetric data

We use a statistical thermodynamic approach and a simple thermodynamic model of hydration to examine the molecular origins of the volumetric properties of solutes. In this model, solute-solvent interactions are treated as a binding reaction. The free energy of hydration of the noninteracting solute species coincides with the free energy of cavity formation, while the free energy of solute-solvent interactions is given by the binding polynomial. By differentiating the relationship for the free energy of hydration with respect to temperature and pressure, one obtains the complete set of equations describing the thermodynamic profile of hydration, including enthalpy, entropy, volume, compressibility, expansibility, and so forth. The model enables one to rigorously define in thermodynamic terms the hydration number and the related concept of hydration shell, which are both widely used as operational definitions in experimental studies. Hydration number, nh, is the effective number of water molecules solvating the solute and represents the derivative of the free energy of hydration with respect to the logarithm of water activity. One traditional way of studying hydration relies on the use of volumetric measurements. However, microscopic interpretation of macroscopic volumetric data is complicated and currently relies on empirical models that are not backed by theory. We use our derived model to link the microscopic determinants of the volumetric properties of a solute and its statistical thermodynamic parameters. In this treatment, the partial molar volume, V degrees, of a solute depends on the cavity volume, hydration number, and the properties of waters of hydration. In contrast, the partial molar isothermal compressibility, K degrees T, and expansibility, E degrees, observables, in addition to the intrinsic compressibility and expansibility of the cavity enclosing the solute, hydration number, and the properties of waters of hydration, contain previously unappreciated relaxation terms that originate from pressure- and temperature-induced perturbation of the equilibrium between the solvated solute species. If significant, the relaxation terms may bring about a new level of nonadditivity to compressibility and expansibility group contributions that goes beyond the overlap of the hydration shells of adjacent groups. We apply our theoretical results to numerical analyses of the volume and compressibility responses to changes in the distribution of solvated species of polar compounds.     

Solute�CSolvent Interactions in Aqueous Glycylglycine�CCuCl2 Solutions: Acoustical and Molecular Dynamics Perspective

Acoustical and molecular dynamics studies were carried out to understand the various interactions present in glycylglycine?CCuCl₂ aqueous solutions. Amongst these interactions, hydrogen bonding and solute?Csolvent interactions have been highlighted in this study. The radial distribution function (RDF) was used to investigate solution structure and hydration parameters. Binding of Cu²⁺ with various polar peptide atoms reveals the nature and degree of binding. The formation of complex clusters between glycylglycine and water molecules increases the relaxation time. The first hydration shell considerably influences the structure of the second shell, facilitating the formation of an ordered hydrogen bonded network. Both experimental and theoretical results have proved to be efficient in analyzing the behavior of molecules and to give a clear idea on molecular interactions in solutions.     

How does the solvent guide the hydration process?

The hydration numbers are investigated of the glycine amino acid in solutions of substances with different effects on the structure of water: urea, monomethylurea, and 1,3-dimethylurea. Glycine loses a half of its hydration water in a 20m urea solution and only a quarter of it in a 20m dimethylurea solution. The constancy of the hydration number of glycine in concentrated dimethylurea solutions is due to the compensatory effect of the interactions in the ternary and binary systems.     

Tryptophan–water interaction in Monellin: Hydration patterns from molecular dynamics simulation

Femtosecond spectroscopy carried out earlier on Monellin and some other systems has given insights into the hydration dynamics of the proteins. In the present work, molecular dynamics simulations have been performed on Monellin to study the hydration dynamics. A method has been described to follow up the molecular events of the protein–water interactions in detail. The time constants of the survival correlation function match well with the reported experimental values. This validates the procedure, adapted here for Monellin, to investigate the hydration dynamics in general.