Rapid NMR screening of RNA secondary structure and binding

RNA folding kinetics control riboswitch sensitivity in vivo

Riboswitches are ligand-responsive gene-regulatory RNA elements that perform key roles in maintaining cellular homeostasis. Understanding how riboswitch sensitivity is controlled is critical to understanding how highly conserved aptamer domains are deployed in a variety of contexts with different sensitivity demands. Here we uncover new roles by which RNA folding dynamics control riboswitch sensitivity in cells. By investigating the Clostridium beijerinckii pfl ZTP riboswitch, we identify multiple mechanistic routes of altering expression platform sequence and structure to slow RNA folding, all of which enhance riboswitch sensitivity. Applying these methods to riboswitches with diverse aptamer architectures that regulate transcription and translation with ON and OFF logic demonstrates the generality of our findings, indicating that any riboswitch that operates in a kinetic regime can be sensitized by slowing expression platform folding. Comparison of the most sensitized versions of these switches to equilibrium aptamer:ligand dissociation constants suggests a limit to the sensitivities achievable by kinetic RNA switches. Our results add to the growing suite of knowledge and approaches that can be used to rationally program cotranscriptional RNA folding for biotechnology applications, and suggest general RNA folding principles for understanding dynamic RNA systems in other areas of biology.

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.

TransUFold: Unlocking the structural complexity of short and long RNA with pseudoknots

The RNA secondary structure is like a blueprint that holds the key to unlocking the mysteries of RNA function and 3D structure. It serves as a crucial foundation for investigating the complex world of RNA, making it an indispensable component of research in this exciting field. However, pseudoknots cannot be accurately predicted by conventional prediction methods based on free energy minimization, which results in a performance bottleneck. To this end, we propose a deep learning-based method called TransUFold to train directly on RNA data annotated with structure information. It employs an encoder-decoder network architecture, named Vision Transformer, to extract long-range interactions in RNA sequences and utilizes convolutions with lateral connections to supplement short-range interactions. Then, a post-processing program is designed to constrain the model's output to produce realistic and effective RNA secondary structures, including pseudoknots. After training TransUFold on benchmark datasets, we outperform other methods in test data on the same family. Additionally, we achieve better results on longer sequences up to 1600 nt, demonstrating the outstanding performance of Vision Transformer in extracting long-range interactions in RNA sequences. Finally, our analysis indicates that TransUFold produces effective pseudoknot structures in long sequences. As more high-quality RNA structures become available, deep learning-based prediction methods like Vision Transformer can exhibit better performance.

Modeling the 3D structure and conformational dynamics of very large RNAs using coarse-grained molecular simulations

We describe a computational approach to building and simulating realistic 3D models of very large RNA molecules (>1000 nucleotides) at a resolution of one "bead" per nucleotide. The method starts with a predicted secondary structure and uses several stages of energy minimization and Brownian dynamics (BD) simulation to build 3D models. A key step in the protocol is the temporary addition of a 4 ^th spatial dimension that allows all predicted helical elements to become disentangled from each other in an effectively automated way. We then use the resulting 3D models as input to Brownian dynamics simulations that include hydrodynamic interactions (HIs) that allow the diffusive properties of the RNA to be modelled as well as enabling its conformational dynamics to be simulated. To validate the dynamics part of the method, we first show that when applied to small RNAs with known 3D structures the BD-HI simulation models accurately reproduce their experimental hydrodynamic radii (Rh). We then apply the modelling and simulation protocol to a variety of RNAs for which experimental Rh values have been reported ranging in size from 85 to 3569 nucleotides. We show that the 3D models, when used in BD-HI simulations, produce hydrodynamic radii that are usually in good agreement with experimental estimates for RNAs that do not contain tertiary contacts that persist even under very low salt conditions. Finally, we show that sampling of the conformational dynamics of large RNAs on timescales of 100 µs is computationally feasible with BD-HI simulations.

Site-Specific Labeling of RNAs with Modified and 19 F-Labeled Nucleotides by Chemo-Enzymatic Synthesis

More than 170 post-transcriptional modifications of RNAs have currently been identified. Detailed biophysical investigations of these modifications have been limited since large RNAs containing these post-transcriptional modifications are difficult to produce. Further, adequate readout of spectroscopic fingerprints are important, necessitating additional labeling procedures beyond the naturally occurring RNA modifications. Here, we report the chemo-enzymatic synthesis of RNA modifications and several structurally similar fluorine-modified analogs further optimizing a recently developed methodology.^[1] This chemo-enzymatic method allows synthesis of also large RNAs. We were able to incorporate 16 modified nucleotides and 6 ¹⁹ F-labeled nucleotides. To showcase the applicability of such modified large RNAs, we incorporated a ¹⁹ F-labeled cytidine into the aptamer domain of the 2'dG sensing riboswitch (2'dG-sw) from Mesoplasma florum, enabling characterizing RNA fold, ligand binding and kinetics. Thanks to the large chemical shift dispersion of ¹⁹ F, we can detect conformational heterogeneity in the apo state of the riboswitch.

Characterizing Fluorescence Properties of Turn-on RNA Aptamers

