Monitoring Hydrogen Exchange During Protein Folding by Fast Pressure Jump NMR Spectroscopy

Connecting Protein Millisecond Conformational Dynamics to Protein Thermal Stability

The stability of protein folded states is crucial for its function, yet the relationship with the protein sequence remains poorly understood. Prior studies have focused on the amino acid composition and thermodynamic couplings within a single folded conformation, overlooking the potential contribution of protein dynamics. Here, we address this gap by systematically analyzing the impact of alanine mutations in the C-terminal β-strand (β5) of ubiquitin, a model protein exhibiting millisecond timescale interconversion between two conformational states differing in the β5 position. Our findings unveil a negative correlation between millisecond dynamics and thermal stability, with alanine substitutions at seemingly flexible C-terminal residues significantly enhancing thermostability. Integrating spectroscopic and computational approaches, we demonstrate that the thermally unfolded state retains a substantial secondary structure but lacks β5 engagement, recapitulating the transition state for millisecond dynamics. Thus, alanine mutations that modulate the stabilities of the folded states with respect to the partially unfolded state impact both the dynamics and stability. Our findings underscore the importance of conformational dynamics with implications for protein engineering and design.

Building a Cost-Efficient High-Pressure Cell for Online High-Field NMR and MRI Using Standard Static Probe Heads: An In Situ Demonstration on Clathrate Hydrate Formation

High-pressure nuclear magnetic resonance (NMR) spectroscopy finds remarkable applications in catalysis, protein biochemistry and biophysics, analytical chemistry, material science, energy, and environmental control but requires expensive probe heads and/or sample cells. This contribution describes the design, construction, and testing of a low-cost 5 mm NMR tube suitable for high-pressure NMR measurements of up to 30 MPa. The sample cell comprises a standard, 5 mm single-crystal sapphire tube that has been fitted to a section of a relatively inexpensive polyether ether ketone (PEEK) HPLC column. PEEK HPLC tubing and connectors enable integration with a gas rig or a standard HPLC pump located outside the stray field of the magnet. The cell is compatible with any 5 mm static NMR probe head, exhibits almost zero background in NMR experiments, and is compatible with any liquid, gas, temperature, or pressure range encountered in HPLC experimentation. A specifically designed transport case enables the safe handling of the pressurized tube outside the probe head. The performance of the setup was evaluated using in situ high-field NMR spectroscopy and MRI performed during the formation of bulk and nanoconfined clathrate hydrates occluding methane, ethane, and hydrogen.

Insights into the Structure of Invisible Conformations of Large Methyl Group Labeled Molecular Machines from High Pressure NMR

Most proteins are highly flexible and can adopt conformations that deviate from the energetically most favorable ground state. Structural information on these lowly populated, alternative conformations is often lacking, despite the functional importance of these states. Here, we study the pathway by which the Dcp1:Dcp2 mRNA decapping complex exchanges between an autoinhibited closed and an open conformation. We make use of methyl Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion (RD) experiments that report on the population of the sparsely populated open conformation as well as on the exchange rate between the two conformations. To obtain volumetric information on the open conformation as well as on the transition state structure we made use of RD measurements at elevated pressures. We found that the open Dcp1:Dcp2 conformation has a lower molecular volume than the closed conformation and that the transition state is close in volume to the closed state. In the presence of ATP the volume change upon opening of the complex increases and the volume of the transition state lies in-between the volumes of the closed and open state. These findings show that ATP has an effect on the volume changes that are associated with the opening-closing pathway of the complex. Our results highlight the strength of pressure dependent NMR methods to obtain insights into structural features of protein conformations that are not directly observable. As our work makes use of methyl groups as NMR probes we conclude that the applied methodology is also applicable to high molecular weight complexes.

SABRE Hyperpolarization with up to 200 bar Parahydrogen in Standard and Quickly Removable Solvents

Parahydrogen (p-H₂)-based techniques are known to drastically enhance NMR signals but are usually limited by p-H₂ supply. This work reports p-H₂-based SABRE hyperpolarization at p-H₂ pressures of hundreds of bar, far beyond the typical ten bar currently reported in the literature. A recently designed high-pressure setup was utilized to compress p-H₂ gas up to 200 bar. The measurements were conducted using a sapphire high-pressure NMR tube and a 43 MHz benchtop NMR spectrometer. In standard methanol solutions, it could be shown that the signal intensities increased with pressure until they eventually reached a plateau. A polarization of about 2%, equal to a molar polarization of 1.2 mmol L^-1, could be achieved for the sample with the highest substrate concentration. While the signal plateaued, the H₂ solubility increased linearly with pressure from 1 to 200 bar, indicating that p-H₂ availability is not the limiting factor in signal enhancement beyond a certain pressure, depending on sample composition. Furthermore, the possibility of using liquefied ethane and compressed CO₂ as removable solvents for hyperpolarization was demonstrated. The use of high pressures together with quickly removable organic/non-organic solvents represents an important breakthrough in the field of hyperpolarization, advancing SABRE as a promising tool for materials science, biophysics, and molecular imaging.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

