Beyond slow two-state protein conformational exchange using CEST: applications to three-state protein interconversion on the millisecond timescale.

Although NMR spectroscopy is routinely used to study the conformational dynamics of biomolecules, robust analyses of the data are challenged in cases where exchange is more complex than two-state, such as when a 'visible' major conformer exchanges with two 'invisible' minor states on the millisecond timescale. It is becoming increasingly clear that chemical exchange saturation transfer (CEST) NMR experiments that were initially developed to study systems undergoing slow interconversion are also sensitive to intermediate-fast timescale biomolecular conformational exchange. Here we investigate the utility of the amide 15N CEST experiment to characterise protein three-state exchange occurring on the millisecond timescale by studying the interconversion between the folded (F) state of the FF domain from human HYPA/FBP11 (WT FF) and two of its folding intermediates I1 and I2. Although 15N CPMG experiments are consistent with the F state interconverting with a single minor state on the millisecond timescale, 15N CEST data clearly establish an exchange process between F and a pair of minor states. A unique three-state exchange model cannot be obtained by analysis of 15N CEST data recorded at a single temperature. However, including the relative sign of the difference in the chemical shifts of the two minor states based on a simple two-state analysis of CEST data recorded at multiple temperatures, results in a robust three-state model in which the F, I1 and I2 states interconvert with each other on the millisecond timescale ( k e x , F I 1 ~ 550 s-1, k e x , F I 2 ~ 1200 s-1, k e x , I 1 I 2 ~ 5000 s-1), with I1 and I2 sparsely populated at ~ 0.15% and ~ 0.35%, respectively, at 15 °C. A computationally demanding grid-search of exchange parameter space is not required to extract the best-fit exchange parameters from the CEST data. The utility of the CEST experiment, thus, extends well beyond studies of conformers in slow exchange on the NMR chemical shift timescale, to include systems with interconversion rates on the order of thousands/second.

Probing excited state 1Hα chemical shifts in intrinsically disordered proteins with a triple resonance-based CEST experiment: Application to a disorder-to-order switch.

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure-function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of 'invisible' excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i + 1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints.

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 µs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways.

Conformational dynamics play critical roles in protein folding, misfolding, function, misfunction, and aggregation. While detecting and studying the different conformational states populated by protein molecules on their free energy surfaces (FESs) remain a challenge, NMR spectroscopy has emerged as an invaluable experimental tool to explore the FES of a protein, as conformational dynamics can be probed at atomic resolution over a wide range of timescales. Here, we use chemical exchange saturation transfer (CEST) to detect "invisible" minor states on the energy landscape of the A39G mutant FF domain that exhibited "two-state" folding kinetics in traditional experiments. Although CEST has mostly been limited to studies of processes with rates between ∼5 to 300 s-1 involving sparse states with populations as low as ∼1%, we show that the line broadening that is often associated with minor state dips in CEST profiles can be exploited to inform on additional conformers, with lifetimes an order of magnitude shorter and populations close to 10-fold smaller than what typically is characterized. Our analysis of CEST profiles that exploits the minor state linewidths of the 71-residue A39G FF domain establishes a folding mechanism that can be described in terms of a four-state exchange process between interconverting states spanning over two orders of magnitude in timescale from ∼100 to ∼15,000 µs. A similar folding scheme is established for the wild-type domain as well. The study shows that the folding of this small domain proceeds through a pair of sparse, partially structured intermediates via two discrete pathways on a volcano-shaped FES.

A CEST NMR experiment to obtain glycine 1Hα chemical shifts in 'invisible' minor states of proteins.

