In-cell RNA structure probing with SHAPE-MaP

RNA infrastructure profiling illuminates transcriptome structure in crowded spaces

RNAs fold into compact structures and undergo protein interactions in cells. These occluded environments can block reagents that probe the underlying RNAs. Probes that can analyze structure in crowded settings can shed light on RNA biology. Here, we employ 2'-OH-reactive probes that are small enough to access folded RNA structure underlying close molecular contacts within cells, providing considerably broader coverage for intracellular RNA structural analysis. The data are analyzed first with well-characterized human ribosomal RNAs and then applied transcriptome-wide to polyadenylated transcripts. The smallest probe acetylimidazole (AcIm) yields 80% greater structural coverage than larger conventional reagent NAIN3, providing enhanced structural information in hundreds of transcripts. The acetyl probe also provides superior signals for identifying m6A modification sites in transcripts, particularly in sites that are inaccessible to a standard probe. Our strategy enables profiling RNA infrastructure, enhancing analysis of transcriptome structure, modification, and intracellular interactions, especially in spatially crowded settings.

The role of structure in regulatory RNA elements

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e., RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Abolished frameshifting for predicted structure-stabilizing SARS-CoV-2 mutants: implications to alternative conformations and their statistical structural analyses

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases, including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition.

RNA structures within Venezuelan equine encephalitis virus E1 alter macrophage replication fitness and contribute to viral emergence

Venezuelan equine encephalitis virus (VEEV) is a mosquito-borne +ssRNA virus belonging to the Togaviridae. VEEV is found throughout Central and South America and is responsible for periodic epidemic/epizootic outbreaks of febrile and encephalitic disease in equines and humans. Endemic/enzootic VEEV is transmitted between Culex mosquitoes and sylvatic rodents, whereas epidemic/epizootic VEEV is transmitted between mosquitoes and equids, which serve as amplification hosts during outbreaks. Epizootic VEEV emergence has been shown to arise from mutation of enzootic VEEV strains. Specifically, epizootic VEEV has been shown to acquire amino acid mutations in the E2 viral glycoprotein that facilitate viral entry and equine amplification. However, the abundance of synonymous mutations which accumulate across the epizootic VEEV genome suggests that other viral determinants such as RNA secondary structure may also play a role in VEEV emergence. In this study we identify novel RNA structures in the E1 gene which specifically alter replication fitness of epizootic VEEV in macrophages but not other cell types. We show that SNPs are conserved within epizootic lineages and that RNA structures are conserved across different lineages. We also identified several novel RNA-binding proteins that are necessary for altered macrophage replication. These results suggest that emergence of VEEV in nature requires multiple mutations across the viral genome, some of which alter cell-type specific replication fitness in an RNA structure-dependent manner.

A novel reporter for helicase activity in translation uncovers DDX3X interactions

DDX3X regulates the translation of a subset of human transcripts containing complex 5' untranslated regions (5' UTRs). In this study, we developed the helicase activity reporter for translation (HART), which uses DDX3X-sensitive 5' UTRs to measure DDX3X-mediated translational activity in cells. To directly measure RNA structure in DDX3X-dependent mRNAs, we used SHAPE-MaP to determine the secondary structures present in DDX3X-sensitive 5' UTRs and then used HART to investigate how sequence alterations influence DDX3X sensitivity. Additionally, we identified residues 38-44 as potential mediators of DDX3X's interaction with the translational machinery. HART revealed that both DDX3X's association with the translational machinery and its helicase activity are required for its function in promoting the translation of DDX3X-sensitive 5' UTRs. These findings suggest DDX3X plays a crucial role in regulating translation through its interaction with the translational machinery during ribosome scanning and establish the HART reporter as a robust, lentivirally encoded, colorimetric measurement of DDX3X-dependent translation in cells.

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.

A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS

A key to understanding the roles of RNA in regulating gene expression is knowing their structures in vivo. One way to obtain this information is through probing the structures of RNA with chemicals. To probe RNA structure directly in cells, membrane-permeable reagents that modify the Watson-Crick (WC) face of unpaired nucleotides can be used. Although dimethyl sulfate (DMS) has led to substantial insight into RNA structure, it has limited nucleotide specificity in vivo, with WC face reactivity only at adenine (A) and cytosine (C) at neutral pH. The reagent 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) was recently shown to modify the WC face of guanine (G) and uracil (U). Although useful at lower concentrations in experiments that measure chemical modifications by reverse transcription (RT) stops, at higher concentrations necessary for detection by mutational profiling (MaP), EDC treatment leads to degradation of RNA. Here, we demonstrate EDC-stimulated degradation of RNA in Gram-negative and Gram-positive bacteria. In an attempt to overcome these limitations, we developed a new carbodiimide reagent, 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide methiodide (ETC), which we show specifically modifies unpaired Gs and Us in vivo without substantial degradation of RNA. We establish ETC as a probe for MaP and optimize the RT conditions and computational analysis in Escherichia coli Importantly, we demonstrate the utility of ETC as a probe for improving RNA structure prediction both alone and with DMS.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

A new mRNA structure prediction based approach to identifying improved signal peptides for bone morphogenetic protein 2

BACKGROUND

