Using relaxation dispersion NMR spectroscopy to determine structures of excited, invisible protein states

The conformational landscape of fold-switcher KaiB is tuned to the circadian rhythm timescale

How can a single protein domain encode a conformational landscape with multiple stably-folded states, and how do those states interconvert? Here, we use real-time and relaxation-dispersion NMR to characterize the conformational landscape of the circadian rhythm protein KaiB from Rhodobacter sphaeroides. Unique among known natural metamorphic proteins, this KaiB variant spontaneously interconverts between two monomeric states: the "Ground" and "Fold-switched" (FS) state. KaiB in its FS state interacts with multiple binding partners, including the central KaiC protein, to regulate circadian rhythms. We find that KaiB itself takes hours to interconvert between the Ground and FS state, underscoring the ability of a single sequence to encode the slow process needed for function. We reveal the rate-limiting step between the Ground and FS state is the cis-trans isomerization of three prolines in the fold-switching region by demonstrating interconversion acceleration by the prolyl isomerase CypA. The interconversion proceeds through a "partially disordered" (PD) state, where the C-terminal half becomes disordered while the N-terminal half remains stably folded. We discovered two additional properties of KaiB's landscape. Firstly, the Ground state experiences cold denaturation: at 4°C, the PD state becomes the majorly populated state. Secondly, the Ground state exchanges with a fourth state, the "Enigma" state, on the millisecond timescale. We combine AlphaFold2-based predictions and NMR chemical shift predictions to predict this "Enigma" state is a beta-strand register shift that eases buried charged residues, and support this structure experimentally. These results provide mechanistic insight in how evolution can design a single sequence that achieves specific timing needed for its function.

Increasing the accuracy of exchange parameters reporting on slow dynamics by performing CEST experiments with 'high' B1 fields

Over the last decade chemical exchange saturation transfer (CEST) NMR methods have emerged as powerful tools to characterize biomolecular conformational dynamics occurring between a visible major state and 'invisible' minor states. The ability of the CEST experiment to detect these minor states, and provide precise exchange parameters, hinges on using appropriate B1 field strengths during the saturation period. Typically, a pair of B1 fields with ω1 (=2πB1) values around the exchange rate kex are chosen. Here we show that the transverse relaxation rate of the minor state resonance (R2,B) also plays a crucial role in determining the B1 fields that lead to the most informative datasets. Using [Formula: see text]  ≥ kex, to guide the choice of B1, instead of kex, leads to data wherefrom substantially more accurate exchange parameters can be derived. The need for higher B1 fields, guided by K, is demonstrated by studying the conformational exchange in two mutants of the 71 residue FF domain with kex ∼ 11 s-1 and ∼ 72 s-1, respectively. In both cases analysis of CEST datasets recorded using B1 field values guided by kex lead to imprecise exchange parameters, whereas using B1 values guided by K resulted in precise site-specific exchange parameters. The conclusions presented here will be valuable while using CEST to study slow processes at sites with large intrinsic relaxation rates, including carbonyl sites in small to medium sized proteins, amide 15N sites in large proteins and when the minor state dips are broadened due to exchange among the minor states.

A unique chaperoning mechanism in class A JDPs recognizes and stabilizes mutant p53

J-domain proteins (JDPs) constitute a large family of molecular chaperones that bind a broad spectrum of substrates, targeting them to Hsp70, thus determining the specificity of and activating the entire chaperone functional cycle. The malfunction of JDPs is therefore inextricably linked to myriad human disorders. Here, we uncover a unique mechanism by which chaperones recognize misfolded clients, present in human class A JDPs. Through a newly identified β-hairpin site, these chaperones detect changes in protein dynamics at the initial stages of misfolding, prior to exposure of hydrophobic regions or large structural rearrangements. The JDPs then sequester misfolding-prone proteins into large oligomeric assemblies, protecting them from aggregation. Through this mechanism, class A JDPs bind destabilized p53 mutants, preventing clearance of these oncoproteins by Hsp70-mediated degradation, thus promoting cancer progression. Removal of the β-hairpin abrogates this protective activity while minimally affecting other chaperoning functions. This suggests the class A JDP β-hairpin as a highly specific target for cancer therapeutics.

Combined NMR and molecular dynamics conformational filter identifies unambiguously dynamic ensembles of Dengue protease NS2B/NS3pro

