Temperature dependence of three-body hydrophobic interactions: potential of mean force, enthalpy, entropy, heat capacity, and nonadditivity

Temperature-dependent three-body hydrophobic interactions are investigated by extensive constant-pressure simulations of methane-like nonpolar solutes in TIP4P model water at six temperatures. A multiple-body hydrophobic interaction is considered to be (i) additive, (ii) cooperative, or (iii) anti-cooperative if its potential of mean force (PMF) is (i) equal to, (ii) smaller than, or (iii) larger than the corresponding pairwise sum of two-methane PMFs. We found that three-methane hydrophobic interactions at the desolvation barrier are anti-cooperative at low to intermediate T, and vary from essentially additive to slightly cooperative at high T. Interactions at the contact minimum are slightly anti-cooperative over a wider temperature range. Enthalpy, entropy, and heat capacity are estimated from the computed PMFs. Contrary to the common expectation that burial of solvent-accessible nonpolar surface area always leads to a decrease in heat capacity, the present results show that the change in heat capacity upon three-methane association is significantly positive at the desolvation barrier and slightly positive at the contact minimum. This suggests that the heat capacity signature of a hydrophobic polymer need not vary uniformly nor monotonically with conformational compactness. Ramifications for protein folding are discussed.     

Global thermodynamics of hydrophobic cavitation, dewetting, and hydration

Pure water experimental and simulation results are combined to predict the thermodynamics of cavity formation, spanning atomic to macroscopic length scales, over the entire ambient liquid temperature range. The resulting cavity equation of state is used to quantify dewetting excess contributions to cavity formation thermodynamics and construct a thermodynamic perturbation theory of hydrophobic hydration. Predictions are compared with large cavity simulations and experimental rare-gas hydration thermodynamics data (for He, Ne, Ar, Kr, Xe, and Rn). Key findings include the strong temperature dependence of the critical length scale for hydrophobic dewetting and the evaluation of fundamental solute-solvent interaction contributions to rare-gas hydration chemical potentials.     

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.     

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.     

Molecular dynamics simulations of pressure effects on hydrophobic interactions.

We report results on the pressure effects on hydrophobic interactions obtained from molecular dynamics simulations of aqueous solutions of methanes in water. A wide range of pressures that is relevant to pressure denaturation of proteins is investigated. The characteristic features of water-mediated interactions between hydrophobic solutes are found to be pressure-dependent. In particular, with increasing pressure we find that (1) the solvent-separated configurations in the solute-solute potential of mean force (PMF) are stabilized with respect to the contact configurations; (2) the desolvation barrier increases monotonically with respect to both contact and solvent-separated configurations; (3) the locations of the minima and the barrier move toward shorter separations; and (4) pressure effects are considerably amplified for larger hydrophobic solutes. Together, these observations lend strong support to the picture of the pressure denaturation process proposed previously by Hummer et al. (Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 1552): with increasing pressure, the transfer of water into protein interior becomes key to the pressure denaturation process, leading to the dissociation of close hydrophobic contacts and subsequent swelling of the hydrophobic protein interior through insertions of water molecules. The pressure dependence of the PMF between larger hydrophobic solutes shows that pressure effects on the interaction between hydrophobic amino acids may be considerably amplified compared to those on the methane-methane PMF.     

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.     

Estimating the temperature dependence of peptide folding entropies and free enthalpies from total energies in molecular dynamics simulations

The temperature dependence of thermodynamic quantities, such as heat capacity, entropy and free enthalpy, may be obtained by using well-known equations that relate these quantities to the enthalpy of the molecular system of interest at a range of temperatures. In turn, the enthalpy of a molecular system can be estimated from molecular dynamics simulations of an appropriate model. To demonstrate this, we have investigated the temperature dependence of the enthalpy, heat capacity, entropy and free enthalpy of a system that consists of a beta-heptapeptide in methanol and have used the statistical mechanics relationships to describe the thermodynamics of the folding/unfolding equilibrium of the peptide. The results illustrate the power of current molecular simulation force fields and techniques in establishing the link between thermodynamic quantities and conformational distributions.     

