SHAPE Directed Discovery of New Functions in Large RNAs

Ligand-binding pockets in RNA and where to find them

RNAs are critical regulators of gene expression, and their functions are often mediated by complex secondary and tertiary structures. Structured regions in RNA can selectively interact with small molecules-via well-defined ligand-binding pockets-to modulate the regulatory repertoire of an RNA. The broad potential to modulate biological function intentionally via RNA-ligand interactions remains unrealized, however, due to challenges in identifying compact RNA motifs with the ability to bind ligands with good physicochemical properties (often termed drug-like). Here, we devise fpocketR, a computational strategy that accurately detects pockets capable of binding drug-like ligands in RNA structures. Remarkably few, roughly 50, of such pockets have ever been visualized. We experimentally confirmed the ligandability of novel pockets detected with fpocketR using a fragment-based approach introduced here, Frag-MaP, that detects ligand-binding sites in cells. Analysis of pockets detected by fpocketR and validated by Frag-MaP reveals dozens of sites able to bind drug-like ligands, supports a model for RNA secondary structural motifs able to bind quality ligands, and creates a broad framework for understanding the RNA ligand-ome.

The Druggable Transcriptome Project: From Chemical Probes to Precision Medicines

RNA presents abundant opportunities as a drug target, offering significant potential for small molecule medicine development. The transcriptome, comprising both coding and noncoding RNAs, is a rich area for therapeutic innovation, yet challenges persist in targeting RNA with small molecules. RNA structure can be predicted with or without experimental data, but discrepancies with the actual biological structure can impede progress. Prioritizing RNA targets supported by genetic or evolutionary evidence enhances success. Further, small molecules must demonstrate binding to RNA in cells, not solely in vitro, to validate both the target and compound. Effective small molecule binders modulate functional sites that influence RNA biology, as binding to nonfunctional sites requires recruiting effector mechanisms, for example degradation, to achieve therapeutic outcomes. Addressing these challenges is critical to unlocking RNA's vast potential for small molecule medicines, and a strategic framework is proposed to navigate this promising field, with a focus on targeting human RNAs.

Ligand-binding pockets in RNA, and where to find them

RNAs are critical regulators of gene expression, and their functions are often mediated by complex secondary and tertiary structures. Structured regions in RNA can selectively interact with small molecules - via well-defined ligand binding pockets - to modulate the regulatory repertoire of an RNA. The broad potential to modulate biological function intentionally via RNA-ligand interactions remains unrealized, however, due to challenges in identifying compact RNA motifs with the ability to bind ligands with good physicochemical properties (often termed drug-like). Here, we devise fpocketR, a computational strategy that accurately detects pockets capable of binding drug-like ligands in RNA structures. Remarkably few, roughly 50, of such pockets have ever been visualized. We experimentally confirmed the ligandability of novel pockets detected with fpocketR using a fragment-based approach introduced here, Frag-MaP, that detects ligand-binding sites in cells. Analysis of pockets detected by fpocketR and validated by Frag-MaP reveals dozens of newly identified sites able to bind drug-like ligands, supports a model for RNA secondary structural motifs able to bind quality ligands, and creates a broad framework for understanding the RNA ligand-ome.

Scaffold-enabled high-resolution cryo-EM structure determination of RNA

Cryo-EM structure determination of protein-free RNAs has remained difficult with most attempts yielding low to moderate resolution and lacking nucleotide-level detail. These difficulties are compounded for small RNAs as cryo-EM is inherently more difficult for lower molecular weight macromolecules. Here we present a strategy for fusing small RNAs to a group II intron that yields high resolution structures of the appended RNA. We demonstrate this technology by determining the structures of the 86-nucleotide (nt) thiamine pyrophosphate (TPP) riboswitch aptamer domain and the recently described 210-nt raiA bacterial non-coding RNA involved in sporulation and biofilm formation. In the case of the TPP riboswitch aptamer domain, the scaffolding approach allowed visualization of the riboswitch ligand binding pocket at 2.5 Å resolution. We also determined the structure of the ligand-free apo state and observe that the aptamer domain of the riboswitch adopts an open Y-shaped conformation in the absence of ligand. Using this scaffold approach, we determined the structure of raiA at 2.5 Å in the core. Our versatile scaffolding strategy enables efficient RNA structure determination for a broad range of small to moderate-sized RNAs, which were previously intractable for high-resolution cryo-EM studies.

Visualizing RNA structure ensembles by single-molecule correlated chemical probing

RNA molecules fold to form complex internal structures. Many of these RNA structures populate ensembles with rheostat-like properties, with each state having a distinct function. Until recently, analysis of RNA structures, especially within cells, was limited to modeling either a single averaged structure or computationally-modeled ensembles. These approaches obscure the intrinsic heterogeneity of many structured RNAs. Single-molecule correlated chemical probing (smCCP) strategies are now making it possible to measure and deconvolute RNA structure ensembles based on efficiently executed chemical probing experiments. Here, we provide an overview of fundamental single-molecule probing principles, review current ensemble deconvolution strategies, and discuss recent applications to diverse biological systems. smCCP is enabling a revolution in understanding how the plasticity of RNA structure is exploited in biological systems to respond to stimuli and alter gene function. The energetics of RNA ensembles are often subtle and a subset can likely be targeted to modulate disease-associated biological processes.