Chemical exchange saturation transfer (CEST) experiments are routinely used to study protein conformational exchange between a 'visible' major state and 'invisible' minor states because they can detect minor states with lifetimes varying from ~ 3 to ~ 100 ms populated to just ~ 0.5%. Consequently several 1H, 15N and 13C CEST experiments have been developed to study exchange and obtain minor state chemical shifts at almost all backbone and sidechain sites in proteins. Conspicuously missing from this extensive set of CEST experiments is a 1H CEST experiment to study exchange at glycine (Gly) 1Hα sites as the existing 1H CEST experiments that have been designed to study dynamics in amide 1H-15N spin systems and methyl 13CH3 groups with three equivalent protons while suppressing 1H-1H NOE induced dips are not suitable for studying exchange in methylene 13CH2 groups with inequivalent protons. Here a Gly 1Hα CEST experiment to obtain the minor state Gly 1Hα chemical shifts is presented. The utility of this experiment is demonstrated on the L99A cavity mutant of T4 Lysozyme (T4L L99A) that undergoes conformational exchange between two compact conformers. The CEST derived minor state Gly 1Hα chemical shifts of T4L L99A are in agreement with those obtained previously using CPMG techniques. The experimental strategy presented here can also be used to obtain methylene proton minor state chemical shifts from protein sidechain and nucleic acid backbone sites.

Confronting the Invisible: Assignment of Protein 1HN Chemical Shifts in Cases of Extreme Broadening.

NMR studies of intrinsically disordered proteins (IDPs) at neutral pH values are hampered by the rapid exchange of backbone amide protons with solvent. Although exchange rates can be modulated by changes in pH, interactions between IDPs that lead to phase separation sometimes only occur at neutral pH values or higher, where backbone amide-based experiments fail. Here we describe a simple NMR experiment for measuring amide proton chemical shifts in cases where 1HN spectra cannot be obtained. The approach uses a weak 1H B1 field, searching for elusive 1HN resonance frequencies that become encoded in the intensities of cross-peaks in three-dimensional 1Hα-detect spectra. Applications to the CAPRIN1 protein in both dilute- and phase-separated states highlight the utility of the method, establishing that accurate 1HN chemical shifts can be obtained even in cases where solvent hydrogen exchange rates are on the order of 1500 s-1.

A Double-Resonance CEST Experiment To Study Multistate Protein Conformational Exchange: An Application to Protein Folding.

Despite the importance of protein dynamics to function, studying exchange between multiple conformational states remains a challenge because sparsely populated states are invisible to conventional techniques. CEST NMR experiments can detect minor states with lifetimes between 5 and 200 ms populated to a level of just ∼1%. However, CEST often cannot provide the exchange mechanism for processes involving three or more states, leaving the role of the detected minor states unknown. Here a double-resonance CEST experiment to determine the kinetics of multistate exchange is presented. The approach that involves irradiating resonances from two minor states simultaneously is used to study the exchange of T4 lysozyme (T4L) between the dominant native state and two minor states, the unfolded state and a second minor state (B), each populated to only ∼4%. Regular CEST does not provide the folding mechanism, but double-resonance CEST clearly shows that T4L can fold directly without going through B.

A Methyl-TROSY-Based 1 H Relaxation Dispersion Experiment for Studies of Conformational Exchange in High Molecular Weight Proteins.

Molecular complexes often sample conformational states that direct them to specific functions. These states can be difficult to observe through traditional biophysical approaches but they can be studied using a variety of different NMR spin relaxation experiments. However, these applications, when focused on moderate to high molecular weight proteins, are complicated by fast relaxing signals that negatively affect the sensitivity and resolution of spectra. Here a methyl 1 H CPMG-based experiment for studies of excited conformational states of protein machines is described that exploits a TROSY-effect to increase signal-to-noise. Complexities from the multiplicity of methyl 1 H transitions are addressed to generate a robust pulse scheme that is applied to a 320 kDa homeostasis protein, p97.

Exchangeable deuterons introduce artifacts in amide 15N CEST experiments used to study protein conformational exchange.