Probing the solvation of the α-helix with extended amide III bands in Raman optical activity

Experimental and theoretical Raman optical activity (ROA) study of α-helical peptides and proteins has suggested that the relative intensity of two extended amide III ROA bands at ∼1340 cm^-1 (I band) and ∼1300 cm^-1 (II band) can be used to monitor the permittivity of the surrounding medium of the α-helix. So far, the ROA intensity ratio, I_I/I_II, has been interpreted from two different viewpoints. The first one is in terms of a direct effect of permittivity around the α-helix. The second one is based on a structural equilibrium of two types of α-helical structures, "hydrated" and "unhydrated" ones. In the present study, temperature- and solvent-dependences of I_I/I_II are measured for highly-α-helical peptides and compared to the theoretical spectra while varying the permittivity or the type of α-helical structure. A fragment method with partial optimization in the normal modes is adopted in density functional theory calculations. The main features of the experimental spectra and a trend of the observed I_I/I_II are well reproduced by the simulations, which leads us to a conclusion that the I_I/I_II is dominantly governed by a direct influence of the permittivity of the environment and just accessorily by the equilibrium of the two types of α-helices. The simulations also opposed the conventional assignments of the I and II bands to "hydrated" and "unhydrated" α-helical structures, respectively. In the case of α-helical proteins, solvent exposure of the α-helix may be monitored by the ROA ratio.

Transition-State Compressibility and Activation Volume of Transient Protein Conformational Fluctuations

Proteins are dynamic entities that intermittently depart from their ground-state structures and undergo conformational transitions as a critical part of their functions. Central to understanding such transitions are the structural rearrangements along the connecting pathway, where the transition state plays a special role. Using NMR relaxation at variable temperature and pressure to measure aromatic ring flips inside a protein core, we obtain information on the structure and thermodynamics of the transition state. We show that the isothermal compressibility coefficient of the transition state is similar to that of short-chain hydrocarbon liquids, implying extensive local unfolding of the protein. Our results further indicate that the required local volume expansions of the protein can occur not only with a net positive activation volume of the protein, as expected from previous studies, but also with zero activation volume by compaction of remote void volume, when averaged over the ensemble of states.

Real-time nuclear magnetic resonance spectroscopy in the study of biomolecular kinetics and dynamics

The review describes the application of nuclear magnetic resonance (NMR) spectroscopy to study kinetics of folding, refolding and aggregation of proteins, RNA and DNA. Time-resolved NMR experiments can be conducted in a reversible or an irreversible manner. In particular, irreversible folding experiments pose large requirements for (i) signal-to-noise due to the time limitations and (ii) synchronising of the refolding steps. Thus, this contribution discusses the application of methods for signal-to-noise increases, including dynamic nuclear polarisation, hyperpolarisation and photo-CIDNP for the study of time-resolved NMR studies. Further, methods are reviewed ranging from pressure and temperature jump, light induction to rapid mixing to induce rapidly non-equilibrium conditions required to initiate folding.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (R_h) of biomolecules in solution. For unfolded proteins, the R_h provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the R_h of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by ¹³C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that R_h values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the R_h of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Experimental estimates of compression heating and decompression cooling in ethylene glycol

The chemical shift difference, Δσ, between the methylene and hydroxyl protons in the high resolution ¹ H nuclear magnetic resonance spectrum of ethylene glycol is shown to be pressure dependent. The equilibrium Δσ values for ethylene glycol are reported as a function of temperature and pressure between ambient conditions, 323 K and 2 kbar, respectively. This surface is used along with Δσ values measured in response to a rapid pressure increase to calculate a temperature rise that is used to infer a temperature change for water that is consistent with theoretical estimates. This work implies that compression heating and decompression cooling are not significant enough to interfere with pressure induced protein folding studies.

Fast electron paramagnetic resonance magic angle spinning simulations using analytical powder averaging techniques

