Complex Ligand-Induced Conformational Changes in tRNAAsp Revealed by Single-Nucleotide Resolution SHAPE Chemistry

2'-OH as a universal handle for studying intracellular RNAs

RNA plays pivotal roles in most cellular processes, serving as both the traditional carrier of genetic information and as a key regulator of cellular functions. The advent of chemical technologies has contributed critically to the analysis of cellular RNA structures, functions, and interactions. Many of these methods and molecules involve the utilization of chemically reactive handles in RNAs, either introduced externally or inherent within the polymer itself. Among these handles, the 2'-hydroxyl (2'-OH) group has emerged as an exceptionally well-suited and general chemical moiety for the modification and profiling of RNAs in intracellular studies. In this review, we provide an overview of the recent advancements in intracellular applications of acylation at the 2'-OH group of RNA. We outline progress made in probing RNA structure and interactomes, controlling RNA function, RNA imaging, and analyzing RNA-small molecule interactions, all achieved in living cells through this simple chemical handle on the biopolymer.

Genomic and transcriptomic profiling of phoenix colonies

Pseudomonas aeruginosa is a Gram-negative bacterium responsible for numerous human infections. Previously, novel antibiotic tolerant variants known as phoenix colonies as well as variants similar to viable but non-culturable (VBNC) colonies were identified in response to high concentrations of aminoglycosides. In this study, the mechanisms behind phoenix colony and VBNC-like colony emergence were further explored using both whole genome sequencing and RNA sequencing. Phoenix colonies were found to have a single nucleotide polymorphism (SNP) in the PA4673 gene, which is predicted to encode a GTP-binding protein. No SNPs were identified within VBNC-like colonies compared to the founder population. RNA sequencing did not detect change in expression of PA4673 but revealed multiple differentially expressed genes that may play a role in phoenix colony emergence. One of these differentially expressed genes, PA3626, encodes for a tRNA pseudouridine synthase which when knocked out led to a complete lack of phoenix colonies. Although not immediately clear whether the identified genes in this study may have interactions which have not yet been recognized, they may contribute to the understanding of how phoenix colonies are able to emerge and survive in the presence of antibiotic exposure.

RNA-Ligand Interactions Quantified by Surface Plasmon Resonance with Reference Subtraction

Structured RNAs bind ligands and are attractive targets for small-molecule drugs. A wide variety of analytical methods have been used to characterize RNA-ligand interactions, but our experience is that most have significant limitations in terms of material requirements and applicability to complex RNAs. Surface plasmon resonance (SPR) potentially overcomes these limitations, but we find that the standard experimental framework measures notable nonspecific electrostatic-mediated interactions, frustrating analysis of weak RNA binders. SPR measurements are typically quantified relative to a non-target reference channel. Here, we show that referencing to a channel containing a non-binding control RNA enables subtraction of nonspecific binding contributions, allowing measurements of accurate and specific binding affinities. We validated this approach for small-molecule binders of two riboswitch RNAs with affinities ranging from nanomolar to millimolar, including low-molecular-mass fragment ligands. SPR implemented with reference subtraction reliably discriminates specific from nonspecific binding, uses RNA and ligand material efficiently, and enables rapid exploration of the ligand-binding landscape for RNA targets.

Subsite Ligand Recognition and Cooperativity in the TPP Riboswitch: Implications for Fragment-Linking in RNA Ligand Discovery

RNA molecules can show high levels of cooperativity in their global folding and interactions with divalent ions. However, cooperativity at individual ligand-RNA interaction sites remains poorly understood. Here, we investigated the binding of thiamine and methylene diphosphonic acid (MDP, a soluble structural analogue of pyrophosphate) to the thiamine pyrophosphate riboswitch. These ligands each bind weakly at proximal subsites, with 10 μM and 1 mM affinities, respectively. The affinity of MDP moderately improves when thiamine or thiamine-like fragments are pre-bound to the RNA. Covalent linking of thiamine and MDP substantially increases riboswitch binding to a notable high affinity of 20 nM. Crystal structures and single-molecule correlated chemical probing revealed favorable induced fit effects upon binding of individual ligands and, unexpectedly, a substantial thermodynamically unfavorable RNA structural rearrangement upon binding of the linked thiamine-MDP ligand. Thus, linking of two ligands of modest affinity, accompanied by an unfavorable structural rearrangement, still yields a potent linked RNA-binding compound. Since complex ligands often bind riboswitches and other RNAs at proximal subsites, principles derived from this work inform and support fragment-linking strategies for identifying small molecules that interact with RNA specifically and with high affinity.

Ions in the Deep Subsurface of Earth, Mars, and Icy Moons: Their Effects in Combination with Temperature and Pressure on tRNA-Ligand Binding

