Genome-Wide Double-Stranded RNA Sequencing Reveals the Functional Significance of Base-Paired RNAs in Arabidopsis

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.     

Cation-dependent cleavage of the duplex form of the subtype-B HIV-1 RNA dimerization initiation site

The crystal structure of subtype-B HIV-1 genomic RNA Dimerization Initiation Site duplex revealed chain cleavage at a specific position resulting in 3′-phosphate and 5′-hydroxyl termini. A crystallographic analysis showed that Ba²⁺, Mn²⁺, Co²⁺ and Zn²⁺ bind specifically on a guanine base close to the cleaved position. The crystal structures also point to a necessary conformational change to induce an ‘in-line’ geometry at the cleavage site. In solution, divalent cations increased the rate of cleavage with pH/pKa compensation, indicating that a cation-bound hydroxide anion is responsible for the cleavage. We propose a ‘Trojan horse’ mechanism, possibly of general interest, wherein a doubly charged cation hosted near the cleavage site as a ‘harmless’ species is further transformed in situ into an ‘aggressive’ species carrying a hydroxide anion.     

Exploring the repertoire of RNA secondary motifs using graph theory; implications for RNA design

Understanding the structural repertoire of RNA is crucial for RNA genomics research. Yet current methods for finding novel RNAs are limited to small or known RNA families. To expand known RNA structural motifs, we develop a two-dimensional graphical representation approach for describing and estimating the size of RNA’s secondary structural repertoire, including naturally occurring and other possible RNA motifs. We employ tree graphs to describe RNA tree motifs and more general (dual) graphs to describe both RNA tree and pseudoknot motifs. Our estimates of RNA’s structural space are vastly smaller than the nucleotide sequence space, suggesting a new avenue for finding novel RNAs. Specifically our survey shows that known RNA trees and pseudoknots represent only a small subset of all possible motifs, implying that some of the ‘missing’ motifs may represent novel RNAs. To help pinpoint RNA-like motifs, we show that the motifs of existing functional RNAs are clustered in a narrow range of topological characteristics. We also illustrate the applications of our approach to the design of novel RNAs and automated comparison of RNA structures; we report several occurrences of RNA motifs within larger RNAs. Thus, our graph theory approach to RNA structures has implications for RNA genomics, structure analysis and design.     

Model-Free RNA Sequence and Structure Alignment Informed by SHAPE Probing Reveals a Conserved Alternate Secondary Structure for 16S rRNA

Discovery and characterization of functional RNA structures remains challenging due to deficiencies in de novo secondary structure modeling. Here we describe a dynamic programming approach for model-free sequence comparison that incorporates high-throughput chemical probing data. Based on SHAPE probing data alone, ribosomal RNAs (rRNAs) from three diverse organisms – the eubacteria E. coli and C. difficile and the archeon H. volcanii – could be aligned with accuracies comparable to alignments based on actual sequence identity. When both base sequence identity and chemical probing reactivities were considered together, accuracies improved further. Derived sequence alignments and chemical probing data from protein-free RNAs were then used as pseudo-free energy constraints to model consensus secondary structures for the 16S and 23S rRNAs. There are critical differences between these experimentally-informed models and currently accepted models, including in the functionally important neck and decoding regions of the 16S rRNA. We infer that the 16S rRNA has evolved to undergo large-scale changes in base pairing as part of ribosome function. As high-quality RNA probing data become widely available, structurally-informed sequence alignment will become broadly useful for de novo motif and function discovery.     

New scoring system to identify RNA G-quadruplex folding

G-quadruplexes (G4s) are non-canonical structures involved in many important cellular processes. To date, the prediction of potential G-quadruplex structures (PG4s) has been based almost exclusively on the sequence of interest agreeing with the algorithm G_x-N-_1–7-G_x-N_1–7-G_x-N_1–7-G_x (where x ≥ 3 and N = A, U, G or C). However, many sequences agreeing with this algorithm do not form G4s and are considered false-positive predictions. Here we show the RNA PG4 candidate in the 3′-untranslated region (UTR) of the TTYH1 gene to be one such false positive. Specifically, G4 folding was observed to be inhibited by the presence of multiple-cytosine tracks, located in the candidate’s genomic context, that adopted a Watson–Crick base-paired structure. Clearly, the neighbouring sequences of a PG4 may influence its folding. The secondary structure of 12 PG4 motifs along with either 15 or 50 nucleotides of their upstream and downstream genomic contexts were evaluated by in-line probing. Data permitted the development of a scoring system for the prediction of PG4s taking into account the effect of the neighbouring sequences. The accuracy of this scoring system was assessed by probing 14 other novel PG4 candidates retrieved in human 5′-UTRs. This new scoring system can be used, in combination with the standard algorithm, to better predict the folding of RNA G4s.     

