A manually curated database of tetrapod mitochondrially encoded tRNA sequences and secondary structures

Background  

Mitochondrial tRNAs have been the subject of study for structural biologists interested in their secondary structure characteristics, evolutionary biologists have researched patterns of compensatory and structural evolution and medical studies have been directed towards understanding the basis of human disease. However, an up to date, manually curated database of mitochondrially encoded tRNAs from higher animals is currently not available.     

PLMItRNA, a database for higher plant mitochondrial tRNAs and tRNA genes.

The PLMItRNA database contains information and multialignments of tRNA genes and molecules detected in higher plant mitochondria. It has been developed from a previous compilation of higher plant mitochondrial tRNA genes [Sagliano,A., Volpicella,M., Gallerani,R. and Ceci,L.R. (1998) Nucleic Acids Res., 26, 154-155] and implemented with data and sequences of tRNA molecules retrieved from the literature. The current version of the database reports information on 171 genes and 16 tRNA molecules from 24 plants. PLMItRNA is accessible via WWW at http://bio-www.ba.cnr.it:8000/srs/     

PLMItRNA, a database for mitochondrial tRNA genes and tRNAs in photosynthetic eukaryotes

The PLMItRNA database for mitochondrial tRNA molecules and genes in VIRIDIPLANTAE: (green plants) [Volpetti,V., Gallerani,R., DeBenedetto,C., Liuni,S., Licciulli,F. and Ceci,L.R. (2000) Nucleic Acids Res., 28, 159-162] has been enlarged to include algae. The database now contains 436 genes and 16 tRNA entries relative to 25 higher plants, eight green algae, four red algae (RHODOPHYTAE:) and two STRAMENOPILES: The PLMItRNA database is accessible via the WWW at http://bio-www.ba.cnr.it:8000/PLMItRNA.     

Identification and characterization of mammalian mitochondrial tRNA nucleotidyltransferases

The CCA-adding enzyme (ATP:tRNA adenylyltransferase or CTP:tRNA cytidylyltransferase (EC )) generates the conserved CCA sequence responsible for the attachment of amino acid at the 3' terminus of tRNA molecules. It was shown that enzymes from various organisms strictly recognize the elbow region of tRNA formed by the conserved D- and T-loops. However, most of the mammalian mitochondrial (mt) tRNAs lack consensus sequences in both D- and T-loops. To characterize the mammalian mt CCA-adding enzymes, we have partially purified the enzyme from bovine liver mitochondria and determined cDNA sequences from human and mouse dbESTs by mass spectrometric analysis. The identified sequences contained typical amino-terminal peptides for mitochondrial protein import and had characteristics of the class II nucleotidyltransferase superfamily that includes eukaryotic and eubacterial CCA-adding enzymes. The human recombinant enzyme was overexpressed in Escherichia coli, and its CCA-adding activity was characterized using several mt tRNAs as substrates. The results clearly show that the human mt CCA-adding enzyme can efficiently repair mt tRNAs that are poor substrates for the E. coli enzyme although both enzymes work equally well on cytoplasmic tRNAs. This suggests that the mammalian mt enzymes have evolved so as to recognize mt tRNAs with unusual structures.     

Downstream elements of mammalian pre-mRNA polyadenylation signals: primary,secondary and higher-order structures

Primary, secondary and higher-order structures of downstream elements of mammalian pre-mRNA polyadenylation signals [poly(A) signals] are re viewed. We have carried out a detailed analysis on our database of 244 human pre-mRNA poly(A) signals in order to characterize elements in their downstream regions. We suggest that the downstream region of the mammalian pre-mRNA poly(A) signal consists of various simple elements located at different distances from each other. Thus, the downstream region is not described by any precise consensus. Searching our database, we found that ~80% of pre-mRNAs with the AAUAAA or AUUAAA core upstream elements contain simple downstream elements, consisting of U-rich and/or 2GU/U tracts, the former occurring ~2-fold more often than the latter. Approximately one-third of the pre-mRNAs analyzed here contain sequences that may form G-quadruplexes. A substantial number of these sequences are located immediately downstream of the poly(A) signal. A possible role of G-rich sequences in the polyadenylation process is discussed. A model of the secondary structure of the SV40 late pre-mRNA poly(A) signal downstream region is presented.     