The dengue protease NS2B/NS3pro has been reported to adopt either an 'open' or a 'closed' conformation. We have developed a conformational filter that combines NMR with MD simulations to identify conformational ensembles that dominate in solution. Experimental values derived from relaxation parameters for the backbone and methyl side chains were compared with the corresponding back-calculated relaxation parameters of different conformational ensembles obtained from free MD simulations. Our results demonstrate a high prevalence for the 'closed' conformational ensemble while the 'open' conformation is absent, indicating that the latter conformation is most probably due to crystal contacts. Conversely, conformational ensembles in which the positioning of the co-factor NS2B results in a 'partially' open conformation, previously described in both MD simulations and X-ray studies, were identified by our conformational filter. Altogether, we believe that our approach allows for unambiguous identification of true conformational ensembles, an essential step for reliable drug discovery.

Discerning conformational dynamics and binding kinetics of GPCRs by 19F NMR

19F NMR provides a way of monitoring conformational dynamics of G-protein coupled receptors (GPCRs) from the perspective of an ensemble. While X-ray crystallography provides exquisitely resolved high-resolution structures of specific states, it generally does not recapitulate the true ensemble of functional states. Fluorine (19F) NMR provides a highly sensitive spectroscopic window into the conformational ensemble, generally permitting the direct quantification of resolvable states. Moreover, straightforward T1- and T2-based relaxation experiments allow for the study of fluctuations within a given state and exchange between states, on timescales spanning nanoseconds to seconds. Conveniently, most biological systems are free of fluorine. Thus, via fluorinated amino acid analogues or thiol-reactive fluorinated tags, F or CF3 reporters can be site specifically incorporated into proteins of interest. In this review, fluorine labeling protocols and 19F NMR experiments will be presented, from the perspective of small molecule NMR (i.e. drug or small molecule interactions with receptors) or macromolecular NMR (i.e. conformational dynamics of receptors and receptor-G-protein complexes).

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.

SPICY: a method for single scan rotating frame relaxometry

T _1ρ is an NMR relaxation mode that is sensitive to low frequency molecular motions, making it an especially valuable tool in biomolecular research. Here, we introduce a new method, SPICY, for measuring T_1ρ relaxation times. In contrast to conventional T_1ρ experiments, in which the sequence is repeated many times to determine the T_1ρ time, the SPICY sequence allows determination of T_1ρ within a single scan, shortening the experiment time remarkably. We demonstrate the method using ¹H T_1ρ relaxation dispersion experiments. Additionally, we combine the sequence with spatial encoding to produce 1D images in a single scan. We show that T_1ρ relaxation times obtained using the single scan approach are in good agreement with those obtained using the traditional experiments.

In-Cell Structural Biology by NMR: The Benefits of the Atomic Scale

In-cell structural biology aims at extracting structural information about proteins or nucleic acids in their native, cellular environment. This emerging field holds great promise and is already providing new facts and outlooks of interest at both fundamental and applied levels. NMR spectroscopy has important contributions on this stage: It brings information on a broad variety of nuclei at the atomic scale, which ensures its great versatility and uniqueness. Here, we detail the methods, the fundamental knowledge, and the applications in biomedical engineering related to in-cell structural biology by NMR. We finally propose a brief overview of the main other techniques in the field (EPR, smFRET, cryo-ET, etc.) to draw some advisable developments for in-cell NMR. In the era of large-scale screenings and deep learning, both accurate and qualitative experimental evidence are as essential as ever to understand the interior life of cells. In-cell structural biology by NMR spectroscopy can generate such a knowledge, and it does so at the atomic scale. This review is meant to deliver comprehensive but accessible information, with advanced technical details and reflections on the methods, the nature of the results, and the future of the field.

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.

Broadband measurement of true transverse relaxation rates in systems with coupled protons: application to the study of conformational exchange

Accurate measurement of transverse relaxation rates in coupled spin systems is important in the study of molecular dynamics, but is severely complicated by the signal modulations caused by scalar couplings in spin echo experiments. The most widely used experiments for measuring transverse relaxation in coupled systems, CPMG and PROJECT, can suppress such modulations, but they also both suppress some relaxation contributions, and average relaxation rates between coupled spins. Here we introduce a new experiment which for the first time allows accurate broadband measurement of transverse relaxation rates of coupled protons, and hence the determination of exchange rate constants in slow exchange from relaxation measurements. The problems encountered with existing methods are illustrated, and the use of the new method is demonstrated for the classic case of hindered amide rotation and for the more challenging problem of exchange between helical enantiomers of a gold(i) complex.

Enzyme dynamics: Looking beyond a single structure

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

Probing conformational transitions towards mutagenic Watson-Crick-like G·T mismatches using off-resonance sugar carbon R1ρ relaxation dispersion