Evaluating Performance of Different RNA Secondary Structure Prediction Programs Using Self-cleaving Ribozymes

Accurate identification of the correct, biologically relevant RNA structures is critical to understanding various aspects of RNA biology since proper folding represents the key to the functionality of all types of RNA molecules and plays pivotal roles in many essential biological processes. Thus, a plethora of approaches have been developed to predict, identify, or solve RNA structures based on various computational, molecular, genetic, chemical, or physicochemical strategies. Purely computational approaches hold distinct advantages over all other strategies in terms of the ease of implementation, time, speed, cost, and throughput, but they strongly underperform in terms of accuracy that significantly limits their broader application. Nonetheless, the advantages of these methods led to a steady development of multiple in silico RNA secondary structure prediction approaches including recent deep learning-based programs. Here, we compared the accuracy of predictions of biologically relevant secondary structures of dozens of self-cleaving ribozyme sequences using seven in silico RNA folding prediction tools with tasks of varying complexity. We found that while many programs performed well in relatively simple tasks, their performance varied significantly in more complex RNA folding problems. However, in general, a modern deep learning method outperformed the other programs in the complex tasks in predicting the RNA secondary structures, at least based on the specific class of sequences tested, suggesting that it may represent the future of RNA structure prediction algorithms.

Chemical diversity of reagents that modify RNA 2'-OH in water: a review

Electrophilic water-soluble compounds have proven versatile in reacting selectively with 2'-OH groups in RNA, enabling structure mapping, probing, caging, labeling, crosslinking, and conjugation of RNAs in vitro and in living cells. While early work focused on one or two types of reagents with limited properties, recent studies have greatly diversified the structure, properties, and applications of these reagents. Here we review the scope of documented RNA hydroxyl-reactive species reported to date, with an eye to the effects of chemical structure on reactivity with RNA and other useful properties. Multiple forms of carbonyl electrophiles are now known to react at the 2'-OH, and recently, sulfonyl and aryl electrophiles have also been documented to form bonds there in high yields as well. In addition to electrophilicity, data also point to significant effects of reagent stability, steric bulk, and chirality on reaction yields and selectivity. Finally, we outline reagent properties and principles that define utility in applications with RNA, with an eye to the design of future reagents.

Characterization of the dual regulation by a c-di-GMP riboswitch Bc1 with a long expression platform from Bacillus thuringiensis

A riboswitch generally regulates the expression of its downstream genes through conformational change in its expression platform (EP) upon ligand binding. The cyclic diguanosine monophosphate (c-di-GMP) class I riboswitch Bc1 is widespread and conserved among Bacillus cereus group species. In this study, we revealed that Bc1 has a long EP with two typical ρ-independent terminator sequences 28 bp apart. The upstream terminator T1 is dominant in vitro, while downstream terminator T2 is more efficient in vivo. Through mutation analysis, we elucidated that Bc1 exerts a rare and incoherent "transcription-translation" dual regulation with T2 playing a crucial role. However, we found that Bc1 did not respond to c-di-GMP under in vitro transcription conditions, and the expressions of downstream genes did not change with fluctuation in intracellular c-di-GMP concentration. To explore this puzzle, we conducted SHAPE-MaP and confirmed the interaction of Bc1 with c-di-GMP. This shows that as c-di-GMP concentration increases, T1 unfolds but T2 remains almost intact and functional. The presence of T2 masks the effect of T1 unwinding, resulting in no response of Bc1 to c-di-GMP. The high Shannon entropy values of EP region imply the potential alternative structures of Bc1. We also found that zinc uptake regulator can specifically bind to the dual terminator coding sequence and slightly trigger the response of Bc1 to c-di-GMP. This work will shed light on the dual-regulation riboswitch and enrich our understanding of the RNA world.IMPORTANCEIn nature, riboswitches are involved in a variety of metabolic regulation, most of which preferentially regulate transcription termination or translation initiation of downstream genes in specific ways. Alternatively, the same or different riboswitches can exist in tandem to enhance regulatory effects or respond to multiple ligands. However, many putative conserved riboswitches have not yet been experimentally validated. Here, we found that the c-di-GMP riboswitch Bc1 with a long EP could form a dual terminator and exhibit non-canonical and incoherent "transcription-translation" dual regulation. Besides, zinc uptake regulator specifically bound to the coding sequence of the Bc1 EP and slightly mediated the action of Bc1. The application of SHAPE-MaP to the dual regulation mechanism of Bc1 may establish the foundation for future studies of such complex untranslated regions in other bacterial genomes.

RNATACs: Multispecific small molecules targeting RNA by induced proximity

