Direct ¹³C-detected NMR experiments for mapping and characterization of hydrogen bonds in RNA

Relaxation Optimized Heteronuclear Experiments for Extending the Size Limit of RNA Nuclear Magnetic Resonance

The application of NMR to large RNAs has been limited by the inability to perform heteronuclear correlation experiments essential for resolving overlapping 1H NMR signals, determining interproton distance restraints and interhelical orientations for structure calculations, and evaluating conformational dynamics. Approaches exploiting 1H-13C correlations that are routinely applied to proteins and small RNAs of ∼60 nucleotides or fewer are impractical for larger RNAs due to rapid dipolar relaxation of protons by their attached carbons. Here we report a 2H-enhanced, 1H-15N correlation approach that enables atom-specific NMR characterization of much larger RNAs. Purine H8 transverse relaxation rates are reduced ∼20-fold with ribose perdeuteration, enabling efficient magnetization transfer via two-bond 1H-15N couplings. We focus on H8-N9 correlation spectra which benefit from favorable N9 chemical shift anisotropy. Chemical shift assignment is enabled by retention of protons at the C1' position, which allow measurement of two-bond H1'-N9 and through-space H1'-H8 correlations with only a minor effect on H8 relaxation. The approach is demonstrated for the 232 nucleotide HIV-1 Rev response element, where chemical shift assignments, 15N-edited nuclear Overhauser effects, and 1H-15N residual dipolar couplings are readily obtained from sensitive, high-resolution spectra. Heteronuclear correlated NMR methods that have been essential for the study of proteins can now be extended to RNAs of at least 78 kDa.

Dissecting the Conformational Heterogeneity of Stem-Loop Substructures of the Fifth Element in the 5'-Untranslated Region of SARS-CoV-2

Throughout the family of coronaviruses, structured RNA elements within the 5' region of the genome are highly conserved. The fifth stem-loop element from SARS-CoV-2 (5_SL5) represents an example of an RNA structural element, repeatedly occurring in coronaviruses. It contains a conserved, repetitive fold within its substructures SL5a and SL5b. We herein report the detailed characterization of the structure and dynamics of elements SL5a and SL5b that are located immediately upstream of the SARS-CoV-2 ORF1a/b start codon. Exploiting the unique ability of solution NMR methods, we show that the structures of both apical loops are modulated by structural differences in the remote parts located in their stem regions. We further integrated our high-resolution models of SL5a/b into the context of full-length 5_SL5 structures by combining different structural biology methods. Finally, we evaluated the impact of the two most common VoC mutations within 5_SL5 with respect to individual base-pair stability.

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.

Shine-Dalgarno Accessibility Governs Ribosome Binding to the Adenine Riboswitch

Translational riboswitches located in the 5' UTR of the messenger RNA (mRNA) regulate translation through variation of the accessibility of the ribosome binding site (RBS). These are the result of conformational changes in the riboswitch RNA governed by ligand binding. Here, we use a combination of single-molecule colocalization techniques (Single-Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) and Single-Molecule Kinetic Analysis of Ribosome Binding (SiM-KARB)) and microscale thermophoresis (MST) to investigate the adenine-sensing riboswitch in Vibrio vulnificus, focusing on the changes of accessibility between the ligand-free and ligand-bound states. We show that both methods faithfully report on the accessibility of the RBS within the riboswitch and that both methods identify an increase in accessibility upon adenine binding. Expanding on the regulatory context, we show the impact of the ribosomal protein S1 on the unwinding of the RNA secondary structure, thereby favoring ribosome binding even for the apo state. The determined rate constants suggest that binding of the ribosome is faster than the time required to change from the ON state to the OFF state, a prerequisite for efficient regulation decision.

Binding of 30S Ribosome Induces Single-stranded Conformation Within and Downstream of the Expression Platform in a Translational Riboswitch

Translational riboswitches are bacterial gene regulatory elements found in the 5'-untranslated region of mRNAs. They operate through a conformational refolding reaction that is triggered by a concentration change of a modulating small molecular ligand. The translation initiation region (TIR) is either released from or incorporated into base pairing interactions through the conformational switch. Hence, initiation of translation is regulated by the accessibility of the Shine-Dalgarno sequence and start codon. Interaction with the 30S ribosome is indispensable for the structural switch between functional OFF and ON states. However, on a molecular level it is still not fully resolved how the ribosome is accommodated near or at the translation initiation region in the context of translational riboswitches. The standby model of translation initiation postulates a binding site where the mRNA enters the ribosome and where it resides until the initiation site becomes unstructured and accessible. We here investigated the adenine-sensing riboswitch from Vibrio vulnificus. By application of a ¹⁹F labelling strategy for NMR spectroscopy that utilizes ligation techniques to synthesize differentially ¹⁹F labelled riboswitch molecules we show that nucleotides directly downstream of the riboswitch domain are first involved in productive interaction with the 30S ribosomal subunit. Upon the concerted action of ligand and the ribosomal protein rS1 the TIR becomes available and subsequently the 30S ribosome can slide towards the TIR. It will be interesting to see whether this is a general feature in translational riboswitches or if riboswitches exist where this region is structured and represent yet another layer of regulation.

