Solvation free energies and hydration structure of N-methyl-p-nitroaniline

Solvation Gibbs energies of N-methyl-p-nitroaniline (MNA) in water and 1-octanol are calculated using the expanded ensemble molecular dynamics method with a force field taken from the literature. The accuracy of the free energy calculations is verified with the experimental Gibbs free energy data and found to reproduce the experimental 1-octanol∕water partition coefficient to within ±0.1 in log unit. To investigate the hydration structure around N-methyl-p-nitroaniline, an independent NVT molecular dynamics simulation was performed at ambient conditions. The local organization of water molecules around the solute MNA molecule was investigated using the radial distribution function (RDF), the coordination number, and the extent of hydrogen bonding. The spatial distribution functions (SDFs) show that the water molecules are distributed above and below the nitrogen atoms parallel to the plane of aromatic ring for both the methylamino and nitro functional groups. It is found that these groups have a significant effect on the hydration of MNA with water molecules forming two weak hydrogen bonds with both the methylamino and nitro groups. The hydration structures around the functional groups in MNA in water are different from those that have been found for methylamine, nitrobenzene, and benzene in aqueous solutions, and these differences together with weak hydrogen bonds explain the lower solubility of MNA in water. The RDFs together with SDFs provide a tool for the understanding the hydration of MNA (and other molecules) and therefore their solubility.     

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: a detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.     

Hydrophobic Effects in Adsorptive Protein Immobilization

In an aqueous environment, noncovalent polar bonding (i.e., through ionic, charge transfer, or hydrogen bonds) may be assumed to be relatively weak because of thestrong chargesolvating and hydrogen-bonding ability of water. It has been suggested, the refore, that apolar (hydrophobic) forces may be the most important single factor providing the driving force for noncovalent intermolecular interactions in an aqueous milieu [1, 2]. However, in aqueous solutions electrostatic interaction may become important when the pertaining charges are shielded by accompanying hydrophobic groups from the quenching effect of water [1]. In this connection it may be pointed out that, on account of the low dielectric constant of the hydrophobic interior of a protein, intramolecular interaction of oppositely charged surface groups can be strong indeed [3].     

Effect of Protic Ionic Liquid and Surfactant Structure on Partitioning of Polyoxyethylene Non‐ionic Surfactants

The partitioning constants and Gibbs free energies of transfer of poly(oxyethylene) n‐alkyl ethers between dodecane and the protic ionic liquids (ILs) ethylammonium nitrate (EAN) and propylammonium nitrate (PAN) are determined. EAN and PAN have a sponge‐like nanostructure that consists of interpenetrating charged and apolar domains. This study reveals that the ILs solvate the hydrophobic and hydrophilic parts of the amphiphiles differently. The ethoxy groups are dissolved in the polar region of both ILs by means of hydrogen bonds. The environment is remarkably water‐like and, as in water, the solubility of the ethoxy groups in EAN decreases on warming, which underscores the critical role of the IL hydrogen‐bond network for solubility. In contrast, amphiphile alkyl chains are not preferentially solvated by the charged or uncharged regions of the ILs. Rather, they experience an average IL composition and, as a result, partitioning from dodecane into the IL increases as the cation alkyl chain is lengthened from ethyl to propyl, because the IL apolar volume fraction increases. Together, these results show that surfactant dissolution in ILs is related to structural compatibility between the head or tail group and the IL nanostructure. Thus, these partitioning studies reveal parameters for the effective molecular design of surfactants in ILs.     

Hydrogen bonding and the cryoprotective properties of glycerol/water mixtures

Molecular dynamics simulations and infrared spectroscopy were used to determine the hydrogen bond patterns of glycerol and its mixtures with water. The ability of glycerol/water mixtures to inhibit ice crystallization is linked to the concentration of glycerol and the hydrogen bonding patterns formed by these solutions. At low glycerol concentrations, sufficient amounts of bulk-like water exist, and at low temperature, these solutions demonstrate crystallization. As the glycerol concentration is increased, the bulk-like water pool is eventually depleted. Water in the first hydration shell becomes concentrated around the polar groups of glycerol, and the alkyl groups of glycerol self-associate. Glycerol-glycerol hydrogen bonds become the dominant interaction in the first hydration shell, and the percolation nature of the water network is disturbed. At glycerol concentrations beyond this point, glycerol/water mixtures remain glassy at low temperatures and the glycerol-water hydrogen bond favors a more linear arrangement. High glycerol concentration mixtures mimic the strong hydrogen bonding pattern seen in ice, yet crystallization does not occur. Hydrogen bond patterns are discussed in terms of hydrogen bond angle distributions and average hydrogen bond number. Shift in infrared frequency of related stretch and bend modes is also reviewed.     