NMR off-resonance R_1ρ relaxation dispersion measurements on base carbon and nitrogen nuclei have revealed that wobble G·T/U mismatches in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-abundance, and mutagenic Watson-Crick-like conformations. As Watson-Crick-like G·T mismatches have base pairing geometries similar to Watson-Crick base pairs, we hypothesized that they would mimic Watson-Crick base pairs with respect to the sugar-backbone conformation as well. Using off-resonance R_1ρ measurements targeting the sugar C3' and C4' nuclei, a structure survey, and molecular dynamics simulations, we show that wobble G·T mismatches adopt sugar-backbone conformations that deviate from the canonical Watson-Crick conformation and that transitions toward tautomeric and anionic Watson-Crick-like G·T mismatches restore the canonical Watson-Crick sugar-backbone. These measurements also reveal kinetic isotope effects for tautomerization in D₂O versus H₂O, which provide experimental evidence in support of a transition state involving proton transfer. The results provide additional evidence in support of mutagenic Watson-Crick-like G·T mismatches, help rule out alternative inverted wobble conformations in the case of anionic G·T^-, and also establish sugar carbons as new non-exchangeable probes of this exchange process.

A Practical Perspective on the Roles of Solution NMR Spectroscopy in Drug Discovery

Solution nuclear magnetic resonance (NMR) spectroscopy is a powerful tool to study structures and dynamics of biomolecules under physiological conditions. As there are numerous NMR-derived methods applicable to probe protein–ligand interactions, NMR has been widely utilized in drug discovery, especially in such steps as hit identification and lead optimization. NMR is frequently used to locate ligand-binding sites on a target protein and to determine ligand binding modes. NMR spectroscopy is also a unique tool in fragment-based drug design (FBDD), as it is able to investigate target-ligand interactions with diverse binding affinities. NMR spectroscopy is able to identify fragments that bind weakly to a target, making it valuable for identifying hits targeting undruggable sites. In this review, we summarize the roles of solution NMR spectroscopy in drug discovery. We describe some methods that are used in identifying fragments, understanding the mechanism of action for a ligand, and monitoring the conformational changes of a target induced by ligand binding. A number of studies have proven that ¹⁹F-NMR is very powerful in screening fragments and detecting protein conformational changes. In-cell NMR will also play important roles in drug discovery by elucidating protein-ligand interactions in living cells.

Exploring exchange processes in proteins by paramagnetic perturbation of NMR spectra

The effect of extrinsic paramagnetic probes on NMR relaxation rates for surface mapping of proteins and other biopolymers is a widely investigated and powerful NMR technique. Here we describe a new application of those probes. It relies on the setting of the relaxation delay to generate magnetization equilibrium and off-equilibrium conditions, in order to tailor the extent of steady state signal recovery with and without the water-soluble nitroxide Tempol. With this approach it is possible to identify signals whose relaxation is affected by exchange processes and, from the relative assignments, to map the protein residues involved in association or conformational interconversion processes on a micro-to-millisecond time scale. This finding is confirmed by the comparison with the results obtained from relaxation dispersion measurements. This simple and convenient method allows preliminary inspection to highlight regions where structural or chemical exchange events are operative, in order to focus on quantitative subsequent determinations by transverse relaxation dispersion experiments or analogous NMR relaxation studies, and/or to gain insights into the predictions of calculations.

Direct evidence for (G)O6···H2-N4(C)+ hydrogen bonding in transient G(syn)-C+ and G(syn)-m5C+ Hoogsteen base pairs in duplex DNA from cytosine amino nitrogen off-resonance R1ρ relaxation dispersion measurements

NMR relaxation dispersion studies have shown that Watson-Crick G-C and A-T base pairs in duplex DNA exist in dynamic equilibrium with their Hoogsteen counterparts. Hoogsteen base pairs form through concurrent rotation of the purine base about the glycosidic bond from an anti to a syn conformation and constriction of the C1'-C1' distance across the base pair by ∼2 Å to allow Hoogsteen type hydrogen bonding. Owing to their unique structure, Hoogsteen base pairs can play important roles in DNA recognition, the accommodation, recognition, and repair of DNA damage, and in DNA replication. NMR relaxation dispersion experiments targeting imino nitrogen and protonated base and sugar carbons have provided insights into many structural features of transient Hoogsteen base pairs, including one of two predicted hydrogen bonds involving (G)N7···H-N3(C)⁺ and (A)N7···H-N3(T). Here, through measurement of cytosine amino (N4) R_1ρ relaxation dispersion, we provide direct evidence for the second (G)O6···H₂-N4(C)⁺ hydrogen bond in G(syn)-C⁺ transient Hoogsteen base pairs. The utility of cytosine N4 R_1ρ relaxation dispersion as a new sensitive probe of transient Hoogsteen base pairs, and cytosine dynamics in general, is further demonstrated by measuring G(syn)-C⁺ Hoogsteen exchange near neutral pH and in the context of the naturally occurring DNA modification 5-methyl cytosine (m⁵C), in DNA samples prepared using chemical synthesis and a ¹⁵N labeled m⁵C phosphoramidite.

