A model for the tertiary structure of mammalian mitochondrial transfer RNAs lacking the entire 'dihydrouridine' loop and stem

The complete mitochondrial genome of the imperiled Bullnose ray Myliobatis freminvillei (Myliobatiformes: Myliobatidae) with comments on its phylogenetic position and claims of diversifying selection affecting protein coding genes in a closely related species

The Bullnose ray Myliobatis freminvillei is a bentho-pelagic eagle ray that inhabits the north Gulf of Mexico and the Western Atlantic Ocean Coast, discontinuously, from Massachusetts, USA to Buenos Aires, Argentina. Myliobatis freminvillei is currently listed as vulnerable by the 2019 IUCN Red List of Threatened Species given that it is often captured as bycatch by artisanal and commercial fisheries, along the coasts of Argentina, Brazil, and Venezuela. This study, for the first time, assembled and characterized the mitochondrial genome of M. fremnvillei. The A+T rich mitochondrial genome of M. fremnvillei is 18,356 bp long and encodes 22 transfer RNA genes (tRNA), 2 ribosomal RNA genes (12S ribosomal RNA and 16S ribosomal RNA), 13 protein coding genes (PCGs), and also contains a non coding control region of 2,617 bp. Nonsynonomous codon usage with a preference for A+T rich codons was observed in all 13 PCGs. Leu (CTA), Ile (ATC), Phe (TTC), Thr (ACA), and Ala (GCC) were the most frequently used codons. Ka/Ks ratios estimated for all 13 PCGs exhibited values < 1, indicating strong purifying selection affecting all these genes. In contrast to the results of a previous study that claimed diversifying selective pressure in two mitochondrial PCGs of Mobula tarapacana, reanalysis of the Ka/Ks values for the same species indicated purifying selection in all 13 PCGs. Of the 22 tRNA genes, all have a cloverleaf secondary structure except tRNA-Ser1 which has a truncated dihydrouridine arm. In the control region, A+T rich microsatellites (n = 42) and short tandem repeats (n = 6) were identified, and the secondary structure of the same region contained numerous hairpin loops. Phylomitogenomic analyses supported the monophyletic status of the order Myliobatiformes and family Myliobatidae. The assembled mitochondrial genome will assist with conservation efforts in Myliobatis fremnvillei.

Structural basis for a degenerate tRNA identity code and the evolution of bimodal specificity in human mitochondrial tRNA recognition

Animal mitochondrial gene expression relies on specific interactions between nuclear-encoded aminoacyl-tRNA synthetases and mitochondria-encoded tRNAs. Their evolution involves an antagonistic interplay between strong mutation pressure on mtRNAs and selection pressure to maintain their essential function. To understand the molecular consequences of this interplay, we analyze the human mitochondrial serylation system, in which one synthetase charges two highly divergent mtRNASer isoacceptors. We present the cryo-EM structure of human mSerRS in complex with mtRNASer(UGA), and perform a structural and functional comparison with the mSerRS-mtRNASer(GCU) complex. We find that despite their common function, mtRNASer(UGA) and mtRNASer(GCU) show no constrain to converge on shared structural or sequence identity motifs for recognition by mSerRS. Instead, mSerRS evolved a bimodal readout mechanism, whereby a single protein surface recognizes degenerate identity features specific to each mtRNASer. Our results show how the mutational erosion of mtRNAs drove a remarkable innovation of intermolecular specificity rules, with multiple evolutionary pathways leading to functionally equivalent outcomes.

The diverse structural modes of tRNA binding and recognition

tRNAs are short noncoding RNAs responsible for decoding mRNA codon triplets, delivering correct amino acids to the ribosome, and mediating polypeptide chain formation. Due to their key roles during translation, tRNAs have a highly conserved shape and large sets of tRNAs are present in all living organisms. Regardless of sequence variability, all tRNAs fold into a relatively rigid three-dimensional L-shaped structure. The conserved tertiary organization of canonical tRNA arises through the formation of two orthogonal helices, consisting of the acceptor and anticodon domains. Both elements fold independently to stabilize the overall structure of tRNAs through intramolecular interactions between the D- and T-arm. During tRNA maturation, different modifying enzymes posttranscriptionally attach chemical groups to specific nucleotides, which not only affect translation elongation rates but also restrict local folding processes and confer local flexibility when required. The characteristic structural features of tRNAs are also employed by various maturation factors and modification enzymes to assure the selection, recognition, and positioning of specific sites within the substrate tRNAs. The cellular functional repertoire of tRNAs continues to extend well beyond their role in translation, partly, due to the expanding pool of tRNA-derived fragments. Here, we aim to summarize the most recent developments in the field to understand how three-dimensional structure affects the canonical and noncanonical functions of tRNA.

Structural basis for shape-selective recognition and aminoacylation of a D-armless human mitochondrial tRNA

Human mitochondrial gene expression relies on the specific recognition and aminoacylation of mitochondrial tRNAs (mtRNAs) by nuclear-encoded mitochondrial aminoacyl-tRNA synthetases (mt-aaRSs). Despite their essential role in cellular energy homeostasis, strong mutation pressure and genetic drift have led to an unparalleled sequence erosion of animal mtRNAs. The structural and functional consequences of this erosion are not understood. Here, we present cryo-EM structures of the human mitochondrial seryl-tRNA synthetase (mSerRS) in complex with mtRNASer(GCU). These structures reveal a unique mechanism of substrate recognition and aminoacylation. The mtRNASer(GCU) is highly degenerated, having lost the entire D-arm, tertiary core, and stable L-shaped fold that define canonical tRNAs. Instead, mtRNASer(GCU) evolved unique structural innovations, including a radically altered T-arm topology that serves as critical identity determinant in an unusual shape-selective readout mechanism by mSerRS. Our results provide a molecular framework to understand the principles of mito-nuclear co-evolution and specialized mechanisms of tRNA recognition in mammalian mitochondrial gene expression.

