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.

Entropy Tug-of-War Determines Solvent Effects in the Liquid-Liquid Phase Separation of a Globular Protein.

Liquid-liquid phase separation (LLPS) plays a key role in the compartmentalization of cells via the formation of biomolecular condensates. Here, we combined atomistic molecular dynamics (MD) simulations and terahertz (THz) spectroscopy to determine the solvent entropy contribution to the formation of condensates of the human eye lens protein γD-Crystallin. The MD simulations reveal an entropy tug-of-war between water molecules that are released from the protein droplets and those that are retained within the condensates, two categories of water molecules that were also assigned spectroscopically. A recently developed THz-calorimetry method enables quantitative comparison of the experimental and computational entropy changes of the released water molecules. The strong correlation mutually validates the two approaches and opens the way to a detailed atomic-level understanding of the different driving forces underlying the LLPS.

Tuning Acid-Base Chemistry at an Electrified Gold/Water Interface.

Acid-base reactions are ubiquitous in solution chemistry, as well as in electrochemistry. However, macroscopic concepts derived in solutions, such as pKa and pH, differ significantly at electrified metal-aqueous interfaces due to specific solvation and applied voltage. Here, we measure the pKa values of an amino acid, glycine, at a gold/water interface under a varying applied voltage by means of spectroscopic titration. With the help of simulations, we propose a general model to understand potential-dependent shifts in pKa values in terms of local hydrophobicity and electric fields. These parameters can be tuned by adjusting the metal surface and applied voltage, respectively, offering promising, but still unexplored, paths to regulate reactivity. Our results change the focus with respect to common interpretations based on, for example, apparent local pH effects and open interesting perspectives for electrochemical reaction steering.

S.O.S: Shape, orientation, and size tune solvation in electrocatalysis.

Current models to understand the reactivity of metal/aqueous interfaces in electrochemistry, e.g., volcano plots, are based on the adsorption free energies of reactants and products, which are often small hydrophobic molecules (such as in CO2 and N2 reduction). Calculations played a major role in the quantification and comprehension of these free energies in terms of the interactions that the reactive species form with the surface. However, solvation free energies also come into play in two ways: (i) by modulating the adsorption free energy together with solute-surface interactions, as the solute has to penetrate the water adlayer in contact with the surface and get partially desolvated (which costs free energy); (ii) by regulating transport across the interface, i.e., the free energy profile from the bulk to the interface, which is strongly non-monotonic due to the unique nature of metal/aqueous interfaces. Here, we use constant potential molecular dynamics to study the solvation contributions, and we uncover huge effects of the shape and orientation (on top of the already known size effect) of small hydrophobic and amphiphilic solutes on their adsorption free energy. We propose a minimal theoretical model, the S.O.S. model, that accounts for size, orientation, and shape effects. These novel aspects are rationalized by recasting the concepts at the base of the Lum-Chandler-Weeks theory of hydrophobic solvation (for small solutes in the so-called volume-dominated regime) into a layer-by-layer form, where the properties of each interfacial region close to the metal are explicitly taken into account.

Hydration water drives the self-assembly of guanosine monophosphate.

Guanosine monophosphate (GMP) is a nucleotide that can self-assemble in aqueous solution under certain conditions. An understanding of the process at the molecular level is an essential step to comprehend the involvement of DNA substructures in transcription and replication, as well as their relationship to genetic diseases such as cancer. We present the temperature-dependent terahertz (1.5-12 THz, 50-400 cm-1) absorptivity spectra of aqueous Na2 GMP solution in comparison with the aqueous solutions of other RNA nucleotides. Distinct absorption features were observed in the spectrum of GMP, which we attribute to the intramolecular modes of the self-assemblies (i.e., G-complexes) that, at 1 M, start to form at 313 K and below. Changes in broad-band features of the terahertz spectrum were also observed, which we associate with the release of hydration water in the temperature-dependent formation of guanine quadruplexes. Using a state-of-the-art THz calorimetry approach correlating spectroscopic to thermodynamic changes, we propose a molecular mechanism of hydrophilic hydration driving GMP self-assembly as a function of temperature. The free energy contribution of hydrophilic hydration is shown as a decisive factor in guanine-quadruplex formation. Our findings spotlight the role of hydration in the formation of macromolecular structures and suggest the potential of hydration tuning for regulating DNA transcription and replication.

