Probing chemical shifts of invisible states of proteins with relaxation dispersion NMR spectroscopy: how well can we do?

Transient interdomain interactions in free USP14 shape its conformational ensemble

The deubiquitinase (DUB) ubiquitin-specific protease 14 (USP14) is a dual domain protein that plays a regulatory role in proteasomal degradation and has been identified as a promising therapeutic target. USP14 comprises a conserved USP domain and a ubiquitin-like (Ubl) domain separated by a 25-residue linker. The enzyme activity of USP14 is autoinhibited in solution, but is enhanced when bound to the proteasome, where the Ubl and USP domains of USP14 bind to the Rpn1 and Rpt1/Rpt2 units, respectively. No structure of full-length USP14 in the absence of proteasome has yet been presented, however, earlier work has described how transient interactions between Ubl and USP domains in USP4 and USP7 regulate DUB activity. To better understand the roles of the Ubl and USP domains in USP14, we studied the Ubl domain alone and in full-length USP14 by nuclear magnetic resonance spectroscopy and used small angle x-ray scattering and molecular modeling to visualize the entire USP14 protein ensemble. Jointly, our results show how transient interdomain interactions between the Ubl and USP domains of USP14 predispose its conformational ensemble for proteasome binding, which may have functional implications for proteasome regulation and may be exploited in the design of future USP14 inhibitors.

Biomolecular NMR spectroscopy in the era of artificial intelligence

Biomolecular nuclear magnetic resonance (NMR) spectroscopy and artificial intelligence (AI) have a burgeoning synergy. Deep learning-based structural predictors have forever changed structural biology, yet these tools currently face limitations in accurately characterizing protein dynamics, allostery, and conformational heterogeneity. We begin by highlighting the unique abilities of biomolecular NMR spectroscopy to complement AI-based structural predictions toward addressing these knowledge gaps. We then highlight the direct integration of deep learning approaches into biomolecular NMR methods. AI-based tools can dramatically improve the acquisition and analysis of NMR spectra, enhancing the accuracy and reliability of NMR measurements, thus streamlining experimental processes. Additionally, deep learning enables the development of novel types of NMR experiments that were previously unattainable, expanding the scope and potential of biomolecular NMR spectroscopy. Ultimately, a combination of AI and NMR promises to further revolutionize structural biology on several levels, advance our understanding of complex biomolecular systems, and accelerate drug discovery efforts.

Intrinsic structural dynamics dictate enzymatic activity and inhibition

Enzymes are known to sample various conformations, many of which are critical for their biological function. However, structural characterizations of enzymes predominantly focus on the most populated conformation. As a result, single-point mutations often produce structures that are similar or essentially identical to those of the wild-type enzyme despite large changes in enzymatic activity. Here, we show for mutants of a histone deacetylase enzyme (HDAC8) that reduced enzymatic activities, reduced inhibitor affinities, and reduced residence times are all captured by the rate constants between intrinsically sampled conformations that, in turn, can be obtained independently by solution NMR spectroscopy. Thus, for the HDAC8 enzyme, the dynamic sampling of conformations dictates both enzymatic activity and inhibitor potency. Our analysis also dissects the functional role of the conformations sampled, where specific conformations distinct from those in available structures are responsible for substrate and inhibitor binding, catalysis, and product dissociation. Precise structures alone often do not adequately explain the effect of missense mutations on enzymatic activity and drug potency. Our findings not only assign functional roles to several conformational states of HDAC8 but they also underscore the paramount role of dynamics, which will have general implications for characterizing missense mutations and designing inhibitors.

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.

CPMG pulse sequence for relaxation dispersion that cancels artifacts independently of spin states

The relaxation dispersion (rd) of nuclear magnetic resonance provides thermodynamic and kinetic information regarding molecules for which the conformations are exchanging in equilibrium. Experiments have often been implemented with Carr-Purcell-Meiboom-Gill (CPMG) pulse sequences for heteronuclear S-spin in SI and SI3 spin systems. One of the most common CPMG sequences contains a sequence called a P-element in the middle to average the different relaxation rates of anti-phase and in-phase coherences; however, its drawback is that the artifacts that can be compensated for are only those in one of the two S-spin doublet magnetization components, for example, SIα or SIβ in an SI spin system. Thus, when the CPMG sequence is followed by a normal heteronuclear single-quantum correlation (HSQC) sequence, the detected signal will also include the other component with accumulated artifacts. To overcome this issue, we developed a new pulse sequence (AFTAC) that can suppress artifacts in both the magnetization components. Its effectiveness was demonstrated by simulations and actual measurements targeting the methyl groups of dimethylsulfoxide and N, N-dimethyltrichloroacetamide. The results demonstrated that the AFTAC sequence sufficiently suppressed the artifacts, despite being followed by an HSQC sequence that detects both components. AFTAC is particularly suitable for the rd measurements of small molecules, which are usually not deuterated and are not subject to transverse relaxation optimized spectroscopy (TROSY) sequences. AFTAC does not require 1H continuous wave irradiation for I-spin decoupling, which is necessary for certain CPMG methods that maintain S-spin in-phase coherence during the CPMG period (Tcpmg). Therefore, AFTAC places less burden on the probe, even with a long Tcpmg. Furthermore, the AFTAC method achieves a similar artifact suppression quality not only in SI but also in SI2 and SI3 spin systems.

Towards autonomous analysis of chemical exchange saturation transfer experiments using deep neural networks