Structural mapping of oligomeric intermediates in an amyloid assembly pathway

Transient oligomers are commonly formed in the early stages of amyloid assembly. Determining the structure(s) of these species and defining their role(s) in assembly is key to devising new routes to control disease. Here, using a combination of chemical kinetics, NMR spectroscopy and other biophysical methods, we identify and structurally characterize the oligomers required for amyloid assembly of the protein ΔN6, a truncation variant of human β₂-microglobulin (β₂m) found in amyloid deposits in the joints of patients with dialysis-related amyloidosis. The results reveal an assembly pathway which is initiated by the formation of head-to-head non-toxic dimers and hexamers en route to amyloid fibrils. Comparison with inhibitory dimers shows that precise subunit organization determines amyloid assembly, while dynamics in the C-terminal strand hint to the initiation of cross-β structure formation. The results provide a detailed structural view of early amyloid assembly involving structured species that are not cytotoxic.

Selective Inhibitors of FKBP51 Employ Conformational Selection of Dynamic Invisible States

The recently discovered SAFit class of inhibitors against the Hsp90 co-chaperone FKBP51 show greater than 10 000-fold selectivity over its closely related paralogue FKBP52. However, the mechanism underlying this selectivity remained unknown. By combining NMR spectroscopy, biophysical and computational methods with mutational analysis, we show that the SAFit molecules bind to a transient pocket in FKBP51. This represents a weakly populated conformation resembling the inhibitor-bound state of FKBP51, suggesting conformational selection rather than induced fit as the major binding mechanism. The inhibitor-bound conformation of FKBP51 is stabilized by an allosteric network of residues located away from the inhibitor-binding site. These residues stabilize the Phe67 side chain in a dynamic outward conformation and are distinct in FKBP52, thus rationalizing the basis for the selectivity of SAFit inhibitors. Our results represent a paradigm for the selective inhibition of transient binding pockets.

An NMR View of Protein Dynamics in Health and Disease

Biological molecules are often highly dynamic, and this flexibility can be critical for function. The large range of sampled timescales and the fact that many of the conformers that are continually explored are only transiently formed and sparsely populated challenge current biophysical approaches. Solution nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for characterizing biomolecular dynamics in detail, even in cases where excursions involve short-lived states. Here, we briefly review a number of NMR experiments for studies of biomolecular dynamics on the microsecond-to-second timescale and focus on applications to protein and nucleic acid systems that clearly illustrate the functional relevance of motion in both health and disease.

NMR analysis of GPCR conformational landscapes and dynamics

Understanding the signal transduction mechanism mediated by the G Protein-Coupled Receptors (GPCRs) in eukaryote cells represents one of the main issues in modern biology. At the molecular level, various biophysical approaches have provided important insights on the functional plasticity of these complex allosteric machines. In this context, X-ray crystal structures published during the last decade represent a major breakthrough in GPCR structural biology, delivering important information on the activation process of these receptors through the description of the three-dimensional organization of their active and inactive states. In complement to crystals and cryo-electronic microscopy structures, information on the probability of existence of different GPCR conformations and the dynamic barriers separating those structural sub-states is required to better understand GPCR function. Among the panel of techniques available, nuclear magnetic resonance (NMR) spectroscopy represents a powerful tool to characterize both conformational landscapes and dynamics. Here, we will outline the potential of NMR to address such biological questions, and we will illustrate the functional insights that NMR has brought in the field of GPCRs in the recent years.

Characterizing Protein-Protein Interactions Using Solution NMR Spectroscopy

In this chapter, we describe how NMR chemical shift titrations can be used to study the interaction between two proteins with emphasis on mapping the interface of the complex and determining the binding affinity from a quantitative analysis of the experimental data. In particular, we discuss the appearance of NMR spectra in different chemical exchange regimes (fast, intermediate, and slow) and how these regimes affect NMR data analysis.

Microsecond Protein Dynamics from Combined Bloch-McConnell and Near-Rotary-Resonance R1p Relaxation-Dispersion MAS NMR

Studying protein dynamics on microsecond-to-millisecond (μs-ms) time scales can provide important insight into protein function. In magic-angle-spinning (MAS) NMR, μs dynamics can be visualized by R 1 ρ rotating-frame relaxation dispersion experiments in different regimes of radio-frequency field strengths: at low RF field strength, isotropic-chemical-shift fluctuation leads to "Bloch-McConnell-type" relaxation dispersion, while when the RF field approaches rotary resonance conditions bond angle fluctuations manifest as increased R 1 ρ rate constants ("Near-Rotary-Resonance Relaxation Dispersion", NERRD). Here we explore the joint analysis of both regimes to gain comprehensive insight into motion in terms of geometric amplitudes, chemical-shift changes, populations and exchange kinetics. We use a numerical simulation procedure to illustrate these effects and the potential of extracting exchange parameters, and apply the methodology to the study of a previously described conformational exchange process in microcrystalline ubiquitin.