Fluorescent RNA aptamers are tools for studying RNA localization and interactions in vivo. The photophysical properties of these in vitro selected RNAs should be characterized prior to cellular imaging experiments. Here, we describe the process of determining the fluorophore affinity, fluorescence enhancement, and fluorescence lifetime(s) of the Mango-III fluorescence turn-on aptamer. Parameters determined through these protocols will aid in establishing conditions for live-cell imaging.

Modified Nucleotides for Chemical and Enzymatic Synthesis of Therapeutic RNA

In recent years, RNA has emerged as a medium with a broad spectrum of therapeutic potential, however, for years, a group of short RNA fragments was studied and considered therapeutic molecules. In nature, RNA plays both functions, with coding and non-coding potential. For RNA, like any other therapeutic, to be used clinically, certain barriers must be crossed. Among them, there are biocompatibility, relatively low toxicity, bioavailability, increased stability, target efficiency and low off-target effects. In the case of RNA, most of these obstacles can be overcome by incorporating modified nucleotides into its structure. This may be achieved by both, in vitro and in vivo biosynthetic methods, as well as chemical synthesis. Some advantages and disadvantages of each approach are summarized here. The wide range of nucleotide analogues has been tested for their utility as monomers for RNA synthesis. Many of them have been successfully implemented, and a lot of pre-clinical and clinical studies involving modified RNA have been carried out. Some of these medications have already been introduced into clinics. After the huge success of RNA-based vaccines that were introduced into widespread use in 2020, and the introduction to the market of some RNA-based drugs, RNA therapeutics containing modified nucleotides appear to be the future of medicine.

The cotranscriptional folding landscape for two cyclic di-nucleotide-sensing riboswitches with highly homologous aptamer domains acting either as ON- or OFF-switches

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 μM to 0.25 μM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

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.

Mapping the conformational landscape of the neutral network of RNA sequences that connect two functional distinctly different ribozymes

During evolution of an RNA world, the development of enzymatic function was essential. Such enzymatic function was linked to RNA sequences capable of adopting specific RNA folds that possess catalytic pockets to promote catalysis. Within this primordial RNA world, initially evolved self-replicating ribozymes presumably mutated to ribozymes with new functions. In 2000, Schultes and Bartel investigated such conversion from one ribozyme to a new ribozyme with distinctly different catalytic functions. Within a neutral network that linked these two prototype ribozymes, a single RNA chain could be identified that exhibited both enzymatic functions. This system serves as a paradigm for an evolutionary system that allows neutral drifts by stepwise mutation from one ribozyme into a different ribozyme without loss of intermittent function. Here, we investigated this complex functional diversification by analyzing several RNA sequences within this neutral network between two ribozymes with class III ligase activity and with self-cleavage reactivity. We utilized rapid RNA sample preparation for NMR spectroscopic studies together with SHAPE analysis and in-line probing to characterize secondary structure changes within the neutral network. Our investigations allowed delineation of the 2 nd structure space and by comparison with the previously determined catalytic function allowed correlation of the structure-function relation of ribozyme function in this neutral network.

Optimization of RNA Aptamer SELEX Methods: Improved Aptamer Transcript 3'-End Homogeneity, PAGE Purification Yield, and Target-Bound Aptamer RNA Recovery

Since its inception in the early 1990s, SELEX remains the gold standard for discovering RNA aptamers specific for proteins and small molecules. The SELEX process has undergone countless modifications and now encompasses a breadth of innovative selection schemes to pare an aptamer library toward target-specific aptamers. Common to all these RNA aptamer SELEX processes are the steps for the preparation of DNA template and in vitro transcription of aptamer RNA. These steps have remained mostly unchanged over the past three decades and would benefit from optimization. We focused on three key areas: improving the homogeneity of in vitro transcribed aptamer RNA, increasing the efficiency of in vitro transcribed aptamer RNA purification by PAGE, and improving the quality of target-bound aptamer RNA recovered during SELEX. Together, these optimizations contribute toward a more efficient SELEX process and are applicable to both protein-based and cell-based RNA aptamer selections.

Overview of Methods for Large-Scale RNA Synthesis

In recent years, it has become clear that RNA molecules are involved in almost all vital cellular processes and pathogenesis of human disorders. The functional diversity of RNA comes from its structural richness. Although composed of only four nucleotides, RNA molecules present a plethora of secondary and tertiary structures critical for intra and intermolecular contacts with other RNAs and ligands (proteins, small metabolites, etc.). In order to fully understand RNA function it is necessary to define its spatial structure. Crystallography, nuclear magnetic resonance and cryogenic electron microscopy have demonstrated considerable success in determining the structures of biologically important RNA molecules. However, these powerful methods require large amounts of sample. Despite their limitations, chemical synthesis and in vitro transcription are usually employed to obtain milligram quantities of RNA for structural studies, delivering simple and effective methods for large-scale production of homogenous samples. The aim of this paper is to provide an overview of methods for large-scale RNA synthesis with emphasis on chemical synthesis and in vitro transcription. We also present our own results of testing the efficiency of these approaches in order to adapt the material acquisition strategy depending on the desired RNA construct.

Characterization of Structure and Dynamics of the Guanidine-II Riboswitch from Escherichia coli by NMR Spectroscopy and Small-Angle X-ray Scattering (SAXS)

