Asymmetry of 13C labeled 3-pyruvate affords improved site specific labeling of RNA for NMR spectroscopy

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Isotope-Labeled RNA Building Blocks for NMR Structure and Dynamics Studies

RNA structural research lags behind that of proteins, preventing a robust understanding of RNA functions. NMR spectroscopy is an apt technique for probing the structures and dynamics of RNA molecules in solution at atomic resolution. Still, RNA analysis by NMR suffers from spectral overlap and line broadening, both of which worsen for larger RNAs. Incorporation of stable isotope labels into RNA has provided several solutions to these challenges. In this review, we summarize the benefits and limitations of various methods used to obtain isotope-labeled RNA building blocks and how they are used to prepare isotope-labeled RNA for NMR structure and dynamics studies.

Quantifying the effects of long-range 13C-13C dipolar coupling on measured relaxation rates in RNA

Selective stable isotope labeling has transformed structural and dynamics analysis of RNA by NMR spectroscopy. These methods can remove ¹³C-¹³C dipolar couplings that complicate ¹³C relaxation analyses. While these phenomena are well documented for sites with adjacent ¹³C nuclei (e.g. ribose C1'), less is known about so-called isolated sites (e.g. adenosine C2). To investigate and quantify the effects of long-range (> 2 Å) ¹³C-¹³C dipolar interactions on RNA dynamics, we simulated adenosine C2 relaxation rates in uniformly [U-¹³C/¹⁵N]-ATP or selectively [2-¹³C]-ATP labeled RNAs. Our simulations predict non-negligible ¹³C-¹³C dipolar contributions from adenosine C4, C5, and C6 to C2 longitudinal (R₁) relaxation rates in [U-¹³C/¹⁵N]-ATP labeled RNAs. Moreover, these contributions increase at higher magnetic fields and molecular weights to introduce discrepancies that exceed 50%. This will become increasingly important at GHz fields. Experimental R₁ measurements in the 61 nucleotide human hepatitis B virus encapsidation signal ε RNA labeled with [U-¹³C/¹⁵N]-ATP or [2-¹³C]-ATP corroborate these simulations. Thus, in the absence of selectively labeled samples, long-range ¹³C-¹³C dipolar contributions must be explicitly taken into account when interpreting adenosine C2 R₁ rates in terms of motional models for large RNAs.

Deleterious effects of carbon-carbon dipolar coupling on RNA NMR dynamics

Many regulatory RNAs undergo dynamic exchanges that are crucial for their biological functions and NMR spectroscopy is a versatile tool for monitoring dynamic motions of biomolecules. Meaningful information on biomolecular dynamics requires an accurate measurement of relaxation parameters such as longitudinal (R₁) rates, transverse (R₂) rates and heteronuclear Overhauser effect (hNOE). However, earlier studies have shown that the large ¹³C-¹³C interactions complicate analysis of the carbon relaxation parameters. To investigate the effect of ¹³C-¹³C interactions on RNA dynamic studies, we performed relaxation measurements on various RNA samples with different labeling patterns and compared these measurements with the computational simulations. For uniformly labeled samples, contributions of the neighboring carbon to R₁ measurements were observed. These contributions increased with increasing magnetic field and overall correlation time ([Formula: see text]) for R₁ rates, necessitating more careful analysis for uniformly labeled large RNAs. In addition, the hNOE measurements were also affected by the adjacent carbon nuclei. Unlike R₁ rates, R_1ρ rates showed relatively good agreement between uniformly- and site-selectively labeled samples, suggesting no dramatic effect from their attached carbon, in agreement with previous observations. Overall, having more accurate rate measurements avoids complex analysis and will be a key for interpreting ¹³C relaxation rates for molecular motion that can provide valuable insights into cellular molecular recognition events.

Benefits of stable isotope labeling in RNA analysis

RNAs are key players in life as they connect the genetic code (DNA) with all cellular processes dominated by proteins. They contain a variety of chemical modifications and many RNAs fold into complex structures. Here, we review recent progress in the analysis of RNA modification and structure on the basis of stable isotope labeling techniques. Mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy are the key tools and many breakthrough developments were made possible by the analysis of stable isotope labeled RNA. Therefore, we discuss current stable isotope labeling techniques such as metabolic labeling, enzymatic labeling and chemical synthesis. RNA structure analysis by NMR is challenging due to two major problems that become even more salient when the size of the RNA increases, namely chemical shift overlaps and line broadening leading to complete signal loss. Several isotope labeling strategies have been developed to provide solutions to these major issues, such as deuteration, segmental isotope labeling or site-specific labeling. Quantification of modified nucleosides in RNA by MS is only possible through the application of stable isotope labeled internal standards. With nucleic acid isotope labeling coupled mass spectrometry (NAIL-MS), it is now possible to analyze the dynamic processes of post-transcriptional RNA modification and demodification. The trend, in both NMR and MS RNA analytics, is without doubt shifting from the analysis of snapshot moments towards the development and application of tools capable of analyzing the dynamics of RNA structure and modification profiles.

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.

Chemical synthesis and NMR spectroscopy of long stable isotope labelled RNA

We showcase the high potential of the 2'-cyanoethoxymethyl (CEM) methodology to synthesize RNAs with naturally occurring modified residues carrying stable isotope (SI) labels for NMR spectroscopic applications. The method was applied to synthesize RNAs with sizes ranging between 60 to 80 nucleotides. The presented approach gives the possibility to selectively modify larger RNAs (>60 nucleotides) with atom-specifically ¹³C/¹⁵N-labelled building blocks. The method harbors the unique potential to address structural as well as dynamic features of these RNAs with NMR spectroscopy but also using other biophysical methods, such as mass spectrometry (MS), or small angle neutron/X-ray scattering (SANS, SAXS).

Applications of NMR to structure determination of RNAs large and small

Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool to investigate the structure and dynamics of RNA, because many biologically important RNAs have conformationally flexible structures, which makes them difficult to crystallize. Functional, independently folded RNA domains, range in size between simple stem-loops of as few as 10-20 nucleotides, to 50-70 nucleotides, the size of tRNA and many small ribozymes, to a few hundred nucleotides, the size of more complex RNA enzymes and of the functional domains of non-coding transcripts. In this review, we discuss new methods for sample preparation, assignment strategies and structure determination for independently folded RNA domains of up to 100 kDa in molecular weight.

In vivo, large-scale preparation of uniformly (15)N- and site-specifically (13)C-labeled homogeneous, recombinant RNA for NMR studies

Knowledge of how ribonucleic acid (RNA) structures fold to form intricate, three-dimensional structures has provided fundamental insights into understanding the biological functions of RNA. Nuclear magnetic resonance (NMR) spectroscopy is a particularly useful high-resolution technique to investigate the dynamic structure of RNA. Effective study of RNA by NMR requires enrichment with isotopes of (13)C or (15)N or both. Here, we present a method to produce milligram quantities of uniformly (15)N- and site-specifically (13)C-labeled RNAs using wild-type K12 and mutant tktA Escherichia coli in combination with a tRNA-scaffold approach. The method includes a double selection protocol to obtain an E. coli clone with consistently high expression of the recombinant tRNA-scaffold. We also present protocols for the purification of the tRNA-scaffold from a total cellular RNA extract and the excision of the RNA of interest from the tRNA-scaffold using DNAzymes. Finally, we showcase NMR applications to demonstrate the benefit of using in vivo site-specifically (13)C-labeled RNA.

(19)F-labeling of the adenine H2-site to study large RNAs by NMR spectroscopy

In comparison to proteins and protein complexes, the size of RNA amenable to NMR studies is limited despite the development of new isotopic labeling strategies including deuteration and ligation of differentially labeled RNAs. Due to the restricted chemical shift dispersion in only four different nucleotides spectral resolution remains limited in larger RNAs. Labeling RNAs with the NMR-active nucleus (19)F has previously been introduced for small RNAs up to 40 nucleotides (nt). In the presented work, we study the natural occurring RNA aptamer domain of the guanine-sensing riboswitch comprising 73 nucleotides from Bacillus subtilis. The work includes protocols for improved in vitro transcription of 2-fluoroadenosine-5'-triphosphat (2F-ATP) using the mutant P266L of the T7 RNA polymerase. Our NMR analysis shows that the secondary and tertiary structure of the riboswitch is fully maintained and that the specific binding of the cognate ligand hypoxanthine is not impaired by the introduction of the (19)F isotope. The thermal stability of the (19)F-labeled riboswitch is not altered compared to the unmodified sequence, but local base pair stabilities, as measured by hydrogen exchange experiments, are modulated. The characteristic change in the chemical shift of the imino resonances detected in a (1)H,(15)N-HSQC allow the identification of Watson-Crick base paired uridine signals and the (19)F resonances can be used as reporters for tertiary and secondary structure transitions, confirming the potential of (19)F-labeling even for sizeable RNAs in the range of 70 nucleotides.

Chemo-enzymatic synthesis of selectively ¹³C/¹⁵N-labeled RNA for NMR structural and dynamics studies

RNAs are an important class of cellular regulatory elements, and they are well characterized by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy in their folded or bound states. However, the apo or unfolded states are more difficult to characterize by either method. Particularly, effective NMR spectroscopy studies of RNAs in the past were hampered by chemical shift overlap of resonances and associated rapid signal loss due to line broadening for RNAs larger than the median size found in the PDB (~25 nt); most functional riboswitches are bigger than this median size. Incorporation of selective site-specific (13)C/(15)N-labeled nucleotides into RNAs promises to overcome this NMR size limitation. Unlike previous isotopic enrichment methods such as phosphoramidite, de novo, uniform-labeling, and selective-biomass approaches, this newer chemical-enzymatic selective method presents a number of advantages for producing labeled nucleotides over these other methods. For example, total chemical synthesis of nucleotides, followed by solid-phase synthesis of RNA using phosphoramidite chemistry, while versatile in incorporating isotope labels into RNA at any desired position, faces problems of low yields (<10%) that drop precipitously for oligonucleotides larger than 50 nt. The alternative method of de novo pyrimidine biosynthesis of NTPs is also a robust technique, with modest yields of up to 45%, but it comes at the cost of using 16 enzymes, expensive substrates, and difficulty in making some needed labeling patterns such as selective labeling of the ribose C1' and C5' and the pyrimidine nucleobase C2, C4, C5, or C6. Biomass-produced, uniformly or selectively labeled NTPs offer a third method, but suffer from low overall yield per labeled input metabolite and isotopic scrambling with only modest suppression of (13)C-(13)C couplings. In contrast to these four methods, our current chemo-enzymatic approach overcomes most of these shortcomings and allows for the synthesis of gram quantities of nucleotides with >80% yields while using a limited number of enzymes, six at most. The unavailability of selectively labeled ribose and base precursors had prevented the effective use of this versatile method until now. Recently, we combined an improved organic synthetic approach that selectively places (13)C/(15)N labels in the pyrimidine nucleobase (either (15)N1, (15)N3, (13)C2, (13)C4, (13)C5, or (13)C6 or any combination) with a very efficient enzymatic method to couple ribose with uracil to produce previously unattainable labeling patterns (Alvarado et al., 2014). Herein we provide detailed steps of both our chemo-enzymatic synthesis of custom nucleotides and their incorporation into RNAs with sizes ranging from 29 to 155 nt and showcase the dramatic improvement in spectral quality of reduced crowding and narrow linewidths. Applications of this selective labeling technology should prove valuable in overcoming two major obstacles, chemical shift overlap of resonances and associated rapid signal loss due to line broadening, that have impeded studying the structure and dynamics of large RNAs such as full-length riboswitches larger than the ~25 nt median size of RNA NMR structures found in the PDB.

Regio-selective chemical-enzymatic synthesis of pyrimidine nucleotides facilitates RNA structure and dynamics studies

Isotope labeling has revolutionized NMR studies of small nucleic acids, but to extend this technology to larger RNAs, site-specific labeling tools to expedite NMR structural and dynamics studies are required. Using enzymes from the pentose phosphate pathway, we coupled chemically synthesized uracil nucleobase with specifically (13) C-labeled ribose to synthesize both UTP and CTP in nearly quantitative yields. This chemoenzymatic method affords a cost-effective preparation of labels that are unattainable by current methods. The methodology generates versatile (13) C and (15) N labeling patterns which, when employed with relaxation-optimized NMR spectroscopy, effectively mitigate problems of rapid relaxation that result in low resolution and sensitivity. The methodology is demonstrated with RNAs of various sizes, complexity, and function: the exon splicing silencer 3 (27 nt), iron responsive element (29 nt), Pro-tRNA (76 nt), and HIV-1 core encapsidation signal (155 nt).

Magnetic resonance access to transiently formed protein complexes

Protein-protein interactions are of utmost importance to an understanding of biological phenomena since non-covalent and therefore reversible couplings between basic proteins leads to the formation of complex regulatory and adaptive molecular systems. Such systems are capable of maintaining their integrity and respond to external stimuli, processes intimately related to living organisms. These interactions, however, span a wide range of dissociation constants, from sub-nanomolar affinities in tight complexes to high-micromolar or even millimolar affinities in weak, transiently formed protein complexes. Herein, we demonstrate how novel NMR and EPR techniques can be used for the characterization of weak protein-protein (ligand) complexes. Applications to intrinsically disordered proteins and transiently formed protein complexes illustrate the potential of these novel techniques to study hitherto unobserved (and unobservable) higher-order structures of proteins.