Overview of Relaxation Dispersion NMR Spectroscopy to Study Protein Dynamics and Protein-Ligand Interactions

Proteins and nucleic acids are central to all biological processes. NMR spectroscopy has proven to be excellent for studying the dynamics of these macromolecules over various timescales. Relaxation rates and heteronuclear nuclear Overhauser-effect values can resolve motion on pico- to nanosecond timescales, residual dipolar couplings provide information on submicro- to millisecond timescales, and even slower dynamics over seconds to hours can be resolved by hydrogen-exchange experiments. Relaxation dispersion experiments are especially valuable because they resolve motion on micro- to millisecond timescales, encompassing biomolecular motions associated with ligand binding, enzymatic catalysis, and domain-domain opening. These experiments provide structural, kinetic, and thermodynamic information on "invisible" excited conformational states. Relaxation dispersion can be applied not only to single biomolecules but also to protein-ligand complexes to study the kinetics and thermodynamics of association and dissociation. We review recent developments in relaxation dispersion methodology, outline the R_1ρ relaxation dispersion experiment, and discuss application to biomolecular interactions. © 2018 by John Wiley & Sons, Inc.

Characterization of Dynamic IDP Complexes by NMR Spectroscopy

NMR spectroscopy has proven to be a key method for studying intrinsically disordered proteins (IDPs). Nonetheless, traditional NMR methods developed for solving structures of ordered protein complexes are insufficient for the full characterization of dynamic IDP complexes, where the energy landscape is broader and more rugged. Furthermore, due to their high sensitivity to environmental changes, NMR studies of IDP complexes must be conducted with extra care and the observed NMR parameters thoroughly evaluated to enable disentanglement of binding events from ensemble distribution changes. In this chapter, written for the non-NMR expert, we start out by outlining sample preparation for IDP complexes, guide through the recording and evaluation of diagnostic ¹H,¹⁵N-HSQC spectra, and delineate more sophisticated NMR strategies to follow for the particular type of complex. The most relevant experiments are then described in terms of aims, needs, pitfalls, analysis, and expected outcomes, with references to recent examples.

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 ¹H 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 ¹H 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 B_o, 2B_o and 3B_o, 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.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

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]).

Computational studies of G protein-coupled receptor complexes: Structure and dynamics

G protein-coupled receptors (GPCRs) are ubiquitously expressed transmembrane proteins associated with a wide range of diseases such as Alzheimer's, Parkinson, schizophrenia, and also implicated in in several abnormal heart conditions. As such, this family of receptors is regarded as excellent drug targets. However, due to the high number of intracellular signaling partners, these receptors have a complex interaction networks and it becomes challenging to modulate their function. Experimentally determined structures give detailed information on the salient structural properties of these signaling complexes but they are far away from providing mechanistic insights into the underlying process. This chapter presents some of the computational tools, namely molecular dynamics, molecular docking, and molecular modeling and related analyses methods that have been used to complement experimental findings.

Characterization of Chemical Exchange Using Relaxation Dispersion of Hyperpolarized Nuclear Spins

Chemical exchange phenomena are ubiquitous in macromolecules, which undergo conformational change or ligand complexation. NMR relaxation dispersion (RD) spectroscopy based on a Carr-Purcell-Meiboom-Gill pulse sequence is widely applied to identify the exchange and measure the lifetime of intermediate states on the millisecond time scale. Advances in hyperpolarization methods improve the applicability of NMR spectroscopy when rapid acquisitions or low concentrations are required, through an increase in signal strength by several orders of magnitude. Here, we demonstrate the measurement of chemical exchange from a single aliquot of a ligand hyperpolarized by dissolution dynamic nuclear polarization (D-DNP). Transverse relaxation rates are measured simultaneously at different pulsing delays by dual-channel ¹⁹F NMR spectroscopy. This two-point measurement is shown to allow the determination of the exchange term in the relaxation rate expression. For the ligand 4-(trifluoromethyl)benzene-1-carboximidamide binding to the protein trypsin, the exchange term is found to be equal within error limits in neutral and acidic environments from D-DNP NMR spectroscopy, corresponding to a pre-equilibrium of trypsin deprotonation. This finding illustrates the capability for determination of binding mechanisms using D-DNP RD. Taking advantage of hyperpolarization, the ligand concentration in the exchange measurements can reach on the order of tens of μM and protein concentration can be below 1 μM, i.e., conditions typically accessible in drug discovery.

