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.     

Slow conformational dynamics at C2'-endo nucleotides in RNA

RNA molecules undergo local conformational dynamics on timescales spanning picoseconds to minutes. Slower local motions have the greater potential to govern RNA folding, ligand recognition, and ribonucleoprotein assembly reactions but are difficult to detect in large RNAs with complex structures. RNA SHAPE chemistry employs acylation of the ribose 2'-hydroxyl position to measure local nucleotide flexibility in RNA and is well-characterized by a mechanism in which each nucleotide samples unreactive (closed) and reactive (open) states. We monitor RNA conformational dynamics over distinct time domains by varying the electrophilicity of the acylating reagent. Select C2'-endo nucleotides are nonreactive toward fast reagents but reactive toward slower SHAPE reagents in both model RNAs and in a large RNA with a tertiary fold. We conclude, first, that the C2'-endo conformation by itself does not govern SHAPE reactivity. However, some C2'-endo nucleotides undergo extraordinarily slow conformational changes, on the order of 10(-4) s(-1). Due to their distinctive local dynamics, C2'-endo nucleotides have the potential to function as rate-determining molecular switches and are likely to play central, currently unexplored, roles in RNA folding and function.     

Time-resolved RNA SHAPE chemistry

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry yields quantitative RNA secondary and tertiary structure information at single nucleotide resolution. SHAPE takes advantage of the discovery that the nucleophilic reactivity of the ribose 2'-hydroxyl group is modulated by local nucleotide flexibility in the RNA backbone. Flexible nucleotides are reactive toward hydroxyl-selective electrophiles, whereas constrained nucleotides are unreactive. Initial versions of SHAPE chemistry, which employ isatoic anhydride derivatives that react on the minute time scale, are emerging as the ideal technology for monitoring equilibrium structures of RNA in a wide variety of biological environments. Here, we extend SHAPE chemistry to a benzoyl cyanide scaffold to make possible facile time-resolved kinetic studies of RNA in approximately 1 s snapshots. We then use SHAPE chemistry to follow the time-dependent folding of an RNase P specificity domain RNA. Tertiary interactions form in two distinct steps with local tertiary contacts forming an order of magnitude faster than long-range interactions. Rate-determining tertiary folding requires minutes despite that no non-native interactions must be disrupted to form the native structure. Instead, overall folding is limited by simultaneous formation of interactions approximately 55 A distant in the RNA. Time-resolved SHAPE holds broad potential for understanding structural biogenesis and the conformational interconversions essential to the functions of complex RNA molecules at single nucleotide resolution.     

Fingerprinting Noncanonical and Tertiary RNA Structures by Differential SHAPE Reactivity

Many RNA structures are composed of simple secondary structure elements linked by a few critical tertiary interactions. SHAPE chemistry has made interrogation of RNA dynamics at single-nucleotide resolution straightforward. However, de novo identification of nucleotides involved in tertiary interactions remains a challenge. Here we show that nucleotides that form noncanonical or tertiary contacts can be detected by comparing information obtained using two SHAPE reagents, N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6). Nucleotides that react preferentially with NMIA exhibit slow local nucleotide dynamics and usually adopt the less common C2'-endo ribose conformation. Experiments and first-principles calculations show that 1M6 reacts preferentially with nucleotides in which one face of the nucleobase allows an unhindered stacking interaction with the reagent. Differential SHAPE reactivities were used to detect noncanonical and tertiary interactions in four RNAs with diverse structures and to identify preformed noncanonical interactions in partially folded RNAs. Differential SHAPE reactivity analysis will enable experimentally concise, large-scale identification of tertiary structure elements and ligand binding sites in complex RNAs and in diverse biological environments.     

A powerful approach for the selection of 2-aminopurine substitution sites to investigate RNA folding

A precise tertiary structure must be adopted to allow the function of many RNAs in cells. Accordingly, increasing resources have been devoted to the elucidation of RNA structures and the folding of RNAs. 2-Aminopurine (2AP), a fluorescent nucleobase analogue, can be substituted in strategic positions of DNA or RNA molecules to act as site-specific probe to monitor folding and folding dynamics of nucleic acids. Recent studies further demonstrated the potential of 2AP modifications in the assessment of folding kinetics during ligand-induced secondary and tertiary RNA structure rearrangements. However, an efficient way to unambiguously identify reliable positions for 2AP sensors is as yet unavailable and would represent a major asset, especially in the absence of crystallographic or NMR structural data for a target molecule. We report evidence of a novel and direct correlation between the 2'-OH flexibility of nucleotides, observed by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probing and the fluorescence response following nucleotide substitutions by 2AP. This correlation leads to a straightforward method, using SHAPE probing with benzoyl cyanide, to select appropriate nucleotide sites for 2AP substitution. This clear correlation is presented for three model RNAs of biological significance: the SAM-II, adenine (addA), and preQ(1) class II (preQ(1)cII) riboswitches.     

