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

Solvation free energy in governing equations for DNA hybridization, protein-ligand binding, and protein folding

This work examines the thermodynamics of model biomolecular interactions using a governing equation that accounts for the participation of bulk water in the equilibria. In the first example, the binding affinities of two DNA duplexes, one of nine and one of 10 base pairs in length, are measured and characterized by isothermal titration calorimetry (ITC) as a function of concentration. The results indicate that the change in solvation free energy that accompanies duplex formation (ΔGS) is large and unfavorable. The duplex with the larger number of G:C pairings yields the largest change in solvation free energy, ΔGS = +460 kcal·mol-1per base pair at 25 °C. A van't Hoff analysis of the data is complicated by the varying degree of intramolecular base stacking within each DNA strand as a function of temperature. A modeling study reveals how the solvation free energy alters the output of a typical ITC experiment and leads to a good, though misleading, fit to the classical equilibrium equation. The same thermodynamic framework is applied to a model protein-ligand interaction, the binding of ribonuclease A with the nucleotide inhibitor 3'-UMP, and to a conformational equilibrium, the change in tertiary structure of α-lactalbumin in molar guanidinium chloride solutions. The ribonuclease study yields a value of ΔGS = +160 kcal·mol-1, whereas the folding equilibrium yields ΔGS ≈ 0, an apparent characteristic of hydrophobic interactions. These examples provide insight on the role of solvation energy in binding equilibria and suggest a pivot in the fundamental application of thermodynamics to solution chemistry.

Unraveling the hydration dynamics of ACC1-13K24 with ATP: From liquid to droplet to amyloid fibril

In order to achieve a comprehensive understanding of protein aggregation processes, an exploration of solvation dynamics, a key yet intricate component of biological phenomena, is mandatory. In the present study, we used Fourier transform infrared spectroscopy and terahertz spectroscopy complemented by atomic force microscopy and kinetic experiments utilizing thioflavin T fluorescence to elucidate the changes in solvation dynamics during liquid-liquid phase separation and subsequent amyloid fibril formation, the latter representing a transition from liquid to solid phase separation. These processes are pivotal in the pathology of neurodegenerative disorders such as Alzheimer's and Parkinson's diseases. We focus on the ACC1-13K24-ATP protein complex, which undergoes fibril formation followed by droplet generation. Our investigation reveals the importance of hydration as a driving force in these processes, offering new insights into the molecular mechanisms at play.

Structural Order in the Hydration Shell of Nonpolar Groups versus that in Bulk Water

The poor solubility of nonpolar compounds in water around room temperature is governed by a large and negative entropy change, whose molecular cause is still debated. Since the Frank and Evans original proposal in 1945, the large and negative entropy change is usually attributed to the formation of ordered structures in the hydration shell of nonpolar groups. However, the existence of such ordered structures has never been proven. The present study is aimed at providing available structural results and thermodynamic arguments disproving the existence of ordered structures in the hydration shell of nonpolar groups.

Raman spectroscopy in the study of amyloid formation and phase separation

Neurodegenerative diseases, such as Alzheimer's and Parkinson's, share a common pathological feature of amyloid structure accumulation. However, the structure-function relationship between these well-ordered, β-sheet-rich, filamentous protein deposits and disease etiology remains to be defined. Recently, an emerging hypothesis has linked phase separation, a process involved in the formation of protein condensates, to amyloid formation, suggesting that liquid protein droplets serve as loci for amyloid initiation. To elucidate how these processes contribute to disease progression, tools that can directly report on protein secondary structural changes are needed. Here, we review recent studies that have demonstrated Raman spectroscopy as a powerful vibrational technique for interrogating amyloid structures; one that offers sensitivity from the global secondary structural level to specific residues. This probe-free technique is further enhanced via coupling to a microscope, which affords structural data with spatial resolution, known as Raman spectral imaging (RSI). In vitro and in cellulo applications of RSI are discussed, highlighting studies of protein droplet aging, cellular internalization of fibrils, and Raman imaging of intracellular water. Collectively, utilization of the myriad Raman spectroscopic methods will contribute to a deeper understanding of protein conformational dynamics in the complex cellular milieu and offer potential clinical diagnostic capabilities for protein misfolding and aggregation processes in disease states.

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.

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.

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.

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.

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.

Studying Local Electrostatics by Terahertz Spectroscopy Using Amines as a Probe

In a previous study[1] we could show that a large amplitude mode of the zwitterion glycine can serve as a sensitive probe for protonation and allows to deduce local pKa values. Here we show that the underlying concept is more general: We present the results of a pH dependent measurement of Terahertz-FTIR (THz-FTIR) spectra of solvated amines, i. e. Diethylamine (DEA), Triethylamine (TEA), and Diisopropylamine (DiPA). We show that amines serve as a sensitive, label free probe for local protonation. Protonation of the amines yield intensity changes which can be quantified by precise THz spectroscopy (30 cm-1 -450 cm-1 ). A detailed analysis allows us to correlate the titration spectra of solvated amines in the THz range with pKa values. This demonstrates the potential of THz spectroscopy to probe the charge state of biomolecules in water in a label free manner.

The explicit role of interfacial hydration during polyethylene glycol induced lipid fusion: a THz spectroscopic investigation

While the phenomenon of excipient mediated membrane fusion has been studied widely, the inherent role of interfacial hydration involved in the process has mostly remained unaddressed. Here we report the experimental validation of the fact that PEG-induced membrane fusion is associated with the dehydration of the membrane(s). We explore the explicit hydration behavior at three different lipids (DOPC, POPC and DPPC) membranes with different aliphatic tails as they undergo fusogenic transition in the presence of PEG of average molecular weight of 4000 using THz-FTIR spectroscopy in the frequency window of 1.5-13.5 THz. Dynamic light scattering and electron microscopic measurements confirm the formation of different intermediate steps of the liposomes during the fusion process: bilayer aggregation, destabilization and finally lipid fusion. We observe that membrane hydration follows a systematic trend with the lipid specificity as the fusion process sets in.

Energy Dissipation into the Solvent during Proton Transfer Occurs via Acoustic Phonons

In photochemistry, rapid energy dissipation into the solvent is mandatory to prevent radiation damages. By optical pump THz spectroscopy, we are able to follow the details of the energy dissipation mechanism upon photoexcitation of the photoacid to the hydrogen-bonded network of water: Based on the frequency maps subsequent to photoexcitation, we propose that energy transfer takes place via propagation of an acoustic phonon. The dissipation into the solvent can be rationalized by an initial first hydration shell response followed by energy dissipation via an acoustic phonon. Surprisingly, for the first 10 ps, the propagation in the water network can be described by a wave packet with a constant group velocity, indicating a long-range correlation. After 300 ps, thermalization in the liquid jet is reached and the THz spectrum reflects a Boltzmann population, corresponding a temperature increase of ΔT = 0.5 °C.

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.

Impact of hydrophobicity on local solvation structures and its connection with the global solubilization thermodynamics of amphiphilic molecules

The relationship between the local solvation structures and global thermodynamics, specifically in the case of amphiphilic molecules, is a complex phenomenon and is not yet fully understood. With the prior knowledge that local solvation structures can impose a significant impact on the overall solvation process, we here combine THz spectroscopic analysis with MD simulations to investigate the impact of the altered hydrophobicity and polarity of amphiphilic solute molecules on the local solvation configurations. We use two water soluble alcohols: ethanol (EtOH) and its fluorinated counterpart, 2,2,2-trifluoroethanol (TFE), as model solutes. Our study is aimed to determine the relative abundance of different hydrogen bonded conformers and to establish a correlation between the spectral signatures (as obtained from THz spectroscopic measurements) and microscopic solute-solvent interactions associated with the local solvation structures (as obtained from MD simulations). Finally, we estimate the possible energetic parameters associated with the alcohol solubilization process. We found that while both the alcohols are completely water soluble, they receive a contrasting solvation energy share in terms of entropy and enthalpy. We understand that these findings are not limited to the specific system studied here but can be broadly extrapolated to other amphiphilic aqueous solutions.

Solvation Plays a Key Role in Antioxidant-Mediated Attenuation of Elevated Creatinine Level: An In Vitro Spectroscopic Investigation

An elevated level of creatinine (CRN) is a mark of kidney ailment, and prolonged retention of such condition could lead to renal failure, associated with severe ischemia. Antioxidants are clinically known to excrete CRN from the body through urine, thereby reducing its level in blood. The molecular mechanism of such an exclusion process is still illusive. As the excretion channel is urine, solvation of the solute is expected to play a pivotal role. Here, we report a detailed time-domain and frequency-domain terahertz (THz) spectroscopic investigation to understand the solvation of CRN in the presence of two model antioxidants, mostly used to treat elevated CRN level: N-Acetyl-l-cysteine (NAC) and ascorbic acid (ASC). FTIR spectroscopy in the mid-infrared region and UV absorption spectroscopy measurements coupled with quantum chemical calculations [at the B3LYP/6-311G++(d,p) level] reveal that both NAC and ASC form HBonded complexes with CRN and rapidly undergo a barrier-less proton transfer process to form creatinium ions. THz measurements provide explicit evidence of the formation of highly solvated complexes compared with bare CRN, which eventually enables its excretion through urine. These observations could provide a foundation for designing more beneficial drugs to resolve kidney diseases..

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.

Caught in the act: real-time observation of the solvent response that promotes excited-state proton transfer in pyranine

Photo-induced excited-state proton transfer (ESPT) reactions are of central importance in many biological and chemical processes. Identifying mechanistic details of the solvent reorganizations that facilitate proton transfer however, is challenging for current experimental and theoretical approaches. Using optical pump THz probe (OPTP) spectroscopy and molecular dynamics simulations, we were able to elucidate the ultrafast changes in the solvation environment for three derivatives of pyranine: the photoacid HPTS, the methoxy derivative MPTS, and the photobase OPTS. Experimentally, we find damped oscillations in the THz signal at short times and our simulations enable their assignment to vibrational energy transfer beatings between the photoexcited chromophore and nearby solvent molecules. The simulations of HPTS reveal strikingly efficient sub-ps energy transfer into a particular solvent mode, that is active near 4 THz, and which can provide the requisite energy required for solvent reorganization promoting proton transfer. Similar oscillations are present in the THz signal for all three derivatives, however the signal is damped rapidly for HPTS (within 0.4 ps) and more slowly for MPTS (within 1.4 ps) and OPTS (within 2.0 ps). For HPTS, we also characterize an additional phonon-like propagation of the proton into the bulk with a 140 ps period and an 83 ps damping time. Thermalization of the solvent occurs on a time scale exceeding 120 ps.

The birth and evolution of solvated electrons in the water

The photo-induced radiolysis of water is an elementary reaction in biology and chemistry, forming solvated electrons, OH radicals, and hydronium cations on fast time scales. Here, we use an optical-pump terahertz-probe spectroscopy setup to trigger the photoionization of water molecules with optical laser pulses at ~400 nm and then time-resolve the transient solvent response with broadband terahertz (THz) fields with a ~90 fs time resolution. We observe three distinct spectral responses. The first is a positive broadband mode that can be attributed to an initial diffuse, delocalized electron with a radius of (22 ± 1) Å, which is short lived (<200 fs) because the absorption is blue-shifting outside of the THz range. The second emerging spectroscopic signature with a lifetime of about 150 ps is attributed to an intermolecular mode associated with a mass rearrangement of solvent molecules due to charge separation of radicals and hydronium cations. After 0.2 ps, we observe a long-lasting THz signature with depleted intensity at 110 cm^-1 that is well reproduced by ab initio molecular dynamics. We interpret this negative band at 110 cm^-1 as the solvent cage characterized by a weakening of the hydrogen bond network in the first and second hydration shells of the cavity occupied by the localized electron.

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.

Electronic properties of amino acids and nucleobases: similarity classes and pairing principles from chemical reactivity indices

We present a new classification scheme for amino acids and nucleobases based on the electronic properties of the individual molecules. Using chemical reactivity indices such as electronegativity, electrophilicity, and chemical hardness, we can identify similarities and differences between each class of amino acids and nucleobases. Notable differences emerge in particular with regard to high, neutral or low electronegativity as well as different combinations of chemical hardness. Our approach allows us to relate these insights to the properties of the side groups in terms of a unique reference scheme. We further show that hydrophobic differences between amino acids are rather negligible in the context of electronic properties. Our classification scheme also rationalizes the occurrence of distinct stable nucleobase pairs and clearly emphasizes certain differences between individual molecules. The stability and abundant occurrence of Watson-Crick nucleobase pairs is further discussed in the context of the minimum electrophilicity principle.

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.