The interactions of ligands with nucleic acids are central to numerous reactions in the biological cell. How such reactions are affected by harsh environmental conditions such as low temperatures, high pressures, and high concentrations of destructive ions is still largely unknown. To elucidate the ions’ role in shaping habitability in extraterrestrial environments and the deep subsurface of Earth with respect to fundamental biochemical processes, we investigated the effect of selected salts (MgCl₂, MgSO₄, and Mg(ClO₄)₂) and high hydrostatic pressure (relevant for the subsurface of that planet) on the complex formation between tRNA and the ligand ThT. The results show that Mg²⁺ salts reduce the binding tendency of ThT to tRNA. This effect is largely due to the interaction of ThT with the salt anions, which leads to a strong decrease in the activity of the ligand. However, at mM concentrations, binding is still favored. The ions alter the thermodynamics of binding, rendering complex formation that is more entropy driven. Remarkably, the pressure favors ligand binding regardless of the type of salt. Although the binding constant is reduced, the harsh conditions in the subsurface of Earth, Mars, and icy moons do not necessarily preclude nucleic acid–ligand interactions of the type studied here.

Mapping the RNA Chaperone Activity of the T. brucei Editosome Using SHAPE Chemical Probing

Mitochondrial pre-mRNAs in African trypanosomes adopt intricately folded, highly stable 2D and 3D structures. The RNA molecules are substrates of a U-nucleotide-specific insertion/deletion-type RNA editing reaction, which is catalyzed by a 0.8 MDa protein complex known as the editosome. RNA binding to the editosome is followed by a chaperone-mediated RNA remodeling reaction. The reaction increases the dynamic of specifically U-nucleotides to lower their base-pairing probability and as a consequence generates a simplified RNA folding landscape that is critical for the progression of the editing reaction cycle. Here we describe a chemical mapping method to quantitatively monitor the chaperone-driven structural changes of pre-edited mRNAs upon editosome binding. The method is known as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). SHAPE is based on the differential electrophilic modification of ribose 2'-hydroxyl groups in structurally constraint (double-stranded) versus structurally unconstrained (single-stranded) nucleotides. Electrophilic anhydrides such as 1-methyl-7-nitroisatoic anhydride are used as probing reagents, and the ribose 2'-modified nucleotides are mapped as abortive cDNA synthesis products. As a result, SHAPE allows the identification of all single-stranded and base-paired regions in a given RNA, and the data are used to compute experimentally derived RNA 2D structures. A side-by-side comparison of the RNA 2D folds in the pre- and post-chaperone states finally maps the chaperone-induced dynamic of the different pre-mRNAs with single-nucleotide resolution.

SHAPE Directed Discovery of New Functions in Large RNAs

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

Aminoacyl-tRNA synthetases as drug targets

Aminoacyl-tRNA synthetases (AARSs) have been considered very attractive drug-targets for decades. This interest probably emerged with the identification of differences in AARSs between prokaryotic and eukaryotic species, which provided a rationale for the development of antimicrobials targeting bacterial AARSs with minimal effect on the homologous human AARSs. Today we know that AARSs are not only attractive, but also valid drug targets as they are housekeeping proteins that: (i) play a fundamental role in protein translation by charging the corresponding amino acid to its cognate tRNA and preventing mistranslation mistakes [1], a critical process during fast growing conditions of microbes; and (ii) present significant differences between microbes and humans that can be used for drug development [2]. Together with the vast amount of available data on both pathogenic and mammalian AARSs, it is expected that, in the future, the numerous reported inhibitors of AARSs will provide the basis to develop new therapeutics for the treatment of human diseases. In this chapter, a detailed summary on the state-of-the-art in drug discovery and drug development for each aminoacyl-tRNA synthetase will be presented.

Aminoglycoside antibiotics can inhibit or activate twister ribozyme cleavage

Interactions between aminoglycoside antibiotics and the twister ribozyme were investigated in this study. An initial screen of 17 RNA-binding antibiotics showed that a number of aminoglycosides inhibit the ribozyme, while a subset of aminoglycosides enhances twister cleavage. Initial kinetic analysis of the twister ribozyme showed a sevenfold inhibition of ribozyme cleavage by paromomycin and a fivefold enhancement of cleavage by sisomicin. Direct binding between the twister ribozyme RNA and paromomycin or sisomicin was measured by microscale thermophoresis. Selective 2'-hydroxyl acylation analysed by primer extension shows that both paromomycin and sisomicin induce distinctive tertiary structure changes to the twister ribozyme. Published crystal structures and mechanistic analysis of the twister ribozyme have deduced a nucleobase-mediated general acid-base catalytic mechanism, in which a conserved guanine plays a key role. Here, we show that paromomycin binding induces a structural transition to the twister ribozyme such that a highly conserved guanine in the active site becomes displaced, leading to inhibition of cleavage. In contrast, sisomicin binding appears to change interactions between P3 and L2, inducing allosteric changes to the active site that enhance twister RNA cleavage. Therefore, we show that small-molecule binding can modulate twister ribozyme activity. These results suggest that aminoglycosides may be used as molecular tools to study this widely distributed ribozyme.

NMR characterization of RNA small molecule interactions