Unraveling the meaning of chemical shifts in protein NMR

Chemical shifts are among the most informative parameters in protein NMR. They provide wealth of information about protein secondary and tertiary structure, protein flexibility, and protein-ligand binding. In this report, we review the progress in interpreting and utilizing protein chemical shifts that has occurred over the past 25years, with a particular focus on the large body of work arising from our group and other Canadian NMR laboratories. More specifically, this review focuses on describing, assessing, and providing some historical context for various chemical shift-based methods to: (1) determine protein secondary and super-secondary structure; (2) derive protein torsion angles; (3) assess protein flexibility; (4) predict residue accessible surface area; (5) refine 3D protein structures; (6) determine 3D protein structures and (7) characterize intrinsically disordered proteins. This review also briefly covers some of the methods that we previously developed to predict chemical shifts from 3D protein structures and/or protein sequence data. It is hoped that this review will help to increase awareness of the considerable utility of NMR chemical shifts in structural biology and facilitate more widespread adoption of chemical-shift based methods by the NMR spectroscopists, structural biologists, protein biophysicists, and biochemists worldwide. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Protein folding by NMR

Protein folding is a highly complex process proceeding through a number of disordered and partially folded nonnative states with various degrees of structural organization. These transiently and sparsely populated species on the protein folding energy landscape play crucial roles in driving folding toward the native conformation, yet some of these nonnative states may also serve as precursors for protein misfolding and aggregation associated with a range of devastating diseases, including neuro-degeneration, diabetes and cancer. Therefore, in vivo protein folding is often reshaped co- and post-translationally through interactions with the ribosome, molecular chaperones and/or other cellular components. Owing to developments in instrumentation and methodology, solution NMR spectroscopy has emerged as the central experimental approach for the detailed characterization of the complex protein folding processes in vitro and in vivo. NMR relaxation dispersion and saturation transfer methods provide the means for a detailed characterization of protein folding kinetics and thermodynamics under native-like conditions, as well as modeling high-resolution structures of weakly populated short-lived conformational states on the protein folding energy landscape. Continuing development of isotope labeling strategies and NMR methods to probe high molecular weight protein assemblies, along with advances of in-cell NMR, have recently allowed protein folding to be studied in the context of ribosome-nascent chain complexes and molecular chaperones, and even inside living cells. Here we review solution NMR approaches to investigate the protein folding energy landscape, and discuss selected applications of NMR methodology to studying protein folding in vitro and in vivo. Together, these examples highlight a vast potential of solution NMR in providing atomistic insights into molecular mechanisms of protein folding and homeostasis in health and disease.

(S)Pinning down protein interactions by NMR

Protein molecules are highly diverse communication platforms and their interaction repertoire stretches from atoms over small molecules such as sugars and lipids to macromolecules. An important route to understanding molecular communication is to quantitatively describe their interactions. These types of analyses determine the amounts and proportions of individual constituents that participate in a reaction as well as their rates of reactions and their thermodynamics. Although many different methods are available, there is currently no single method able to quantitatively capture and describe all types of protein reactions, which can span orders of magnitudes in affinities, reaction rates, and lifetimes of states. As the more versatile technique, solution NMR spectroscopy offers a remarkable catalogue of methods that can be successfully applied to the quantitative as well as qualitative descriptions of protein interactions. In this review we provide an easy-access approach to NMR for the non-NMR specialist and describe how and when solution state NMR spectroscopy is the method of choice for addressing protein ligand interaction. We describe very briefly the theoretical background and illustrate simple protein-ligand interactions as well as typical strategies for measuring binding constants using NMR spectroscopy. Finally, this review provides examples of caveats of the method as well as the options to improve the outcome of an NMR analysis of a protein interaction reaction.

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Determining the structures, kinetics, thermodynamics and mechanisms that underlie conformational exchange processes in proteins remains extremely difficult. Only in favourable cases is it possible to provide atomic-level descriptions of sparsely populated and transiently formed alternative conformations. Here we benchmark the ability of enhanced-sampling molecular dynamics simulations to determine the free energy landscape of the L99A cavity mutant of T4 lysozyme. We find that the simulations capture key properties previously measured by NMR relaxation dispersion methods including the structure of a minor conformation, the kinetics and thermodynamics of conformational exchange, and the effect of mutations. We discover a new tunnel that involves the transient exposure towards the solvent of an internal cavity, and show it to be relevant for ligand escape. Together, our results provide a comprehensive view of the structural landscape of a protein, and point forward to studies of conformational exchange in systems that are less characterized experimentally.

DOI:http://dx.doi.org/10.7554/eLife.17505.001

