An exchange-free measure of 15N transverse relaxation: an NMR spectroscopy application to the study of a folding intermediate with pervasive chemical exchange

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.

Poly(ADP-ribosyl)ation enhances nucleosome dynamics and organizes DNA damage repair components within biomolecular condensates

Nucleosomes, the basic structural units of chromatin, hinder recruitment and activity of various DNA repair proteins, necessitating modifications that enhance DNA accessibility. Poly(ADP-ribosyl)ation (PARylation) of proteins near damage sites is an essential initiation step in several DNA-repair pathways; however, its effects on nucleosome structural dynamics and organization are unclear. Using NMR, cryoelectron microscopy (cryo-EM), and biochemical assays, we show that PARylation enhances motions of the histone H3 tail and DNA, leaving the configuration of the core intact while also stimulating nuclease digestion and ligation of nicked nucleosomal DNA by LIG3. PARylation disrupted interactions between nucleosomes, preventing self-association. Addition of LIG3 and XRCC1 to PARylated nucleosomes generated condensates that selectively partition DNA repair-associated proteins in a PAR- and phosphorylation-dependent manner in vitro. Our results establish that PARylation influences nucleosomes across different length scales, extending from the atom-level motions of histone tails to the mesoscale formation of condensates with selective compositions.

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.

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.

Correlating histone acetylation with nucleosome core particle dynamics and function

Epigenetic modifications of chromatin play a critical role in regulating the fidelity of the genetic code and in controlling the translation of genetic information into the protein components of the cell. One key posttranslational modification is acetylation of histone lysine residues. Molecular dynamics simulations, and to a smaller extent experiment, have established that lysine acetylation increases the dynamics of histone tails. However, a systematic, atomic resolution experimental investigation of how this epigenetic mark, focusing on one histone at a time, influences the structural dynamics of the nucleosome beyond the tails, and how this translates into accessibility of protein factors such as ligases and nucleases, has yet to be performed. Herein, using NMR spectroscopy of nucleosome core particles (NCPs), we evaluate the effects of acetylation of each histone on tail and core dynamics. We show that for histones H2B, H3, and H4, the histone core particle dynamics are little changed, even though the tails have increased amplitude motions. In contrast, significant increases to H2A dynamics are observed upon acetylation of this histone, with the docking domain and L1 loop particularly affected, correlating with increased susceptibility of NCPs to nuclease digestion and more robust ligation of nicked DNA. Dynamic light scattering experiments establish that acetylation decreases inter-NCP interactions in a histone-dependent manner and facilitates the development of a thermodynamic model for NCP stacking. Our data show that different acetylation patterns result in nuanced changes to NCP dynamics, modulating interactions with other protein factors, and ultimately controlling biological output.

Surface electrostatics dictate RNA-binding protein CAPRIN1 condensate concentration and hydrodynamic properties

Biomolecular condensates concentrate proteins, nucleic acids, and small molecules and play an essential role in many biological processes. Their formation is tuned by a balance between energetically favorable and unfavorable contacts, with charge-charge interactions playing a central role in some systems. The positively charged intrinsically disordered carboxy-terminal region of the RNA-binding protein CAPRIN1 is one such example, phase separating upon addition of negatively charged ATP or high concentrations of sodium chloride (NaCl). Using solution NMR spectroscopy, we measured residue-specific near-surface electrostatic potentials (ϕENS) of CAPRIN1 along its NaCl-induced phase separation trajectory to compare with those obtained using ATP. In both cases, electrostatic shielding decreases ϕENS values, yet surface potentials of CAPRIN1 in the two condensates can be different, depending on the amount of NaCl or ATP added. Our results establish that even small differences in ϕENS can significantly affect the level of protein enrichment and the mechanical properties of the condensed phase, leading, potentially, to the regulation of biological processes.

NMR Methods to Study the Dynamics of SH2 Domain-Phosphopeptide Complexes

Nuclear magnetic resonance (NMR) spectroscopy is the method of choice for studying the dynamics of biological macromolecules in solution. By exploiting the intricate interplay between the effects of protein motion (both overall rotational diffusion and internal mobility) and nuclear spin relaxation, NMR allows molecular motion to be probed at atomic resolution over a wide range of timescales, including picosecond (bond vibrations and methyl-group rotations), nanosecond (loop motions and rotational diffusion), and microsecond-millisecond (ligand binding, allostery). In this chapter, we describe different NMR pulse schemes (R1, R1ρ, heteronuclear NOE, and CPMG relaxation dispersion) to characterize the dynamics of SH2 domains. As an example, we use the N-SH2 domain of protein tyrosine phosphatase SHP2 in complex with two phosphopeptides derived from immune checkpoint receptor PD-1 (ITIM and ITSM).

Mechanism of Crosstalk between the LSD1 Demethylase and HDAC1 Deacetylase in the CoREST Complex

