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.     

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.     

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.     

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.     

Selective 2'-hydroxyl acylation analyzed by protection from exoribonuclease (RNase-detected SHAPE) for direct analysis of covalent adducts and of nucleotide flexibility in RNA

RNA SHAPE chemistry yields quantitative, single-nucleotide resolution structural information based on the reaction of the 2'-hydroxyl group of conformationally flexible nucleotides with electrophilic SHAPE reagents. However, SHAPE technology has been limited by the requirement that sites of RNA modification be detected by primer extension. Primer extension results in loss of information at both the 5' and 3' ends of an RNA and requires multiple experimental steps. Here we describe RNase-detected SHAPE that uses a processive, 3'→5' exoribonuclease, RNase R, to detect covalent adducts in 5'-end-labeled RNA in a one-tube experiment. RNase R degrades RNA but stops quantitatively three and four nucleotides 3' of a nucleotide containing a covalent adduct at the ribose 2'-hydroxyl or the pairing face of a nucleobase, respectively. We illustrate this technology by characterizing ligand-induced folding for the aptamer domain of the Escherichia coli thiamine pyrophosphate riboswitch RNA. RNase-detected SHAPE is a facile, two-day approach that can be used to analyze diverse covalent adducts in any RNA molecule, including short RNAs not amenable to analysis by primer extension and RNAs with functionally important structures at their 5' or 3' ends.     

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.     

RNA flexibility in the dimerization domain of a gamma retrovirus

Retroviruses are the causative agents of serious diseases, such as acquired immunodeficiency syndromes and several cancers, and are also useful gene therapy vectors. Retroviruses contain two sense-strand RNA genomes, which become linked at their 5' ends to form an RNA dimer. Understanding the molecular basis for dimerization may yield new approaches for controlling viral infectivity. Because this RNA domain is highly conserved within retrovirus groups, it has not been possible to define a consensus structure for the 5' dimerization domain by comparative sequence analysis. Here, we defined a 170-nucleotide minimal dimerization active sequence (MiDAS) for a representative gamma retrovirus, the Moloney murine sarcoma virus, by stringent competitive dimerization. We then analyzed the structure at every nucleotide in the MiDAS monomeric starting state with quantitative selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry. Notably, SHAPE analysis demonstrated that the RNA monomer contains an extensive flexible domain spanning 50 nucleotides. These findings support a structural model in which RNA flexibility directly facilitates retroviral genome dimerization by reducing the energetic cost of disrupting pre-existing base pairings in the monomer.     

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.     

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.     

Following the dynamics of changes in solvent accessibility of 16 S and 23 S rRNA during ribosomal subunit association using synchrotron-generated hydroxyl radicals

We have probed the structure and dynamics of ribosomal RNA in the Escherichia coli ribosome using equilibrium and time-resolved hydroxyl radical (OH) RNA footprinting to explore changes in the solvent-accessible surface of the rRNA with single-nucleotide resolution. The goal of these studies is to better understand the structural transitions that accompany association of the 30 S and 50 S subunits and to build a foundation for the quantitative analysis of ribosome structural dynamics during translation. Clear portraits of the subunit interface surfaces for 16 S and 23 S rRNA were obtained by constructing difference maps between the OH protection maps of the free subunits and that of the associated ribosome. In addition to inter-subunit contacts consistent with the crystal structure, additional OH protections are evident in regions at or near the subunit interface that reflect association-induced conformational changes. Comparison of these data with the comparable difference maps of the solvent-accessible surface of the rRNA calculated for the Thermus thermophilus X-ray crystal structures shows extensive agreement but also distinct differences. As a prelude to time-resolved OH footprinting studies, the reactivity profiles obtained using Fe(II)EDTA and X-ray generated OH were comprehensively compared. The reactivity patterns are similar except for a small number of nucleotides that have decreased reactivity to OH generated from Fe(II)EDTA compared to X-rays. These nucleotides are generally close to ribosomal proteins, which can quench diffusing radicals by virtue of side-chain oxidation. Synchrotron X-ray OH footprinting was used to monitor the kinetics of association of the 30 S and 50 S subunits. The rates individually measured for the inter-subunit contacts are comparable within experimental error. The application of this approach to the study of ribosome dynamics during the translation cycle is discussed.     

Molecular sensing by the aptamer domain of the FMN riboswitch: a general model for ligand binding by conformational selection

Understanding the nature of the free state of riboswitch aptamers is important for illuminating common themes in gene regulation by riboswitches. Prior evidence indicated the flavin mononucleotide (FMN)-binding riboswitch aptamer adopted a ‘bound-like’ structure in absence of FMN, suggesting only local conformational changes upon ligand binding. In the scope of pinpointing the general nature of such changes at the nucleotide level, we performed SHAPE mapping experiments using the aptamer domain of two phylogenetic variants, both in absence and in presence of FMN. We also solved the crystal structures of one of these domains both free (3.3 Å resolution) and bound to FMN (2.95 Å resolution). Our comparative study reveals that structural rearrangements occurring upon binding are restricted to a few of the joining regions that form the binding pocket in both RNAs. This type of binding event with minimal structural perturbations is reminiscent of binding events by conformational selection encountered in other riboswitches and various RNAs.     

SHAMS: Combining chemical modification of RNA with mass spectrometry to examine polypurine tract-containing RNA/DNA hybrids

Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) has gained popularity as a facile method of examining RNA structure both in vitro and in vivo, exploiting accessibility of the ribose 2′-OH to acylation by N-methylisatoic anhydride (NMIA) in unpaired or flexible configurations. Subsequent primer extension terminates at the site of chemical modification, and these products are fractionated by high-resolution gel electrophoresis. When applying SHAPE to investigate structural features associated with the wild-type and analog-substituted polypurine tract (PPT)–containing RNA/DNA hybrids, their size (20–25 base pairs) rendered primer extension impractical. As an alternative method of detection, we reasoned that chemical modification could be combined with tandem mass spectrometry, relying on the mass increment of RNA fragments containing the NMIA adduct (M_r = 133 Da). Using this approach, we demonstrate both specific modification of the HIV-1 PPT RNA primer and variations in its acylation pattern induced by replacing template nucleotides with a non-hydrogen-bonding thymine isostere. Our selective 2′-hydroxyl acylation analyzed by mass spectrometry strategy (SHAMS) should find utility when examining the structure of small RNA fragments or RNA/DNA hybrids where primer extension cannot be performed.     

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.     

Simplified RNA secondary structure mapping by automation of SHAPE data analysis

SHAPE (Selective 2'-hydroxyl acylation analysed by primer extension) technology has emerged as one of the leading methods of determining RNA secondary structure at the nucleotide level. A significant bottleneck in using SHAPE is the complex and time-consuming data processing that is required. We present here a modified data collection method and a series of algorithms, embodied in a program entitled Fast Analysis of SHAPE traces (FAST), which significantly reduces processing time. We have used this method to resolve the secondary structure of the first ~900 nt of the hepatitis C virus (HCV) genome, including the entire core gene. We have also demonstrated the ability of SHAPE/FAST to detect the binding of a small molecule inhibitor to the HCV internal ribosomal entry site (IRES). In conclusion, FAST allows for high-throughput data processing to match the current high-throughput generation of data possible with SHAPE, reducing the barrier to determining the structure of RNAs of interest.