Macromolecules often exchange between functional states on timescales that can be accessed with NMR spectroscopy and many NMR tools have been developed to characterise the kinetics and thermodynamics of the exchange processes, as well as the structure of the conformers that are involved. However, analysis of the NMR data that report on exchanging macromolecules often hinges on complex least-squares fitting procedures as well as human experience and intuition, which, in some cases, limits the widespread use of the methods. The applications of deep neural networks (DNNs) and artificial intelligence have increased significantly in the sciences, and recently, specifically, within the field of biomolecular NMR, where DNNs are now available for tasks such as the reconstruction of sparsely sampled spectra, peak picking, and virtual decoupling. Here we present a DNN for the analysis of chemical exchange saturation transfer (CEST) data reporting on two- or three-site chemical exchange involving sparse state lifetimes of between approximately 3-60 ms, the range most frequently observed via experiment. The work presented here focuses on the ¹H CEST class of methods that are further complicated, in relation to applications to other nuclei, by anti-phase features. The developed DNNs accurately predict the chemical shifts of nuclei in the exchanging species directly from anti-phase ¹H^N CEST profiles, along with an uncertainty associated with the predictions. The performance of the DNN was quantitatively assessed using both synthetic and experimental anti-phase CEST profiles. The assessments show that the DNN accurately determines chemical shifts and their associated uncertainties. The DNNs developed here do not contain any parameters for the end-user to adjust and the method therefore allows for autonomous analysis of complex NMR data that report on conformational exchange.

NMR Relaxation Dispersion Methods for the Structural and Dynamic Analysis of Quickly Interconverting, Low-Populated Conformational Substates

Most biomolecular processes involve proteins shuttling among different conformational states, particularly from highly populated ground states to the lowly populated excited states that determine the interconversion rates and biological function, and which are invisible to most structural biology techniques. These structural transitions are rare and relatively fast: happen in the millisecond-microsecond timescale (ms-μs). NMR spectroscopy can access these timescales via relaxation dispersion techniques (RD-NMR). The exchange parameters extracted from RD-NMR experiments provide pivotal information on these otherwise invisible states that reports on key properties of the high free energy, reactive regions of the protein's energy landscape, including the mechanisms of folding/unfolding and of the interconversion between active and inactive states. Here, we describe a simple, step-by-step protocol to carry out RD-NMR experiments on proteins to detect the existence of such conformational substates and characterize their structural properties (chemical shifts).

Two β-Lactamase Variants with Reduced Clavulanic Acid Inhibition Display Different Millisecond Dynamics

The β-lactamase of Mycobacterium tuberculosis, BlaC, is susceptible to inhibition by clavulanic acid. The ability of this enzyme to escape inhibition through mutation was probed using error-prone PCR combined with functional screening in Escherichia coli. The variant that was found to confer the most inhibitor resistance, K234R, as well as variant G132N that was found previously were characterized using X-ray crystallography and nuclear magnetic resonance (NMR) relaxation experiments to probe structural and dynamic properties. The G132N mutant exists in solution in two almost equally populated conformations that exchange with a rate of ca. 88 s-1. The conformational change affects a broad region of the enzyme. The crystal structure reveals that the Asn132 side chain forces the peptide bond between Ser104 and Ile105 in a cis-conformation. The crystal structure suggests multiple conformations for several side chains (e.g., Ser104 and Ser130) and a short loop (positions 214 to 216). In the K234R mutant, the active-site dynamics are significantly diminished with respect to the wild-type enzyme. These results show that multiple evolutionary routes are available to increase inhibitor resistance in BlaC and that active-site dynamics on the millisecond time scale are not required for catalytic function.

Quantitative Exchange NMR-Based Analysis of Huntingtin-SH3 Interactions Suggests an Allosteric Mechanism of Inhibition of Huntingtin Aggregation

Huntingtin polypeptides (htt^ex1), encoded by exon 1 of the htt gene and containing an expanded polyglutamine tract, form fibrils that accumulate within neuronal inclusion bodies, resulting in the fatal neurodegenerative condition known as Huntington's disease. Htt^ex1 comprises three regions: a 16-residue N-terminal amphiphilic domain (NT), a polyglutamine tract of variable length (Q_n), and a polyproline-rich domain containing two polyproline tracts. The NT region of htt^ex1 undergoes prenucleation transient oligomerization on the sub-millisecond time scale, resulting in a productive tetramer that promotes self-association and nucleation of the polyglutamine tracts. Here we show that binding of Fyn SH3, a small intracellular proline-binding domain, to the first polyproline tract of htt^ex1Q₃₅ inhibits fibril formation by both NMR and a thioflavin T fluorescence assay. The interaction of Fyn SH3 with htt^ex1Q₇ was investigated using NMR experiments designed to probe kinetics and equilibria at atomic resolution, including relaxation dispersion, and concentration-dependent exchange-induced chemical shifts and transverse relaxation in the rotating frame. Sub-millisecond exchange between four species is demonstrated: two major states comprising free (P) and SH3-bound (PL) monomeric htt^ex1Q₇, and two sparsely populated dimers in which either both subunits (P₂L₂) or only a single subunit (P₂L) is bound to SH3. Binding of SH3 increases the helical propensity of the NT domain, resulting in a 25-fold stabilization of the P₂L₂ dimer relative to the unliganded P₂ dimer. The P₂L₂ dimer, in contrast to P₂, does not undergo any detectable oligomerization to a tetramer, thereby explaining the allosteric inhibition of htt^ex1 fibril formation by Fyn SH3.

A weakened interface in the P182L variant of HSP27 associated with severe Charcot-Marie-Tooth neuropathy causes aberrant binding to interacting proteins

HSP27 is a human molecular chaperone that forms large, dynamic oligomers and functions in many aspects of cellular homeostasis. Mutations in HSP27 cause Charcot-Marie-Tooth (CMT) disease, the most common inherited disorder of the peripheral nervous system. A particularly severe form of CMT disease is triggered by the P182L mutation in the highly conserved IxI/V motif of the disordered C-terminal region, which interacts weakly with the structured core domain of HSP27. Here, we observed that the P182L mutation disrupts the chaperone activity and significantly increases the size of HSP27 oligomers formed in vivo, including in motor neurons differentiated from CMT patient-derived stem cells. Using NMR spectroscopy, we determined that the P182L mutation decreases the affinity of the HSP27 IxI/V motif for its own core domain, leaving this binding site more accessible for other IxI/V-containing proteins. We identified multiple IxI/V-bearing proteins that bind with higher affinity to the P182L variant due to the increased availability of the IxI/V-binding site. Our results provide a mechanistic basis for the impact of the P182L mutation on HSP27 and suggest that the IxI/V motif plays an important, regulatory role in modulating protein-protein interactions.