Theoretical basis for stabilizing messenger RNA through secondary structure design

RNA hydrolysis presents problems in manufacturing, long-term storage, world-wide delivery and in vivo stability of messenger RNA (mRNA)-based vaccines and therapeutics. A largely unexplored strategy to reduce mRNA hydrolysis is to redesign RNAs to form double-stranded regions, which are protected from in-line cleavage and enzymatic degradation, while coding for the same proteins. The amount of stabilization that this strategy can deliver and the most effective algorithmic approach to achieve stabilization remain poorly understood. Here, we present simple calculations for estimating RNA stability against hydrolysis, and a model that links the average unpaired probability of an mRNA, or AUP, to its overall hydrolysis rate. To characterize the stabilization achievable through structure design, we compare AUP optimization by conventional mRNA design methods to results from more computationally sophisticated algorithms and crowdsourcing through the OpenVaccine challenge on the Eterna platform. We find that rational design on Eterna and the more sophisticated algorithms lead to constructs with low AUP, which we term ‘superfolder’ mRNAs. These designs exhibit a wide diversity of sequence and structure features that may be desirable for translation, biophysical size, and immunogenicity. Furthermore, their folding is robust to temperature, computer modeling method, choice of flanking untranslated regions, and changes in target protein sequence, as illustrated by rapid redesign of superfolder mRNAs for B.1.351, P.1 and B.1.1.7 variants of the prefusion-stabilized SARS-CoV-2 spike protein. Increases in in vitro mRNA half-life by at least two-fold appear immediately achievable.     

'In-line attack' conformational effect plays a modest role in an enzyme-catalyzed RNA cleavage: a free energy simulation study

Since the proposal of ‘in-line attack’ conformation as a possibly important intermediate in RNA cleavage, its structure has been captured in various protein and RNA enzymes; these structures strengthen the belief that this conformation plays an essential role in the catalysis of RNA cleavage. As generally discussed, this intermediate structure can be involved in energy barrier reduction in two possible ways, e.g. through either conformational effect or electrostatic effect. In order to quantitatively elucidate the contribution of conformational effect in this type of enzyme catalysis, free energy simulations were performed on the RNA structures both in a splicing endonuclease complex and in the aqueous solution. Our free energy simulation results revealed that the ‘in-line attack’ conformational effect plays a modest role in facilitating the reaction rate enhancement (~12-fold) compared with the overall 10¹²-fold rate increase. The close agreement between the present computational estimation and an experimental measurement on the spontaneous RNA cleavage in an in vitro evolved ATP aptamer motives us to realize that the conformation distribution of an enzyme substrate prior to rather than after its binding determines the upper bound of the rate enhancement ability through the conformational strategy.     

Fingerprinting the junctions of RNA structure by an open-paddlewheel diruthenium compound

RNA function is determined by its structural organization. The RNA structure consists of the combination of distinct secondary structure motifs connected by junctions that play an essential role in RNA folding. Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) probing is an established methodology to analyze the secondary structure of long RNA molecules in solution, which provides accurate data about unpaired nucleotides. However, the residues located at the junctions of RNA structures usually remain undetected. Here we report an RNA probing method based on the use of a novel open-paddlewheel diruthenium (OPW-Ru) compound [Ru₂Cl₂(µ-DPhF)₃(DMSO)] (DPhF = N,N′-diphenylformamidinate). This compound has four potential coordination sites in a singular disposition to establish covalent bonds with substrates. As a proof of concept, we have analyzed the reactivity of OPW-Ru toward RNA using two viral internal ribosome entry site (IRES) elements whose function depends on the structural organization of the molecule. Our study suggests that the compound OPW-Ru preferentially attacks at positions located one or two nucleotides away from junctions or bulges of the RNA structure. The OPW-Ru fingerprinting data differ from that obtained by other chemical reagents and provides new information about RNA structure features.     

The HIV-2 Rev-response element: determining secondary structure and defining folding intermediates

Interaction between the viral protein Rev and the RNA motifs known as Rev response elements (RREs) is required for transport of unspliced and partially spliced human immunodeficiency virus (HIV)-1 and HIV-2 RNAs from the nucleus to the cytoplasm during the later stages of virus replication. A more detailed understanding of these nucleoprotein complexes and the host factors with which they interact should accelerate the development of new antiviral drugs targeting cis-acting RNA regulatory signals. In this communication, the secondary structures of the HIV-2 RRE and two RNA folding precursors have been identified using the SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) chemical probing methodology together with a novel mathematical approach for determining the secondary structures of RNA conformers present in a mixture. A complementary chemical probing technique was also used to support these secondary structure models, to confirm that the RRE2 RNA undergoes a folding transition and to obtain information about the relative positioning of RRE2 substructures in three dimensions. Our analysis collectively suggests that the HIV-2 RRE undergoes two conformational transitions before assuming the energetically most favorable conformer. The 3D models for the HIV-2 RRE and folding intermediates are also presented, wherein the Rev-binding stem–loops (IIB and I) are located coaxially in the former, which is in agreement with previous models for HIV-1 Rev-RRE binding.     

