Methyl-Based NMR Spectroscopy Methods for Uncovering Structural Dynamics in Large Proteins and Protein Complexes

Synthesis, NMR kinetics and dynamic structure of a 17-mer heptaloop RNA hairpin carrying a 3-N-methyluridine nucleotide residue in the loop region

A 17-mer RNA hairpin (5'GGGAGUXAGCGGCUCCC3') carrying 3-N-methyluridine (m3U) at position X (m3U7-RNA), designed to represent the anticodon stem-loop (ACSL) region of tRNAs to study an open loop state (O-state), was synthesized, purified by HPLC, and characterized by MALDI-ToF_MS and NMR methods. 1H-NMR data revealed primary (P-state in 56.1%), secondary (S-state in 43.9%) and tertiary (∼5-6%) ACSL conformations. Exchange rate constant (kex) for interconversion between P and S states is 112 sec-1 (<Δω ∼454 rad/sec), confirming a slow exchange regime between the two states. Forward (kPS) and backward (kSP) rate constants are 49.166 sec-1 and 62.792 sec-1, respectively, leading to a longer life-time (20.339 msec) for the P-state and a shorter life-time (15.926 msec) for the S-state. In accordance with conformational populations determined by 1H-NMR, dynamics of the P/S/tertiary states of m3U7-RNA and its wild-type counterpart (wt-RNA) were studied by three independent MD production simulations. Cluster analysis revealed that wt-RNA reflects the structural characteristics of the ACSL region of tRNAs. The P-state of m3U7-RNA was found to be structurally similar to wt-RNA but lacks an intraloop H-bond between m3U7 and C10 (U33 and nt36 in tRNAs). In the S-state of m3U7-RNA, m3U7 flips out of the loop region. O-state loop conformations of m3U7-RNA were also clustered (4.8%), wherein the loop nucleotides m3U7.A8.G9.C10.G11 stack one after another. We propose that the O-state of m3U7-RNA is the most suitable conformation that makes the loop accessible for complementary nucleotides and for non-enzymatic primordial replication of small circular RNAs.Communicated by Ramaswamy H. Sarma.

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

Cell-free synthesis of proteins with selectively 13C-labelled methyl groups from inexpensive precursors

The novel eCell system maintains the activity of the entire repertoire of metabolic Escherichia coli enzymes in cell-free protein synthesis. We show that this can be harnessed to produce proteins with selectively 13 C-labelled amino acids from inexpensive 13 C-labelled precursors. The system is demonstrated with selective 13 C labelling of methyl groups in the proteins ubiquitin and peptidyl-prolyl cis-trans isomerase B. Starting from 3-13 C-pyruvate, 13 C-HSQC cross-peaks are obtained devoid of one-bond 13 C-13 C scalar couplings. Starting from 2-13 C-methyl-acetolactate, single methyl groups of valine and leucine are labelled. Labelling efficiencies are 70 % or higher, and the method allows us to produce perdeuterated proteins with protonated methyl groups in a residue-selective manner. The system uses the isotope-labelled precursors sparingly and is readily scalable.

AI-designed NMR spectroscopy RF pulses for fast acquisition at high and ultra-high magnetic fields

Nuclear magnetic resonance (NMR) spectroscopy is a powerful high-resolution tool for characterizing biomacromolecular structure, dynamics, and interactions. However, the lengthy longitudinal relaxation of the nuclear spins significantly extends the total experimental time, especially at high and ultra-high magnetic field strengths. Although longitudinal relaxation-enhanced techniques have sped up data acquisition, their application has been limited by the chemical shift dispersion. Here we combined an evolutionary algorithm and artificial intelligence to design ¹H and ¹⁵N radio frequency (RF) pulses with variable phase and amplitude that cover significantly broader bandwidths and allow for rapid data acquisition. We re-engineered the basic transverse relaxation optimized spectroscopy experiment and showed that the RF shapes enhance the spectral sensitivity of well-folded proteins up to 180 kDa molecular weight. These RF shapes can be tailored to re-design triple-resonance experiments for accelerating NMR spectroscopy of biomacromolecules at high fields.

Biochemical methods to map and quantify allosteric motions in human glucokinase

