Temperature dependence of hydrophobic hydration dynamics: from retardation to acceleration

Bridging the Gap in Cryopreservation Mechanism: Unraveling the Interplay between Structure, Dynamics, and Thermodynamics in Cryoprotectant Aqueous Solutions

Cryoprotectants play a crucial role in preserving biological material, ensuring their viability during storage and facilitating crucial applications such as the conservation of medical compounds, tissues, and organs for transplantation. However, the precise mechanism by which cryoprotectants modulate the thermodynamic properties of water to impede the formation and growth of ice crystals, thus preventing long-term damage, remains elusive. This is evident in the use of empirically optimized recipes for mixtures that typically contain DMSO, glycerol, and various sugar constituents. Here, we use terahertz calorimetry, Overhauser nuclear polarization, and molecular dynamics simulations to show that DMSO exhibits a robust structuring effect on water around its methyl groups, reaching a maximum at a DMSO mole fraction of XDMSO = 0.33. In contrast, glycerol exerts a smaller water-structuring effect, even at higher concentrations (Scheme 1). These results potentially suggest that the wrapped water around DMSO's methyl group, which can be evicted upon ligand binding, may render DMSO a more surface-active cryoprotectant than glycerol, while glycerol may participate more as a viscogen that acts on the entire sample. These findings shed light on the molecular intricacies of cryoprotectant solvation behavior and have potentially significant implications for optimizing cryopreservation protocols.

The influence of force fields on the structure and dynamics of water confined in ZIF-8 from atomistic simulations

The complexity of modeling flexible crystals, such as ZIF-8, mainly stems from the handling of intramolecular interactions. Numerous force fields have been proposed in the literature to describe the interactions between atoms in ZIF-8. We employ seven force fields to examine the structure and dynamic behavior of water molecules confined in ZIF-8, with the aim of investigating the impact of force fields on simulation results. Various structural characterization methods consistently indicate that the choice of different force fields has quantitative effects but no qualitative effects on the structural characteristics of confined water. Additionally, the force fields do not impact the qualitative description of the diffusion mechanism. Both mean-square displacement and van Hove autocorrelation function reveal two characteristic movements of water molecules diffusing in ZIF-8: a short-time intra-cavity hopping process and a long-time inter-cavity hopping process. However, the framework flexibility is found to play a crucial role in determining the order of spatial arrangement and local structure, self-diffusion coefficient and reorientational dynamics of confined water. Specifically, the DREIDING force field gives rise to an unrealistic stiff framework, enhancing the order of spatial arrangement and diminishing the local ordered structure of confined water. Meanwhile, it results in much slower translational and reorientational dynamics. Hence, the general DREIDING force field cannot be considered for providing a quantitative description of the water structure and dynamics.

Effects of polarizability and charge transfer on water dynamics and the underlying activation energies

A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy-entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics.

Structure and dynamics of a water/methanol mixture confined in zeolitic imidazolate framework ZIF-8 from atomistic simulations

A classical atomistic simulation study is reported for the microscopic structure and dynamics of a water/methanol mixture confined in flexible nanoporous zeolitic imidazolate framework ZIF-8. Both the radial density distribution and vivid two-dimensional density profile demonstrate that methanol molecules can roughly be viewed as "embedded" between two layers of water molecules to form a "sandwich" structure. The reason for the formation of such a specific structure is explained based on the hydrogen-bonding state and the strength of various hydrogen bonds. The investigation of guest molecular diffusion shows that the self-diffusion coefficient of confined water is generally one to two orders of magnitude smaller than that of bulk water. In addition, the dependence of the self-diffusion coefficient on loading is non-monotonic: the self-diffusion coefficient firstly shows a significant increase and then decreases at higher loading. Moreover, both the structure and dynamics of the hydrogen bond (HB) network of confined water molecules are investigated in a spatially resolved manner. The results indicate that both the HB structure and dynamics of water molecules near the ZIF-8 surface deviate significantly from those of bulk water. However, while water molecules located at the pore center are relatively similar to bulk water molecules with respect to the HB structure, they exhibit strong slowdown in HB dynamics when compared with bulk water. This simulation study elucidates in detail the structural and dynamical properties of a water/methanol mixture in nanoscopic ZIF-8 confinement, which is expected to provide a deep insight into the role of porous fillers, such as ZIF-8, in improving the performance of the dehydration of alcohols via pervaporation and other related processes.

Hydration and Seamless Integration of Hydrogen Peroxide in Water

Raman multivariate curve resolution is used to decompose the vibrational spectra of aqueous hydrogen peroxide (H2O2) into pure water, dilute H2O2, and concentrated H2O2 spectral components. The dilute spectra reveal four sub-bands in the OH stretch region, assigned to the OH stretch and Fermi resonant bend overtone of H2O2, and two nonequivalent OH groups on water molecules that donate a hydrogen bond to H2O2. At high concentrations, a spectral component resembling pure H2O2 emerges. Our results further demonstrate that H2O2 perturbs the structure of water significantly less than either methanol or sodium chloride of the same concentration, as evidenced by comparing the hydration-shell spectra of tert-butyl alcohol dissolved in the three aqueous solutions.

Revealing the Frank-Evans "Iceberg" Structures within the Solvation Layer around Hydrophobic Solutes

Using computer simulations, the structural properties of solvation water of three model hydrophobic molecules, methane and two fullerenes (C60 and C80), were studied. Systems were simulated at temperatures in the range of 250-298 K. By analyzing both the local ordering of the molecules of water in the solvation layers and the structure of hydrogen bond network, it is shown that in the solvation layer of hydrophobic molecules, ordered aggregates consisting of water molecules are formed. Even though it is difficult to define the exact structure of these aggregates, their existence alone is clearly noticeable. Moreover, these aggregates become more pronounced with the decrease of temperature. The existence of the ordered aggregates around the hydrophobic solutes complies with the concept of "icebergs" proposed by Frank and Evans.

