Influence of nucleotide identity on ribose 2′-hydroxyl reactivity in RNA

Hydroxyl-selective electrophiles, including N-methylisatoic anhydride (NMIA) and 1-methyl-7-nitroisatoic anhydride (1M7), are broadly useful for RNA structure analysis because they react preferentially with the ribose 2′-OH group at conformationally unconstrained or flexible nucleotides. Each nucleotide in an RNA has the potential to form an adduct with these reagents to yield a comprehensive, nucleotide-resolution, view of RNA structure. However, it is possible that factors other than local structure modulate reactivity. To evaluate the influence of base identity on the intrinsic reactivity of each nucleotide, we analyze NMIA and 1M7 reactivity using four distinct RNAs, under both native and denaturing conditions. We show that guanosine and adenosine residues have identical intrinsic 2′-hydroxyl reactivities at pH 8.0 and are 1.4 and 1.7 times more reactive than uridine and cytidine, respectively. These subtle, but statistically significant, differences do not impact the ability of selective 2′-hydroxyl acylation analyzed by primer extension-based (SHAPE) methods to establish an RNA secondary structure or monitor RNA folding in solution because base-specific influences are much smaller than the reactivity differences between paired and unpaired nucleotides.     

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.     

End-to-end crosstalk within the hepatitis C virus genome mediates the conformational switch of the 3′X-tail region

The hepatitis C virus (HCV) RNA genome contains multiple structurally conserved domains that make long-distance RNA–RNA contacts important in the establishment of viral infection. Microarray antisense oligonucelotide assays, improved dimethyl sulfate probing methods and 2′ acylation chemistry (selective 2’-hydroxyl acylation and primer extension, SHAPE) showed the folding of the genomic RNA 3′ end to be regulated by the internal ribosome entry site (IRES) element via direct RNA–RNA interactions. The essential cis-acting replicating element (CRE) and the 3′X-tail region adopted different 3D conformations in the presence and absence of the genomic RNA 5′ terminus. Further, the structural transition in the 3′X-tail from the replication-competent conformer (consisting of three stem-loops) to the dimerizable form (with two stem-loops), was found to depend on the presence of both the IRES and the CRE elements. Complex interplay between the IRES, the CRE and the 3′X-tail region would therefore appear to occur. The preservation of this RNA–RNA interacting network, and the maintenance of the proper balance between different contacts, may play a crucial role in the switch between different steps of the HCV cycle.     

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

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

Pyrvinium pamoate changes alternative splicing of the serotonin receptor 2C by influencing its RNA structure

The serotonin receptor 2C plays a central role in mood and appetite control. It undergoes pre-mRNA editing as well as alternative splicing. The RNA editing suggests that the pre-mRNA forms a stable secondary structure in vivo. To identify substances that promote alternative exons inclusion, we set up a high-throughput screen and identified pyrvinium pamoate as a drug-promoting exon inclusion without editing. Circular dichroism spectroscopy indicates that pyrvinium pamoate binds directly to the pre-mRNA and changes its structure. SHAPE (selective 2′-hydroxyl acylation analysed by primer extension) assays show that part of the regulated 5′-splice site forms intramolecular base pairs that are removed by this structural change, which likely allows splice site recognition and exon inclusion. Genome-wide analyses show that pyrvinium pamoate regulates >300 alternative exons that form secondary structures enriched in A–U base pairs. Our data demonstrate that alternative splicing of structured pre-mRNAs can be regulated by small molecules that directly bind to the RNA, which is reminiscent to an RNA riboswitch.     

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE): quantitative RNA structure analysis at single nucleotide resolution

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) interrogates local backbone flexibility in RNA at single-nucleotide resolution under diverse solution environments. Flexible RNA nucleotides preferentially sample local conformations that enhance the nucleophilic reactivity of 2'-hydroxyl groups toward electrophiles, such as N-methylisatoic anhydride (NMIA). Modified sites are detected as stops in an optimized primer extension reaction, followed by electrophoretic fragment separation. SHAPE chemistry scores local nucleotide flexibility at all four ribonucleotides in a single experiment and discriminates between base-paired versus unconstrained or flexible residues with a dynamic range of 20-fold or greater. Quantitative SHAPE reactivity information can be used to establish the secondary structure of an RNA, to improve the accuracy of structure prediction algorithms, to monitor structural differences between related RNAs or a single RNA in different states, and to detect ligand binding sites. SHAPE chemistry rarely needs significant optimization and requires two days to complete for an RNA of 100-200 nucleotides.     

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

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