From Local Hydration Motifs in Aqueous Acetic Acid Solutions to Macroscopic Mixing Thermodynamics: A Quantitative Connection from THz-Calorimetry.

We report the results of THz measurements (30-440 cm-1) of aqueous acetic acid solutions over the full mixing range (XAcAc = 0-1). We recorded spectroscopic observables as a function of the acetic acid concentration in the frequency range of the intermolecular stretch at 150 cm-1 and of the librational modes at 350-440 cm-1. This allowed us to unravel changes in hydrophobic and hydrophilic hydration motifs, respectively. By means of a novel THz-calorimetry approach, we quantitatively correlated these changes in local hydration motifs to excess mixing entropy and enthalpy. We find that ΔHmix is determined by both hydrophobic and hydrophilic solvation contributions. In contrast, ΔSmix is governed by hydrophobic cavity formation. Our results further suggest that acetic acid-water mixtures are systems at the edge of phase separation due to endothermic contributions from both hydrophilic and hydrophobic solvation in a large portion of the miscibility range. Our work establishes a quantitative relationship between the balance of local hydrophobic and hydrophilic solvation motifs and the macroscopic mixing thermodynamic properties.

Hydration makes a difference! How to tune protein complexes between liquid-liquid and liquid-solid phase separation.

Understanding how protein rich condensates formed upon liquid-liquid phase separation (LLPS) evolve into solid aggregates is of fundamental importance for several medical applications, since these are suspected to be hot-spots for many neurotoxic diseases. This requires developing experimental approaches to observe in real-time both LLPS and liquid-solid phase separation (LSPS), and to unravel the delicate balance of protein and water interactions dictating the free energy differences between the two. We present a vibrational THz spectroscopy approach that allows doing so from the point of view of hydration water. We focus on a cellular prion protein of high medical relevance, which we can drive to undergo either LLPS or LSPS with few mutations. We find that it is a subtle balance of hydrophobic and hydrophilic solvation contributions that allows tuning between LLPS and LSPS. Hydrophobic hydration provides an entropic driving force to phase separation, through the release of hydration water into the bulk. Water hydrating hydrophilic groups provides an enthalpic driving force to keep the condensates in a liquid state. As a result, when we modify the protein by a few mutations to be less hydrophilic, we shift from LLPS to LSPS. This molecular understanding paves the way for a rational design of proteins.

Local solvation structures govern the mixing thermodynamics of glycerol-water solutions.

Glycerol is a major cryoprotective agent and is widely used to promote protein stabilization. By a combined experimental and theoretical study, we show that global thermodynamic mixing properties of glycerol and water are dictated by local solvation motifs. We identify three hydration water populations, i.e., bulk water, bound water (water hydrogen bonded to the hydrophilic groups of glycerol) and cavity wrap water (water hydrating the hydrophobic moieties). Here, we show that for glycerol experimental observables in the THz regime allow quantification of the abundance of bound water and its partial contribution to the mixing thermodynamics. Specifically, we uncover a 1 : 1 connection between the population of bound waters and the mixing enthalpy, which is further corroborated by the simulation results. Therefore, the changes in global thermodynamic quantity - mixing enthalpy - are rationalized at the molecular level in terms of changes in the local hydrophilic hydration population as a function of glycerol mole fraction in the full miscibility range. This offers opportunities to rationally design polyol water, as well as other aqueous mixtures to optimize technological applications by tuning mixing enthalpy and entropy based on spectroscopic screening.

Real-time measure of solvation free energy changes upon liquid-liquid phase separation of α-elastin.