Proteins are the workhorses of cells, where they perform a wide range of roles. To do so, they adopt specific three-dimensional structures that enable them to interact with other molecules as necessary. Often a protein needs to be able to shift between different states with distinct structures as it goes about its job. To fully understand how a protein works, it is important to be able to characterize these different structures and how the protein changes between them.

Many of the experimental techniques used to study protein structure rely on isolating the individual structural forms of a protein. Since many structures only exist briefly, this can be very difficult. To complement experimental results, computer simulations allow researchers to model how atoms behave within a molecule. However, a number of factors limit how well these models represent what happens experimentally, such as the accuracy of the physical description used for the modeling.

Wang et al. set out to test and benchmark how well computer simulations could model changes in structure for a protein called T4 lysozyme, which has been studied extensively using experimental techniques. T4 lysozyme exists in two different states that have distinct structures. By comparing existing detailed experimental measurements with the results of their simulations, Wang et al. found that the simulations could capture key aspects of how T4 lysozyme changes its shape.

The simulations described the structure of the protein in both states and accurately determined the relative proportion of molecules that are found in each state. They could also determine how long it takes for a molecule to change its shape from one state to the other. The findings allowed Wang et al. to describe in fine detail – down to the level of individual atoms – how the protein changes its shape and how mutations in the protein affect its ability to do so. A key question for future studies is whether these insights can be extended to other proteins that are less well characterized experimentally than T4 lysozyme.

DOI:http://dx.doi.org/10.7554/eLife.17505.002

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

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.

Mechanisms of amyloid formation revealed by solution NMR

Amyloid fibrils are proteinaceous elongated aggregates involved in more than fifty human diseases. Recent advances in electron microscopy and solid state NMR have allowed the characterization of fibril structures to different extents of refinement. However, structural details about the mechanism of fibril formation remain relatively poorly defined. This is mainly due to the complex, heterogeneous and transient nature of the species responsible for assembly; properties that make them difficult to detect and characterize in structural detail using biophysical techniques. The ability of solution NMR spectroscopy to investigate exchange between multiple protein states, to characterize transient and low-population species, and to study high molecular weight assemblies, render NMR an invaluable technique for studies of amyloid assembly. In this article we review state-of-the-art solution NMR methods for investigations of: (a) protein dynamics that lead to the formation of aggregation-prone species; (b) amyloidogenic intrinsically disordered proteins; and (c) protein-protein interactions on pathway to fibril formation. Together, these topics highlight the power and potential of NMR to provide atomic level information about the molecular mechanisms of one of the most fascinating problems in structural biology.

The role of slow and fast protein motions in allosteric interactions

Allostery is fundamentally thermodynamic in nature. Long-range communication in proteins may be mediated not only by changes in the mean conformation with enthalpic contribution but also by changes in dynamic fluctuations with entropic contribution. The important role of protein motions in mediating allosteric interactions has been established by NMR spectroscopy. By using CAP as a model system, we have shown how changes in protein structure and internal dynamics can allosterically regulate protein function and activity. The results indicate that changes in conformational entropy can give rise to binding enhancement, binding inhibition, or have no effect in the expected affinity, depending on the magnitude and sign of enthalpy-entropy compensation. Moreover, allosteric interactions can be regulated by the modulation a low-populated conformation states that serve as on-pathway intermediates for ligand binding. Taken together, the interplay between fast internal motions, which are intimately related to conformational entropy, and slow internal motions, which are related to poorly populated conformational states, can regulate protein activity in a way that cannot be predicted on the basis of the protein's ground-state structure.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

A robust approach to correct for pronounced errors in temperature measurements by controlling radiation damping feedback fields in solution NMR

Accurate temperature measurement is a requisite for obtaining reliable thermodynamic and kinetic information in all NMR experiments. A widely used method to calibrate sample temperature depends on a secondary standard with temperature-dependent chemical shifts to report the true sample temperature, such as the hydroxyl proton in neat methanol or neat ethylene glycol. The temperature-dependent chemical shift of the hydroxyl protons arises from the sensitivity of the hydrogen-bond network to small changes in temperature. The frequency separation between the alkyl and the hydroxyl protons are then converted to sample temperature. Temperature measurements by this method, however, have been reported to be inconsistent and incorrect in modern NMR, particularly for spectrometers equipped with cryogenically-cooled probes. Such errors make it difficult or even impossible to study chemical exchange and molecular dynamics or to compare data acquired on different instruments, as is frequently done in biomolecular NMR. In this work, we identify the physical origins for such errors to be unequal amount of dynamical frequency shifts on the alkyl and the hydroxyl protons induced by strong radiation damping (RD) feedback fields. Common methods used to circumvent RD may not suppress such errors. A simple, easy-to-implement solution was demonstrated that neutralizes the RD effect on the frequency separation by a "selective crushing recovery" pulse sequence to equalize the transverse magnetization of both spin species. Experiments using cryoprobes at 500 MHz and 800 MHz demonstrated that this approach can effectively reduce the errors in temperature measurements from about ±4.0 K to within ±0.4 K in general.

