The RNA transport element of the murine musD retrotransposon requires long-range intramolecular interactions for function

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

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.

RNA-mediated translation regulation in viral genomes: computational advances in the recognition of sequences and structures

RNA structures are widely distributed across all life forms. The global conformation of these structures is defined by a variety of constituent structural units such as helices, hairpin loops, kissing-loop motifs and pseudoknots, which often behave in a modular way. Their ubiquitous distribution is associated with a variety of functions in biological processes. The location of these structures in the genomes of RNA viruses is often coordinated with specific processes in the viral life cycle, where the presence of the structure acts as a checkpoint for deciding the eventual fate of the process. These structures have been found to adopt complex conformations and exert their effects by interacting with ribosomes, multiple host translation factors and small RNA molecules like miRNA. A number of such RNA structures have also been shown to regulate translation in viruses at the level of initiation, elongation or termination. The role of various computational studies in the preliminary identification of such sequences and/or structures and subsequent functional analysis has not been fully appreciated. This review aims to summarize the processes in which viral RNA structures have been found to play an active role in translational regulation, their global conformational features and the bioinformatics/computational tools available for the identification and prediction of these structures.

KRAB-zinc finger protein gene expansion in response to active retrotransposons in the murine lineage

The Krüppel-associated box zinc finger protein (KRAB-ZFP) family diversified in mammals. The majority of human KRAB-ZFPs bind transposable elements (TEs), however, since most TEs are inactive in humans it is unclear whether KRAB-ZFPs emerged to suppress TEs. We demonstrate that many recently emerged murine KRAB-ZFPs also bind to TEs, including the active ETn, IAP, and L1 families. Using a CRISPR/Cas9-based engineering approach, we genetically deleted five large clusters of KRAB-ZFPs and demonstrate that target TEs are de-repressed, unleashing TE-encoded enhancers. Homozygous knockout mice lacking one of two KRAB-ZFP gene clusters on chromosome 2 and chromosome 4 were nonetheless viable. In pedigrees of chromosome 4 cluster KRAB-ZFP mutants, we identified numerous novel ETn insertions with a modest increase in mutants. Our data strongly support the current model that recent waves of retrotransposon activity drove the expansion of KRAB-ZFP genes in mice and that many KRAB-ZFPs play a redundant role restricting TE activity.

New windows into retroviral RNA structures

Background

The multiple roles of both viral and cellular RNAs have become increasingly apparent in recent years, and techniques to model them have become significantly more powerful, enabling faster and more accurate visualization of RNA structures.

Main body

Techniques such as SHAPE (selective 2’OH acylation analysed by primer extension) have revolutionized the field, and have been used to examine RNAs belonging to many and diverse retroviruses. Secondary structure probing reagents such as these have been aided by the development of faster methods of analysis either via capillary or next-generation sequencing, allowing the analysis of entire genomes, and of retroviral RNA structures within virions. Techniques to model the three-dimensional structures of these large RNAs have also recently developed.

Conclusions

The flexibility of retroviral RNAs, both structural and functional, is clear from the results of these new experimental techniques. Retroviral RNA structures and structural changes control many stages of the lifecycle, and both the RNA structures themselves and their interactions with ligands are potential new drug targets. In addition, our growing understanding of retroviral RNA structures is aiding our knowledge of cellular RNA form and function.

Probing the Structures of Viral RNA Regulatory Elements with SHAPE and Related Methodologies

Viral RNAs were selected by evolution to possess maximum functionality in a minimal sequence. Depending on the classification of the virus and the type of RNA in question, viral RNAs must alternately be replicated, spliced, transcribed, transported from the nucleus into the cytoplasm, translated and/or packaged into nascent virions, and in most cases, provide the sequence and structural determinants to facilitate these processes. One consequence of this compact multifunctionality is that viral RNA structures can be exquisitely complex, often involving intermolecular interactions with RNA or protein, intramolecular interactions between sequence segments separated by several thousands of nucleotides, or specialized motifs such as pseudoknots or kissing loops. The fluidity of viral RNA structure can also present a challenge when attempting to characterize it, as genomic RNAs especially are likely to sample numerous conformations at various stages of the virus life cycle. Here we review advances in chemoenzymatic structure probing that have made it possible to address such challenges with respect to cis-acting elements, full-length viral genomes and long non-coding RNAs that play a major role in regulating viral gene expression.