Protein molecules sample different conformations in solution and characterizing these conformations is crucial to understanding protein function. 15N CEST experiments are now routinely used to study slow conformational exchange of protein molecules between a 'visible' major state and 'invisible' minor states. These experiments have also been adapted to measure the solvent exchange rates of amide protons by exploiting the one bond deuterium isotope effect on the amide 15N chemical shifts. However at moderately high temperatures (~ 50 °C) that are sometimes required to populate protein minor conformers to levels (~ 1%) that can be detected by CEST experiments solvent H/D exchange can lead to 'dips' in low B115N CEST profiles that can be wrongly assigned to the conformational exchange process being characterized. This is demonstrated in the case of ~ 18 kDa T4 Lysozyme (T4L) at 50 °C and the ~ 11 kDa E. coli hibernation promoting factor (HPF) at 52 °C. This problem is trivially solved by eliminating the exchangeable deuterons in the solvent by using either an external D2O lock or by using a small amount (~ 1-3%) of a molecule like d6-DMSO that does not contain exchangeable deuterons to lock the spectrometer.

Measuring Diffusion Constants of Invisible Protein Conformers by Triple-Quantum 1 H CPMG Relaxation Dispersion.

Proteins are not locked in a single structure but often interconvert with other conformers that are critical for function. When such conformers are sparsely populated and transiently formed they become invisible to routine biophysical methods, however they can be studied in detail by NMR spin-relaxation experiments. Few experiments are available in the NMR toolkit, however, for characterizing the hydrodynamic properties of invisible states. Herein we describe a CPMG-based experiment for measuring translational diffusion constants of invisible states using a pulsed-field gradient approach that exploits methyl 1 H triple-quantum coherences. An example, involving diffusion of a sparsely populated and hence invisible unfolded protein ensemble is presented, without the need for the addition of denaturants that tend to destroy weak interactions that can be involved in stabilizing residual structure in the unfolded state.

A methyl 1H double quantum CPMG experiment to study protein conformational exchange.

Protein conformational changes play crucial roles in enabling function. The Carr-Purcell-Meiboom-Gill (CPMG) experiment forms the basis for studying such dynamics when they involve the interconversion between highly populated and sparsely formed states, the latter having lifetimes ranging from ~ 0.5 to ~ 5 ms. Among the suite of experiments that have been developed are those that exploit methyl group probes by recording methyl 1H single quantum (Tugarinov and Kay in J Am Chem Soc 129:9514-9521, 2007) and triple quantum (Yuwen et al. in Angew Chem Int Ed Engl 55:11490-11494, 2016) relaxation dispersion profiles. Here we build upon these by developing a third experiment in which methyl 1H double quantum coherences evolve during a CPMG relaxation element. By fitting single, double, and triple quantum datasets, akin to recording the single quantum dataset at static magnetic fields of Bo, 2Bo and 3Bo, we show that accurate exchange values can be obtained even in cases where exchange rates exceed 10,000 s-1. The utility of the double quantum experiment is demonstrated with a pair of cavity mutants of T4 lysozyme (T4L) with ground and excited states interchanged and with exchange rates differing by fourfold (~ 900 s-1 and ~ 3600 s-1), as well as with a fast-folding domain where the unfolded state lifetime is ~ 80 µs.

Atomic resolution mechanism of ligand binding to a solvent inaccessible cavity in T4 lysozyme.

Ligand binding sites in proteins are often localized to deeply buried cavities, inaccessible to bulk solvent. Yet, in many cases binding of cognate ligands occurs rapidly. An intriguing system is presented by the L99A cavity mutant of T4 Lysozyme (T4L L99A) that rapidly binds benzene (~106 M-1s-1). Although the protein has long served as a model system for protein thermodynamics and crystal structures of both free and benzene-bound T4L L99A are available, the kinetic pathways by which benzene reaches its solvent-inaccessible binding cavity remain elusive. The current work, using extensive molecular dynamics simulation, achieves this by capturing the complete process of spontaneous recognition of benzene by T4L L99A at atomistic resolution. A series of multi-microsecond unbiased molecular dynamics simulation trajectories unequivocally reveal how benzene, starting in bulk solvent, diffuses to the protein and spontaneously reaches the solvent inaccessible cavity of T4L L99A. The simulated and high-resolution X-ray derived bound structures are in excellent agreement. A robust four-state Markov model, developed using cumulative 60 µs trajectories, identifies and quantifies multiple ligand binding pathways with low activation barriers. Interestingly, none of these identified binding pathways required large conformational changes for ligand access to the buried cavity. Rather, these involve transient but crucial opening of a channel to the cavity via subtle displacements in the positions of key helices (helix4/helix6, helix7/helix9) leading to rapid binding. Free energy simulations further elucidate that these channel-opening events would have been unfavorable in wild type T4L. Taken together and via integrating with results from experiments, these simulations provide unprecedented mechanistic insights into the complete ligand recognition process in a buried cavity. By illustrating the power of subtle helix movements in opening up multiple pathways for ligand access, this work offers an alternate view of ligand recognition in a solvent-inaccessible cavity, contrary to the common perception of a single dominant pathway for ligand binding.