The Dihydrouridine landscape from tRNA to mRNA: a perspective on synthesis, structural impact and function

The universal dihydrouridine (D) epitranscriptomic mark results from a reduction of uridine by the Dus family of NADPH-dependent reductases and is typically found within the eponym D-loop of tRNAs. Despite its apparent simplicity, D is structurally unique, with the potential to deeply affect the RNA backbone and many, if not all, RNA-connected processes. The first landscape of its occupancy within the tRNAome was reported 20 years ago. Its potential biological significance was highlighted by observations ranging from a strong bias in its ecological distribution to the predictive nature of Dus enzymes overexpression for worse cancer patient outcomes. The exquisite specificity of the Dus enzymes revealed by a structure-function analyses and accumulating clues that the D distribution may expand beyond tRNAs recently led to the development of new high-resolution mapping methods, including Rho-seq that established the presence of D within mRNAs and led to the demonstration of its critical physiological relevance.

Complete chemical structures of human mitochondrial tRNAs

Mitochondria generate most cellular energy via oxidative phosphorylation. Twenty-two species of mitochondrial (mt-)tRNAs encoded in mtDNA translate essential subunits of the respiratory chain complexes. mt-tRNAs contain post-transcriptional modifications introduced by nuclear-encoded tRNA-modifying enzymes. They are required for deciphering genetic code accurately, as well as stabilizing tRNA. Loss of tRNA modifications frequently results in severe pathological consequences. Here, we perform a comprehensive analysis of post-transcriptional modifications of all human mt-tRNAs, including 14 previously-uncharacterized species. In total, we find 18 kinds of RNA modifications at 137 positions (8.7% in 1575 nucleobases) in 22 species of human mt-tRNAs. An up-to-date list of 34 genes responsible for mt-tRNA modifications are provided. We identify two genes required for queuosine (Q) formation in mt-tRNAs. Our results provide insight into the molecular mechanisms underlying the decoding system and could help to elucidate the molecular pathogenesis of human mitochondrial diseases caused by aberrant tRNA modifications.

The MRPP1/MRPP2 complex is a tRNA-maturation platform in human mitochondria

Mitochondrial polycistronic transcripts are extensively processed to give rise to functional mRNAs, rRNAs and tRNAs; starting with the release of tRNA elements through 5′-processing by RNase P (MRPP1/2/3-complex) and 3′-processing by RNase Z (ELAC2). Here, we show using in vitro experiments that MRPP1/2 is not only a component of the mitochondrial RNase P but that it retains the tRNA product from the 5′-processing step and significantly enhances the efficiency of ELAC2-catalyzed 3′-processing for 17 of the 22 tRNAs encoded in the human mitochondrial genome. Furthermore, MRPP1/2 retains the tRNA product after ELAC2 processing and presents the nascent tRNA to the mitochondrial CCA-adding enzyme. Thus, in addition to being an essential component of the RNase P reaction, MRPP1/2 serves as a processing platform for several down-stream tRNA maturation steps in human mitochondria. These findings are of fundamental importance for our molecular understanding of disease-related mutations in MRPP1/2, ELAC2 and mitochondrial tRNA genes.

Losing the stem-loop structure from metazoan mitochondrial tRNAs and co-evolution of interacting factors

Conventional tRNAs have highly conserved sequences, four-armed cloverleaf secondary structures, and L-shaped tertiary structures. However, metazoan mitochondrial tRNAs contain several exceptional structures. Almost all tRNAs^Ser for AGY/N codons lack the D-arm. Furthermore, in some nematodes, no four-armed cloverleaf-type tRNAs are present: two tRNAs^Ser without the D-arm and 20 tRNAs without the T-arm are found. Previously, we showed that in nematode mitochondria, an extra elongation factor Tu (EF-Tu) has evolved to support interaction with tRNAs lacking the T-arm, which interact with C-terminal domain 3 in conventional EF-Tu. Recent mitochondrial genome analyses have suggested that in metazoan lineages other than nematodes, tRNAs without the T-arm are present. Furthermore, even more simplified tRNAs are predicted in some lineages. In this review, we discuss mitochondrial tRNAs with divergent structures, as well as protein factors, including EF-Tu, that support the function of truncated metazoan mitochondrial tRNAs.

Idiosyncrasies in decoding mitochondrial genomes

Mitochondria originate from the α-proteobacterial domain of life. Since this unique event occurred, mitochondrial genomes of protozoans, fungi, plants and metazoans have highly derived and diverged away from the common ancestral DNA. These resulting genomes highly differ from one another, but all present-day mitochondrial DNAs have a very reduced coding capacity. Strikingly however, ATP production coupled to electron transport and translation of mitochondrial proteins are the two common functions retained in all mitochondrial DNAs. Paradoxically, most components essential for these two functions are now expressed from nuclear genes. Understanding how mitochondrial translation evolved in various eukaryotic models is essential to acquire new knowledge of mitochondrial genome expression. In this review, we provide a thorough analysis of the idiosyncrasies of mitochondrial translation as they occur between organisms. We address this by looking at mitochondrial codon usage and tRNA content. Then, we look at the aminoacyl-tRNA-forming enzymes in terms of peculiarities, dual origin, and alternate function(s). Finally we give examples of the atypical structural properties of mitochondrial tRNAs found in some organisms and the resulting adaptive tRNA-protein partnership.