Probing the Interaction of Huntingtin Exon-1 Polypeptides with the Chaperonin Nanomachine GroEL

Huntington's disease arises from polyQ expansion within the exon-1 region of huntingtin (htt^ex1 ), resulting in an aggregation-prone protein that accumulates in neuronal inclusion bodies. We investigate the interaction of various htt^ex1 constructs with the bacterial analog (GroEL) of the human chaperonin Hsp60. Using fluorescence spectroscopy and electron and atomic force microscopy, we show that GroEL inhibits fibril formation. The binding kinetics of htt^ex1 constructs with intact GroEL and a mini-chaperone comprising the apical domain is characterized by relaxation-based NMR measurements. The lifetimes of the complexes range from 100 to 400 μs with equilibrium dissociation constants (K_D ) of ∼1-2 mM. The binding interface is formed by the N-terminal amphiphilic region of htt^ex1 (which adopts a partially helical conformation) and the H and I helices of the GroEL apical domain. Sequestration of monomeric htt^ex1 by GroEL likely increases the critical concentration required for fibrillization.

Backbone chemical shift assignments for the SARS-CoV-2 non-structural protein Nsp9: intermediate (ms - μs) dynamics in the C-terminal helix at the dimer interface

The Betacoronavirus SARS-CoV-2 non-structural protein Nsp9 is a 113-residue protein that is essential for viral replication, and consequently, a potential target for the development of therapeutics against COVID19 infections. To capture insights into the dynamics of the protein's backbone in solution and accelerate the identification and mapping of ligand-binding surfaces through chemical shift perturbation studies, the backbone ¹H, ¹³C, and ¹⁵N NMR chemical shifts for Nsp9 have been extensively assigned. These assignments were assisted by the preparation of an ~ 70% deuterated sample and residue-specific, ¹⁵N-labelled samples (V, L, M, F, and K). A major feature of the assignments was the "missing" amide resonances for N96-L106 in the ¹H-¹⁵N HSQC spectrum, a region that comprises almost the complete C-terminal α-helix that forms a major part of the homodimer interface in the crystal structure of SARS-CoV-2 Nsp9, suggesting this region either undergoes intermediate motion in the ms to μs timescale and/or is heterogenous. These "missing" amide resonances do not unambiguously appear in the ¹H-¹⁵N HSQC spectrum of SARS-CoV-2 Nsp9 collected at a concentration of 0.0007 mM. At this concentration, at the detection limit, native mass spectrometry indicates the protein is exclusively in the monomeric state, suggesting the intermediate motion in the C-terminal of Nsp9 may be due to intramolecular dynamics. Perhaps this intermediate ms to μs timescale dynamics is the physical basis for a previously suggested "fluidity" of the C-terminal helix that may be responsible for homophilic (Nsp9-Nsp9) and postulated heterophilic (Nsp9-Unknown) protein-protein interactions.

RING NMR dynamics: software for analysis of multiple NMR relaxation experiments

Molecular motions are fundamental to the existence of life, and NMR spectroscopy remains one of the most useful and powerful methods to measure their rates and molecular characteristics. Multiple experimental methods are available for measuring the NMR relaxation properties and these can require different methods for extracting model parameters. We present here a new software application, RING NMR Dynamics, that is designed to support analysis of multiple relaxation types. The initial release of RING NMR Dynamics supports the analysis of exponential decay experiments such as T₁ and T₂, as well as CEST and R₂ and R_1ρ relaxation dispersion. The software runs on multiple operating systems in both a command line mode and a user-friendly GUI that allows visualizing and simulating relaxation data. Interaction with another program, NMRFx Analyst, allows drilling down from the derived relaxation parameters to the raw spectral data.

Fusidic acid resistance through changes in the dynamics of the drug target

Antibiotic resistance in clinically important bacteria can be mediated by target protection mechanisms, whereby a protein binds to the drug target and protects it from the inhibitory effects of the antibiotic. The most prevalent source of clinical resistance to the antibiotic fusidic acid (FA) is expression of the FusB family of proteins that bind to the drug target (Elongation factor G [EF-G]) and promote dissociation of EF-G from FA-stalled ribosome complexes. FusB binding causes changes in both the structure and conformational flexibility of EF-G, but which of these changes drives FA resistance was not understood. We present here detailed characterization of changes in the conformational flexibility of EF-G in response to FusB binding and show that these changes are responsible for conferring FA resistance. Binding of FusB to EF-G causes a significant change in the dynamics of domain III of EF-GC3 that leads to an increase in a minor, more disordered state of EF-G domain III. This is sufficient to overcome the steric block of transmission of conformational changes within EF-G by which FA prevents release of EF-G from the ribosome. This study has identified an antibiotic resistance mechanism mediated by allosteric effects on the dynamics of the drug target.

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 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 ¹H, ¹⁵N and ¹³C 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 ¹H CEST experiment to study exchange at glycine (Gly) ¹H^α sites as the existing ¹H CEST experiments that have been designed to study dynamics in amide ¹H-¹⁵N spin systems and methyl ¹³CH₃ groups with three equivalent protons while suppressing ¹H-¹H NOE induced dips are not suitable for studying exchange in methylene ¹³CH₂ groups with inequivalent protons. Here a Gly ¹H^α CEST experiment to obtain the minor state Gly ¹H^α 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 ¹H^α 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.

A distal regulatory region of a class I human histone deacetylase

Histone deacetylases (HDACs) are key enzymes in epigenetics and important drug targets in cancer biology. Whilst it has been established that HDACs regulate many cellular processes, far less is known about the regulation of these enzymes themselves. Here, we show that HDAC8 is allosterically regulated by shifts in populations between exchanging states. An inactive state is identified, which is stabilised by a range of mutations and resembles a sparsely-populated state in equilibrium with active HDAC8. Computational models show that the inactive and active states differ by small changes in a regulatory region that extends up to 28 Å from the active site. The regulatory allosteric region identified here in HDAC8 corresponds to regions in other class I HDACs known to bind regulators, thus suggesting a general mechanism. The presented results pave the way for the development of allosteric HDAC inhibitors and regulators to improve the therapy for several disease states.

