Search tips
Search criteria 


Logo of ccLink to Publisher's site
Cell Cycle. 2011 October 15; 10(20): 3566–3570.
Published online 2011 October 15. doi:  10.4161/cc.10.20.17857
PMCID: PMC3266182

Identification and modeling of a phosphatase-like domain in a tRNA 2′-O-ribosyl phosphate transferase Rit1p


Cytoplasmic initiator tRNAs from plants and fungi are excluded from participating in translational elongation by the presence of a unique 2′-phosphoribosyl modification of purine 64, introduced posttranscriptionally by the enzyme Rit1p. Members of the Rit1p family show no obvious similarity to other proteins or domains, there is no structural information available to guide experimental analyses, and the mechanism of action of this enzyme remains a mystery. Using protein fold recognition, we identified a phosphatase-like domain in the C-terminal part of Rit1p. A comparative model of the C-terminal domain was constructed and used to predict the function of conserved residues and to propose the mechanism of action of Rit1p. The model will facilitate experimental analyses of Rit1p and its interactions with the initiator tRNA substrate.

Key words: fold recognition, homology modeling, tRNA modification, Rit1p, bioinformatics


tRNA's role in decoding the genome is critical to the accuracy and efficiency of protein synthesis. Thus, the understanding of tRNA gene expression, folding of the transcripts, posttranscriptional processing, subcellular localization and their interactions with other macromolecules is essential for the full understanding of protein translation.14 A key step in the initiation of protein synthesis involves the recognition of the correct starting codon on the mRNA. The protein-synthesizing apparatus of the cell must discriminate between different AUG codons encoding the same amino acid methionine—the start codon and those in the internal coding region of the mRNA. A specialized initiator tRNA (tRNAiMet in Eukaryota and tRNAfmet in Prokaryota and organelles) exclusively recognizes the start codon located in the ribosomal P site during the initiation of translation, whereas the elongator tRNA (tRNAmMet) is used to recognize internal methionine codons in the ribosomal A site.5 In bacteria, mitochondria and chloroplasts, the Met of Met-tRNAfmet is formylated, which blocks the α-amino group from forming a peptide bond. On the other hand, the eukaryotic cytoplasmic Met-tRNAiMet is not formylated, but possesses a number of unique sequence features not found or rare in elongator tRNAs. These features include the absence of residue 17, presence of an additional residue A20 and the replacement of T54 and pyrimidine 60 in the TΨC loop (found in virtually all elongator tRNAs) by A54 and A60. As revealed by the crystal structure of tRNAiMet from Saccharomyces cerevisiae,6 these features cluster together in space to form a unique substructure stabilized by an extensive and intricate network of hydrogen bonds.

Cytoplasmic tRNAsiMet from fungi and plants (but not from vertebrates) carry a posttranscriptional modification that links a 5′-phosphoribosyl group by a glycosidic bond to the O2′ of ribose of purine 64 in the T stem.7 It was shown that the demodification of A64 in yeast tRNAiMet or G64 in wheat germ tRNAsiMet allows these tRNAs to read not only the start AUG codons, but also the internal AUG codons.8 Thus, this posttranscriptional modification excludes the initiator tRNA from participating in translational elongation.

Thus far, over 100 posttrascriptionally modified nucleosides have been reported in different types of RNAs from diverse organisms.9 The largest number of modified nucleosides with the greatest structural diversity is found in tRNAs.10 These modifications are introduced by a large number of enzymes, which may act alone or in complexes and pathways.1114 While many of these enzymes have been extensively characterized, very little is known about others (for details, see the collection of reviews in books, refs. 1517). For example, the RIT1 gene encoding the 2′-O-ribosyl phosphate transferase (Rit1p) was identified in yeast more than 10 y ago.18 However, the mechanism of action of this enzyme remains a mystery. In particular, no structural information is available for Rit1p to guide experimental analyses and determine the molecular basis of interactions with tRNAiMet.