Apoptotic regulation and tRNA

Apoptotic regulation is critical to organismal homeostasis and protection against many human disease processes such as cancer. Significant research efforts over the past several decades have illuminated many signaling molecules and effecter proteins responsible for this form of programmed cell death. Recent evidence suggests that transfer RNA (tRNA) regulates apoptotic sensitivity at the level of cytochrome c-mediated apoptosome formation. This finding unexpectedly places tRNA at the nexus of cellular biosynthesis and survival. Here we review the current understanding of both the apoptotic machinery and tRNA biology. We describe the evidence linking tRNA and cytochrome c in depth, and speculate on the implications of this link in cell biology.

Unique features of animal mitochondrial translation systems. The non-universal genetic code, unusual features of the translational apparatus and their relevance to human mitochondrial diseases

In animal mitochondria, several codons are non-universal and their meanings differ depending on the species. In addition, the tRNA structures that decipher codons are sometimes unusually truncated. These features seem to be related to the shortening of mitochondrial (mt) genomes, which occurred during the evolution of mitochondria. These organelles probably originated from the endosymbiosis of an aerobic eubacterium into an ancestral eukaryote. It is plausible that these events brought about the various characteristic features of animal mt translation systems, such as genetic code variations, unusually truncated tRNA and rRNA structures, unilateral tRNA recognition mechanisms by aminoacyl-tRNA synthetases, elongation factors and ribosomes, and compensation for RNA deficits by enlarged proteins. In this article, we discuss molecular mechanisms for these phenomena. Finally, we describe human mt diseases that are caused by modification defects in mt tRNAs.

Diversity and similarity in the tRNA world: overall view and case study on malaria-related tRNAs

Transfer RNAs (tRNAs) are ancient macromolecules that have evolved under various environmental pressures as adaptors in translation in all forms of life but also towards alternative structures and functions. The present knowledge on both "canonical" and "deviating" signature motifs retrieved from vertical and horizontal sequence comparisons is briefly reviewed. Novel characteristics, proper to tRNAs from a given translation system, are revealed by a case study on the nuclear and organellar tRNA sets from malaria-related organisms. Unprecedented distinctive features for Plasmodium falciparum apicoplastic tRNAs appear, which provide novel routes to be explored towards anti-malarial drugs. The ongoing high-throughput sequencing programs are expected to allow for further horizontal comparisons and to reveal other signatures of either full or restricted sets of tRNAs.

Core flexibility of a truncated metazoan mitochondrial tRNA

Secondary and tertiary structures of tRNAs are remarkably preserved from bacteria to humans, the notable exception being the mitochondrial (m) tRNAs of metazoans, which often deviate substantially from the canonical cloverleaf (secondary) or ‘L’-shaped (tertiary) structure. Many metazoan mtRNAs lack either the TψC (T) or dihydrouridine (D) loops of the canonical cloverleaf, which are known to confer structural rigidity to the folded structure. Thus, the absence of canonical TψC–D interactions likely results in greater dispersion of anticodon-acceptor interstem angle than for canonical tRNAs. To test this hypothesis, we have assessed the dispersion of the anticodon-acceptor angle for bovine mtRNA^Ser(AGY), which lacks the canonical D arm and is thus incapable of forming stabilizing interarm interactions. Using the method of transient electric birefringence (TEB), and by changing the helical torsion angle between a core mtRNA bend and a second bend of known angle/rigidity, we have demonstrated that the core of mtRNA^Ser(AGY) has substantially greater flexibility than its well-characterized canonical counterpart, yeast cytoplasmic tRNA^Phe. These results suggest that increased flexibility, in addition to a more open interstem angle, would allow both noncanonical and canonical mtRNAs to utilize the same protein synthetic apparatus.

Structure-specific tRNA determinants for editing a mischarged amino acid

Alanyl-tRNA synthetase efficiently aminoacylates tRNAAla and an RNA minihelix that comprises just one domain of the two-domain L-shaped tRNA structure. It also clears mischarged tRNAAla using a specialized domain in its C-terminal half. In contrast to full-length tRNAAla, minihelixAla was robustly mischarged and could not be edited. Addition in trans of the missing anticodon-containing domain did not activate editing of mischarged minihelixAla. To understand these differences between minihelixAla and tRNAAla, several chimeric full tRNAs were constructed. These had the acceptor stem of a non-cognate tRNA replaced with the stem of tRNAAla. The chimeric tRNAs collectively introduced multiple sequence changes in all parts but the acceptor stem. However, although the acceptor stem in isolation (as the minihelix) lacked determinants for editing, alanyl-tRNA synthetase effectively cleared a mischarged amino acid from each chimeric tRNA. Thus, a covalently continuous two-domain structure per se, not sequence, is a major determinant for clearance of errors of aminoacylation by alanyl-tRNA synthetase. Because errors of aminoacylation are known to be deleterious to cell growth, structure-specific determinants constitute a powerful selective pressure to retain the format of the two-domain L-shaped tRNA.

Dual mode recognition of two isoacceptor tRNAs by mammalian mitochondrial seryl-tRNA synthetase

