|
|
||||||||
1 Department of Applied Chemistry, Faculty of Engineering, Ehime University, Bunkyo 3, Matsuyama 790-8577, Japan
2 Department of Chemistry and Biotechnology, Graduate School of Engineering, University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan
3 Department of Molecular Biology, School of Life Science, Tokyo University of Pharmacy and Life Science, Horinouchi 1432-1, Hachioji 192-0392, Japan
4 Venture Business Laboratory, Ehime University, Bunkyo 3, Matsuyama 790-8577, Japan
5 Cell-free Science and Technology Research Center, Ehime University, Matsuyama, 790-8577, Japan
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Recent computational research has reported that bacterial TrmD proteins and some RNA methyltransferases such as SpoU (TrmH) (Persson et al. 1997; CavailIé et al. 1999; Hori et al. 2002, 2003) might share a common evolutionary origin and form a single "SPOUT" (SpoU-TrmD) super-family (Anantharaman et al. 2002). Work in this laboratory over a number of years has helped delineate the properties of tRNA (Gm18) methyltransferase [SpoU (trmH)] (Hori et al. 1989, 1998, 2002, 2003; Matsumoto et al. 1990; Nureki et al. 2004; Watanabe et al. 2005). TrmH catalyzes methyl-transfer from AdoMet to the 2'-OH of ribose of the conserved G18 (Hori et al. 1998). Recently, we determined the crystal structure of T. thermophilus TrmH (Nureki et al. 2004), and three groups independently reported the structures of TrmD proteins from Haemophilus influenzae (Ahn et al. 2003), E. coli (Elkins et al. 2003) and Aquifex aeolicus (Liu et al. 2003). These structural studies reveal that both TrmH and TrmD proteins have a deep trefoil knot structure (Nureki et al. 2002). Thus, structural similarities of these proteins were established.
To investigate the origin of SPOUT proteins, we have studied the RNA modification enzymes from A. aeolicus (Hori et al. 2003), an organism which was isolated from a hot spring in Yellowstone National Park, and which can grow at nearly 95 °C. The 16S rRNA gene of Aquifex has been analyzed from the perspective of molecular evolution, and it has been suggested that this bacterium may be the earliest diverging eubacterium (Burggraf et al. 1992). The complete genome of this organism was determined in 1998 (Deckert et al. 1998). Because the m1G37 modification prevents frameshift errors during protein synthesis on the ribosome (Björk et al. 1989; Hagervall et al. 1990; 1993; Farabaugh & Björk 1999), this modification is essentially required in all living organisms (Björk et al. 2001). In fact, trmD genes have been reported from all of the eubacterial genomes analyzed. In this paper, we demonstrate the substrate tRNA recognition mechanism of A. aeolicus TrmD, which catalyzes the methyl-transfer to G37 not only in a G36G37 sequence but also in an A36G37 sequence. Further, the molecular evolution of this modification is discussed.
| Results |
|---|
|
|
|---|
To characterize A. aeolicus TrmD, we amplified a DNA fragment containing the ORF by PCR and expressed it as a recombinant protein in E. coli. The yield of protein from this expression system in the E. coli BL21(DE3) strain was 100 µg per 1 L culture. In order to improve the protein yield, we exchanged the DNA fragment between SacI and XhoI sites for an artificial DNA fragment containing codons preferentially used in E. coli (see Experimental procedures). The resultant plasmid was introduced into E. coli BL21(DE3)-Codonplus-RIL strain. This alteration resulted in a remarkable fivefold improvement in the yield of recombinant protein.
We purified the TrmD protein to homogeneity. The purified protein is shown on a 15% SDS-polyacrylamide gel in Fig. 1.
|
To investigate the substrate RNA specificity of A. aeolicus TrmD, we analyzed methyl-transfer to various tRNA transcripts (Fig. 2). tRNA transcript (0.3 A260 units) was incubated with 300 ng of the purified TrmD protein (5 pmol of the dimer enzyme) and 38 µM [methyl-14C]-AdoMet for 10 min at 60 °C. The RNA was loaded on to a 10% polyacrylamide gel (7 M urea). The gel was stained with methylene blue (Fig. 3, left) to visualize RNA and dried. The 14C-incorporation into the RNA was monitored using a Fuji Photo Film BAS2000 imaging analyzer (Fig. 3, right). The kinetic parameters measured at 55 °C are given in Table 1. As expected, A. aeolicus TrmD methylated tRNA transcripts possessing a G36G37 sequence in their anti-codon loop, a result consistent with the recognition mechanism of E. coli TrmD (Holmes et al. 1992, 1995; Redlak et al. 1997) (Fig. 3, lanes 17, and Table 1). We also tested tRNA transcripts possessing an A36G37 sequence (yeast tRNAPhe (GAA), Haloferax volcanii tRNALeu (CAA), etc.) (see Fig. 2). Surprisingly, methyl-transfer to these tRNAs was observed (Fig. 3, lanes 812, and Table 1). The methyl-transfer activity to these tRNAs was considerably stronger compared to that to tRNAs possessing a G36G37 sequence (Fig. 3, lanes 17, and Table 1). In Table 1, E. coli tRNALeu, yeast tRNAPhe, H. volcanii tRNATrp and tRNATyr were methylated slower than the other tRNA transcripts. These tRNA transcripts have many AU pairs and/or mismatch pairs (GU and AC pairs) in the stem structure as compared to the other transcripts. Because we measured the kinetic parameters at 55 °C, the stability of the stem structure may influence the methyl acceptance activities of these transcripts. We also tested yeast tRNAAsp (GUC) and H. volcanii tRNAGlu (UUC) possessing a C36G37 sequence, yeast mutant tRNAPhe (GAU) possessing an U36G37 sequence, and E. coli tRNASer (UGA) possessing an A36A37 sequence (Fig. 2). These tRNA species were not methylated at all (Fig. 3, lanes 1316, and Table 1). From these results, we conclude that A. aeolicus TrmD catalyzes the methyl-transfer to tRNA possessing a purine36G37 sequence.
|
|
|
In order to check the modified nucleotide, we first performed the 5'-mononucleotide analysis by two-dimensional thin layer chromatography (Keith 1995). Yeast tRNAPhe was methylated with [methyl-14C]-AdoMet, and then completely digested with nuclease P1. The 5'-mononucleotides were separated by two-dimensional TLC and the mobility of the labeled spot was compared to those of cold markers: unlabeled pA, pU, pG, and pC were separated on the same thin-layer plate and marked by UV (254 nm) irradiation. As a result, [14C]-pm1G could be identified as a single labeled spot (data not shown).
To identify the modified nucleoside and the position of the methylated site precisely, we employed LC/MS using ESI/iontrap mass spectrometry (Figs 4 and 5). In these experiments, we used yeast tRNAPhe transcript, because the elution profile of the yeast tRNAPhe fragments by RNase T1 digestion was investigated in our previous studies (Hori et al. 2002; Okamoto et al. 2004). Fig. 4 shows the result of the nucleoside analysis. Yeast tRNAPhe transcript was methylated by the TrmD protein with cold AdoMet at 60 °C for 4 h and then completely digested with nuclease P1 and bacterial alkaline phosphatase. The m1G nucleoside eluted at 25 min (Fig. 4A), and the mass signals of the protonated molecule of m1G nucleoside (m/z = 298), the fragment ion derived from the base moiety (m/z = 166) and the sodium adduct (m/z = 320) were clearly identified (Fig. 4B).
|
|
Methyl group acceptance activities of the truncated RNA molecules and replacements of the anti-codon stem sequence
To search the additional identity elements recognized by A. aeolicus TrmD, we prepared nine truncated RNA molecules of yeast tRNAPhe (Fig. 6). The methyl acceptance activities of these truncated RNA molecules were tested at 37 and 55 °C. Because both experiments showed the same tendency, the results at 55 °C are shown in Fig. 6. Unexpectedly, the micro helix RNA corresponding to the anti-codon arm was efficiently methylated (Fig. 6B). The kinetic study showed that the Km value for this micro helix was very large compared to that for the full-length yeast tRNAPhe (Table 2). However, because Vmax value for this micro helix was increased, the initial velocity of the micro helix was around 80% of the full-length transcript. When the anti-codon stem structure was disrupted, the methyl group acceptance activity was completely lost (Fig. 6C). This result clearly showed that the anti-codon stem structure is one of the additional identity elements. Furthermore, it is also clear that the A. aeolicus TrmD does not require the L-shaped tRNA structure. Moreover, the kinetic study showed that the deletion of the D-arm structure increased the Vmax value (transcript D, E, and G in Table 2). Indeed, these transcripts were methylated faster than the full-length yeast tRNAPhe transcript (Fig. 6 and Table 2). Because the existence of the nucleotides corresponding to positions 1925 decreased the Vmax value (transcript F in Table 2), this region seems to affect the methyl group acceptance activity. In the case of E. coli TrmD, it has been reported that the deletions of both T- and D-arms decreased the methyl group acceptance activity severely: the initial velocity to the mutant transcript is 1/100 of that of the full-length transcript (Redlak et al. 1997). This transcript corresponds to the transcript E in the current study (Fig. 6 and Table 2). Thus, the effect of the D-stem structure seems to differ between A. aeolicus and E. coli TrmD proteins.
|
|
|
To address the structural differences of eubacterial TrmD proteins, we performed phylogenetic tree analysis. Because the phylogenetic tree in Fig. 8 does not have the root, the tree does not represent the precise evolutionary process followed by TrmD proteins: it simply shows the resemblances of amino acid sequences. This analysis demonstrates that A. aeolicus TrmD is considerably different from the other TrmD proteins of eubacteria. For example, the value of the local bootstrap probability between A. aeolicus and Thermotoga maritima TrmD proteins was calculated as 67, even though T. maritima TrmD is the most similar protein to A. aeolicus TrmD (Fig. 8). As found in Fig. 8, the A. aeolicus TrmD protein was grouped with T. martima and M. capricolum TrmD proteins. In small subunit rRNA based molecular phylogenetic analyses, T. martima and A. aeolicus are placed as the most basal species of eubacteria (Woese et al. 1990; Burggraf et al. 1992). However, the M. capricolum is a member of Firmicutes, so that its phylogenetic position is different from the result of the small subunit rRNA analysis. This result suggests that the evolutionary rate of the M. capricolum TrmD protein may be accelerated. Furthermore, recently, another approach to determine the relative branching order of bacterial divisions has been reported: Griffiths & Gupta (2001) used signature sequences in the Hsp60, Hsp70, CTP synthase and alanyl-tRNA synthetase to determine the phylogenetic placement of Fibrobacter succinogenes. In their work, A. aeolicus was classified in the Cytophaga-Flavobacteria-Bacteroides group based on their signature sequence analyses. Thus, their results suggest that A. aeolicus is not an earliest diverged eubacterium. More recently, the relationship of the species in the Aquificae has been investigated by inserts and deletions (so called "indels") and signature sequences (Griffiths & Gupta 2006).
|
Distribution of A36 recognition by TrmD
Recently, Brulé et al. (2004) reported that E. coli TrmD catalyzes methyl-transfer to transcripts of yeast tRNAPhe, tRNAAsp A36 mutant and human mitochondria tRNAPro A36 mutant, although the initial velocities of the methyl-transfer to these transcripts were slow compared to the E. coli tRNALeu1 transcript harboring the sequence G36G37. Both our own and their kinetic studies of tRNA transcripts showed that the substrate specificity of A. aeolicus TrmD is rather broad, because A. aeolicus TrmD efficiently methylated all A36 tRNA transcripts as shown in Fig. 3. However, the recognition of A36 is the common feature between E. coli and A. aeolicus TrmD proteins. Furthermore, it has also been reported that E. coli TrmD did not methylate yeast tRNAAsp C36 and U36 mutants (Brulé et al. 2004). E. coli TrmD has this property in common with TrmD of A. aeolicus as demonstrated in this study. Thus, E. coli TrmD recognizes purine36G37 sequence in the same way as A. aeolicus TrmD, suggesting that the purine36G37 recognition by TrmD may be common to the eubacterial domain.
| Discussion |
|---|
|
|
|---|
In the eubacteria domain, tRNA with an A36G37 sequence is very rare (Sprinzl et al. 1998). In fact, the A. aeolicus genome codes 42 tRNA genes except for tRNASec and tmRNA; however, tRNA with an A36G37 sequence is not found (Deckert et al. 1998). Thus, A. aeolicus TrmD does not need to methylate tRNA species with an A36G37 sequence in the living cell. However, it is noteworthy that several tRNA genes with an A36G37 sequence can be found in eubacterial genomes (Sprinzl et al. 1998). Furthermore, it has been reported that tRNAPhe (GAA) from Mycoplasma sp. (Kid) has an A36m1G37 modification (Kimball et al. 1974). It therefore seems that there is widespread distribution of the A36 specificity of TrmD proteins.
A. aeolicus and E. coli TrmD proteins do not recognize a pyrimidine36G37 sequence. In general, it is difficult to distinguish the structures of purine and pyrimidine precisely by hydrogen bonds from amino acid residues. Therefore, the size and/or hydrophobicity of the purine ring at position 36 should be a positive determinant for bacterial TrmD proteins. A. aeolicus may diverge at the earliest period in eubacterial evolution and the amino acid sequence of A. aeolcius TrmD is considerably different from those of other eubacterial TrmD proteins. Therefore, this determinant may be a remnant derived from the common ancestor of eubacterial TrmD proteins. Furthermore, the experiments using the truncated RNA molecules showed that the RNA recognition mechanism of the A. aeolicus TrmD is considerably different from that of E. coli TrmD: the A. aeolicus TrmD methylated the truncated RNA molecules efficiently and the D-stem structure of the tRNA inhibited the activity. In contrast, it has been reported that E. coli TrmD requires the T- and D-arm structures for the efficient methyl acceptance activity (Holmes et al. 1992; Redlak et al. 1997). These differences may be caused by the differences between the protein folds of A. aeolicus and E. coli TrmD proteins, because Liu et al. (2003) have reported the crystal structure of the A. aeolicus TrmD as a new fold, in which the trefoil knot is not presented.
Recently, we have reported the crystal structure of TrmH (Nureki et al. 2004), which is classified into the SPOUT protein group together with TrmD (Anantharaman et al. 2002). The catalytic domain of the SPOUT proteins is classified as class IV fold, which has a deep trefoil knot structure. We proposed a docking model of TrmH with tRNA complex (Nureki et al. 2004). Our model suggests that the interaction of TrmH and tRNA is quite different from that of TrmD and tRNA, even though the protein folds resemble each other. In order to understand the relationship of these two modification enzymes functionally, structurally and evolutionally, further studies will be necessary.
| Experimental procedures |
|---|
|
|
|---|
[Methyl-14C]-AdoMet (1.95 GBq/mmol) and [methyl-3H]-AdoMet (2.47 TBq/mmol) were purchased from ICN. Cold AdoMet was obtained from Sigma, USA. DE52 was a product of Whatman, UK. CM-Toyopearl 650M came from Tosoh, Japan. DNA oligomers were bought from Invitrogen, USA, and T7 RNA polymerase was from Toyobo, Japan. Other chemical reagents were of analytical grade.
Construction of TrmD expression system in E. coli
The trmD gene was amplified by PCR from A. aeolicus genomic DNA using the following primers: Nde1trmDF, 5'-CGG GCA TAT GAG TTC TAA TCC CCT CCG CTT CTT CG-3'; BamH1trmDR, 5'-CGG GGG ATC CTT AAA GAA GGT GCT TGT GC-3'. The underlined sequences show the sites of restriction enzymes, NdeI and BamHI, respectively. The amplified DNA fragment was cloned into the pET30a expression vector (Novagen, USA) utilizing NdeI and BamHI sites. In order to increase the yield of the expressed protein, the DNA sequences between SacI and XhoI sites were replaced by an artificial DNA sequence composed of codons optimized for translation in E. coli. Synthetic DNA oligomers (DC18) were phosphorylated, annealed, and ligated into the vector described above. DNA oligomers used to introduce this alteration were as follows: DC1, 5'-CCT GTC TGG ACA CCA CAA ACT GAT CGA ACT GTG GAA GCT GTG GCA CAG G-3'; DC2, 5'-ATC GAA AAC ACG GTA AAA AAG CGT CCG GAT CTG ATC CCG AAA GAC CTG A-3'; DC3, 5'-CAG AAC TGG AAA AAG ACA TTC TGA ATA GTA TCC TGT CTG GCA AGT CAT T-3'; DC4, 5'-CAA GGA ATG GCT GAA AGA GCA CAA GCA CCT TCT TTA ATA AC-3'; DC5, 5'-TCG AGT TAT TAA AGA AGG TGC TTG TGC TCT TTC AGC CAT TCC TTG AAT GAC TTG CC-3'; DC6, 5'-AGA CAG GAT ACT ATT CAG AAT GTC TTT TTC CAG TTC TGT CAG GTC TTT CGG-3'; DC7, 5'-GAT CAG ATC CGG ACG CTT TTT TAC CGT GTT TTC GAT CCT GTG CCA C-3'; DC8, 5'-AGC TTC CAC AGT TCG ATC AGT TTG TGG TGT CCA GAC AGG AGC T-3'. The resulting plasmid was named pET30trmD and introduced into E. coli BL21(DE3)-Codonplus-RIL strain (Stratagene, USA) for expression.
Expression and purification of TrmD
The expression of TrmD in E. coli was induced by addition of IPTG to a final concentration of 1 mM. After 4 h, the cells were collected by centrifugation at 6500 x g for 20 min, frozen in liquid nitrogen, and stored at 80 °C until required. The cells (5 g) were suspended in 25 mL of the buffer A (50 mM Tris-HCl (pH 7.6), 5 mM MgCl2, 6 mM 2-mercaptoethanol, 50 mM KCl), and disrupted with an ultrasonic disruptor model UD-200 (Tomy, Japan). The cell debris was removed by centrifugation at 8000 g for 20 min. The supernatant fraction was then heated at 70 °C for 30 min. The denatured proteins were removed by centrifugation at 8000 g for 20 min. The supernatant fraction was applied on to a DE52 column (column volume, 10 mL). The flow-through fractions were collected and dialyzed against buffer B (50 mM HEPES-KOH (pH 6.8), 6 mM 2-mercaptoethanol). The dialyzed sample was loaded on to a CM-Toyoperal 650M column (column volume, 10 mL) and TrmD protein eluted with buffer B containing 200 mM KCl. The eluted sample was dialyzed against the buffer A and concentrated with Centriprep YM-10 centrifugal filter devices (Millipore). Glycerol was added to the purified protein to a final concentration of 50% v/v and stored at 30 °C. Quantification of protein was performed with a Bio-Rad protein assay kit using bovine serum albumin as the standard.
Preparation of tRNA transcripts
The yeast tRNAPhe, E. coli tRNALeu (CAG), tRNASer (UGA) and A. aeolcius tRNA transcripts were prepared as described previously (Hori et al. 2002, 2003). Briefly, yeast tRNAPhe, E. coli tRNALeu (CAG) and tRNASer (UGA) genes with T7 RNA polymerase promoter were amplified from plasmids by PCR and then run-off transcriptions were carried out. The template DNAs for A. aeolcius tRNA transcripts were prepared by annealing of two primers and in vitro DNA polymerase reaction: one primer contains the T7 RNA polymerase promoter sequence. For example, in the case of A. aeolicus tRNALeu (CAG), the following primers were used: AALeu (CAG) F, 5'-TAA TAC GAC TCA CTA TAG CGG GTG TGG CGG AAC TGG CAG ACG CGC CGG ACT CAG GAT CCG-3'; AALeu (CAG) R, 5'-GG TGC GGG TGG CGG GAG TCG AAC CCG CAC GCC CTT GCG GGC ACC GGA TCC TGA GTC CGG-3'. The template DNA was purified by phenol-chloroform extraction, recovered by ethanol precipitation and used for the run-off transcription. The primers for yeast tRNAAsp, tRNAPhe (A36U mutant), and H. volcanii tRNA transcripts were newly synthesized and used for in vitro transcription. All transcripts were purified by 10% polyacrylamide gel (7 M urea) electrophoresis.
Measurements of the enzymatic activities
The standard assay used during the enzyme purification measured the transfer of the methyl group from [methyl-14C]-AdoMet to A. aeolicus tRNAPro(GGG) transcript: 300 ng of the protein, 13.5 µM transcript, and 38 µM [methyl-14C]-AdoMet in 30 µL of buffer C (50 mM Tris-HCl (pH 7.6), 5 mM MgCl2, 6 mM 2-mercaptoethanol, 50 mM KCl) were incubated for 5 min at 60 °C. An aliquot (20 µL) was then used for the conventional filter assay. RNA transcripts were prepared as reported previously (Hori et al. 2003). In order to visualize the methyl-transfer reaction, we used 10% polyacrylamide gel electrophoresis and an imaging analyzer system. Briefly, tRNA (0.3 A260 units) was incubated with 300 ng (5 pmol of dimer enzyme) of the purified TrmD and 38 µM [methyl-14C]-AdoMet for 10 min at 60 °C in 30 µL of buffer C, extracted with phenol, collected by ethanol precipitation, and then loaded on to a 10% polyacrylamide gel (7 M urea). The gel was stained with methylene blue and dried. The 14C-methyl incorporation into tRNAs was monitored with a Fuji Photo Film BAS2000 imaging analyzer. Apparent kinetic parameters, Km and Vmax, were determined by a LineweaverBurk plot of the methyl-transfer reaction as measured by the [methyl-3H]-AdoMet by the filter assay: 1 µg of the protein, 38 µM AdoMet, various concentration of the transcript in 50 µL of buffer C (50 mM Tris-HCl (pH 7.6), 5 mM MgCl2, 6 mM 2-mercaptoethanol, 50 mM KCl) were incubated for 5 min at 55 °C. RNA concentrations in the assay were changed to the methyl acceptance activity of the transcript. RNA concentrations in the first measurement were 0.5, 1, 2, 3, 5, 8, and 12 µM. Based on this first measurement, we changed the RNA concentrations. When the estimated Km value was relatively small (< 2 µM), RNA concentrations were changed to be 0.2, 0.3, 0.5, 1, 2, 3, 5, 8, 12, and 15 µM. When the estimated Km value was relatively large (> 2 µM), RNA concentrations were changed to be 0.5, 1, 2, 3, 5, 10, 15, 20 and 30 µM.
Mass spectrometry
Yeast tRNAPhe (50 µg) was methylated with 20 µg enzyme and cold AdoMet for 4 h at 60 °C in 100 µL of buffer C. The RNA was extracted with phenol, recovered by ethanol precipitation, and loaded on to a 10% polyacrylamide gel (7 M urea) electrophoresis. The RNA was visualized by UV (254 nm) irradiation on a thin layer plate (Funacell P-254, Japan), excised, and extracted with 400 µL of gel elution buffer (0.5 M ammonium acetate, 10 mM MgCl2, 1 mM EDTA, and 0.1% SDS). The extracted sample was passed through a Steradisc 13 filter unit (0.2 µM, Kurabo, Japan), and the RNA was recovered by ethanol precipitation. An LCQ ion-trap (IT) mass spectrometer (ThermoFinnigan) equipped with an electrospray ionization (ESI) source and a HP1100 liquid chromatography system (Agilent) was used to analyze the total nucleosides and RNase T1-digested RNA. For nucleoside analysis, the sample (5 µg) was completely digested with nuclease P1 (3 µg) and bacterial alkaline phosphatase (0.1 unit) in 20 µL of reaction mixture containing 50 mM HEPES-KOH (pH 7.5) at 37 °C for 3 h. The hydrolysates were subjected to mass spectrometric analysis as previously described (Hori et al. 2002; Kaneko et al. 2003; Okamoto et al. 2004). For fragment analysis, methylated or unmethylated control tRNA (5 µg each) was digested with RNase T1 (2.5 unit) in 20 µL of reaction mixture containing 25 mM NH4OAc (pH 5.3) at 37 °C for 1 h and subjected to mass spectrometric analysis. Oligo-nucleotides produced by RNase T1 digestion were analyzed in negative ion form by LC/MS as previously described (Hori et al. 2002; Kaneko et al. 2003; Okamoto et al. 2004).
Amino acid sequence alignment and phylogenetic analyses
The amino acid sequences of eubacterial trmD proteins were retrieved from GENBANK. The amino acid sequences were aligned with CLUSTAL X (Thompson et al. 1997) using the default settings. The well aligned regions (163 sites) were selected with GBLOCKS (Castresana 2000).
For Maximum likelihood (ML) analysis (no site-by-site rate variations), the ML pair-wise distances of all combinations were estimated using the PROTML routine in MOLPHY 2.3b (Adachi & Hasegawa 1996), with the D (pair-wise distance calculation) option. A Neighbor-Joining (NJ) tree was then constructed using the NJDIST routine in MOLPHY. The topology of the NJ tree was then used as the initial tree for the ML tree search by the NNI search routine of PROTML.
For reconstruction of the NJ tree (Saitou & Nei 1987), the alpha parameter of gamma distribution of site-by-site variation of evolutionary rate was estimated with Tree-Puzzle 5.0 (Schmidt et al. 2002) with JTT model. The pair-wise distances of all combinations were estimated with PROTDIST in PHYLIP 3.6 (Felsenstein 1997) using the site-by-site variation model. The NJ tree was then reconstructed with NEIGHBOR in PHYLIP 3.6. For bootstrap analysis, the 1000 resamplings were performed with SEQBOOT in PHYLIP 3.6. The PROTDIST was used for estimation of pair-wise distances, and NEIGHBOR was used for reconstruction of NJ trees of resamplings. The bootstrap probability of each node was calculated with CONSENSE of PHYLIP 3.6.
For Maximum parsimony (MP) analysis, the PAUP 4.01b10 was used. The MP tree was estimated by heuristic search under the default condition except use of the option ó the random step-wise addition (10 replicates) of taxon. For bootstrap analysis, 100 resamplings were used.
| Acknowledgements |
|---|
| Footnotes |
|---|
aPresent address: Institute of Environmental Microbiology, Kyowa Kako Co. Ltd, Tadao 2-15-5, Machida 194-0035, Japan
| References |
|---|
|
|
|---|
Ahn, H.J., Kim, H.W., Yoon, H.J., Lee, B., Suh, S.S. & Yang, J.K. (2003) Crystal structure of tRNA (m1G37) methyltransferase: insights into tRNA recognition. EMBO J. 22, 25932603.[CrossRef][Medline]
Anantharaman, V., Koonin, E.V. & Aravind, L. (2002) SPOUT: a class of methyltransferases that includes spoU and trmD RNA methylase superfamilies, and novel superfamilies of predicted prokaryotic RNA methylases. J. Mol. Microbiol. Biotechnol. 4, 7175.[Medline]
Behm-Ansmant, I., Urban, A. & Ma, X., Yu Y.T., Motorin, Y. & Branlant, C. (2003) The Saccharomyces cerevisiae U2 snRNA: pseudouridine-synthase Pus7p is a novel multisite-multisubstrate RNA: Psi-synthase also acting on tRNAs. RNA 9, 13711382.
Björk, G.R. (1992) The role of modified nucleosides in tRNA interactions. In: Transfer RNA in Protein Synthesis (eds D.L. Hatfield, B.J. Lee & R.M. Pirtle), pp. 2385, Boca Raton, FL: CRC Press.
Björk, G.R. (1995) Biosynthesis and function of modified nucleosides. In: tRNA: Structure, Biosynthesis and Function (eds D. Söll & U.L. RajBhandary), pp. 165205. Washington, DC: ASM Press.
Björk, G.R., Wikstrom, P.M. & Byström, A.S. (1989) Prevention of translational frameshifting by the modified nucleoside 1-methylguanosine. Science 244, 986989.
Björk, G.R., Kerstin, J., Nilsson, K., Johansson, M.J.O., Byström, A.S. & Persson, O.P. (2001) A primordial tRNA modification required for the evolution of life? EMBO J. 20, 231239.[CrossRef][Medline]
BruIé H., Elliott, M., Redlak, M., Zehner, Z.E. & Holmes, W.M. (2004) Isolation and characterization of the human tRNA-(N1G37) methyltransferase (TRM5) and comparison to the Escherichia coli TrmD protein. Biochemistry 43, 92439255.[CrossRef][Medline]
Burggraf, S., Olsen, G.J., Stetter, K.O. & Woese, C.R. (1992) A phylogenetic analysis of Aquifex pyrophilus. Syst. Appl. Microbiol. 15, 353356.
Byström, A.S. & Björk, G.R. (1982) Chromosomal location and cloning of the gene (trmD) responsible for the synthesis of tRNA (m1G) methyltransferase in Escherichia coli K-12. Mol. Gen. Genet. 188, 440446.[CrossRef][Medline]
Byström, A.S., Hjaimarsson, K.J., Wikstrom, P.M. & Björk, G.R. (1983) The nucleotide sequence of an Escherichia coli operon containing genes for the tRNA (m1G) methyltransferase, the ribosomal proteins S16 and L19 and a 21-K polypeptide. EMBO J. 2, 899905.[Medline]
Castresana, J. (2000) Election of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 4, 540552.
CavailIé, J., Chetouani, F. & Bachellerie, J.-P. (1999) The yeast Saccharomyces cerevisiae YDL112w ORF encodes the putative 2'-O-ribose methyltransferase catalyzing the formation of Gm18 in tRNAs. RNA 5, 6681.[Abstract]
Christian, T., Evilia, C., Williams, S. & Hou, Y.-M. (2004) Distinct origins of tRNA (m1G37) methyltransferase. J. Mol. Biol. 339, 707719.[CrossRef][Medline]
Deckert, G., Warren, P.V., Gaasterland, T., Young, W.G., Lenox, A.L., Graham, D.E., Overbeek, R., Snead, M.A., Keller, M., Aujay, M., Huber, R., Feldman, R.A., Short, J.M., Olsen, G.J. & Swanson, R.V. (1998) The complete genome of the hyperthermophilic bacterium Aquifex aeolicus. Nature 392, 353358.[CrossRef][Medline]
Elkins, P.A., Watts, J.M., Zalacain, M., van Thiel, A., Vitazka, P.R., Redlak, M., Andraos-Selim, C., Rastinejad, F. & Holmes, W.M. (2003) Insights into catalysis by a knotted TrmD tRNA methyltransferase. J. Mol. Biol. 333, 931949.[CrossRef][Medline]
Farabaugh, P.J. & Björk, G.R. (1999) How translational accuracy influences reading frame maintenance. EMBO J. 18, 14271434.[CrossRef][Medline]
Felsenstein, J. (1997) An alternating least squares approach to inferring phylogenies from pairwise distances. Syst. Biol. 46, 101111.[CrossRef][Medline]
Garcia, G.R. & Goodenough-Lashhua, D.M. (1998) Appendix 3: general properties of RNA-modifing and ñediting enzymes. In: Modification and Editing of RNA (eds H. Grosjean & R. Benne), pp. 555560. Washington, DC: ASM Press.
Griffiths, E. & Gupta, R.S. (2001) The use of signature sequences in different proteins to determine the relative branching order of bacterial divisions: evidence that Fibrobacter diverged at a similar time to Chlamydia and the Cytophaga-Flavobacterium-Bacteroides division. Microbiology 147, 26112622.
Griffiths, E. & Gupta, R.S. (2006) Molecular signatures in protein sequences that are characteristics of the phylum Aquificae. Int. J. Syst. Evol. Microbiol. 56, 99107.[CrossRef][Medline]
Grosjean, H., Droogmans, L., Giegé, R. & Uhlenbeck, O.C. (1990) Guanosine modifications in runoff transcripts of synthetic transfer RNA-Phe genes microinjected into Xenopus oocytes. Biochim. Biophys. Acta 1050, 267273.[Medline]
Hagervall, T.G., Ericson, J.U., Esberg, K.B., Li, J.N. & Björk, G.R. (1990) Role of tRNA modification in translational fidelity. Biochim. Biophys. Acta 1050, 263266.[Medline]
Hagervall, T.G., Tuohy, T.M.F., Atkins, J.F. & Björk, G.R. (1993) Deficiency of 1-methylguanosine in tRNA from Salmonella typhimurium induces frameshifting by quadruplet translocation. J. Mol. Biol. 232, 756765.[CrossRef][Medline]
Hjalmarsson, K.J., Byström, A.S. & Björk, G.R. (1983) Purification and characterization of transfer RNA (guanine-1) methyltransferase from Escherichia coli. J. Biol. Chem. 258, 13431351.
Holmes, W.M., Andraos-Selim, C., Roberts, I. & Wahab, S.Z. (1992) Structural requirements for tRNA methylation. Action of Escherichia coli tRNA (guanosine-1) methyltransferase on tRNALeu1 structural variants. J. Biol. Chem. 267, 1344013445.
Holmes, W.M., Andraos-Selim, C. & Redlak, M. (1995) tRNAm1G methyltransferase interactions: touching bases with structure. Biochimie 77, 6265.[Medline]
Hori, H., Saneyoshi, M., Kumagai, I., Miura, K. & Watanabe, K. (1989) Effects of modification of 4-thiouridine in E. coli tRNA (fMet) on its methyl acceptor activity by thermostable Gm-methylases. J. Biochem. 106, 798802.
Hori, H., Yamazaki, N., Matsumoto, T., Watanabe, Y., Ueda, T., Nishikawa, K., Kumagai, I. & Watanabe, K. (1998) Substrate recognition of tRNA (Guanosine-2'-)-methyltransferase from Thermus thermophilus HB27. J. Biol. Chem. 273, 2572125727.
Hori, H., Suzuki, T., Sugawara, K., Inoue, Y., Shibata, T., Kuramitsu, S., Yokoyama, S., Oshima, T. & Watanabe, K. (2002) Identification and characterization of tRNA (Gm18) methyltransferase from Thermus thermophilus HB8: domain structure and conserved amino acid sequence motifs. Genes Cells 7, 259272.[Abstract]
Hori, H., Kubota, S., Watanabe, K., Kim, J.M., Ogasawara, T., Sawasaki, T. & Endo, Y. (2003) Aquifex aeolicus tRNA (Gm18) methyltransferase has unique substrate specificity. tRNA recognition mechanism of the enzyme. J. Biol. Chem. 278, 2508125090.
Kaneko, T., Suzuki, T., Kapushoc, S.T., Rubio, M.A., Ghazvini, J., Watanabe, K., Simpson, L. & Suzuki, T. (2003) Wobble modification differences and subcellular localization of tRNAs in Leishmania tarentolae: implication for tRNA sorting mechanism. EMBO J. 22, 657667.[CrossRef][Medline]
Keith, G. (1995) Mobilities of modified ribonucleotides on two-dimensional cellulose thin-layer chromatography. Biochimie 77, 142144.[Medline]
Kimball, M.E., Szeto, K.S. & Söll, D. (1974) The nucleotide sequence of phenylalanine tRNA from Mycoplasma sp. (Kid). Nucleic Acids Res. 1, 17211732.
Lecointe, F., Simos, G., Sauer, A., Hurt, E.C., Motorin, Y. & Grosjean, H. (1998) Characterization of yeast protein Deg1 as pseudouridine synthase (Pus3) catalyzing the formation of psi 38 and psi 39 in tRNA anticodon loop. J. Biol. Chem. 273, 13161323.
Li, J. & Björk, G.R. (1999) Structural alterations of the tRNA (m1G37) methyltransferase from Salmonella typhimurium affect tRNA substrate specificity. RNA 5, 395408.[Abstract]
Liu, J., Wang, W., Shin, D.H., Yokota, H., Kim, R. & Kim, S.-H. (2003) Crystal structure of tRNA (m1G37) methyltransferase from Aquifex aeolicus at 2.6 Å resolution: a novel methyltransferase fold. Proteins 53, 326328.[CrossRef][Medline]
Matsumoto, T., Nishikawa, K., Hori, H., Ohta, T., Miura, K. & Watanabe, K. (1990) Recognition sites of tRNA by a thermostable tRNA (guanosine-2'-)-methyltransferase from Thermus thermophilus HB27. J. Biochem. 107, 331338.
McCloskey, J.A. & Crain, P.F. (1998) The RNA modification databaseó1998. Nucleic Acids Res. 26, 196197.
Motorin, Y. & Grosjean, H. (1998) Appendix 1: chemical structure and classification of posttranscriptionally modified nucleosides in RNA. In: Modification and Editing of RNA (eds H. Grosjean & R. Benne), pp. 543550. Washington, DC: ASM Press.
Motorin, Y. & Grosjean, H. (1999) Multisite-specific tRNA: m5C-methyltransferase (Trm4) in yeast Saccharomyces cerevisiae: identification of the gene and substrate specificity of the enzyme. RNA 5, 11051118.[Abstract]
Motorin, Y., Keith, G., Simon, C., Foiret, D., Simos, G., Hurt, E. & Grosjean, H. (1998) The yeast tRNA: pseudouridine synthase Pus1p displays a multisite substrate specificity. RNA 4, 856859.[Abstract]
Nureki, O., Shirouzu, M., Hashimoto, K., Ishitani, R., Terada, T., Tamakoshi, M., Oshima, T., Chijimatsu, M., Takio, K., Vassylyev, D.G., Shibata, T., Inoue, Y., Kuramitsu, S. & Yokoyama, S. (2002) An enzyme with a deep trefoil knot for the active-site architecture. Acta Crystallogr. D Biol. Crystallogr. 58, 11291137.[CrossRef][Medline]
Nureki, O., Watanabe, K., Fukai, S., Ishii, R., Endo, Y., Hori, H. & Yokoyama, S. (2004) Deep knot structure for construction of active site and cofactor binding site of tRNA modification enzyme. Structure 12, 593602.[Medline]
Okamoto, H., Watanabe, K., Ikeuchi, Y., Suzuki, T., Endo, Y. & Hori, H. (2004) Substrate tRNA recognition mechanism of tRNA (m7G46) methyltransferase from Aquifex aeolicus. J. Biol. Chem. 279, 4915149159.
Persson, B.C., Jager, G. & Gustafsson, C. (1997) The spoU gene of Escherichia coli, the fourth gene of the spoT operon, is essential for tRNA (Gm18) 2'-O-methyltransferase activity. Nucleic Acids Res. 25, 40934097.
Pintard, L., Lecointe, F., Bujnicki, J.M., Bonnerot, C., Grosjean, H. & Lapeyre, B. (2002) Trm7p catalyses the formation of two 2'-O-methylriboses in yeast tRNA anticodon loop. EMBO J. 21, 18111820.[CrossRef][Medline]
Redlak, M., Andraos-Selim, C., Giegé. R., Florentz, C. & Holmes, W.M. (1997) Interaction of tRNA with tRNA (guanosine-1) methyltransferase: binding specificity determinants involve the dinucleotide G36pG37 and tertiary structure. Biochemistry 36, 86998709.[CrossRef][Medline]
Rozenski, J., McCloskey, J.A. & Crain, P.F. (1999) The RNA Modification Database: 1999 update. Nucleic Acids Res. 27, 196197.
Saitou, N. & Nei, M. (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406425.[Abstract]
Schmidt, H.A., Strimmer, K., Vingron, M. & von Haeseler, A. (2002) TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics 3, 502504.
Sprinzl, M., Horn, C., Brown, M., Ioudovitch, A. & Steinberg, S. (1998) Compilation of tRNA sequences and sequences of tRNA genes. Nucleic Acids Res. 26, 148153.
Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F. & Higgins, D.G. (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25, 48764882.
Watanabe, K., Nureki, O., Fukai, S., Ishii, R., Okamoto, H., Yokoyama, S., Endo, Y. & Hori, H. (2005) Roles of conserved amino acid sequence motifs in the SpoU (TrmH) RNA methyltransferase family. J. Biol. Chem. 280, 1036810377.
Woese, C.R., Kandler, O. & Wheelis, M.L. (1990) Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc. Natl. Acad. Sci. USA 87, 45764579.
Received: 25 March 2006
Accepted: 27 August 2006
This article has been cited by other articles:
![]() |
T. Toyooka, T. Awai, T. Kanai, T. Imanaka, and H. Hori Stabilization of tRNA (m1G37) methyltransferase [TrmD] from Aquifex aeolicus by an intersubunit disulfide bond formation Genes Cells, August 1, 2008; 13(8): 807 - 816. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lee, G. Kramer, D. E. Graham, and D. R. Appling Yeast Mitochondrial Initiator tRNA Is Methylated at Guanosine 37 by the Trm5-encoded tRNA (Guanine-N1-)-methyltransferase J. Biol. Chem., September 21, 2007; 282(38): 27744 - 27753. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP |