RNase J1 endonuclease activity as a probe of RNA secondary structure

Reliable determination of RNA secondary structure depends on both computer algorithms and experimental probing of nucleotides in single- or double-stranded conformation. Here we describe the exploitation of the endonucleolytic activity of the Bacillus subtilis enzyme RNase J1 as a probe of RNA structure. RNase J1 cleaves in single-stranded regions and, in vitro at least, the enzyme has relatively relaxed nucleotide specificity. We confirmed the feasibility of the approach on an RNA of known structure, B. subtilis tRNA^Thr. We then used RNase J1 to solve the secondary structure of the 5′ end of the hbs mRNA. Finally, we showed that RNase J1 can also be used in footprinting experiments by probing the interaction between the 30S ribosomal subunit and the Shine–Dalgarno element of the hbs 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.     

Structural insights into the targeting of mRNA GU-rich elements by the three RRMs of CELF1

The CUG-BP, Elav-like family (CELF) of RNA-binding proteins control gene expression at a number of different levels by regulating pre-mRNA splicing, deadenylation and mRNA stability. We present structural insights into the binding selectivity of CELF member 1 (CELF1) for GU-rich mRNA target sequences of the general form 5′-UGUN_xUGUN_yUGU and identify a high affinity interaction (K_d ∼ 100 nM for x = 2 and y = 4) with simultaneous binding of all three RNA recognition motifs within a single 15-nt binding element. RNA substrates spin-labelled at either the 3′ or 5′ terminus result in differential nuclear magnetic resonance paramagnetic relaxation enhancement effects, which are consistent with a non-sequential 2-1-3 arrangement of the three RNA recognition motifs on UGU sites in a 5′ to 3′ orientation along the RNA target. We further demonstrate that CELF1 binds to dispersed single-stranded UGU sites at the base of an RNA hairpin providing a structural rationale for recognition of CUG expansion repeats and splice site junctions in the regulation of alternative splicing.     

Characterization of a novel eukaryal nick-sealing RNA ligase from Naegleria gruberi

The proteome of the amoebo-flagellate protozoan Naegleria gruberi is rich in candidate RNA repair enzymes, including 15 putative RNA ligases, one of which, NgrRnl, is a eukaryal homolog of Deinococcus radiodurans RNA ligase, DraRnl. Here we report that purified recombinant NgrRnl seals nicked 3′-OH/5′-PO₄ duplexes in which the 3′-OH strand is RNA. It does so via the “classic” ligase pathway, entailing reaction with ATP to form a covalent NgrRnl–AMP intermediate, transfer of AMP to the nick 5′-PO₄, and attack of the RNA 3′-OH on the adenylylated nick to form a 3′–5′ phosphodiester. Unlike members of the four known families of ATP-dependent RNA ligases, NgrRnl lacks a carboxy-terminal appendage to its nucleotidyltransferase domain. Instead, it contains a defining amino-terminal domain that we show is important for 3′-OH/5′-PO₄ nick-sealing and ligase adenylylation, but dispensable for phosphodiester synthesis at a preadenylylated nick. We propose that NgrRnl, DraRnl, and their homologs from diverse bacteria, viruses, and unicellular eukarya comprise a new “Rnl5 family” of nick-sealing ligases with a signature domain organization.     

The 1.19 A X-ray structure of 2'-O-Me(CGCGCG)(2) duplex shows dehydrated RNA with 2-methyl-2,4-pentanediol in the minor groove