Riboswitches are regulatory RNA elements that undergo functionally important allosteric conformational switching upon binding of specific ligands. The here investigated guanidine‐II riboswitch binds the small cation, guanidinium, and forms a kissing loop‐loop interaction between its P1 and P2 hairpins. We investigated the structural changes to support previous studies regarding the binding mechanism. Using NMR spectroscopy, we confirmed the structure as observed in crystal structures and we characterized the kissing loop interaction upon addition of Mg²⁺ and ligand for the riboswitch aptamer from Escherichia coli. We further investigated closely related mutant constructs providing further insight into functional differences between the two (different) hairpins P1 and P2. Formation of intermolecular interactions were probed by small‐angle X‐ray scattering (SAXS) and NMR DOSY data. All data are consistent and show the formation of oligomeric states of the riboswitch induced by Mg²⁺ and ligand binding.

NMR chemical shift assignments of RNA oligonucleotides to expand the RNA chemical shift database

RNAs play myriad functional and regulatory roles in the cell. Despite their significance, three-dimensional structure elucidation of RNA molecules lags significantly behind that of proteins. NMR-based studies are often rate-limited by the assignment of chemical shifts. Automation of the chemical shift assignment process can greatly facilitate structural studies, however, accurate chemical shift predictions rely on a robust and complete chemical shift database for training. We searched the Biological Magnetic Resonance Data Bank (BMRB) to identify sequences that had no (or limited) chemical shift information. Here, we report the chemical shift assignments for 12 RNA hairpins designed specifically to help populate the BMRB.

Production of Structured RNA Fragments by In Vitro Transcription and HPLC Purification

The understanding of the functional importance of RNA has increased enormously in the last decades. This has required research on the RNA molecules themselves, with the concomitant need for obtaining purified RNA samples, such as for structural studies by NMR or other methods. The main method to create labeled and unlabeled RNA, T7 in vitro transcription, suffers from sequence-dependent yield and often low homogeneity for short constructs (<100 nt) and requires laborious purification. Additionally, the design of structured RNA fragments mimicking the structure of a larger biological RNA is often not straightforward. Secondary structure simulations can be used to make reliable predictions about the folding of a particular RNA fragment. In this article, we describe how to design an RNA construct of interest from a larger sequence, and we combine several previously published improvements of the in vitro transcription method, such as the use of 2'-methoxy modifications and dimethyl sulfoxide or the use of tandem repeats, to increase yield and purity of in vitro-transcribed RNA. Together with a high-performance liquid chromatography (HPLC) purification procedure using both reversed-phase ion-pairing and ion-exchange HPLC, we provide a robust protocol to obtain highly pure RNA of short to intermediate length in large quantities. The protocol optimizes yield, especially for RNA starting with nucleotides other than G. At the same time, it is simplified, and the required time is reduced. The protocols described here constitute a versatile pipeline for the production of purified RNA samples and are suitable for users with little experience in liquid chromatography. © 2021 The Authors. Basic Protocol 1: RNA construct design Basic Protocol 2: DNA template production and in vitro transcription Alternate Protocol: Tandem transcription and RNase H cleavage Basic Protocol 3: Reversed-phase ion-pairing HPLC purification Basic Protocol 4: Ion-exchange HPLC purification.

New Frontiers for Site-Directed RNA Editing: Harnessing Endogenous ADARs

RNA editing activity can be exploited for the restoration of disease-causing nonsense and missense mutations and as a tool to manipulate the transcriptome in a simple and programmable way. The general concept is called site-directed RNA editing and has high potential for translation into the clinics. Due to its different mode of action RNA editing may well complement gene editing and other gene therapy options. In this method chapter, we particularly highlight RNA editing strategies that harness endogenous ADARs. Such strategies circumvent the delivery and expression of engineered editases and are notably precise and simple. This is particularly true if endogenous ADARs are recruited with chemically modified antisense oligonucleotides, an approach we call RESTORE (recruiting endogenous ADAR to specific transcripts for oligonucleotide-mediated RNA editing). To foster the research and development of RESTORE we now report a detailed protocol for the procedure of editing reactions, and a protocol for the generation of partly chemically modified RESTORE ASOs with a combination of in vitro transcription and ligation.

A Tale of Two Transitions: The Unfolding Mechanism of the prfA RNA Thermosensor

RNA thermosensors (RNATs), found in the 5' untranslated region (UTR) of some bacterial messenger RNAs (mRNAs), control the translation of the downstream gene in a temperature-dependent manner. In Listeria monocytogenes, the expression of a key transcription factor, PrfA, is mediated by an RNAT in its 5' UTR. PrfA functions as a master regulator of virulence in L. monocytogenes, controlling the expression of many virulence factors. The temperature-regulated expression of PrfA by its RNAT element serves as a signal of successful host invasion for the bacteria. Structurally, the prfA RNAT bears little resemblance to known families of RNATs, and prior studies demonstrated that the prfA RNAT is highly responsive over a narrow temperature range. Herein, we have undertaken a comprehensive mutational and thermodynamic analysis to ascertain the molecular determinants of temperature sensitivity. We provide evidence to support the idea that the prfA RNAT unfolding is different from that of cssA, a well-characterized RNAT, suggesting that these RNATs function via distinct mechanisms. Our data show that the unfolding of the prfA RNAT occurs in two distinct events and that the internal loops play an important role in mediating the cooperativity of RNAT unfolding. We further demonstrated that regions distal to the ribosome binding site (RBS) not only contribute to RNAT structural stability but also impact translation of the downstream message. Our collective results provide insight connecting the thermal stability of the prfA RNAT structure, unfolding energetics, and translational control.