Allosteric regulation of protein function is ubiquitous in biology. Allostery originates from ligand-mediated alterations in polypeptide structure and/or dynamics, which produce a cooperative kinetic or thermodynamic response to changing ligand concentrations. Establishing a mechanistic description of individual allosteric events requires both mapping the relevant changes in protein structure and quantifying the rates of differential conformational dynamics in the absence and presence of effectors. In this chapter, we describe three biochemical approaches to understand the dynamic and structural signatures of protein allostery using the well-established cooperative enzyme glucokinase as a case study. The combined application of pulsed proteolysis, biomolecular nuclear magnetic resonance spectroscopy and hydrogen-deuterium exchange mass spectrometry offers complementary information that can used to establish molecular models for allosteric proteins, especially when differential protein dynamics are involved.

Slow methyl axes motions in perdeuterated villin headpiece subdomain probed by cross-correlated NMR relaxation measurements

Protein methyl groups can participate in multiple motional modes on different time scales. Sub-nanosecond to nano-second time scale motions of methyl axes are particularly challenging to detect for small proteins in solutions. In this work we employ NMR relaxation interference between the methyl H-H/H-C dipole-dipole interactions [Sun&Tugarinov, J. Magn. Reason. 2012] to characterize methyl axes motions as a function of temperature in a small model protein villin headpiece subdomain (HP36), in which all non-exchangeable protons are deuterated with the exception of methyl groups of leucine and valine residues. The data points to the existence of slow motional modes of methyl axes on sub-nanosecond to nanosecond time scales. Further, at high temperatures for which the overall tumbling of the protein is on the order of 2 ns, we observe a coupling between the slow internal motion and the overall molecular tumbling, based on the anomalous order parameters and their temperature-dependent trends. The addition of 28%(w/w) glycerol-d8 increases the viscosity of the solvent and separates the timescales of internal and overall tumbling, thus permitting for another view of the necessity of the coupling assumption for these sites at high temperatures.

Motional clustering in supra-τc conformational exchange influences NOE cross-relaxation rate

Biomolecular spin relaxation processes, such as the NOE, are commonly modeled by rotational τ_c-tumbling combined with fast motions on the sub-τ_c timescale. Motions on the supra-τ_c timescale, in contrast, are considered to be completely decorrelated to the molecular tumbling and therefore invisible. Here, we show how supra-τ_c dynamics can nonetheless influence the NOE build-up between methyl groups. This effect arises because supra-τ_c motions can cluster the fast-motion ensembles into discrete states, affecting distance averaging as well as the fast-motion order parameter and hence the cross-relaxation rate. We present a computational approach to estimate methyl-methyl cross-relaxation rates from extensive (>100×τ_c) all-atom molecular dynamics (MD) trajectories on the example of the 723-residue protein Malate Synthase G. The approach uses Markov state models (MSMs) to resolve transitions between metastable states and thus to discriminate between sub-τ_c and supra-τ_c conformational exchange. We find that supra-τ_c exchange typically increases NOESY cross-peak intensities. The methods described in this work extend the theory of modeling sub-μs dynamics in spin relaxation and thus contribute to a quantitative estimation of NOE cross-relaxation rates from MD simulations, eventually leading to increased precision in structural and functional studies of large proteins.

NMR Assignment of Methyl Groups in Immobilized Proteins Using Multiple-Bond 13C Homonuclear Transfers, Proton Detection, and Very Fast MAS

In nuclear magnetic resonance spectroscopy of proteins, methyl protons play a particular role as extremely sensitive reporters on dynamics, allosteric effects, and protein–protein interactions, accessible even in high-molecular-weight systems approaching 1 MDa. The notorious issue of their chemical shift assignment is addressed here by a joint use of solid-state 1H-detected methods at very fast (nearly 100 kHz) magic-angle spinning, partial deuteration, and high-magnetic fields. The suitability of a series of RF schemes is evaluated for the efficient coherence transfer across entire 13C side chains of methyl-containing residues, which is key for establishing connection between methyl and backbone 1H resonances. The performance of ten methods for recoupling of either isotropic 13C–13C scalar or anisotropic dipolar interactions (five variants of TOBSY, FLOPSY, DIPSI, WALTZ, RFDR, and DREAM) is evaluated experimentally at two state-of-the-art magic-angle spinning (55 and 94.5 kHz) and static magnetic field conditions (18.8 and 23.5 T). Model isotopically labeled compounds (alanine and Met-Leu-Phe tripeptide) and ILV-methyl and amide-selectively protonated, and otherwise deuterated chicken α-spectrin SH3 protein are used as convenient reference systems. Spin dynamics simulations in SIMPSON are performed to determine optimal parameters of these RF schemes, up to recently experimentally attained spinning frequencies (200 kHz) and B0 field strengths (28.2 T). The concept of linearization of 13C side chain by appropriate isotope labeling is revisited and showed to significantly increase sensitivity of methyl-to-backbone correlations. A resolution enhancement provided by 4D spectroscopy with non-uniform (sparse) sampling is demonstrated to remove ambiguities in simultaneous resonance assignment of methyl proton and carbon chemical shifts.