The transcriptional corepressor complex CoREST is one of seven histone deacetylase complexes that regulate the genome through controlling chromatin acetylation. The CoREST complex is unique in containing both histone demethylase and deacetylase enzymes, LSD1 and HDAC1, held together by the RCOR1 scaffold protein. To date, it has been assumed that the enzymes function independently within the complex. Now, we report the assembly of the ternary complex. Using both structural and functional studies, we show that the activity of the two enzymes is closely coupled and that the complex can exist in at least two distinct states with different kinetics. Electron microscopy of the complex reveals a bi-lobed structure with LSD1 and HDAC1 enzymes at opposite ends of the complex. The structure of CoREST in complex with a nucleosome reveals a mode of chromatin engagement that contrasts with previous models.

•

The activities of LSD1 and HDAC1 are closely coupled in the CoREST complex

•

Both LSD1 and HDAC1 exist in two different kinetic states

•

CoREST has a bi-lobed, flexible structure with the two enzymes located at opposite ends

•

CoREST interacts with methylated nucleosomes via LSD1, but not HDAC1 or RCOR1

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.

Multiquantum Chemical Exchange Saturation Transfer NMR to Quantify Symmetrical Exchange: Application to Rotational Dynamics of the Guanidinium Group in Arginine Side Chains

Chemical exchange saturation transfer (CEST) NMR experiments have emerged as a powerful tool for characterizing dynamics in proteins. We show here that the CEST approach can be extended to systems with symmetrical exchange, where the NMR signals of all exchanging species are severely broadened. To achieve this, multiquantum CEST (MQ-CEST) is introduced, where the CEST pulse is applied to a longitudinal multispin order density element and the CEST profiles are encoded onto nonbroadened nuclei. The MQ-CEST approach is demonstrated on the restricted rotation of guanidinium groups in arginine residues within proteins. These groups and their dynamics are essential for many enzymes and for noncovalent interactions through the formation of hydrogen bonds, salt-bridges, and π-stacking interactions, and their rate of rotation is highly indicative of the extent of interactions formed. The MQ-CEST method is successfully applied to guanidinium groups in the 19 kDa L99A mutant of T4 lysozyme.

Activation Loop Dynamics Are Coupled to Core Motions in Extracellular Signal-Regulated Kinase-2

The activation loop segment in protein kinases is a common site for regulatory phosphorylation. In extracellular signal-regulated kinase 2 (ERK2), dual phosphorylation and conformational rearrangement of the activation loop accompany enzyme activation. X-ray structures show the active conformation to be stabilized by multiple ion pair interactions between phosphorylated threonine and tyrosine residues in the loop and six arginine residues in the kinase core. Despite the extensive salt bridge network, nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiments show that the phosphorylated activation loop is conformationally mobile on a microsecond to millisecond time scale. The dynamics of the loop match those of previously reported global exchange within the kinase core region and surrounding the catalytic site that have been found to facilitate productive nucleotide binding. Mutations in the core region that alter these global motions also alter the dynamics of the activation loop. Conversely, mutations in the activation loop perturb the global exchange within the kinase core. Together, these findings provide evidence for coupling between motions in the activation loop and those surrounding the catalytic site in the active state of the kinase. Thus, the activation loop segment in dual-phosphorylated ERK2 is not held statically in the active X-ray conformation but instead undergoes exchange between conformers separated by a small energetic barrier, serving as part of a dynamic allosteric network controlling nucleotide binding and catalytic function.

β-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.

Determining Binding Kinetics of Intrinsically Disordered Proteins by NMR Spectroscopy

The unique structural flexibility of intrinsically disordered proteins (IDPs) is central to their diverse functions in cellular processes. Protein-protein interactions involving IDPs are frequently transient and dynamic in nature. Nuclear magnetic resonance (NMR) spectroscopy is an especially powerful tool for characterizing the structural propensities, dynamics, and interactions of IDPs. Here we describe applications of the Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiment in combination with NMR titrations to characterize the kinetics and mechanisms of interactions between intrinsically disordered proteins and their targets. We illustrate the method with reference to interactions between the activation domain of the human T-cell leukemia virus type-I (HTLV-1) basic leucine zipper protein (HBZ) and its cellular binding partner, the KIX domain of the transcriptional coactivator CBP.

Gcn4-Mediator Specificity Is Mediated by a Large and Dynamic Fuzzy Protein-Protein Complex

Transcription activation domains (ADs) are inherently disordered proteins that often target multiple coactivator complexes, but the specificity of these interactions is not understood. Efficient transcription activation by yeast Gcn4 requires its tandem ADs and four activator-binding domains (ABDs) on its target, the Mediator subunit Med15. Multiple ABDs are a common feature of coactivator complexes. We find that the large Gcn4-Med15 complex is heterogeneous and contains nearly all possible AD-ABD interactions. Gcn4-Med15 forms via a dynamic fuzzy protein-protein interface, where ADs bind the ABDs in multiple orientations via hydrophobic regions that gain helicity. This combinatorial mechanism allows individual low-affinity and specificity interactions to generate a biologically functional, specific, and higher affinity complex despite lacking a defined protein-protein interface. This binding strategy is likely representative of many activators that target multiple coactivators, as it allows great flexibility in combinations of activators that can cooperate to regulate genes with variable coactivator requirements.