Heats of solution of several freons in water from 5 to 45°C

An extensive set of measurements of the heats of solution of a number of fluorine-containing gases (CCl₂F₂, CClF₃, CBrF₃, CF₄) have been performed from 5 to 45°C. The temperature dependence allows accurate determination of the heat capacity change in the solution process as a function of temperature. To a first approximation these heat capacity results agree with the predictions of the simple two energy state model for water molecules in the first solvation shell. The further extension of the general applicability of this model, originally developed to account for the thermodynamic properties of solvation of hydrocarbon and other small apolar gases in water, suggests that the unique thermodynamic properties of hydrophobic solvation are largely determined by the water molecules in the first solvation shell.     

Comparison of thermodynamic properties of coarse-grained and atomic-level simulation models.

Thermodynamic data are often used to calibrate or test amomic-level (AL) force fields for molecular dynamics (MD) simulations. In contrast, the majority of coarse-grained (CG) force fields do not rely extensively on thermodynamic quantities. Recently, a CG force field for lipids, hydrocarbons, ions, and water, in which approximately four non-hydrogen atoms are mapped onto one interaction site, has been proposed and applied to study various aspects of lipid systems. To date, no extensive investigation of its capability to describe salvation thermodynamics has been undertaken. In the present study, a detailed picture of vaporization, solvation, and phase-partitioning thermodynamics for liquid hydrocarbons and water was obtained at CG and AL resolutions, in order to compare the two types or models and evaluate their ability to describe thermodynamic properties in the temperature range between 263 and 343 K. Both CG and AL models capture the experimental dependence of the thermodynamic properties on the temperature, albeit a systematically weaker dependence is found for the CG model. Moreover, deviations are found for solvation thermodynamics and for the corresponding enthalpy-entropy compensation for the CG model. Particularly water/oil repulsion seems to be overestimated. However, the results suggest that the thermodynamic properties considered should be reproducible by a CG model provided it is reparametrized on the basis of these liquid-phase properties.     

A statistical model for instable thermodynamical systems

A generic model is presented for statistical systems which display thermodynamic features in contrast to our everyday experience, such as infinite and negative heat capacities. Such system are instable in terms of classical equilibrium thermodynamics. Using our statistical model, we are able to investigate states of instable systems which are undefined in the framework of equilibrium thermodynamics. We show that a region of negative heat capacity in the adiabatic environment, leads to a first order like phase transition when the system is coupled to a heat reservoir. This phase transition takes place without a phase coexistence. Nevertheless, all intermediate states are stable due to fluctuations. When two instable system are brought in thermal contact, the temperature of the composed system is lower than the minimum temperature of the individual systems. Generally, the equilibrium states of instable system cannot be simply decomposed into equilibrium states of the individual systems. The properties of instable system depend on the environment, ensemble equivalence is broken.     

Thermodynamics of Micellization from Heat‐Capacity Measurements

Differential scanning calorimetry (DSC), the most important technique for studying the thermodynamics of structural transitions of biological macromolecules, is seldom used in quantitative thermodynamic studies of surfactant micellization/demicellization. The reason for this could be ascribed to an insufficient understanding of the temperature dependence of the heat capacity of surfactant solutions (DSC data) in terms of thermodynamics, which leads to problems with the design of experiments and interpretation of the output signals. We address these issues by careful design of DSC experiments performed with solutions of ionic and nonionic surfactants at various surfactant concentrations, and individual and global mass‐action model analysis of the obtained DSC data. Our approach leads to reliable thermodynamic parameters of micellization for all types of surfactants, comparable with those obtained by using isothermal titration calorimetry (ITC). In summary, we demonstrate that DSC can be successfully used as an independent method to obtain temperature‐dependent thermodynamic parameters for micellization.     

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.     

Molecular simulations suggest protein salt bridges are uniquely suited to life at high temperatures