Bacterial Transcriptional Regulators: A Road Map for Functional, Structural, and Biophysical Characterization

The different niches through which bacteria move during their life cycle require a fast response to the many environmental queues they encounter. The sensing of these stimuli and their correct response is driven primarily by transcriptional regulators. This kind of protein is involved in sensing a wide array of chemical species, a process that ultimately leads to the regulation of gene transcription. The allosteric-coupling mechanism of sensing and regulation is a central aspect of biological systems and has become an important field of research during the last decades. In this review, we summarize the state-of-the-art techniques applied to unravel these complex mechanisms. We introduce a roadmap that may serve for experimental design, depending on the answers we seek and the initial information we have about the system of study. We also provide information on databases containing available structural information on each family of transcriptional regulators. Finally, we discuss the recent results of research about the allosteric mechanisms of sensing and regulation involving many transcriptional regulators of interest, highlighting multipronged strategies and novel experimental techniques. The aim of the experiments discussed here was to provide a better understanding at a molecular level of how bacteria adapt to the different environmental threats they face.

Analysis of conformational exchange processes using methyl-TROSY-based Hahn echo measurements of quadruple-quantum relaxation

Transverse nuclear spin relaxation is a sensitive probe of chemical exchange on timescales on the order of microseconds to milliseconds. Here we present an experiment for the simultaneous measurement of the relaxation rates of two quadruple-quantum transitions in 13 CH3 -labelled methyl groups. These coherences are protected against relaxation by intra-methyl dipolar interactions and so have unexpectedly long lifetimes within perdeuterated biomacromolecules. However, these coherences also have an order of magnitude higher sensitivity to chemical exchange broadening than lower order coherences and therefore provide ideal probes of dynamic processes. We show that analysis of the static magnetic field dependence of zero-, double- and quadruple-quantum Hahn echo relaxation rates provides a robust indication of chemical exchange and can determine the signed relative magnitudes of proton and carbon chemical shift differences between ground and excited states. We also demonstrate that this analysis can be combined with established Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements, providing improved precision in parameter estimates, particularly in the determination of 1 H chemical shift differences.

Observation of Sub-Microsecond Protein Methyl-Side Chain Dynamics by Nanoparticle-Assisted NMR Spin Relaxation

Amino-acid side-chain properties in proteins are key determinants of protein function. NMR spin relaxation of side chains is an important source of information about local protein dynamics and flexibility. However, traditional solution NMR relaxation methods are most sensitive to sub-nanosecond dynamics lacking information on slower ns-μs time-scale motions. Nanoparticle-assisted NMR spin relaxation (NASR) of methyl-side chains is introduced here as a window into these ns-μs dynamics. NASR utilizes the transient and nonspecific interactions between folded proteins and slowly tumbling spherical nanoparticles (NPs), whereby the increase of the relaxation rates reflects motions on time scales from ps all the way to the overall tumbling correlation time of the NPs ranging from hundreds of ns to μs. The observed motional amplitude of each methyl group can then be expressed by a model-free NASR S² order parameter. The method is demonstrated for ²H-relaxation of CH₂D methyl moieties and cross-correlated relaxation of CH₃ groups for proteins Im7 and ubiquitin in the presence of anionic silica-nanoparticles. Both types of relaxation experiments, dominated by either quadrupolar or dipolar interactions, yield highly consistent results. Im7 shows additional dynamics on the intermediate time scales taking place in a functionally important loop, whereas ubiquitin visits the majority of its conformational substates on the sub-ns time scale. These experimental observations are in good agreement with 4-10 μs all-atom molecular dynamics trajectories. NASR probes side-chain dynamics on a much wider range of motional time scales than previously possible, thereby providing new insights into the interplay between protein structure, dynamics, and molecular interactions that govern protein function.

Allosteric and dynamic control of RNA-dependent RNA polymerase function and fidelity