Carboxylic ester hydrolases: Classification and database derived from their primary,secondary, and tertiary structures

We classified the carboxylic ester hydrolases (CEHs) into families and clans by use of multiple sequence alignments, secondary structure analysis, and tertiary structure superpositions. Our work for the first time has fully established their systematic structural classification. Family members have similar primary, secondary, and tertiary structures, and their active sites and reaction mechanisms are conserved. Families may be gathered into clans by their having similar secondary and tertiary structures, even though primary structures of members of different families are not similar. CEHs were gathered from public databases by use of Basic Local Alignment Search Tool (BLAST) and divided into 91 families, with 36 families being grouped into five clans. Members of one clan have standard α/β‐hydrolase folds, while those of other two clans have similar folds but with different sequences of their β‐strands. The other two clans have members with six‐bladed β‐propeller and three‐α‐helix bundle tertiary structures. Those families not in clans have a large variety of structures or have no members with known structures. At the time of writing, the 91 families contained 321,830 primary structures and 1378 tertiary structures. From these data, we constructed an accessible database: CASTLE (CArboxylic eSTer hydroLasEs, http://www.castle.cbe.iastate.edu ).     

MamMiBase: a mitochondrial genome database for mammalian phylogenetic studies

MamMiBase, the mammalian mitochondrial genome database, is a relational database of complete mitochondrial genome sequences of mammalian species. The database is useful for phylogenetic analysis, since it allows a ready retrieval of nucleotide and aminoacid individual alignments, in three different formats (NEXUS for PAUP program, for MEGA program and for PHYLIP program) of the 13 protein coding mitochondrial genes. The user may download the sequences that are useful for him/her based on their parameters values, such as sequence length, p-distances, base content, transition transversion ratio, gamma, which are also given by MamMiBase. A simple phylogenetic tree (neighbor-joining tree with Jukes Cantor distance) is also available for download, useful for parameter calculations and other simple tasks. AVAILABILITY: MamMiBase is available at http://www.mammibase.lncc.br     

The primary structure of yeast mitochondrial tyrosine tRNA

The mitochondrial tyrosine tRNA from Saccharomyces cerevisiae has been sequenced. It has two interesting structural features: (i) it lacks two semi-invariant purine residues in the D-loop which are involved in tertiary interactions in the yeast cytoplasmic tRNAPhe; (ii) it has a large variable loop and therefore resembles procaryotic tRNAsTyr rather than eucaryotic cytoplasmic ones.     

PLMItRNA,a database on the heterogeneous genetic origin of mitochondrial tRNA genes and tRNAs in photosynthetic eukaryotes

The updated version of PLMItRNA reports information and multialignments on 609 genes and 34 tRNA molecules active in the mitochondria of Viridiplantae (27 Embryophyta and 10 Chlorophyta), and photosynthetic algae (one Cryptophyta, four Rhodophyta and two Stramenopiles). Colour-code based tables reporting the different genetic origin of identified genes allow hyper-textual link to single entries. Promoter sequences identified for tRNA genes in the mitochondrial genomes of Angiospermae are also reported. The PLMItRNA database is accessible at http://bighost.area.ba.cnr.it/PLMItRNA/.     

Nuclear and mitochondrial revertants of a yeast mitochondrial tRNA mutant

Summary We isolated revertants capable of respiration from the respiratory deficient yeast mutant, FF1210-6C/ 170, which displays greatly decreased mitochondrial protein synthesis due to a single base substitution at the penultimate base of the tRNA^Asp gene on mitochondrial (mt) DNA. Three classical types of revertant were identified: (1) same-site revertants; (2) intragenic revertants which restore the base pairing in the acceptor stem of the mitochondrial tRNA^Asp; and (3) extragenic suppressors located in nuclear DNA. In addition a fourth type of revertant was identified in which the mutant tRNA^Asp is amplified due to the maintenance of both the original mutant mtDNA and a modified form of the mutant mtDNA in which only a small region around the tRNA^Asp gene is retained and amplified. The latter form resembles the mtDNA in vegetative petite (rho ^-) strains which normally segregates rapidly from the wild-type mtDNA. Each revertant type was characterized genetically and by both DNA sequence analysis of the mitochondrial tRNA^Asp gene and analysis of the quantity and size of RNA containing the tRNA^Asp sequence. These results indicate that the mitochondrial tRNA^Asp of the mutant retains a low level of activity and that the presence of the terminal base pair in tRNA^Asp is a determinant of both tRNA^Asp function and the maintenance of wild-type levels of tRNA^Asp.     

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.     

Predicting sub-cellular localization of tRNA synthetases from their primary structures

Since endo-symbiotic events occur, all genes of mitochondrial aminoacyl tRNA synthetase (AARS) were lost or transferred from ancestral mitochondrial genome into the nucleus. The canonical pattern is that both cytosolic and mitochondrial AARSs coexist in the nuclear genome. In the present scenario all mitochondrial AARSs are nucleus-encoded, synthesized on cytosolic ribosomes and post-translationally imported from the cytosol into the mitochondria in eukaryotic cell. The site-based discrimination between similar types of enzymes is very challenging because they have almost same physico-chemical properties. It is very important to predict the sub-cellular location of AARSs, to understand the mitochondrial protein synthesis. We have analyzed and optimized the distinguishable patterns between cytosolic and mitochondrial AARSs. Firstly, support vector machines (SVM)-based modules have been developed using amino acid and dipeptide compositions and achieved Mathews correlation coefficient (MCC) of 0.82 and 0.73, respectively. Secondly, we have developed SVM modules using position-specific scoring matrix and achieved the maximum MCC of 0.78. Thirdly, we developed SVM modules using N-terminal, intermediate residues, C-terminal and split amino acid composition (SAAC) and achieved MCC of 0.82, 0.70, 0.39 and 0.86, respectively. Finally, a SVM module was developed using selected attributes of split amino acid composition (SA-SAAC) approach and achieved MCC of 0.92 with an accuracy of 96.00%. All modules were trained and tested on a non-redundant data set and evaluated using fivefold cross-validation technique. On the independent data sets, SA-SAAC based prediction model achieved MCC of 0.95 with an accuracy of 97.77%. The web-server 'MARSpred' based on above study is available at http://www.imtech.res.in/raghava/marspred/.     

Dimerization of a pathogenic human mitochondrial tRNA

Mutations of human mitochondrial transfer RNA (tRNA) are implicated in a variety of multisystemic diseases. The most prevalent pathogenic mitochondrial mutation is the A3243G substitution within the gene for tRNA(Leu(UUR)). Here we describe the pronounced structural change promoted by this mutation. The A3243G mutation induces the formation of a tRNA dimer that strongly self-associates under physiological conditions. The dimerization interface in the mutant tRNA is a self-complementary hexanucleotide in the D-stem, a particularly weak structural element within tRNA(Leu(UUR)). Aminoacylation of the A3243G mutant is significantly attenuated, and mutational studies indicate that dimerization is partially responsible for the observed loss of function. The disruption of a conserved tertiary structural contact also contributes to the functional defect. The pathogenic mutation is proposed to interfere with the cellular function of human mitochondrial tRNA(Leu(UUR)) by destabilizing the native structure and facilitating the formation of a dimeric complex with low biological activity.     

Bleomycin cleaves DNA depending on DNA primary, secondary, and tertiary structures

The cleavage by bleomycin-Fe(II) complex in the presence of dithiothreitol was investigated by using 3'- or 5'-end-labeled DNA containing the region of the bacteriophage G4 origin of complementary strand synthesis as substrates. Bleomycin cleaved single-stranded DNA substrates preferentially at inverted repeat sequences, which potentially form stem-and-loop structures, in addition to the primary sequence specificity previously reported. DNA sequences preferentially cleaved in the double-stranded substrate were resistant when they lay outside the stem regions. These results suggest the formation of three predicted stem-and-loop structures and other possible secondary structures near the replication origin. Changes of the degree of bleomycin-induced DNA cleavage in a NaCl concentration between 0 and 50 mM suggest that a subtle change of ionic conditions within the double helix, or of DNA conformation, or of both, may occur at 0-50 mM NaCl. Bleomycin appears to be a useful reagent for analyzing secondary and tertiary structures of DNA.