19 F NMR-Based Fragment Screening for 14 Different Biologically Active RNAs and 10 DNA and Protein Counter-Screens

We report here the nuclear magnetic resonance ¹⁹F screening of 14 RNA targets with different secondary and tertiary structure to systematically assess the druggability of RNAs. Our RNA targets include representative bacterial riboswitches that naturally bind with nanomolar affinity and high specificity to cellular metabolites of low molecular weight. Based on counter‐screens against five DNAs and five proteins, we can show that RNA can be specifically targeted. To demonstrate the quality of the initial fragment library that has been designed for easy follow‐up chemistry, we further show how to increase binding affinity from an initial fragment hit by chemistry that links the identified fragment to the intercalator acridine. Thus, we achieve low‐micromolar binding affinity without losing binding specificity between two different terminator structures.

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.

The conformational landscape of transcription intermediates involved in the regulation of the ZMP-sensing riboswitch from Thermosinus carboxydivorans

Recently, prokaryotic riboswitches have been identified that regulate transcription in response to change of the concentration of secondary messengers. The ZMP (5-Aminoimidazole-4-carboxamide ribonucleotide (AICAR))-sensing riboswitch from Thermosinus carboxydivorans is a transcriptional ON-switch that is involved in purine and carbon-1 metabolic cycles. Its aptamer domain includes the pfl motif, which features a pseudoknot, impeding rho-independent terminator formation upon stabilization by ZMP interaction. We herein investigate the conformational landscape of transcriptional intermediates including the expression platform of this riboswitch and characterize the formation and unfolding of the important pseudoknot structure in the context of increasing length of RNA transcripts. NMR spectroscopic data show that even surprisingly short pre-terminator stems are able to disrupt ligand binding and thus metabolite sensing. We further show that the pseudoknot structure, a prerequisite for ligand binding, is preformed in transcription intermediates up to a certain length. Our results describe the conformational changes of 13 transcription intermediates of increasing length to delineate the change in structure as mRNA is elongated during transcription. We thus determine the length of the key transcription intermediate to which addition of a single nucleotide leads to a drastic drop in ZMP affinity.

A robust and versatile method for production and purification of large-scale RNA samples for structural biology

Recent findings in genome-wide transcriptomics revealed that RNAs are involved in almost every biological process, across all domains of life. The characterization of native RNAs of unknown function and structure is particularly challenging due to their typical low abundance in the cell and the inherent sensitivity toward ubiquitous RNA degrading enzymes. Therefore, robust in vitro synthesis and extensive work-up methods are often needed to obtain samples amenable for biochemical, biophysical, and structural studies. Here, we present a protocol that combines the most recent advances in T7 in vitro transcription methodology with reverse phase ion pairing and ion exchange HPLC purification of RNAs for the production of yield-optimized large-scale samples. The method is easy to follow, robust and suitable for users with little or no experience within the field of biochemistry or chromatography. The complete execution of this method, for example, for production of isotopically labeled NMR samples, can be performed in less than a week.

Base-pair conformational switch modulates miR-34a targeting of Sirt1 mRNA

MicroRNAs (miRNAs) regulate the levels of translation of messenger RNAs (mRNAs). At present, the major parameter that can explain the selection of the target mRNA and the efficiency of translation repression is the base pairing between the 'seed' region of the miRNA and its counterpart mRNA¹. Here we use R_1ρ relaxation-dispersion nuclear magnetic resonance² and molecular simulations³ to reveal a dynamic switch-based on the rearrangement of a single base pair in the miRNA-mRNA duplex-that elongates a weak five-base-pair seed to a complete seven-base-pair seed. This switch also causes coaxial stacking of the seed and supplementary helix fitting into human Argonaute 2 protein (Ago2), reminiscent of an active state in prokaryotic Ago^4,5. Stabilizing this transient state leads to enhanced repression of the target mRNA in cells, revealing the importance of this miRNA-mRNA structure. Our observations tie together previous findings regarding the stepwise miRNA targeting process from an initial 'screening' state to an 'active' state, and unveil the role of the RNA duplex beyond the seed in Ago2.

Challenges to optimizing RNA nanostructures for large scale production and controlled therapeutic properties

Nucleic acids have been utilized to construct an expansive collection of nanoarchitectures varying in design, physicochemical properties, cellular processing and biomedical applications. However, the broader therapeutic adaptation of nucleic acid nanoassemblies in general, and RNA-based nanoparticles in particular, have faced several challenges in moving towards (pre)clinical settings. For one, the large-batch synthesis of nucleic acids is still under development, with multi-stranded and chemically modified assemblies requiring greater production capacity while maintaining consistent medical-grade outputs. Furthermore, the unknown immunostimulation by these nanomaterials poses additional challenges, necessary to be overcome for optimizing future development of clinically approved RNA nanoparticles.

Genetic Code Expansion Facilitates Position-Selective Labeling of RNA for Biophysical Studies

Nature relies on reading and synthesizing the genetic code with high fidelity. Nucleic acid building blocks that are orthogonal to the canonical A-T and G-C base-pairs are therefore uniquely suitable to facilitate position-specific labeling of nucleic acids. Here, we employ the orthogonal kappa-xanthosine-base-pair for in vitro transcription of labeled RNA. We devised an improved synthetic route to obtain the phosphoramidite of the deoxy-version of the kappa nucleoside in solid phase synthesis. From this DNA template, we demonstrate the reliable incorporation of xanthosine during in vitro transcription. Using NMR spectroscopy, we show that xanthosine introduces only minor structural changes in an RNA helix. We furthermore synthesized a clickable 7-deaza-xanthosine, which allows to site-specifically modify transcribed RNA molecules with fluorophores or other labels.

A Modified In Vitro Transcription Approach to Improve RNA Synthesis and Ribozyme Cleavage Efficiency

RNA elements such as catalytic RNA, riboswitch, microRNA, and long non-coding RNA perform a major role in cellular processes. A complete understanding of cellular processes is impossible without knowing the structure-function relationship of participating RNA molecules that ultimately requires large quantities of pure RNAs. Thus, structural/functional analyses of emerging RNAs necessitate revised protocols for improved RNA quantity and quality. Here we present a modified in vitro transcription protocol to enhance ribozyme cleaving efficiency and RNA yield by working on two variables, i.e., incubation temperature and limiting GTPs. Following an improved RNA synthesis, the target RNA is purified from transcription mixture components through denaturing size-exclusion chromatography. The protocol confirms that cyclic elevated incubation temperatures during transcription and increased concentrations of GTPs improve the production rate of RNA. Our modified in vitro transcription method improves the ribozyme cleaving efficiency and targets RNA yield by four- to fivefold that can benefit almost any RNA-related study from protein-RNA interaction analysis to crystallography.

Conformational switch in the ribosomal protein S1 guides unfolding of structured RNAs for translation initiation

Initiation of bacterial translation requires that the ribosome-binding site in mRNAs adopts single-stranded conformations. In Gram-negative bacteria the ribosomal protein S1 (rS1) is a key player in resolving of structured elements in mRNAs. However, the exact mechanism of how rS1 unfolds persistent secondary structures in the translation initiation region (TIR) is still unknown. Here, we show by NMR spectroscopy that Vibrio vulnificus rS1 displays a unique architecture of its mRNA-binding domains, where domains D3 and D4 provide the mRNA-binding platform and cover the nucleotide binding length of the full-length rS1. D5 significantly increases rS1’s chaperone activity, although it displays structural heterogeneity both in isolation and in presence of the other domains, albeit to varying degrees. The heterogeneity is induced by the switch between the two equilibrium conformations and is triggered by an order-to-order transition of two mutually exclusive secondary structures (β-strand-to-α-helix) of the ‘AERERI’ sequence. The conformational switching is exploited for melting of structured 5′-UTR’s, as the conformational heterogeneity of D5 can compensate the entropic penalty of complex formation. Our data thus provides a detailed understanding of the intricate coupling of protein and RNA folding dynamics enabling translation initiation of structured mRNAs.

Isothermal Titration Calorimetry Measurements of Riboswitch-Ligand Interactions

One of the many ways by which bacteria control gene expression is through cis-acting regulatory mRNA elements called riboswitches. By specifically binding to small molecules or metabolites and pairing the binding event to an RNA structure change, riboswitches link a metabolic input to a transcriptional or translational output. For over a decade, isothermal titration calorimetry (ITC) has been used to investigate how riboswitches interact with small molecules. We present methods for assaying RNA-ligand interactions using ITC and analyzing resulting data to estimate thermodynamic parameters associated with binding.

Chemo-Enzymatic Synthesis of Position-Specifically Modified RNA for Biophysical Studies including Light Control and NMR Spectroscopy

The investigation of non-coding RNAs requires RNAs containing modifications at every possible position within the oligonucleotide. Here, we present the chemo-enzymatic RNA synthesis containing photoactivatable or ¹³ C,¹⁵ N-labelled nucleosides. All four ribonucleotides containing ortho-nitrophenylethyl (NPE) photocages, photoswitchable azobenzene C-nucleotides and ¹³ C,¹⁵ N-labelled nucleotides were incorporated position-specifically in high yields. We applied this approach for the synthesis of light-inducible 2'dG-sensing riboswitch variants and detected ligand-induced structural reorganization upon irradiation by NMR spectroscopy. This chemo-enzymatic method opens the possibility to incorporate a wide range of modifications at any desired position of RNAs of any lengths beyond the limits of solid-phase synthesis.

Accurate Measurement of Residual Dipolar Couplings in Large RNAs by Variable Flip Angle NMR