QuShape: Rapid,accurate, and best-practices quantification of nucleic acid probing information,resolved by capillary electrophoresis

Chemical probing of RNA and DNA structure is a widely used and highly informative approach for examining nucleic acid structure and for evaluating interactions with protein and small-molecule ligands. Use of capillary electrophoresis to analyze chemical probing experiments yields hundreds of nucleotides of information per experiment and can be performed on automated instruments. Extraction of the information from capillary electrophoresis electropherograms is a computationally intensive multistep analytical process, and no current software provides rapid, automated, and accurate data analysis. To overcome this bottleneck, we developed a platform-independent, user-friendly software package, QuShape, that yields quantitatively accurate nucleotide reactivity information with minimal user supervision. QuShape incorporates newly developed algorithms for signal decay correction, alignment of time-varying signals within and across capillaries and relative to the RNA nucleotide sequence, and signal scaling across channels or experiments. An analysis-by-reference option enables multiple, related experiments to be fully analyzed in minutes. We illustrate the usefulness and robustness of QuShape by analysis of RNA SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) experiments.     

Identification of acylation products in SHAPE chemistry

SHAPE chemistry (selective 2′-hydroxyl acylation analyzed by primer extension) has been developed to specifically target flexible nucleotides (often unpaired nucleotides) independently to their purine or pyrimidine nature for RNA secondary structure determination. However, to the best of our knowledge, the structure of 2′-O-acylation products has never been confirmed by NMR or X-ray data. We have realized the acylation reactions between cNMP and NMIA under SHAPE chemistry conditions and identified the acylation products using standard NMR spectroscopy and LC–MS/MS experiments. For cAMP and cGMP, the major acylation product is the 2′-O-acylated compound (>99%). A trace amount of N-acylated cAMP has also been identified by LC–UV–MS². While for cCMP, the isolated acylation products are composed of 96% of 2′-O-acylated, 4% of N,O-diacylated, and trace amount of N-acylated compounds. In addition, the characterization of the major 2′-O-acylated compound by NMR showed slight differences in the conformation of the acylated sugar between the three cyclic nucleotides. This interesting result should be useful to explain some unexpected reactivity of the SHAPE chemistry.     

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.     

Visualizing the global secondary structure of a viral RNA genome with cryo-electron microscopy

The lifecycle, and therefore the virulence, of single-stranded (ss)-RNA viruses is regulated not only by their particular protein gene products, but also by the secondary and tertiary structure of their genomes. The secondary structure of the entire genomic RNA of satellite tobacco mosaic virus (STMV) was recently determined by selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE). The SHAPE analysis suggested a single highly extended secondary structure with much less branching than occurs in the ensemble of structures predicted by purely thermodynamic algorithms. Here we examine the solution-equilibrated STMV genome by direct visualization with cryo-electron microscopy (cryo-EM), using an RNA of similar length transcribed from the yeast genome as a control. The cryo-EM data reveal an ensemble of branching patterns that are collectively consistent with the SHAPE-derived secondary structure model. Thus, our results both elucidate the statistical nature of the secondary structure of large ss-RNAs and give visual support for modern RNA structure determination methods. Additionally, this work introduces cryo-EM as a means to distinguish between competing secondary structure models if the models differ significantly in terms of the number and/or length of branches. Furthermore, with the latest advances in cryo-EM technology, we suggest the possibility of developing methods that incorporate restraints from cryo-EM into the next generation of algorithms for the determination of RNA secondary and tertiary structures.     

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

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

Synergistic SHAPE/Single-Molecule Deconvolution of RNA Conformation under Physiological Conditions