β-Lactamase of Mycobacterium tuberculosis Shows Dynamics in the Active Site That Increase upon Inhibitor Binding

The Mycobacterium tuberculosis β-lactamase BlaC is a broad-spectrum β-lactamase that can convert a range of β-lactam antibiotics. Enzymes with low specificity are expected to exhibit active-site flexibility. To probe the motions in BlaC, we studied the dynamic behavior in solution using nuclear magnetic resonance (NMR) spectroscopy. ¹⁵N relaxation experiments show that BlaC is mostly rigid on the pico- to nanosecond timescale. Saturation transfer experiments indicate that also on the high-millisecond timescale BlaC is not dynamic. Using relaxation dispersion experiments, clear evidence was obtained for dynamics in the low-millisecond range, with an exchange rate of ca. 860 s^-1 The dynamic amide groups are localized in the active site. Upon formation of an adduct with the inhibitor avibactam, extensive line broadening occurs, indicating an increase in magnitude of the active-site dynamics. Furthermore, the rate of the motions increases significantly. Upon reaction with the inhibitor clavulanic acid, similar line broadening is accompanied by duplication of NMR signals, indicative of at least one additional, slower exchange process (exchange rate, k _ex, of <100 s^-1), while for this inhibitor also loss of pico- to nanosecond timescale rigidity is observed for some amides in the α domain. Possible sources of the observed dynamics, such as motions in the omega loop and rearrangements of active-site residues, are discussed. The increase in dynamics upon ligand binding argues against a model of inhibitor binding through conformational selection. Rather, the induced dynamics may serve to maximize the likelihood of sampling the optimal conformation for hydrolysis of the bound ligand.

Methyl TROSY spectroscopy: A versatile NMR approach to study challenging biological systems

A major goal in structural biology is to unravel how molecular machines function in detail. To that end, solution-state NMR spectroscopy is ideally suited as it is able to study biological assemblies in a near natural environment. Based on methyl TROSY methods, it is now possible to record high-quality data on complexes that are far over 100 kDa in molecular weight. In this review, we discuss the theoretical background of methyl TROSY spectroscopy, the information that can be extracted from methyl TROSY spectra and approaches that can be used to assign methyl resonances in large complexes. In addition, we touch upon insights that have been obtained for a number of challenging biological systems, including the 20S proteasome, the RNA exosome, molecular chaperones and G-protein-coupled receptors. We anticipate that methyl TROSY methods will be increasingly important in modern structural biology approaches, where information regarding static structures is complemented with insights into conformational changes and dynamic intermolecular interactions.

Dynamic design: manipulation of millisecond timescale motions on the energy landscape of cyclophilin A

Proteins need to interconvert between many conformations in order to function, many of which are formed transiently, and sparsely populated. Particularly when the lifetimes of these states approach the millisecond timescale, identifying the relevant structures and the mechanism by which they interconvert remains a tremendous challenge. Here we introduce a novel combination of accelerated MD (aMD) simulations and Markov state modelling (MSM) to explore these 'excited' conformational states. Applying this to the highly dynamic protein CypA, a protein involved in immune response and associated with HIV infection, we identify five principally populated conformational states and the atomistic mechanism by which they interconvert. A rational design strategy predicted that the mutant D66A should stabilise the minor conformations and substantially alter the dynamics, whereas the similar mutant H70A should leave the landscape broadly unchanged. These predictions are confirmed using CPMG and R1ρ solution state NMR measurements. By efficiently exploring functionally relevant, but sparsely populated conformations with millisecond lifetimes in silico, our aMD/MSM method has tremendous promise for the design of dynamic protein free energy landscapes for both protein engineering and drug discovery.

Characterizing Post-Translational Modifications and Their Effects on Protein Conformation Using NMR Spectroscopy

The diversity of the cellular proteome substantially exceeds the number of genes coded by the DNA of an organism because one or more residues in a majority of eukaryotic proteins are post-translationally modified (PTM) by the covalent conjugation of specific chemical groups. We now know that PTMs alter protein conformation and function in ways that are not entirely understood at the molecular level. NMR spectroscopy has been particularly successful as an analytical tool in elucidating the themes underlying the structural role of PTMs. In this Perspective, we focus on the NMR-based characterization of three abundant PTMs: phosphorylation, acetylation, and glycosylation. We detail NMR methods that have found success in detecting these modifications at a site-specific level. We also highlight NMR studies that have mapped the conformational changes ensuing from these PTMs as well as evaluated their relation to function. The NMR toolbox is expanding rapidly with experiments available to probe not only the average structure of biomolecules but also how this structure changes with time on time scales ranging from picoseconds to seconds. The atomic resolution insights into the biomolecular structure, dynamics, and mechanism accessible from NMR spectroscopy ensure that NMR will continue to be at the forefront of research in the structural biology of PTMs.

Revisiting 1HN CPMG relaxation dispersion experiments: a simple modification can eliminate large artifacts

Carr-Purcell-Meiboom-Gill relaxation dispersion experiments are commonly used to probe biomolecular dynamics on the millisecond timescale. The simplest experiment involves using backbone ¹⁵N spins as probes of motion and pulse sequences are now available for providing accurate dispersion profiles in this case. In contrast, ¹H-based experiments recorded on fully protonated samples are less common because of difficulties associated with homonuclear scalar couplings that can result in transfer of magnetization between coupled spins, leading to significant artifacts. Herein we examine a version of the ¹H^N CPMG experiment that has been used in our laboratory where a pair of CPMG pulse trains comprising non-selective, high power ¹H refocusing pulses sandwich an amide selective pulse that serves to refocus scalar-coupled evolution by the end of the train. The origin of the artifacts in our original scheme is explained and a new, significantly improved sequence is presented. The utility of the new experiment is demonstrated by obtaining flat ¹H^N dispersion profiles in a protonated protein system that is not expected to undergo millisecond timescale dynamics, and subsequently by measuring profiles on a cavity mutant of T4 lysozyme that exchanges between a pair of distinct states, establishing that high quality data can be generated even for fully protonated samples.