The crystal and molecular structure of 2′-O-Me(CGCGCG)₂ has been determined at 1.19 Å resolution, at 100 K, using synchrotron radiation. The structure in space group P3₂12 is a half-turn right-handed helix that includes two 2-methyl-2,4-pentanediol (MPD) molecules bound in the minor groove. The structure deviates from A-form RNA. The duplex is overwound with an average value of 9.7 bp per turn, characterised as having a C3′-endo sugar pucker, very low base pair rise and high helical twist and inclination angles. The structure includes 65 ordered water molecules. Only a single row of water molecules is observed in the minor groove due to the presence of hydrophobic 2′-O-methyl groups. As many as five magnesium ions are located in the structure. Two are in the major groove and interact with O⁶ and N⁷ of guanosine and N⁴ of cytidine residues through their hydration spheres. This work provides the first example of molecular interactions of nucleic acids with MPD, which was used as a precipitant, cryo-solvent and resolution enhancing agent. The two MPD molecules intrude into the hydration network in the minor groove, each forming hydrogen bonds between their secondary hydroxyl group and exo-amino functions of guanosine residues. Comparison of the 2′-O-Me(CGCGCG)₂ structure in the P3₂12 and P6₁22 crystals delineates stability of the water network within the minor groove to dehydration by MPD and is of interest for evaluating factors governing small molecule binding to RNA. Intrusion of MPD into the minor groove of 2′-O-Me(CGCGCG)₂ is discussed with respect to RNA dehydration, a prerequisite of Z-RNA formation.     

Pathology-related substitutions in human mitochondrial tRNA(Ile) reduce precursor 3' end processing efficiency in vitro

The human mitochondrial genome encodes 22 tRNAs interspersed among the two rRNAs and 11 mRNAs, often without spacers, suggesting that tRNAs must be efficiently excised. Numerous maternally transmitted diseases and syndromes arise from mutations in mitochondrial tRNAs, likely due to defect(s) in tRNA metabolism. We have systematically explored the effect of pathogenic mutations on tRNA^Ile precursor 3′ end maturation in vitro by 3′-tRNase. Strikingly, four pathogenic tRNA^Ile mutations reduce 3′-tRNase processing efficiency (V_max / K_M) to ~10-fold below that of wild-type, principally due to lower V_max. The structural impact of mutations was sought by secondary structure probing and wild-type tRNA^Ile precursor was found to fold into a canonical cloverleaf. Among the mutant tRNA^Ile precursors with the greatest 3′ end processing deficiencies, only G4309A displays a secondary structure substantially different from wild-type, with changes in the T domain proximal to the substitution. Reduced efficiency of tRNA^Ile precursor 3′ end processing, in one case associated with structural perturbations, could thus contribute to human mitochondrial diseases caused by mutant tRNAs.     

New scoring system to identify RNA G-quadruplex folding

G-quadruplexes (G4s) are non-canonical structures involved in many important cellular processes. To date, the prediction of potential G-quadruplex structures (PG4s) has been based almost exclusively on the sequence of interest agreeing with the algorithm G_x-N-_1–7-G_x-N_1–7-G_x-N_1–7-G_x (where x ≥ 3 and N = A, U, G or C). However, many sequences agreeing with this algorithm do not form G4s and are considered false-positive predictions. Here we show the RNA PG4 candidate in the 3′-untranslated region (UTR) of the TTYH1 gene to be one such false positive. Specifically, G4 folding was observed to be inhibited by the presence of multiple-cytosine tracks, located in the candidate’s genomic context, that adopted a Watson–Crick base-paired structure. Clearly, the neighbouring sequences of a PG4 may influence its folding. The secondary structure of 12 PG4 motifs along with either 15 or 50 nucleotides of their upstream and downstream genomic contexts were evaluated by in-line probing. Data permitted the development of a scoring system for the prediction of PG4s taking into account the effect of the neighbouring sequences. The accuracy of this scoring system was assessed by probing 14 other novel PG4 candidates retrieved in human 5′-UTRs. This new scoring system can be used, in combination with the standard algorithm, to better predict the folding of RNA G4s.     

The secondary structure of the 5' end of the FIV genome reveals a long-range interaction between R/U5 and gag sequences, and a large, stable stem-loop

Feline immunodeficiency virus (FIV) is a lentivirus that infects cats and is related to human immunodeficiency virus (HIV). Although it is a common worldwide infection, and has potential uses as a human gene therapy vector and as a nonprimate model for HIV infection, little detail is known of the viral life cycle. Previous experiments have shown that its packaging signal includes two or more regions within the first 511 nucleotides of the genomic RNA. We have undertaken a secondary structural analysis of this RNA by minimal free-energy structural prediction, biochemical mapping, and phylogenetic analysis, and show that it contains five conserved stem–loops and a conserved long-range interaction between heptanucleotide sequences 5′-CCCUGUC-3′ in R/U5 and 5′-GACAGGG-3′ in gag. This long-range interaction is similar to that seen in primate lentiviruses where it is thought to be functionally important. Along with strains that infect domestic cats, this heptanucleotide interaction can also occur in species-specific FIV strains that infect pumas, lions, and Pallas' cats where the heptanucleotide sequences involved vary. We have analyzed spliced and genomic FIV RNAs and see little structural change or sequence conservation within single-stranded regions of the 5′ UTR that are important for viral packaging, suggesting that FIV may employ a cotranslational packaging mechanism.     