RNA-targeting small molecules (rSMs) have become an attractive modality to tackle traditionally undruggable proteins and expand the druggable space. Among many innovative concepts, RNA-targeting chimeras (RNATACs) represent a new class of multispecific, induced proximity small molecules that act by chemically bringing RNA targets into proximity with an endogenous RNA effector, such as a ribonuclease (RNase). Depending on the RNA effector, RNATACs can alter the stability, localization, translation, or splicing of the target RNA. Although still in its infancy, this new modality has the potential for broad applications in the future to treat diseases with high unmet need. In this review, we discuss potential advantages of RNATACs, recent progress in the field, and challenges to this cutting-edge technology.

Scaffold-enabled high-resolution cryo-EM structure determination of RNA

Cryo-EM structure determination of protein-free RNAs has remained difficult with most attempts yielding low to moderate resolution and lacking nucleotide-level detail. These difficulties are compounded for small RNAs as cryo-EM is inherently more difficult for lower molecular weight macromolecules. Here we present a strategy for fusing small RNAs to a group II intron that yields high resolution structures of the appended RNA, which we demonstrate with the 86-nucleotide thiamine pyrophosphate (TPP) riboswitch, and visualizing the riboswitch ligand binding pocket at 2.5 Å resolution. We also determined the structure of the ligand-free apo state and observe that the aptamer domain of the riboswitch undergoes a large-scale conformational change upon ligand binding, illustrating how small molecule binding to an RNA can induce large effects on gene expression. This study both sets a new standard for cryo-EM riboswitch visualization and offers a versatile strategy applicable to a broad range of small to moderate-sized RNAs, which were previously intractable for high-resolution cryo-EM studies.

RNAvigate: efficient exploration of RNA chemical probing datasets

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA, including local nucleotide flexibility, RNA secondary structure, protein and ligand binding, through-space interaction networks, and multistate structural ensembles. Deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multilayered relationships. Current platforms lack the broad accessibility, flexibility and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library that automatically parses 21 standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs 18 plot types. RNAvigate enables efficient exploration of nuanced relationships between multiple layers of RNA structure information and across multiple experimental conditions. Compatibility with Jupyter notebooks enables nonburdensome, reproducible, transparent and organized sharing of multistep analyses and data visualization strategies. RNAvigate simplifies and accelerates discovery and characterization of RNA-centric functions in biology.

Minimization of the E. coli ribosome, aided and optimized by community science

The ribosome is a ribonucleoprotein complex found in all domains of life. Its role is to catalyze protein synthesis, the messenger RNA (mRNA)-templated formation of amide bonds between α-amino acid monomers. Amide bond formation occurs within a highly conserved region of the large ribosomal subunit known as the peptidyl transferase center (PTC). Here we describe the step-wise design and characterization of mini-PTC 1.1, a 284-nucleotide RNA that recapitulates many essential features of the Escherichia coli PTC. Mini-PTC 1.1 folds into a PTC-like structure under physiological conditions, even in the absence of r-proteins, and engages small molecule analogs of A- and P-site tRNAs. The sequence of mini-PTC 1.1 differs from the wild type E. coli ribosome at 12 nucleotides that were installed by a cohort of citizen scientists using the on-line video game Eterna. These base changes improve both the secondary structure and tertiary folding of mini-PTC 1.1 as well as its ability to bind small molecule substrate analogs. Here, the combined input from Eterna citizen-scientists and RNA structural analysis provides a robust workflow for the design of a minimal PTC that recapitulates many features of an intact ribosome.

Small molecule approaches to targeting RNA

The development of innovative methodologies to identify RNA binders has attracted enormous attention in chemical biology and drug discovery. Although antibiotics targeting bacterial ribosomal RNA have been on the market for decades, the renewed interest in RNA targeting reflects the need to better understand complex intracellular processes involving RNA. In this context, small molecules are privileged tools used to explore the biological functions of RNA and to validate RNAs as therapeutic targets, and they eventually are to become new drugs. Despite recent progress, the rational design of specific RNA binders requires a better understanding of the interactions which occur with the RNA target to reach the desired biological response. In this Review, we discuss the challenges to approaching this underexplored chemical space, together with recent strategies to bind, interact and affect biologically relevant RNAs.

A compact regulatory RNA element in mouse Hsp70 mRNA

Hsp70 (70 kDa heat shock protein) performs molecular chaperone functions by assisting the folding of newly synthesized and misfolded proteins, thereby counteracting various cell stresses and preventing multiple diseases, including neurodegenerative disorders and cancers. It is well established that, immediately after heat shock, Hsp70 gene expression is mediated by a canonical mechanism of cap-dependent translation. However, the molecular mechanism of Hsp70 expression during heat shock remains elusive. Intriguingly, the 5' end of Hsp70 messenger RNA (mRNA) appears to form a compact structure with the potential to regulate protein expression in a cap-independent manner. Here, we determined the minimal length of the mHsp70 5'-terminal mRNA sequence that is required for RNA folding into a highly compact structure. This span of this RNA element was mapped and the secondary structure characterized by chemical probing, resulting in a secondary structural model that includes multiple stable stems, including one containing the canonical start codon. All of these components, including a short stretch of the 5' open reading frame (ORF), were shown to be vital for RNA folding. This work provides a structural basis for future investigations on the role of translational regulatory structures in the 5' untranslated region and ORF sequences of Hsp70 during heat shock.