Characterising side chains in large proteins by protonless 13C-detected NMR spectroscopy

Side chains cover protein surfaces and are fundamental to processes as diverse as substrate recognition, protein folding and enzyme catalysis. However, characterisation of side-chain motions has so far been restricted to small proteins and methyl-bearing side chains. Here we present a class of methods, based on ¹³C-detected NMR spectroscopy, to more generally quantify motions and interactions of side chains in medium-to-large proteins. A single, uniformly isotopically labelled sample is sufficient to characterise the side chains of six different amino acid types. Side-chain conformational dynamics on the millisecond time-scale can be quantified by incorporating chemical exchange saturation transfer (CEST) into the presented methods, whilst long-range ¹³C-¹³C scalar couplings reporting on nanosecond to millisecond motions can be quantified in proteins as large as 80 kDa. The presented class of methods promises characterisation of side-chain behaviour at a level that has so far been reserved for the protein backbone.

Analysis of side-chain motions by NMR has so far been restricted to small proteins and methyl-bearing side chains. Here, the authors present NMR methods based on ¹³C direct detection of highly deuterated protein samples that yield sharp and well-resolved signals and allow the characterisation of side-chain conformational dynamics of six different amino acid types in medium-to-large proteins.

Chemical Exchange

The phenomenon of chemical or conformational exchange in NMR spectroscopy has enabled detailed characterization of time-dependent aspects of biomolecular function, including folding, molecular recognition, allostery, and catalysis, on timescales from microsecond to second. Importantly, NMR methods based on a variety of spin relaxation parameters have been developed that provide quantitative information on interconversion kinetics, thermodynamic properties, and structural features of molecular states populated to a fraction of a percent at equilibrium and otherwise unobservable by other NMR approaches. The ongoing development of more sophisticated experimental techniques and the necessity to apply these methods to larger and more complex molecular systems engenders a corresponding need for theoretical advances describing such techniques and facilitating data analysis in applications. This review surveys current aspects of the theory of chemical exchange, as utilized in ZZ-exchange; Hahn and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo; and R_1ρ, chemical exchange saturation transfer (CEST), and dark state saturation transfer (DEST) spin-locking experiments. The review emphasizes theoretical results for kinetic topologies with more than two interconverting states, both to obtain compact analytical forms suitable for data analysis and to establish conditions for distinguishability between alternative kinetic schemes.

Characterization of intrinsically disordered proteins and their dynamic complexes: From in vitro to cell-like environments

Over the last two decades, it has become increasingly clear that a large fraction of the human proteome is intrinsically disordered or contains disordered segments of significant length. These intrinsically disordered proteins (IDPs) play important regulatory roles throughout biology, underlining the importance of understanding their conformational behavior and interaction mechanisms at the molecular level. Here we review recent progress in the NMR characterization of the structure and dynamics of IDPs in various functional states and environments. We describe the complementarity of different NMR parameters for quantifying the conformational propensities of IDPs in their isolated and phosphorylated states, and we discuss the challenges associated with obtaining structural models of dynamic protein-protein complexes involving IDPs. In addition, we review recent progress in understanding the conformational behavior of IDPs in cell-like environments such as in the presence of crowding agents, in membrane-less organelles and in the complex environment of the human cell.

Probing Protein-Protein Interactions Using Asymmetric Labeling and Carbonyl-Carbon Selective Heteronuclear NMR Spectroscopy

Protein-protein interactions (PPIs) regulate a plethora of cellular processes and NMR spectroscopy has been a leading technique for characterizing them at the atomic resolution. Technically, however, PPIs characterization has been challenging due to multiple samples required to characterize the hot spots at the protein interface. In this paper, we review our recently developed methods that greatly simplify PPI studies, which minimize the number of samples required to fully characterize residues involved in the protein-protein binding interface. This original strategy combines asymmetric labeling of two binding partners and the carbonyl-carbon label selective (CCLS) pulse sequence element implemented into the heteronuclear single quantum correlation (¹H-15N HSQC) spectra. The CCLS scheme removes signals of the J-coupled 15N⁻13C resonances and records simultaneously two individual amide fingerprints for each binding partner. We show the application to the measurements of chemical shift correlations, residual dipolar couplings (RDCs), and paramagnetic relaxation enhancements (PRE). These experiments open an avenue for further modifications of existing experiments facilitating the NMR analysis of PPIs.

Characterization of the interaction interface and conformational dynamics of human TGIF1 homeodomain upon the binding of consensus DNA

The TG interacting factor-1 homeodomain (TGIF1-HD) binds with the consensus DNA motif 5'-TGTCA-3' in gene promoters through its three-amino acid loop extension (TALE) type homeodomain, and then recruits co-regulators to regulate gene expression. Although the solution NMR structure of human TGIF1-HD has been reported previously, little is known about its DNA binding mechanism. NMR titrations have been extensively used to study mechanisms of ligand binding to target proteins; however, an intermediate exchange occurred predominantly between TGIF1-HD in the free and bound states when titrated with the consensus DNA, which resulted in poor-quality NMR spectra and precluded further exploration of its interaction interface and conformational dynamics. Here, the helix α3 of TGIF1-HD was speculated as the specific DNA binding interface by hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments, and subsequently confirmed by chemical exchange saturation transfer (CEST) spectroscopy. In addition, simultaneous conformational changes in other regions, including α1 and α2, were induced by DNA binding, explaining the observation of chemical shift perturbations from extensive residues besides those located in α3. Further, low-populated DNA-bound TGIF1-HD undergoing a slow exchange at a rate of 130.2 ± 3.6 s^-1 was derived from the analysis of the CEST data, and two residues, R220 and R221, located in the middle of α3 were identified to be crucial for DNA binding. Our study provides structural and dynamic insights into the mechanisms of TGIF1-HD recognition of extensive promoter DNA.

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.

Removal of slow-pulsing artifacts in in-phase 15N relaxation dispersion experiments using broadband 1H decoupling