Simulations describing the spin physics underpinning nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopy play an important role in the design of new experiments. When experiments are performed in the solid state, samples are commonly composed of powders or glasses, with molecules oriented at a large number of angles with respect to the laboratory frame. These powder angles must be represented in simulations to account for anisotropic interactions. Numerical techniques are typically used to accurately compute such powder averages. A large number of Euler angles are usually required, leading to lengthy simulation times. This is particularly true in broad spectra, such as those observed in EPR. The combination of the traditionally separate techniques of EPR and magic angle spinning (MAS) NMR could play an important role in future electron detected experiments, combined with dynamic nuclear polarization, which will allow for exceptional detection sensitivity of NMR spin coherences. Here, we present a method of reducing the required number of Euler angles in magnetic resonance simulations by analytically performing the powder average over one of the Euler angles in the static and MAS cases for the TEMPO nitroxide radical in a 7 T field. In the static case, this leads to a 97.5% reduction in simulation time over the fully numerical case and reproduces the expected spinning sideband manifold when simulated with a MAS frequency of 150 kHz. This technique is applicable to more traditional NMR experiments as well, such as those involving quadrupolar nuclei or multiple dimensions.

Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics

As theory and experiment have shown, protein dehydration is a major contributor to protein folding. Dehydration upon folding can be characterized directly by all-atom simulations of fast pressure drops, which create desolvated pockets inside the nascent hydrophobic core. Here, we study pressure-drop refolding of three λ-repressor fragment (λ_6-85) mutants computationally and experimentally. The three mutants report on tertiary structure formation via different fluorescent helix-helix contact pairs. All-atom simulations of pressure drops capture refolding and unfolding of all three mutants by a similar mechanism, thus validating the nonperturbative nature of the fluorescent contact probes. Analysis of simulated interprobe distances shows that the α-helix 1-3 pair distance displays a slower characteristic time scale than the 1-2 or 3-2 pair distance. To see whether slow packing of α-helices 1 and 3 is reflected in the rate-limiting folding step, fast pressure-drop relaxation experiments captured refolding on a millisecond time scale. These experiments reveal that refolding monitored by 1-3 contact formation indeed is much slower than when monitored by 1-2 or 3-2 contact formation. Unlike the case of the two-state folder [three-α-helix bundle (α₃D)], whose drying and core formation proceed in concert, λ_6-85 repeatedly dries and rewets different local tertiary contacts before finally forming a solvent-excluded core, explaining the non-two-state behavior observed during refolding in molecular dynamics simulations. This work demonstrates that proteins can explore desolvated pockets and dry globular states numerous times before reaching the native conformation.

Interrupted Pressure-Jump NMR Experiments Reveal Resonances of On-Pathway Protein Folding Intermediate

Previous pressure-jump NMR experiments on a pressure-sensitized double mutant of ubiquitin showed evidence that its folding occurs via two parallel, comparably efficient pathways: a single barrier and a two-barrier pathway. An interrupted folding NMR experiment is introduced, where for a brief period the pressure is dropped to atmospheric conditions (1 bar), followed by a jump back to high pressure for signal detection. Conventional, forward sampling of the indirect dimension during the low-pressure period correlates the ¹⁵N or ¹³C' chemical shifts of the unfolded protein at 1 bar to the ¹H frequencies of both the unfolded and folded proteins at high pressure. Remarkably, sampling the data of the same experiment in the reverse direction yields the frequencies of proteins present at the end of the low-pressure interval, which include unfolded, intermediate, and folded species. Although the folding intermediate ¹⁵N shifts differ strongly from natively folded protein, its ¹³C' chemical shifts, which are more sensitive probes for secondary structure, closely match those of the folded protein and indicate that the folding intermediate must have a structure that is quite similar to the native state.

Practical aspects of high-pressure NMR spectroscopy and its applications in protein biophysics and structural biology

Pressure and temperature are the two fundamental variables of thermodynamics. Temperature and chemical perturbation are central experimental tools for the exploration of macromolecular structure and dynamics. Though it has long been recognized that hydrostatic pressure offers a complementary and often unique view of macromolecular structure, stability and dynamics, it has not been employed nearly as much. For solution NMR applications the limited use of high-pressure is undoubtedly traced to difficulties of employing pressure in the context of modern multinuclear and multidimensional NMR. Limitations in pressure tolerant NMR sample cells have been overcome and enable detailed studies of macromolecular energy landscapes, hydration, dynamics and function. Here we review the practical considerations for studies of biological macromolecules at elevated pressure, with a particular emphasis on applications in protein biophysics and structural biology.

Equilibrium and Kinetic Unfolding of GB1: Stabilization of the Native State by Pressure