Exciting discoveries of naturally occurring ligand-sensing and disease-linked noncoding RNAs have promoted significant interests in understanding RNA-small molecule interactions. NMR spectroscopy is a powerful tool for characterizing intermolecular interactions. In this review, we describe protocols and approaches for applying NMR spectroscopy to investigate interactions between RNA and small molecules. We review protocols for RNA sample preparation, methods for identifying RNA-binding small molecules, approaches for mapping RNA-small molecule interactions, determining complex structures, and characterizing binding kinetics. We hope this review will provide a guideline to streamline NMR applications in studying RNA-small molecule interactions, facilitating both basic mechanistic understandings of RNA functions and translational efforts in developing RNA-targeted therapeutics.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

RNA Structural Differentiation: Opportunities with Pattern Recognition

Our awareness and appreciation of the many regulatory roles of RNA have dramatically increased in the past decade. This understanding, in addition to the impact of RNA in many disease states, has renewed interest in developing selective RNA-targeted small molecule probes. However, the fundamental guiding principles in RNA molecular recognition that could accelerate these efforts remain elusive. While high-resolution structural characterization can provide invaluable insight, examples of well-characterized RNA structures, not to mention small molecule:RNA complexes, remain limited. This Perspective provides an overview of the current techniques used to understand RNA molecular recognition when high-resolution structural information is unavailable. We will place particular emphasis on a new method, pattern recognition of RNA with small molecules (PRRSM), that provides rapid insight into critical components of RNA recognition and differentiation by small molecules as well as into RNA structural features.

Drugging tRNA aminoacylation

Inhibition of tRNA aminoacylation has proven to be an effective antimicrobial strategy, impeding an essential step of protein synthesis. Mupirocin, the well-known selective inhibitor of bacterial isoleucyl-tRNA synthetase, is one of three aminoacylation inhibitors now approved for human or animal use. However, design of novel aminoacylation inhibitors is complicated by the steadfast requirement to avoid off-target inhibition of protein synthesis in human cells. Here we review available data regarding known aminoacylation inhibitors as well as key amino-acid residues in aminoacyl-tRNA synthetases (aaRSs) and nucleotides in tRNA that determine the specificity and strength of the aaRS-tRNA interaction. Unlike most ligand-protein interactions, the aaRS-tRNA recognition interaction represents coevolution of both the tRNA and aaRS structures to conserve the specificity of aminoacylation. This property means that many determinants of tRNA recognition in pathogens have diverged from those of humans-a phenomenon that provides a valuable source of data for antimicrobial drug development.

Multiplex Aptamer Discovery through Apta-Seq and Its Application to ATP Aptamers Derived from Human-Genomic SELEX

Laboratory-evolved RNAs bind a wide variety of targets and serve highly diverse functions, including as diagnostic and therapeutic aptamers. The majority of aptamers have been identified using in vitro selection (SELEX), a molecular evolution technique based on selecting target-binding RNAs from highly diverse pools through serial rounds of enrichment and amplification. In vitro selection typically yields multiple distinct motifs of highly variable abundance and target-binding affinities. The discovery of new aptamers is often limited by the difficulty of characterizing the selected motifs, because testing of individual sequences tends to be a tedious process. To facilitate the discovery of new aptamers within in vitro selected pools, we developed Apta-Seq, a multiplex analysis based on quantitative, ligand-dependent 2' acylation of solvent-accessible regions of the selected RNA pools, followed by reverse transcription (SHAPE) and deep sequencing. The method reveals, in a single sequencing experiment, the identity, structural features, and target dissociation constants for aptamers present in the selected pool. Application of Apta-Seq to a human genomic pool enriched for ATP-binding RNAs yielded three new aptamers, which together with previously identified human aptamers suggest that ligand-binding RNAs may be common in mammals.

Thermodynamic and Kinetic Analyses of Iron Response Element (IRE)-mRNA Binding to Iron Regulatory Protein, IRP1

Comparison of kinetic and thermodynamic properties of IRP1 (iron regulatory protein1) binding to FRT (ferritin) and ACO2 (aconitase2) IRE-RNAs, with or without Mn²⁺, revealed differences specific to each IRE-RNA. Conserved among animal mRNAs, IRE-RNA structures are noncoding and bind Fe²⁺ to regulate biosynthesis rates of the encoded, iron homeostatic proteins. IRP1 protein binds IRE-RNA, inhibiting mRNA activity; Fe²⁺ decreases IRE-mRNA/IRP1 binding, increasing encoded protein synthesis. Here, we observed heat, 5 °C to 30 °C, increased IRP1 binding to IRE-RNA 4-fold (FRT IRE-RNA) or 3-fold (ACO2 IRE-RNA), which was enthalpy driven and entropy favorable. Mn²⁺ (50 µM, 25 °C) increased IRE-RNA/IRP1 binding (K_d) 12-fold (FRT IRE-RNA) or 6-fold (ACO2 IRE-RNA); enthalpic contributions decreased ~61% (FRT) or ~32% (ACO2), and entropic contributions increased ~39% (FRT) or ~68% (ACO2). IRE-RNA/IRP1 binding changed activation energies: FRT IRE-RNA 47.0 ± 2.5 kJ/mol, ACO2 IRE-RNA 35.0 ± 2.0 kJ/mol. Mn²⁺ (50 µM) decreased the activation energy of RNA-IRP1 binding for both IRE-RNAs. The observations suggest decreased RNA hydrogen bonding and changed RNA conformation upon IRP1 binding and illustrate how small, conserved, sequence differences among IRE-mRNAs selectively influence thermodynamic and kinetic selectivity of the protein/RNA interactions.