Novel Biochemical Tools for Probing HIV RNA Structure

Functional analysis of viral RNA requires knowledge of secondary structure arrangements and tertiary base interactions. Thus, high-throughput and comprehensive methods for assessing RNA structure are highly desirable. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) has proven highly useful for modeling the secondary structures of HIV and other retroviral RNAs in recent years. This technology is not without its limitations however, as SHAPE data can be severely compromised when the RNA under study is structurally heterogeneous. In addition, the method reveals little information regarding the three-dimensional (3D) organization of an RNA. This chapter outlines four detailed SHAPE-related methodologies that circumvent these limitations. "Ensemble" and "in-gel" variations of SHAPE permit structural analysis of individual conformers within structurally heterogeneous mixtures of RNA, while probing strategies that utilize "through-space" cleavage reagents such as methidiumpropyl-EDTA (MPE) and peptides appended with an ATCUN (amino terminal copper/nickel binding motif) can provide insight into 3D organization. Combinational application of these techniques provides a formidable arsenal for exploring the structures of HIV RNAs and associated nucleoprotein complexes.

A small stem-loop structure of the Ebola virus trailer is essential for replication and interacts with heat-shock protein A8

Ebola virus (EBOV) is a single-stranded negative-sense RNA virus belonging to the Filoviridae family. The leader and trailer non-coding regions of the EBOV genome likely regulate its transcription, replication, and progeny genome packaging. We investigated the cis-acting RNA signals involved in RNA–RNA and RNA–protein interactions that regulate replication of eGFP-encoding EBOV minigenomic RNA and identified heat shock cognate protein family A (HSC70) member 8 (HSPA8) as an EBOV trailer-interacting host protein. Mutational analysis of the trailer HSPA8 binding motif revealed that this interaction is essential for EBOV minigenome replication. Selective 2′-hydroxyl acylation analyzed by primer extension analysis of the secondary structure of the EBOV minigenomic RNA indicates formation of a small stem-loop composed of the HSPA8 motif, a 3′ stem-loop (nucleotides 1868–1890) that is similar to a previously identified structure in the replicative intermediate (RI) RNA and a panhandle domain involving a trailer-to-leader interaction. Results of minigenome assays and an EBOV reverse genetic system rescue support a role for both the panhandle domain and HSPA8 motif 1 in virus replication.

Characterizing 3D RNA structure by single molecule FRET

The importance of elucidating the three dimensional structures of RNA molecules is becoming increasingly clear. However, traditional protein structural techniques such as NMR and X-ray crystallography have several important drawbacks when probing long RNA molecules. Single molecule Förster resonance energy transfer (smFRET) has emerged as a useful alternative as it allows native sequences to be probed in physiological conditions and allows multiple conformations to be probed simultaneously. This review serves to describe the method of generating a three dimensional RNA structure from smFRET data from the biochemical probing of the secondary structure to the computational refinement of the final model.

Automated RNA 3D Structure Prediction with RNAComposer

RNAs adopt specific structures to perform their activities and these are critical to virtually all RNA-mediated processes. Because of difficulties in experimentally assessing structures of large RNAs using NMR, X-ray crystallography, or cryo-microscopy, there is currently great demand for new high-resolution 3D structure prediction methods. Recently we reported on RNAComposer, a knowledge-based method for the fully automated RNA 3D structure prediction from a user-defined secondary structure. RNAComposer method is especially suited for structural biology users. Since our initial report in 2012, both servers, freely available at http://rnacomposer.ibch.poznan.pl and http://rnacomposer.cs.put.poznan.pl have been often visited. Therefore this chapter provides guidance for using RNAComposer and discusses points that should be considered when predicting 3D RNA structure. An application example presents current scope and limitations of RNAComposer.