Structural Analysis of the Hepatitis B Virus RNA Encapsidation Signal ε by Selective 2'-Hydroxyl Acylation Analyzed by Primer Extension (SHAPE)

RNA structure is crucial for RNA function, including in viral cis-elements such as the hepatitis B virus (HBV) RNA encapsidation signal ε. Interacting with the viral polymerase ε mediates packaging of the pregenomic (pg) RNA into capsids, initiation of reverse transcription, and it affects the mRNA functions of pgRNA. As free RNA, the 61-nucleotide (nt) ε sequence adopts a bipartite stem-loop structure with a central bulge and an apical loop. Due to stable Watson-Crick base pairing, this was already predicted by early RNA folding programs and confirmed by classical enzymatic and chemical structure probing. A newer, high-resolution probing technique exploits the selective acylation of solvent-accessible 2'-hydroxyls in the RNA backbone by electrophilic compounds such as 2-methylnicotinic acid imidazolide (NAI), followed by mapping of the modified sites by primer extension. This SHAPE principle has meanwhile been extended to numerous applications. Here we provide a basic protocol for NAI-based SHAPE of isolated HBV ε RNA which already provided insights into the impact of mutations, and preliminarily, of polymerase binding on the RNA structural dynamics. While the focus is on NAI modification, we also briefly cover target RNA preparation by in vitro transcription, primer extension using a radiolabeled primer, and analysis of the resulting cDNAs by denaturing polyacrylamide gelelectrophoresis (PAGE). Given the high tolerance of SHAPE chemistry to different conditions, including applicability in live cells, we expect this technique to greatly facilitate deciphering the conformational dynamics underlying the various functions of the ε element, especially in concert with the recently solved three-dimensional structure of the free RNA.

RNAvigate: Efficient exploration of RNA chemical probing datasets

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA including local nucleotide flexibility, RNA secondary structure, protein- and ligand-binding, through-space interaction networks, and multi-state structural ensembles. Performing these experiments, by themselves, does not directly lead to biological insight. Instead, deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multi-layered relationships. Current platforms lack the broad accessibility, flexibility, and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library. RNAvigate currently automatically parses twenty-one standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs eighteen plot types. These features enable efficient exploration of nuanced relationships between chemical probing data, RNA structure, and motif annotations across multiple experimental samples. Compatibility with Jupyter Notebooks enables non-burdensome, reproducible, transparent and organized sharing of multi-step analyses and data visualization strategies. RNAvigate simplifies examination of multi-layered RNA structure information and accelerates discovery and characterization of RNA-centric functions in biology.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a pivotal role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N1-methylguanine DMS modifications to enable high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than current two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are straightforward to perform and will broadly facilitate improved RNA structural analysis in living cells.

Analysis of lncRNAs in Lupinus mutabilis (Tarwi) and Their Potential Role in Drought Response

Lupinus mutabilis is a legume with high agronomic potential and available transcriptomic data for which lncRNAs have not been studied. Therefore, our objective was to identify, characterize, and validate the drought-responsive lncRNAs in L. mutabilis. To achieve this, we used a multilevel approach based on lncRNA prediction, annotation, subcellular location, thermodynamic characterization, structural conservation, and validation. Thus, 590 lncRNAs were identified by at least two algorithms of lncRNA identification. Annotation with the PLncDB database showed 571 lncRNAs unique to tarwi and 19 lncRNAs with homology in 28 botanical families including Solanaceae (19), Fabaceae (17), Brassicaceae (17), Rutaceae (17), Rosaceae (16), and Malvaceae (16), among others. In total, 12 lncRNAs had homology in more than 40 species. A total of 67% of lncRNAs were located in the cytoplasm and 33% in exosomes. Thermodynamic characterization of S03 showed a stable secondary structure with -105.67 kcal/mol. This structure included three regions, with a multibranch loop containing a hairpin with a SECIS-like element. Evaluation of the structural conservation by CROSSalign revealed partial similarities between L. mutabilis (S03) and S. lycopersicum (Solyc04r022210.1). RT-PCR validation demonstrated that S03 was upregulated in a drought-tolerant accession of L. mutabilis. Finally, these results highlighted the importance of lncRNAs in tarwi improvement under drought conditions.

Mutation signature filtering enables high-fidelity RNA structure probing at all four nucleobases with DMS