Osmium tetroxide as a probe of RNA structure

Structured RNAs have a central role in cellular function. The capability of structured RNAs to adopt fixed architectural structures or undergo dynamic conformational changes contributes to their diverse role in the regulation of gene expression. Although numerous biophysical and biochemical tools have been developed to study structured RNAs, there is a continuing need for the development of new methods for the investigation of RNA structures, especially methods that allow RNA structure to be studied in solution close to its native cellular conditions. Here we use osmium tetroxide (OsO₄) as a chemical probe of RNA structure. In this method, we have used fluorescence-based sequencing technologies to detect OsO₄ modified RNA. We characterized the requirements for OsO₄ modification of RNA by investigating three known structured RNAs: the M-box, glycine riboswitch RNAs, and tRNA^asp Our results show that OsO₄ predominantly modifies RNA at uracils that are conformationally exposed on the surface of the RNA. We also show that changes in OsO₄ reactivity at flexible positions in the RNA correlate with ligand-driven conformational changes in the RNA structure. Osmium tetroxide modification of RNA will provide insights into the structural features of RNAs that are relevant to their underlying biological functions.

Challenge of mimicking the influences of the cellular environment on RNA structure by PEG-induced macromolecular crowding

There are large differences between the cellular environment and the conditions widely used to study RNA in vitro. SHAPE RNA structure probing in Escherichia coli cells has shown that the cellular environment stabilizes both long-range and local tertiary interactions in the adenine riboswitch aptamer domain. Synthetic crowding agents are widely used to understand the forces that stabilize RNA structure and in efforts to recapitulate the cellular environment under simplified experimental conditions. Here, we studied the structure and ligand binding ability of the adenine riboswitch in the presence of the macromolecular crowding agent, polyethylene glycol (PEG). Ethylene glycol and low-molecular mass PEGs destabilized RNA structure and caused the riboswitch to sample secondary structures different from those observed in simple buffered solutions or in cells. In the presence of larger PEGs, longer-range loop-loop interactions were more similar to those in cells than in buffer alone, consistent with prior work showing that larger PEGs stabilize compact RNA states. Ligand affinity was weakened by low-molecular mass PEGs but increased with high-molecular mass PEGs, indicating that PEG cosolvents exert complex chemical and steric effects on RNA structure. Regardless of polymer size, however, nucleotide-resolution structural characteristics observed in cells were not recapitulated in PEG solutions. Our results reveal that the cellular environment is difficult to recapitulate in vitro; mimicking the cellular state will likely require a combination of crowding agents and other chemical species.

Selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistries exploit small electrophilic reagents that react with 2'-hydroxyl groups to interrogate RNA structure at single-nucleotide resolution. Mutational profiling (MaP) identifies modified residues by using reverse transcriptase to misread a SHAPE-modified nucleotide and then counting the resulting mutations by massively parallel sequencing. The SHAPE-MaP approach measures the structure of large and transcriptome-wide systems as accurately as can be done for simple model RNAs. This protocol describes the experimental steps, implemented over 3 d, that are required to perform SHAPE probing and to construct multiplexed SHAPE-MaP libraries suitable for deep sequencing. Automated processing of MaP sequencing data is accomplished using two software packages. ShapeMapper converts raw sequencing files into mutational profiles, creates SHAPE reactivity plots and provides useful troubleshooting information. SuperFold uses these data to model RNA secondary structures, identify regions with well-defined structures and visualize probable and alternative helices, often in under 1 d. SHAPE-MaP can be used to make nucleotide-resolution biophysical measurements of individual RNA motifs, rare components of complex RNA ensembles and entire transcriptomes.

tRNAs as antibiotic targets

Transfer RNAs (tRNAs) are central players in the protein translation machinery and as such are prominent targets for a large number of natural and synthetic antibiotics. This review focuses on the role of tRNAs in bacterial antibiosis. We will discuss examples of antibiotics that target multiple stages in tRNA biology from tRNA biogenesis and modification, mature tRNAs, aminoacylation of tRNA as well as prevention of proper tRNA function by small molecules binding to the ribosome. Finally, the role of deacylated tRNAs in the bacterial “stringent response” mechanism that can lead to bacteria displaying antibiotic persistence phenotypes will be discussed.

Rapid kinetics of iron responsive element (IRE) RNA/iron regulatory protein 1 and IRE-RNA/eIF4F complexes respond differently to metal ions