Signal peptide (SP) engineering has proven able to improve production of many proteins yet is a laborious process that still relies on trial and error. mRNA structure around the translational start site is important in translation initiation and has rarely been considered in this context, with recent improvements in in silico mRNA structure potentially rendering it a useful predictive tool for SP selection. Here we attempt to create a method to systematically screen candidate signal peptide sequences in silico based on both their nucleotide and amino acid sequences. Several recently released computational tools were used to predict signal peptide activity (SignalP), localization target (DeepLoc) and predicted mRNA structure (MXFold2). The method was tested with Bone Morphogenetic Protein 2 (BMP2), an osteogenic growth factor used clinically for bone regeneration. It was hoped more effective BMP2 SPs could improve BMP2-based gene therapies and reduce the cost of recombinant BMP2 production.

RESULTS

Amino acid sequence analysis indicated 2,611 SPs from the TGF-β superfamily were predicted to function when attached to BMP2. mRNA structure prediction indicated structures at the translational start site were likely highly variable. The five sequences with the most accessible translational start sites, a codon optimized BMP2 SP variant and the well-established hIL2 SP sequence were taken forward to in vitro testing. The top five candidates showed non-significant improvements in BMP2 secretion in HEK293T cells. All showed reductions in secretion versus the native sequence in C2C12 cells, with several showing large and significant decreases. None of the tested sequences were able to increase alkaline phosphatase activity above background in C2C12s. The codon optimized control sequence and hIL2 SP showed reasonable activity in HEK293T but very poor activity in C2C12.

CONCLUSIONS

These results support the use of peptide sequence based in silico tools for basic predictions around signal peptide activity in a synthetic biology context. However, mRNA structure prediction requires improvement before it can produce reliable predictions for this application. The poor activity of the codon optimized BMP2 SP variant in C2C12 emphasizes the importance of codon choice, mRNA structure, and cellular context for SP activity.

RNA structures within Venezuelan equine encephalitis virus E1 alter macrophage replication fitness and contribute to viral emergence

Venezuelan equine encephalitis virus (VEEV) is a mosquito-borne +ssRNA virus belonging to the Togaviridae. VEEV is found throughout Central and South America and is responsible for periodic epidemic/epizootic outbreaks of febrile and encephalitic disease in equines and humans. Endemic/enzootic VEEV is transmitted between Culex mosquitoes and sylvatic rodents, whereas epidemic/epizootic VEEV is transmitted between mosquitoes and equids, which serve as amplification hosts during outbreaks. Epizootic VEEV emergence has been shown to arise from mutation of enzootic VEEV strains. Specifically, epizootic VEEV has been shown to acquire amino acid mutations in the E2 viral glycoprotein that facilitate viral entry and equine amplification. However, the abundance of synonymous mutations which accumulate across the epizootic VEEV genome suggests that other viral determinants such as RNA secondary structure may also play a role in VEEV emergence. In this study we identify novel RNA structures in the E1 gene which specifically alter replication fitness of epizootic VEEV in macrophages but not other cell types. We show that SNPs are conserved within epizootic lineages and that RNA structures are conserved across different lineages. We also identified several novel RNA-binding proteins that are necessary for altered macrophage replication. These results suggest that emergence of VEEV in nature requires multiple mutations across the viral genome, some of which alter cell-type specific replication fitness in an RNA structure-dependent manner.

Abolished frameshifting for predicted structure-stabilizing SARS-CoV-2 mutants: Implications to alternative conformations and their statistical structural analyses

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition.

Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

Dissection of protein and RNA regions required for SPEN binding to XIST A-repeat RNA

XIST noncoding RNA promotes the initiation of X chromosome silencing by recruiting the protein SPEN to one X chromosome in female mammals. The SPEN protein is also called SHARP (SMRT and HDAC-associated repressor protein) and MINT (Msx-2 interacting nuclear target) in humans. SPEN recruits N-CoR2 and HDAC3 to initiate histone deacetylation on the X chromosome, leading to the formation of repressive chromatin marks and silencing gene expression. We dissected the contributions of different RNA and protein regions to the formation of a human XIST-SPEN complex in vitro and identified novel sequence and structure determinants that may contribute to X chromosome silencing initiation. Binding of SPEN to XIST RNA requires RRM 4 of the protein, in contrast to the requirement of RRM 3 and RRM 4 for specific binding to SRA RNA. Measurements of SPEN binding to full-length, dimeric, trimeric, or other truncated versions of the A-repeat region revealed that high-affinity binding of XIST to SPEN in vitro requires a minimum of four A-repeat segments. SPEN binding to XIST A-repeat RNA changes the accessibility of the RNA at specific nucleotide sequences, as indicated by changes in RNA reactivity through chemical structure probing. Based on computational modeling, we found that inter-repeat duplexes formed by multiple A-repeats can present an unpaired adenosine in the context of a double-stranded region of RNA. The presence of this specific combination of sequence and structural motifs correlates with high-affinity SPEN binding in vitro. These data provide new information on the molecular basis of the XIST and SPEN interaction.

Analysis of Telomerase RNA Structure in Physcomitrium patens Indicates Functionally Relevant Transitions Between OPEN and CLOSED Conformations