Femtomole SHAPE reveals regulatory structures in the authentic XMRV RNA genome

Higher-order structure influences critical functions in nearly all noncoding and coding RNAs. Most single-nucleotide resolution RNA structure determination technologies cannot be used to analyze RNA from scarce biological samples, like viral genomes. To make quantitative RNA structure analysis applicable to a much wider array of RNA structure-function problems, we developed and applied high-sensitivity selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) to structural analysis of authentic genomic RNA of the xenotropic murine leukemia virus-related virus (XMRV). For analysis of fluorescently labeled cDNAs generated in high-sensitivity SHAPE experiments, we developed a two-color capillary electrophoresis approach with zeptomole molecular detection limits and subfemtomole sensitivity for complete SHAPE experiments involving hundreds of individual RNA structure measurements. High-sensitivity SHAPE data correlated closely (R = 0.89) with data obtained by conventional capillary electrophoresis. Using high-sensitivity SHAPE, we determined the dimeric structure of the XMRV packaging domain, examined dynamic interactions between the packaging domain RNA and viral nucleocapsid protein inside virion particles, and identified the packaging signal for this virus. Despite extensive sequence differences between XMRV and the intensively studied Moloney murine leukemia virus, architectures of the regulatory domains are similar and reveal common principles of gammaretrovirus RNA genome packaging.     

A deoxyribozyme that forms a three-helix-junction complex with its RNA substrates and has general RNA branch-forming activity

We recently used in vitro selection to identify 7S11, a deoxyribozyme that synthesizes 2',5'-branched RNA. The 7S11 DNA enzyme mediates the nucleophilic attack of an adenosine 2'-hydroxyl group at a 5'-triphosphate, forming 2',5'-branched RNA in a reaction that resembles the first step of in vivo RNA splicing. Here, we describe 7S11 characterization experiments that have two important implications for nucleic acid chemistry and biochemistry. First, on the basis of a comprehensive analysis of its substrate sequence requirements, 7S11 is shown to be generally applicable for the synthesis of a wide range of 2',5'-branched RNAs. Strict substrate sequence requirements are found at the two RNA nucleotides that directly form the branched linkage, and these requirements correspond to those nucleotides found most commonly at these two positions in natural spliced RNAs. Outside of these two nucleotides, most substrate sequences are tolerated with useful ligation activity, although rates and yields vary. Because chemical synthesis approaches to branched RNA are extremely limited in scope, the deoxyribozyme-based route using 7S11 will enable many experiments that require branched RNA. Second, comprehensive nucleotide covariation experiments demonstrate that 7S11 and its RNA substrates adopt a three-helix-junction structure in which the branch-site nucleotide is located at the intersection of the three helices. Because many natural ribozymes have multi-helix junctions, 7S11 is an interesting model system for catalytic nucleic acids.     

Dynamics of large elongated RNA by NMR carbon relaxation

We present an NMR strategy for characterizing picosecond-to-nanosecond internal motions in uniformly 13C/15N-labeled RNAs that combines measurements of R1, R1rho, and heteronuclear 13C{1H} NOEs for protonated base (C2, C5, C6, and C8) and sugar (C1') carbons with a domain elongation strategy for decoupling internal from overall motions and residual dipolar coupling (RDC) measurements for determining the average RNA global conformation and orientation of the principal axis of the axially symmetric rotational diffusion. TROSY-detected pulse sequences are presented for the accurate measurement of nucleobase carbon R1 and R1rho rates in large RNAs. The relaxation data is analyzed using a model free formalism which takes into account the very high anisotropy of overall rotational diffusion (Dratio approximately 4.7), asymmetry of the nucleobase CSAs and noncollinearity of C-C, C-H dipolar and CSA interactions under the assumption that all interaction tensors for a given carbon experience identical isotropic internal motions. The approach is demonstrated and validated on an elongated HIV-1 TAR RNA (taum approximately 18 ns) both in free form and bound to the ligand argininamide (ARG). Results show that, while ARG binding reduces the amplitude of collective helix motions and local mobility at the binding pocket, it leads to a drastic increase in the local mobility of "spacer" bulge residues linking the two helices which undergo virtually unrestricted internal motions (S2 approximately 0.2) in the ARG bound state. Our results establish the ability to quantitatively study the dynamics of RNAs which are significantly larger and more anisotropic than customarily studied by NMR carbon relaxation.     

Deoxyribozymes with 2'-5' RNA ligase activity

In vitro selection was used to identify deoxyribozymes that ligate two RNA substrates. In the ligation reaction, a 2'-5' RNA phosphodiester linkage is created from a 2',3'-cyclic phosphate and a 5'-hydroxyl group. The new Mg(2+)-dependent deoxyribozymes provide 50-60% yield of ligated RNA in overnight incubations at pH 7.5 and 37 degrees C, and they afford 40-50% yield in 1 h at pH 9.0 and 37 degrees C. Various RNA substrate sequences may be joined by simple Watson-Crick covaration of the DNA binding arms that interact with the two RNA substrates. The current deoxyribozymes have some RNA substrate sequence requirements at the nucleotides immediately surrounding the ligation junction (either UAUA GGAA or UAUN GGAA, where the arrow denotes the ligation site and N equals any nucleotide). One of the new deoxyribozymes was used to prepare by ligation the Tetrahymena group I intron RNA P4-P6 domain, a representative structured RNA. Nondenaturing gel electrophoresis revealed that a 2'-5' linkage between nucleotides A233 and G234 of P4-P6 does not disrupt its Mg(2+)-dependent folding (DeltaDeltaG degrees ' < 0.2 kcal/mol). This demonstrates that a 2'-5' linkage does not necessarily interfere with structure in a folded RNA. Therefore, these non-native linkages may be acceptable in modified RNAs when structure/function relationships are investigated. Deoxyribozymes that ligate RNA should be particularly useful for preparing site-specifically modified RNAs for studies of RNA structure, folding, and catalysis.     