Structural RNA domains are widely involved in the regulation of biological functions, such as gene expression, gene modification, and gene repair. Activity of these dynamic regions depends sensitively on the global fold of the RNA, in particular, on the binding affinity of individual conformations to effector molecules in solution. Consequently, both the 1) structure and 2) conformational dynamics of noncoding RNAs prove to be essential in understanding the coupling that results in biological function. Toward this end, we recently reported observation of three conformational states in the metal-induced folding pathway of the tRNA-like structure domain of Brome Mosaic Virus, via single-molecule fluorescence resonance energy transfer studies. We report herein selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE)-directed structure predictions as a function of metal ion concentrations ([Mⁿ⁺]) to confirm the three-state folding model, as well as test 2° structure models from the literature. Specifically, SHAPE reactivity data mapped onto literature models agrees well with the secondary structures observed at 0–10 mM [Mg²⁺], with only minor discrepancies in the E hairpin domain at low [Mg²⁺]. SHAPE probing and SHAPE-directed structure predictions further confirm the stepwise unfolding pathway previously observed in our single-molecule studies. Of special relevance, this means that reduction in the metal-ion concentration unfolds the 3′ pseudoknot interaction before unfolding the long-range stem interaction. This work highlights the synergistic power of combining 1) single-molecule Förster resonance energy transfer and 2) SHAPE-directed structure-probing studies for detailed analysis of multiple RNA conformational states. In particular, single-molecule guided deconvolution of the SHAPE reactivities permits 2° structure predictions of isolated RNA conformations, thereby substantially improving on traditional limitations associated with current structure prediction algorithms.     

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.     

Recruitment of the 40S Ribosome Subunit to the 3′-Untranslated Region (UTR) of a Viral mRNA,via the eIF4 Complex,Facilitates Cap-independent Translation

Barley yellow dwarf virus mRNA, which lacks both cap and poly(A) tail, has a translation element (3′-BTE) in its 3′-UTR essential for efficient translation initiation at the 5′-proximal AUG. This mechanism requires eukaryotic initiation factor 4G (eIF4G), subunit of heterodimer eIF4F (plant eIF4F lacks eIF4A), and 3′-BTE-5′-UTR interaction. Using fluorescence anisotropy, SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) analysis, and toeprinting, we found that (i) 40S subunits bind to BTE (K_d = 350 ± 30 nm), (ii) the helicase complex eIF4F-eIF4A-eIF4B-ATP increases 40S subunit binding (K_d = 120 ± 10 nm) to the conserved stem-loop I of the 3′-BTE by exposing more unpaired bases, and (iii) long distance base pairing transfers this complex to the 5′-end of the mRNA, where translation initiates. Although 3′-5′ interactions have been recognized as important in mRNA translation, barley yellow dwarf virus employs a novel mechanism utilizing the 3′-UTR as the primary site of ribosome recruitment.     

RecJ-like protein from Pyrococcus furiosus has 3′–5′ exonuclease activity on RNA: implications for proofreading of 3′-mismatched RNA primers in DNA replication

Replicative DNA polymerases require an RNA primer for leading and lagging strand DNA synthesis, and primase is responsible for the de novo synthesis of this RNA primer. However, the archaeal primase from Pyrococcus furiosus (Pfu) frequently incorporates mismatched nucleoside monophosphate, which stops RNA synthesis. Pfu DNA polymerase (PolB) cannot elongate the resulting 3′-mismatched RNA primer because it cannot remove the 3′-mismatched ribonucleotide. This study demonstrates the potential role of a RecJ-like protein from P. furiosus (PfRecJ) in proofreading 3′-mismatched ribonucleotides. PfRecJ hydrolyzes single-stranded RNA and the RNA strand of RNA/DNA hybrids in the 3′–5′ direction, and the kinetic parameters (K_m and K_cat) of PfRecJ during RNA strand digestion are consistent with a role in proofreading 3′-mismatched RNA primers. Replication protein A, the single-stranded DNA–binding protein, stimulates the removal of 3′-mismatched ribonucleotides of the RNA strand in RNA/DNA hybrids, and Pfu DNA polymerase can extend the 3′-mismatched RNA primer after the 3′-mismatched ribonucleotide is removed by PfRecJ. Finally, we reconstituted the primer-proofreading reaction of a 3′-mismatched ribonucleotide RNA/DNA hybrid using PfRecJ, replication protein A, Proliferating cell nuclear antigen (PCNA) and PolB. Given that PfRecJ is associated with the GINS complex, a central nexus in archaeal DNA replication fork, we speculate that PfRecJ proofreads the RNA primer in vivo.