Direct Characterization of Hydrophobic Hydration during Cold and Pressure Denaturation

Unfolding under Pressure: An NMR Perspective

This review aims to analyse the role of solution nuclear magnetic resonance spectroscopy in pressure-induced in vitro studies of protein unfolding. Although this transition has been neglected for many years because of technical difficulties, it provides important information about the forces that keep protein structure together. We first analyse what pressure unfolding is, then provide a critical overview of how NMR spectroscopy has contributed to the field and evaluate the observables used in these studies. Finally, we discuss the commonalities and differences between pressure-, cold- and heat-induced unfolding. We conclude that, despite specific peculiarities, in both cold and pressure denaturation the important contribution of the state of hydration of nonpolar side chains is a major factor that determines the pressure dependence of the conformational stability of proteins.

The Cold Thermal Response of an Amyloid Oligomer Differs from Typical Globular Protein Cold Denaturation

In contrast with the general behavior of folded proteins, the cold thermal response of amyloid assemblies is difficult to elicit with simple models. We exploit exhaustive simulations to evaluate the thermal response of a barrel-shaped model amyloid oligomer, with a distinct hydrophobic core akin to that of folded proteins. Cumulative thermal data over the range of 210-483 K indicate a sharp inflection and rise in structural stability as the temperature is decreased below the melting temperature of the water model. This is not commensurate with the equilibrium free energy profile obtained with core packing as the order parameter. However, energetic analyses and the size of their fluctuations indicate the crucial role of hydration in mediating structural transitions, beyond the expected temperature-dependent hydrophobic effect. Structural ordering of the hydration layer over bulk water is maximized at the transition and vanishes at high temperatures. This is a first direct demonstration of the microscopic influence of hydration water on the low-temperature response of an amyloid assembly close to the cryo-regime.

How Hydrophobic Hydration Destabilizes Surfactant Micelles at Low Temperature: A Coarse-Grained Simulation Study

Micelles are self-assembled aggregates of amphiphilic surfactant molecules that are important in a variety of applications, including drug delivery, detergency, and catalysis. It is known that the micellization process is driven by the same physiochemical forces that drive protein folding, aggregation, and biological membrane self-assembly. Nevertheless, the molecular details of how micelle stability changes in water at low temperature are not fully clear. We develop and use a coarse-grained model to investigate how the interplay between nonionic surfactants and the surrounding water at the nanoscale affects the stability of micelles at high and low temperatures. Simulations of preformed C₁₂E₅ micelles in explicit water at a range of temperatures reveal the existence of two distinct surfactant conformations within the micelle, a bent structure and an extended structure, the latter being more prevalent at low temperature. Favorable interactions of the surfactant with more ordered solvation water stabilizes the extended configuration, allowing nanoscale wetting of the dry, hydrophobic core of the micelle, leading to the micelle breaking. Taken together, our coarse-grained simulations unravel how energetic and structural changes of the surfactant and the surrounding water destabilize micelles at low temperature, which is a direct consequence of the weakened hydrophobicity. Our approach thus provides an effective mean for extracting the molecular-level changes during hydrophobicity-driven destabilization of molecular self-assembly, which is important in a wide range of fields, including biology, polymer science, and nanotechnology.

Interplay between Conformational Heterogeneity and Hydration in the Folding Landscape of a Designed Three-Helix Bundle

Water is known to play a critical role in protein folding and stability. Here we develop and employ a coarse-grained (CG) model to directly explore the role of water in shaping the conformational landscape explored during protein folding. For this purpose, we simulate a designed sequence with binary patterning of neutral and hydrophobic residues, which is capable of folding to a three-helix bundle in explicit water. We find two folded states of this sequence, with rotation of the helices occurring to trade between hydrophobic packing and water expulsion from the core. This work provides insight into the role of water and hydrophobicity in generating competing folded states for a protein.

Signatures of protein thermal denaturation and local hydrophobicity in domain specific hydration behavior: a comparative molecular dynamics study

We investigate, using atomistic molecular dynamics simulations, the association of surface hydration accompanying local unfolding in the mesophilic protein Yfh1 under a series of thermal conditions spanning its cold and heat denaturation temperatures. The results are benchmarked against the thermally stable protein, Ubq, and behavior at the maximum stability temperature. Local unfolding in Yfh1, predominantly in the beta sheet regions, is in qualitative agreement with recent solution NMR studies; the corresponding Ubq unfolding is not observed. Interestingly, all domains, except for the beta sheet domains of Yfh1, show increased effective surface hydrophobicity with increase in temperature, as reflected by the density fluctuations of the hydration layer. Velocity autocorrelation functions (VACF) of oxygen atoms of water within the hydration layers and the corresponding vibrational density of states (VDOS) are used to characterize alteration in dynamical behavior accompanying the temperature dependent local unfolding. Enhanced caging effects accompanying transverse oscillations of the water molecules are found to occur with the increase in temperature preferentially for the beta sheet domains of Yfh1. Helical domains of both proteins exhibit similar trends in VDOS with changes in temperature. This work demonstrates the existence of key signatures of the local onset of protein thermal denaturation in solvent dynamical behavior.

Contribution of Water to Pressure and Cold Denaturation of Proteins

The mechanisms of cold and pressure denaturation of proteins are matter of debate and are commonly understood as due to water-mediated interactions. Here, we study several cases of proteins, with or without a unique native state, with or without hydrophilic residues, by means of a coarse-grain protein model in explicit solvent. We show, using Monte Carlo simulations, that taking into account how water at the protein interface changes its hydrogen bond properties and its density fluctuations is enough to predict protein stability regions with elliptic shapes in the temperature-pressure plane, consistent with previous theories. Our results clearly identify the different mechanisms with which water participates to denaturation and open the perspective to develop advanced computational design tools for protein engineering.