Aminoglycoside microarrays to explore interactions of antibiotics with RNAs and proteins

RNA is an important target for drug discovery efforts. Several clinically used aminoglycoside antibiotics bind to bacterial rRNA and inhibit protein synthesis. Aminoglycosides, however, are losing efficacy due to their inherent toxicity and the increase in antibiotic resistance. Targeting of other RNAs is also becoming more attractive thanks to the discovery of new potential RNA drug targets through genome sequencing and biochemical efforts. Identification of new compounds that target RNA is therefore urgent, and we report here on the development of rapid screening methods to probe binding of low molecular weight ligands to proteins and RNAs. A series of aminoglycosides has been immobilized onto glass microscope slides, and binding to proteins and RNAs has been detected by fluorescence. Construction and analysis of the arrays is completed by standard DNA genechip technology. Binding of immobilized aminoglycosides to proteins that are models for study of aminoglycoside toxicity (DNA polymerase and phospholipase C), small RNA oligonucleotide mimics of aminoglycoside binding sites in the ribosome (rRNA A-site mimics), and a large (approximately 400 nucleotide) group I ribozyme RNA is detected. The ability to screen large RNAs alleviates many complications associated with binding experiments that use isolated truncated regions from larger RNAs. These studies lay the foundation for rapid identification of small organic ligands from combinatorial libraries that exhibit strong and selective RNA binding while displaying decreased affinity to toxicity-causing proteins.     

Molecular dynamics study of the ribosomal A-site

Many aminoglycosidic antibiotics target the A-site of 16S RNA in the small ribosomal subunit and affect the fidelity of protein translation in bacteria. Upon binding, aminoglycosides displace two adenines (A1492 and A1493 for E. coli numbering) that are involved in tRNA anticodon loop recognition. The major difference in the aminoglycosidic binding site between the prokaryota and eukaryota is an adenine into guanine substitution in the position 1408. This mutation likely affects the dynamics of near A1492 and A1493 and hinders the binding of aminoglycosides to eukaryotic ribosomes. With multiple 20 ns long all-atom molecular dynamics simulations, we study the flexibility of a 22 nucleotide RNA fragment which mimics the aminoglycosidic binding site. Simulations are carried out for both native and A1408G mutated RNA as well as for their complexes with aminoglycosidic representative paromomycin. We observe intra- and extrahelical configurations of A1492 and A1493, which differ between the prokaryotic and the mutated structure. We obtain configurations of the A-site that are also observed in the NMR and crystal structures. Our studies show the differences in the internal mobility of the A-site, as well as that in ion and water density distributions inside of the binding cleft, between the prokaryotic and mutated RNA. We also compare the performance of two force field parameters for RNA, Amber and Charmm.     

Probing conformational variations at the ATPase site of the RNA helicase DbpA by high-field electron-nuclear double resonance spectroscopy