13C Direct Detected NMR for Challenging Systems

Thanks to recent improvements in NMR spectrometer hardware and pulse sequence design, modern 13C NMR has become a useful tool for biomolecular applications. The complete assignment of a protein can be accomplished by using 13C detected multinuclear experiments and it can provide unique information relevant for the study of a variety of different biomolecules including paramagnetic proteins and intrinsically disordered proteins. A wide range of NMR observables can be measured, concurring to the structural and dynamic characterization of a protein in isolation, as part of a larger complex, or even inside a living cell. We present the different properties of 13C with respect to 1H, which provide the rationale for the experiments developed and their application, the technical aspects that need to be faced, and the many experimental variants designed to address different cases. Application areas where these experiments successfully complement proton NMR are also described.

NMR Spectroscopy of Large Functional RNAs: From Sample Preparation to Low-Gamma Detection

NMR spectroscopy is a potent method for the structural and biophysical characterization of RNAs. The application of NMR spectroscopy is restricted in RNA size and most often requires isotope-labeled or even selectively labeled RNAs. Additionally, new NMR pulse sequences, such as the heteronuclear-detected NMR experiments, are introduced. We herein provide detailed protocols for the preparation of isotope-labeled RNA for NMR spectroscopy via in vitro transcription. This protocol covers all steps, from the preparation of DNA template to the transcription of milligram RNA quantities. Moreover, we present a protocol for a chemo-enzymatic approach to introduce a single modified nucleotide at any position of any RNA. Regarding NMR methodology, we share protocols for the implementation of a suite of heteronuclear-detected NMR experiments including ¹³ C-detected experiments for ribose assignment and amino groups, the CN-spin filter heteronuclear single quantum coherence (HSQC) for imino groups and the ¹⁵ N-detected band-selective excitation short transient transverse-relaxation-optimized spectroscopy (BEST-TROSY) experiment. © 2020 The Authors. Basic Protocol 1: Preparation of isotope-labeled RNA samples with in vitro transcription using T7 RNAP, DEAE chromatography, and RP-HPLC purification Alternate Protocol 1: Purification of isotope-labeled RNA from in vitro transcription with preparative PAGE Alternate Protocol 2: Purification of isotope-labeled RNA samples from in vitro transcription via centrifugal concentration Support Protocol 1: Preparation of DNA template from plasmid Support Protocol 2: Preparation of PCR DNA as template Support Protocol 3: Preparation of T7 RNA Polymerase (T7 RNAP) Support Protocol 4: Preparation of yeast inorganic pyrophosphatase (YIPP) Basic Protocol 2: Preparation of site-specific labeled RNAs using a chemo-enzymatic synthesis Support Protocol 5: Synthesis of modified nucleoside 3',5'-bisphosphates Support Protocol 6: Preparation of T4 RNA Ligase 2 Support Protocol 7: Setup of NMR spectrometer for heteronuclear-detected NMR experiments Support Protocol 8: IPAP and DIPAP for homonuclear decoupling Basic Protocol 3: ¹³ C-detected 3D (H)CC-TOCSY, (H)CPC, and (H)CPC-CCH-TOCSY experiments for ribose assignment Basic Protocol 4: ¹³ C-detected 2D CN-spin filter HSQC experiment Basic Protocol 5: ¹³ C-detected C(N)H-HDQC experiment for the detection of amino groups Support Protocol 9: ¹³ C-detected CN-HSQC experiment for amino groups Basic Protocol 6: ¹³ C-detected "amino"-NOESY experiment Basic Protocol 7: ¹⁵ N-detected BEST-TROSY experiment.

More than Proton Detection-New Avenues for NMR Spectroscopy of RNA