Animal mitochondrial translation systems contain two serine tRNAs, corresponding to the codons AGY (Y = U and C) and UCN (N = U, C, A, and G), each possessing an unusual secondary structure; tRNA(GCU)(Ser) (for AGY) lacks the entire D arm, whereas tRNA(UGA)(Ser) (for UCN) has an unusual cloverleaf configuration. We previously demonstrated that a single bovine mitochondrial seryl-tRNA synthetase (mt SerRS) recognizes these topologically distinct isoacceptors having no common sequence or structure. Recombinant mt SerRS clearly footprinted at the TPsiC loop of each isoacceptor, and kinetic studies revealed that mt SerRS specifically recognized the TPsiC loop sequence in each isoacceptor. However, in the case of tRNA(UGA)(Ser), TPsiC loop-D loop interaction was further required for recognition, suggesting that mt SerRS recognizes the two substrates by distinct mechanisms. mt SerRS could slightly but significantly misacylate mitochondrial tRNA(Gln), which has the same TPsiC loop sequence as tRNA(UGA)(Ser), implying that the fidelity of mitochondrial translation is maintained by kinetic discrimination of tRNAs in the network of aminoacyl-tRNA synthetases.

Translation ability of mitochondrial tRNAsSer with unusual secondary structures in an in vitro translation system of bovine mitochondria

BACKGROUND

Metazoan mitochondrial (mt) tRNAs are structurally quite different from the canonical cloverleaf secondary structure. The mammalian mt tRNASerGCU for AGY codons (Y = C or U) lacks the entire D arm, whereas tRNASerUGA for UCN codons (N = A, G, C or U) has an extended anti-codon stem. It has been a long-standing problem to prove experimentally how these tRNAsSer work in the mt translation system.

RESULTS

To solve the above-mentioned problem, we examined their translational abilities in an in vitro bovine mitochondrial translation system using transcripts of altered tRNASer analogues derived from bovine mitochondria. Both tRNASer analogues had almost the same ability to form ternary complexes with mt EF-Tu and GTP. The D-arm-lacking tRNASer GCU analogue had considerably lower translational activity than the tRNASerUGA analogue and produced mostly short oligopeptides, up to a tetramer. In addition, tRNASerGCU analogue was disfavoured by the ribosome when other tRNAs capable of decoding the cognate codon were available.

CONCLUSION

Both mt tRNASerGCU and tRNASerUGA analogues with unusual secondary structure were found to be capable of translation on the ribosome. However, the tRNASerGCU analogue has some molecular disadvantage on the ribosome, which probably derives from the lack of a D arm.

Search for characteristic structural features of mammalian mitochondrial tRNAs

A number of mitochondrial (mt) tRNAs have strong structural deviations from the classical tRNA cloverleaf secondary structure and from the conventional L-shaped tertiary structure. As a consequence, there is a general trend to consider all mitochondrial tRNAs as "bizarre" tRNAs. Here, a large sequence comparison of the 22 tRNA genes within 31 fully sequenced mammalian mt genomes has been performed to define the structural characteristics of this specific group of tRNAs. Vertical alignments define the degree of conservation/variability of primary sequences and secondary structures and search for potential tertiary interactions within each of the 22 families. Further horizontal alignments ascertain that, with the exception of serine-specific tRNAs, mammalian mt tRNAs do fold into cloverleaf structures with mostly classical features. However, deviations exist and concern large variations in size of the D- and T-loops. The predominant absence of the conserved nucleotides G18G19 and T54T55C56, respectively in these loops, suggests that classical tertiary interactions between both domains do not take place. Classification of the tRNA sequences according to their genomic origin (G-rich or G-poor DNA strand) highlight specific features such as richness/poorness in mismatches or G-T pairs in stems and extremely low G-content or C-content in the D- and T-loops. The resulting 22 "typical" mammalian mitochondrial sequences built up a phylogenetic basis for experimental structural and functional investigations. Moreover, they are expected to help in the evaluation of the possible impacts of those point mutations detected in human mitochondrial tRNA genes and correlated with pathologies.

Search for differences in post-transcriptional modification patterns of mitochondrial DNA-encoded wild-type and mutant human tRNALys and tRNALeu(UUR)

Post-transcriptional modifications are characteristic features of tRNAs and have been shown in a number of cases to influence both their structural and functional properties, including structure stabilization, amino-acylation and codon recognition. We have developed an approach which allows the investigation of the post-transcriptional modification patterns of human mitochondrial wild-type and mutant tRNAs at both the qualitative and the quantitative levels. Specific tRNA species are long-term labeled in vivo with [32P]orthophosphate, isolated in a highly selective way, enzymatically digested to mononucleotides and then subjected to two-dimensional thin layer chromatographic analysis. The wild-type tRNALysand the corresponding tRNALyscarrying the A8344G mutation associated with the MERRF (Myoclonic Epilepsy with Ragged Red Fibers) syndrome exhibit the same modified nucleotides at the same molar concentrations. By contrast, a quantitatively different modification pattern was observed between the wild-type tRNALeu(UUR)and its counterpart carrying the A3243G mutation associated with the MELAS (Mitochondrial Myopathy, Encephalopathy with Lactic Acidosis and Stroke-like episodes) syndrome, the latter exhibiting a 50% decrease in m2G content. Complementary sequencing of tRNALeu(UUR)has allowed the localization of this modification at position 10 within the D-stem of the tRNA. The decreased level of this modification may have important implications for understanding the molecular mechanism underlying the MELAS-associated mitochondrial dysfunction.

Determination of the angle between the acceptor and anticodon stems of a truncated mitochondrial tRNA

