tRNA Secondary Structures Research Articles

The Mitochondrial Sequences of Heptathela hangzhouensis and Ornithoctonus huwena Reveal Unique Gene Arrangements and Atypical tRNAs

We have sequenced the complete mitochondrial genomes of the spiders Heptathela hangzhouensis and Ornithoctonus huwena. Both genomes encode 13 protein-coding genes, 22 tRNA genes, and 2 ribosomal RNA genes. H. hangzhouensis, a species of the suborder Mesothelae and a representative of the most basal clade of Araneae, possesses a gene order identical to that of Limulus polyphemus of Xiphosura. On the other hand, O. huwena, a representative of suborder Opisthothelae, infraorder Mygalomorphae, was found to have seven tRNA genes positioned differently from those of Limulus. The rrnL-trnL1-nad1 arrangement shared by the araneomorph families Salticidae, Nesticidae, and Linyphiidae and the mygalomorph family Theraphosidae is a putative synapomorphy joining the mygalomorph with the araneomorph. Between the two species examined, base compositions also differ significantly. The lengths of most protein-coding genes in H. hangzhouensis and O. huwena mtDNA are either identical to or slightly shorter than their Limulus counterparts. Usage of initiation and termination codons in these protein-coding genes seems to follow patterns conserved among most arthropod and some other metazoan mitochondrial genomes. The sequences of the 3' ends of rrnS and rrnL in the two species are similar to those reported for Limulus, and the entire genes are shortened by about 100-250 nucleotides with respect to Limulus. The lengths of most tRNA genes from the two species are distinctly shorter than those of Limulus and the sequences reveal unusual inferred tRNA secondary structures. Our finding provides new molecular evidence supporting that the suborder Mesothelae is basal to opisthothelids.

Corpus based learning of stochastic, context-free grammars combined with Hidden Markov Models for tRNA modelling

In this paper, a new method for modelling tRNA secondary structures is presented. This method is based on the combination of stochastic context-free grammars (SCFG) and Hidden Markov Models (HMM). HMM are used to capture the local relations in the loops of the molecule (nonstructured regions) and SCFG are used to capture the long term relations between nucleotides of the arms (structured regions). Given annotated public databases, the HMM and SCFG models are learned by means of automatic inductive learning methods. Two SCFG learning methods have been explored. Both of them take advantage of the structural information associated with the training sequences: one of them is based on a stochastic version of the Sakakibara algorithm and the other one is based on a Corpus based algorithm. A final model is then obtained by merging of the HMM of the nonstructured regions and the SCFG of the structured regions. Finally, the performed experiments on the tRNA sequence corpus and the non-tRNA sequence corpus give significant results. Comparative experiments with another published method are also presented.

Aminoacylation properties of pathology-related human mitochondrial tRNA(Lys) variants.

In vitro transcription has proven to be a successful tool for preparation of functional RNAs, especially in the tRNA field, in which, despite the absence of post-transcriptional modifications, transcripts are correctly folded and functionally active. Human mitochondrial (mt) tRNA(Lys) deviates from this principle and folds into various inactive conformations, due to the absence of the post-transcriptional modification m(1)A9 which hinders base-pairing with U64 in the native tRNA. Unavailability of a functional transcript is a serious drawback for structure/function investigations as well as in deciphering the molecular mechanisms by which point mutations in the mt tRNA(Lys) gene cause severe human disorders. Here, we show that an engineered in vitro transcribed "pseudo-WT" tRNA(Lys) variant is efficiently recognized by lysyl-tRNA synthetase and can substitute for the WT tRNA as a valuable reference molecule. This has been exploited in a systematic analysis of the effects on aminoacylation of nine pathology-related mutations described so far. The sole mutation located in a loop of the tRNA secondary structure, A8344G, does not affect aminoacylation efficiency. Out of eight mutations located in helical domains converting canonical Watson-Crick pairs into G-U pairs or C.A mismatches, six have no effect on aminoacylation (A8296G, U8316C, G8342A, U8356C, U8362G, G8363A), and two lead to drastic decreases (5000- to 7000-fold) in lysylation efficiencies (G8313A and G8328A). This screening, allowing for analysis of the primary impact level of all mutations affecting one tRNA under comparable conditions, indicates distinct molecular origins for different disorders.

ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences.