GPCRs: Lipid-Dependent Membrane Receptors That Act as Drug Targets

G protein-coupled receptors (GPCRs) are the largest class of molecules involved in signal transduction across cell membranes and represent major targets in the development of novel drug candidates in all clinical areas. Although there have been some recent leads, structural information on GPCRs is relatively rare due to the difficulty associated with crystallization. A specific reason for this is the intrinsic flexibility displayed by GPCRs, which is necessary for their functional diversity. Since GPCRs are integral membrane proteins, interaction of membrane lipids with them constitutes an important area of research in GPCR biology. In particular, membrane cholesterol has been reported to have a modulatory role in the function of a number of GPCRs. The role of membrane cholesterol in GPCR function is discussed with specific example of the serotonin1A receptor. Recent results show that GPCRs are characterized with structural motifs that preferentially associate with cholesterol. An emerging and important concept is oligomerization of GPCRs and its role in GPCR function and signaling. The role of membrane cholesterol in GPCR oligomerization is highlighted. Future research in GPCR biology would offer novel insight in basic biology and provide new avenues for drug discovery.

In-cell NMR and EPR spectroscopy of biomacromolecules

The dream of cell biologists is to be able to watch biological macromolecules perform their duties in the intracellular environment of live cells. Ideally, the observation of both the location and the conformation of these macromolecules with biophysical techniques is desired. The development of many fluorescence techniques, including superresolution fluorescence microscopy, has significantly enhanced our ability to spot proteins and other molecules in the crowded cellular environment. However, the observation of their structure and conformational changes while they attend their business is still very challenging. In principle, NMR and EPR spectroscopy can be used to investigate the conformation and dynamics of biological macromolecules in living cells. The development of in-cell magnetic resonance techniques has demonstrated the feasibility of this approach. Herein we review the different techniques with a focus on liquid-state in-cell NMR spectroscopy, provide an overview of applications, and discuss the challenges that lie ahead.

Measuring hydrogen exchange rates in invisible protein excited states

Hydrogen exchange rates have become a valuable probe for studying the relationship between dynamics and structure and for dissecting the mechanism by which proteins fold to their native conformation. Typically measured rates correspond to averages over all protein states from which hydrogen exchange can occur. Here we describe a new NMR experiment based on chemical exchange saturation transfer that provides an avenue for obtaining uncontaminated, per-residue amide hydrogen exchange rates for interconverting native and invisible states so long as they can be separated on the basis of distinct (15)N chemical shifts. The approach is applied to the folding reaction of the Fyn SH3 domain that exchanges between a highly populated, NMR-visible native state and a conformationally excited, NMR-invisible state, corresponding to the unfolded ensemble. Excellent agreement between experimentally derived hydrogen exchange rates of the excited state at a pair of pHs is obtained, taking into account the expected dependence of exchange on pH. Extracted rates for the unfolded ensemble have been used to test hydrogen exchange predictions based on the primary protein sequence that are used in many analyses of solvent exchange rates, with a Pearson correlation coefficient of 0.84 obtained.

Measurement of histidine pKa values and tautomer populations in invisible protein states

The histidine imidazole side chain plays a critical role in protein function and stability. Its importance for catalysis is underscored by the fact that histidines are localized to active sites in ∼ 50% of all enzymes. NMR spectroscopy has become an important tool for studies of histidine side chains through the measurement of site-specific pK(a)s and tautomer populations. To date, such studies have been confined to observable protein ground states; however, a complete understanding of the role of histidine electrostatics in protein function and stability requires that similar investigations be extended to rare, transiently formed conformers that populate the energy landscape, yet are often "invisible" in standard NMR spectra. Here we present NMR experiments and a simple strategy for studies of such conformationally excited states based on measurement of histidine (13)Cγ, (13)Cδ2 chemical shifts and (1)Hε-(13)Cε one-bond scalar couplings. The methodology is first validated and then used to obtain pKa values and tautomer distributions for histidine residues of an invisible on-pathway folding intermediate of the colicin E7 immunity protein. Our results imply that the side chains of H40 and H47 are exposed in the intermediate state and undergo significant conformational rearrangements during folding to the native structure. Further, the pKa values explain the pH-dependent stability differences between native and intermediate states over the pH range 5.5-6.5 and they suggest that imidazole deprotonation is not a barrier to the folding of this protein.