Measuring the signs of the methyl 1H chemical shift differences between major and 'invisible' minor protein conformational states using methyl 1H multi-quantum spectroscopy.

Carr-Purcell-Meiboom-Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the µs to millisecond (ms) timescale between a visible major state and 'invisible' minor states. The exchange rate(s) ([Formula: see text]), population(s) of the minor state(s) and the absolute value of the chemical shift difference [Formula: see text] (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of [Formula: see text] that is required to reconstruct the spectrum of the 'invisible' minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl [Formula: see text] CPMG experiment (Yuwen in Angew Chem 55:11490-11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of [Formula: see text] is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~ 18 kDa cavity mutant of T4 Lysozyme ([Formula: see text]) and (2) the [Formula: see text] kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus ([Formula: see text]).

CPMG Experiments for Protein Minor Conformer Structure Determination.

CPMG relaxation dispersion NMR experiments have emerged as a powerful method to characterize protein minor states that are in exchange with a visible dominant conformation, and have lifetimes between ~0.5 and 5 milliseconds (ms) and populations greater than 0.5%. The structure of the minor state can, in favorable cases, be determined from the parameters provided by the CPMG relaxation dispersion experiments. Here, we go through the intricacies of setting up these powerful CPMG experiments.

Probing conformational dynamics in biomolecules via chemical exchange saturation transfer: a primer.

Although Chemical Exchange Saturation Transfer (CEST) type NMR experiments have been used to study chemical exchange processes in molecules since the early 1960s, there has been renewed interest in the past several years in using this approach to study biomolecular conformational dynamics. The methodology is particularly powerful for the study of sparsely populated, transiently formed conformers that are recalcitrant to investigation using traditional biophysical tools, and it is complementary to relaxation dispersion and magnetization transfer experiments that have traditionally been used to study chemical exchange processes. Here we discuss the concepts behind the CEST experiment, focusing on practical aspects as well, we review some of the pulse sequences that have been developed to characterize protein and RNA conformational dynamics, and we discuss a number of examples where the CEST methodology has provided important insights into the role of dynamics in biomolecular function.

Simultaneous NMR characterisation of multiple minima in the free energy landscape of an RNA UUCG tetraloop.

RNA molecules in solution tend to undergo structural fluctuations of relatively large amplitude and to populate a range of different conformations some of which with low populations. It is still very challenging, however, to characterise the structures of these low populated states and to understand their functional roles. In the present study, we address this problem by using NMR residual dipolar couplings (RDCs) as structural restraints in replica-averaged metadynamics (RAM) simulations. By applying this approach to a 14-mer RNA hairpin containing the prototypical UUCG tetraloop motif, we show that it is possible to construct the free energy landscape of this RNA molecule. This free energy landscapes reveals the surprisingly rich dynamics of the UUCG tetraloop and identifies the multiple substates that exist in equilibrium owing to thermal fluctuations. The approach that we present is general and can be applied to the study of the free energy landscapes of other RNA or RNA-protein systems.

Enhancing the Sensitivity of CPMG Relaxation Dispersion to Conformational Exchange Processes by Multiple-Quantum Spectroscopy.

A triple-quantum (1) H Carr-Purcell-Meiboom-Gill NMR relaxation dispersion experiment is presented that uses methyl group probes as reporters of conformational exchange in highly deuterated, methyl-protonated proteins. Significantly larger dispersion profiles, by as much as a factor of nine, can be obtained relative to single-quantum approaches, thus offering very significant advantages in applications involving interconverting conformers with only small changes in structure or in studies of rare states that are at very low populations. Applications to a number of protein systems are presented where the utility of the method, including its improved sensitivity to chemical exchange processes, is established.

Atomistic picture of conformational exchange in a T4 lysozyme cavity mutant: an experiment-guided molecular dynamics study.

Despite the importance of dynamics to protein function there is little information about the states that are formed as the protein explores its conformational landscape or about the mechanism by which transitions between the different states occur. Here we used a combined NMR spin relaxation and unbiased molecular dynamics (MD) approach to investigate the exchange process by which a cavity in an L99A mutant of T4 lysozyme (T4L 99A) interconverts between an empty and occupied form that involves repositioning of an aromatic residue, Phe114. Although structures of the end-states of the exchange process are available, insight into the mechanism by which the transition takes place cannot be obtained from experiment and the timescales involved are too slow to address using brute force MD. Using spin relaxation NMR methods, we have identified a triple-mutant of T4L that undergoes the same exchange process as T4L L99A but where the minor state lifetime has decreased significantly so that the spontaneous conformational transition to the major state can be studied using all-atom MD simulations. The simulation trajectories were analyzed using Markov state models and the energy landscape so obtained is in good agreement with expectations based on NMR studies. Notably there is no large-scale perturbation of the structure during the transition, multiple intermediates are formed between the two similar exchanging conformations and the free energy barrier between these two well-folded, compact forms is small (6kBT), only slightly larger than for processes considered to be barrierless.

Probing the free energy landscape of the fast-folding gpW protein by relaxation dispersion NMR.

The topographic features of the free energy landscapes that govern the thermodynamics and kinetics of conformational transitions in proteins, which in turn are integral for function, are not well understood. This reflects the experimental challenges associated with characterizing these multidimensional surfaces, even for small proteins. Here we focus on a 62-residue protein, gpW, that folds very rapidly into a native structure with an α/ß topology in which α-helices are at the N- and C-terminal ends of the molecule with a central ß-hairpin positioned orthogonally to the helices. Using relaxation dispersion NMR spectroscopy to probe the conformational fluctuations in gpW at 1 °C, we found that the native state interconverts with a transiently formed, sparsely populated second state with a lifetime of 250 µs, consistent with the global folding-unfolding rate under these conditions. In this low-populated state, the ß-hairpin is unfolded whereas the α-helices remain predominantly formed. Our results argue for a hierarchical stability of secondary structural elements and demonstrate the existence of a complex free energy landscape even in this small, fast-folding single-domain protein.

Viscosity-dependent kinetics of protein conformational exchange: microviscosity effects and the need for a small viscogen.

Conformational rearrangements are critical to a variety of biological processes including protein folding and misfolding, ligand binding, enzyme catalysis, and signal transduction. Viscosity-dependent kinetics measurements can provide crucial insights into the dynamics of protein conformational exchange by highlighting the relative importance of frictional forces derived from either solvent or from internal protein interactions in activating the exchange reaction. Here, we analyze the kinetics of interconversion between the native and intermediate states of the four helix bundle FF domain recorded in solutions containing the viscogens glycerol or bovine serum albumin (BSA), using the viscosity measured from the translational diffusion of probes of different sizes. In the large viscogen BSA, we demonstrate that vastly different internal friction values are obtained using the different viscosity measures, leading to conflicting interpretations of the role of solvent friction in the interconversion. We show that this can be a consequence of the small effective hydrodynamic radius of the protein conformational transition and differences between solution micro- and macroscopic viscosities that are germane in this case. In general, correct values of internal friction can only be obtained by carrying out measurements using small viscogens.