Understanding of the molecular mechanisms of protein function requires detailed insight into the conformational landscape accessible to the protein. Conformational changes can be crucial for biological processes, such as ligand binding, protein folding, and catalysis. NMR spectroscopy is exquisitely sensitive to such dynamic changes in protein conformations. In particular, Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments are a powerful tool to investigate protein dynamics on a millisecond time scale. CPMG experiments that probe the chemical shift modulation of 15N in-phase magnetization are particularly attractive, due to their high sensitivity. These experiments require high power 1H decoupling during the CPMG period to keep the 15N magnetization in-phase. Recently, an improved version of the in-phase 15N-CPMG experiment was introduced, offering greater ease of use by employing a single 1H decoupling power for all CPMG pulsing rates. In these experiments however, incomplete decoupling of off-resonance amide 1H spins introduces an artefactual dispersion of relaxation rates, the so-called slow-pulsing artifact. Here, we analyze the slow-pulsing artifact in detail and demonstrate that it can be suppressed through the use of composite pulse decoupling (CPD). We report the performances of various CPD schemes and show that CPD decoupling based on the 90x-240y-90x element results in high-quality dispersion curves free of artifacts, even for amides with high 1H offset.

Studying sparsely populated conformational states in RNA combining chemical synthesis and solution NMR spectroscopy

Using chemical synthesis and solution NMR spectroscopy, RNA structural ensembles including a major ground state and minor populated excited states can be studied at atomic resolution. In this work, atom-specific ¹³C labeled RNA building blocks - a 5-¹³C-uridine and a 2,8-¹³C₂-adenosine building block - are used to introduce isolated ¹³C-¹H-spin topologies into a target RNA to probe such structural ensembles via NMR spectroscopy. First, the 5-¹³C-uridine 2'-O-TBDMS-phosphoramidite building block was introduced into a 21 nucleotide (nt) tP5c stem construct of the tP5abc subdomain of the Tetrahymena group I ribozyme. Then, the 2,8-¹³C₂-adenosine 2'-O-TBDMS-phosphoramidite building block was incorporated into a 9 kDa and a 15 kD construct derived from the epsilon (ε) RNA element of the duck Hepatitis B virus. The 2,8-¹³C₂-adenosine resonances of the 9 kDa 28 nt sequence could be mapped to the full-length 53 nt construct. The isolated NMR active nuclei pairs were used to probe for low populated excited states (<10%) via ¹³C-Carr-Purcell-Meiboom-Gill (CPMG)-relaxation dispersion NMR spectroscopy. The ¹³C-CPMG relaxation dispersion experiment recapitulated a secondary structure switching event in the P5c hairpin of the group I intron construct previously revealed by ¹⁵N relaxation dispersion experiments. In the ε-HBV RNA an unfolding event occurring on the millisecond time scale was found in the upper stem in-line with earlier observations. This unpaired conformational state is presumed to be important for the binding of the epsilon reverse transcriptase (RT) enzyme. Thus, a full description of an RNA's folding landscape helps to obtain a deeper understanding of its function, as these high energy conformational states often represent functionally important intermediates involved in (un)folding or ribozyme catalysis.

Methyl NMR spectroscopy: Measurement of dynamics in viral RNA-directed RNA polymerases

Measurement of nuclear spin relaxation provides a powerful approach to access information about biomolecular conformational dynamics over several orders of magnitude in timescale. In several cases this knowledge in combination with spatial information from three-dimensional structures yields unique insight into protein stability and the kinetics and thermodynamics of their interactions and function. However, due to intrinsic difficulties in studying large systems using solution state nuclear magnetic resonance (NMR) approaches, until recently these measurements were limited to small-to-medium-sized systems. However, the development of a wide range of novel strategies that allow the selective isotope labeling of methyl groups in proteins have allowed the exploitation of the unique relaxation properties of this spin-system. This has in turn enabled the extension of NMR approaches to high molecular weight proteins including a variety of enzymes and their complexes. Here, we recount our experiences in obtaining assignments of the methyl resonances for two representative members of a class of RNA-directed RNA polymerases (RdRps) encoded by bacteriophages of the Cystoviridae family. We demonstrate the utility of these methyl probes, limited in number for one case and more numerous for the other, to investigate the conformational dynamics of RdRps on the fast (ps-ns) and slow (μs-ms) timescales.

Tyr51: Key Determinant of the Low Thermostability of the Colwellia psychrerythraea Cold-Shock Protein

Cold-shock proteins (Csps) are expressed at lower-than-optimum temperatures, and they function as RNA chaperones; however, no structural studies on psychrophilic Csps have been reported. Here, we aimed to investigate the structure and dynamics of the Csp of psychrophile Colwellia psychrerythraea 34H, ( Cp-Csp). Although Cp-Csp shares sequence homology, common folding patterns, and motifs, including a five β-stranded barrel, with its thermophilic counterparts, its thermostability (37 °C) was markedly lower than those of other Csps. Cp-Csp binds heptathymidine with an affinity of 10^-7 M, thereby increasing its thermostability to 50 °C. Nuclear magnetic resonance spectroscopic analysis of the Cp-Csp structure and backbone dynamics revealed a flexible structure with only one salt bridge and 10 residues in the hydrophobic cavity. Notably, Cp-Csp contains Tyr51 instead of the conserved Phe in the hydrophobic core, and its phenolic hydroxyl group projects toward the surface. The Y51F mutation increased the stability of hydrophobic packing and may have allowed for the formation of a K3-E21 salt bridge, thereby increasing its thermostability to 43 °C. Cp-Csp exhibited conformational exchanges in its ribonucleoprotein motifs 1 and 2 (754 and 642 s^-1), and heptathymidine binding markedly decreased these motions. Cp-Csp lacks salt bridges and has longer flexible loops and a less compact hydrophobic cavity resulting from Tyr51 compared to mesophilic and thermophilic Csps. These might explain the low thermostability of Cp-Csp. The conformational flexibility of Cp-Csp facilitates its accommodation of nucleic acids at low temperatures in polar oceans and its function as an RNA chaperone for cold adaptation.

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