The molecular mechanism of stabilization of proteins by TMAO and its ability to counteract the effects of urea

Trimethylamine n-oxide (TMAO) is a naturally occurring osmolyte that stabilizes proteins and offsets the destabilizing effects of urea. To investigate the molecular mechanism of these effects, we have studied the thermodynamics of interaction between TMAO and protein functional groups. The solubilities of a homologous series of cyclic dipeptides were measured by differential refractive index and the dissolution heats were determined calorimetrically as a function of TMAO concentration at 25 degrees C. The transfer free energy of the amide unit (-CONH-) from water to 1 M TMAO is large and positive, indicating an unfavorable interaction between the TMAO solution and the amide unit. This unfavorable interaction is enthalpic in origin. The interaction between TMAO and apolar groups is slightly favorable. The transfer free energy of apolar groups from water to TMAO consists of favorable enthalpic and unfavorable entropic contributions. This is in contrast to the contributions for the interaction between urea and apolar groups. Molecular dynamics simulations were performed to provide a structural framework for the interpretation of these results. The simulations show enhancement of water structure by TMAO in the form of a slight increase in the number of hydrogen bonds per water molecule, stronger water hydrogen bonds, and long-range spatial ordering of the solvent. These findings suggest that TMAO stabilizes proteins via enhancement of water structure, such that interactions with the amide unit are discouraged.     

Partial molar enthalpies of solution as indicators of interactions in mixtures of 2-butoxyethanol and 2-butanol with water

Calorimetric measurements have been made of differential enthalpies of solution of both components in the binary system 2-butoxyethanol-water and of 2-butanol in the system 2-butanol-water as a function of composition at three different temperatures. The heat capacity changes for dissolution were calculated from the temperature variation of the solution enthalpies. Drastic changes of the solution properties are seen with increasing solute concentration in water-rich solutions. In the 2-butoxyethanol-water system, which could be studied over the whole composition range, four different regions can be identified. At extreme dilution in water, the solute is fully hydrated with a primary hydration layer of monolayer thickness involved in long-range secondary hydration. In dilute solutions the primary hydration layer is unchanged but the secondary hydration diminishes with increasing solute concentration. In semi-dilute solution the primary hydration layer breaks down and the particular hydrophobic characteristics of hydrocarbon groups in aqueous solution disappear. At higher solute content the mixtures show no hydrophobic character but the behavior of regular mixtures of polar solutes.     

Estimation of the hydration of polar groups of α-amino acids by differential scanning calorimetry

The heat capacity of hydration of zwitterions derived from aliphatic amino acids depends linearly on the surface area of the amino acid side radicals accessible to water molecules with the slopeb = 2.35±0.11 J mol^–1 K^–1 Å^–2 at 298 K. The linear correlation between hydration heat capacities of zwitterions of aliphatic amino acids and the corresponding aliphatic alcohols with a coefficient of approximately unity confirms the assumption that hydrophobic hydration does not depend on the nature of the surrounding groups. Using the assumption that the hydration of hydrocarbon radicals is independent of the neighboring groups, theb value has been used to calculate the contributions of polar groups. The contributions of OH, COON, and CONH groups of the side radicals in polar amino acids in the zwitterion form are close to zero; in the case of organic nonionic molecules, these contributions are negative. The increments for polar groups obtained for the zwitterions can be used for the calculation of the heat capacities of proteins and polypeptides incorporating charged amino acid residues. The difference between hydrophilic and hydrophobic hydration mechanisms is manifested not only as different magnitudes and signs of heat capacities and temperature coefficients but also in the fact that the neighboring polar (charged) groups have an effect on hydrophilic hydration but have no effect on hydrophobic hydration.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 2237–2242, September, 1996.     

Interaction of urea with amino acids: implications for urea-induced protein denaturation