The RNA helicase DbpA promotes RNA remodeling coupled to ATP hydrolysis. It is unique because of its specificity to hairpin 92 of 23S rRNA (HP92). Although DbpA kinetic pathways leading to ATP hydrolysis and RNA unwinding have been recently elucidated, the molecular (atomic) basis for the coupling of ATP hydrolysis to RNA remodeling remains unclear. This is, in part, due to the lack of detailed structural information on the ATPase site in the presence and absence of RNA in solution. We used high-field pulse ENDOR (electron-nuclear double resonance) spectroscopy to detect and analyze fine conformational changes in the protein's ATPase site in solution. Specifically, we substituted the essential Mg(2+) cofactor in the ATPase active site for paramagnetic Mn(2+) and determined its close environment with different nucleotides (ADP, ATP, and the ATP analogues ATPγS and AMPPnP) in complex with single- and double-stranded RNA. We monitored the Mn(2+) interactions with the nucleotide phosphates through the (31)P hyperfine couplings and the coordination by protein residues through (13)C hyperfine coupling from (13)C-enriched DbpA. We observed that the nucleotide binding site of DbpA adopts different conformational states upon binding of different nucleotides. The ENDOR spectra revealed a clear distinction between hydrolyzable and nonhydrolyzable nucleotides prior to RNA binding. Furthermore, both the (13)C and the (31)P ENDOR spectra were found to be highly sensitive to changes in the local environment of the Mn(2+) ion induced by the hydrolysis. More specifically, ATPγS was efficiently hydrolyzed upon binding of RNA, similar to ATP. Importantly, the Mn(2+) cofactor remains bound to a single protein side chain and to one or two nucleotide phosphates in all complexes, whereas the remaining metal coordination positions are occupied by water. The conformational changes in the protein's ATPase active site associated with the different DbpA states occur in remote coordination shells of the Mn(2+) ion. Finally, a competitive Mn(2+) binding site was found for single-stranded RNA construct.     

2'-mercaptonucleotide interference reveals regions of close packing within folded RNA molecules

The 2'-hydroxyl group makes essential contributions to RNA structure and function. As an approach to assess the ability of a mercapto group to serve as a functional analogue for the 2'-hydroxyl group, we synthesized 2'-mercaptonucleotides for use in nucleotide analogue interference mapping. To correlate the observed interference effects with tertiary structure, we used the independently folding DeltaC209 P4-P6 domain from the Tetrahymena group I intron. We generated populations of DeltaC209 P4-P6 molecules containing 2'-mercaptonucleotides located randomly throughout the domain and separated the folded molecules from the unfolded molecules by nondenaturing gel electrophoresis. Iodine-induced cleavage of the RNA molecules revealed the sites at which 2'-mercaptonucleotides interfere with folding. These interferences cluster in the most densely packed regions of the tertiary structure, occurring only at sites that lack the space and flexibility to accommodate a sulfur atom. Interference mapping with 2'-mercaptonucleotides therefore provides a method by which to identify structurally rigid and densely packed regions within folded RNA molecules.     

Probing RNA hairpins with cobalt(III)hexammine and electrospray ionization mass spectrometry

In this work, electrospray ionization mass spectrometry (ESI MS) was employed to study the interactions of cobalt(III) hexammine, Co(NH3)6(3+), with five RNA hairpins representing the 790 loop of 16S ribosomal RNA and 1920 loop of 23S ribosomal RNA. The RNAs varied in mismatch identity (G.U versus A.C) and level of base modification (pseudouridine versus uridine). Co(NH3)6(3+) binding was observed with the four RNA hairpins that contained a G.U wobble pair in the stem region. ESI MS revealed 1:1 and 1:2 complex formation with all RNAs. Weaker binding was observed with the fifth RNA hairpin that contained an A.C wobble pair in the stem region. The effects of pH on Co(NH3)6(3+) binding were also examined.     

Synthesis of 2'-N-methylamino-2'-deoxyguanosine and 2'-N,N-dimethylamino-2'-deoxyguanosine and their incorporation into RNA by phosphoramidite chemistry

The 2'-hydroxyl groups within RNA contribute in essential ways to RNA structure and function. Previously, we designed an atomic mutation cycle (AMC) that uses ribonucleoside analogues bearing different C-2'-substituents, including -OCH(3), -NH(2), -NHMe, and -NMe(2), to identify hydroxyl groups within RNA that donate functionally significant hydrogen bonds. To enable AMC analysis of the nucleophilic guanosine cofactor in the Tetrahymena ribozyme reaction and at other guanosines whose 2'-hydroxyl groups impart critical functional contributions, we describe here the syntheses of 2'-methylamino-2'-deoxyguanosine (G(NHMe)) and 2'-N,N-dimethylamino-2'-deoxyguanosine (G(NMe(2))) and their corresponding phosphoramidites. The key step in obtaining the nucleosides involved S(N)2 displacement of 2'-β-triflate from an appropriate guanosine derivative by methylamine or dimethylamine. We readily obtained the G(NMe(2)) phosphoramidite and incorporated it into RNA. However, the G(NHMe) phosphoramidite posed a significantly greater challenge due to lack of a suitable -2'-NHMe protecting group. After testing several strategies, we established that allyloxycarbonyl (Alloc) provided suitable protection for 2'-N-methylamino group during the phosphoramidite synthesis and the subsequent RNA synthesis. This work enables AMC analysis of guanosine's 2'-hydroxyl group within RNA.