Interaction Between the a3 Region of Factor VIII and the TIL'E' Domains of the von Willebrand Factor

The von Willebrand factor (VWF) and coagulation factor VIII (FVIII) are intricately involved in hemostasis. A tight, noncovalent complex between VWF and FVIII prolongs the half-life of FVIII in plasma, and failure to form this complex leads to rapid clearance of FVIII and bleeding diatheses such as hemophilia A and von Willebrand disease (VWD) type 2N. High-resolution insight into the complex between VWF and FVIII has so far been strikingly lacking. This is particularly the case for the flexible a3 region of FVIII, which is imperative for high-affinity binding. Here, a structural and biophysical characterization of the interaction between VWF and FVIII is presented with focus on two of the domains that have been proven pivotal for mediating the interaction, namely the a3 region of FVIII and the TIL'E' domains of VWF. Binding between the FVIII a3 region and VWF TIL'E' was here observed using NMR spectroscopy, where chemical shift changes were localized to two β-sheet regions on the edge of TIL'E' upon FVIII a3 region binding. Isothermal titration calorimetry and NMR spectroscopy were used to characterize the interaction between FVIII and TIL'E' as well as mutants of TIL'E', which further highlights the importance of the β-sheet region of TIL'E' for high-affinity binding. Overall, the results presented provide new insight into the role the FVIII a3 region plays for complex formation between VWF and FVIII and the β-sheet region of TIL'E' is shown to be important for FVIII binding. Thus, the results pave the way for further high-resolution insights into this imperative complex.

Frequency selective coherence transfer NMR spectroscopy to study the structural dynamics of high molecular weight proteins

Multidimensional nuclear magnetic resonance (NMR) spectroscopy has enabled detailed characterizations of protein structures and dynamics that are closely linked to functions. However, it leads to a large sensitivity loss in applications to high molecular weight proteins, which is caused by spin relaxation during the frequency discrimination period in the indirect dimension. Here, we describe a selective coherence transfer scheme, which enables us to selectively observe ¹H nuclei bonded to ¹⁵N or ¹³C nuclei with specified resonance frequencies. By utilizing this scheme, we achieved a 2.5- to 6-fold increase in signal height per unit of time with this scheme by avoiding the relaxation loss in the indirect dimension, as compared to the conventional two-dimensional heteronuclear correlation spectroscopy. We also demonstrated the effectiveness of this approach with applications to the membrane protein KirBac1.1, and characterized the functionally relevant conformational exchange process in both detergent micelles and a reconstituted membrane environment, corresponding to the apparent molecular masses of 220 kDa and 300 kDa, respectively.

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.

Arginine Side-Chain Hydrogen Exchange: Quantifying Arginine Side-Chain Interactions in Solution

The rate with which labile backbone hydrogen atoms in proteins exchange with the solvent has long been used to probe protein interactions in aqueous solutions. Arginine, an essential amino acid found in many interaction interfaces, is capable of an impressive range of interactions via its guanidinium group. The hydrogen exchange rate of the guanidinium hydrogens therefore becomes an important measure to quantify side-chain interactions. Herein we present an NMR method to quantify the hydrogen exchange rates of arginine side-chain 1 Hϵ protons and thus present a method to gauge the strength of arginine side-chain interactions. The method employs 13 C-detection and the one-bond deuterium isotope shift observed for 15 Nϵ to generate two exchanging species in 1 H2 O/2 H2 O mixtures. An application to the protein T4 Lysozyme is shown, where protection factors calculated from the obtained exchange rates correlate well with the interactions observed in the crystal structure. The methodology presented provides an important step towards characterising interactions of arginine side-chains in enzymes, in phase separation, and in protein interaction interfaces in general.

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.

Ab Initio Prediction of NMR Spin Relaxation Parameters from Molecular Dynamics Simulations