Influence of crowding on hydrophobic hydration-shell structure

The influence of molecular crowding on water structure, and the associated crossover behavior, is quantified using Raman multivariate curve resolution (Raman-MCR) hydration-shell vibrational spectroscopy of aqueous tert-butyl alcohol, 2-butyl alcohol and 2-butoxyethanol solutions of variable concentration and temperature. Changes in the hydration-shell OH stretch band shape and mean frequency are used to identify the temperature at which the hydration-shell crosses over from a more ordered to less ordered structure, relative to pure water. The influence of crowding on the crossover is found to depend on solute size and shape in a way that is correlated with the corresponding infinitely dilute hydration-shell structure (and the corresponding first hydration-shell spectra are invariably very similar to pure water). Analysis of the results using a Muller-like two-state equilibrium between more ordered and less ordered hydration-shell structures implies that crossover temperature changes are dictated primarily by enthalpic stabilization of the more ordered hydration-shell structures.

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.

CO2 Diffusion in Various Carbonated Beverages: A Molecular Dynamics Study

Carbonated beverages are widely enjoyed in spare time, yet there remain many physical and chemical processes clouded at the molecular level. In this report, we employ molecular dynamics simulations to estimate the diffusion coefficients of CO₂ and the molecular origin of its variations in three model systems with characteristic features of champagnes, sugar-based cola drinks, and club sodas. The computed diffusion coefficients of CO₂ are in good agreement with experimental data. Analyses of hydrogen bonding and the solvent's structural and dynamic properties reveal that the change of CO₂ diffusion coefficient is closely associated with the diffusional behavior of the solvent water itself, as a result of changes in the number and strength of hydrogen bonding interactions among the species and solvent.

Aqueous TMAO solutions as seen by theoretical THz spectroscopy: hydrophilic versus hydrophobic water

All THz resonances of aqueous TMAO solutions are computed and assigned based on ab initio molecular dynamics simulations.

Ab Initio Simulations of Water Dynamics in Aqueous TMAO Solutions: Temperature and Concentration Effects

We use ab initio molecular dynamics simulation to study the effect of hydrophobic groups on the dynamics of water molecules in aqueous solutions of trimethylamine N-oxide (TMAO). We show that hydrophobic groups induce a moderate (<2-fold) slowdown of water reorientation and hydrogen-bond dynamics in dilute solutions, but that this slowdown rapidly increases with solute concentration. In addition, the slowdown factor is found to vary very little with temperature, thus suggesting an entropic origin. All of these results are in quantitative agreement with prior classical molecular dynamics simulations and with the previously suggested excluded-volume model. The hydrophilic TMAO headgroup is found to affect water dynamics more strongly than the hydrophobic moiety, and the magnitude of this slowdown is very sensitive to the strength of the water-solute hydrogen-bond.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Hydration structure of osmolyte TMAO: concentration/pressure-induced response

The effect of solute concentration/pressure on the TMAO hydration structure was studied to understand its protective action under abiotic stressors.

Fast and slow dynamics and the local structure of liquid and supercooled water next to a hydrophobic amino acid

We study, through molecular dynamics simulations, the structure and orientational dynamics of water next to a blocked hydrophobic amino acid, valine (Val), above and below the freezing point of water.

Comparative study of hydration shell dynamics around a hyperactive antifreeze protein and around ubiquitin

The hydration layer surrounding a protein plays an essential role in its biochemical function and consists of a heterogeneous ensemble of water molecules with different local environments and different dynamics. What determines the degree of dynamical heterogeneity within the hydration shell and how this changes with temperature remains unclear. Here, we combine molecular dynamics simulations and analytic modeling to study the hydration shell structure and dynamics of a typical globular protein, ubiquitin, and of the spruce budworm hyperactive antifreeze protein over the 230-300 K temperature range. Our results show that the average perturbation induced by both proteins on the reorientation dynamics of water remains moderate and changes weakly with temperature. The dynamical heterogeneity arises mostly from the distribution of protein surface topographies and is little affected by temperature. The ice-binding face of the antifreeze protein induces a short-ranged enhancement of water structure and a greater slowdown of water reorientation dynamics than the non-ice-binding faces whose effect is similar to that of ubiquitin. However, the hydration shell of the ice-binding face remains less tetrahedral than the bulk and is not "ice-like". We finally show that the hydrogen bonds between water and the ice-binding threonine residues are particularly strong due to a steric confinement effect, thereby contributing to the strong binding of the antifreeze protein on ice crystals.

Characterization of the Local Structure in Liquid Water by Various Order Parameters

A wide range of geometric order parameters have been suggested to characterize the local structure of liquid water and its tetrahedral arrangement, but their respective merits have remained elusive. Here, we consider a series of popular order parameters and analyze molecular dynamics simulations of water, in the bulk and in the hydration shell of a hydrophobic solute, at 298 and 260 K. We show that these parameters are weakly correlated and probe different distortions, for example the angular versus radial disorders. We first combine these complementary descriptions to analyze the structural rearrangements leading to the density maximum in liquid water. Our results reveal no sign of a heterogeneous mixture and show that the density maximum arises from the depletion in interstitial water molecules upon cooling. In the hydration shell of the hydrophobic moiety of propanol, the order parameters suggest that the water local structure is similar to that in the bulk, with only a very weak depletion in ordered configurations, thus confirming the absence of any iceberg-type structure. Finally, we show that the main structural fluctuations that affect water reorientation dynamics in the bulk are angular distortions, which we explain by the jump hydrogen-bond exchange mechanism.