The "complex" RNA post-transcriptional element of a "simple" retrovirus

Replication of retroviruses and transposition of endogenous retroelements exploits a unique mechanism of post-transcriptional regulation as a means of exporting their incompletely-spliced mRNAs (which serve as both the genomic RNA and the template for protein synthesis). Following discovery of the Rev response element (RRE) that mediates nucleocytoplasmic export of the full-length and singly-spliced human immunodeficiency virus type 1 (HIV-1) genome, equivalent cis-acting regulatory elements have been characterized for both complex and simple retroviruses and retroelements, together with the obligate viral and host proteins with which they interact. The exception to this is the gammaretrovirus family of simple retroviruses, exemplified by reticuloendotheliosis virus (REV), murine leukemia virus (MLV) and xenotropic MLV-related retrovirus (XMRV). In this commentary, we discuss our recent data that reported structural and functional data on the MLV/XMRV post-transcriptional regulatory element (designated the PTE). The PTE was characterized by a highly-structured region of multiple stem-loops (SL1 – SL7) overlapping the pro and 5′ portion of the pol open reading frames, comprising a bipartite export signal whose structures are separated by ∼1400 nt. In addition, structural probing suggested that SL3 nucleotides were involved in pseudoknot formation. These data, when compared with RNA transport elements of complex retroviruses (HIV) and simple murine retrotransposons (musD), collectively present an emerging picture that long-range tertiary interactions are critical mediators of their biological function.

Evolutionary conservation of a molecular machinery for export and expression of mRNAs with retained introns

Intron retention is one of the least studied forms of alternative splicing. Through the use of retrovirus and other model systems, it was established many years ago that mRNAs with retained introns are subject to restriction both at the level of nucleocytoplasmic export and cytoplasmic expression. It was also demonstrated that specific cis-acting elements in the mRNA could serve to bypass these restrictions. Here we show that one of these elements, the constitutive transport element (CTE), first identified in the retrovirus MPMV and subsequently in the human NXF1 gene, is a highly conserved element. Using GERP analysis, CTEs with strong primary sequence homology, predicted to display identical secondary structure, were identified in NXF genes from >30 mammalian species. CTEs were also identified in the predicted NXF1 genes of zebrafish and coelacanths. The CTE from the zebrafish NXF1 was shown to function efficiently to achieve expression of mRNA with a retained intron in human cells in conjunction with zebrafish Nxf1 and cofactor Nxt proteins. This demonstrates that all essential functional components for expression of mRNA with retained introns have been conserved from fish to man.

Insights into secondary and tertiary interactions of dengue virus RNA by SHAPE

Dengue virus (DENV) is a single-stranded positive-sense RNA virus belonging to the Flaviviridae family. The DENV RNA genome contains multiple cis-acting elements that continue to unravel their essential role in managing viral molecular processes. Attempts have been made to predict the secondary structure of DENV RNA using a variety of computational tools. Nevertheless, a greater degree of accuracy is achieved when these methods are complemented with structure probing experimentation. This chapter outlines detailed methodology for the structural study of DENV subgenomic minigenome RNA by applying high-throughput selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). High-throughput SHAPE combines a novel chemical probing technology with reverse transcription, capillary electrophoresis, and secondary structure prediction software to rapidly and reproducibly determine the structure of RNAs from several hundred to several thousand nucleotides at single-nucleotide resolution. This methodology investigates local structure for all positions in a sequence-independent manner and as such it is particularly useful in predicting RNA secondary and tertiary interactions.

Gammaretrovirus mRNA expression is mediated by a novel, bipartite post-transcriptional regulatory element