Structural domains within the 3' untranslated region of Turnip crinkle virus

The genomes of positive-strand RNA viruses undergo conformational shifts that complicate efforts to equate structures with function. We have initiated a detailed analysis of secondary and tertiary elements within the 3′ end of Turnip crinkle virus (TCV) that are required for viral accumulation in vivo. MPGAfold, a massively parallel genetic algorithm, suggested the presence of five hairpins (H4a, H4b, and previously identified hairpins H4, H5, and Pr) and one H-type pseudoknot (Ψ₃) within the 3′-terminal 194 nucleotides (nt). In vivo compensatory mutagenesis analyses confirmed the existence of H4a, H4b, Ψ₃ and a second pseudoknot (Ψ₂) previously identified in a TCV satellite RNA. In-line structure probing of the 194-nt fragment supported the coexistence of H4, H4a, H4b, Ψ₃ and a pseudoknot that connects H5 and the 3′ end (Ψ₁). Stepwise replacements of TCV elements with the comparable elements from Cardamine chlorotic fleck virus indicated that the complete 142-nt 3′ end, and subsets containing Ψ₃, H4a, and H4b or Ψ₃, H4a, H4b, H5, and Ψ₂, form functional domains for virus accumulation in vivo. A new 3-D molecular modeling protocol (RNA2D3D) predicted that H4a, H4b, H5, Ψ₃, and Ψ₂ are capable of simultaneous existence and bears some resemblance to a tRNA. The related Japanese iris necrotic ring virus does not have comparable domains. These results provide a framework for determining how interconnected elements participate in processes that require 3′ untranslated region sequences such as translation and replication.     

High salt solution structure of a left-handed RNA double helix

Right-handed RNA duplexes of (CG)_n sequence undergo salt-induced helicity reversal, forming left-handed RNA double helices (Z-RNA). In contrast to the thoroughly studied Z-DNA, no Z-RNA structure of natural origin is known. Here we report the NMR structure of a half-turn, left-handed RNA helix (CGCGCG)₂ determined in 6 M NaClO₄. This is the first nucleic acid motif determined at such high salt. Sequential assignments of non-exchangeable proton resonances of the Z-form were based on the hitherto unreported NOE connectivity path [H6_(n)-H5′/H5″_(n)-H8_(n+1)-H1′_(n+1)-H6_(n+2)] found for left-handed helices. Z-RNA structure shows several conformational features significantly different from Z-DNA. Intra-strand but no inter-strand base stacking was observed for both CpG and GpC steps. Helical twist angles for CpG steps have small positive values (4–7°), whereas GpC steps have large negative values (−61°). In the full-turn model of Z-RNA (12.4 bp per turn), base pairs are much closer to the helix axis than in Z-DNA, thus both the very deep, narrow minor groove with buried cytidine 2′-OH groups, and the major groove are well defined. The 2′-OH group of cytidines plays a crucial role in the Z-RNA structure and its formation; 2′-O-methylation of cytidine, but not of guanosine residues prohibits A to Z helicity reversal.     

CvfA Protein and Polynucleotide Phosphorylase Act in an Opposing Manner to Regulate Staphylococcus aureus Virulence