Targeting a Dark Excited State of HIV-1 Nucleocapsid by Antiretroviral Thioesters Revealed by NMR Spectroscopy

HIV-1 nucleocapsid (NCp7) is a two Cys₂ HisCys zinc knuckle (N-Zn and C-Zn) protein that plays a key role in viral replication. NCp7 conformational dynamics is characterized by NMR relaxation dispersion and chemical exchange saturation transfer measurements. While the N-Zn knuckle is conformationally stable, the C-Zn knuckle interconverts on the millisecond timescale between the major state, in which the zinc is coordinated by three cysteines and a histidine, and two folded minor species (with populations around 1 %) in which one of the coordination bonds (Cys413-Sγ-Zn or His421-Nϵ2-Zn) is hydrolyzed. These findings explain why antiretroviral thioesters specifically disrupt the C-Zn knuckle by initial acylation of Cys413, and show that transient, sparsely-populated ("dark"), excited states of proteins can present effective targets for rational drug design.

Effect of amino acid mutations on the conformational dynamics of amyloidogenic immunoglobulin light-chains: A combined NMR and in silico study

The conformational dynamics of a pathogenic κ4 human immunoglobulin light-chain variable domain, SMA, associated with AL amyloidosis, were investigated by ¹⁵N relaxation dispersion NMR spectroscopy. Compared to a homologous light-chain, LEN, which differs from SMA at eight positions but is non-amyloidogenic in vivo, we find that multiple residues in SMA clustered around the N-terminus and CDR loops experience considerable conformational exchange broadening caused by millisecond timescale protein motions, consistent with a destabilized dimer interface. To evaluate the contribution of each amino acid substitution to shaping the dynamic conformational landscape of SMA, NMR studies were performed for each SMA-like point mutant of LEN followed by in silico analysis for a subset of these proteins. These studies show that a combination of only three mutations located within or directly adjacent to CDR3 loop at the dimer interface, which remarkably include both destabilizing (Q89H and Y96Q) and stabilizing (T94H) mutations, largely accounts for the differences in conformational flexibility between LEN and SMA. Collectively, our studies indicate that a correct combination of stabilizing and destabilizing mutations is key for immunoglobulin light-chains populating unfolded intermediates that result in amyloid formation, and underscore the complex nature of correlations between light-chain conformational flexibility, thermodynamic stability and amyloidogenicity.

Measurement of protein backbone 13CO and 15N relaxation dispersion at high resolution

Peak overlap in crowded regions of two-dimensional spectra prevents characterization of dynamics for many sites of interest in globular and intrinsically disordered proteins. We present new three-dimensional pulse sequences for measurement of Carr-Purcell-Meiboom-Gill relaxation dispersions at backbone nitrogen and carbonyl positions. To alleviate increase in the measurement time associated with the additional spectral dimension, we use non-uniform sampling in combination with two distinct methods of spectrum reconstruction: compressed sensing and co-processing with multi-dimensional decomposition. The new methodology was validated using disordered protein CD79A from B-cell receptor and an SH3 domain from Abp1p in exchange between its free form and bound to a peptide from the protein Ark1p. We show that, while providing much better resolution, the 3D NUS experiments give the similar accuracy and precision of the dynamic parameters to ones obtained using traditional 2D experiments. Furthermore, we show that jackknife resampling of the spectra yields robust estimates of peak intensities errors, eliminating the need for recording duplicate data points.

Structure and Dynamics of an Intrinsically Disordered Protein Region That Partially Folds upon Binding by Chemical-Exchange NMR

Many intrinsically disordered proteins (IDPs) and protein regions (IDRs) engage in transient, yet specific, interactions with a variety of protein partners. Often, if not always, interactions with a protein partner lead to partial folding of the IDR. Characterizing the conformational space of such complexes is challenging: in solution-state NMR, signals of the IDR in the interacting region become broad, weak, and often invisible, while X-ray crystallography only provides information on fully ordered regions. There is thus a need for a simple method to characterize both fully and partially ordered regions in the bound state of IDPs. Here, we introduce an approach based on monitoring chemical exchange by NMR to investigate the state of an IDR that folds upon binding through the observation of the free state of the protein. Structural constraints for the bound state are obtained from chemical shifts, and site-specific dynamics of the bound state are characterized by relaxation rates. The conformation of the interacting part of the IDR was determined and subsequently docked onto the structure of the folded partner. We apply the method to investigate the interaction between the disordered C-terminal region of Artemis and the DNA binding domain of Ligase IV. We show that we can accurately reproduce the structure of the core of the complex determined by X-ray crystallography and identify a broader interface. The method is widely applicable to the biophysical investigation of complexes of disordered proteins and folded proteins.

Biophysical characterization of the calmodulin-like domain of Plasmodium falciparum calcium dependent protein kinase 3

Calcium dependent protein kinases are unique to plants and certain parasites and comprise an N-terminal segment and a kinase domain that is regulated by a C-terminal calcium binding domain. Since the proteins are not found in man they are potential drug targets. We have characterized the calcium binding lobes of the regulatory domain of calcium dependent protein kinase 3 from the malaria parasite Plasmodium falciparum. Despite being structurally similar, the two lobes differ in several other regards. While the monomeric N-terminal lobe changes its structure in response to calcium binding and shows global dynamics on the sub-millisecond time-scale both in its apo and calcium bound states, the C-terminal lobe could not be prepared calcium-free and forms dimers in solution. If our results can be generalized to the full-length protein, they suggest that the C-terminal lobe is calcium bound even at basal levels and that activation is caused by the structural reorganization associated with binding of a single calcium ion to the N-terminal lobe.