NMR approaches using nucleotide-specific deuterium labeling schemes have enabled structural studies of biologically relevant RNAs of increasing size and complexity. Although local structure is well-determined using these methods, definition of global structural features, including relative orientations of independent helices, remains a challenge. Residual dipolar couplings, a potential source of orientation information, have not been obtainable for large RNAs due to poor sensitivity resulting from rapid heteronuclear signal decay. Here we report a novel multiple quantum NMR method for RDC determination that employs flip angle variation rather than a coupling evolution period. The accuracy of the method and its utility for establishing interhelical orientations are demonstrated for a 36-nucleotide RNA, for which comparative data could be obtained. Applied to a 78 kDa Rev response element from the HIV-1 virus, which has an effective rotational correlation time of ca. 160 ns, the method yields sensitivity gains of an order of magnitude or greater over existing approaches. Solution-state access to structural organization in RNAs of at least 230 nucleotides is now possible.

HIV-1 Matrix Protein Interactions with tRNA: Implications for Membrane Targeting

The N-terminally myristoylated matrix (MA) domain of the HIV-1 Gag polyprotein promotes virus assembly by targeting Gag to the inner leaflet of the plasma membrane. Recent studies indicate that, prior to membrane binding, MA associates with cytoplasmic tRNAs (including tRNA^Lys3), and in vitro studies of tRNA-dependent MA interactions with model membranes have led to proposals that competitive tRNA interactions contribute to membrane discrimination. We have characterized interactions between native, mutant, and unmyristylated (myr-) MA proteins and recombinant tRNA^Lys3 by NMR spectroscopy and isothermal titration calorimetry. NMR experiments confirm that tRNA^Lys3 interacts with a patch of basic residues that are also important for binding to the plasma membrane marker, phosphatidylinositol-4,5-bisphosphate [PI(4,5)P₂]. Unexpectedly, the affinity of MA for tRNA^Lys3 (K_d = 0.63 ± 0.03 μM) is approximately 1 order of magnitude greater than its affinity for PI(4,5)P₂-enriched liposomes (K_d(apparent) = 10.2 ± 2.1 μM), and NMR studies indicate that tRNA^Lys3 binding blocks MA association with liposomes, including those enriched with PI(4,5)P₂, phosphatidylserine, and cholesterol. However, the affinity of MA for tRNA^Lys3 is diminished by mutations or sample conditions that promote myristate exposure. Since Gag-Gag interactions are known to promote myristate exposure, our findings support virus assembly models in which membrane targeting and genome binding are mechanistically coupled.

A guide to large-scale RNA sample preparation

RNA is becoming more important as an increasing number of functions, both regulatory and enzymatic, are being discovered on a daily basis. As the RNA boom has just begun, most techniques are still in development and changes occur frequently. To understand RNA functions, revealing the structure of RNA is of utmost importance, which requires sample preparation. We review the latest methods to produce and purify a variation of RNA molecules for different purposes with the main focus on structural biology and biophysics. We present a guide aimed at identifying the most suitable method for your RNA and your biological question and highlighting the advantages of different methods. Graphical abstract In this review we present different methods for large-scale production and purification of RNAs for structural and biophysical studies.

Life times of metastable states guide regulatory signaling in transcriptional riboswitches

Transcriptional riboswitches modulate downstream gene expression by a tight coupling of ligand-dependent RNA folding kinetics with the rate of transcription. RNA folding pathways leading to functional ON and OFF regulation involve the formation of metastable states within well-defined sequence intervals during transcription. The kinetic requirements for the formation and preservation of these metastable states in the context of transcription remain unresolved. Here, we reversibly trap the previously defined regulatory relevant metastable intermediate of the Mesoplasma florum 2′-deoxyguanosine (2′dG)-sensing riboswitch using a photocaging-ligation approach, and monitor folding to its native state by real-time NMR in both presence and absence of ligand. We further determine transcription rates for two different bacterial RNA polymerases. Our results reveal that the riboswitch functions only at transcription rates typical for bacterial polymerases (10–50 nt s⁻¹) and that gene expression is modulated by 40–50% only, while subtle differences in folding rates guide population ratios within the structural ensemble to a specific regulatory outcome.

Riboswitches are RNA-based regulatory elements, which regulate downstream gene expression by binding of small molecular weight ligands. Here the authors demonstrate the molecular mechanism of a transcriptional riboswitch that integrates changes in transcription rates, metabolite concentration, and kinetic on- and off-rates of ligand binding.

A T7 RNA Polymerase Mutant Enhances the Yield of 5'-Thienoguanosine-Initiated RNAs

Spectroscopic methods, which are used to establish RNA structure-function relationships, require strategies for post-synthetic, site-specific incorporation of chemical probes into target RNAs. For RNAs larger than 50 nt, the enzymatic incorporation of a nucleoside or nucleotide monophosphate guanosine analogue (G analogue) at their 5'-end is routinely achieved by T7 RNA polymerase (T7RNAP)-mediated in vitro transcription (IVT) of the appropriate DNA template containing a GTP-initiating class III Φ6.5 promoter. However, when high G analogue:GTP ratios are used to bias G analogue incorporation at the 5'-end, RNA yield is compromised. Here, we show that the use of a T7RNAP P266L mutant in IVT with 10:1 thienoguanosine (^th G):GTP increased the percent incorporation and yield of 5'-^th G-initiated precursor tRNA for a net ≈threefold gain compared to IVT with wild-type T7RNAP. We also demonstrated that a one-pot multienzyme approach, consisting of transcription by T7RNAP P266L and post-transcriptional cleanup by polyphosphatase and an exonuclease, led to essentially near-homogeneous 5'-^th G-modified transcripts. This approach should be of broad utility in preparing 5'-modified RNAs.

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.