Telomerase RNA (TR) conformation determines its function as a template for telomere synthesis and as a scaffold for the assembly of the telomerase nucleoprotein complex. Experimental analyses of TR secondary structure using DMS-Map Seq and SHAPE-Map Seq techniques show its CLOSED conformation as the consensus structure where the template region cannot perform its function. Our data show that the apparent discrepancy between experimental results and predicted TR functional conformation, mostly ignored in published studies, can be explained using data analysis based on single-molecule structure prediction from individual sequencing reads by the recently established DaVinci method. This method results in several clusters of secondary structures reflecting the structural dynamics of TR, possibly related to its multiple functional states. Interestingly, the presumed active (OPEN) conformation of TR corresponds to a minor fraction of TR under in vivo conditions. Therefore, structural polymorphism and dynamic TR transitions between CLOSED and OPEN conformations may be involved in telomerase activity regulation as a switch that functions independently of total TR transcript levels.

Late consolidation of rRNA structure during co-transcriptional assembly in E. coli by time-resolved DMS footprinting

The production of new ribosomes requires proper folding of the rRNA and the addition of more than 50 ribosomal proteins. The structures of some assembly intermediates have been determined by cryo-electron microscopy, yet these structures do not provide information on the folding dynamics of the rRNA. To visualize the changes in rRNA structure during ribosome assembly in E. coli cells, transcripts were pulse-labeled with 4-thiouridine and the structure of newly made rRNA probed at various times by dimethyl sulfate modification and mutational profiling sequencing (4U-DMS-MaPseq). The in-cell DMS modification patterns revealed that many long-range rRNA tertiary interactions and protein binding sites through the 16S and 23S rRNA remain partially unfolded 1.5 min after transcription. By contrast, the active sites were continually shielded from DMS modification, suggesting that these critical regions are guarded by cellular factors throughout assembly. Later, bases near the peptidyl tRNA site exhibited specific rearrangements consistent with the binding and release of assembly factors. Time-dependent structure-probing in cells suggests that many tertiary interactions throughout the new ribosomal subunits remain mobile or unfolded until the late stages of subunit maturation.

RNA Folding, Mutation, and Detection

The structure of RNA molecules is absolutely critical to their functions in a biological system. RNA structure is dynamic and changes in response to cellular needs. Within the last few decades, there has been an increased interest in studying the structure of RNA molecules and how they change to support the needs of the cell in different conditions. Selective 2'-hydroxyl acylation-based mutational profiling using high-throughput sequencing is a powerful method to predict the secondary structure of RNA molecules both in vivo and in immunopurified samples. Selective 2'-hydroxyl acylation-based mutational profiling using high-throughput sequencing works by adding bulky groups onto accessible "flexible" bases in an RNA molecule that are not involved in any base-pairing or RNA-protein interactions. When the RNA is reverse transcribed into cDNA, the bulky groups are incorporated as base mutations, which can be compared to an unmodified control to identify the locations of flexible bases. The comparison of sequence data between modified and unmodified samples allows the computer software program (developed to generate reactivity profiles) to generate RNA secondary structure models. These models can be compared in a variety of conditions to determine how specific stimuli influence RNA secondary structures.

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.

A Molecular Orchestration of Plant Translation under Abiotic Stress

The complexities of translational strategies make this stage of implementing genetic information one of the most challenging to comprehend and, simultaneously, perhaps the most engaging. It is evident that this diverse range of strategies results not only from a long evolutionary history, but is also of paramount importance for refining gene expression and metabolic modulation. This notion is particularly accurate for organisms that predominantly exhibit biochemical and physiological reactions with a lack of behavioural ones. Plants are a group of organisms that exhibit such features. Addressing unfavourable environmental conditions plays a pivotal role in plant physiology. This is particularly evident with the changing conditions of global warming and the irrevocable loss or depletion of natural ecosystems. In conceptual terms, the plant response to abiotic stress comprises a set of elaborate and intricate strategies. This is influenced by a range of abiotic factors that cause stressful conditions, and molecular genetic mechanisms that fine-tune metabolic pathways allowing the plant organism to overcome non-standard and non-optimal conditions. This review aims to focus on the current state of the art in the field of translational regulation in plants under abiotic stress conditions. Different regulatory elements and patterns are being assessed chronologically. We deem it important to focus on significant high-performance techniques for studying the genetic information dynamics during the translation phase.

RNA Infrastructure Profiling Illuminates Transcriptome Structure in Crowded Spaces

RNAs can fold into compact three-dimensional structures, and most RNAs undergo protein interactions in the cell. These compact and occluded environments can block the ability of structure-probing agents to provide useful data about the folding and modification of the underlying RNA. The development of probes that can analyze structure in crowded settings, and differentiate the proximity of interactions, can shed new light on RNA biology. To this end, here we employ 2'-OH-reactive probes that are small enough to access folded RNA structure underlying many close molecular contacts within cells, providing considerably broader coverage for intracellular RNA structural analysis. We compare reverse transcriptase stops in RNA-Seq data from probes of small and standard size to assess RNA-protein proximity and evaluate solvent-exposed tunnels adjacent to RNA. The data are analyzed first with structurally characterized complexes (human 18S and 28S RNA), and then applied transcriptome-wide to polyadenylated transcripts in HEK293 cells. In our transcriptome profile, the smallest probe acetylimidazole (AcIm) yields 80% greater structural coverage than larger conventional reagent NAIN3, providing enhanced structural information in hundreds of transcripts. We further show that acetyl probes provide superior signals for identifying m6A modification sites in transcripts, and provide information regarding methylation sites that are inaccessible to a larger standard probe. RNA infrastructure profiling (RISP) enables enhanced analysis of transcriptome structure, modification, and interactions in living cells, especially in spatially crowded settings.