NMR spectroscopy allows an all-atom view on pressure-induced protein folding, separate detection of different folding states, determination of their population, and the measurement of the folding kinetics at equilibrium. Here, we studied the folding of protein GB1 at pH 2 in a temperature and pressure dependent way. We find that the midpoints of temperature-induced unfolding increase with higher pressure. NMR relaxation dispersion experiments disclosed that the unfolding kinetics slow down at elevated pressure while the folding kinetics stay virtually the same. Therefore, pressure is stabilizing the native state of GB1. These findings extend the knowledge of the influence of pressure on protein folding kinetics, where so far typically a destabilization by increased activation volumes of folding was observed. Our findings thus point toward an exceptional section in the pressure-temperature phase diagram of protein unfolding. The stabilization of the native state could potentially be caused by a shift of p K_a values of glutamates and aspartates in favor of the negatively charged state as judged from pH sensitive chemical shifts.

Monitoring 15N Chemical Shifts During Protein Folding by Pressure-Jump NMR

Pressure-jump hardware permits direct observation of protein NMR spectra during a cyclically repeated protein folding process. For a two-state folding protein, the change in resonance frequency will occur nearly instantaneously when the protein clears the transition state barrier, resulting in a monoexponential change of the ensemble-averaged chemical shift. However, protein folding pathways can be more complex and contain metastable intermediates. With a pseudo-3D NMR experiment that utilizes stroboscopic observation, we measure the ensemble-averaged chemical shifts, including those of exchange-broadened intermediates, during the folding process. Such measurements for a pressure-sensitized mutant of ubiquitin show an on-pathway kinetic intermediate whose ¹⁵N chemical shifts differ most from the natively folded protein for strands β5, its preceding turn, and the two strands that pair with β5 in the native structure.

Study of protein folding under native conditions by rapidly switching the hydrostatic pressure inside an NMR sample cell

Development of specialized instrumentation enables rapid switching of the hydrostatic pressure inside an operating NMR spectrometer. This technology allows observation of protein signals during the repeated folding process. Applied to ubiquitin, a previously extensively studied model of protein folding, the methodology reveals an initially highly dynamic state that deviates relatively little from random coil behavior but also provides evidence for numerous repeatedly failed folding events, previously only observed in computer simulations. Above room temperature, direct NMR evidence shows a ∼50% fraction of proteins folding through an on-pathway kinetic intermediate, thereby revealing two equally efficient parallel folding pathways.

In general, small proteins rapidly fold on the timescale of milliseconds or less. For proteins with a substantial volume difference between the folded and unfolded states, their thermodynamic equilibrium can be altered by varying the hydrostatic pressure. Using a pressure-sensitized mutant of ubiquitin, we demonstrate that rapidly switching the pressure within an NMR sample cell enables study of the unfolded protein under native conditions and, vice versa, study of the native protein under denaturing conditions. This approach makes it possible to record 2D and 3D NMR spectra of the unfolded protein at atmospheric pressure, providing residue-specific information on the folding process. ¹⁵N and ¹³C chemical shifts measured immediately after dropping the pressure from 2.5 kbar (favoring unfolding) to 1 bar (native) are close to the random-coil chemical shifts observed for a large, disordered peptide fragment of the protein. However, ¹⁵N relaxation data show evidence for rapid exchange, on a ∼100-μs timescale, between the unfolded state and unstable, structured states that can be considered as failed folding events. The NMR data also provide direct evidence for parallel folding pathways, with approximately one-half of the protein molecules efficiently folding through an on-pathway kinetic intermediate, whereas the other half fold in a single step. At protein concentrations above ∼300 μM, oligomeric off-pathway intermediates compete with folding of the native state.

Propensity for cis-Proline Formation in Unfolded Proteins

In unfolded proteins, peptide bonds involving Pro residues exist in equilibrium between the minor cis and major trans conformations. Folded proteins predominantly contain trans-Pro bonds, and slow cis-trans Pro isomerization in the unfolded state is often found to be a rate-limiting step in protein folding. Moreover, kinases and phosphatases that act upon Ser/Thr-Pro motifs exhibit preferential recognition of either the cis- or trans-Pro conformer. Here, NMR spectra obtained at both atmospheric and high pressures indicate that the population of cis-Pro falls well below previous estimates, an effect attributed to the use of short peptides with charged termini in most prior model studies. For the intrinsically disordered protein α-synuclein, cis-Pro populations at all of its five X-Pro bonds are less than 5 %, with only modest ionic strength dependence and no detectable effect of the previously demonstrated interaction between the N- and C-terminal halves of the protein. Comparison to small peptides with the same amino-acid sequence indicates that peptides, particularly those with unblocked, oppositely charged amino and carboxyl end groups, strongly overestimate the amount of cis-Pro.