A computer program, ARAGORN, identifies tRNA and tmRNA genes. The program employs heuristic algorithms to predict tRNA secondary structure, based on homology with recognized tRNA consensus sequences and ability to form a base-paired cloverleaf. tmRNA genes are identified using a modified version of the BRUCE program. ARAGORN achieves a detection sensitivity of 99% from a set of 1290 eubacterial, eukaryotic and archaeal tRNA genes and detects all complete tmRNA sequences in the tmRNA database, improving on the performance of the BRUCE program. Recently discovered tmRNA genes in the chloroplasts of two species from the 'green' algae lineage are detected. The output of the program reports the proposed tRNA secondary structure and, for tmRNA genes, the secondary structure of the tRNA domain, the tmRNA gene sequence, the tag peptide and a list of organisms with matching tmRNA peptide tags.

BRUCE: a program for the detection of transfer-messenger RNA genes in nucleotide sequences

A computer program, BRUCE, was developed for the identification of transfer-messenger RNA (tmRNA) genes. The program employs heuristic algorithms to search for a tRNA(Ala)-like secondary structure surrounding a short sequence encoding the tag peptide. In the 57 completely sequenced bacterial genomes where tmRNA genes have been reported previously, BRUCE identified all with no false positives. In addition, BRUCE found 99 of the 100 tmRNAs identified previously in other bacteria, red chloroplasts and cyanelles. The output of the program reports the proposed tRNA secondary structure, the tmRNA gene sequence and the tag peptide.

Hybridization of antisense oligonucleotides with yeast tRNAPhe: factors determining the efficiency of interaction

With the object of elucidating the factors that affect the rate and efficiency of oligonucleotide hybridization with complementary regions in natural RNA, the interaction of tRNAPhe and in vitro transcript of this tRNA with synthetic oligonucleotides matching the 32—61 region comprising the variable loop and anticodon and TΨC stems of the molecule was studied in detail. Gel retardation was used to assay thermodynamic and kinetic parameters of tRNA interaction with oligonucleotides. Oligonucleotide hybridization rates do not depend on the oligonucleotide length and are determined by the stability of RNA structure within the binding site. The rate constant for hybridization correlates with the number of base pairs in the RNA opened during binding with the oligonucleotide. Factors stabilizing the tRNA secondary structure (the presence of Mg2+ ions, low temperature, hypermodified bases) decrease significantly both the rate and efficiency of hybridization. Analysis of the tRNAPhe structure in the presence of an oligonucleotide shows that binding of oligonucleotides with the 45—61 region of tRNAPhe results only in unfolding of the TΨC hairpin and does not disturb the secondary structure of the rest of the molecule.

Complete 5' and 3' end maturation of group II intron-containing tRNA precursors.

Higher plant chloroplasts provide the only experimentally validated example of functional tRNA genes that are disrupted by group II introns. Here, precursor transcripts for tRNA(Gly)(UCC), tRNA(Val)(UAC), and tRNA(Ala)(UGC) were investigated for processing of 5' leader and 3' trailer sequences in vivo. Use of intron-specific primer pairs and inclusion of a barley chloroplast splicing mutant specifically allowed us to evaluate the potential effect of intervening sequences that disrupt tRNA secondary and tertiary structures. The data suggest that (1) neither integrity of the dihydrouridine nor the anticodon domain is required for the nucleotidyltransferase-mediated addition of 3'-terminal CCA; (2) interruption of these two structural elements by group II introns does not interfere with nucleotide-specific 5' maturation by RNase P; (3) processing intermediates of chloroplast tRNAs can be 3' polyadenylated; and (4) plastid DNA-encoded proteins are not required for 3' and 5' maturation of plastid tRNAs.