Determinants of DDX3X sensitivity uncovered using a helicase activity in translation reporter

DDX3X regulates the translation of a subset of human transcripts containing complex 5' untranslated regions (5' UTRs). In this study we developed the helicase activity reporter for translation (HART) which uses DDX3X-sensitive 5' UTRs to measure DDX3X mediated translational activity in cells. To dissect the structural underpinnings of DDX3X dependent translation, we first used SHAPE-MaP to determine the secondary structures present in DDX3X-sensitive 5' UTRs and then employed HART to investigate how their perturbation impacts DDX3X-sensitivity. Additionally, we identified residues 38-44 as potential mediators of DDX3X's interaction with the translational machinery. HART revealed that both DDX3X's association with the ribosome complex as well as its helicase activity are required for its function in promoting the translation of DDX3X-sensitive 5' UTRs. These findings suggest DDX3X plays a crucial role regulating translation through its interaction with the translational machinery during ribosome scanning, and establish the HART reporter as a robust, lentivirally encoded measurement of DDX3X-dependent translation in cells.

Pervasive downstream RNA hairpins dynamically dictate start-codon selection

Translational reprogramming allows organisms to adapt to changing conditions. Upstream start codons (uAUGs), which are prevalently present in mRNAs, have crucial roles in regulating translation by providing alternative translation start sites1-4. However, what determines this selective initiation of translation between conditions remains unclear. Here, by integrating transcriptome-wide translational and structural analyses during pattern-triggered immunity in Arabidopsis, we found that transcripts with immune-induced translation are enriched with upstream open reading frames (uORFs). Without infection, these uORFs are selectively translated owing to hairpins immediately downstream of uAUGs, presumably by slowing and engaging the scanning preinitiation complex. Modelling using deep learning provides unbiased support for these recognizable double-stranded RNA structures downstream of uAUGs (which we term uAUG-ds) being responsible for the selective translation of uAUGs, and allows the prediction and rational design of translating uAUG-ds. We found that uAUG-ds-mediated regulation can be generalized to human cells. Moreover, uAUG-ds-mediated start-codon selection is dynamically regulated. After immune challenge in plants, induced RNA helicases that are homologous to Ded1p in yeast and DDX3X in humans resolve these structures, allowing ribosomes to bypass uAUGs to translate downstream defence proteins. This study shows that mRNA structures dynamically regulate start-codon selection. The prevalence of this RNA structural feature and the conservation of RNA helicases across kingdoms suggest that mRNA structural remodelling is a general feature of translational reprogramming.

Identifying the structures of individual RNA isoforms inside cells

Nano-DMS-MaP focuses in on the structures of individual RNA isoforms, enabling direct examination of the structural diversity of different RNAs inside cells.

In vivo secondary structural analysis of Influenza A virus genomic RNA

Influenza A virus (IAV) is a respiratory virus that causes epidemics and pandemics. Knowledge of IAV RNA secondary structure in vivo is crucial for a better understanding of virus biology. Moreover, it is a fundament for the development of new RNA-targeting antivirals. Chemical RNA mapping using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) coupled with Mutational Profiling (MaP) allows for the thorough examination of secondary structures in low-abundance RNAs in their biological context. So far, the method has been used for analyzing the RNA secondary structures of several viruses including SARS-CoV-2 in virio and in cellulo. Here, we used SHAPE-MaP and dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) for genome-wide secondary structure analysis of viral RNA (vRNA) of the pandemic influenza A/California/04/2009 (H1N1) strain in both in virio and in cellulo environments. Experimental data allowed the prediction of the secondary structures of all eight vRNA segments in virio and, for the first time, the structures of vRNA5, 7, and 8 in cellulo. We conducted a comprehensive structural analysis of the proposed vRNA structures to reveal the motifs predicted with the highest accuracy. We also performed a base-pairs conservation analysis of the predicted vRNA structures and revealed many highly conserved vRNA motifs among the IAVs. The structural motifs presented herein are potential candidates for new IAV antiviral strategies.

Structural analysis of MALAT1 long noncoding RNA in cells and in evolution

Although not canonically polyadenylated, the long noncoding RNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) is stabilized by a highly conserved 76-nt triple helix structure on its 3' end. The entire MALAT1 transcript is over 8000 nt long in humans. The strongest structural conservation signal in MALAT1 (as measured by covariation of base pairs) is in the triple helix structure. Primary sequence analysis of covariation alone does not reveal the degree of structural conservation of the entire full-length transcript, however. Furthermore, RNA structure is often context dependent; RNA binding proteins that are differentially expressed in different cell types may alter structure. We investigate here the in-cell and cell-free structures of the full-length human and green monkey (Chlorocebus sabaeus) MALAT1 transcripts in multiple tissue-derived cell lines using SHAPE chemical probing. Our data reveal levels of uniform structural conservation in different cell lines, in cells and cell-free, and even between species, despite significant differences in primary sequence. The uniformity of the structural conservation across the entire transcript suggests that, despite seeing covariation signals only in the triple helix junction of the lncRNA, the rest of the transcript's structure is remarkably conserved, at least in primates and across multiple cell types and conditions.