NMR monitoring of the SELEX process to confirm enrichment of structured RNA

RNA aptamers are RNA molecules that bind to a target molecule with high affinity and specificity using uniquely-folded tertiary structures. RNA aptamers are selected from an RNA pool typically comprising up to 1015 different sequences generated by iterative steps of selection and amplification known as Systematic Evolution of Ligands by EXponential enrichment (SELEX). Over several rounds of SELEX, the diversity of the RNA pool decreases and the aptamers are enriched. Hence, monitoring of the enrichment of these RNA pools is critical for the successful selection of aptamers, and several methods for monitoring them have been developed. In this study, we measured one-dimensional imino proton NMR spectra of RNA pools during SELEX. The spectrum of the initial RNA pool indicates that the RNAs adopt tertiary structures. The structural diversity of the RNA pools was shown to depend highly on the design of the primer-binding sequence. Furthermore, we demonstrate that enrichment of RNA aptamers can be monitored using NMR. The RNA pools can be recovered from the NMR tube after measurement of NMR spectra. We also can monitor target binding in the NMR tubes. Thus, we propose using NMR to monitor the enrichment of structured aptamers during the SELEX process.

Integrated structural biology to unravel molecular mechanisms of protein-RNA recognition

Recent advances in RNA sequencing technologies have greatly expanded our knowledge of the RNA landscape in cells, often with spatiotemporal resolution. These techniques identified many new (often non-coding) RNA molecules. Large-scale studies have also discovered novel RNA binding proteins (RBPs), which exhibit single or multiple RNA binding domains (RBDs) for recognition of specific sequence or structured motifs in RNA. Starting from these large-scale approaches it is crucial to unravel the molecular principles of protein-RNA recognition in ribonucleoprotein complexes (RNPs) to understand the underlying mechanisms of gene regulation. Structural biology and biophysical studies at highest possible resolution are key to elucidate molecular mechanisms of RNA recognition by RBPs and how conformational dynamics, weak interactions and cooperative binding contribute to the formation of specific, context-dependent RNPs. While large compact RNPs can be well studied by X-ray crystallography and cryo-EM, analysis of dynamics and weak interaction necessitates the use of solution methods to capture these properties. Here, we illustrate methods to study the structure and conformational dynamics of protein-RNA complexes in solution starting from the identification of interaction partners in a given RNP. Biophysical and biochemical techniques support the characterization of a protein-RNA complex and identify regions relevant in structural analysis. Nuclear magnetic resonance (NMR) is a powerful tool to gain information on folding, stability and dynamics of RNAs and characterize RNPs in solution. It provides crucial information that is complementary to the static pictures derived from other techniques. NMR can be readily combined with other solution techniques, such as small angle X-ray and/or neutron scattering (SAXS/SANS), electron paramagnetic resonance (EPR), and Förster resonance energy transfer (FRET), which provide information about overall shapes, internal domain arrangements and dynamics. Principles of protein-RNA recognition and current approaches are reviewed and illustrated with recent studies.

NMR Structural Profiling of Transcriptional Intermediates Reveals Riboswitch Regulation by Metastable RNA Conformations

Gene repression induced by the formation of transcriptional terminators represents a prime example for the coupling of RNA synthesis, folding, and regulation. In this context, mapping the changes in available conformational space of transcription intermediates during RNA synthesis is important to understand riboswitch function. A majority of riboswitches, an important class of small metabolite-sensing regulatory RNAs, act as transcriptional regulators, but the dependence of ligand binding and the subsequent allosteric conformational switch on mRNA transcript length has not yet been investigated. We show a strict fine-tuning of binding and sequence-dependent alterations of conformational space by structural analysis of all relevant transcription intermediates at single-nucleotide resolution for the I-A type 2'dG-sensing riboswitch from Mesoplasma florum by NMR spectroscopy. Our results provide a general framework to dissect the coupling of synthesis and folding essential for riboswitch function, revealing the importance of metastable states for RNA-based gene regulation.

Transcriptional start site heterogeneity modulates the structure and function of the HIV-1 genome