Metal ion binding was previously shown to destabilize IRE-RNA/IRP1 equilibria and enhanced IRE-RNA/eIF4F equilibria. In order to understand the relative importance of kinetics and stability, we now report rapid rates of protein/RNA complex assembly and dissociation for two IRE-RNAs with IRP1, and quantitatively different metal ion response kinetics that coincide with the different iron responses in vivo. k_on, for FRT IRE-RNA binding to IRP1 was eight times faster than ACO2 IRE-RNA. Mn²⁺ decreased k_on and increased k_off for IRP1 binding to both FRT and ACO2 IRE-RNA, with a larger effect for FRT IRE-RNA. In order to further understand IRE-mRNA regulation in terms of kinetics and stability, eIF4F kinetics with FRT IRE-RNA were determined. k_on for eIF4F binding to FRT IRE-RNA in the absence of metal ions was 5-times slower than the IRP1 binding to FRT IRE-RNA. Mn²⁺ increased the association rate for eIF4F binding to FRT IRE-RNA, so that at 50 µM Mn²⁺ eIF4F bound more than 3-times faster than IRP1. IRP1/IRE-RNA complex has a much shorter life-time than the eIF4F/IRE-RNA complex, which suggests that both rate of assembly and stability of the complexes are important, and that allows this regulatory system to respond rapidly to change in cellular iron.

Safer one-pot synthesis of the 'SHAPE' reagent 1-methyl-7-nitroisatoic anhydride (1m7)

Estimating the reactivity of 2'-hydroxyl groups along an RNA chain of interest aids in the modeling of the folded RNA structure; flexible loops tend to be reactive, whereas duplex regions are generally not. Among the most useful reagents for probing 2'-hydroxyl reactivity is 1-methyl-7-nitroisatoic anhydride (1m7), but the absence of a reliable, inexpensive source has prevented widespread adoption. An existing protocol for the conversion of an inexpensive precursor 4-nitroisatoic anhydride (4NIA) recommends the use of NaH in dimethylformamide (DMF), a reagent combination that most molecular biology labs are not equipped to handle, and that does not scale safely in any case. Here we describe a safer, one-pot method for bulk conversion of 4NIA to 1m7 that reduces costs and bypasses the use of NaH. We show that 1m7 produced by this method is free of side products and can be used to probe RNA structure in vitro.

The cellular environment stabilizes adenine riboswitch RNA structure

There are large differences between the intracellular environment and the conditions widely used to study RNA structure and function in vitro. To assess the effects of the crowded cellular environment on RNA, we examined the structure and ligand binding function of the adenine riboswitch aptamer domain in healthy, growing Escherichia coli cells at single-nucleotide resolution on the minute time scale using SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension). The ligand-bound aptamer structure is essentially the same in cells and in buffer at 1 mM Mg(2+), the approximate Mg(2+) concentration we measured in cells. In contrast, the in-cell conformation of the ligand-free aptamer is much more similar to the fully folded ligand-bound state. Even adding high Mg(2+) concentrations to the buffer used for in vitro analyses did not yield the conformation observed for the free aptamer in cells. The cellular environment thus stabilizes the aptamer significantly more than does Mg(2+) alone. Our results show that the intracellular environment has a large effect on RNA structure that ultimately favors highly organized conformations.

Riboswitch control of aminoglycoside antibiotic resistance

The majority of riboswitches are regulatory RNAs that regulate gene expression by binding small-molecule metabolites. Here we report the discovery of an aminoglycoside-binding riboswitch that is widely distributed among antibiotic-resistant bacterial pathogens. This riboswitch is present in the leader RNA of the resistance genes that encode the aminoglycoside acetyl transferase (AAC) and aminoglycoside adenyl transferase (AAD) enzymes that confer resistance to aminoglycoside antibiotics through modification of the drugs. We show that expression of the AAC and AAD resistance genes is regulated by aminoglycoside binding to a secondary structure in their 5' leader RNA. Reporter gene expression, direct measurements of drug RNA binding, chemical probing, and UV crosslinking combined with mutational analysis demonstrate that the leader RNA functions as an aminoglycoside-sensing riboswitch in which drug binding to the leader RNA leads to the induction of aminoglycosides antibiotic resistance.

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.

Kinetics of tRNA folding monitored by aminoacylation

We describe a strategy for tracking Mg²⁺-initiated folding of ³²P-labeled tRNA molecules to their native structures based on the capacity for aminoacylation by the cognate aminoacyl-tRNA synthetase enzyme. The approach directly links folding to function, paralleling a common strategy used to study the folding of catalytic RNAs. Incubation of unfolded tRNA with magnesium ions, followed by the addition of aminoacyl-tRNA synthetase and further incubation, yields a rapid burst of aminoacyl-tRNA formation corresponding to the prefolded tRNA fraction. A subsequent slower increase in product formation monitors continued folding in the presence of the enzyme. Further analysis reveals the presence of a parallel fraction of tRNA that folds more rapidly than the majority of the population. The application of the approach to study the influence of post-transcriptional modifications in folding of Escherichia coli tRNA₁(Gln) reveals that the modified bases increase the folding rate but do not affect either the equilibrium between properly folded and misfolded states or the folding pathway. This assay allows the use of ³²P-labeled tRNA in integrated studies combining folding, post-transcriptional processing, and aminoacylation reactions.

Exploring RNA structural codes with SHAPE chemistry