We previously identified CvfA (SA1129) as a Staphylococcus aureus virulence factor using a silkworm infection model. S. aureus cvfA-deleted mutants exhibit decreased expression of the agr locus encoding a positive regulator of hemolysin genes and decreased hemolysin production. CvfA protein hydrolyzes a 2′,3′-cyclic phosphodiester bond at the RNA 3′ terminus, producing RNA with a 3′-phosphate (3′-phosphorylated RNA, RNA with a 3′-phosphate). Here, we report that the cvfA-deleted mutant phenotype (decreased agr expression and hemolysin production) was suppressed by disrupting pnpA-encoding polynucleotide phosphorylase (PNPase) with 3′- to 5′-exonuclease activity. The suppression was blocked by introducing a pnpA-encoding PNPase with exonuclease activity but not by a pnpA-encoding mutant PNPase without exonuclease activity. Therefore, loss of PNPase exonuclease activity suppressed the cvfA-deleted mutant phenotype. Purified PNPase efficiently degraded RNA with 2′,3′-cyclic phosphate at the 3′ terminus (2′,3′-cyclic RNA), but it inefficiently degraded 3′-phosphorylated RNA. These findings indicate that 3′-phosphorylated RNA production from 2′,3′-cyclic RNA by CvfA prevents RNA degradation by PNPase and contributes to the expression of agr and hemolysin genes. We speculate that in the cvfA-deleted mutant, 2′,3′-cyclic RNA is not converted to the 3′-phosphorylated form and is efficiently degraded by PNPase, resulting in the loss of RNA essential for expressing agr and hemolysin genes, whereas in the cvfA/pnpA double-disrupted mutant, 2′,3′-cyclic RNA is not degraded by PNPase, leading to hemolysin production. These findings suggest that CvfA and PNPase competitively regulate RNA degradation essential for S. aureus virulence.     

Specificity and Evolutionary Conservation of the Escherichia coli RNA Pyrophosphohydrolase RppH

Bacterial RNA degradation often begins with conversion of the 5′-terminal triphosphate to a monophosphate by the RNA pyrophosphohydrolase RppH, an event that triggers rapid ribonucleolytic attack. Besides its role as the master regulator of 5′-end-dependent mRNA decay, RppH is important for the ability of pathogenic bacteria to invade host cells, yet little is known about how it chooses its targets. Here, we show that Escherichia coli RppH (EcRppH) requires at least two unpaired nucleotides at the RNA 5′ end and prefers three or more such nucleotides. It can tolerate any nucleotide at the first three positions but has a modest preference for A at the 5′ terminus and either a G or A at the second position. Mutational analysis has identified EcRppH residues crucial for substrate recognition or catalysis. The promiscuity of EcRppH differentiates it from its Bacillus subtilis counterpart, which has a strict RNA sequence requirement. EcRppH orthologs likely to share its relaxed sequence specificity are widespread in all classes of Proteobacteria, except Deltaproteobacteria, and in flowering plants. By contrast, the phylogenetic range of recognizable B. subtilis RppH orthologs appears to be restricted to the order Bacillales. These findings help to explain the selective influence of RppH on bacterial mRNA decay and show that RppH-dependent degradation has diversified significantly during the course of evolution.     

Biochemical analysis of pistol self-cleaving ribozymes

Pistol RNAs are members of a distinct class of self-cleaving ribozymes that was recently discovered by using a bioinformatics search strategy. Several hundred pistol ribozymes share a consensus sequence including 10 highly conserved nucleotides and many other modestly conserved nucleotides associated with specific secondary structure features, including three base-paired stems and a pseudoknot. A representative pistol ribozyme from the bacterium Lysinibacillus sphaericus was found to promote RNA strand scission with a rate constant of ∼10 min⁻¹ under physiological Mg²⁺ and pH conditions. The reaction proceeds via the nucleophilic attack of a 2′-oxygen atom on the adjacent phosphorus center, and thus adheres to the same general catalytic mechanism of internal phosphoester transfer as found with all other classes of natural self-cleaving ribozymes discovered to date. Analyses of the kinetic characteristics and the metal ion requirements of the cleavage reaction reveal that members of this ribozyme class likely use several catalytic strategies to promote the rapid cleavage of RNA.     

Imp3 unfolds stem structures in pre-rRNA and U3 snoRNA to form a duplex essential for small subunit processing