A series of explicit solvent molecular dynamics simulations has been performed to investigate the temperature dependence of salt bridge interactions between two freely diffusing amino acids. The simulations, performed at 25, 50, 75, and 100 degrees C, allow a large number of distinct association and dissociation events to be directly observed, without the imposition of additional forces to drive association. Analysis of contact frequencies for atom pairs demonstrates that the number of salt bridge contacts between the two molecules is unaffected by temperature, whereas the numbers of hydrophobic and polar contacts are greatly diminished. A second, independent set of simulations-using rigid, prototypical molecule types-allows the differing temperature dependences of hydrophobic, polar, and salt bridge interactions to be unambiguously examined. In the prototype molecule simulations, the salt bridge interaction is found to substantially increase in stability at 100 degrees C relative to 25 degrees C. This difference in behavior between flexible amino acids and rigid prototype molecules is perhaps a direct manifestation of the effects of conformational entropy on association thermodynamics. Overall, the results demonstrate that salt bridge interactions are extremely resilient to temperature increases and, as such, are uniquely suited to promoting protein stability at high temperatures.     

Gaussian excitations model for glass-former dynamics and thermodynamics

We describe a model for the thermodynamics and dynamics of glass-forming liquids in terms of excitations from an ideal glass state to a Gaussian manifold of configurationally excited states. The quantitative fit of this three parameter model to the experimental data on excess entropy and heat capacity shows that "fragile" behavior, indicated by a sharply rising excess heat capacity as the glass transition is approached from above, occurs in anticipation of a first-order transition--usually hidden below the glass transition--to a "strong" liquid state of low excess entropy. The distinction between fragile and strong behavior of glass formers is traced back to an order of magnitude difference in the Gaussian width of their excitation energies. Simple relations connect the excess heat capacity to the Gaussian width parameter, and the liquid-liquid transition temperature, and strong, testable, predictions concerning the distinct properties of energy landscape for fragile liquids are made. The dynamic model relates relaxation to a hierarchical sequence of excitation events each involving the probability of accumulating sufficient kinetic energy on a separate excitable unit. Super-Arrhenius behavior of the relaxation rates, and the known correlation of kinetic with thermodynamic fragility, both follow from the way the rugged landscape induces fluctuations in the partitioning of energy between vibrational and configurational manifolds. A relation is derived in which the configurational heat capacity, rather than the configurational entropy of the Adam-Gibbs equation, controls the temperature dependence of the relaxation times, and this gives a comparable account of the experimental observations without postulating a divergent length scale. The familiar coincidence of zero mobility and Kauzmann temperatures is obtained as an approximate extrapolation of the theoretical equations. The comparison of the fits to excess thermodynamic properties of laboratory glass formers, and to configurational thermodynamics from simulations, reveals that the major portion of the excitation entropy responsible for fragile behavior resides in the low-frequency vibrational density of states. The thermodynamic transition predicted for fragile liquids emerges from beneath the glass transition in case of laboratory water and the unusual heat capacity behavior observed for this much studied liquid can be closely reproduced by the model.     

Two-Gaussian excitations model for the glass transition

We develop a modified "two-state" model with Gaussian widths for the site energies of both ground and excited states, consistent with expectations for a disordered system. The thermodynamic properties of the system are analyzed in configuration space and found to bridge the gap between simple two-state models ("logarithmic" model in configuration space) and the random energy model ("Gaussian" model in configuration space). The Kauzmann singularity given by the random energy model remains for very fragile liquids but is suppressed or eliminated for stronger liquids. The sharp form of constant-volume heat capacity found by recent simulations for binary mixed Lennard-Jones and soft-sphere systems is reproduced by the model, as is the excess entropy and heat capacity of a variety of laboratory systems, strong and fragile. The ideal glass in all cases has a narrow Gaussian, almost invariant among molecular and atomic glassformers, while the excited-state Gaussian depends on the system and its width plays a role in the thermodynamic fragility. The model predicts the possibility of first-order phase transitions for fragile liquids. The analysis of laboratory data for toluene and o-terphenyl indicates that fragile liquids resolve the Kauzmann paradox by a first-order transition from supercooled liquid to ideal-glass state at a temperature between T(g) and Kauzmann temperature extrapolated from experimental data. We stress the importance of the temperature dependence of the energy landscape, predicted by the fluctuation-dissipation theorem, in analyzing the liquid thermodynamics.     

