RNA structure has critical roles in processes ranging from ligand sensing to the regulation of translation, polyadenylation and splicing. However, a lack of genome-wide in vivo RNA structural data has limited our understanding of how RNA structure regulates gene expression in living cells. Here we present a high-throughput, genome-wide in vivo RNA structure probing method, structure-seq, in which dimethyl sulphate methylation of unprotected adenines and cytosines is identified by next-generation sequencing. Application of this method to Arabidopsis thaliana seedlings yielded the first in vivo genome-wide RNA structure map at nucleotide resolution for any organism, with quantitative structural information across more than 10,000 transcripts. Our analysis reveals a three-nucleotide periodic repeat pattern in the structure of coding regions, as well as a less-structured region immediately upstream of the start codon, and shows that these features are strongly correlated with translation efficiency. We also find patterns of strong and weak secondary structure at sites of alternative polyadenylation, as well as strong secondary structure at 5' splice sites that correlates with unspliced events. Notably, in vivo structures of messenger RNAs annotated for stress responses are poorly predicted in silico, whereas mRNA structures of genes related to cell function maintenance are well predicted. Global comparison of several structural features between these two categories shows that the mRNAs associated with stress responses tend to have more single-strandedness, longer maximal loop length and higher free energy per nucleotide, features that may allow these RNAs to undergo conformational changes in response to environmental conditions. Structure-seq allows the RNA structurome and its biological roles to be interrogated on a genome-wide scale and should be applicable to any organism.

Many protein enzymes use general acid-base catalysis as a way to increase reaction rates. The amino acid histidine is optimized for this function because it has a pK(a) (where K(a) is the acid dissociation constant) near physiological pH. The RNA enzyme (ribozyme) from hepatitis delta virus catalyzes self-cleavage of a phosphodiester bond. Reactivity-pH profiles in monovalent or divalent cations, as well as distance to the leaving-group oxygen, implicate cytosine 75 (C75) of the ribozyme as the general acid and ribozyme-bound hydrated metal hydroxide as the general base in the self-cleavage reaction. Moreover, C75 has a pK(a) perturbed to neutrality, making it "histidine-like." Anticooperative interaction is observed between protonated C75 and a metal ion, which serves to modulate the pK(a) of C75. General acid-base catalysis expands the catalytic repertoire of RNA and may provide improved rate acceleration.

The hepatitis delta virus (HDV) ribozyme and HDV-like ribozymes are self-cleaving RNAs found throughout all kingdoms of life. These RNAs fold into a double-nested pseudoknot structure and cleave RNA, yielding 2',3'-cyclic phosphate and 5'-hydroxyl termini. The active site nucleotide C75 has a pK(a) shifted >2 pH units toward neutrality and has been implicated as a general acid/base in the cleavage reaction. An active site Mg(2+) ion that helps activate the 2'-hydroxyl for nucleophilic attack has been characterized biochemically; however, this ion has not been visualized in any previous structures. To create a snapshot of the ribozyme in a state poised for catalysis, we have crystallized and determined the structure of the HDV ribozyme bound to an inhibitor RNA containing a deoxynucleotide at the cleavage site. This structure includes the wild-type C75 nucleotide and Mg(2+) ions, both of which are required for maximal ribozyme activity. This structure suggests that the position of C75 does not change during the cleavage reaction. A partially hydrated Mg(2+) ion is also found within the active site where it interacts with a newly resolved G.U reverse wobble. Although the inhibitor exhibits crystallographic disorder, we modeled the ribozyme-substrate complex using the conformation of the inhibitor strand observed in the hammerhead ribozyme. This model suggests that the pro-R(P) oxygen of the scissile phosphate and the 2'-hydroxyl nucleophile are inner-sphere ligands to the active site Mg(2+) ion. Thus, the HDV ribozyme may use a combination of metal ion Lewis acid and nucleobase general acid strategies to effect RNA cleavage.

Previous studies have shown that the translation level of in vitro transcribed messenger RNA (mRNA) is enhanced when its uridines are replaced with pseudouridines; however, the reason for this enhancement has not been identified. Here, we demonstrate that in vitro transcripts containing uridine activate RNA-dependent protein kinase (PKR), which then phosphorylates translation initiation factor 2-alpha (eIF-2α), and inhibits translation. In contrast, in vitro transcribed mRNAs containing pseudouridine activate PKR to a lesser degree, and translation of pseudouridine-containing mRNAs is not repressed. RNA pull-down assays demonstrate that mRNA containing uridine is bound by PKR more efficiently than mRNA with pseudouridine. Finally, the role of PKR is validated by showing that pseudouridine- and uridine-containing RNAs were translated equally in PKR knockout cells. These results indicate that the enhanced translation of mRNAs containing pseudouridine, compared to those containing uridine, is mediated by decreased activation of PKR.

The human double-stranded RNA- (dsRNA) activated protein kinase (PKR) has a dsRNA-binding domain (dsRBD) that contains two tandem copies of the dsRNA-binding motif (dsRBM). The minimal-length polypeptide required to bind dsRNA contains both dsRBMs, as determined by mobility-shift and filter-binding assays. Mobility-shift experiments indicate binding requires a minimum of 16 base pairs of dsRNA, while a minimal-length site for saturation of longer RNAs is 11 base pairs. Bulge defects in the helix disfavor binding, and single-stranded tails do not strongly influence the dsRNA length requirement. These polypeptides do not bind an RNA-DNA hybrid duplex or dsDNA as judged by either mobility-shift or competition experiments, suggesting 2'-OH contacts on both strands of the duplex stabilize binding. Related experiments on chimeric duplexes in which specific sets of 2'-OHs are substituted with 2'-H or 2'-OCH3 reveal that the 2'-OHs required for binding are located along the entire 11 basepair site. These results are supported by Fe(II) EDTA footprinting experiments that show protein-dependent protection of the minor groove of dsRNA. The dependence of dsRNA-protein binding on salt concentration suggests that only one ionic contact is made between the protein and dsRNA phosphate backbone and that at physiological salt concentrations 90% of the free energy of binding is nonelectrostatic. Thus, the specificity of PKR for dsRNA over RNA-DNA hybrids and dsDNA is largely due to molecular recognition of a network of 2'-OHs involving both strands of dsRNA and present along the entire 11 base-pair site.