Mitochondrial genome organization and vertebrate phylogenetics

With the advent of DNA sequencing techniques the organization of the vertebrate mitochondrial genome shows variation between higher taxonomic levels. The most conserved gene order is found in placental mammals, turtles, fishes, some lizards and Xenopus. Birds, other species of lizards, crocodilians, marsupial mammals, snakes, tuatara, lamprey, and some other amphibians and one species of fish have gene orders that are less conserved. The most probable mechanism for new gene rearrangements seems to be tandem duplication and multiple deletion events, always associated with tRNA sequences. Some new rearrangements seem to be typical of monophyletic groups and the use of data from these groups may be useful for answering phylogenetic questions involving vertebrate higher taxonomic levels. Other features such as the secondary structure of tRNA, and the start and stop codons of protein-coding genes may also be useful in comparisons of vertebrate mitochondrial genomes.

Mitochondrial sequence evolution in spiders: intraspecific variation in tRNAs lacking the TPsiC Arm.

Analyses of mitochondrial DNA sequences from three species of Habronattus jumping spiders (Chelicerata: Arachnida: Araneae) reveal unusual inferred tRNA secondary structures and gene arrangements, providing new information on tRNA evolution within chelicerate arthropods. Sequences from the protein-coding genes NADH dehydrogenase subunit 1 (ND1), cytochrome oxidase subunit I (COI), and subunit II (COII) were obtained, along with tRNA, tRNA, and large-subunit ribosomal RNA (16S) sequences; these revealed several peculiar features. First, inferred secondary structures of tRNA and, likely, tRNA, lack the TPsiC arm and the variable arm and therefore do not form standard cloverleaf structures. In place of these arms is a 5-6-nt T arm-variable loop (TV) replacement loop such as that originally described from nematode mitochondrial tRNAs. Intraspecific variation occurs in the acceptor stem sequences in both tRNAs. Second, while the proposed secondary structure of the 3' end of 16S is similar to that reported for insects, the sequence at the 5' end is extremely divergent, and the entire gene is truncated about 300 nt with respect to Drosophila yakuba. Third, initiation codons appear to consist of ATY (ATT and ATC) and TTG for ND1 and COII, respectively. Finally, Habronattus shares the same ND1-tRNA-16S gene arrangement as insects and crustaceans, thus illustrating variation in a tRNA gene arrangement previously proposed as a character distinguishing chelicerates from insects and crustaceans.

Mytilus mitochondrial DNA contains a functional gene for a tRNASer(UCN) with a dihydrouridine arm-replacement loop and a pseudo-tRNASer(UCN) gene.

A 2500-nucleotide pair (ntp) sequence of F-type mitochondrial (mt) DNA of the Pacific Rim mussel Mytilus californianus (class Bivalvia, phylum Mollusca) that contains two complete (ND2 and ND3) and two partial (COI and COIII) protein genes and nine tRNA genes is presented. Seven of the encoded tRNAs (Ala, Arg, His, Met(AUA), Pro, Ser(UCN), and Trp) have the potential to fold into the orthodox four-armed tRNA secondary structure, while two [tRNASer(AGN) and a second tRNASer(UCN)] will fold only into tRNAs with a dihydrouridine (DHU) arm-replacement loop. Comparison of these mt-tRNA gene sequences with previously published, corresponding M. edulis F-type mtDNA indicates that similarity between the four-armed tRNASer(UCN) genes is only 63.8% compared with an average of 92.1% (range 86.2-98. 5%) for the remaining eight tRNA genes. Northern blot analysis indicated that mature tRNAs encoded by the DHU arm-replacement loop-containing tRNASer(UCN), tRNASer(AGN), tRNAMet(AUA), tRNATrp, and tRNAPro genes occur in M. californianus mitochondria, strengthening the view that all of these genes are functional. However, Northern blot and 5' RACE (rapid amplification of cDNA ends) analyses indicated that the four-armed tRNASer(UCN) gene is transcribed into a stable RNA that includes the downstream COI sequence and is not processed into a mature tRNA. On the basis of these observations the M. californianus and M. edulis four-armed tRNASer(UCN) sequences are interpreted as pseudo-tRNASer(UCN) genes.