All RNA viruses encode an RNA-dependent RNA polymerase (RdRp) responsible for genome replication. It is now recognized that enzymes in general, and RdRps specifically, are dynamic macromolecular machines such that their moving parts, including active site loops, play direct functional roles. While X-ray crystallography has provided deep insight into structural elements important for RdRp function, this methodology generally provides only static snapshots, and so is limited in its ability to report on dynamic fluctuations away from the lowest energy conformation. Nuclear magnetic resonance (NMR), molecular dynamics (MD) simulations and other biophysical techniques have brought new insight into RdRp function by their ability to characterize the trajectories, kinetics and thermodynamics of conformational motions. In particular, these methodologies have identified coordinated motions among conserved structural motifs necessary for nucleotide selection and incorporation. Disruption of these motions through amino acid substitutions or inhibitor binding impairs RdRp function. Understanding and re-engineering these motions thus provides exciting new avenues for anti-viral strategies. This chapter outlines the basics of these methodologies, summarizes the dynamic motions observed in different RdRps important for nucleotide selection and incorporation, and illustrates how this information can be leveraged towards rational vaccine strain development and anti-viral drug design.

Methodological advancements for characterising protein side chains by NMR spectroscopy

The surface of proteins is covered by side chains of polar amino acids that are imperative for modulating protein functionality through the formation of noncovalent intermolecular interactions. However, despite their tremendous importance, the unique structures of protein side chains require tailored approaches for investigation by nuclear magnetic resonance spectroscopy and so have traditionally been understudied compared with the protein backbone. Here, we review substantial recent methodological advancements within nuclear magnetic resonance spectroscopy to address this issue. Specifically, we consider advancements that provide new insight into methyl-bearing side chains, show the potential of using non-natural amino acids and reveal the actions of charged side chains. Combined, the new methods promise unprecedented characterisations of side chains that will further elucidate protein function.

Simultaneous measurement of 1HC/N-R2's for rapid acquisition of backbone and sidechain paramagnetic relaxation enhancements (PREs) in proteins

Paramagnetic relaxation enhancements (PREs) are routinely used to provide long-range distance restraints for the determination of protein structures, to resolve protein dynamics, ligand-protein binding sites, and lowly populated species, using Nuclear Magnetic Resonance Spectroscopy (NMR). Here, we propose a simultaneous ¹H-¹⁵ N, ¹H-¹³C SESAME based pulse scheme for the rapid acquisition of ¹H^C/N-R₂ relaxation rates for the determination of backbone and sidechain PREs of proteins. The ¹H^N-R₂ rates from the traditional and our approach on Ubiquitin (UBQ) are well correlated (R² = 0.99), revealing their potential to be used quantitatively. Comparison of the S57C UBQ calculated and experimental PREs provided backbone and side chain Q factors of 0.23 and 0.24, respectively, well-fitted to the UBQ NMR structure, showing that our approach can be used to acquire accurate PRE rates from the functionally important sites of proteins but in at least half the time as traditional methods.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

The Role of Conformational Dynamics in the Recognition and Regulation of Ubiquitination

The ubiquitination pathway is central to many cell signaling and regulatory events. One of the intriguing aspects of the pathway is the combinatorial sophistication of substrate recognition and ubiquitin chain building determinations. The abundant structural and biological data portray several characteristic protein folds among E2 and E3 proteins, and the understanding of the combinatorial complexity that enables interaction with much of the human proteome is a major goal to developing targeted and selective manipulation of the pathway. With the commonality of some folds, there are likely other aspects that can provide differentiation and recognition. These aspects involve allosteric effects and conformational dynamics that can direct recognition and chain building processes. In this review, we will describe the current state of the knowledge for conformational dynamics across a wide timescale, address the limitations of present approaches, and illustrate the potential to make new advances in connecting dynamics with ubiquitination regulation.

The dynamic nature of the Mre11-Rad50 DNA break repair complex

The Mre11-Rad50-Nbs1/Xrs2 protein complex plays a pivotal role in the detection and repair of DNA double strand breaks. Through traditional and emerging structural biology techniques, various functional structural states of this complex have been visualized; however, relatively little is known about the transitions between these states. Indeed, it is these structural transitions that are important for Mre11-Rad50-mediated DNA unwinding at a break and the activation of downstream repair signaling events. Here, we present a brief overview of the current understanding of the structure of the core Mre11-Rad50 complex. We then highlight our recent studies emphasizing the contributions of solution state NMR spectroscopy and other biophysical techniques in providing insight into the structures and dynamics associated with Mre11-Rad50 functions.

Recent advances in solid-state relaxation dispersion techniques