RNA is the central conduit for gene expression. This role depends on an ability to encode information at two levels: in its linear sequence and in the complex structures RNA can form by folding back on itself. Understanding the global structure-function interrelationships mediated by RNA remains a great challenge in molecular and structural biology. In this Account, we discuss evolving work in our laboratory focused on creating facile, generic, quantitative, accurate, and highly informative approaches for understanding RNA structure in biologically important environments. The core innovation derives from our discovery that the nucleophilic reactivity of the ribose 2'-hydroxyl in RNA is gated by local nucleotide flexibility. The 2'-hydroxyl is reactive at conformationally flexible positions but is unreactive at nucleotides constrained by base pairing. Sites of modification in RNA can be detected efficiently either using primer extension or by protection from exoribonucleolytic degradation. This technology is now called SHAPE, for selective 2'-hydroxyl acylation analyzed by primer extension (or protection from exoribonuclease). SHAPE reactivities are largely independent of nucleotide identity but correlate closely with model-free measurements of molecular order. The simple SHAPE reaction is thus a robust, nucleotide-resolution, biophysical measurement of RNA structure. SHAPE can be used to provide an experimental correction to RNA folding algorithms and, in favorable cases, yield kilobase-scale secondary structure predictions with high accuracies. SHAPE chemistry is based on very simple reactive carbonyl centers that can be varied to yield slow- and fast-reacting reagents. Differential SHAPE reactivities can be used to detect specific RNA positions with slow local nucleotide dynamics. These positions, which are often in the C2'-endo conformation, have the potential to function as molecular timers that regulate RNA folding and function. In addition, fast-reacting SHAPE reagents can be used to visualize RNA structural biogenesis and RNA-protein assembly reactions in one second snapshots in very straightforward experiments. The application of SHAPE to challenging problems in biology has revealed surprises in well-studied systems. New regions have been identified that are likely to have critical functional roles on the basis of their high levels of RNA structure. For example, SHAPE analysis of large RNAs, such as authentic viral RNA genomes, suggests that RNA structure organizes regulatory motifs and regulates splicing, protein folding, genome recombination, and ribonucleoprotein assembly. SHAPE has also revealed limitations to the hierarchical model for RNA folding. Continued development and application of SHAPE technologies will advance our understanding of the many ways in which the genetic code is expressed through the underlying structure of RNA.

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.

Nonhierarchical ribonucleoprotein assembly suggests a strain-propagation model for protein-facilitated RNA folding

Proteins play diverse and critical roles in cellular ribonucleoproteins (RNPs) including promoting formation of and stabilizing active RNA conformations. Yet, the conformational changes required to convert large RNAs into active RNPs have proven difficult to characterize fully. Here we use high-resolution approaches to monitor both local nucleotide flexibility and solvent accessibility for nearly all nucleotides in the bI3 group I intron RNP in four assembly states: the free RNA, maturase-bound RNA, Mrs1-bound RNA, and the complete six-component holocomplex. The free RNA is misfolded relative to the secondary structure required for splicing. The maturase and Mrs1 proteins each stabilized long-range tertiary interactions, but neither protein alone induced folding into the functional secondary structure. In contrast, simultaneous binding by both proteins results in large secondary structure rearrangements in the RNA and yielded the catalytically active group I intron structure. Secondary and tertiary folding of the RNA component of the bI3 RNP are thus not independent: RNA folding is strongly nonhierarchical. These results emphasize that protein-mediated stabilization of RNA tertiary interactions functions to pull the secondary structure into an energetically disfavored, but functional, conformation and emphasize a new role for facilitator proteins in RNP assembly.

SHAPE-directed RNA secondary structure prediction

The diverse functional roles of RNA are determined by its underlying structure. Accurate and comprehensive knowledge of RNA structure would inform a broader understanding of RNA biology and facilitate exploiting RNA as a biotechnological tool and therapeutic target. Determining the pattern of base pairing, or secondary structure, of RNA is a first step in these endeavors. Advances in experimental, computational, and comparative analysis approaches for analyzing secondary structure have yielded accurate structures for many small RNAs, but only a few large (>500 nts) RNAs. In addition, most current methods for determining a secondary structure require considerable effort, analytical expertise, and technical ingenuity. In this review, we outline an efficient strategy for developing accurate secondary structure models for RNAs of arbitrary length. This approach melds structural information obtained using SHAPE chemistry with structure prediction using nearest-neighbor rules and the dynamic programming algorithm implemented in the RNAstructure program. Prediction accuracies reach >or=95% for RNAs on the kilobase scale. This approach facilitates both development of new models and refinement of existing RNA structure models, which we illustrate using the Gag-Pol frameshift element in an HIV-1 M-group genome. Most promisingly, integrated experimental and computational refinement brings closer the ultimate goal of efficiently and accurately establishing the secondary structure for any RNA sequence.

Advances in RNA structure analysis by chemical probing

RNA is arguably the most versatile biological macromolecule because of its ability both to encode and to manipulate genetic information. The diverse roles of RNA depend on its ability to fold back on itself to form biologically functional structures that bind small molecule and large protein ligands, to change conformation, and to affect the cellular regulatory state. These features of RNA biology can be structurally interrogated using chemical mapping experiments. The usefulness and applications of RNA chemical probing technologies have expanded dramatically over the past five years because of several critical advances. These innovations include new sequence-independent RNA chemistries, algorithmic tools for high-throughput analysis of complex data sets composed of thousands of measurements, new approaches for interpreting chemical probing data for both secondary and tertiary structure prediction, facile methods for following time-dependent processes, and the willingness of individual research groups to tackle increasingly bold problems in RNA structural biology.