Complete suboptimal folding of RNA and the stability of secondary structures.

An algorithm is presented for generating rigorously all suboptimal secondary structures between the minimum free energy and an arbitrary upper limit. The algorithm is particularly fast in the vicinity of the minimum free energy. This enables the efficient approximation of statistical quantities, such as the partition function or measures for structural diversity. The density of states at low energies and its associated structures are crucial in assessing from a thermodynamic point of view how well-defined the ground state is. We demonstrate this by exploring the role of base modification in tRNA secondary structures, both at the level of individual sequences from Escherichia coli and by comparing artificially generated ensembles of modified and unmodified sequences with the same tRNA structure. The two major conclusions are that (1) base modification considerably sharpens the definition of the ground state structure by constraining energetically adjacent structures to be similar to the ground state, and (2) sequences whose ground state structure is thermodynamically well defined show a significant tendency to buffer single point mutations. This can have evolutionary implications, since selection pressure to improve the definition of ground states with biological function may result in increased neutrality.

Complete suboptimal folding of RNA and the stability of secondary structures

An algorithm is presented for generating rigorously all suboptimal secondary structures between the minimum free energy and an arbitrary upper limit. The algorithm performs particularly well in the vicinity of the minimum free energy. This enables the efficient approximation of statistical quantities, such as the partition function or measures for structural diversity. The density of states at low energies and its associated structures are crucial in assessing from a thermodynamic point of view how well defined the ground state is. We demonstrate this by exploring the role of base modification in tRNA secondary structures, both at the level of individual sequences from E. coli and by comparing artificially generated ensembles of modified and unmodified sequences with the same tRNA structure. The two major conclusions are that (1) base modification considerably sharpens the definition of the ground state structure by constraining energetically adjacent structures to be similar to the ground state, and (2) sequences whose ground state structure is thermodynamically well defined show a significant tendency to buffer single point mutations. This can have evolutionary implications, since selection pressure to improve the definition of ground states with biological function may result in increased neutrality.

In vitro selection of RNAs aminoacylated by Escherichia coli leucyl-tRNA synthetase

To investigate systematically the RNA sequences necessary for aminoacylation by Escherichia coli leucyl-tRNA synthetase, RNAs with leucylation activity were isolated by in vitro selection from a library of tRNA Leu variants possessing randomized sequences in the D-loop, the variable arm, and the T-loop. After two rounds of selection, most of the selected variants showed the following features: (1) the tertiary interaction between nucleotides at positions 15 and 48 was A15-U48; (2) the continuous G18G19 sequence, which is invariant in canonical tRNAs, appeared at the fixed position in the D-loop; and (3) the nucleotide at position 20a in the D-loop was A. These selected nucleotides and their positions, concentrating on the hinge region of tRNA, were identical to those of native tRNA Leu. In contrast, although the long variable arm is the most characteristic of the tRNA Leu structure, the primary and secondary structures were not correlated with the leucylation activity. These findings indicate that A15-U48, A20a, and G18G19 located at specific positions are involved in the tertiary folding of leucine-accepting tRNA molecules. With increases in the selection cycle, the D-loop sequence and the secondary structure of the variable arm became similar to those of tRNA Leu, suggesting that tRNA Leu represents an optimized RNA sequence for leucylation.

A disease-associated G5703A mutation in human mitochondrial DNA causes a conformational change and a marked decrease in steady-state levels of mitochondrial tRNA(Asn).