¹H-¹⁵N NMR spin relaxation and relaxation dispersion experiments can reveal the time scale and extent of protein motions across the ps-ms range, where the ps-ns dynamics revealed by fundamental quantities R₁, R₂, and heteronuclear NOE can be well-sampled by molecular dynamics simulations (MD). Although the principles of relaxation prediction from simulations are well-established, numerous NMR-MD comparisons have hitherto focused upon the aspect of order parameters S² due to common artifacts in the prediction of transient dynamics. We therefore summarize here all necessary components and highlight existing and proposed solutions, such as the inclusion of quantum mechanical zero-point vibrational corrections and separate MD convergence of global and local motions in coarse-grained and all-atom force fields, respectively. For the accuracy of the MD prediction to be tested, two model proteins GB3 and Ubiquitin are used to validate five atomistic force fields against published NMR data supplemented by the coarse-grained force field MARTINI+EN. In Amber and CHARMM-type force fields, quantitative agreement was achieved for structured elements with minimum adjustment of global parameters. Deviations from experiment occur in flexible loops and termini, indicating differences in both the extent and time scale of backbone motions. The lack of systematic patterns and water model dependence suggests that modeling of the local environment limits prediction accuracy. Nevertheless, qualitative accuracy in a 2 μs CHARMM36m Stam2 VHS domain simulation demonstrates the potential of MD-based interpretation in combination with NMR-measured dynamics, increasing the utility of spin relaxation in integrative structural biology.

Comprehensive analysis of NMR data using advanced line shape fitting