Chemical probing experiments have transformed RNA structure analysis, enabling high-throughput measurement of base-pairing in living cells. Dimethyl sulfate (DMS) is one of the most widely used structure probing reagents and has played a prominent role in enabling next-generation single-molecule probing analyses. However, DMS has traditionally only been able to probe adenine and cytosine nucleobases. We previously showed that, using appropriate conditions, DMS can also be used to interrogate base-pairing of uracil and guanines in vitro at reduced accuracy. However, DMS remained unable to informatively probe guanines in cells. Here, we develop an improved DMS mutational profiling (MaP) strategy that leverages the unique mutational signature of N 1 -methylguanine DMS modifications to enable robust, high-fidelity structure probing at all four nucleotides, including in cells. Using information theory, we show that four-base DMS reactivities convey greater structural information than comparable two-base DMS and SHAPE probing strategies. Four-base DMS experiments further enable improved direct base-pair detection by single-molecule PAIR analysis, and ultimately support RNA structure modeling at superior accuracy. Four-base DMS probing experiments are easily performed and will broadly facilitate improved RNA structural analysis in living cells.

Structure-first identification of RNA elements that regulate dengue virus genome architecture and replication

The genomes of RNA viruses encode the information required for replication in host cells both in their linear sequence and in complex higher-order structures. A subset of these RNA genome structures show clear sequence conservation, and have been extensively described for well-characterized viruses. However, the extent to which viral RNA genomes contain functional structural elements-unable to be detected by sequence alone-that nonetheless are critical to viral fitness is largely unknown. Here, we devise a structure-first experimental strategy and use it to identify 22 structure-similar motifs across the coding sequences of the RNA genomes for the four dengue virus serotypes. At least 10 of these motifs modulate viral fitness, revealing a significant unnoticed extent of RNA structure-mediated regulation within viral coding sequences. These viral RNA structures promote a compact global genome architecture, interact with proteins, and regulate the viral replication cycle. These motifs are also thus constrained at the levels of both RNA structure and protein sequence and are potential resistance-refractory targets for antivirals and live-attenuated vaccines. Structure-first identification of conserved RNA structure enables efficient discovery of pervasive RNA-mediated regulation in viral genomes and, likely, other cellular RNAs.

Single-Molecule Correlated Chemical Probing: A Revolution in RNA Structure Analysis

RNA molecules convey biological information both in their linear sequence and in their base-paired secondary and tertiary structures. Chemical probing experiments, which involve treating an RNA with a reagent that modifies conformationally dynamic nucleotides, have broadly enabled examination of short- and long-range RNA structure in diverse contexts, including in living cells. For decades, chemical probing experiments have been interpreted in a per-nucleotide way, such that the reactivity measured at each nucleotide reports the average structure at a position over all RNA molecules within a sample. However, there are numerous important cases where per-nucleotide chemical probing falls short, including for RNAs that are bound by proteins, RNAs that form complex higher order structures, and RNAs that sample multiple conformations.Recent experimental and computational innovations have started a revolution in RNA structure analysis by transforming chemical probing into a massively parallel, single-molecule experiment. Enabled by a specialized reverse transcription strategy called mutational profiling (MaP), multiple chemical modification events can be measured within individual RNA molecules. Nucleotides that communicate structurally through direct base pairing or large-scale folding-unfolding transitions will react with chemical probes in a correlated manner, thereby revealing structural complexity hidden to conventional approaches. These single-molecule correlated chemical probing (smCCP) experiments can be interpreted to directly identify nucleotides that base pair (the PAIR-MaP strategy) and to reveal long-range, through-space structural communication (RING-MaP). Correlated probing can also define the thermodynamic populations of complex RNA ensembles (DANCE-MaP). Complex RNA-protein networks can be interrogated by cross-linking proteins to RNA and measuring correlations between cross-linked positions (RNP-MaP).smCCP thus visualizes RNA secondary and higher-order structure with unprecedented accuracy, defining novel structures, RNA-protein interaction networks, time-resolved dynamics, and allosteric structural switches. These strategies are not mutually exclusive; in favorable cases, multiple levels of RNA structure ─ base pairing, through-space structural communication, and equilibrium ensembles ─ can be resolved concurrently. The physical experimentation required for smCCP is profoundly simple, and experiments are readily performed in cells on RNAs of any size, including large noncoding RNAs and mRNAs. Single-molecule correlated chemical probing is paving the way for a new generation of biophysical studies on RNA in living systems.

RSCanner: rapid assessment and visualization of RNA structure content

MOTIVATION

The increasing availability of RNA structural information that spans many kilobases of transcript sequence imposes a need for tools that can rapidly screen, identify, and prioritize structural modules of interest.

RESULTS

We describe RNA Structural Content Scanner (RSCanner), an automated tool that scans RNA transcripts for regions that contain high levels of secondary structure and then classifies each region for its relative propensity to adopt stable or dynamic structures. RSCanner then generates an intuitive heatmap enabling users to rapidly pinpoint regions likely to contain a high or low density of discrete RNA structures, thereby informing downstream functional or structural investigation.

AVAILABILITY AND IMPLEMENTATION