We submitted the sequence of S. cerevisiae Rit1p (NCBI gene identification number 6323939) to the GeneSilico metaserver19 to predict the protein structure by fold-recognition methods. This analysis revealed that Rit1p is composed of two domains (roughly residues 1–340 and 341–512). The N-terminal domain (NTD) exhibited no evident similarity to any known protein structure or family. However, the C-terminal domain (CTD) of Rit1p was found to exhibit significant similarity to structures of dual-specificity phosphatases (DUSPs). Nearly all individual fold-recognition methods reported DUSP structures as the most preferred matches (e.g., HHSEARCH: 2 hcm_A, score 97.88, and 1wrm_A, score 97.28, with 100 being ideal match, and mGen-Threader: 1wrm_A, score 0.001, with 0 being ideal match). Based on these primary predictions, the consensus predictor PCONS unequivocally identified the DUSP fold as the only reasonable template for modeling: the top 10 matches exhibited scores in the range of 1.94–2.65, which indicate a highly confident prediction (see Methods for explanation of scores). The similarity between Rit1p and DUSPs was originally described in a B.Sc. thesis20 of the first author of this article, and it has been independently identified using HHsearch21 by Arcady Mushegian.22

DUSPs are a heterogeneous group of protein phosphatases that can dephosphorylate both phosphotyrosine and phosphoserine or phosphothreonine residues, also within the same substrate. They have been implicated as major modulators of signaling pathways critical for cell growth and differentiation (reviewed in refs. 2325). They are often composed of multiple domains, with the conserved catalytic domains fused to other domains with various folds and functions.23 There are also known cases of DUSP-related proteins, such as MKSTYX, that lack the conserved residues in the active site and appear enzymatically inactive. While the evolutionary relationship and the resulting structural similarity between the CTD of Rit1p and DUSPs appears evident beyond any reasonable doubt, the question arises whether Rit1p possesses the DUSP-like active site and may exhibit a similar activity.

Based on the alignments from the metaserver, we have built a comparative model of the CTD (Fig. 1) using the “FRankenstein's monster” modeling approach.26 Details of the modeling procedure and an explanation of the model assessment scores are presented in the Materials and Methods section. The best-scoring structure was evaluated as potentially “extremely good” by the PROQ method (predicted LGscore: 4.317), and the MetaMQAP method27 predicted that the model's root mean square deviation to the (currently unknown) true structure is around 3.4 Å. These values indicate that our structural prediction is likely to be correct and support the prediction that Rit1 has a common evolutionary ancestor with phosphatases. The model can, therefore, be used to make functional inferences at the level of individual residues (however, not with atomic precision) in connection with the analysis of sequence conservation in the Rit1p family.

Figure 1
A model of S. cerevisiae Rit1 CTD. Coordinates are available from The protein backbone is shown as a ribbon. The substrate phosphate group is shown in light gray and residues predicted to be important for catalysis ...

A search of the non-redundant database at the NCBI revealed that Rit1 homologs are present in green plants, fungi and Giardia lamblia, and that all members of this family possess both domains. Analysis of the multiple sequence alignment (Fig. 2) reveals that the NTD is more conserved and contains several sequence motifs rich in invariant or nearly invariant residues, while the CTD contains two conserved motifs. One of the conserved motifs (hhhx-CxxGxDhS, where h indicates a hydrophobic residue and x a non-conserved residue) in the CTD corresponds to the active site motif of dual-specificity phosphatases. Interestingly, most (but not all) Rit1 family members possess a conserved Cys residue (C445 in S. cerevisiae Rit1p) that corresponds to the catalytic Cys residue of DUSPs, which is used to form a cysteinyl-phosphate intermediate in the phosphate group transfer reaction.28 A catalytic Asp residue that fulfills a role of general acid in the reaction performed by DUSPs has no homologous counterpart in the Rit1 family. However, a spatially equivalent position in the putative phosphatebinding pocket is occupied by a conserved Asp residue (D450 in S. cerevisiae Rit1p). This spatial equivalence, which became apparent only upon inspection of the model and was not obvious from the sequence alignment, suggests that the CTD of Rit1 proteins possesses a putative active site similar to those of DUSPs.

Figure 2
Multiple alignment of representative sequences from the Rit1 family and the known dual specificity phosphatase structure (2 hcm) identified by the fold-recognition analysis as the best template for modeling of the CTD. Secondary structures predicted for ...