NMR spectroscopy is uniquely suited for atomic resolution studies of biomolecules such as proteins, nucleic acids and metabolites, since detailed information on structure and dynamics are encoded in positions and line shapes of peaks in NMR spectra. Unfortunately, accurate determination of these parameters is often complicated and time consuming, in part due to the need for different software at the various analysis steps and for validating the results. Here, we present an integrated, cross-platform and open-source software that is significantly more versatile than the typical line shape fitting application. The software is a completely redesigned version of PINT ( https://pint-nmr.github.io/PINT/ ). It features a graphical user interface and includes functionality for peak picking, editing of peak lists and line shape fitting. In addition, the obtained peak intensities can be used directly to extract, for instance, relaxation rates, heteronuclear NOE values and exchange parameters. In contrast to most available software the entire process from spectral visualization to preparation of publication-ready figures is done solely using PINT and often within minutes, thereby, increasing productivity for users of all experience levels. Unique to the software are also the outstanding tools for evaluating the quality of the fitting results and extensive, but easy-to-use, customization of the fitting protocol and graphical output. In this communication, we describe the features of the new version of PINT and benchmark its performance.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Measurement of 15N longitudinal relaxation rates in 15NH4+ spin systems to characterise rotational correlation times and chemical exchange

Many chemical and biological processes rely on the movement of monovalent cations and an understanding of such processes can therefore only be achieved by characterising the dynamics of the involved ions. It has recently been shown that 15N-ammonium can be used as a proxy for potassium to probe potassium binding in bio-molecules such as DNA quadruplexes and enzymes. Moreover, equations have been derived to describe the time-evolution of 15N-based spin density operator elements of 15NH4+ spin systems. Herein NMR pulse sequences are derived to select specific spin density matrix elements of the 15NH4+ spin system and to measure their longitudinal relaxation in order to characterise the rotational correlation time of the 15NH4+ ion as well as report on chemical exchange events of the 15NH4+ ion. Applications to 15NH4+ in acidic aqueous solutions are used to cross-validate the developed pulse sequence while measurements of spin-relaxation rates of 15NH4+ bound to a 41kDa domain of the bacterial Hsp70 homologue DnaK are presented to show the general applicability of the derived pulse sequence. The rotational correlation time obtained for 15N-ammonium bound to DnaK is similar to the correlation time that describes the rotation about the threefold axis of a methyl group. The methodology presented here provides, together with the previous theoretical framework, an important step towards characterising the motional properties of cations in macromolecular systems.

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.

Characterizing Active Site Conformational Heterogeneity along the Trajectory of an Enzymatic Phosphoryl Transfer Reaction

States along the phosphoryl transfer reaction catalyzed by the nucleoside monophosphate kinase UmpK were captured and changes in the conformational heterogeneity of conserved active site arginine side‐chains were quantified by NMR spin‐relaxation methods. In addition to apo and ligand‐bound UmpK, a transition state analog (TSA) complex was utilized to evaluate the extent to which active site conformational entropy contributes to the transition state free energy. The catalytically essential arginine side‐chain guanidino groups were found to be remarkably rigid in the TSA complex, indicating that the enzyme has evolved to restrict the conformational freedom along its reaction path over the energy landscape, which in turn allows the phosphoryl transfer to occur selectively by avoiding side reactions.

Spectral density mapping at multiple magnetic fields suitable for (13)C NMR relaxation studies

Standard spectral density mapping protocols, well suited for the analysis of (15)N relaxation rates, introduce significant systematic errors when applied to (13)C relaxation data, especially if the dynamics is dominated by motions with short correlation times (small molecules, dynamic residues of macromolecules). A possibility to improve the accuracy by employing cross-correlated relaxation rates and on measurements taken at several magnetic fields has been examined. A suite of protocols for analyzing such data has been developed and their performance tested. Applicability of the proposed protocols is documented in two case studies, spectral density mapping of a uniformly labeled RNA hairpin and of a selectively labeled disaccharide exhibiting highly anisotropic tumbling. Combination of auto- and cross-correlated relaxation data acquired at three magnetic fields was applied in the former case in order to separate effects of fast motions and conformational or chemical exchange. An approach using auto-correlated relaxation rates acquired at five magnetic fields, applicable to anisotropically moving molecules, was used in the latter case. The results were compared with a more advanced analysis of data obtained by interpolation of auto-correlated relaxation rates measured at seven magnetic fields, and with the spectral density mapping of cross-correlated relaxation rates. The results showed that sufficiently accurate values of auto- and cross-correlated spectral density functions at zero and (13)C frequencies can be obtained from data acquired at three magnetic fields for uniformly (13)C-labeled molecules with a moderate anisotropy of the rotational diffusion tensor. Analysis of auto-correlated relaxation rates at five magnetic fields represents an alternative for molecules undergoing highly anisotropic motions.

Relaxation dispersion NMR spectroscopy for the study of protein allostery

Allosteric transmission of information between distant sites in biological macromolecules often involves collective transitions between active and inactive conformations. Nuclear magnetic resonance (NMR) spectroscopy can yield detailed information on these dynamics. In particular, relaxation dispersion techniques provide structural, dynamic, and mechanistic information on conformational transitions occurring on the millisecond to microsecond timescales. In this review, we provide an overview of the theory and analysis of Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments and briefly describe their application to the study of allosteric dynamics in the homeodomain from the PBX transcription factor (PBX-HD). CPMG NMR data show that local folding (helix/coil) transitions in one part of PBX-HD help to communicate information between two distant binding sites. Furthermore, the combination of CPMG and other spin relaxation data show that this region can also undergo local misfolding, reminiscent of conformational ensemble models of allostery.

Heteronuclear transverse and longitudinal relaxation in AX4 spin systems: application to (15)N relaxations in (15)NH4(+)

The equations that describe the time-evolution of transverse and longitudinal (15)N magnetisations in tetrahedral ammonium ions, (15)NH4(+), are derived from the Bloch-Wangsness-Redfield density operator relaxation theory. It is assumed that the relaxation of the spin-states is dominated by (1) the intra-molecular (15)N-(1)H and (1)H-(1)H dipole-dipole interactions and (2) interactions of the ammonium protons with remote spins, which also include the contribution to the relaxations that arise from the exchange of the ammonium protons with the bulk solvent. The dipole-dipole cross-correlated relaxation mechanisms between each of the (15)N-(1)H and (1)H-(1)H interactions are explicitly taken into account in the derivations. An application to (15)N-ammonium bound to a 41kDa domain of the protein DnaK is presented, where a comparison between experiments and simulations show that the ammonium ion rotates rapidly within its binding site with a local correlation time shorter than approximately 1ns. The theoretical framework provided here forms the basis for further investigations of dynamics of AX4 spin systems, with ammonium ions in solution and bound to proteins of particular interest.

Chemical exchange in biomacromolecules: past, present, and future

The perspective reviews quantitative investigations of chemical exchange phenomena in proteins and other biological macromolecules using NMR spectroscopy, particularly relaxation dispersion methods. The emphasis is on techniques and applications that quantify the populations, interconversion kinetics, and structural features of sparsely populated conformational states in equilibrium with a highly populated ground state. Applications to folding, molecular recognition, catalysis, and allostery by proteins and nucleic acids are highlighted.

Allosteric coupling between the intracellular coupling helix 4 and regulatory sites of the first nucleotide-binding domain of CFTR

Cystic fibrosis is caused by mutations in CFTR (cystic fibrosis transmembrane conductance regulator), leading to folding and processing defects and to chloride channel gating misfunction. CFTR is regulated by ATP binding to its cytoplasmic nucleotide-binding domains, NBD1 and NBD2, and by phosphorylation of the NBD1 regulatory insert (RI) and the regulatory extension (RE)/R region. These regulatory effects are transmitted to the rest of the channel via NBD interactions with intracellular domain coupling helices (CL), particularly CL4. Using a sensitive method for detecting inter-residue correlations between chemical shift changes in NMR spectra, an allosteric network was revealed within NBD1, with a construct lacking RI. The CL4-binding site couples to the RI-deletion site and the C-terminal residues of NBD1 that precede the R region in full-length CFTR. Titration of CL4 peptide into NBD1 perturbs the conformational ensemble in these sites with similar titration patterns observed in F508del, the major CF-causing mutant, and in suppressor mutants F494N, V510D and Q637R NBD1, as well as in a CL4-NBD1 fusion construct. Reciprocally, the C-terminal mutation, Q637R, perturbs dynamics in these three sites. This allosteric network suggests a mechanism synthesizing diverse regulatory NBD1 interactions and provides biophysical evidence for the allosteric coupling required for CFTR function.

Enhanced accuracy of kinetic information from CT-CPMG experiments by transverse rotating-frame spectroscopy

Micro-to-millisecond motions of proteins transmit pivotal signals for protein function. A powerful technique for the measurement of these motions is nuclear magnetic resonance spectroscopy. One of the most widely used methodologies for this purpose is the constant-time Carr-Purcell-Meiboom-Gill (CT-CPMG) relaxation dispersion experiment where kinetic and structural information can be obtained at atomic resolution. Extraction of accurate kinetics determined from CT-CPMG data requires refocusing frequencies that are much larger than the nuclei's exchange rate between states. We investigated the effect when fast processes are probed by CT-CPMG experiments via simulation and show that if the intrinsic relaxation rate (R(CT-CPMG)(2,0)) is not known a priori the extraction of accurate kinetics is hindered. Errors on the order of 50 % in the exchange rate are attained when processes become fast, but are minimized to 5 % with a priori (CT-CPMG)(2,0)) information. To alleviate this shortcoming, we developed an experimental scheme probing (CT-CPMG)(2,0)) with large amplitude spin-lock fields, which specifically contains the intrinsic proton longitudinal Eigenrelaxation rate. Our approach was validated with ubiquitin and the Oscillatoria agardhii agglutinin (OAA). For OAA, an underestimation of 66 % in the kinetic rates was observed if (CT-CPMG)(2,0)) is not included during the analysis of CT-CPMG data and result in incorrect kinetics and imprecise amplitude information. This was overcome by combining CT-CPMG with (CT-CPMG)(2,0)) measured with a high power R1ρ experiment. In addition, the measurement of (CT-CPMG)(2,0)) removes the ambiguities in choosing between different models that describe CT-CPMG data.

Conformational stabilization of ubiquitin yields potent and selective inhibitors of USP7

Protein conformation and function are often inextricably linked, such that the states a protein adopts define its enzymatic activity or its affinity for various partners. Here we combine computational design with macromolecular display to isolate functional conformations of ubiquitin that tightly bind the catalytic core of the oncogenic ubiquitin-specific protease 7 (USP7) deubiquitinase. Structural and biochemical characterization of these ubiquitin variants suggest that remodeled backbone conformations and core packing poise these molecules for stronger interactions, leading to potent and specific inhibition of enzymatic activity. A ubiquitin variant expressed in human tumor cell lines binds and inhibits endogenous USP7, thereby enhancing Mdm2 proteasomal turnover and stabilizing p53. In sum, we have developed an approach to rationally target macromolecular libraries toward the remodeling of protein conformation, shown that engineering of ubiquitin conformation can greatly increase its interaction with deubiquitinases and developed powerful tools to probe the cellular role of USP7.

Is buffer a good proxy for a crowded cell-like environment? A comparative NMR study of calmodulin side-chain dynamics in buffer and E. coli lysate

Biophysical studies of protein structure and dynamics are typically performed in a highly controlled manner involving only the protein(s) of interest. Comparatively fewer such studies have been carried out in the context of a cellular environment that typically involves many biomolecules, ions and metabolites. Recently, solution NMR spectroscopy, focusing primarily on backbone amide groups as reporters, has emerged as a powerful technique for investigating protein structure and dynamics in vivo and in crowded “cell-like” environments. Here we extend these studies through a comparative analysis of Ile, Leu, Val and Met methyl side-chain motions in apo, Ca²⁺-bound and Ca²⁺, peptide-bound calmodulin dissolved in aqueous buffer or in E. coli lysate. Deuterium spin relaxation experiments, sensitive to pico- to nano-second time-scale processes and Carr-Purcell-Meiboom-Gill relaxation dispersion experiments, reporting on millisecond dynamics, have been recorded. Both similarities and differences in motional properties are noted for calmodulin dissolved in buffer or in lysate. These results emphasize that while significant insights can be obtained through detailed “test-tube” studies, experiments performed under conditions that are “cell-like” are critical for obtaining a comprehensive understanding of protein motion in vivo and therefore for elucidating the relation between motion and function.

Local folding and misfolding in the PBX homeodomain from a three-state analysis of CPMG relaxation dispersion NMR data

NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments represent a powerful approach for characterizing protein internal motions and for gaining insight into fundamental biological processes such as protein folding, catalysis, and allostery. In most cases, CPMG data are analyzed assuming that the protein exchanges between two different conformational states. Systems exchanging among more than two states are far more challenging to characterize by CPMG NMR. For example, in the case of three-state exchange in the fast time scale regime, it is difficult to uniquely connect the parameters extracted from CPMG analyses with the physical parameters of most interest, intercoversion rates, populations, and chemical shift differences for exchanging states. We have developed a grid search selection procedure that allows these physical parameters to be uniquely determined from CPMG data, based on additional information, which in this study comprises ligand-induced chemical shift perturbations. We applied this approach to the PBX homeodomain (PBX-HD), a three-helix protein with a C-terminal extension that folds into a fourth helix upon binding to DNA. We recently showed that the C-terminal extension transiently folds, even in the absence DNA, in a process that is likely tied to the cooperative binding of PBX-HD to DNA and other homeodomains. Using the grid search selection procedure, we found that PBX-HD undergoes exchange between three different conformational states, a major form in which the C-terminal extension is unfolded, the previously identified state in which the C-terminal extension forms a fourth helix, and an additional state in which the C-terminal extension is misfolded.

Increasing the exchange time-scale that can be probed by CPMG relaxation dispersion NMR

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy has emerged as a valuable tool to characterize conformational exchange between major and minor states in a large variety of biomolecules. The window of exchange that is amenable for study, corresponding to rates on the order of 2000 s(-1) or less, is limiting, however. Here we show that a combined analysis of both amide (15)N and (1)H(N) CPMG profiles and major state exchange induced (15)N chemical shift changes leads to significant increases in the exchange time scale for which accurate exchange parameters and chemical shift differences between the interconverting states can be obtained. The utility of the approach is illustrated with examples involving a pair of protein systems that are in the moderately fast exchange regime. In these cases the analysis of dispersion profiles alone is not sufficient to obtain robust measures of exchange parameters and chemical shift differences. Inclusion of major state exchange induced (15)N chemical shift changes measured in ((15)N-(1)H(N)) HMQC and HSQC data sets in addition to the (15)N and (1)H(N) dispersion profiles in the analysis "breaks" the correlation in parameters, allowing accurate values to be obtained. The approach is straightforward to implement and makes use of HMQC/HSQC data sets that are recorded as a matter of routine to obtain chemical shifts of the excited state. It promises to increase the range of exchanging systems involving low populated, transiently formed excited states that can be studied by relaxation dispersion NMR.

Active site dynamics in NADH oxidase from Thermus thermophilus studied by NMR spin relaxation

We have characterized the backbone dynamics of NADH oxidase from Thermus thermophilus (NOX) using a recently-developed suite of NMR experiments designed to isolate exchange broadening, together with (15)N R (1), R (1ρ ), and {(1)H}-(15)N steady-state NOE relaxation measurements performed at 11.7 and 18.8 T. NOX is a 54 kDa homodimeric enzyme that belongs to a family of structurally homologous flavin reductases and nitroreductases with many potential biotechnology applications. Prior studies have suggested that flexibility is involved in the catalytic mechanism of the enzyme. The active site residue W47 was previously identified as being particularly important, as its level of solvent exposure correlates with enzyme activity, and it was observed to undergo "gating" motions in computer simulations. The NMR data are consistent with these findings. Signals from W47 are dynamically broadened beyond detection and several other residues in the active site have significant R ( ex ) contributions to transverse relaxation rates. In addition, the backbone of S193, whose side chain hydroxyl proton hydrogen bonds directly with the FMN cofactor, exhibits extensive mobility on the ns-ps timescale. We hypothesize that these motions may facilitate structural rearrangements of the active site that allow NOX to accept both FMN and FAD as cofactors.

Recent developments in (15)N NMR relaxation studies that probe protein backbone dynamics

Nuclear Magnetic Resonance (NMR) relaxation is a powerful technique that provides information about internal dynamics associated with configurational energetics in proteins, as well as site-specific information involved in conformational equilibria. In particular, (15)N relaxation is a useful probe to characterize overall and internal backbone dynamics of proteins because the relaxation mainly reflects reorientational motion of the N-H bond vector. Over the past 20 years, experiments and protocols for analysis of (15)N R (1), R 2, and the heteronuclear (15)N-{(1)H} NOE data have been well established. The development of these methods has kept pace with the increase in the available static-magnetic field strength, providing dynamic parameters optimized from data fitting at multiple field strengths. Using these methodological advances, correlation times for global tumbling and order parameters and correlation times for internal motions of many proteins have been determined. More recently, transverse relaxation dispersion experiments have extended the range of NMR relaxation studies to the milli- to microsecond time scale, and have provided quantitative information about functional conformational exchange in proteins. Here, we present an overview of recent advances in (15)N relaxation experiments to characterize protein backbone dynamics.

Mechanism for recognition of polyubiquitin chains: balancing affinity through interplay between multivalent binding and dynamics

RAP80 plays a key role in signal transduction in the DNA damage response by recruiting proteins to DNA damage foci by binding K63-polyubiquitin chains with two tandem ubiquitin-interacting motifs (tUIM). It is generally recognized that the typically weak interaction between ubiquitin (Ub) and various recognition motifs is intensified by themes such as tandem recognition motifs and Ub polymerization to achieve biological relevance. However, it remains an intricate problem to develop a detailed molecular mechanism to describe the process that leads to amplification of the Ub signal. A battery of solution-state NMR methods and molecular dynamics simulations were used to demonstrate that RAP80-tUIM employs mono- and multivalent interactions with polyUb chains to achieve enhanced affinity in comparison to monoUb interactions for signal amplification. The enhanced affinity is balanced by unfavorable entropic effects that include partial quenching of rapid reorientation between individual UIM domains and individual Ub domains in the bound state. For the RAP80-tUIM-polyUb interaction, increases in affinity with increasing chain length are a result of increased numbers of mono- and multivalent binding sites in the longer polyUb chains. The mono- and multivalent interactions are characterized by intrinsically weak binding and fast off-rates; these weak interactions with fast kinetics may be an important factor underlying the transient nature of protein-protein interactions that comprise DNA damage foci.

Protein flexibility and conformational entropy in ligand design targeting the carbohydrate recognition domain of galectin-3

Rational drug design is predicated on knowledge of the three-dimensional structure of the protein-ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. (15)N and (2)H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein-carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.

Simple tests for the validation of multiple field spin relaxation data

(15)N spin relaxation data is widely used to extract detailed dynamic information regarding bond vectors such as the amide N-H bond of the protein backbone. Analysis is typically carried using the Lipari-Szabo model-free approach. Even though the original model-free equation can be determined from single field R (1), R (2) and NOE, over-determination of more complex motional models is dependent on the recording of multiple field datasets. This is especially important for the characterization of conformational exchange which affects R (2) in a field dependent manner. However, severe artifacts can be introduced if inconsistencies arise between experimental setups with different magnets (or samples). Here, we propose the use of simple tests as validation tools for the assessment of consistency between different datasets recorded at multiple magnetic fields. Synthetic data are used to show the effects of inconsistencies on the proposed tests. Moreover, an analysis of data currently deposited in the BMRB is performed. Finally, two cases from our laboratory are presented. These tests are implemented in the open-source program relax, and we propose their use as a routine check-up for assessment of multiple field dataset consistency prior to any analysis such as model-free calculations. We believe this will aid in the extraction of higher quality dynamics information from (15)N spin relaxation data.

Selective characterization of microsecond motions in proteins by NMR relaxation

The three-dimensional structures of macromolecules fluctuate over a wide range of time-scales. Separating the individual dynamic processes according to frequency is of importance in relating protein motions to biological function and stability. We present here a general NMR method for the specific characterization of microsecond motions at backbone positions in proteins even in the presence of other dynamics such as large-amplitude nanosecond motions and millisecond chemical exchange processes. The method is based on measurement of relaxation rates of four bilinear coherences and relies on the ability of strong continuous radio frequency fields to quench millisecond chemical exchange. The utility of the methodology is demonstrated and validated through two specific examples focusing on the thermo-stable proteins, ubiquitin and protein L, where it is found that small-amplitude microsecond dynamics are more pervasive than previously thought. Specifically, these motions are localized to alpha helices, loop regions, and regions along the rim of beta sheets in both of the proteins examined. A third example focuses on a 28 kDa ternary complex of the chaperone Chz1 and the histones H2A.Z/H2B, where it is established that pervasive microsecond motions are localized to a region of the chaperone that is important for stabilizing the complex. It is further shown that these motions can be well separated from extensive millisecond dynamics that are also present and that derive from exchange of Chz1 between bound and free states. The methodology is straightforward to implement, and data recorded at only a single static magnetic field are required.

Deuterium spin probes of backbone order in proteins: 2H NMR relaxation study of deuterated carbon alpha sites

(2)H spin relaxation NMR experiments to study the dynamics of deuterated backbone alpha-positions, D(alpha), are developed. To date, solution-state (2)H relaxation measurements in proteins have been confined to side-chain deuterons-primarily (13)CH(2)D or (13)CHD(2) methyl groups. It is shown that quantification of (2)H relaxation rates at D(alpha) backbone positions and the derivation of associated order parameters of C(alpha)-D(alpha) bond vector motions in small [U-(15)N,(13)C,(2)H]-labeled proteins is feasible with reasonable accuracy. The utility of the developed methodology is demonstrated on a pair of proteins-ubiquitin (8.5 kDa) at 10, 27, and 40 degrees C, and a variant of GB1 (6.5 kDa) at 22 degrees C. In both proteins, the D(alpha)-derived parameters of the global rotational diffusion tensor are in good agreement with those obtained from (15)N relaxation rates. Semiquantitative solution-state NMR measurements yield an average value of the quadrupolar coupling constant, QCC, for D(alpha) sites in proteins equal to 174 kHz. Using a uniform value of QCC for all D(alpha) sites, we show that C(alpha)-D(alpha) bond vectors are motionally distinct from the backbone amide N-H bond vectors, with (2)H-derived squared order parameters of C(alpha)-D(alpha) bond vector motions, S(2)(CalphaDalpha), on average slightly higher than their N-H amides counterparts, S(2)(NH). For ubiquitin, the (2)H-derived backbone mobility compares well with that found in a 1-mus molecular dynamics simulation.