Experimental and computational determination of tRNA dynamics

As the molecular representation of the genetic code, tRNA plays a central role in the translational machinery where it interacts with several proteins and other RNAs during the course of protein synthesis. These interactions exploit the dynamic flexibility of tRNA. In this minireview, we discuss the effects of modified bases, ions, and proteins on tRNA structure and dynamics and the challenges of observing its motions over the cycle of translation.

Time-resolved RNA SHAPE chemistry: quantitative RNA structure analysis in one-second snapshots and at single-nucleotide resolution

RNA selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry exploits the discovery that conformationally dynamic nucleotides preferentially adopt configurations that facilitate reaction between the 2'-OH group and a hydroxyl-selective electrophile, such as benzoyl cyanide (BzCN), to form a 2'-O-adduct. BzCN is ideally suited for quantitative, time-resolved analysis of RNA folding and ribonucleoprotein (RNP) assembly mechanisms because this reagent both reacts with flexible RNA nucleotides and also undergoes auto-inactivating hydrolysis with a half-life of 0.25 s at 37 degrees C. RNA folding is initiated by addition of Mg(2+) or protein, or other change in solution conditions, and nucleotide resolution structural images are obtained by adding aliquots of the evolving reaction to BzCN and then 'waiting' for 1 second. Sites of the 2'-O-adduct formation are subsequently scored as stops to primer extension using reverse transcriptase. This time-resolved SHAPE protocol makes it possible to obtain 1-second structural snapshots in time-resolved kinetic studies for RNAs of arbitrary length and complexity in a straightforward and concise experiment.

Direct Fe2+ sensing by iron-responsive messenger RNA:repressor complexes weakens binding

Fe(2+) is now shown to weaken binding between ferritin and mitochondrial aconitase messenger RNA noncoding regulatory structures ((iron-responsive element) (IRE)-RNAs) and the regulatory proteins (IRPs), which adds a direct role of iron to regulation that can complement the well known regulatory protein modification and degradative pathways related to iron-induced mRNA translation. We observe that the K(d) value increases 17-fold in 5'-untranslated region IRE-RNA:repressor complexes; Fe(2+), is studied in the absence of O(2). Other metal ions, Mn(2+) and Mg(2+) have similar effects to Fe(2+) but the required Mg(2+) concentration is 100 times greater than for Fe(2+) or Mn(2+). Metal ions also weaken ethidium bromide binding to IRE-RNA with no effect on IRP fluorescence, using Mn(2+) as an O(2)-resistant surrogate for Fe(2+), indicating that metal ions bound IRE-RNA but not IRP: Fe(2+) decreases IRP repressor complex stability of ferritin IRE-RNA 5-10 times compared with 2-5 times for mitochondrial aconitase IRE-RNA, over the same concentration range, suggesting that differences among IRE-RNA structures contribute to the differences in the iron responses observed in vivo. The results show the IRE-RNA:repressor complex literally responds to Fe(2+), selectively for each IRE-mRNA.

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.

Native-like RNA tertiary structures using a sequence-encoded cleavage agent and refinement by discrete molecular dynamics

The difficulty of analyzing higher order RNA structure, especially for folding intermediates and for RNAs whose functions require domains that are conformationally flexible, emphasizes the need for new approaches for modeling RNA tertiary structure accurately. Here, we report a concise approach that makes use of facile RNA structure probing experiments that are then interpreted using a computational algorithm, carefully tailored to optimize both the resolution and refinement speed for the resulting structures, without requiring user intervention. The RNA secondary structure is first established using SHAPE chemistry. We then use a sequence-directed cleavage agent, which can be placed arbitrarily in many helical motifs, to obtain high quality inter-residue distances. We interpret this in-solution chemical information using a fast, coarse grained, discrete molecular dynamics engine in which each RNA nucleotide is represented by pseudoatoms for the phosphate, ribose, and nucleobase groups. By this approach, we refine base paired positions in yeast tRNA(Asp) to 4 A rmsd without any preexisting information or assumptions about secondary or tertiary structures. This blended experimental and computational approach has the potential to yield native-like models for the diverse universe of functionally important RNAs whose structures cannot be characterized by conventional structural methods.

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.

High-throughput SHAPE and hydroxyl radical analysis of RNA structure and ribonucleoprotein assembly

RNA folds to form complex structures vital to many cellular functions. Proteins facilitate RNA folding at both the secondary and tertiary structure levels. An absolute prerequisite for understanding RNA folding and ribonucleoprotein (RNP) assembly reactions is a complete understanding of the RNA structure at each stage of the folding or assembly process. Here we provide a guide for comprehensive and high-throughput analysis of RNA secondary and tertiary structure using SHAPE and hydroxyl radical footprinting. As an example of the strong and sometimes surprising conclusions that can emerge from high-throughput analysis of RNA folding and RNP assembly, we summarize the structure of the bI3 group I intron RNA in four distinct states. Dramatic structural rearrangements occur in both secondary and tertiary structure as the RNA folds from the free state to the active, six-component, RNP complex. As high-throughput and high-resolution approaches are applied broadly to large protein-RNA complexes, other proteins previously viewed as making simple contributions to RNA folding are also likely to be found to exert multifaceted, long-range, cooperative, and nonadditive effects on RNA folding. These protein-induced contributions add another level of control, and potential regulatory function, in RNP complexes.