Salt-dependent heat capacity changes for RNA duplex formation

Two short non-self-complementary RNA duplexes have been probed extensively to assess whether the heat capacity change (DeltaCP) for RNA duplex formation is sensitive to the solution ionic strength. Isothermal titration calorimetry data for duplex formation were collected as a function of temperature to obtain the DeltaCP values under each condition. Ionic strength was varied by repeating the measurements in the presence of 0.1-1.5 M added NaCl. In both cases, the DeltaCP for duplex formation was shown to have a significant dependence on the solution ionic strength, varying linearly with the log of the added NaCl. In each case, the DeltaCP became more negative at high ionic strength. Heat capacity changes have been much better studied in the area of protein folding where DeltaCP is not strongly dependent upon ionic strength but instead depends primarily on the burial of hydrophobic surface. Whereas hydrophobic burial is almost certainly involved in RNA duplex formation as well, the polyanionic nature of the backbone and the associated ion condensation effects are almost certainly a contributor to the observed ionic strength dependence of DeltaCPs for duplex formation. These data also hint at the possibility that nearest neighbor effects may be evident in the heat capacity changes associated with nucleic acid folding.     

Pressure and temperature dependence of hydrophobic hydration: volumetric, compressibility, and thermodynamic signatures

The combined effect of pressure and temperature on hydrophobic hydration of a nonpolar methanelike solute is investigated by extensive simulations in the TIP4P model of water. Using test-particle insertion techniques, free energies of hydration under a range of pressures from 1 to 3000 atm are computed at eight temperatures ranging from 278.15 to 368.15 K. Corresponding enthalpy, entropy, and heat capacity accompanying the hydration process are estimated from the temperature dependence of the free energies. Partial molar and excess volumes calculated using pressure derivatives of the simulated free energies are consistent with those determined by direct volume simulations; but direct volume determination offers more reliable estimates for compressibility. At 298.15 K, partial molar and excess isothermal compressibilities of methane are negative at 1 atm. Partial molar and excess adiabatic (isentropic) compressibilities are estimated to be also negative under the same conditions. But partial molar and excess isothermal compressibilities are positive at high pressures, with a crossover from negative to positive compressibility at approximately 100-1000 atm. This trend is consistent with experiments on aliphatic amino acids and pressure-unfolded states of proteins. For the range of pressures simulated, hydration heat capacity exhibits little pressure dependence, also in apparent agreement with experiment. When pressure is raised at constant room temperature, hydration free energy increases while its entropic component remains essentially constant. Thus, the increasing unfavorability of hydration under raised pressure is seen as largely an enthalpic effect. Ramifications of the findings of the authors for biopolymer conformational transitions are discussed.     

Estimation of the heat capacity of individual substances on the basis of experimental enthalpy increments

A method of determination of the high-temperature (T > 298.15 K) heat capacity function of individual substances by differentiation of enthalpy increment measurements with respect to temperature has been suggested. The approach is based on approximation of enthalpy increment measurements by a polynomial without free term, which makes it possible to strictly meet the boundary conditions of differentiation determined by the properties of the enthalpy increment function. In combination with the Shomate method, the approach enables the optimal choice of the form of temperature dependence of heat capacity that is consistent best of all with the concepts of classical thermodynamics and provides the best reproducibility of enthalpy increment measurements. When low-temperature heat capacity measurements for individual substances are not available, i.e., the Shomate method is inapplicable, the suggested strategy makes it possible to obtain the heat capacity change function in the temperature range of enthalpy increment measurements and to estimate the heat capacity of the substance at standard temperature.