Our prediction of a DUSP domain in Rit1p with a partially conserved DUSP-like active site suggests that this domain may be directly involved in the phosphotransferase activity of Rit1p, based on the same mechanism used by DUSPs. However, the structure of the remaining N-terminal part of Rit1p remains unknown. The structurally uncharacterized NTD is actually more conserved than the DUSP-like CTD. Therefore, it cannot be excluded that the transfer of the phosphoribosyl group is catalyzed by the conserved NTD, while the CTD performs some secondary function. There are known cases of inactivation of the enzymatic function of a protein and utilization of the enzyme-like domain for non-enzymatic roles. Among RNA modification enzymes, relevant examples are the Gcd10p subunit of tRNA:m1A58 MTase Trm6p,29,30 Sen15 and Sen54 subunits of the tRNA splicing endonuclease31 or the N-terminal domain of rRNA MTase RsmC.32 In these proteins, “degenerated” enzyme-like domains are typically used for RNA binding. The partial conservation of the DUSP-like active site suggests that the CTD of Rit1p may be involved in interactions with the phosphate group, but this interaction could be limited to binding of the phosphate backbone of the tRNA substrate. Alternatively, the NTD and the CTD may fulfill distinct enzymatic functions; a number of RNA modification enzymes that catalyze more than one reaction have been characterized, exemplified by the trifunctional capping enzyme in the mimivirus, which is composed of a metal-dependent RNA triphosphatase module fused to a guanylyltransferase and methyltransferase domains.33 Such can be the case with Rit1p one of its domains could be involved in an activity that yet remains to be discovered.

The ultimate test for our predictions may be provided by experimental analysis, which is out of scope of this article. The multiple sequence alignment, which was used as the basis for our modeling of the CTD, suggests a number of conserved residues in both domains that may be targeted by site-directed mutagenesis for functional studies of Rit1p. The charged or hydrophilic character of these conserved residues suggests involvement in recognition and binding of the tRNA substrate or in other aspects of catalysis. Potential targets for further analyses include invariant residues in the NTD (e.g., S48, D52, N72, R74, H97, R107, N109, D136, T141, D280, D281 and E283) and highly conserved C445, D450, S452 and R502 in the CTD. These experiments may provide clues as to whether Rit1p CTD is enzymatically active or plays only a supporting role by, e.g., binding the substrate RNA.

Materials and Methods

Sequence database searches and alignment.

Searches of the non-redundant version of current sequence database (nr) were performed using PSI-BLAST,34 with the e-value threshold of 1e-3. The multiple sequence alignment of Rit1p and proteins identified in the database was calculated using MUSCLE35 with default parameters and refined by hand to ensure that no unwarranted gaps had been introduced within α-helices and β-strands.

Protein fold prediction.

Protein structure predictions were performed via the GeneSilico metaserver gateway19 for secondary structure prediction and fold recognition. Structural predictions were performed both for the full-length Rit1p sequence and for the individual domains. It is important to indicate that different FR servers use completely different scoring systems with different scales (e.g., Z-scores, e-values, percent values, etc.). Moreover, the meaning of scores changes over time and may not be the same as reported in original publications describing the methods, as servers are modified and databases grow continuously. The comparable reliability thresholds for a number of servers are estimated, e.g., by the Livebench benchmark36,37 conducted by Leszek Rychlewski and coworkers. The GeneSilico metaserver uses the PCONS38 consensus method to compare different predictions with each other and to assign its overall score. According to Livebench-8 analysis, a PCONS score of 0.67 corresponds to the average score for their 8th incorrect predictions, which amount to 5% of incorrect predictions for all targets in the Livebench test set.39 In general, scores > 1 indicate confident predictions with low chance of false positives.

Protein structure modeling.