The promoter in HIV type 1 (HIV-1) proviral DNA contains three sequential guanosines at the U3-R boundary that have been proposed to function as sites for transcription initiation. Here we show that all three sites are used in cells infected with HIV-1 and that viral RNAs containing a single 5' capped guanosine (^Cap1G) are specifically selected for packaging in virions, consistent with a recent report [Masuda et al. (2015) Sci Rep 5:17680]. In addition, we now show that transcripts that begin with two or three capped guanosines (^Cap2G or ^Cap3G) are enriched on polysomes, indicating that RNAs synthesized from different transcription start sites have different functions in viral replication. Because genomes are selected for packaging as dimers, we examined the in vitro monomer-dimer equilibrium properties of ^Cap1G, ^Cap2G, and ^Cap3G 5'-leader RNAs in the NL4-3 strain of HIV-1. Strikingly, under physiological-like ionic conditions in which the ^Cap1G 5'-leader RNA adopts a dimeric structure, the ^Cap2G and ^Cap3G 5'-leader RNAs exist predominantly as monomers. Mutagenesis studies designed to probe for base-pairing interactions suggest that the additional guanosines of the 2G and 3G RNAs remodel the base of the PolyA hairpin, resulting in enhanced sequestration of dimer-promoting residues and stabilization of the monomer. Our studies suggest a mechanism through which the structure, function, and fate of the viral genome can be modulated by the transcriptionally controlled presence or absence of a single 5' guanosine.

NMR detection of intermolecular interaction sites in the dimeric 5'-leader of the HIV-1 genome

HIV type-1 (HIV-1) contains a pseudodiploid RNA genome that is selected for packaging and maintained in virions as a noncovalently linked dimer. Genome dimerization is mediated by conserved elements within the 5'-leader of the RNA, including a palindromic dimer initiation signal (DIS) that has been proposed to form kissing hairpin and/or extended duplex intermolecular contacts. Here, we have applied a ²H-edited NMR approach to directly probe for intermolecular interactions in the full-length, dimeric HIV-1 5'-leader (688 nucleotides; 230 kDa). The interface is extensive and includes DIS:DIS base pairing in an extended duplex state as well as intermolecular pairing between elements of the upstream Unique-5' (U5) sequence and those near the gag start site (AUG). Other pseudopalindromic regions of the leader, including the transcription activation (TAR), polyadenylation (PolyA), and primer binding (PBS) elements, do not participate in intermolecular base pairing. Using a ²H-edited one-dimensional NMR approach, we also show that the extended interface structure forms on a time scale similar to that of overall RNA dimerization. Our studies indicate that a kissing dimer-mediated structure, if formed, exists only transiently and readily converts to the extended interface structure, even in the absence of the HIV-1 nucleocapsid protein or other RNA chaperones.

A general method for rapid and cost-efficient large-scale production of 5' capped RNA

The eukaryotic mRNA 5' cap structure is indispensible for pre-mRNA processing, mRNA export, translation initiation, and mRNA stability. Despite this importance, structural and biophysical studies that involve capped RNA are challenging and rare due to the lack of a general method to prepare mRNA in sufficient quantities. Here, we show that the vaccinia capping enzyme can be used to produce capped RNA in the amounts that are required for large-scale structural studies. We have therefore designed an efficient expression and purification protocol for the vaccinia capping enzyme. Using this approach, the reaction scale can be increased in a cost-efficient manner, where the yields of the capped RNA solely depend on the amount of available uncapped RNA target. Using a large number of RNA substrates, we show that the efficiency of the capping reaction is largely independent of the sequence, length, and secondary structure of the RNA, which makes our approach generally applicable. We demonstrate that the capped RNA can be directly used for quantitative biophysical studies, including fluorescence anisotropy and high-resolution NMR spectroscopy. In combination with (13)C-methyl-labeled S-adenosyl methionine, the methyl groups in the RNA can be labeled for methyl TROSY NMR spectroscopy. Finally, we show that our approach can produce both cap-0 and cap-1 RNA in high amounts. In summary, we here introduce a general and straightforward method that opens new means for structural and functional studies of proteins and enzymes in complex with capped RNA.

Chemo-enzymatic labeling for rapid assignment of RNA molecules

Even though Nuclear Magnetic Resonance (NMR) spectroscopy is one of the few techniques capable of determining atomic resolution structures of RNA, it is constrained by two major problems of chemical shift overlap of resonances and rapid signal loss due to line broadening. Emerging tools to tackle these problems include synthesis of atom specifically labeled or chemically modified nucleotides. Herein we review the synthesis of these nucleotides, the design and production of appropriate RNA samples, and the application and analysis of the NMR experiments that take advantage of these labels.

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

In RNA secondary structure determination, it is essential to determine whether a nucleotide is base-paired and not. Base-pairing of nucleotides is mediated by hydrogen bonds. The NMR characterization of hydrogen bonds relies on experiments correlating the NMR resonances of exchangeable protons and can be best performed for structured parts of the RNA, where labile hydrogen atoms are protected from solvent exchange. Functionally important regions in RNA, however, frequently reveal increased dynamic disorder which often leads to NMR signals of exchangeable protons that are broadened beyond (1)H detection. Here, we develop (13)C direct detected experiments to observe all nucleotides in RNA irrespective of whether they are involved in hydrogen bonds or not. Exploiting the self-decoupling of scalar couplings due to the exchange process, the hydrogen bonding behavior of the hydrogen bond donor of each individual nucleotide can be determined. Furthermore, the adaption of HNN-COSY experiments for (13)C direct detection allows correlations of donor-acceptor pairs and the localization of hydrogen-bond acceptor nucleotides. The proposed (13)C direct detected experiments therefore provide information about molecular sites not amenable by conventional proton-detected methods. Such information makes the RNA secondary structure determination by NMR more accurate and helps to validate secondary structure predictions based on bioinformatics.