Significant departures from the canonical (cloverleaf) secondary structure of transfer (t)RNAs can be found among the mitochondrial (m)tRNAs of higher metazoans; these mtRNAs thus pose a challenge to the concept of an invariant, L-shaped tertiary conformation for all tRNAs. For bovine mtRNASer(AGY), which lacks the entire "dihydrouridine" (dhU) arm, two distinct tertiary models have been proposed: the first model preserves the L-shaped conformation at the expense of overall size; the second model preserves the absolute distance between the 3' terminus and the anticodon loop, while allowing the acceptor-anticodon interstem angle to vary. We have tested the central predictions of these two models by performing a series of transient electric birefringence measurements on bovine mtRNASer(AGY) constructs in which the aminoacyl-acceptor and anticodon stems were each extended by approximately 70 bp. This mtRNA species is particularly amenable to analysis, since the native bovine (heart) mtRNA is completely unmodified outside of the anticodon loop. For magnesium ion concentrations above 1 mM, the interstem angle for the extended mtRNA, 120(+/-5) degrees, is approximately 50% larger than the corresponding angle for yeast tRNAPhe (70-80 degrees) under the same ionic conditions. Furthermore, the interstem angles of the two tRNAs exhibit strikingly different responses to the addition of magnesium ions: the interstem angle for yeast tRNAPhe is reduced by nearly 50 % upon addition of 2 mM magnesium ions, whereas the angle for mtRNASer(AGY) increases by about 10%. Our data thus support a central prediction of the second model; namely, that truncated mtRNAs will possess more open interstem angles. In addition, we demonstrate that birefringence amplitude data can be used to provide model-independent estimates for the interstem angles.

Higher-order structure and thermal instability of bovine mitochondrial tRNASerUGA investigated by proton NMR spectroscopy

Although mammalian mitochondrial serine-specific tRNA with the anticodon UGA (tRNASerUGA) appears to possess an almost normal cloverleaf secondary structure, it exhibits an extraordinarily low melting temperature (tm). An in vitro tRNASerUGA transcript without modified nucleosides had an even lower tm and slightly less hyperchromicity, but its tertiary structure was apparently very similar to that of the native counterpart judging from its aminoacylation activity and the body of experimental evidence so far obtained for canonical tRNAs. The transcript was therefore used to investigate the higher-order structure and thermal instability of tRNASerUGA. 1H-NMR analysis of the transcript showed that it takes a nearly L-shaped tertiary structure with similar tertiary base-pairings to those found in yeast tRNAPhe, which is representative of canonical tRNAs. However, magnesium ion titration revealed that Mg2+ affected the chemical shifts of the tRNASerUGA transcript differently than those of canonical tRNAs so far studied; the former was less sensitive toward Mg2+, especially in the D-arm region. This observation was confirmed by NMR analysis with paramagnetic manganese ion titration. Hill plots derived from the CD spectral changes caused by titration with Mg2+ suggested that the tRNASerUGA transcript had fewer Mg2+ binding sites than those of yeast tRNAPhe as well as its transcript, a finding that was consistent with the NMR data. We thus surmise that the thermal instability of both the transcript and tRNASerUGA itself originated from a reduction in the number of the divalent ion binding sites within the tRNA molecule. These results suggest a new type of thermal instability for mitochondrial tRNA.

Effect of a mutation in the anticodon of human mitochondrial tRNAPro on its post-transcriptional modification pattern

Although the gene sequences of all 22 tRNAs encoded in the human mitochondrial genome are known, little information exists about their sequences at the RNA level. This becomes a crucial limitation when searching for a molecular understanding of the growing number of maternally inherited human diseases correlated with point mutations in tRNA genes. Here we describe the sequence of human mt-tRNAPropurified from placenta. It shows absence of editing events in this tRNA and highlights the presence of eight post-transcriptional modifications. These include T54, never found so far in an animal mt-tRNA, and m1G37, a modification known to have fundamental functional properties in a number of canonical tRNAs. Occurrence of m1G37 was further investigated in an analysis of the substrate properties of in vitro transcripts of human mt-tRNAProtowards pure Escherichia coli methylguanosine transferase. This enzyme properly methylates G37 in mt-tRNA and is sensitive to the presence of a second G at position 36, neighboring the target nucleotide for methylation. Since mutation of nt 36 was shown to be correlated with myopathy, the potential consequences of non-modification or under-modification of mt-tRNA nucleotides in expression of the particular myopathy and of mitochondrial diseases in general are discussed.

Structural requirements for enzymatic formation of threonylcarbamoyladenosine (t6A) in tRNA: an in vivo study with Xenopus laevis oocytes