Homology modeling of the catalytic core was performed using the “FRankenstein's monster” approach (see ref. 26 and 40 for a detailed description). Briefly, preliminary models were built with MODELLER41 based on the sequence alignment between Rit1p and the template structure 2 hcm obtained from various fold-recognition servers queried via the GeneSilico metaserver. The preliminary models were scored by MetaMQAP to predict their accuracy at the level of individual residues,27 and a hybrid model was generated by merging fragments with consensus alignment with those non-consensus fragments that exhibited the best MetaMQAP scores. Additional evaluations of protein structure quality were performed with PROQ.42 The hybrid model was refined with the REFINER method, which uses a reduced representation of the protein chain and a statistical potential of mean force to describe intramolecular interactions.43 REFINER is a real-space version of a lattice-based algorithm CABS44 that we had earlier successfully combined with the “FRankenstein's monster” method in CASP6.45 The folding was performed with restraints on predicted secondary structure. Models generated during the simulation had their full-atom representation rebuilt and were scored using PROQ and MetaMQAP.


This analysis was supported by Polish Ministry of Science and Higher Education (MNiSW, grants number N301 010 31/0219 and Iuventus Plus 0067/P01/2010/70 to A.C.). J.M.B. was supported by the European Research Council (ERC, StG grant RNA+p = 123D) and by the “Ideas for Poland” fellowship from the Foundation for Polish Science.


fold recognition
expectation value
C-terminal domain
protein tyrosine phosphatase

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.