The p53 endoplasmic reticulum stress-response pathway evolved in humans but not in mice via PERK-regulated p53 mRNA structures

Cellular stress conditions activate p53-dependent pathways to counteract the inflicted damage. To achieve the required functional diversity, p53 is subjected to numerous post-translational modifications and the expression of isoforms. Little is yet known how p53 has evolved to respond to different stress pathways. The p53 isoform p53/47 (p47 or ΔNp53) is linked to aging and neural degeneration and is expressed in human cells via an alternative cap-independent translation initiation from the 2nd in-frame AUG at codon 40 (+118) during endoplasmic reticulum (ER) stress. Despite an AUG codon in the same location, the mouse p53 mRNA does not express the corresponding isoform in either human or mouse-derived cells. High-throughput in-cell RNA structure probing shows that p47 expression is attributed to PERK kinase-dependent structural alterations in the human p53 mRNA, independently of eIF2α. These structural changes do not take place in murine p53 mRNA. Surprisingly, PERK response elements required for the p47 expression are located downstream of the 2nd AUG. The data show that the human p53 mRNA has evolved to respond to PERK-mediated regulation of mRNA structures in order to control p47 expression. The findings highlight how p53 mRNA co-evolved with the function of the encoded protein to specify p53-activities under different cellular conditions.

Strategies for targeting RNA with small molecule drugs

INTRODUCTION

Historically, therapeutic treatment of disease has been restricted to targeting proteins. Of the approximately 20,000 translated human proteins, approximately 1600 are associated with diseases. Strikingly, less than 15% of disease-associated proteins are predicted or known to be 'druggable.' While the concept and narrative of protein druggability continue to evolve with the development of novel technological and pharmacological advances, most of the human proteome remains undrugged. Recent genomic studies indicate that less than 2% of the human genome encodes for proteins, and while as much as 75% of the genome is transcribed, RNA has largely been ignored as a druggable target for therapeutic interventions.

AREAS COVERED

This review delineates the theory and techniques involved in the development of small molecule inhibitors of RNAs from brute force, high-throughput screening technologies to de novo molecular design using computational machine and deep learning. We will also highlight the potential pitfalls and limitations of targeting RNA with small molecules.

EXPERT OPINION

Although significant advances have recently been made in developing systems to identify small molecule inhibitors of RNAs, many challenges remain. Focusing on RNA structure and ligand binding sites may help bring drugging RNA in line with traditional protein drug targeting.

Chemical Probing of RNA Structure In Vivo Using SHAPE-MaP and DMS-MaP

The functional roles of RNAs are often regulated by their structure. Selective 2' hydroxyl acylation analyzed by primer extension (SHAPE) and dimethyl sulfate (DMS) base reactivity can be employed to investigate the flexibility of nucleotides and correlate it to the constraints imparted by base-pairing and/or protein-binding. In vivo, a multitude of proteins could bind an RNA molecule, regulating its structure and function. Hence, to obtain a more comprehensive view of the RNA structure-function relationship in vivo, it may be required to characterize both the RNA structure and the RNA-protein interaction network. In this chapter, we describe methods for characterizing the in vivo nucleotide flexibility of RNA in cells by SHAPE-MaP (SHAPE by Mutational Profiling) and DMS-MaP. In another chapter, we will discuss the characterization of RNA-protein interaction network by RNP-MaP.

Synthetic biology tools to promote the folding and function of RNA aptamers in mammalian cells

RNA aptamers are structured RNAs that can bind to diverse ligands, including proteins, metabolites, and other small molecules. RNA aptamers are widely used as in vitro affinity reagents. However, RNA aptamers have not been highly successful as bioactive intracellular molecules that can bind target molecules and influence cellular processes. We describe how poor RNA aptamer expression and especially poor RNA aptamer folding have limited the use of RNA aptamers in RNA synthetic biology applications. We discuss innovative new approaches that promote RNA aptamer folding in living cells and how these approaches have improved the function of aptamers in mammalian cells. These new approaches are making RNA aptamer-based synthetic biology and RNA aptamer therapeutic applications much more achievable.

Dissecting Functional Biological Interactions Using Modular RNA Nanoparticles

Nucleic acid nanoparticles (NANPs) are an exciting and innovative technology in the context of both basic and biomedical research. Made of DNA, RNA, or their chemical analogs, NANPs are programmed for carrying out specific functions within human cells. NANPs are at the forefront of preventing, detecting, and treating disease. Their nucleic acid composition lends them biocompatibility that provides their cargo with enhanced opportunity for coordinated delivery. Of course, the NANP system of targeting specific cells and tissues is not without its disadvantages. Accumulation of NANPs outside of the target tissue and the potential for off-target effects of NANP-mediated cargo delivery present challenges to research and medical professionals and these challenges must be effectively addressed to provide safe treatment to patients. Importantly, development of NANPs with regulated biological activities and immunorecognition becomes a promising route for developing versatile nucleic acid therapeutics. In a basic research context, NANPs can assist investigators in fine-tuning the structure-function relationship of final formulations and in this review, we explore the practical applications of NANPs in laboratory and clinical settings and discuss how we can use established nucleic acid research techniques to design effective NANPs.

Examples of Structural Motifs in Viral Genomes and Approaches for RNA Structure Characterization