Post-transcriptional regulatory mechanisms of several complex and simple retroviruses and retroelements have been elucidated, with the exception of the gammaretrovirus family. We found that, similar to the other retroviruses, gag gene expression of MuLV and XMRV depends on post-transcriptional regulation mediated via an RNA sequence overlapping the pro-pol open reading frame, termed the Post-Transcriptional Element (PTE). PTE function can be replaced by heterologous RNA export elements, e.g. CTE of simian type D retroviruses. Alternatively, Gag particle production is achieved using an RNA/codon optimized gag gene. PTE function is transferable and can replace HIV Rev-RRE-regulated expression of HIV gag. Analysis of PTE by SHAPE revealed a highly structured RNA comprising seven stem-loop structures, with the 5′ and 3′ stem-loops forming an essential bipartite signal. MuLV and XMRV PTE share 98% identity and have highly similar RNA structures, with changes mostly located to single-stranded regions. PTE identification strongly suggests that all retroviruses and retroelements share common strategies of post-transcriptional gene regulation to produce Gag. Expression depends on complex RNA structures embedded within retroviral mRNA, in coding regions or the 3′ untranslated region. These specific structures serve as recognition signals for either cellular or viral proteins.

Tackling structures of long noncoding RNAs

RNAs are important catalytic machines and regulators at every level of gene expression. A new class of RNAs has emerged called long non-coding RNAs, providing new insights into evolution, development and disease. Long non-coding RNAs (lncRNAs) predominantly found in higher eukaryotes, have been implicated in the regulation of transcription factors, chromatin-remodeling, hormone receptors and many other processes. The structural versatility of RNA allows it to perform various functions, ranging from precise protein recognition to catalysis and metabolite sensing. While major housekeeping RNA molecules have long been the focus of structural studies, lncRNAs remain the least characterized class, both structurally and functionally. Here, we review common methodologies used to tackle RNA structure, emphasizing their potential application to lncRNAs. When considering the complexity of lncRNAs and lack of knowledge of their structure, chemical probing appears to be an indispensable tool, with few restrictions in terms of size, quantity and heterogeneity of the RNA molecule. Probing is not constrained to in vitro analysis and can be adapted to high-throughput sequencing platforms. Significant efforts have been applied to develop new in vivo chemical probing reagents, new library construction protocols for sequencing platforms and improved RNA prediction software based on the experimental evidence.

In-gel probing of individual RNA conformers within a mixed population reveals a dimerization structural switch in the HIV-1 leader