Accurate SHAPE-directed RNA structure determination

Almost all RNAs can fold to form extensive base-paired secondary structures. Many of these structures then modulate numerous fundamental elements of gene expression. Deducing these structure-function relationships requires that it be possible to predict RNA secondary structures accurately. However, RNA secondary structure prediction for large RNAs, such that a single predicted structure for a single sequence reliably represents the correct structure, has remained an unsolved problem. Here, we demonstrate that quantitative, nucleotide-resolution information from a SHAPE experiment can be interpreted as a pseudo-free energy change term and used to determine RNA secondary structure with high accuracy. Free energy minimization, by using SHAPE pseudo-free energies, in conjunction with nearest neighbor parameters, predicts the secondary structure of deproteinized Escherichia coli 16S rRNA (>1,300 nt) and a set of smaller RNAs (75-155 nt) with accuracies of up to 96-100%, which are comparable to the best accuracies achievable by comparative sequence analysis.

The UA_handle: a versatile submotif in stable RNA architectures

Stable RNAs are modular and hierarchical 3D architectures taking advantage of recurrent structural motifs to form extensive non-covalent tertiary interactions. Sequence and atomic structure analysis has revealed a novel submotif involving a minimal set of five nucleotides, termed the UA_handle motif (5'XU/AN(n)X3'). It consists of a U:A Watson-Crick: Hoogsteen trans base pair stacked over a classic Watson-Crick base pair, and a bulge of one or more nucleotides that can act as a handle for making different types of long-range interactions. This motif is one of the most versatile building blocks identified in stable RNAs. It enters into the composition of numerous recurrent motifs of greater structural complexity such as the T-loop, the 11-nt receptor, the UAA/GAN and the G-ribo motifs. Several structural principles pertaining to RNA motifs are derived from our analysis. A limited set of basic submotifs can account for the formation of most structural motifs uncovered in ribosomal and stable RNAs. Structural motifs can act as structural scaffoldings and be functionally and topologically equivalent despite sequence and structural differences. The sequence network resulting from the structural relationships shared by these RNA motifs can be used as a proto-language for assisting prediction and rational design of RNA tertiary structures.

ShapeFinder: a software system for high-throughput quantitative analysis of nucleic acid reactivity information resolved by capillary electrophoresis

Analysis of the long-range architecture of RNA is a challenging experimental and computational problem. Local nucleotide flexibility, which directly reports underlying base pairing and tertiary interactions in an RNA, can be comprehensively assessed at single nucleotide resolution using high-throughput selective 2'-hydroxyl acylation analyzed by primer extension (hSHAPE). hSHAPE resolves structure-sensitive chemical modification information by high-resolution capillary electrophoresis and typically yields quantitative nucleotide flexibility information for 300-650 nucleotides (nt) per experiment. The electropherograms generated in hSHAPE experiments provide a wealth of structural information; however, significant algorithmic analysis steps are required to generate quantitative and interpretable data. We have developed a set of software tools called ShapeFinder to make possible rapid analysis of raw sequencer data from hSHAPE, and most other classes of nucleic acid reactivity experiments. The algorithms in ShapeFinder (1) convert measured fluorescence intensity to quantitative cDNA fragment amounts, (2) correct for signal decay over read lengths extending to 600 nts or more, (3) align reactivity data to the known RNA sequence, and (4) quantify per nucleotide reactivities using whole-channel Gaussian integration. The algorithms and user interface tools implemented in ShapeFinder create new opportunities for tackling ambitious problems involving high-throughput analysis of structure-function relationships in large RNAs.

Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA

The functions of most RNA molecules are critically dependent on the distinct local dynamics that characterize secondary structure and tertiary interactions and on structural changes that occur upon binding by proteins and small molecule ligands. Measurements of RNA dynamics at nucleotide resolution set the foundation for understanding the roles of individual residues in folding, catalysis, and ligand recognition. In favorable cases, local order in small RNAs can be quantitatively analyzed by NMR in terms of a generalized order parameter, S2. Alternatively, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry measures local nucleotide flexibility in RNAs of any size using structure-sensitive reagents that acylate the 2'-hydroxyl position. In this work, we compare per-residue RNA dynamics, analyzed by both S2 and SHAPE, for three RNAs: the HIV-1 TAR element, the U1A protein binding site, and the Tetrahymena telomerase stem loop 4. We find a very strong correlation between the two measurements: nucleotides with high SHAPE reactivities consistently have low S2 values. We conclude that SHAPE chemistry quantitatively reports local nucleotide dynamics and can be used with confidence to analyze dynamics in large RNAs, RNA-protein complexes, and RNAs in vivo.