The molecular mechanism of urea-induced protein denaturation is not yet fully understood. Mainly two opposing mechanisms are controversially discussed, according to which either hydrophobic, or polar interactions are the dominant driving force. To resolve this question, we have investigated the interactions between urea and all 20 amino acids by comprehensive molecular dynamics simulations of 22 tripeptides. Calculation of atomic contact frequencies between the amino acids and solvent molecules revealed a clear profile of solvation preferences by either water or urea. Almost all amino acids showed preference for contacts with urea molecules, whereas charged and polar amino acids were found to have slight preferences for contact with water molecules. Particularly strong preference for contacts to urea were seen for aromatic and apolar side-chains, as well as for the protein backbone of all amino acids. Further, protein-urea hydrogen bonds were found to be significantly weaker than protein-water or water-water hydrogen bonds. Our results suggest that hydrophobic interactions are the dominant driving force, while hydrogen bonds between urea and the protein backbone contribute markedly to the overall energetics by avoiding unfavorable unsatisfied hydrogen bond sites on the backbone. In summary, we suggest a combined mechanism that unifies the two current and seemingly opposing views.     

Hydration heat capacities of hydrophobic solutes at 248–373 K

A new equation is suggested to define the temperature dependence of the Gibbs energy of hydration of hydrophobic substances: ΔG ⁰ = b ₀ + b ₁ T + b ₂lnT. According to this equation, the hydration heat capacity is in inverse proportion to temperature. Consistent values of hydration heat capacity of nonpolar solutes have been obtained for different temperatures using data on solubility and dissolution enthalpy. The contributions of the hydrocarbon radicals and OH group to the heat capacity of hydration of the compounds were found for the temperature range 248–373 K. The hydration heat capacity of the hydroxyl group has a weak dependence on temperature and increases by only 12 J/(mol·K) in the specified temperature interval. Changes in the hydration entropy of hydrophobic and OH groups are calculated for the temperature increasing from 248 K to 373 K.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

Local heterogeneous dynamics of water around lysozyme: a computer simulation study

Water present near the surface of a protein exhibits dynamic properties different from that of water in the pure bulk state. In this work, we have carried out atomistic molecular dynamics simulation of an aqueous solution of hen egg-white lysozyme. Attempts have been made to explore the correlation between the local heterogeneous mobility of water around the protein segments and the rigidity of the hydration layers with the microscopic dynamics of hydrogen bonds formed by water molecules with the protein residues. The kinetics of breaking and reformation of hydrogen bonds involving the surface water molecules have been calculated. It is found that the reformations of broken hydrogen bonds are more frequent for the hydration layers of those segments of the protein that are more rigid. The calculation of the low-frequency vibrational modes of hydration layer water molecules reveals that the protein influences the transverse and longitudinal degrees of freedom of water around it in a differential manner. These findings can be verified by appropriate experimental studies.     

Microheterogeneity in aqueous-organic solutions: Heat capacities,volumes and expansibilities of some alcohols,aminoalcohol and tertiary amines in water

The heat capacities per unit volume and the densities of aqueous solutions of 2-propanol, neopentanol, tert-amylalcohol, 2-amino-2-methylpropanol, triethylamine and diethylmethylamine were measured, in many cases as a function of temperature, over the whole mole fraction or solubility range. Apparent and partial molal heat capacities, volumes and expansibilities were derived. The concentration dependence of these functions suggest the existence of transitions in some of these systems, in the water-rich region, qualitatively similar to micellization. The large relaxation contribution observed with some of the thermodynamic functions of hydrophobic alcohols and amines suggests a reinforcement of hydrophobic hydration due to strong hydrogen-bonding interactions of the polar groups with water.     

Temperature dependence of dimerization and dewetting of large-scale hydrophobes: a molecular dynamics study

We studied by molecular dynamics simulations the temperature dependence of hydrophobic association and drying transition of large-scale solutes. Similar to the behavior of small solutes, we found the association process to be characterized by a large negative heat capacity change. The origin of this large change in heat capacity is the high fragility of hydrogen bonds between water molecules at the interface with hydrophobic solutes; an increase in temperature breaks more hydrogen bonds at the interface than in the bulk. With increasing temperature, both entropy and enthalpy changes for association strongly decrease, while the change in free energy weakly varies, exhibiting a small minimum at high temperatures. At around T=Ts=360 K, the change in entropy is zero, a behavior similar to the solvation of small nonpolar solutes. Unexpectedly, we find that at Ts, there is still a substantial orientational ordering of the interfacial water molecules relative to the bulk. Nevertheless, at this point, the change in entropy vanishes due to a compensating contribution of translational entropy. Thus, at Ts, there is rotational order and translational disorder of the interfacial water relative to bulk water. In addition, we studied the temperature dependence of the drying-wetting transition. By calculating the contact angle of water on the hydrophobic surface at different temperatures, we compared the critical distance observed in the simulations with the critical distance predicted by macroscopic theory. Although the deviations of the predicted from the observed values are very small (8-23%), there seems to be an increase in the deviations with an increase in temperature. We suggest that these deviations emerge due to increased fluctuations, characterizing finite systems, as the temperature increases.     