Biological condensates are known to retain a large fraction of water to remain in a liquid and reversible state. Local solvation contributions from water hydrating hydrophilic and hydrophobic protein surfaces were proposed to play a prominent role for the formation of condensates through liquid-liquid phase separation (LLPS). However, although the total free energy is accessible by calorimetry, the partial solvent contributions to the free energy changes upon LLPS remained experimentally inaccessible so far. Here, we show that the recently developed THz calorimetry approach allows to quantify local hydration enthalpy and entropy changes upon LLPS of α-elastin in real time, directly from experimental THz spectroscopy data. We find that hydrophobic solvation dominates the entropic solvation term, whereas hydrophilic solvation mainly contributes to the enthalpy. Both terms are in the order of hundreds of kJ/mol, which is more than one order of magnitude larger than the total free energy changes at play during LLPS. However, since we show that entropy/enthalpy mostly compensates, a small entropy/enthalpy imbalance is sufficient to tune LLPS. Theoretically, a balance was proposed before. Here we present experimental evidence based on our spectroscopic approach. We finally show that LLPS can be steered by inducing small changes of solvation entropy/enthalpy compensation via concentration or temperature in α-elastin.

Liquid-Liquid Phase Separation? Ask the Water!

Water is more than an inert spectator during liquid-liquid phase separation (LLPS), the reversible compartmentalization of protein solutions into a protein-rich and a dilute phase. We show that LLPS is driven by changes in hydration entropy and enthalpy. Tuning LLPS by adjusting experimental parameters, e.g., addition of co-solutes, is a major goal for biological and medical applications. This requires a general model to quantify thermodynamic driving forces. Here, we develop such a model based on the measured amplitudes of characteristic THz-features of two hydration populations: "Cavity-wrap" water hydrating hydrophobic patches is released during LLPS leading to an increase in entropy. "Bound" water hydrating hydrophilic patches is retained since it is enthalpically favorable. We introduce a THz-phase diagram mapping these spectroscopic/thermodynamic changes. This provides not only a precise understanding of hydrophobic and hydrophilic hydration driving forces as a function of temperature and concentration but also a rational means to tune LLPS.

On the Trail of Molecular Hydrophilicity and Hydrophobicity at Aqueous Interfaces.

Uncovering microscopic hydrophilicity and hydrophobicity at heterogeneous aqueous interfaces is essential as it dictates physico/chemical properties such as wetting, the electrical double layer, and reactivity. Several molecular and spectroscopic descriptors were proposed, but a major limitation is the lack of connections between them. Here, we combine density functional theory-based MD simulations (DFT-MD) and SFG spectroscopy to explore how interfacial water responds in contact with self-assembled monolayers (SAM) of tunable hydrophilicity. We introduce a microscopic metric to track the transition from hydrophobic to hydrophilic interfaces. This metric combines the H/V descriptor, a structural descriptor based on the preferential orientation within the water network in the topmost binding interfacial layer (BIL) and spectroscopic fingerprints of H-bonded and dangling OH groups of water carried by BIL-resolved SFG spectra. This metric builds a bridge between molecular descriptors of hydrophilicity/hydrophobicity and spectroscopically measured quantities and provides a recipe to quantitatively or qualitatively interpret experimental SFG signals.

Spectroscopic Fingerprints of Cavity Formation and Solute Insertion as a Measure of Hydration Entropic Loss and Enthalpic Gain.

Hydration free energies are dictated by a subtle balance of hydrophobic and hydrophilic interactions. We present here a spectroscopic approach, which gives direct access to the two main contributions: Using THz-spectroscopy to probe the frequency range of the intermolecular stretch (150-200 cm-1 ) and the hindered rotations (450-600 cm-1 ), the local contributions due to cavity formation and hydrophilic interactions can be traced back. We show that via THz calorimetry these fingerprints can be correlated 1 : 1 with the group specific solvation entropy and enthalpy. This allows to deduce separately the hydrophobic (i.e. cavity formation) and hydrophilic contributions to thermodynamics, as shown for hydrated alcohols as a case study. Accompanying molecular dynamics simulations quantitatively support our experimental results. In the future our approach will allow to dissect hydration contributions in inhomogeneous mixtures and under non-equilibrium conditions.

Unexpectedly High Capacitance of the Metal Nanoparticle/Water Interface: Molecular-Level Insights into the Electrical Double Layer.