The relationship between conserved structural motifs and their biological function in the virus replication cycle is the interest of many researchers around the world. RNA structure is closely related to RNA function. Therefore, technological progress in high-throughput approaches for RNA structure analysis and the development of new ones are very important. In this mini review, we discuss a few perspectives on the structural elements of viral genomes and some methods used for RNA structure prediction and characterization. Based on the recent literature, we describe several examples of studies concerning the viral genomes, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza A virus (IAV). Herein, we emphasize that a better understanding of viral genome architecture allows for the discovery of the structure-function relationship, and as a result, the discovery of new potential antiviral therapeutics.

Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.

ScanFold 2.0: a rapid approach for identifying potential structured RNA targets in genomes and transcriptomes

A major limiting factor in target discovery for both basic research and therapeutic intervention is the identification of structural and/or functional RNA elements in genomes and transcriptomes. This was the impetus for the original ScanFold algorithm, which provides maps of local RNA structural stability, evidence of sequence-ordered (potentially evolved) structure, and unique model structures comprised of recurring base pairs with the greatest structural bias. A key step in quantifying this propensity for ordered structure is the prediction of secondary structural stability for randomized sequences which, in the original implementation of ScanFold, is explicitly evaluated. This slow process has limited the rapid identification of ordered structures in large genomes/transcriptomes, which we seek to overcome in this current work introducing ScanFold 2.0. In this revised version of ScanFold, we no longer explicitly evaluate randomized sequence folding energy, but rather estimate it using a machine learning approach. For high randomization numbers, this can increase prediction speeds over 100-fold compared to ScanFold 1.0, allowing for the analysis of large sequences, as well as the use of additional folding algorithms that may be computationally expensive. In the testing of ScanFold 2.0, we re-evaluate the Zika, HIV, and SARS-CoV-2 genomes and compare both the consistency of results and the time of each run to ScanFold 1.0. We also re-evaluate the SARS-CoV-2 genome to assess the quality of ScanFold 2.0 predictions vs several biochemical structure probing datasets and compare the results to those of the original ScanFold program.

Cooperation and competition by RNA-binding proteins in cancer

Post-transcriptional regulation of gene expression plays a major role in determining the cellular proteome in health and disease. Post-transcriptional control mechanisms are disrupted in many cancers, contributing to multiple processes of tumorigenesis. RNA-binding proteins (RBPs), the main post-transcriptional regulators, often show altered expression and activity in cancer cells. Dysregulation of RBPs contributes to many cancer phenotypes, functioning in complex regulatory networks with other cellular players such as non-coding RNAs, signaling mediators and transcription factors to alter the expression of oncogenes and tumor suppressor genes. RBPs often function combinatorially, based on their binding to target sequences/structures on shared mRNA targets, to regulate the expression of cancer-related genes. This gives rise to cooperativity and competition between RBPs in mRNA binding and resultant functional outcomes in post-transcriptional processes such as mRNA splicing, stability, export and translation. Cooperation and competition is also observed in the case of interaction of RBPs and microRNAs with mRNA targets. RNA structural change is a common mechanism mediating the cooperative/competitive interplay between RBPs and between RBPs and microRNAs. RNA modifications, leading to changes in RNA structure, add a new dimension to cooperative/competitive binding of RBPs to mRNAs, further expanding the RBP regulatory landscape. Therefore, cooperative/competitive interplay between RBPs is a major determinant of the RBP interactome and post-transcriptional regulation of gene expression in cancer cells.

An overview of structural approaches to study therapeutic RNAs

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

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.

Rotavirus RNA chaperone mediates global transcriptome-wide increase in RNA backbone flexibility

Due to genome segmentation, rotaviruses must co-package eleven distinct genomic RNAs. The packaging is mediated by virus-encoded RNA chaperones, such as the rotavirus NSP2 protein. While the activities of distinct RNA chaperones are well studied on smaller RNAs, little is known about their global effect on the entire viral transcriptome. Here, we used Selective 2'-hydroxyl Acylation Analyzed by Primer Extension and Mutational Profiling (SHAPE-MaP) to examine the secondary structure of the rotavirus transcriptome in the presence of increasing amounts of NSP2. SHAPE-MaP data reveals that despite the well-documented helix-unwinding activity of NSP2 in vitro, its incubation with cognate rotavirus transcripts does not induce a significant change in the SHAPE reactivities. However, a quantitative analysis of mutation rates measured by mutational profiling reveals a global 5-fold rate increase in the presence of NSP2. We demonstrate that the normalization procedure used in deriving SHAPE reactivities from mutation rates can mask an important global effect of an RNA chaperone. Analysis of the mutation rates reveals a larger effect on stems rather than loops. Together, these data provide the first experimentally derived secondary structure model of the rotavirus transcriptome and reveal that NSP2 acts by globally increasing RNA backbone flexibility in a concentration-dependent manner.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

Current Status of Epitranscriptomic Marks Affecting lncRNA Structures and Functions