Ribonucleic acid oligonucleotides (RNAs) play pivotal roles in cellular function (riboswitches), chemical biology applications (SELEX-derived aptamers), cell biology and biomedical applications (transcriptomics). Furthermore, a growing number of RNA forms (long non-coding RNAs, circular RNAs) but also RNA modifications are identified, showing the ever increasing functional diversity of RNAs. To describe and understand this functional diversity, structural studies of RNA are increasingly important. However, they are often more challenging than protein structural studies as RNAs are substantially more dynamic and their function is often linked to their structural transitions between alternative conformations. NMR is a prime technique to characterize these structural dynamics with atomic resolution. To extend the NMR size limitation and to characterize large RNAs and their complexes above 200 nucleotides, new NMR techniques have been developed. This Minireview reports on the development of NMR methods that utilize detection on low-γ nuclei (heteronuclei like 13 C or 15 N with lower gyromagnetic ratio than 1 H) to obtain unique structural and dynamic information for large RNA molecules in solution. Experiments involve through-bond correlations of nucleobases and the phosphodiester backbone of RNA for chemical shift assignment and make information on hydrogen bonding uniquely accessible. Previously unobservable NMR resonances of amino groups in RNA nucleobases are now detected in experiments involving conformational exchange-resistant double-quantum 1 H coherences, detected by 13 C NMR spectroscopy. Furthermore, 13 C and 15 N chemical shifts provide valuable information on conformations. All the covered aspects point to the advantages of low-γ nuclei detection experiments in RNA.

Novel 13 C-detected NMR Experiments for the Precise Detection of RNA Structure

Up to now, NMR spectroscopic investigations of RNA have utilized imino proton resonances as reporters for base pairing and RNA structure. The nucleobase amino groups are often neglected, since most of their resonances are broadened beyond detection due to rotational motion around the C-NH₂ bond. Here, we present ¹³ C-detected NMR experiments for the characterization of all RNA amino groups irrespective of their motional behavior. We have developed a C(N)H-HDQC experiment that enables the observation of a complete set of sharp amino resonances through the detection of proton-NH₂ double quantum coherences. Further, we present an "amino"-NOESY experiment to detect NOEs to amino protons, which are undetectable by any other conventional NOESY experiment. Together, these experiments allow the exploration of additional chemical shift information and inter-residual proton distances important for high-resolution RNA secondary and tertiary structure determination.

Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme

In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Ligand-modulated folding of the full-length adenine riboswitch probed by NMR and single-molecule FRET spectroscopy

The full-length translation-regulating add adenine riboswitch (Asw) from Vibrio vulnificus has a more complex conformational space than its isolated aptamer domain. In addition to the predicted apo (apoA) and holo conformation that feature the conserved three-way junctional purine riboswitch aptamer, it adopts a second apo (apoB) conformation with a fundamentally different secondary structure. Here, we characterized the ligand-dependent conformational dynamics of the full-length add Asw by NMR and by single-molecule FRET (smFRET) spectroscopy. Both methods revealed an adenine-induced secondary structure switch from the apoB-form to the apoA-form that involves no tertiary structural interactions between aptamer and expression platform. This strongly suggests that the add Asw triggers translation by capturing the apoA-form secondary structure in the holo state. Intriguingly, NMR indicated a homogenous, docked aptamer kissing loop fold for apoA and holo, while smFRET showed persistent aptamer kissing loop docking dynamics between comparably stable, undocked and docked substates of the apoA and the holo conformation. Unraveling the folding of large junctional riboswitches thus requires the integration of complementary solution structural techniques such as NMR and smFRET.

Evaluation of 15N-detected H-N correlation experiments on increasingly large RNAs

Recently, ¹⁵N-detected multidimensional NMR experiments have been introduced for the investigation of proteins. Utilization of the slow transverse relaxation of nitrogen nuclei in a ¹⁵N-TROSY experiment allowed recording of high quality spectra for high molecular weight proteins, even in the absence of deuteration. Here, we demonstrate the applicability of three ¹⁵N-detected H-N correlation experiments (TROSY, BEST-TROSY and HSQC) to RNA. With the newly established ¹⁵N-detected BEST-TROSY experiment, which proves to be the most sensitive ¹⁵N-detected H-N correlation experiment, spectra for five RNA molecules ranging in size from 5 to 100 kDa were recorded. These spectra yielded high resolution in the ¹⁵N-dimension even for larger RNAs since the increase in line width with molecular weight is more pronounced in the ¹H- than in the ¹⁵N-dimension. Further, we could experimentally validate the difference in relaxation behavior of imino groups in AU and GC base pairs. Additionally, we showed that ¹⁵N-detected experiments theoretically should benefit from sensitivity and resolution advantages at higher static fields but that the latter is obscured by exchange dynamics within the RNAs.