Definitive secondary structural mapping of RNAs in vitro can be complicated by the presence of more than one structural conformer or multimerization of some of the molecules. Until now, probing a single structure of conformationally flexible RNA molecules has typically relied on introducing stabilizing mutations or adjusting buffer conditions or RNA concentration. Here, we present an in-gel SHAPE (selective 2'OH acylation analysed by primer extension) approach, where a mixed structural population of RNA molecules is separated by non-denaturing gel electrophoresis and the conformers are individually probed within the gel matrix. Validation of the technique using a well-characterized RNA stem-loop structure, the HIV-1 trans-activation response element, showed that authentic structure was maintained and that the method was accurate and highly reproducible. To further demonstrate the utility of in-gel SHAPE, we separated and examined monomeric and dimeric species of the HIV-1 packaging signal RNA. Extensive differences in acylation sensitivity were seen between monomer and dimer. The results support a recently proposed structural switch model of RNA genomic dimerization and packaging, and demonstrate the discriminatory power of in-gel SHAPE.

RNA secondary structure prediction using high-throughput SHAPE

Understanding the function of RNA involved in biological processes requires a thorough knowledge of RNA structure. Toward this end, the methodology dubbed "high-throughput selective 2' hydroxyl acylation analyzed by primer extension", or SHAPE, allows prediction of RNA secondary structure with single nucleotide resolution. This approach utilizes chemical probing agents that preferentially acylate single stranded or flexible regions of RNA in aqueous solution. Sites of chemical modification are detected by reverse transcription of the modified RNA, and the products of this reaction are fractionated by automated capillary electrophoresis (CE). Since reverse transcriptase pauses at those RNA nucleotides modified by the SHAPE reagents, the resulting cDNA library indirectly maps those ribonucleotides that are single stranded in the context of the folded RNA. Using ShapeFinder software, the electropherograms produced by automated CE are processed and converted into nucleotide reactivity tables that are themselves converted into pseudo-energy constraints used in the RNAStructure (v5.3) prediction algorithm. The two-dimensional RNA structures obtained by combining SHAPE probing with in silico RNA secondary structure prediction have been found to be far more accurate than structures obtained using either method alone.

Secondary structure of bacteriophage T4 gene 60 mRNA: implications for translational bypassing

Translational bypassing is a unique phenomenon of bacteriophage T4 gene 60 mRNA wherein the bacterial ribosome produces a single polypeptide chain from a discontinuous open reading frame (ORF). Upon reaching the 50-nucleotide untranslated region, or coding gap, the ribosome either dissociates or bypasses the interruption to continue translating the remainder of the ORF, generating a subunit of a type II DNA topoisomerase. Mutational and computational analyses have suggested that a compact structure, including a stable hairpin, forms in the coding gap to induce bypassing, yet direct evidence is lacking. Here we have probed the secondary structure of gene 60 mRNA with both Tb³⁺ ions and the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reagent 1M7 under conditions where bypassing is observed. The resulting experimentally informed secondary structure models strongly support the presence of the predicted coding gap hairpin and highlight the benefits of using Tb³⁺ as a second, complementary probing reagent. Contrary to several previously proposed models, however, the rest of the coding gap is highly reactive with both probing reagents, suggesting that it forms only a short stem-loop. Mutational analyses coupled with functional assays reveal that two possible base-pairings of the coding gap with other regions of the mRNA are not required for bypassing. Such structural autonomy of the coding gap is consistent with its recently discovered role as a mobile genetic element inserted into gene 60 mRNA to inhibit cleavage by homing endonuclease MobA.

Postage for the messenger: designating routes for nuclear mRNA export

Transcription of mRNA occurs in the nucleus, making the translocation of mRNA across the nuclear envelope (NE) boundary a critical determinant of proper gene expression and cell survival. A major mRNA export route occurs via the NXF1-dependent pathway through the nuclear pore complexes (NPCs) embedded in the NE. However, recent findings have discovered new evidence supporting the existence of multiple mechanisms for crossing the NE, including both NPC-mediated and NE budding-mediated pathways. An analysis of the trans-acting factors and cis components that define these pathways reveals shared elements as well as mechanistic differences. We review here the current understanding of the mechanisms that characterize each pathway and highlight the determinants that influence mRNA transport fate.

Rise of the RNA machines: exploring the structure of long non-coding RNAs

Novel, profound and unexpected roles of long non-coding RNAs (lncRNAs) are emerging in critical aspects of gene regulation. Thousands of lncRNAs have been recently discovered in a wide range of mammalian systems, related to development, epigenetics, cancer, brain function and hereditary disease. The structural biology of these lncRNAs presents a brave new RNA world, which may contain a diverse zoo of new architectures and mechanisms. While structural studies of lncRNAs are in their infancy, we describe existing structural data for lncRNAs, as well as crystallographic studies of other RNA machines and their implications for lncRNAs. We also discuss the importance of dynamics in RNA machine mechanism. Determining commonalities between lncRNA systems will help elucidate the evolution and mechanistic role of lncRNAs in disease, creating a structural framework necessary to pursue lncRNA-based therapeutics.

Probing Retroviral and Retrotransposon Genome Structures: The "SHAPE" of Things to Come

Understanding the nuances of RNA structure as they pertain to biological function remains a formidable challenge for retrovirus research and development of RNA-based therapeutics, an area of particular importance with respect to combating HIV infection. Although a variety of chemical and enzymatic RNA probing techniques have been successfully employed for more than 30 years, they primarily interrogate small (100–500 nt) RNAs that have been removed from their biological context, potentially eliminating long-range tertiary interactions (such as kissing loops and pseudoknots) that may play a critical regulatory role. Selective 2′ hydroxyl acylation analyzed by primer extension (SHAPE), pioneered recently by Merino and colleagues, represents a facile, user-friendly technology capable of interrogating RNA structure with a single reagent and, combined with automated capillary electrophoresis, can analyze an entire 10,000-nucleotide RNA genome in a matter of weeks. Despite these obvious advantages, SHAPE essentially provides a nucleotide “connectivity map,” conversion of which into a 3-D structure requires a variety of complementary approaches. This paper summarizes contributions from SHAPE towards our understanding of the structure of retroviral genomes, modifications to which technology that have been developed to address some of its limitations, and future challenges.

Automated 3D structure composition for large RNAs

Understanding the numerous functions that RNAs play in living cells depends critically on knowledge of their three-dimensional structure. Due to the difficulties in experimentally assessing structures of large RNAs, there is currently great demand for new high-resolution structure prediction methods. We present the novel method for the fully automated prediction of RNA 3D structures from a user-defined secondary structure. The concept is founded on the machine translation system. The translation engine operates on the RNA FRABASE database tailored to the dictionary relating the RNA secondary structure and tertiary structure elements. The translation algorithm is very fast. Initial 3D structure is composed in a range of seconds on a single processor. The method assures the prediction of large RNA 3D structures of high quality. Our approach needs neither structural templates nor RNA sequence alignment, required for comparative methods. This enables the building of unresolved yet native and artificial RNA structures. The method is implemented in a publicly available, user-friendly server RNAComposer. It works in an interactive mode and a batch mode. The batch mode is designed for large-scale modelling and accepts atomic distance restraints. Presently, the server is set to build RNA structures of up to 500 residues.

The mechanisms of RNA SHAPE chemistry

The biological functions of RNA are ultimately governed by the local environment at each nucleotide. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry is a powerful approach for measuring nucleotide structure and dynamics in diverse biological environments. SHAPE reagents acylate the 2'-hydroxyl group at flexible nucleotides because unconstrained nucleotides preferentially sample rare conformations that enhance the nucleophilicity of the 2'-hydroxyl. The critical corollary is that some constrained nucleotides must be poised for efficient reaction at the 2'-hydroxyl group. To identify such nucleotides, we performed SHAPE on intact crystals of the Escherichia coli ribosome, monitored the reactivity of 1490 nucleotides in 16S rRNA, and examined those nucleotides that were hyper-reactive toward SHAPE and had well-defined crystallographic conformations. Analysis of these conformations revealed that 2'-hydroxyl reactivity is broadly facilitated by general base catalysis involving multiple RNA functional groups and by two specific orientations of the bridging 3'-phosphate group. Nucleotide analog studies confirmed the contributions of these mechanisms to SHAPE reactivity. These results provide a strong mechanistic explanation for the relationship between SHAPE reactivity and local RNA dynamics and will facilitate interpretation of SHAPE information in the many technologies that make use of this chemistry.

Correlating SHAPE signatures with three-dimensional RNA structures

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a facile technique for quantitative analysis of RNA secondary structure. In general, low SHAPE signal values indicate Watson-Crick base-pairing, and high values indicate positions that are single-stranded within the RNA structure. However, the relationship of SHAPE signals to structural properties such as non-Watson-Crick base-pairing or stacking has thus far not been thoroughly investigated. Here, we present results of SHAPE experiments performed on several RNAs with published three-dimensional structures. This strategy allows us to analyze the results in terms of correlations between chemical reactivities and structural properties of the respective nucleotide, such as different types of base-pairing, stacking, and phosphate-backbone interactions. We find that the RNA SHAPE signal is strongly correlated with cis-Watson-Crick/Watson-Crick base-pairing and is to a remarkable degree not dependent on other structural properties with the exception of stacking. We subsequently generated probabilistic models that estimate the likelihood that a residue with a given SHAPE score participates in base-pairing. We show that several models that take SHAPE scores of adjacent residues into account perform better in predicting base-pairing compared with individual SHAPE scores. This underscores the context sensitivity of SHAPE and provides a framework for an improved interpretation of the response of RNA to chemical modification.