RSCanner is freely available as both R script and R Markdown files, along with full documentation and test data (https://github.com/pylelab/RSCanner).

Integrating experimental data with molecular simulations to investigate RNA structural dynamics

Conformational dynamics is crucial for ribonucleic acid (RNA) function. Techniques such as nuclear magnetic resonance, cryo-electron microscopy, small- and wide-angle X-ray scattering, chemical probing, single-molecule Förster resonance energy transfer, or even thermal or mechanical denaturation experiments probe RNA dynamics at different time and space resolutions. Their combination with accurate atomistic molecular dynamics (MD) simulations paves the way for quantitative and detailed studies of RNA dynamics. First, experiments provide a quantitative validation tool for MD simulations. Second, available data can be used to refine simulated structural ensembles to match experiments. Finally, comparison with experiments allows for improving MD force fields that are transferable to new systems for which data is not available. Here we review the recent literature and provide our perspective on this field.

Probing RNA Conformations Using a Polymer-Electrolyte Solid-State Nanopore

Nanopore systems have emerged as a leading platform for the analysis of biomolecular complexes with single-molecule resolution. The conformation of biomolecules, such as RNA, is highly dependent on the electrolyte composition, but solid-state nanopore systems often require high salt concentration to operate, precluding analysis of macromolecular conformations under physiologically relevant conditions. Here, we report the implementation of a polymer-electrolyte solid-state nanopore system based on alkali metal halide salts dissolved in 50% w/v poly(ethylene) glycol (PEG) to augment the performance of our system. We show that polymer-electrolyte bath governs the translocation dynamics of the analyte which correlates with the physical properties of the salt used in the bath. This allowed us to identify CsBr as the optimal salt to complement PEG to generate the largest signal enhancement. Harnessing the effects of the polymer-electrolyte, we probed the conformations of the Chikungunya virus (CHIKV) RNA genome fragments under physiologically relevant conditions. Our system was able to fingerprint CHIKV RNA fragments ranging from ∼300 to ∼2000 nt length and subsequently distinguish conformations between the co-transcriptionally folded and the natively refolded ∼2000 nt CHIKV RNA. We envision that the polymer-electrolyte solid-state nanopore system will further enable structural and conformational analyses of individual biomolecules under physiologically relevant conditions.

Global 5'-UTR RNA structure regulates translation of a SERPINA1 mRNA

SERPINA1 mRNAs encode the protease inhibitor α-1-antitrypsin and are regulated through post-transcriptional mechanisms. α-1-antitrypsin deficiency leads to chronic obstructive pulmonary disease (COPD) and liver cirrhosis, and specific variants in the 5'-untranslated region (5'-UTR) are associated with COPD. The NM_000295.4 transcript is well expressed and translated in lung and blood and features an extended 5'-UTR that does not contain a competing upstream open reading frame (uORF). We show that the 5'-UTR of NM_000295.4 folds into a well-defined multi-helix structural domain. We systematically destabilized mRNA structure across the NM_000295.4 5'-UTR, and measured changes in (SHAPE quantified) RNA structure and cap-dependent translation relative to a native-sequence reporter. Surprisingly, despite destabilizing local RNA structure, most mutations either had no effect on or decreased translation. Most structure-destabilizing mutations retained native, global 5'-UTR structure. However, those mutations that disrupted the helix that anchors the 5'-UTR domain yielded three groups of non-native structures. Two of these non-native structure groups refolded to create a stable helix near the translation initiation site that decreases translation. Thus, in contrast to the conventional model that RNA structure in 5'-UTRs primarily inhibits translation, complex folding of the NM_000295.4 5'-UTR creates a translation-optimized message by promoting accessibility at the translation initiation site.

Cell Compartment-Specific Folding of Ty1 Long Terminal Repeat Retrotransposon RNA Genome

The structural transitions RNAs undergo during trafficking are not well understood. Here, we used the well-developed yeast Ty1 retrotransposon to provide the first structural model of genome (g) RNA in the nucleus from a retrovirus-like transposon. Through a detailed comparison of nuclear Ty1 gRNA structure with those established in the cytoplasm, virus-like particles (VLPs), and those synthesized in vitro, we detected Ty1 gRNA structural alterations that occur during retrotransposition. Full-length Ty1 gRNA serves as the mRNA for Gag and Gag-Pol proteins and as the genome that is reverse transcribed within VLPs. We show that about 60% of base pairs predicted for the nuclear Ty1 gRNA appear in the cytoplasm, and active translation does not account for such structural differences. Most of the shared base pairs are represented by short-range interactions, whereas the long-distance pairings seem unique for each compartment. Highly structured motifs tend to be preserved after nuclear export of Ty1 gRNA. In addition, our study highlights the important role of Ty1 Gag in mediating critical RNA-RNA interactions required for retrotransposition.

ShapeSorter: a fully probabilistic method for detecting conserved RNA structure features supported by SHAPE evidence