Rational guide RNA engineering for small-molecule control of CRISPR/Cas9 and gene editing

It is important to control CRISPR/Cas9 when sufficient editing is obtained. In the current study, rational engineering of guide RNAs (gRNAs) is performed to develop small-molecule-responsive CRISPR/Cas9. For our purpose, the sequence of gRNAs are modified to introduce ligand binding sites based on the rational design of ligand–RNA pairs. Using short target sequences, we demonstrate that the engineered RNA provides an excellent scaffold for binding small molecule ligands. Although the ‘stem–loop 1’ variants of gRNA induced variable cleavage activity for different target sequences, all ‘stem–loop 3’ variants are well tolerated for CRISPR/Cas9. We further demonstrate that this specific ligand–RNA interaction can be utilized for functional control of CRISPR/Cas9 in vitro and in human cells. Moreover, chemogenetic control of gene editing in human cells transfected with all-in-one plasmids encoding Cas9 and designer gRNAs is demonstrated. The strategy may become a general approach for generating switchable RNA or DNA for controlling other biological processes.     

Functional specialization of monocot DCL3 and DCL5 proteins through the evolution of the PAZ domain

Monocot DICER-LIKE3 (DCL3) and DCL5 produce distinct 24-nt small interfering RNAs (siRNAs), heterochromatic siRNAs (hc-siRNAs) and phased secondary siRNAs (phasiRNAs), respectively. The former small RNAs are linked to silencing of transposable elements and heterochromatic repeats, and the latter to reproductive processes. It is assumed that these DCLs evolved from an ancient ‘eudicot-type’ DCL3 ancestor, which may have produced both types of siRNAs. However, how functional differentiation was achieved after gene duplication remains elusive. Here, we find that monocot DCL3 and DCL5 exhibit biochemically distinct preferences for 5′ phosphates and 3′ overhangs, consistent with the structural properties of their in vivo double-stranded RNA substrates. Importantly, these distinct substrate specificities are determined by the PAZ domains of DCL3 and DCL5, which have accumulated mutations during the course of evolution. These data explain the mechanism by which these DCLs cleave their cognate substrates from a fixed end, ensuring the production of functional siRNAs. Our study also indicates how plants have diversified and optimized RNA silencing mechanisms during evolution.     

The separation between the 5′-3′ ends in long RNA molecules is short and nearly constant

RNA molecules play different roles in coding, decoding and gene expression regulation. Such roles are often associated to the RNA secondary or tertiary structures. The folding dynamics lead to multiple secondary structures of long RNA molecules, since an RNA molecule might fold into multiple distinct native states. Despite an ensemble of different structures, it has been theoretically proposed that the separation between the 5′ and 3′ ends of long single-stranded RNA molecules (ssRNA) remains constant, independent of their base content and length. Here, we present the first experimental measurements of the end-to-end separation in long ssRNA molecules. To determine this separation, we use single molecule Fluorescence Resonance Energy Transfer of fluorescently end-labeled ssRNA molecules ranging from 500 to 5500 nucleotides in length, obtained from two viruses and a fungus. We found that the end-to-end separation is indeed short, within 5–9 nm. It is remarkable that the separation of the ends of all RNA molecules studied remains small and similar, despite the origin, length and differences in their secondary structure. This implies that the ssRNA molecules are ‘effectively circularized’ something that might be a general feature of RNAs, and could result in fine-tuning for translation and gene expression regulation.     

Targeted inhibition of the hepatitis C internal ribosomal entry site genomic RNA with oligonucleotide conjugates

Hepatitis C is a major public health concern, with an estimated 170 million people infected worldwide and an urgent need for new drug development. An attractive therapeutic approach is to prevent the ‘cap-independent’ translation initiation of the viral proteins by interfering with both the structure and function of the hepatitis C viral internal ribosomal entry site (HCV IRES). Towards this goal, we report the design, synthesis and purification of novel bi-functional molecules containing DNA or RNA antisenses attached to functional groups performing RNA hydrolysis. These 5′ or 3′-coupled conjugates bind the HCV IRES with affinity and specificity and elicit targeted hydrolysis of the viral genomic RNA after short (1 h) incubation at low (500 nM) concentration at 37°C in vitro. Additional secondary cleavage sites are induced and their mapping within the RNA structure indicates that functional domains IIIb-e are excised from the IRES that, based on cryo-EM studies, becomes incapable of binding the small ribosomal subunit and initiation factor 3 (eIF3). All these molecules inhibit, in a dose-dependent manner, the ‘IRES-dependent’ translation in vitro. The 5′-coupled imidazole conjugate reduces viral protein synthesis by half at a 300 nM concentration (IC50), corresponding to a 4-fold increase of activity when compared to the naked oligonucleotide. These new conjugates are now being tested for activity on infected hepatic cell lines.     