The non-uniform early structural response of globular proteins to cold denaturing conditions: a case study with Yfh1

The mechanism of cold denaturation in proteins is often incompletely understood due to limitations in accessing the denatured states at extremely low temperatures. Using atomistic molecular dynamics simulations, we have compared early (nanosecond timescale) structural and solvation properties of yeast frataxin (Yfh1) at its temperature of maximum stability, 292 K (Ts), and the experimentally observed temperature of complete unfolding, 268 K (Tc). Within the simulated timescales, discernible "global" level structural loss at Tc is correlated with a distinct increase in surface hydration. However, the hydration and the unfolding events do not occur uniformly over the entire protein surface, but are sensitive to local structural propensity and hydrophobicity. Calculated infrared absorption spectra in the amide-I region of the whole protein show a distinct red shift at Tc in comparison to Ts. Domain specific calculations of IR spectra indicate that the red shift primarily arises from the beta strands. This is commensurate with a marked increase in solvent accessible surface area per residue for the beta-sheets at Tc. Detailed analyses of structure and dynamics of hydration water around the hydrophobic residues of the beta-sheets show a more bulk water like behavior at Tc due to preferential disruption of the hydrophobic effects around these domains. Our results indicate that in this protein, the surface exposed beta-sheet domains are more susceptible to cold denaturing conditions, in qualitative agreement with solution NMR experimental results.

Observation of solvent penetration during cold denaturation of E. coli phosphofructokinase-2

Phosphofructokinase-2 is a dimeric enzyme that undergoes cold denaturation following a highly cooperative N2 2I mechanism with dimer dissociation and formation of an expanded monomeric intermediate. Here, we use intrinsic fluorescence of a tryptophan located at the dimer interface to show that dimer dissociation occurs slowly, over several hours. We then use hydrogen-deuterium exchange mass spectrometry experiments, performed by taking time points over the cold denaturation process, to measure amide exchange throughout the protein during approach to the cold denatured state. As expected, a peptide corresponding to the dimer interface became more solvent exposed over time at 3°C; unexpectedly, amide exchange increased throughout the protein over time at 3°C. The rate of increase in amide exchange over time at 3°C was the same for each region and equaled the rate of dimer dissociation measured by tryptophan fluorescence, suggesting that dimer dissociation and formation of the cold denatured intermediate occur without appreciable buildup of folded monomer. The observation that throughout the protein amide exchange increases as phosphofructokinase-2 cold denatures provides experimental evidence for theoretical predictions that cold denaturation primarily occurs by solvent penetration into the hydrophobic core of proteins in a sequence-independent manner.

Single-molecule spectroscopy of cold denaturation and the temperature-induced collapse of unfolded proteins

Recent Förster resonance energy transfer (FRET) experiments show that heat-unfolded states of proteins become more compact with increasing temperature. At the same time, NMR results indicate that cold-denatured proteins are more expanded than heat-denatured proteins. To clarify the connection between these observations, we investigated the unfolded state of yeast frataxin, whose cold denaturation occurs at temperatures above 273 K, with single-molecule FRET. This method allows the unfolded state dimensions to be probed not only in the cold- and heat-denatured range but also in between, i.e., in the presence of folded protein, and can thus be used to link the two regimes directly. The results show a continuous compaction of unfolded frataxin from 274 to 320 K, with a slight re-expansion at higher temperatures. Cold- and heat-denatured states are thus essentially two sides of the same coin, and their behavior can be understood within the framework of the overall temperature dependence of the unfolded state dimensions.

Effects of sequence and solvation on the temperature-pressure conformational landscape of proteinlike heteropolymers

We study the role of sequence and solvation in shaping the temperature-pressure (T, P) conformational landscape of model heteropolymers with a coarse-grained model. We design foldable primarily hydrophobic sequences with fixed polar content in water at physiological conditions, which demonstrate (T, P) dependence of conformational stability similar to biological proteins. Inherent helicity emerges as a result of local polar-polar interactions in the sequences that mimic biological α helices. The helical propensity is reduced upon solvation and remains unaltered at cold T and high P, which is driven by the T-P induced changes of the hydration shell. Consequently, at nonphysiological conditions the weakening of hydrophobic interactions facilitates population of non-native, helical, compact conformations stabilized through direct nonlocal interactions between polar residues.

Collapse of a hydrophobic polymer in a mixture of denaturants

The solvent quality of an aqueous mixture of two good solvents, urea and guanidinium chloride (GdmCl), for a hydrophobic polymer was investigated using atomistic molecular dynamics simulations. A counterintuitive collapse of the polymer was found, suggesting that mixing the two denaturants reduces the solvent quality. This cononsolvency of the polymer in the urea + GdmCl mixture is found to be caused by the preferential adsorption of urea on the polymer. The polymer collapses as a result of indirect long-range interactions between monomers resulting from the presence of urea clouds surrounding them. Surprisingly, urea behaves as the better solvent in the mixture not because there exists a stronger affinity of the polymer for urea. Instead, attractive interactions between two unlike denaturant molecules combined with the direct dispersion interactions of the polymer with both denaturants determine the solvent quality of the mixture.

Unifying microscopic mechanism for pressure and cold denaturations of proteins

We study the stability of globular proteins as a function of temperature and pressure through NPT simulations of a coarse-grained model. We reproduce the elliptical stability of proteins and highlight a unifying microscopic mechanism for pressure and cold denaturations. The mechanism involves the solvation of nonpolar residues with a thin layer of water. These solvated states have lower volume and lower hydrogen-bond energy compared to other conformations of nonpolar solutes. Hence, these solvated states are favorable at high pressure and low temperature, and they facilitate protein unfolding under these thermodynamical conditions.