There is an increased interest in the determination of RNA structures in vivo as it is now possible to probe them in a high-throughput manner, e.g. using SHAPE protocols. By now, there exist a range of computational methods that integrate experimental SHAPE-probing evidence into computational RNA secondary structure prediction. The state-of-the-art in this field is currently provided by computational methods that employ the minimum-free energy strategy for prediction RNA secondary structures with SHAPE-probing evidence. These methods, however, rely on the assumption that transcripts in vivo fold into the thermodynamically most stable configuration and ignore evolutionary evidence for conserved RNA structure features. We here present a new computational method, ShapeSorter, that predicts RNA structure features without employing the thermodynamic strategy. Instead, ShapeSorter employs a fully probabilistic framework to identify RNA structure features that are supported by evolutionary and SHAPE-probing evidence. Our method can capture RNA structure heterogeneity, pseudo-knotted RNA structures as well as transient and mutually exclusive RNA structure features. Moreover, it estimates P-values for the predicted RNA structure features which allows for easy filtering and ranking. We investigate the merits of our method in a comprehensive performance benchmarking and conclude that ShapeSorter has a significantly superior performance for predicting base-pairs than the existing state-of-the-art methods.

Probing RNA structures and functions by solvent accessibility: an overview from experimental and computational perspectives

Characterizing RNA structures and functions have mostly been focused on 2D, secondary and 3D, tertiary structures. Recent advances in experimental and computational techniques for probing or predicting RNA solvent accessibility make this 1D representation of tertiary structures an increasingly attractive feature to explore. Here, we provide a survey of these recent developments, which indicate the emergence of solvent accessibility as a simple 1D property, adding to secondary and tertiary structures for investigating complex structure-function relations of RNAs.

Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP

7SK is a conserved noncoding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFbresponsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription and establishes the power of DANCE-MaP to define RNA dynamics in cells.

The In Vivo and In Vitro Architecture of the Hepatitis C Virus RNA Genome Uncovers Functional RNA Secondary and Tertiary Structures

Hepatitis C virus (HCV) is a positive-strand RNA virus that remains one of the main contributors to chronic liver disease worldwide. Studies over the last 30 years have demonstrated that HCV contains a highly structured RNA genome and many of these structures play essential roles in the HCV life cycle. Despite the importance of riboregulation in this virus, most of the HCV RNA genome remains functionally unstudied. Here, we report a complete secondary structure map of the HCV RNA genome in vivo, which was studied in parallel with the secondary structure of the same RNA obtained in vitro. Our results show that HCV is folded extensively in the cellular context. By performing comprehensive structural analyses on both in vivo data and in vitro data, we identify compact and conserved secondary and tertiary structures throughout the genome. Genetic and evolutionary functional analyses demonstrate that many of these elements play important roles in the virus life cycle. In addition to providing a comprehensive map of RNA structures and riboregulatory elements in HCV, this work provides a resource for future studies aimed at identifying therapeutic targets and conducting further mechanistic studies on this important human pathogen. IMPORTANCE HCV has one of the most highly structured RNA genomes studied to date, and it is a valuable model system for studying the role of RNA structure in protein-coding genes. While previous studies have identified individual cases of regulatory RNA structures within the HCV genome, the full-length structure of the HCV genome has not been determined in vivo. Here, we present the complete secondary structure map of HCV determined both in cells and from corresponding transcripts generated in vitro. In addition to providing a comprehensive atlas of functional secondary structural elements throughout the genomic RNA, we identified a novel set of tertiary interactions and demonstrated their functional importance. In terms of broader implications, the pipeline developed in this study can be applied to other long RNAs, such as long noncoding RNAs. In addition, the RNA structural motifs characterized in this study broaden the repertoire of known riboregulatory elements.

The epitranscriptome toolbox

In the last decade, the notion that mRNA modifications are involved in regulation of gene expression was demonstrated in thousands of studies. To date, new technologies and methods allow accurate identification, transcriptome-wide mapping, and functional characterization of a growing number of RNA modifications, providing important insights into the biology of these marks. Most of the methods and approaches were developed for studying m⁶A, the most prevalent internal mRNA modification. However, unique properties of other RNA modifications stimulated the development of additional approaches. In this technical primer, we will discuss the available tools and approaches for detecting and studying different RNA modifications.

The multiple molecular dimensions of long noncoding RNAs that regulate gene expression and tumorigenesis

PURPOSE OF REVIEW

LncRNAs are emerging as key regulators of gene expression and they ensure homeostasis during cell differentiation and development, replication, and adaptation to the environment. Because of their key central role in regulating the biology of living cells, it is crucial to characterize how lncRNAs function at the genetic, transcriptomic, and mechanistic level.

RECENT FINDINGS

The low endogenous abundance and high molecular complexity of lncRNAs pose unique challenges for their characterization but new methodological advances in biochemistry, biophysics and cell biology have recently made it possible to characterize an increasing number of these transcripts, including oncogenic and tumor suppressor lncRNAs. These recent studies specifically address important issues that had remained controversial, such as the selectivity of lncRNA mechanisms of action, the functional importance of lncRNA sequences, secondary and tertiary structures, and the specificity of lncRNA interactions with proteins.

SUMMARY