Secondary structure prediction and structure-specific sequence analysis of single-stranded DNA

DNA sequence analysis by oligonucleotide binding is often affected by interference with the secondary structure of the target DNA. Here we describe an approach that improves DNA secondary structure prediction by combining enzymatic probing of DNA by structure-specific 5′-nucleases with an energy minimization algorithm that utilizes the 5′-nuclease cleavage sites as constraints. The method can identify structural differences between two DNA molecules caused by minor sequence variations such as a single nucleotide mutation. It also demonstrates the existence of long-range interactions between DNA regions separated by >300 nt and the formation of multiple alternative structures by a 244 nt DNA molecule. The differences in the secondary structure of DNA molecules revealed by 5′-nuclease probing were used to design structure-specific probes for mutation discrimination that target the regions of structural, rather than sequence, differences. We also demonstrate the performance of structure-specific ‘bridge’ probes complementary to non-contiguous regions of the target molecule. The structure-specific probes do not require the high stringency binding conditions necessary for methods based on mismatch formation and permit mutation detection at temperatures from 4 to 37°C. Structure-specific sequence analysis is applied for mutation detection in the Mycobacterium tuberculosis katG gene and for genotyping of the hepatitis C virus.     

Direct structural analysis of modified RNA by fluorescent in-line probing

Chemical probing is a common method for the structural characterization of RNA. Typically, RNA is radioactively end-labelled, subjected to probing conditions, and the cleavage fragment pattern is analysed by gel electrophoresis. In recent years, many chemical modifications, like fluorophores, were introduced into RNA, but methods are lacking that detect the influence of the modification on the RNA structure with single-nucleotide resolution. Here, we first demonstrate that a 5'-terminal (32)P label can be replaced by a dye label for in-line probing of riboswitch RNAs. Next, we show that small, highly structured FRET-labelled Diels-Alderase ribozymes can be directly probed, using the internal or terminal FRET dyes as reporters. The probing patterns indeed reveal whether or not the attachment of the dyes influences the structure. The existence of two dye labels in typical FRET constructs is found to be beneficial, as 'duplexing' allows observation of the complete RNA on a single gel. Structural information can be derived from the probing gels by deconvolution of the superimposed band patterns. Finally, we use fluorescent in-line probing to experimentally validate the structural consequences of photocaging, unambiguously demonstrating the intentional destruction of selected elements of secondary or tertiary structure.     

‘RNA walk’ a novel approach to study RNA–RNA interactions between a small RNA and its target

In this study we describe a novel method to investigate the RNA–RNA interactions between a small RNA and its target that we termed ‘RNA walk’. The method is based on UV-induced AMT cross-linking in vivo followed by affinity selection of the hybrid molecules and mapping the intermolecular adducts by RT–PCR or real-time PCR. Domains carrying the cross-linked adducts fail to efficiently amplify by PCR compared with non-cross-linked domains. This method was calibrated and used to study the interaction between a special tRNA-like molecule (sRNA-85) that is part of the trypanosome signal recognition particle (SRP) complex and the ribosome. Four contact sites between sRNA-85 and rRNA were identified by ‘RNA walk’ and were further fine-mapped by primer extension. Two of the contact sites are expected; one contact site mimics the interaction of the mammalian Alu domain of SRP with the ribosome and the other contact sites include a canonical tRNA interaction. The two other cross-linked sites could not be predicted. We propose that ‘RNA walk, is a generic method to map target RNA small RNAs interactions in vivo.     

A novel SHAPE reagent enables the analysis of RNA structure in living cells with unprecedented accuracy

Due to the mounting evidence that RNA structure plays a critical role in regulating almost any physiological as well as pathological process, being able to accurately define the folding of RNA molecules within living cells has become a crucial need. We introduce here 2-aminopyridine-3-carboxylic acid imidazolide (2A3), as a general probe for the interrogation of RNA structures in vivo. 2A3 shows moderate improvements with respect to the state-of-the-art selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) reagent NAI on naked RNA under in vitro conditions, but it significantly outperforms NAI when probing RNA structure in vivo, particularly in bacteria, underlining its increased ability to permeate biological membranes. When used as a restraint to drive RNA structure prediction, data derived by SHAPE-MaP with 2A3 yields more accurate predictions than NAI-derived data. Due to its extreme efficiency and accuracy, we can anticipate that 2A3 will rapidly take over conventional SHAPE reagents for probing RNA structures both in vitro and in vivo.