We have investigated the specificity of the eukaryotic enzymatic machinery that transforms adenosine at position 37 (3' adjacent to anticodon) of several tRNAs into threonylcarbamoyladenosine (t6A37). To this end, 28 variants of yeast initiator tRNAMet and yeast tRNAVal, devoid of modified nucleotide, were produced by in vitro transcription with T7 polymerase of the corresponding synthetic tRNA genes and microinjected into the cytoplasm of Xenopus laevis oocytes. Threonylcarbamoyl incorporation was analyzed in tRNA transcripts mutated in the anticodon loop by substitution, deletion, or Insertion of nucleotides, or in the overall 3D structure of the tRNA by altering critical tertiary interactions. Specifically, we tested the effects of altering ribonucleotides in the anticodon loop, changes of the loop size, perturbations of the overall tRNA 3D structure due to mutations disruptive of the tertiary base pairs, and truncated tRNAs. The results indicate that, in addition to the targeted A37, only U36 was absolutely required. However, A38 in the anticodon loop considerably facilitates the quantitative conversion of A37 into t6A37 catalyzed by the enzymes present in X. laevis. The anticodon positions 34 and 35 were absolutely "neutral" and can accept any of the four canonical nucleotides A, U, C, or G. The anticodon loop size may vary from six to eight nucleotides, and the anticodon stem may have one mismatch pair of the type AxC or GxU at location 30-40 without affecting the efficiency of t6A37 formation and still t6A37 is efficiently formed. Although threonylcarbamoylation of A37 occurred with tRNA having limited perturbations of 3D structure, the overall L-shaped architecture of the tRNA substrate was required for efficient enzymatic conversion of A37 to t6A37. These results favor the idea that unique enzymatic machinery located in the oocyte cytoplasm catalyzes the formation of t6A37 in all U36A37-containing tRNAs (anticodon NNU). Microinjection of the yeast tRNAMeti into the cytoplasm of X. laevis oocytes also revealed the enzymatic activities for several other nucleotide modifications, respectively m1Gg, m2G10, m(2)2G26, m7G46, D47, m5C48/49, and m1A58.

Usefulness of functional and structural solution data for the modeling of tRNA-like structures

Structures of large RNAs are not easily solved by X-ray crystallography or by NMR spectroscopy. This paper reviews the alternate methodology based on enzymatic and chemical mapping data collected on RNAs combined with graphical modeling for the construction of three-dimensional models. The different steps that lead to the establishment of the models are critically discussed. It is shown how the correctness of an RNA model can be strengthened by establishing correlations between the structure and the functionality of the molecule and its variants. Finally, the predictive potential of a model is discussed The approach is illustrated by results obtained on plant viral tRNA-like structures, and particularly on that of brome mosaic virus (BMV) RNA.

Yeast tRNA(Met) recognition by methionyl-tRNA synthetase requires determinants from the primary, secondary and tertiary structure: a review

The primordial role of the CAU anticodon in methionine identity of the tRNA has been established by others nearly a decade ago in Escherichia coli and yeast tRNA(Met). We show here that the CAU triplet alone is unable to confer methionine acceptance to a tRNA. This requires the contribution of the discriminatory base A73 and the non-anticodon bases of the anticodon loop. To better understand the functional communication between the anticodon and the active site, we analysed the binding and aminoacylation of tRNA(Met) based anticodon and acceptor-stem minihelices and of tRNA(Met) chimeras where the central core region of yeast tRNA(Met) is replaced by that of unusual mitochondrial forms lacking either a D-stem or a T-stem. These studies suggest that the high selectivity of the anticodon bases in tRNA(Met) implies the L-conformation of the tRNA and the presence of a D-stem. The importance of a L-structure for recognition of tRNA(Met) was also deduced from mutations of tertiary interactions known to play a general role in tRNA(Met) folding.

Nucleotide sequence and gene organization of the starfish Asterina pectinifera mitochondrial genome

The 16,260-bp mitochondrial DNA (mtDNA) from the starfish Asterina pectinifera has been sequenced. The genes for 13 proteins, two rRNAs and 22 tRNAs are organized in an extremely economical fashion, similar to those of other animal mtDNAs, with some of the genes overlapping each other. The gene organization is the same as that for another echinoderm, sea urchin, except for the inversion of a 4.6-kb segment that contains genes for two proteins, 13 tRNAs and the 16S rRNA. Judging from the organization of the protein coding genes, mammalian mtDNAs resemble the sea urchin mtDNA more than that of the starfish. The region around the 3' end of the 12S rRNA gene of the starfish shows a high similarity with those for vertebrates. This region encodes a possible stem and loop structure; similar potential structures occur in this region of vertebrate mtDNAs and also in nonmitochondrial small subunit rRNA. A similar stem and loop structure is also found at the 3' end of the 16S rRNA genes in A. pectinifera, in another starfish Pisaster ochraceus, in vertebrates and in Drosophila, but not in sea urchins. The full sequence data confirm the presumption that AGA/AGG, AUA and AAA codons, respectively, code for serine, isoleucine, and asparagine in the starfish mitochondria, and that AGA/AGG codons are read by tRNA(GCUSer), which possesses a truncated dihydrouridine arm, that was previously suggested from a partial mtDNA sequence. The structural characteristics of tRNAs and possible mechanisms for the change in the mitochondrial genetic code are also discussed.

Influence of modified nucleosides on tRNA structure as probed by two plant nucleases

Two new enzymatic probes have been used for structural investigations of native and unmodified transcript tRNA molecules. Both probes were single-strand-specific nucleases isolated from higher plants. The results obtained after enzymatic hydrolysis of tRNAs support the earlier hypothesis that posttranscriptional modifications in tRNA help to stabilize its structure and make it more rigid.

The mitochondrial ribosomal RNA genes of the nematodes Caenorhabditis elegans and Ascaris suum: consensus secondary-structure models and conserved nucleotide sets for phylogenetic analysis

The small- and large-subunit mitochondrial ribosomal RNA genes (mt-s-rRNA and mt-l-rRNA) of the nematode worms Caenorhabditis elegans and Ascaris suum encode the smallest rRNAs so far reported for metazoa. These size reductions correlate with the previously described, smaller, structurally anomalous mt-tRNAs of C. elegans and A. suum. Using primer extension analysis, the 5' end nucleotides of the mt-s-rRNA and mt-l-rRNA genes were determined to be adjacent to the 3' end nucleotides of the tRNA(Glu) and tRNA(His) genes, respectively. Detailed, consensus secondary-structure models were constructed for the mt-s-rRNA genes and the 3' 64% of mt-l-rRNA genes of the two nematodes. The mt-s-rRNA secondary-structure model bears a remarkable resemblance to the previously defined universal core structure of E. coli 16S rRNA: most of the nucleotides that have been classified as variable or semiconserved in the E. coli model appear to have been eliminated from the C. elegans and A. suum sequences. Also, the secondary structure model constructed for the 3' 64% of the mt-l-rRNA is similar to the corresponding portion of the previously defined E. coli 23S rRNA core secondary structure. The proposed C. elegans/A. suum mt-s-rRNA and mt-l-rRNA models include all of the secondary-structure element-forming sequences that in E. coli rRNAs contain nucleotides important for A-site and P-site (but not E-site) interactions with tRNAs. Sets of apparently homologous sequences within the mt-s-rRNA and mt-l-rRNA core structures, derived by alignment of the C. elegans and A. suum mt-rRNAs to the corresponding mt-rRNAs of other eukaryotes, and E. coli rRNAs were used in maximum-likelihood analyses. The patterns of divergence of metazoan phyla obtained show considerable agreement with the most prevalent metazoan divergence patterns derived from more classical, morphological, and developmental data.

Fitting the structurally diverse animal mitochondrial tRNAs(Ser) to common three-dimensional constraints

We propose three-dimensional models for animal mitochondrial (amt) tRNAs lacking the D-domain based on consideration of universal constraints on tRNA to maintain functionality. The available tRNA sequences are classified into two groups, and distinct models are proposed for both classes derived from common structural features. The distance between the anticodon and the acceptor stem is comparable in the models and corresponds to that observed in conventional tRNAs. This fact averts the problem of how a shorter mitochondrial tRNA could function within the context of a protein synthesis machinery suited to full-sized tRNAs. In the models, the angle which defines the relationship between the helical domains composed of the acceptor/T-stem and the anticodon/D-stem is greater than in conventional tRNAs. These structures resemble more a "boomerang" than an "L". However, even in the boomerang model, the inner surface of tRNA would be sufficiently uncluttered to avoid steric clashes when two tRNA molecules cohabit the ribosome.

Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics

Mitochondrial DNA sequences are often used to construct molecular phylogenetic trees among closely related animals. In order to examine the usefulness of mtDNA sequences for deep-branch phylogenetics, genes in previously reported mtDNA sequences were analyzed among several animals that diverged 20-600 million years ago. Unambiguous alignment was achieved for stem-forming regions of mitochondrial tRNA genes by virtue of their conservative secondary structures. Sequences derived from stem parts of the mitochondrial tRNA genes appeared to accumulate much variation linearly for a long period of time: nearly 100 Myr for transition differences and more than 350 Myr for transversion differences. This characteristic could be attributed, in part, to the structural variability of mitochondrial tRNAs, which have fewer restrictions on their tertiary structure than do nonmitochondrial tRNAs. The tRNA sequence data served to reconstruct a well-established phylogeny of the animals with 100% bootstrap probabilities by both maximum parsimony and neighbor-joining methods. By contrast, mitochondrial protein genes coding for cytochrome b and cytochrome oxidase subunit I did not reconstruct the established phylogeny or did so only weakly, although a variety of fractions of the protein gene sequences were subjected to tree-building. This discouraging phylogenetic performance of mitochondrial protein genes, especially with respect to branches originating over 300 Myr ago, was not simply due to high randomness in the data. It may have been due to the relative susceptibility of the protein genes to natural selection as compared with the stem parts of mitochondrial tRNA genes. On the basis of these results, it is proposed that mitochondrial tRNA genes may be useful in resolving deep branches in animal phylogenies with divergences that occurred some hundreds of Myr ago. For this purpose, we designed a set of primers with which mtDNA fragments encompassing clustered tRNA genes were successfully amplified from various vertebrates by the polymerase chain reaction.

Non-canonical substrates of aminoacyl-tRNA synthetases: the tRNA-like structure of brome mosaic virus genomic RNA

A 3-D model of the tyrosylable tRNA-like domain of the genomic brome mosaic virus RNAs was built by computer modelling based on solution probing of the molecule with different chemical and enzymatic reagents. This model encompasses four major structural domains, including two peculiar substructures oriented perpendicularly and mimicking a tRNA structure, and a fifth domain which makes the connection with the rest of the viral RNA. After recalling the different steps that led to the present structural knowledge of the BMV tRNA-like domain, we review its novel structural features revealed by the modelling and that did not appear in older versions of 3-D models of this structure. These features comprise additional base-pairs, hairpin loops, new tertiary long-range interactions, and a second pseudoknot. The main goal of this paper is to strengthen the validity of the model by establishing correlations between the putative 3-D conformation and the functional properties of the domain. For that, we show how the present structural model rationalises mutagenic and footprinting data that have established the importance of specific regions of the RNA for its recognition and aminoacylation by yeast tyrosyl-tRNA synthetase. We discuss further how the model corroborates mutational analyses performed to understand recognition of this RNA domain by the (ATP,CTP):tRNA nucleotidyl-transferase and by the viral replicase. The published mutants of the BMV tRNA-like domain fall into two classes. In one class, the mutants leave unchanged the overall architecture of the molecule, thereby affecting functions directly. In the second class, the overall architecture of the mutants is perturbed, and thus functions are affected indirectly.

The TYMV tRNA-like structure

The genomic RNA from turnip yellow mosaic virus presents a 3'-end functionally and structurally related to tRNAs. This report summarizes our knowledge about the peculiar structure of the tRNA-like domain and its interaction with tRNA specific proteins, like RNAse P, tRNA nucleotidyl-transferase, aminoacyl-tRNA synthetases, and elongation factors. It discusses also the biological role of this structure in the viral life cycle. A brief survey of our knowledge of other tRNA mimicries in biological systems, as well as their relevance for understanding canonical tRNA, will also be presented.

The translation system of mammalian mitochondria

Oligoribonucleotides and mRNA were used to define properties of the bovine mitoribosomal mRNA binding site. The RNA binding domain on the 28 S subunit spans approx. 80 nucleotides of the template, based on ribosome protection experiments, but the major interaction with the ribosome occurs over a 30 nucleotide stretch. The binding site for E. coli IF3 is conserved in bovine mitoribosomes, but mitochondrial factors appear essential for proper interaction of mRNA with mitoribosomes. The small subunit of bovine mitoribosomes contains a high-affinity binding site for guanyl nucleotides, further indication of specialized mechanisms for initiation complex formation and function of mammalian mitochondrial ribosomes.

Nobel lecture. Enzymatic cleavage of RNA by RNA

The discovery and characterization of the catalytic RNA subunit of the enzyme ribonuclease P of Escherichia coli is described.

RNA pseudoknots. Stability and loop size requirements

The effects of ionic conditions, loop size and loop sequence on the formation of pseudoknots by RNA oligonucleotides have been investigated using biochemical and biophysical methods. An oligonucleotide with the sequence 5' GCGAUUUCUGACCGCUUUUUUGUCAG 3' and oligonucleotides with variations in the sequences of the two loop regions, denoted by bold face type, were studied. Each sequence with the potential to form a pseudoknot can also form two stable hairpins. The pseudoknot structure is stabilized relative to the hairpins by addition of Mg2+. Even in the presence of Mg2+, the pseudoknots formed by the sequences investigated are only marginally more stable (1.5 to 2 kcal mol-1 in free energy at 37 degrees C) than either of the constituent hairpins. The pseudoknot structure is the stable conformation in the presence of Mg2+ when the first loop region is at least three nucleotides and the second is at least four nucleotides. Further deletion of nucleotides from the loop regions stabilizes possible hairpin structures relative to the pseudoknot and equilibria among secondary and tertiary structures result.

Secondary structure of the central region of bacteriophage MS2 RNA. Conservation and biological significance

The RNA of the Escherichia coli RNA phages is highly structured with 75% of the nucleotides estimated to take part in base-pairing. We have used enzymatic and chemical sensitivity of nucleotides, phylogenetic sequence comparison and the phenotypes of constructed mutants to develop a secondary structure model for the central region (900 nucleotides) of the group I phage MS2. The RNA folds into a number of, mostly irregular, helices and is further condensed by several long-distance interactions. There is substantial conservation of helices between the related groups I and II, attesting to the relevance of discrete RNA folding. In general, the secondary structure is thought to be needed to prevent annealing of plus and minus strand and to confer protection against RNase. Superimposed, however, are features required to regulate translation and replication. The MS2 RNA section studied here contains three translational start sites, as well as the binding sites for the coat protein and the replicase enzyme. Considering the density of helices along the RNA, it is not unexpected to find that all these sites lie in helical regions. This fact, however, does not mean that these sites are recognized as secondary structure elements by their interaction partners. This holds true only for the coat protein binding site. The other four sites function in the unfolded state and the stability of the helix in which they are contained serves to negatively control their accessibility. Mutations that stabilize helices containing ribosomal binding sites reduce their efficiency and vice versa. Comparison of homologous helices in different phage RNAs indicates that base substitutions have occurred in such a way that the thermodynamic stability of the helix is maintained. The evolution of individual helices shows several distinct size-reduction patterns. We have observed codon deletions from loop areas and shortening of hairpins by base-pair deletions from either the bottom, the middle or the top of stem structures. Evidence for the coaxial stacking of some helical segments is discussed.

Genomic organisation of nuclear tRNAGly and tRNALeu genes in Trypanosoma brucei

We have isolated a 0.3-kb HaeIII restriction fragment from Trypanosoma brucei which contains two tRNA genes. Secondary structure models predict that the two genes identified encode tRNA molecules which specify glycine (anticodon UCC) and leucine (anticodon CAG). The two genes are separated by 86 nucleotides, transcribed in the same direction and contain features of conventional RNA polymerase III transcription units. Southern blot analysis indicates the presence of multicopy tRNa gene families in T. brucei.

Computer modeling from solution data of spinach chloroplast and of Xenopus laevis somatic and oocyte 5 S rRNAs

Detailed atomic models of a eubacterial 5 S rRNA (spinach chloroplast 5 S rRNA) and of a eukaryotic 5 S rRNA (somatic and oocyte 5 S rRNA from Xenopus laevis) were built using computer graphic. Both models integrate stereochemical constraints and experimental data on the accessibility of bases and phosphates towards several structure-specific probes. The base sequence was first inserted on to three-dimensional structural fragments picked up in a specially devised databank. The fragments were modified and assembled interactively on an Evans & Sutherland PS330. Modeling was finalized by stereochemical and energy refinement. In spite of some uncertainty in the relative spatial orientation of the substructures, the broad features of the models can be generalized and several conclusions can be reached: (1) both models adopt a distorted Y-shape structure, with helices B and D not far from colinearity; (2) no tertiary interactions exist between loop c and region d or loop e; (3) the internal loops, in particular region d, contain several non-canonical base-pairs of A.A, U.U and A.G types; (4) invariant residues appear to be more important for protein or RNA binding than for maintaining the tertiary structure. The models are corroborated by footprinting experiments with ribosomal proteins and by the analysis of various mutants. Such models help to clarify the structure-function relationship of 5 S rRNA and are useful for designing site-directed mutagenesis experiments.