The electrical double-layer plays a key role in important interfacial electrochemical processes from catalysis to energy storage and corrosion. Therefore, understanding its structure is crucial for the progress of sustainable technologies. We extract new physico-chemical information on the capacitance and structure of the electrical double-layer of platinum and gold nanoparticles at the molecular level, employing single nanoparticle electrochemistry. The charge storage ability of the solid/liquid interface is larger by one order-of-magnitude than predicted by the traditional mean-field models of the double-layer such as the Gouy-Chapman-Stern model. Performing molecular dynamics simulations, we investigate the possible relationship between the measured high capacitance and adsorption strength of the water adlayer formed at the metal surface. These insights may launch the active tuning of solid-solvent and solvent-solvent interactions as an innovative design strategy to transform energy technologies towards superior performance and sustainability.

The role of hydrophobic hydration in the free energy of chemical reactions at the gold/water interface: Size and position effects.

Metal/water interfaces catalyze a large variety of chemical reactions, which often involve small hydrophobic molecules. In the present theoretical study, we show that hydrophobic hydration at the Au(100)/water interface actively contributes to the reaction free energy by up to several hundreds of meV. This occurs either in adsorption/desorption reaction steps, where the vertical distance from the surface changes in going from reactants to products, or in addition and elimination reaction steps, where two small reactants merge into a larger product and vice versa. We find that size and position effects cannot be captured by treating them as independent variables. Instead, their simultaneous evaluation allows us to map the important contributions, and we provide examples of their combinations for which interfacial reactions can be either favored or disfavored. By taking a N2 and a CO2 reduction pathway as test cases, we show that explicitly considering hydrophobic effects is important for the selectivity and rate of these relevant interfacial processes.

Stripping away ion hydration shells in electrical double-layer formation: Water networks matter.

The double layer at the solid/electrolyte interface is a key concept in electrochemistry. Here, we present an experimental study combined with simulations, which provides a molecular picture of the double-layer formation under applied voltage. By THz spectroscopy we are able to follow the stripping away of the cation/anion hydration shells for an NaCl electrolyte at the Au surface when decreasing/increasing the bias potential. While Na+ is attracted toward the electrode at the smallest applied negative potentials, stripping of the Cl- hydration shell is observed only at higher potential values. These phenomena are directly measured by THz spectroscopy with ultrabright synchrotron light as a source and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations.

Cation enrichment in the ion atmosphere is promoted by local hydration of DNA.

Electrostatic interactions are central to the structure and function of nucleic acids, including their folding, condensation, and interaction with proteins and other charged molecules. These interactions are profoundly affected by ions surrounding nucleic acids, the constituents of the so-called ion atmosphere. Here, we report precise Fourier Transform-Terahertz/Far-Infrared (FT-THz/FIR) measurements in the frequency range 30-500 cm-1 for a 24-bp DNA solvated in a series of alkali halide (NaCl, NaF, KCl, CsCl, and CsF) electrolyte solutions which are sensitive to changes in the ion atmosphere. Cation excess in the ion atmosphere is detected experimentally by observation of cation modes of Na+, K+, and Cs+ in the frequency range between 70-90 cm-1. Based on MD simulations, we propose that the magnitude of cation excess (which is salt specific) depends on the ability of the electrolyte to perturb the water network at the DNA interface: In the NaF atmosphere, the ions reduce the strength of interactions between water and the DNA more than in case of a NaCl electrolyte. Here, we explicitly take into account the solvent contribution to the chemical potential in the ion atmosphere: A decrease in the number of bound water molecules in the hydration layer of DNA is correlated with enhanced density fluctuations, which decrease the free energy cost of ion-hydration, thus promoting further ion accumulation within the DNA atmosphere. We propose that taking into account the local solvation is crucial for understanding the ion atmosphere.

Local Mutations Can Serve as a Game Changer for Global Protein Solvent Interaction.

Although it is well-known that limited local mutations of enzymes, such as matrix metalloproteinases (MMPs), may change enzyme activity by orders of magnitude as well as its stability, the completely rational design of proteins is still challenging. These local changes alter the electrostatic potential and thus local electrostatic fields, which impacts the dynamics of water molecules close the protein surface. Here we show by a combined computational design, experimental, and molecular dynamics (MD) study that local mutations have not only a local but also a global effect on the solvent: In the specific case of the matrix metalloprotease MMP14, we found that the nature of local mutations, coupled with surface morphology, have the ability to influence large patches of the water hydrogen-bonding network at the protein surface, which is correlated with stability. The solvent contribution can be experimentally probed via terahertz (THz) spectroscopy, thus opening the door to the exciting perspective of rational protein design in which a systematic tuning of hydration water properties allows manipulation of protein stability and enzymatic activity.