We introduced mitochondrial DNA (mtDNA) from a patient with a mitochondrial myopathy into established mtDNA-less human osteosarcoma cells. The resulting transmitochondrial cybrid lines, containing either exclusively wild-type or mutated (G5703A transition in the tRNA[Asn] gene) mtDNA, were characterized and analyzed for oxidative phosphorylation function and steady-state levels of different RNA species. Functional studies showed that the G5703A mutation severely impairs oxidative phosphorylation function and mitochondrial protein synthesis. We detected a marked reduction in tRNA(Asn) steady-state levels which was not associated with an accumulation of intermediate transcripts containing tRNA(Asn) sequences or decreased transcription. Native polyacrylamide gel electrophoresis showed that the residual tRNA(Asn) fraction in mutant cybrids had an altered conformation, suggesting that the mutation destabilized the tRNA(Asn) secondary or tertiary structure. Our results suggest that the G5703 mutation causes a conformational change in the tRNA(Asn) which may impair aminoacylation. This alteration leads to a severe reduction in the functional tRNA(Asn) pool by increasing its in vivo degradation by mitochondrial RNases.

Expression of bovine mitochondrial tRNASer GCU derivatives in Escherichia coli.

By replacing a stretch of five A-U base pairs in the acceptor stem with G-C pairs, mitochondrial tRNA-SerGCU lacking a D arm could be expressed in Escherichia coli cells in considerable amounts. The expressed tRNA with no modified nucleoside was serylated in vitro with the mitochondrial enzyme. The tRNASerGCU derivatives carrying identity elements for alanine tRNA and the related anticodons were expressed. However, this expression event did not affect cell growth, probably because the expression started from the late log phase, which suggests that these mitochondrial tRNA derivatives are not involved in E.coli gene expression systems. Although there are some restrictions in the secondary structure of tRNAs that can be expressed by this method, it could prove useful for preparing large amounts of heterologous tRNAs in vivo.

Evolutionary analysis of sea urchin mitochondrial tRNAs: folding of the molecules as suggested by the non-random occurrence of nucleotides

Comparative analyses of the mitochondrial tRNA sequences of the sea urchins Arbacia lixula, Paracentrotus lividus and Strongylocentrotus purpuratus revealed that conserved nucleotides may be involved in determining the typical L-shaped spatial conformation of tRNAs. These results shed light on the specific tertiary interactions that allow the folding of the atypical mitochondrial tRNAs into a functional form. A consensus mitochondrial tRNA secondary structure was derived. It shows the presence of nucleotides virtually conserved only in these organisms that represent a sort of molecular signature in sea urchins and suggests a possible physiological role. Finally, we speculate that the non-canonical structure of animal tRNAs, as well as the deviations from the universality of the genetic code, may be due to the reduction in size of the metazoan mitochondrial genome, with the concomitant acquisition of new functions by the mitochondrial tRNAs.

RNA editing of tRNA(Phe) and tRNA(Cys) in mitochondria of Oenothera berteriana is initiated in precursor molecules.

We have analyzed the role of RNA editing in the correction of mismatched base pairs in tRNA secondary structures in mitochondria of the flowering plant Oenothera berteriana. Comparison of genomic and cDNA sequences from unprocessed primary transcripts of the newly characterized genes for tRNA(Cys), tRNA(Asn) and tRNA(Ile) and the previously described gene for tRNA(Phe) revealed single nucleotide discrepancies in the tRNA(Cys) and tRNA(Phe) sequences. While the change in the anticodon stem of tRNA(Cys) alters a C-T to a T-T mismatch, the nucleotide transition in the tRNA(Phe) restores a conventional T-A Watson-Crick base pair, replacing a C-A mismatch in the acceptor stem. Since both nucleotide alterations are conversions from genomic cytidines to thymidines in the cDNA (uridines in the tRNAs), they are attributed to RNA editing, which is observed in nearly all mRNAs from plant mitochondria.

RNA sequence analysis using covariance models.