1. McFarlane RJ, Whitehall SK. tRNA genes in eukaryotic genome organization and reorganization. Cell Cycle. 2009;8:3102–3106. doi: 10.4161/cc.8.19.9625. [PubMed] [Cross Ref]
2. Huynh LN, Thangavel M, Chen T, Cottrell R, Mitchell JM, PraetoriusIbba M. Linking tRNA localization with activation of nutritional stress responses. Cell Cycle. 2010;9:3112–3118. doi: 10.4161/cc.9.15.12525. [PubMed] [Cross Ref]
3. Kantidakis T, White RJ. A feedback loop between mTOR and tRNA expression? Cell Cycle. 2010;9:2934–2935. doi: 10.4161/cc.9.15.12737. [PubMed] [Cross Ref]
4. Lopez Sanchez MI, Mercer TR, Davies SM, Shearwood AM, Nygard KK, Richman TR, et al. RNA processing in human mitochondria. Cell Cycle. 2011;10 [PubMed]
5. Kozak M. Comparison of initiation of protein synthesis in procaryotes, eucaryotes and organelles. Microbiol Rev. 1983;47:1–45. [PMC free article] [PubMed]
6. Basavappa R, Sigler PB. The 3 A crystal structure of yeast initiator tRNA: functional implications in initiator/elongator discrimination. EMBO J. 1991;10:3105–3111. [PubMed]
7. Desgrès J, Keith G, Kuo KC, Gehrke CW. Presence of phosphorylated Oribosyladenosine in Tpsistem of yeast methionine initiator tRNA. Nucleic Acids Res. 1989;17:865–882. doi: 10.1093/nar/17.3.865. [PMC free article] [PubMed] [Cross Ref]
8. Kiesewetter S, Ott G, Sprinzl M. The role of modified purine 64 in initiator/elongator discrimination of tRNA(iMet) from yeast and wheat germ. Nucleic Acids Res. 1990;18:4677–4682. doi: 10.1093/nar/18.16.4677. [PMC free article] [PubMed] [Cross Ref]
9. Cantara WA, Crain PF, Rozenski J, McCloskey JA, Harris KA, Zhang X, et al. The RNA Modification Database, RNAMDB: 2011 update. Nucleic Acids Res. 2011;39:195–201. doi: 10.1093/nar/gkq1028. [PMC free article] [PubMed] [Cross Ref]
10. Jühling F, Morl M, Hartmann RK, Sprinzl M, Stadler PF, Putz J. tRNAdb 2009: compilation of tRNA sequences and tRNA genes. Nucleic Acids Res. 2009;37:159–162. doi: 10.1093/nar/gkn772. [PMC free article] [PubMed] [Cross Ref]
11. Dunin-Horkawicz S, Czerwoniec A, Gajda MJ, Feder M, Grosjean H, Bujnicki JM. MODOMICS: a database of RNA modification pathways. Nucleic Acids Res. 2006;34:145–149. doi: 10.1093/nar/gkj084. [PMC free article] [PubMed] [Cross Ref]
12. Czerwoniec A, Dunin-Horkawicz S, Purta E, Kaminska KH, Kasprzak JM, Bujnicki JM, et al. MODOMICS: a database of RNA modification pathways. 2008 update. Nucleic Acids Res. 2009;37:118–121. doi: 10.1093/nar/gkn710. [PMC free article] [PubMed] [Cross Ref]
13. Noma A, Kirino Y, Ikeuchi Y, Suzuki T. Biosynthesis of wybutosine, a hypermodified nucleoside in eukaryotic phenylalanine tRNA. EMBO J. 2006;25:2142–2154. doi: 10.1038/sj.emboj.7601105. [PubMed] [Cross Ref]
14. Rossmanith W, Holzmann J. Processing mitochondrial (t)RNAs: new enzyme, old job. Cell Cycle. 2009;8:1650–1653. doi: 10.4161/cc.8.11.8502. [PubMed] [Cross Ref]
15. Grosjean H. Fine-tuning of RNA functions by modification and editing. Berlin-Heidelberg: Springer-Verlag; 2005.
16. Grosjean H. DNA and RNA Modification Enzymes: Structure, Mechanism, Function and Evolution. 2009:682.
17. Smith HC, editor. RNA and DNA editing. Hoboken, New Jersey: John Wiley & Sons, Inc; 2008.
18. Åström SU, Bystrom AS. Rit1, a tRNA backbone-modifying enzyme that mediates initiator and elongator tRNA discrimination. Cell. 1994;79:535–546. doi: 10.1016/0092-8674(94)90262-3. [PubMed] [Cross Ref]
19. Kurowski MA, Bujnicki JM. GeneSilico protein structure prediction meta-server. Nucleic Acids Res. 2003;31:3305–3307. doi: 10.1093/nar/gkg557. [PMC free article] [PubMed] [Cross Ref]
20. Czerwoniec A. Structure prediction and phylogenetic analysis of the Rit1p phosphoribosyltransferase in Polish. Faculty of Biology. Poznan: Adam Mickiewicz University; 2006. p. 37.
21. Söding J. Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005;21:951–960. doi: 10.1093/bioinformatics/bti125. [PubMed] [Cross Ref]
22. Mushegian AR, editor. Foundations of Comparative Genomics. Academic Press; 2007. Finding Sequence Similarities; pp. 5–23.
23. Patterson KI, Brummer T, O'Brien PM, Daly RJ. Dual-specificity phosphatases: critical regulators with diverse cellular targets. Biochem J. 2009;418:475–489. [PubMed]
24. Cerignoli F, Rahmouni S, Ronai Z, Mustelin T. Regulation of MAP kinases by the VHR dual-specific phosphatase: implications for cell growth and differentiation. Cell Cycle. 2006;5:2210–2215. doi: 10.4161/cc.5.19.3267. [PubMed] [Cross Ref]
25. St Clair S, Manfredi JJ. The dual specificity phosphatase Cdc25C is a direct target for transcriptional repression by the tumor suppressor p53. Cell Cycle. 2006;5:709–713. doi: 10.4161/cc.5.7.2628. [PubMed] [Cross Ref]
26. Kosinski J, Cymerman IA, Feder M, Kurowski MA, Sasin JM, Bujnicki JMA. “FRankenstein's monster” approach to comparative modeling: merging the finest fragments of Fold-Recognition models and iterative model refinement aided by 3D structure evaluation. Proteins. 2003;53:369–379. doi: 10.1002/prot.10545. [PubMed] [Cross Ref]
27. Pawlowski M, Gajda MJ, Matlak R, Bujnicki JM. MetaMQAP: a meta-server for the quality assessment of protein models. BMC Bioinformatics. 2008;9:403. doi: 10.1186/1471-2105-9-403. [PMC free article] [PubMed] [Cross Ref]
28. Zhang ZY, Dixon JE. Protein tyrosine phosphatases: mechanism of catalysis and substrate specificity. Adv Enzymol Relat Areas Mol Biol. 1994;68:1–36. [PubMed]
29. Bujnicki JM. In silico analysis of the tRNA:m1A58 methyltransferase family: homology-based fold prediction and identification of new members from Eubacteria and Archaea. FEBS Lett. 2001;507:123–127. doi: 10.1016/S0014-5793(01)02962-3. [PubMed] [Cross Ref]
30. Ozanick SG, Bujnicki JM, Sem DS, Anderson JT. Conserved amino acids in each subunit of the heteroligomeric tRNA m1A58 Mtase from Saccharomyces cerevisiae contribute to tRNA binding. Nucleic Acids Res. 2007 [PMC free article] [PubMed]
31. Bujnicki JM, Rychlewski L. Prediction of a common fold for all four subunits of the yeast tRNA splicing endonuclease: implications for the evolution of the EndA/Sen family. FEBS Lett. 2000;486:328–329. doi: 10.1016/S0014-5793(00)02322-X. [PubMed] [Cross Ref]
32. Sunita S, Purta E, Durawa M, Tkaczuk KL, Swaathi J, Bujnicki JM, et al. Functional specialization of domains tandemly duplicated within 16S rRNA methyltransferase RsmC. Nucleic Acids Res. 2007;35:4264–4274. doi: 10.1093/nar/gkm411. [PMC free article] [PubMed] [Cross Ref]
33. Benarroch D, Smith P, Shuman S. Characterization of a trifunctional mimivirus mRNA capping enzyme and crystal structure of the RNA triphosphatase domain. Structure. 2008;16:501–512. doi: 10.1016/j.str.2008.01.009. [PubMed] [Cross Ref]
34. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–3402. doi: 10.1093/nar/25.17.3389. [PMC free article] [PubMed] [Cross Ref]
35. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–1797. doi: 10.1093/nar/gkh340. [PMC free article] [PubMed] [Cross Ref]
36. Bujnicki JM, Elofsson A, Fischer D, Rychlewski L. LiveBench-1: continuous benchmarking of protein structure prediction servers. Protein Sci. 2001;10:352–361. doi: 10.1110/ps.40501. [PubMed] [Cross Ref]
37. Bujnicki JM, Elofsson A, Fischer D, Rychlewski L. LiveBench-2: large-scale automated evaluation of protein structure prediction servers. Proteins. 2001:184–191. doi: 10.1002/prot.10039. [PubMed] [Cross Ref]
38. Lundström J, Rychlewski L, Bujnicki J, Elofsson A. Pcons: a neural-network-based consensus predictor that improves fold recognition. Protein Sci. 2001;10:2354–2362. doi: 10.1110/ps.08501. [PubMed] [Cross Ref]
39. Rychlewski L, Fischer D. LiveBench-8: the large-scale, continuous assessment of automated protein structure prediction. Protein Sci. 2005;14:240–245. doi: 10.1110/ps.04888805. [PubMed] [Cross Ref]
40. Kosinski J, Gajda MJ, Cymerman IA, Kurowski MA, Pawlowski M, Boniecki M, et al. FRankenstein becomes a cyborg: the automatic recombination and realignment of fold recognition models in CASP6. Proteins. 2005;61:106–113. doi: 10.1002/prot.20726. [PubMed] [Cross Ref]
41. Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234:779–815. doi: 10.1006/jmbi.1993.1626. [PubMed] [Cross Ref]
42. Wallner B, Elofsson A. Can correct protein models be identified? Protein Sci. 2003;12:1073–1086. doi: 10.1110/ps.0236803. [PubMed] [Cross Ref]
43. Boniecki M, Rotkiewicz P, Skolnick J, Kolinski A. Protein fragment reconstruction using various modeling techniques. J Comput Aided Mol Des. 2003;17:725–738. doi: 10.1023/B:JCAM.0000017486.83645.a0. [PubMed] [Cross Ref]
44. Kolinski A. Protein modeling and structure prediction with a reduced representation. Acta Biochim Pol. 2004;51:349–371. [PubMed]
45. Kolinski A, Bujnicki JM. Generalized protein structure prediction based on combination of fold-recognition with de novo folding and evaluation of models. Proteins. 2005;61:84–90. doi: 10.1002/prot.20723. [PubMed] [Cross Ref]

Articles from Cell Cycle are provided here courtesy of Landes Bioscience