Molecular Fingerprints of Hydrophobicity at Aqueous Interfaces from Theory and Vibrational Spectroscopies.

Hydrophobicity/hydrophilicity of aqueous interfaces at the molecular level results from a subtle balance in the water-water and water-surface interactions. This is characterized here via density functional theory-molecular dynamics (DFT-MD) coupled with vibrational sum frequency generation (SFG) and THz-IR absorption spectroscopies. We show that water at the interface with a series of weakly interacting materials is organized into a two-dimensional hydrogen-bonded network (2D-HB-network), which is also found above some macroscopically hydrophilic silica and alumina surfaces. These results are rationalized through a descriptor that measures the number of "vertical" and "horizontal" hydrogen bonds formed by interfacial water, quantifying the competition between water-surface and water-water interactions. The 2D-HB-network is directly revealed by THz-IR absorption spectroscopy, while the competition of water-water and water-surface interactions is quantified from SFG markers. The combination of SFG and THz-IR spectroscopies is thus found to be a compelling tool to characterize the finest details of molecular hydrophobicity at aqueous interfaces.

Size dependence of hydrophobic hydration at electrified gold/water interfaces.

Hydrophobic hydration at metal/water interfaces actively contributes to the energetics of electrochemical reactions, e.g. [Formula: see text] and [Formula: see text] reduction, where small hydrophobic molecules are involved. In this work, constant applied potential molecular dynamics is employed to study hydrophobic hydration at a gold/water interface. We propose an adaptation of the Lum-Chandler-Weeks (LCW) theory to describe the free energy of hydrophobic hydration at the interface as a function of solute size and applied voltage. Based on this model we are able to predict the free energy cost of cavity formation at the interface directly from the free energy cost in the bulk plus an interface-dependent correction term. The interfacial water network contributes significantly to the free energy, yielding a preference for outer-sphere adsorption at the gold surface for ideal hydrophobes. We predict an accumulation of small hydrophobic solutes of sizes comparable to CO or [Formula: see text], while the free energy cost to hydrate larger hydrophobes, above 2.5-Å radius, is shown to be greater at the interface than in the bulk. Interestingly, the transition from the volume dominated to the surface dominated regimes predicted by the LCW theory in the bulk is also found to take place for hydrophobes at the Au/water interface but occurs at smaller cavity radii. By applying the adapted LCW theory to a simple model addition reaction, we illustrate some implications of our findings for electrochemical reactions.

The key role of solvent in condensation: Mapping water in liquid-liquid phase-separated FUS.

Formation of biomolecular condensates through liquid-liquid phase separation (LLPS) has emerged as a pervasive principle in cell biology, allowing compartmentalization and spatiotemporal regulation of dynamic cellular processes. Proteins that form condensates under physiological conditions often contain intrinsically disordered regions with low-complexity domains. Among them, the RNA-binding proteins FUS and TDP-43 have been a focus of intense investigation because aberrant condensation and aggregation of these proteins is linked to neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal dementia. LLPS occurs when protein-rich condensates form surrounded by a dilute aqueous solution. LLPS is per se entropically unfavorable. Energetically favorable multivalent protein-protein interactions are one important aspect to offset entropic costs. Another proposed aspect is the release of entropically unfavorable preordered hydration water into the bulk. We used attenuated total reflection spectroscopy in the terahertz frequency range to characterize the changes in the hydrogen bonding network accompanying the FUS enrichment in liquid-liquid phase-separated droplets to provide experimental evidence for the key role of the solvent as a thermodynamic driving force. The FUS concentration inside LLPS droplets was determined to be increased to 2.0 mM independent of the initial protein concentration (5 or 10 µM solutions) by fluorescence measurements. With terahertz spectroscopy, we revealed a dewetting of hydrophobic side chains in phase-separated FUS. Thus, the release of entropically unfavorable water populations into the bulk goes hand in hand with enthalpically favorable protein-protein interaction. Both changes are energetically favorable, and our study shows that both contribute to the thermodynamic driving force in phase separation.