Long non-coding RNAs (lncRNAs) belong to a class of non-protein-coding RNAs with their lengths longer than 200 nucleotides. Most of the mammalian genome is transcribed as RNA, yet only a small percent of the transcribed RNA corresponds to exons of protein-coding genes. Thus, the number of lncRNAs is predicted to be several times higher than that of protein-coding genes. Because of sheer number of lncRNAs, it is often difficult to elucidate the functions of all lncRNAs, especially those arising from their relationship to their binding partners, such as DNA, RNA, and proteins. Due to their binding to other macromolecules, it has become evident that the structures of lncRNAs influence their functions. In this regard, the recent development of epitranscriptomics (the field of study to investigate RNA modifications) has become important to further elucidate the structures and functions of lncRNAs. In this review, the current status of lncRNA structures and functions influenced by epitranscriptomic marks is discussed.

A disease-linked lncRNA mutation in RNase MRP inhibits ribosome synthesis

RMRP encodes a non-coding RNA forming the core of the RNase MRP ribonucleoprotein complex. Mutations cause Cartilage Hair Hypoplasia (CHH), characterized by skeletal abnormalities and impaired T cell activation. Yeast RNase MRP cleaves a specific site in the pre-ribosomal RNA (pre-rRNA) during ribosome synthesis. CRISPR-mediated disruption of RMRP in human cells lines caused growth arrest, with pre-rRNA accumulation. Here, we analyzed disease-relevant primary cells, showing that mutations in RMRP impair mouse T cell activation and delay pre-rRNA processing. Patient-derived human fibroblasts with CHH-linked mutations showed similar pre-rRNA processing delay. Human cells engineered with the most common CHH mutation (70AG in RMRP) show specifically impaired pre-rRNA processing, resulting in reduced mature rRNA and a reduced ratio of cytosolic to mitochondrial ribosomes. Moreover, the 70AG mutation caused a reduction in intact RNase MRP complexes. Together, these results indicate that CHH is a ribosomopathy.

Determination of RNA Structure with In Vitro SHAPE Experiments

The structure of an RNA molecule is often critical for its correct functioning, post-transcriptional regulation, and/or translation. Predicting RNA secondary structures with in silico tools is relatively straightforward with the large array of software and webservers available. However, for long RNAs and RNA at high temperatures, in silico predictions are less accurate and require experimental validation. To this end, a variety of structural probing reagents are commonly used, both for in vitro and in vivo mapping of RNA structure. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments make use of a nonbase-specific modifying reagent, acylating all conformationally flexible (mainly single-stranded or unpaired) nucleotides and have been employed both for in vitro and in vivo modification of RNA. Here, we describe a protocol for an in vitro SHAPE experiment, starting from in vitro transcribed RNA using a T7 polymerase system. This RNA is folded and subsequently modified in vitro with the SHAPE-reagent N-methyl isatoic anhydride (NMIA). Primer extension employing a radioactive ³²P-labeled primer that binds the RNA downstream of the structure of interest will generate cDNA until the reverse transcriptase enzyme is halted by the introduced SHAPE modifications. Denaturing acrylamide gel electrophoresis of the pool of ³²P-labeled cDNAs and the corresponding sequencing ladders, followed by autoradiography, will expose these stops in reverse transcription (RT) and will therefore enable to identify single-stranded nucleotides in the RNA of interest. These RT stops and NMIA-modification efficiencies can be quantified with ImageJ software and can be used to validate or increase the accuracy of RNA secondary structure predictions.

An ultra low-input method for global RNA structure probing uncovers Regnase-1-mediated regulation in macrophages

To enable diverse functions and precise regulation, an RNA sequence often folds into complex yet distinct structures in different cellular states. Probing RNA in its native environment is essential to uncovering RNA structures of biological contexts. However, current methods generally require large amounts of input RNA and are challenging for physiologically relevant use. Here, we report smartSHAPE, a new RNA structure probing method that requires very low amounts of RNA input due to the largely reduced artefact of probing signals and increased efficiency of library construction. Using smartSHAPE, we showcased the profiling of the RNA structure landscape of mouse intestinal macrophages upon inflammation, and provided evidence that RNA conformational changes regulate immune responses. These results demonstrate that smartSHAPE can greatly expand the scope of RNA structure-based investigations in practical biological systems, and also provide a research paradigm for the study of post-transcriptional regulation.

Alternative poly-adenylation modulates α1-antitrypsin expression in chronic obstructive pulmonary disease

α1-anti-trypsin (A1AT), encoded by SERPINA1, is a neutrophil elastase inhibitor that controls the inflammatory response in the lung. Severe A1AT deficiency increases risk for Chronic Obstructive Pulmonary Disease (COPD), however, the role of A1AT in COPD in non-deficient individuals is not well known. We identify a 2.1-fold increase (p = 2.5x10-6) in the use of a distal poly-adenylation site in primary lung tissue RNA-seq in 82 COPD cases when compared to 64 controls and replicate this in an independent study of 376 COPD and 267 controls. This alternative polyadenylation event involves two sites, a proximal and distal site, 61 and 1683 nucleotides downstream of the A1AT stop codon. To characterize this event, we measured the distal ratio in human primary tissue short read RNA-seq data and corroborated our results with long read RNA-seq data. Integrating these results with 3' end RNA-seq and nanoluciferase reporter assay experiments we show that use of the distal site yields mRNA transcripts with over 50-fold decreased translation efficiency and A1AT expression. We identified seven RNA binding proteins using enhanced CrossLinking and ImmunoPrecipitation precipitation (eCLIP) with one or more binding sites in the SERPINA1 3' UTR. We combined these data with measurements of the distal ratio in shRNA knockdown experiments, nuclear and cytoplasmic fractionation, and chemical RNA structure probing. We identify Quaking Homolog (QKI) as a modulator of SERPINA1 mRNA translation and confirm the role of QKI in SERPINA1 translation with luciferase reporter assays. Analysis of single-cell RNA-seq showed differences in the distribution of the SERPINA1 distal ratio among hepatocytes, macrophages, αβ-Tcells and plasma cells in the liver. Alveolar Type 1,2, dendritic cells and macrophages also vary in their distal ratio in the lung. Our work reveals a complex post-transcriptional mechanism that regulates alternative polyadenylation and A1AT expression in COPD.