Kinetic Insights into the Binding between the nSH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The interaction between CrkII and cAbl is implicated in diverse cellular processes. This interaction starts with the binding of the N-terminal Src homology 3 (nSH3) domain of CrkII to the proline-rich motifs of cAbl (PRMs^cAbl). Despite its critical importance, the detailed binding mechanism between the nSH3 domain and PRMs remains elusive. In this study, we used nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiment to study the binding kinetics between the nSH3 domain of CrkII and PRMs^cAbl. Our results highlight that the nSH3 domain binds to three PRMs^cAbl with very high on- and off-rate constants, indicating the transient nature of the binding. To further characterize the binding transition state, we conducted the Eyring and linear free energy relationship analyses using temperature-dependent kinetic data. These data indicate that the binding transition state of the nSH3 domain and PRM is accompanied by small activation enthalpy, owing to partial desolvation of the transition state. These results also highlight the similarity between the transition and free states, in terms of structure and energetics. Although the binding of the nSH3 domain and PRM displays the features consistent with a diffusion-limited process within our experimental conditions, further tests are necessary to determine if the binding is a true diffusion-limited process.

When fast is better: protein folding fundamentals and mechanisms from ultrafast approaches

Protein folding research stalled for decades because conventional experiments indicated that proteins fold slowly and in single strokes, whereas theory predicted a complex interplay between dynamics and energetics resulting in myriad microscopic pathways. Ultrafast kinetic methods turned the field upside down by providing the means to probe fundamental aspects of folding, test theoretical predictions and benchmark simulations. Accordingly, experimentalists could measure the timescales for all relevant folding motions, determine the folding speed limit and confirm that folding barriers are entropic bottlenecks. Moreover, a catalogue of proteins that fold extremely fast (microseconds) could be identified. Such fast-folding proteins cross shallow free energy barriers or fold downhill, and thus unfold with minimal co-operativity (gradually). A new generation of thermodynamic methods has exploited this property to map folding landscapes, interaction networks and mechanisms at nearly atomic resolution. In parallel, modern molecular dynamics simulations have finally reached the timescales required to watch fast-folding proteins fold and unfold in silico. All of these findings have buttressed the fundamentals of protein folding predicted by theory, and are now offering the first glimpses at the underlying mechanisms. Fast folding appears to also have functional implications as recent results connect downhill folding with intrinsically disordered proteins, their complex binding modes and ability to moonlight. These connections suggest that the coupling between downhill (un)folding and binding enables such protein domains to operate analogically as conformational rheostats.

Exploration of the Conformational Dynamics of Major Histocompatibility Complex Molecules

Major histocompatibility complex (MHC) molecules are loaded with a wide variety of self- and non-self-peptides in their binding grooves and present these to T cell receptors (TCRs) in order to activate the adaptive immune system. A large number of crystal structures of different MHC alleles with different bound peptides have been determined, and they have been found to be quite similar to one another regardless of the bound peptide sequence. The structures do not change markedly even when forming complexes with TCRs. Nonetheless, the degree of TCR activation does differ markedly depending on the peptide presented by the MHC. Recent structural studies in solution rather than as crystals have suggested that the conformational dynamics of MHC molecules may be responsible for the MHC stability differences. Furthermore, it was shown that the conformational dynamics of MHC molecules is important for peptide loading and presentation to TCR. Here, we describe the static and dynamic structures of MHC molecules and appropriate methods to analyze them. We focus particularly on nuclear magnetic resonance (NMR), one of the most powerful tools to study dynamic properties of proteins. The number of such studies in the literature is limited, but in this review, we show that NMR is valuable for elucidating the structural dynamics of MHC molecules.

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.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

(13)C-NMR studies on disulfide bond isomerization in bovine pancreatic trypsin inhibitor (BPTI)

Conformational isomerization of disulfide bonds is associated with the dynamics and thus the functional aspects of proteins. However, our understanding of the isomerization is limited by experimental difficulties in probing it. We explored the disulfide conformational isomerization of the Cys14-Cys38 disulfide bond in bovine pancreatic trypsin inhibitor (BPTI), by performing an NMR line-shape analysis of its Cys carbon peaks. In this approach, 1D (13)C spectra were recorded at small temperature intervals for BPTI samples selectively labeled with site-specifically (13)C-enriched Cys, and the recorded peaks were displayed in the order of the temperature after the spectral scales were normalized to a carbon peak. Over the profile of the line-shape, exchange broadening that altered with temperature was manifested for the carbon peaks of Cys14 and Cys38. The Cys14-Cys38 disulfide bond reportedly exists in equilibrium between a high-populated (M) and two low-populated states (m c14 and m c38). Consistent with the three-site exchange model, biphasic exchange broadening arising from the two processes was observed for the peak of the Cys14 α-carbon. As the exchange broadening is maximized when the exchange rate equals the chemical shift difference in Hz between equilibrating sites, semi-quantitative information that was useful for establishing conditions for (13)C relaxation dispersion experiments was obtained through the carbon line-shape profile. With respect to the m c38 isomerization, the (1)H-(13)C signals at the β-position of the minor state were resolved from the major peaks and detected by exchange experiments at a low temperature.

Evaluating the influence of initial magnetization conditions on extracted exchange parameters in NMR relaxation experiments: applications to CPMG and CEST

Transient excursions of native protein states to functionally relevant higher energy conformations often occur on the μs-ms timescale. NMR spectroscopy has emerged as an important tool to probe such processes using techniques such as Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion and Chemical Exchange Saturation Transfer (CEST). The extraction of kinetic and structural parameters from these measurements is predicated upon mathematical modeling of the resulting relaxation profiles, which in turn relies on knowledge of the initial magnetization conditions at the start of the CPMG/CEST relaxation elements in these experiments. Most fitting programs simply assume initial magnetization conditions that are given by equilibrium populations, which may be incorrect in certain implementations of experiments. In this study we have quantified the systematic errors in extracted parameters that are generated from analyses of CPMG and CEST experiments using incorrect initial boundary conditions. We find that the errors in exchange rates (k ex ) and populations (p E ) are typically small (<10 %) and thus can be safely ignored in most cases. However, errors become larger and cannot be fully neglected (20-40 %) as k ex falls near the lower limit of each method or when short CPMG/CEST relaxation elements are used in these experiments. The source of the errors can be rationalized and their magnitude given by a simple functional form. Despite the fact that errors tend to be small, it is recommended that the correct boundary conditions be implemented in fitting programs so as to obtain as robust estimates of exchange parameters as possible.