Recent developments in methods for RNA sequencing using in vitro 32P-labeling

A variety of approaches that utilize in vitro 32P-labeling of RNA and of oligonucleotides in the sequence analysis of RNAs are described. These include 1) methods for 5'- and 3'- end labeling of RNAs; 2) end labeling and sequencing of oligonucleotides present in complete T1 RNase or pancreatic RNase digests of RNA; 3) use of random endonucleases, such as nuclease P1, for terminal sequence analysis of end labeled RNAs; and 4) use of base specific enzymes or chemical reagents in the sequence analysis of end-labeled RNAs. Also described is an approach to RNA sequencing, applied so far to tRNAs, which is based on partial and random alkaline cleavage of an RNA to generate a series of overlapping oligonucleotide fragments, all containing the original 3'-end of the RNA. Analysis of the 5'- end group of each of these oligonucleotides (following 5'-end labeling with 32P) provides the sequence of most of the tRNA. The above methods have been used to derive the sequences of several tRNAs, the ribosomal 5S and 5 x 8S RNAs, a viroid RNA, and large segments of both prokaryotic and eukaryotic ribosomal and messenger RNAs.     

The relationship of thermodynamic stability at a G x U recognition site to tRNA aminoacylation specificity

The G x U pair at the third position in the acceptor helix of Escherichia coli tRNA(Ala) is critical for aminoacylation. The features that allow G x U recognition are likely to include direct interaction of alanyl-tRNA synthetase with distinctive atomic groups and indirect recognition of the structural and stability information encoded in the sequence of G x U and its immediate context. The present work investigates the thermodynamic stability and acceptor activity for a comprehensive set of variant RNAs with substitutions of the G x U pair of E. coli tRNA(Ala). The four RNAs with Watson-Crick substitutions had a lower acceptor activity and a higher stability relative to the G x U RNA. On the other hand, the RNAs with mispair substitutions had a lower stability, but either a higher or a lower acceptor activity. Thus, the entire set of variant RNAs does not exhibit a correlation between thermodynamic stability of the free, unbound tRNA and its acceptor activity. The substantial acceptor activity of tRNAs with particular mispair substitutions may be explained by their ability to assume the conformational preferences of alanyl-tRNA synthetase. Moreover, the G x U pair may provide a point of deformability for the substrate tRNA to adapt to the enzyme's active site.     

Substrate recognition by a yeast 2'-phosphotransferase involved in tRNA splicing and by its Escherichia coli homolog.

The final step of tRNA splicing in Saccharomyces cerevisiae requires 2'-phosphotransferase (Tpt1) to transfer the 2'-phosphate from ligated tRNA to NAD, producing mature tRNA and ADP ribose-1' '-2' '-cyclic phosphate. To address how Tpt1 protein recognizes substrate RNAs, we measured the steady-state kinetic parameters of Tpt1 protein with 2'-phosphorylated ligated tRNA and a variety of related substrates. Tpt1 protein has a high apparent affinity for ligated tRNA (K(m,RNA), 0.35 nM) and a low turnover rate (k(cat), 0.3 min(-1)). Tpt1 protein recognizes both tRNA and the internal 2'-phosphate of RNAs. Steady-state kinetic analysis reveals that as RNAs lose structure and length, K(m,RNA) and k(cat) both increase commensurately. For a 2'-phosphorylated octadecamer derived from the anticodon stem-loop of ligated tRNA, K(m,RNA) and k(cat) are 5- and 8-fold higher, respectively, than for ligated tRNA, whereas for a simple substrate like pApA(p)pA, K(m,RNA) and k(cat) are 430- and 150-fold higher, respectively. Tpt1 is not detectably active on a trimer with a terminal 5'- or 3'-phosphate and is very inefficient at removal of a terminal 2'-phosphate unless there is an adjacent 3'-phosphate or phosphodiester. The K(m,NAD) for Tpt1 is substrate dependent: K(m,NAD) is 10 microM with ligated tRNA, 200 microM with pApA(p)pA, and 600 microM with pApApA(p). Preliminary analysis of KptA, a functional Tpt1 protein homologue from Escherichia coli, reveals that KptA protein is strikingly similar to yeast Tpt1 in its kinetic parameters, although E. coli is not known to have a 2'-phosphorylated RNA substrate.     

Conversion of mammalian tRNA 3' processing endoribonuclease to four-base-recognizing RNA cutters.

The spermidine-dependent, sequence-specific endoribonuclease (RNase 65) activities in mammalian cell extracts require both protein and 3' truncated tRNA, species of which direct their substrate sequence specificity. Computer analysis for searching possible base pairing between substrate RNAs and their corresponding 3' truncated tRNA, suggested a unified model for substrate recognition mechanism, in which a four-nucleotide (nt) sequence in the target tRNAs 1 nt upstream of their cleavage site, base pairs with the 5' terminal 4 nt sequence of their corresponding 3' truncated tRNA. This model was supported by experiments with several RNA substrates containing a substituted nucleotide in the target 4 nt sequence. In this model, the tRNA substrates and their corresponding 3' truncated tRNA form a complex resembling a 5' processed tRNA precursor containing a 3' trailer, suggesting that the protein component of RNase 65 is identical to tRNA 3' processing endoribonuclease (3' tRNase). Actually, 3' tRNase purified from pig liver cleaved the target RNAs at the expected sites only in the presence of their corresponding 3' truncated tRNA. These results show that the 3' tRNase can be converted to 4 nt specific RNA cutters using the 3' truncated tRNAs.     

The tRNA-like structure at the 3' terminus of turnip yellow mosaic virus RNA. Differences and similarities with canonical tRNA.

The 3' terminus of TYMV RNA, which possesses tRNA-like properties, has been studied. A 3' terminal fragment of 112 nucleotides was obtained by cleavage with RNase H after hybridization of a synthetic oligodeoxynucleotide to the viral RNA. The accessibility of cytidine and adenosine residues was probed with chemical modification. Enzymatic digestion studies were performed with RNase T1, nuclease S1 and the double-strand specific RNase from the venom of the cobra Naja naja oxiana. A model is proposed for the secondary structure of the 3' terminal region of TYMV RNA comprising 86 nucleotides. The main feature of this secondary structure is the absence of a conventional acceptor stem as present in canonical tRNA. However, the terminal 42 nucleotides can be folded in a tertiary structure which bears strong resemblance with the acceptor arm of canonical tRNA. Comparison of this region of TYMV RNA with that of other RNAs from both the tymovirus group and the tobamovirus group gives support to our proposal for such a three-dimensional arrangement. The consequences for the recognition by TYMV RNA of tRNA-specific enzymes is discussed.     

Quantum chemical studies of structures and binding in noncanonical RNA base pairs: the trans Watson-Crick:Watson-Crick family

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6-31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from -8.19 kcal/mol to -21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32-2.19 A. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.     

Structural consensus of RNA templates replicated by Q beta replicase.

Q beta replicase replicates a variety of enzyme-specific small RNAs in addition to the phage genomic RNA. The sequence analysis has revealed that all these RNAs are potentially capable of forming a consensus secondary structure element. It represents a stalk which is formed by the 5'-GGG ... and ... CCCA-3' complementary stretches at the termini of the replicating RNA molecules and adjacent 5'- and 3'-hairpins, which may form a stacking with the stalk. The structure found is rather similar to the analogous structure in the tRNA molecule. The genomic RNA of the Q beta phage and other related phages can also form a similar structural element.     

Presence of queuine in Drosophila melanogaster: correlation of free pool with queuosine content of tRNA and effect of mutations in pteridine metabolism.

Queuine, a modified form of 7-deazaguanine present in certain transfer RNAs, is shown to occur in Drosophila melanogaster adults in a free form and its concentration varies as a function of age, nutrition and genotype. In several, but not all mutant strains, the concentrations of queuine and the Q(+) (queuine-containing) form of tRNATyr are correlated. The bioassay employs L-M cells which respond to the presence of queuine by an increase in their Q(+)tRNAAsp that is accompanied by a decrease in the Q(-)tRNAAsp isoacceptors. The increase in Q(+)tRNATyr in Drosophila that occurs on a yeast diet is accompanied by an increase in queuine. Similarly the increase of Q(+)tRNAs with age also is accompanied by an increase in free queuine. In two mutants, brown and sepia, these correlations were either diminished or failed to occur. Indeed, the extract of both mutants inhibited the response of the L-M cells to authentic queuine. When the pteridines that occur at abnormally high levels in sepia were used at 1 x 10(-6)M, the inhibition of the L-M cell assay occurred in the order biopterin greater than pterin greater than sepiapterin. These pteridines were also inhibitory for the purified guanine:tRNA transglycosylase from rabbit but the relative effectiveness then was pterin greater than biopterin greater than sepiapterin. Pterin was competitive with guanine in the enzyme reaction with Ki = 0.9 x 10(-7)M. Also when an extract of sepia was chromatographed on Sephadex G-50, the pteridine-containing fractions only were inhibitory toward the L-M cell assay or the enzyme assay. These results indicate that free queuine occurs in Drosophila but also that certain pteridines may interfere with the incorporation of queuine into RNA.     

The human pseudouridine synthase PUS7 recognizes RNA with an extended multi-domain binding surface

The human pseudouridine synthase PUS7 is a versatile RNA modification enzyme targeting many RNAs thereby playing a critical role in development and brain function. Whereas all target RNAs of PUS7 share a consensus sequence, additional recognition elements are likely required, and the structural basis for RNA binding by PUS7 is unknown. Here, we characterize the structure–function relationship of human PUS7 reporting its X-ray crystal structure at 2.26 Å resolution. Compared to its bacterial homolog, human PUS7 possesses two additional subdomains, and structural modeling studies suggest that these subdomains contribute to tRNA recognition through increased interactions along the tRNA substrate. Consistent with our modeling, we find that all structural elements of tRNA are required for productive interaction with PUS7 as the consensus sequence of target RNA alone is not sufficient for pseudouridylation by human PUS7. Moreover, PUS7 binds several, non-modifiable RNAs with medium affinity which likely enables PUS7 to screen for productive RNA substrates. Following tRNA modification, the product tRNA has a significantly lower affinity for PUS7 facilitating its dissociation. Taken together our studies suggest a combination of structure-specific and sequence-specific RNA recognition by PUS7 and provide mechanistic insight into its function.     

Structure of RNAs replicated by the DNA-dependent T7 RNA polymerase.

The DNA-dependent RNA polymerase of bacteriophage T7 efficiently and specifically replicates two structurally related RNAs, termed X and Y RNAs. Replication of both RNAs involves synthesis of complementary strands initiated with pppC and pppG. RNAs transcribed from DNA template containing the established sequences of X and Y RNAs were efficiently replicated by T7 RNA polymerase. Both RNAs possess palindromic sequences with a dual axis of symmetry, permitting formation of hairpin-, dumbbell-, or cloverleaf-type structures. The template must consist of RNA and not DNA sequence, and the terminal unpaired dinucleotides of the RNA are necessary for replication. Nucleotidyl transferase activity of E. coli adenylates the unpaired CCOH dinucleotide at the 3' end of a C strand of X RNA. This feature, as well as the length (64 nucleotides) and compact structure of X and Y RNAs, suggests that they may resemble tRNA molecules and tRNA-like structures at the 3' termini of many plant viral RNA genomes.     

tRNA-like properties of tobacco rattle virus RNA.

The 3' terminal forty nucleotides of tobraviral RNAs readily fold into a tertiary structure, resembling that of tymo- and tobamoviral RNAs. The latter RNAs possess a tRNA-like structure at their 3' end that is recognized by a number of tRNA-specific enzymes (Rietveld et al. (1984), EMBO J. 3, 2613-2619). Characteristic for their aminoacyl acceptor arm is the presence of a so-called pseudoknot which we now also find in a corresponding position at the 3' terminus of TRV RNA2 (PSG strain). The nucleotide sequences of all tobraviral RNAs analysed so far indicate that they all possess a similar 3' terminal structure. A domain resembling the anticodon arm of canonical tRNA is not readily recognizable. TRV RNA2 can be adenylated with CTP, ATP; tRNA nucleotidyl transferase and ATP. It is unable, however, to accept any of the twenty common amino acids when incubated with ATP and aminoacyl-tRNA synthetases from wheat germ or yeast. We conclude that TRV RNA contains a tRNA-like structure, which, in contrast to the tymo- and tobamoviral tRNA-like structures, cannot be aminoacylated. It is unlikely therefore, that aminoacylation of plant viral RNAs with a tRNA-like structure is a prerequisite for viral RNA replication.     

Quantum Chemical Studies of Structures and Binding in Noncanonical RNA Base pairs: The Trans Watson-Crick:Watson-Crick Family

Abstract

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6–31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from ?8.19 kcal/mol to ?21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32–2.19 Å. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.     

Plant small nuclear RNAs. Nucleolar U3 snRNA is present in plants: partial characterization

Nuclei, isolated from a number of plant species by either of two independent, newly developed methods, regularly contained a common set of low-molecular-mass RNAs. Partial characterization of these RNAs, based on cell fractionation, polyacrylamide gel electrophoretic and chemical sequencing techniques, as well as comparison with literature data, revealed that, in addition to tRNA, 5S RNA and 5.8S RNA, plant nuclei contain two families of low-molecular-mass RNAs, that are counterparts of vertebrate U1 and U5 RNAs respectively, and three individual low-molecular-mass RNA species. One of these may be related to vertebrate U6 RNA. The two others are true eukaryotic U2 and U3 RNAs, respectively, on the basis of the following lines of evidence obtained from analyses of broad bean nuclear RNAs. The 3'-end portion (121 nucleotides sequenced) of broad bean U2 RNA shows a nearly perfect sequence homology with that of authentic pea U2 RNA. Broad bean U3 RNA is localized in the nucleolus and its 3'-end portion (164 nucleotides sequenced) (a) shows sequence homology with that of both rat U3 RNA (48%) and Dictyostelium D2 RNA (39%), (b) has a secondary structure which fits perfectly that proposed for both rat U3 RNA and Dictyostelium D2 RNA, and (c) contains the specific sequence which, in a model based on the primary structure of rat U3 RNA, is supposed to be involved in the processing of eukaryotic 32S pre-ribosomal RNA. This is the first report on the occurrence in plants of nucleolar U3 RNA.     

In vivo incorporation of the intrinsic photolabel 4-thiouridine into Escherichia coli RNAs.

The in vivo incorporation of the photoactivable uridine analogue 4-thiouridine into the RNAs of an Escherichia coli K12 pyrD strain has been optimised. s4Urd uptake in RNAs appears to be strikingly dependent upon the age of the preculture, i.e. the number of generations the cells have undergone immediately before dilution in the thiolation medium. Conditions have been set up where efficient RNA thiolation occurs in cells growing exponentially at 50 to 70% the rate of the control. The substitution level s4U/U is maximal after growth for 9 to 10 generations in the thiolation medium and reaches 17 +/- 3% in tRNA and bulk RNA. Most of ribosomal derived ribonucleoproteins, 65 +/- 5%, sediment as 70S ribosomes (s4U/U = 7 +/- 2%) on a high Mg2+ sucrose gradient. The thiolated RNAs were characterized by their migration on a thiol-specific affinity electrophoretic gel.     

Structural studies on tRNA acceptor stem microhelices: exchange of the discriminator base A73 for G in human tRNALeu switches the acceptor specificity from leucine to serine possibly by decreasing the stability of the terminal G1-C72 base pair.

Correct recognition of transfer RNAs (tRNAs) by aminoacyl-tRNA synthetases (aaRS) is crucial to the maintenance of translational fidelity. The discriminator base A73 in human tRNALeuis critical for its specific recognition by the aaRS. Exchanging A73 for G abolishes leucine acceptance and converts it into a serine acceptor in vitro . Two RNA microhelices of 24 nt length that correspond to the tRNALeuacceptor stem and differ only in the discriminator base were synthesized: a wild-type tRNALeumicrohelix, where nt 21 corresponds to the discriminator base position 73, and an A21G mutant microhelix. To investigate whether different identities of both tRNAs are caused by conformational differences, NMR and UV melting experiments were performed on both microhelices. Two-dimentional NOESY spectra showed both microhelices to exhibit the same overall conformation at their 3'-CCA ends. Thermodynamic analysis and melting behaviour of the base-paired imino protons observed by NMR spectroscopy suggest that the A21G (A73G in tRNA) exchange results in a decrease of melting transition cooperativity and a destabilization of the terminal G1-C20 (G1-C72 in tRNA) base pair. Furthermore, the fact that this 3'-terminal imino proton is more solvent-exposed at physiological temperature might be another indication for the importance of the stability of the terminal base pair for specific tRNA recognition.     

The Q-base of asparaginyl-tRNA is dispensable for efficient -1 ribosomal frameshifting in eukaryotes

The frameshift signal of the avian coronavirus infectious bronchitis virus (IBV) contains two cis-acting signals essential for efficient frameshifting, a heptameric slippery sequence (UUUAAAC) and an RNA pseudoknot structure located downstream. The frameshift takes place at the slippery sequence with the two ribosome-bound tRNAs slipping back simultaneously by one nucleotide from the zero phase (U UUA AAC) to the -1 phase (UUU AAA). Asparaginyl-tRNA, which decodes the A-site codon AAC, has the modified base Q at the wobble position of the anticodon (5' QUU 3') and it has been speculated that Q may be required for frameshifting. To test this, we measured frameshifting in cos cells that had been passaged in growth medium containing calf serum or horse serum. Growth in horse serum, which contains no free queuine, eliminates Q from the cellular tRNA population upon repeated passage. Over ten cell passages, however, we found no significant difference in frameshift efficiency between the cell types, arguing against a role for Q in frameshifting. We confirmed that the cells cultured in horse serum were devoid of Q by purifying tRNAs and assessing their Q-content by tRNA transglycosylase assays and coupled HPLC-mass spectroscopy. Supplementation of the growth medium of cells grown either on horse serum or calf serum with free queuine had no effect on frameshifting either. These findings were recapitulated in an in vitro system using rabbit reticulocyte lysates that had been largely depleted of endogenous tRNAs and resupplemented with Q-free or Q-containing tRNA populations. Thus Q-base is not required for frameshifting at the IBV signal and some other explanation is required to account for the slipperiness of eukaryotic asparaginyl-tRNA.