This review describes two rotating-frame (R_1ρ) relaxation dispersion methods, namely the Bloch-McConnell Relaxation Dispersion and the Near-rotary Resonance Relaxation Dispersion, which enable the study of microsecond time-scale conformational fluctuations in the solid state using magic-angle-spinning nuclear magnetic resonance spectroscopy. The goal is to provide the reader with key ideas, experimental descriptions, and practical considerations associated with R_1ρ measurements that are needed for analyzing relaxation dispersion and quantifying conformational exchange. While the focus is on protein motion, many presented concepts can be equally well adapted to study the microsecond time-scale dynamics of other bio- (e.g. lipids, polysaccharides, nucleic acids), organic (e.g. pharmaceutical compounds), or inorganic molecules (e.g., metal organic frameworks). This article summarizes the essential contributions made by recent theoretical and experimental solid-state NMR studies to our understanding of protein motion. Here we discuss recent advances in fast MAS applications that enable the observation and atomic level characterization of sparsely populated conformational states which are otherwise inaccessible for other experimental methods. Such high-energy states are often associated with protein functions such as molecular recognition, ligand binding, or enzymatic catalysis, as well as with disease-related properties such as misfolding and amyloid formation.

Automated assignment of methyl NMR spectra from large proteins

As structural biology trends towards larger and more complex biomolecular targets, a detailed understanding of their interactions and underlying structures and dynamics is required. The development of methyl-TROSY has enabled NMR spectroscopy to provide atomic-resolution insight into the mechanisms of large molecular assemblies in solution. However, the applicability of methyl-TROSY has been hindered by the laborious and time-consuming resonance assignment process, typically performed with domain fragmentation, site-directed mutagenesis, and analysis of NOE data in the context of a crystal structure. In response, several structure-based automatic methyl assignment strategies have been developed over the past decade. Here, we present a comprehensive analysis of all available methods and compare their input data requirements, algorithmic strategies, and reported performance. In general, the methods fall into two categories: those that primarily rely on inter-methyl NOEs, and those that utilize methyl PRE- and PCS-based restraints. We discuss their advantages and limitations, and highlight the potential benefits from standardizing and combining different methods.

Mutation of Conserved Mre11 Residues Alter Protein Dynamics to Separate Nuclease Functions

Naked and protein-blocked DNA ends occur naturally during immune cell development, meiosis, and at telomeres as well as from aborted topoisomerase reactions, collapsed replication forks, and other stressors. Damaged DNA ends are dangerous in cells and if left unrepaired can lead to genomic rearrangement, loss of genetic information, and eventually cancer. Mre11 is part of the Mre11-Rad50-Nbs1 complex that recognizes DNA double-strand breaks and has exonuclease and endonuclease activities that help to initiate the repair processes to resolve these broken DNA ends. In fact, these activities are crucial for proper DNA damage repair pathway choice. Here, using Pyrococcus furiosus Mre11, we question how two Mre11 separation-of-function mutants, one previously described but the second first described here, maintain endonuclease activity in the absence of exonuclease activity. To start, we performed solution-state NMR experiments to assign the side-chain methyl groups of the 64-kDa Mre11 nuclease and capping domains, which allowed us to describe the structural differences between Mre11 bound to exo- and endonuclease substrates. Then, through biochemical and biophysical characterization, including NMR structural and dynamics studies, we compared the two mutants and determined that both affect the dynamic features and double-stranded DNA binding properties of Mre11, but in different ways. In total, our results illuminate the structural and dynamic landscape of Mre11 nuclease function.

Automatic structure-based NMR methyl resonance assignment in large proteins

Isotopically labeled methyl groups provide NMR probes in large, otherwise deuterated proteins. However, the resonance assignment constitutes a bottleneck for broader applicability of methyl-based NMR. Here, we present the automated MethylFLYA method for the assignment of methyl groups that is based on methyl-methyl nuclear Overhauser effect spectroscopy (NOESY) peak lists. MethylFLYA is applied to five proteins (28–358 kDa) comprising a total of 708 isotope-labeled methyl groups, of which 612 contribute NOESY cross peaks. MethylFLYA confidently assigns 488 methyl groups, i.e. 80% of those with NOESY data. Of these, 459 agree with the reference, 6 were different, and 23 were without reference assignment. MethylFLYA assigns significantly more methyl groups than alternative algorithms, has an average error rate of 1%, modest runtimes of 0.4–1.2 h, and can handle arbitrary isotope labeling patterns and data from other types of NMR spectra.

The structures and dynamics of large proteins can be studied with methyl-based NMR but peak assignment is still challenging. Here the authors present MethylFLYA that allows automated assignment of methyl groups and apply it to five proteins with molecular weights in the range from 28 to 358 kDa.