We describe a general approach to several RNA sequence analysis problems using probabilistic models that flexibly describe the secondary structure and primary sequence consensus of an RNA sequence family. We call these models 'covariance models'. A covariance model of tRNA sequences is an extremely sensitive and discriminative tool for searching for additional tRNAs and tRNA-related sequences in sequence databases. A model can be built automatically from an existing sequence alignment. We also describe an algorithm for learning a model and hence a consensus secondary structure from initially unaligned example sequences and no prior structural information. Models trained on unaligned tRNA examples correctly predict tRNA secondary structure and produce high-quality multiple alignments. The approach may be applied to any family of small RNA sequences.

Computer simulation of tRNA secondary structure folding.

Computer simulation results of folding linear RNA molecules into secondary structures are presented. The structure is formed by two interacting processes: the RNA molecular chain growth (beginning from an initial length, L0), and the structuring (secondary structure sequential growth in the region of the existing molecular chain, based on the local free energy minimization by sequential addition of elementary substructures--stems). It was found that the final secondary structure formation is greatly influenced by the 'structuring period' T (the ratio of the molecular chain growth rate to the structuring rate), and the direction of RNA synthesis. The computer simulation has been performed for 219 and 906 tRNA genes from two published catalogues, on the whole two-dimensional domain (T,L0) parameters, by using four known free-energy models. Minimum stem length and molecular chain growth direction have been also varied. The calculated secondary structures have been compared to the natural tRNA structures given in the catalogues, and the region of best coincidence for the model parameters has been determined. It has been proved that, on average, > 86% of the paired bases of natural tRNA structures appear in the folding simulation.

Correlations between fluorine-19 nuclear magnetic resonance chemical shift and the secondary and tertiary structure of 5-fluorouracil-substituted tRNA.

To complete assignment of the 19F nuclear magnetic resonance (NMR) spectrum of 5-fluorouracil-substituted Escherichia coli tRNA(Val), resonances from 5-fluorouracil residues involved in tertiary interactions have been identified. Because these assignments could not be made directly by the base-replacement method used to assign 5-fluorouracil residues in loop and stem regions of the tRNA, alternative assignment strategies were employed. FU54 and FU55 were identified by 19F homonuclear Overhauser experiments and were then assigned by comparison of their 19F NMR spectra with those of 5-fluorouracil-labeled yeast tRNA(Phe) mutants having FU54 replaced by adenine and FU55 replaced by cytosine. FU8 and FU12, were assigned from the 19F NMR spectrum of the tRNA(Val) mutant in which the base triple G9-C23-G12 substituted for the wild-type A9-A23-FU12. Although replacement of the conserved U8 (FU8) with A or C disrupts the tertiary structure of tRNA(Val), it has only a small effect on the catalytic turnover number of valyl-tRNA synthetase, while reducing the affinity of the tRNA for enzyme. Analysis of the 19F chemical shift assignments of all 14 resonances in the spectrum of 5-fluorouracil-substituted tRNAVal indicated a strong correlation to tRNA secondary and tertiary structure. 5-Fluorouracil residues in loop regions gave rise to peaks in the central region of the spectrum, 4.4 to 4.9 parts per million (p.p.m.) downfield from free 5-fluorouracil. However, the signal from FU59, in the T-loop of tRNA(Val), was shifted more than 1 p.p.m. downfield, to 5.9 p.p.m., presumably because of the involvement of this fluorouracil in the tertiary interactions between the T and D-loops. The 19F chemical shift moved upfield, to the 2.0 to 2.8 p.p.m. range, when fluorouracil was base-paired with adenine in helical stems. This upfield shift was less pronounced for the fluorine of the FU7.A66 base-pair, located at the base of the acceptor stem, an indication that FU7 is only partially stacked on the adjacent G49 in the continuous acceptor stem/T-stem helix. An unanticipated finding was that the 19F resonances of 5-fluorouracil residues wobble base-paired with guanine were shifted 4 to 5 p.p.m. downfield of those from fluorouracil residues paired with A. In the 19F NMR spectra of all fluorinated tRNAs studied, the farthest downfield peak corresponded to FU55, which replaced the conserved pseudouridine normally found at this position.