Eukaryotic ribosome biogenesis requires rapid hybridization between the U3 snoRNA and the pre-rRNA to direct cleavages at the A₀, A₁, and A₂ sites in pre-rRNA that liberate the small subunit precursor. The bases involved in hybridization of one of the three duplexes that U3 makes with pre-rRNA, designated the U3-18S duplex, are buried in conserved structures: box A/A′ stem–loop in U3 snoRNA and helix 1 (H1) in the 18S region of the pre-rRNA. These conserved structures must be unfolded to permit the necessary hybridization. Previously, we reported that Imp3 and Imp4 promote U3-18S hybridization in vitro, but the mechanism by which these proteins facilitate U3-18S duplex formation remained unclear. Here, we directly addressed this question by probing base accessibility with chemical modification and backbone accessibility with ribonuclease activity of U3 and pre-rRNA fragments that mimic the secondary structure observed in vivo. Our results demonstrate that U3-18S hybridization requires only Imp3. Binding to each RNA by Imp3 provides sufficient energy to unfold both the 18S H1 and the U3 box A/A′ stem structures. The Imp3 unfolding activity also increases accessibility at the U3-dependent A₀ and A₁ sites, perhaps signaling cleavage at these sites to generate the 5′ mature end of 18S. Imp4 destabilizes the U3-18S duplex to aid U3 release, thus differentiating the roles of these proteins. Protein-dependent unfolding of these structures may serve as a switch to block U3-pre-rRNA interactions until recruitment of Imp3, thereby preventing premature and inaccurate U3-dependent pre-rRNA cleavage and folding events in eukaryotic ribosome biogenesis.     

Identification of cis-Acting Nucleotides and a Structural Feature in West Nile Virus 3′-Terminus RNA That Facilitate Viral Minus Strand RNA Synthesis

The 3′-terminal nucleotides (nt) of West Nile virus (WNV) genomic RNA form a penultimate 16-nt small stem-loop (SSL) and an 80-nt terminal stem-loop (SL). These RNA structures are conserved in divergent flavivirus genomes. A previous in vitro study using truncated WNV 3′ RNA structures predicted a putative tertiary interaction between the 5′ side of the 3′-terminal SL and the loop of the SSL. Although substitution or deletion of the 3′ G (nt 87) within the SSL loop, which forms the only G-C pair in the predicted tertiary interaction, in a WNV infectious clone was lethal, a finding consistent with the involvement in a functionally relevant pseudoknot interaction, extensive mutagenesis of nucleotides in the terminal SL did not identify a cis-acting pairing partner for this SSL 3′ G. However, both the sequence and the structural context of two adjacent base pairs flanked by symmetrical internal loops in the 3′-terminal SL were shown to be required for efficient viral RNA replication. Nuclear magnetic resonance analysis confirmed the predicted SSL and SL structures but not the tertiary interaction. The SSL was previously reported to contain one of three eEF1A binding sites, and G87 in the SSL loop was shown to be involved in eEF1A binding. The nucleotides at the bottom part of the 3′-terminal SL switch between 3′ RNA-RNA and 3′-5′ RNA-RNA interactions. The data suggest that interaction of the 3′ SL RNA with eEF1A at three sites and a unique metastable structural feature may participate in regulating structural changes in the 3′-terminal SL.     

Crystal Structure of the Dengue Virus Methyltransferase Bound to a 5′-Capped Octameric RNA

The N-terminal domain of the flavivirus NS5 protein functions as a methyltransferase (MTase). It sequentially methylates the N7 and 2′-O positions of the viral RNA cap structure (GpppA→^7meGpppA→^7meGpppA_2′-O-me). The same NS5 domain could also have a guanylyltransferase activity (GTP+ppA-RNA→GpppA). The mechanism by which this protein domain catalyzes these three distinct functions is currently unknown. Here we report the crystallographic structure of DENV-3 MTase in complex with a 5′-capped RNA octamer (G_pppAGAACCUG) at a resolution of 2.9 Å. Two RNA octamers arranged as kissing loops are encircled by four MTase monomers around a 2-fold non-crystallography symmetry axis. Only two of the four monomers make direct contact with the 5′ end of RNA. The RNA structure is stabilised by the formation of several intra and intermolecular base stacking and non-canonical base pairs. The structure may represent the product of guanylylation of the viral genome prior to the subsequent methylation events that require repositioning of the RNA substrate to reach to the methyl-donor sites. The crystal structure provides a structural explanation for the observed trans-complementation of MTases with different methylation defects.