To Knot or Not to Knot: Multiple Conformations of the SARS-CoV-2 Frameshifting RNA Element

The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has a profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure, has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribosome during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention.

To knot or not to knot: Multiple conformations of the SARS-CoV-2 frameshifting RNA element

The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribosome during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention.

The RNA structurome in the asexual blood stages of malaria pathogen plasmodium falciparum

Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and non-coding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis that combined high-throughput nuclease mapping of RNA structures by duplex RNA-seq, SHAPE-directed RNA structure validation, immunoaffinity purification and characterization of antisense RNAs collectively measured differentially base-paired RNA regions throughout the parasite's asexual RBC cycle. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, approximately 71% of the genes with secondary structures are found to be protein coding mRNAs. The mapping pattern of these base-paired RNAs corresponded to all regions of mRNAs, including the 5' UTR, CDS and 3' UTR as well as the start and stop codons. Histone family genes which are known to form secondary structures in their mRNAs and transcripts from genes which are important for transcriptional and post-transcriptional control, such as the unique plant-like transcription factor family, ApiAP2, DNA-/RNA-binding protein, Alba3 and proteins important for RBC invasion and malaria cytoadherence also showed strong accumulation of duplex RNA reads in various asexual stages in P. falciparum. Intriguingly, our study determined stage-specific, dynamic relationships between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs.

In Vivo RNA Chemical Footprinting Analysis in Archaea

RNA structural conformation and dynamics govern the functional properties of all RNA/RNP. Accordingly, defining changes of RNA structure and dynamics in various conditions may provide detailed insight into how RNA structural properties regulate the function of RNA/RNP. Traditional chemical footprinting analysis using chemical modifiers allows to sample the dynamics and conformation landscape of diverse RNA/RNP. However, many chemical modifiers are limited in their capacity to provide unbiased information reflecting the in vivo RNA/RNP structural landscape. In the recent years, the development of selective-2'-hydroxyl acylation analyzed by primer extension (SHAPE) methodology that uses powerful new chemical modifiers has significantly improved in vitro and in vivo structural probing of secondary and tertiary interactions of diverse RNA species at the single nucleotide level.Although the original discovery of Archaea as an independent domain of life is intimately linked to the technological development of RNA analysis, our understanding of in vivo RNA structural conformation and dynamics in this domain of life remains scarce.This protocol describes the in vivo use of SHAPE chemistry in two evolutionary divergent model Archaea, Sulfolobus acidocaldarius and Haloferax volcanii.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

SHAPE Directed Discovery of New Functions in Large RNAs

RNA lies upstream of nearly all biology and functions as the central conduit of information exchange in all cells. RNA molecules encode information both in their primary sequences and in complex structures that form when an RNA folds back on itself. From the time of discovery of mRNA in the late 1950s until quite recently, we had only a rudimentary understanding of RNA structure across vast regions of most messenger and noncoding RNAs. This deficit is now rapidly being addressed, especially by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, mutational profiling (MaP), and closely related platform technologies that, collectively, create chemical microscopes for RNA. These technologies make it possible to interrogate RNA structure, quantitatively, at nucleotide resolution, and at large scales, for entire mRNAs, noncoding RNAs, and viral RNA genomes. By applying comprehensive structure probing to diverse problems, we and others are showing that control of biological function mediated by RNA structure is ubiquitous across prokaryotic and eukaryotic organisms.Work over the past decade using SHAPE-based analyses has clarified key principles. First, the method of RNA structure probing matters. SHAPE-MaP, with its direct and one-step readout that probes nearly every nucleotide by reaction at the 2'-hydroxyl, gives a more detailed and accurate readout than alternatives. Second, comprehensive chemical probing is essential. Focusing on fragments of large RNAs or using meta-gene or statistical analyses to compensate for sparse data sets misses critical features and often yields structure models with poor predictive power. Finally, every RNA has its own internal structural personality. There are myriad ways in which RNA structure modulates sequence accessibility, protein binding, translation, splice-site choice, phase separation, and other fundamental biological processes. In essentially every instance where we have applied rigorous and quantitative SHAPE technologies to study RNA structure-function interrelationships, new insights regarding biological regulatory mechanisms have emerged. RNA elements with more complex higher-order structures appear more likely to contain high-information-content clefts and pockets that bind small molecules, broadly informing a vigorous field of RNA-targeted drug discovery.The broad implications of this collective work are twofold. First, it is long past time to abandon depiction of large RNAs as simple noodle-like or gently flowing molecules. Instead, we need to emphasize that nearly all RNAs are punctuated with distinctive internal structures, a subset of which modulate function in profound ways. Second, structure probing should be an integral component of any effort that seeks to understand the functional nexuses and biological roles of large RNAs.