Influence of graft interface polarity on hydration/dehydration of grafted thermoresponsive polymer brushes and steroid separation using all-aqueous chromatography

We have prepared poly( N-isopropylacrylamide) (PIPAAm) brush-grafted surfaces with varied temperature-responsive hydrophobic properties through surface-initiated atom transfer radical polymerization (ATRP). These temperature-responsive surfaces were characterized by chromatographic analysis using modified silica beads as a chromatographic stationary phase in aqueous mobile phase. Mixed silane self-assembled monolayers (SAMs) comprising ATRP initiator and silanes with various terminal functional groups were formed on the silica bead surfaces. IPAAm was then polymerized by ATRP using the CuCl/CuCl2/Me6TREN catalyst system in 2-propanol at 25 degrees C for 16 h. The chromatographic retention behavior of steroids on the resulting PIPAAm brushes made on more polar silane components was distinct from that on more apolar silane interfaces. Retention times for steroids on PIPAAm mixed apolar silane graft interfaces were significantly longer than those on analogous polar silane interfaces due to enhanced dehydration of PIPAAm brushes on apolar silane-grafted surfaces. Changes in retention factor, k', on polar silane PIPAAm-grafted interfaces were relatively large compared to that on apolar PIPAAm grafted interfaces due to larger hydration/dehydration alterations of grafted PIPAAm brushes on the former surfaces. Applied step-temperature gradients from 50 to 10 degrees C show that PIPAAm brushes on polar silane interfaces tend to hydrate more, leading to shorter retention times. In conclusion, the polarity of the grafted interface significantly influences the grafted PIPAAm brush hydration/dehydration characteristics and subsequently also the temperature-modulated separation of bioactive compounds in all-aqueous chromatography.     

Water's hydrogen bonds in the hydrophobic effect: a simple model

We propose a simple analytical model to account for water's hydrogen bonds in the hydrophobic effect. It is based on computing a mean-field partition function for a water molecule in the first solvation shell around a solute molecule. The model treats the orientational restrictions from hydrogen bonding, and utilizes quantities that can be obtained from bulk water simulations. We illustrate the principles in a 2-dimensional Mercedes-Benz-like model. Our model gives good predictions for the heat capacity of hydrophobic solvation, reproduces the solvation energies and entropies at different temperatures with only one fitting parameter, and accounts for the solute size dependence of the hydrophobic effect. Our model supports the view that water's hydrogen bonding propensity determines the temperature dependence of the hydrophobic effect. It explains the puzzling experimental observation that dissolving a nonpolar solute in hot water has positive entropy.     

Hydrophobic forces and hydrogen bonds in the adhesion between retinoid-coated surfaces

Interactions between hydrophobic chains of lipid monolayers and interactions between hydrophilic headgroups of lipid bilayers (with or without a molecular recognition step) are now well documented, especially for commonly used lipids. Here, we report force measurements between a new class of fluorinated lipid layers whose headgroups (synthetic ligands of retinoid receptors) display a very unusual polar/apolar character and can interact via a combination of hydrophobic forces and hydrogen bonds. Although these two interactions produce adhesion and are therefore not easily distinguishable, we show that it is possible to extract both contributions unambiguously. Experiments are performed both in pure water, where the adhesion is a combination of hydrophobic forces and hydrogen bonds, and in Tris buffer, where the hydrophobic effect is the dominant short-range attractive force. The contribution of hydrophobic forces scaled down to molecular interactions is deduced from force versus distance profiles, and the same value is found independently in pure water and Tris buffer, about 1 kBT. We also show that retinoid lipid layers attract each other through a very long-range (100 nm) exponential force, which is insensitive to the pH and the salinity. The origin of this long-range attraction is discussed on the basis of previously proposed mechanisms.     

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.