These recent achievements, coupled to population-wide medical and genomic approaches that connect lncRNAs with human diseases and to recent advances in RNA-targeted drug development, open unprecedented new perspectives for exploiting lncRNAs as pharmacological targets or biomarkers to monitor and cure cancer, in addition to metabolic, developmental and cardiovascular diseases.

Distinct MUNC lncRNA structural domains regulate transcription of different promyogenic factors

Many lncRNAs have been discovered using transcriptomic data; however, it is unclear what fraction of lncRNAs is functional and what structural properties affect their phenotype. MUNC lncRNA (also known as ^DRReRNA) acts as an enhancer RNA for the Myod1 gene in cis and stimulates the expression of other promyogenic genes in trans by recruiting the cohesin complex. Here, experimental probing of the RNA structure revealed that MUNC contains multiple structural domains not detected by prediction algorithms in the absence of experimental information. We show that these specific and structurally distinct domains are required for induction of promyogenic genes, for binding genomic sites and gene expression regulation, and for binding the cohesin complex. Myod1 induction and cohesin interaction comprise only a subset of MUNC phenotype. Our study reveals unexpectedly complex, structure-driven functions for the MUNC lncRNA and emphasizes the importance of experimentally determined structures for understanding structure-function relationships in lncRNAs.

Direct detection of RNA modifications and structure using single-molecule nanopore sequencing

Modifications are present on many classes of RNA, including tRNA, rRNA, and mRNA. These modifications modulate diverse biological processes such as genetic recoding and mRNA export and folding. In addition, modifications can be introduced to RNA molecules using chemical probing strategies that reveal RNA structure and dynamics. Many methods exist to detect RNA modifications by short-read sequencing; however, limitations on read length inherent to short-read-based methods dissociate modifications from their native context, preventing single-molecule modification analysis. Here, we demonstrate direct RNA nanopore sequencing to detect endogenous and exogenous RNA modifications on long RNAs at the single-molecule level. We detect endogenous 2'-O-methyl and base modifications across E. coli and S. cerevisiae ribosomal RNAs as shifts in current signal and dwell times distally through interactions with the helicase motor protein. We further use the 2'-hydroxyl reactive SHAPE reagent acetylimidazole to probe RNA structure at the single-molecule level with readout by direct nanopore sequencing.

Stephenson et al. employ direct RNA nanopore sequencing to detect endogenous and exogenous modifications on single RNA molecules. The authors demonstrate detection of endogenous 2'-O-methylation (Nm) on native ribosomal RNAs, confirming known modification patterns. They describe the development of nanoSHAPE, a method that involves exogenously labeling RNA with a small-adduct-generating chemical probe that can reveal RNA structure using long-read sequencing.

Identifying proximal RNA interactions from cDNA-encoded crosslinks with ShapeJumper

SHAPE-JuMP is a concise strategy for identifying close-in-space interactions in RNA molecules. Nucleotides in close three-dimensional proximity are crosslinked with a bi-reactive reagent that covalently links the 2'-hydroxyl groups of the ribose moieties. The identities of crosslinked nucleotides are determined using an engineered reverse transcriptase that jumps across crosslinked sites, resulting in a deletion in the cDNA that is detected using massively parallel sequencing. Here we introduce ShapeJumper, a bioinformatics pipeline to process SHAPE-JuMP sequencing data and to accurately identify through-space interactions, as observed in complex JuMP datasets. ShapeJumper identifies proximal interactions with near-nucleotide resolution using an alignment strategy that is optimized to tolerate the unique non-templated reverse-transcription profile of the engineered crosslink-traversing reverse-transcriptase. JuMP-inspired strategies are now poised to replace adapter-ligation for detecting RNA-RNA interactions in most crosslinking experiments.

Piercing the fog of the RNA structure-ome

Machine learning is poised to transform RNA structure and function discovery

Direct Mapping of Higher-Order RNA Interactions by SHAPE-JuMP

Higher-order structure governs function for many RNAs. However, discerning this structure for large RNA molecules in solution is an unresolved challenge. Here, we present SHAPE-JuMP (selective 2'-hydroxyl acylation analyzed by primer extension and juxtaposed merged pairs) to interrogate through-space RNA tertiary interactions. A bifunctional small molecule is used to chemically link proximal nucleotides in an RNA structure. The RNA cross-link site is then encoded into complementary DNA (cDNA) in a single, direct step using an engineered reverse transcriptase that "jumps" across cross-linked nucleotides. The resulting cDNAs contain a deletion relative to the native RNA sequence, which can be detected by sequencing, that indicates the sites of cross-linked nucleotides. SHAPE-JuMP measures RNA tertiary structure proximity concisely across large RNA molecules at nanometer resolution. SHAPE-JuMP is especially effective at measuring interactions in multihelix junctions and loop-to-helix packing, enables modeling of the global fold for RNAs up to several hundred nucleotides in length, facilitates ranking of structural models by consistency with through-space restraints, and is poised to enable solution-phase structural interrogation and modeling of complex RNAs.