MODOMICS: a database of RNA modifications and related information. 2023 update.

The MODOMICS database was updated with recent data and now includes new data types related to RNA modifications. Changes to the database include an expanded modification catalog, encompassing both natural and synthetic residues identified in RNA structures. This addition aids in representing RNA sequences from the RCSB PDB database more effectively. To manage the increased number of modifications, adjustments to the nomenclature system were made. Updates in the RNA sequences section include the addition of new sequences and the reintroduction of sequence alignments for tRNAs and rRNAs. The protein section was updated and connected to structures from the RCSB PDB database and predictions by AlphaFold. MODOMICS now includes a data annotation system, with 'Evidence' and 'Estimated Reliability' features, offering clarity on data support and accuracy. This system is open to all MODOMICS entries, enhancing the accuracy of RNA modification data representation. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes.

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

MODOMICS: a database of RNA modification pathways. 2021 update.

The MODOMICS database has been, since 2006, a manually curated and centralized resource, storing and distributing comprehensive information about modified ribonucleosides. Originally, it only contained data on the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. Over the years, prompted by the accumulation of new knowledge and new types of data, it has been updated with new information and functionalities. In this new release, we have created a catalog of RNA modifications linked to human diseases, e.g., due to mutations in genes encoding modification enzymes. MODOMICS has been linked extensively to RCSB Protein Data Bank, and sequences of experimentally determined RNA structures with modified residues have been added. This expansion was accompanied by including nucleotide 5'-monophosphate residues. We redesigned the web interface and upgraded the database backend. In addition, a search engine for chemically similar modified residues has been included that can be queried by SMILES codes or by drawing chemical molecules. Finally, previously available datasets of modified residues, biosynthetic pathways, and RNA-modifying enzymes have been updated. Overall, we provide users with a new, enhanced, and restyled tool for research on RNA modification. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

MODOMICS: a database of RNA modification pathways. 2017 update.

MODOMICS is a database of RNA modifications that provides comprehensive information concerning the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. In the current database version, we included the following new features and data: extended mass spectrometry and liquid chromatography data for modified nucleosides; links between human tRNA sequences and MINTbase - a framework for the interactive exploration of mitochondrial and nuclear tRNA fragments; new, machine-friendly system of unified abbreviations for modified nucleoside names; sets of modified tRNA sequences for two bacterial species, updated collection of mammalian tRNA modifications, 19 newly identified modified ribonucleosides and 66 functionally characterized proteins involved in RNA modification. Data from MODOMICS have been linked to the RNAcentral database of RNA sequences. MODOMICS is available at http://modomics.genesilico.pl.

S6:S18 ribosomal protein complex interacts with a structural motif present in its own mRNA.

Prokaryotic ribosomal protein genes are typically grouped within highly conserved operons. In many cases, one or more of the encoded proteins not only bind to a specific site in the ribosomal RNA, but also to a motif localized within their own mRNA, and thereby regulate expression of the operon. In this study, we computationally predicted an RNA motif present in many bacterial phyla within the 5' untranslated region of operons encoding ribosomal proteins S6 and S18. We demonstrated that the S6:S18 complex binds to this motif, which we hereafter refer to as the S6:S18 complex-binding motif (S6S18CBM). This motif is a conserved CCG sequence presented in a bulge flanked by a stem and a hairpin structure. A similar structure containing a CCG trinucleotide forms the S6:S18 complex binding site in 16S ribosomal RNA. We have constructed a 3D structural model of a S6:S18 complex with S6S18CBM, which suggests that the CCG trinucleotide in a specific structural context may be specifically recognized by the S18 protein. This prediction was supported by site-directed mutagenesis of both RNA and protein components. These results provide a molecular basis for understanding protein-RNA recognition and suggest that the S6S18CBM is involved in an auto-regulatory mechanism.

MODOMICS: a database of RNA modification pathways--2013 update.

MODOMICS is a database of RNA modifications that provides comprehensive information concerning the chemical structures of modified ribonucleosides, their biosynthetic pathways, RNA-modifying enzymes and location of modified residues in RNA sequences. In the current database version, accessible at http://modomics.genesilico.pl, we included new features: a census of human and yeast snoRNAs involved in RNA-guided RNA modification, a new section covering the 5'-end capping process, and a catalogue of 'building blocks' for chemical synthesis of a large variety of modified nucleosides. The MODOMICS collections of RNA modifications, RNA-modifying enzymes and modified RNAs have been also updated. A number of newly identified modified ribonucleosides and more than one hundred functionally and structurally characterized proteins from various organisms have been added. In the RNA sequences section, snRNAs and snoRNAs with experimentally mapped modified nucleosides have been added and the current collection of rRNA and tRNA sequences has been substantially enlarged. To facilitate literature searches, each record in MODOMICS has been cross-referenced to other databases and to selected key publications. New options for database searching and querying have been implemented, including a BLAST search of protein sequences and a PARALIGN search of the collected nucleic acid sequences.

RNA methyltransferases involved in 5' cap biosynthesis.

In eukaryotes and viruses that infect them, the 5' end of mRNA molecules, and also many other functionally important RNAs, are modified to form a so-called cap structure that is important for interactions of these RNAs with many nuclear and cytoplasmic proteins. The RNA cap has multiple roles in gene expression, including enhancement of RNA stability, splicing, nucleocytoplasmic transport, and translation initiation. Apart from guanosine addition to the 5' end in the most typical cap structure common to transcripts produced by RNA polymerase II (in particular mRNA), essentially all cap modifications are due to methylation. The complexity of the cap structure and its formation can range from just a single methylation of the unprocessed 5' end of the primary transcript, as in mammalian U6 and 7SK, mouse B2, and plant U3 RNAs, to an elaborate m(7)Gpppm(6,6)AmpAmpCmpm(3)Um structure at the 5' end of processed RNA in trypanosomes, which are formed by as many as 8 methylation reactions. While all enzymes responsible for methylation of the cap structure characterized to date were found to belong to the same evolutionarily related and structurally similar Rossmann Fold Methyltransferase superfamily, that uses the same methyl group donor, S-adenosylmethionine; the enzymes also exhibit interesting differences that are responsible for their distinct functions. This review focuses on the evolutionary classification of enzymes responsible for cap methylation in RNA, with a focus on the sequence relationships and structural similarities and dissimilarities that provide the basis for understanding the mechanism of biosynthesis of different caps in cellular and viral RNAs. Particular attention is paid to the similarities and differences between methyltransferases from human cells and from human pathogens that may be helpful in the development of antiviral and antiparasitic drugs.

Computational modeling of RNA 3D structures, with the aid of experimental restraints.

In addition to mRNAs whose primary function is transmission of genetic information from DNA to proteins, numerous other classes of RNA molecules exist, which are involved in a variety of functions, such as catalyzing biochemical reactions or performing regulatory roles. In analogy to proteins, the function of RNAs depends on their structure and dynamics, which are largely determined by the ribonucleotide sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore, the majority of known RNAs remain structurally uncharacterized. To address this problem, computational structure prediction methods were developed that simulate either the physical process of RNA structure formation ("Greek science" approach) or utilize information derived from known structures of other RNA molecules ("Babylonian science" approach). All computational methods suffer from various limitations that make them generally unreliable for structure prediction of long RNA sequences. However, in many cases, the limitations of computational and experimental methods can be overcome by combining these two complementary approaches with each other. In this work, we review computational approaches for RNA structure prediction, with emphasis on implementations (particular programs) that can utilize restraints derived from experimental analyses. We also list experimental approaches, whose results can be relatively easily used by computational methods. Finally, we describe case studies where computational and experimental analyses were successfully combined to determine RNA structures that would remain out of reach for each of these approaches applied separately.

2'-O-ribose methylation of cap2 in human: function and evolution in a horizontally mobile family.

The 5' cap of human messenger RNA consists of an inverted 7-methylguanosine linked to the first transcribed nucleotide by a unique 5'-5' triphosphate bond followed by 2'-O-ribose methylation of the first and often the second transcribed nucleotides, likely serving to modify efficiency of transcript processing, translation and stability. We report the validation of a human enzyme that methylates the ribose of the second transcribed nucleotide encoded by FTSJD1, henceforth renamed HMTR2 to reflect function. Purified recombinant hMTr2 protein transfers a methyl group from S-adenosylmethionine to the 2'-O-ribose of the second nucleotide of messenger RNA and small nuclear RNA. Neither N(7) methylation of the guanosine cap nor 2'-O-ribose methylation of the first transcribed nucleotide are required for hMTr2, but the presence of cap1 methylation increases hMTr2 activity. The hMTr2 protein is distributed throughout the nucleus and cytosol, in contrast to the nuclear hMTr1. The details of how and why specific transcripts undergo modification with these ribose methylations remains to be elucidated. The 2'-O-ribose RNA cap methyltransferases are present in varying combinations in most eukaryotic and many viral genomes. With the capping enzymes in hand their biological purpose can be ascertained.

Are there double knots in proteins? Prediction and in vitro verification based on TrmD-Tm1570 fusion from C. nitroreducens.

We have been aware of the existence of knotted proteins for over 30 years-but it is hard to predict what is the most complicated knot that can be formed in proteins. Here, we show new and the most complex knotted topologies recorded to date-double trefoil knots (31 #31). We found five domain arrangements (architectures) that result in a doubly knotted structure in almost a thousand proteins. The double knot topology is found in knotted membrane proteins from the CaCA family, that function as ion transporters, in the group of carbonic anhydrases that catalyze the hydration of carbon dioxide, and in the proteins from the SPOUT superfamily that gathers 31 knotted methyltransferases with the active site-forming knot. For each family, we predict the presence of a double knot using AlphaFold and RoseTTaFold structure prediction. In the case of the TrmD-Tm1570 protein, which is a member of SPOUT superfamily, we show that it folds in vitro and is biologically active. Our results show that this protein forms a homodimeric structure and retains the ability to modify tRNA, which is the function of the single-domain TrmD protein. However, how the protein folds and is degraded remains unknown.

Insights into the structure, function and evolution of the radical-SAM 23S rRNA methyltransferase Cfr that confers antibiotic resistance in bacteria.

The Cfr methyltransferase confers combined resistance to five classes of antibiotics that bind to the peptidyl tranferase center of bacterial ribosomes by catalyzing methylation of the C-8 position of 23S rRNA nucleotide A2503. The same nucleotide is targeted by the housekeeping methyltransferase RlmN that methylates the C-2 position. Database searches with the Cfr sequence have revealed a large group of closely related sequences from all domains of life that contain the conserved CX(3)CX(2)C motif characteristic of radical S-adenosyl-l-methionine (SAM) enzymes. Phylogenetic analysis of the Cfr/RlmN family suggests that the RlmN subfamily is likely the ancestral form, whereas the Cfr subfamily arose via duplication and horizontal gene transfer. A structural model of Cfr has been calculated and used as a guide for alanine mutagenesis studies that corroborate the model-based predictions of a 4Fe-4S cluster, a SAM molecule coordinated to the iron-sulfur cluster (SAM1) and a SAM molecule that is the putative methyl group donor (SAM2). All mutations at predicted functional sites affect Cfr activity significantly as assayed by antibiotic susceptibility testing and primer extension analysis. The investigation has identified essential amino acids and Cfr variants with altered reaction mechanisms and represents a first step towards understanding the structural basis of Cfr activity.

LEF1/beta-catenin complex regulates transcription of the Cav3.1 calcium channel gene (Cacna1g) in thalamic neurons of the adult brain.

beta-Catenin, together with LEF1/TCF transcription factors, activates genes involved in the proliferation and differentiation of neuronal precursor cells. In mature neurons, beta-catenin participates in dendritogenesis and synaptic function as a component of the cadherin cell adhesion complex. However, the transcriptional activity of beta-catenin in these cells remains elusive. In the present study, we found that in the adult mouse brain, beta-catenin and LEF1 accumulate in the nuclei of neurons specifically in the thalamus. The particular electrophysiological properties of thalamic neurons depend on T-type calcium channels. Cav3.1 is the predominant T-type channel subunit in the thalamus, and we hypothesized that the Cacna1g gene encoding Cav3.1 is a target of the LEF1/beta-catenin complex. We demonstrated that the expression of Cacna1g is high in the thalamus and is further increased in thalamic neurons treated in vitro with LiCl or WNT3A, activators of beta-catenin. Luciferase reporter assays confirmed that the Cacna1G promoter is activated by LEF1 and beta-catenin, and footprinting analysis revealed four LEF1 binding sites in the proximal region of this promoter. Chromatin immunoprecipitation demonstrated that the Cacna1g proximal promoter is occupied by beta-catenin in vivo in the thalamus, but not in the hippocampus. Moreover, WNT3A stimulation enhanced T-type current in cultured thalamic neurons. Together, our data indicate that the LEF1/beta-catenin complex regulates transcription of Cacna1g and uncover a novel function for beta-catenin in mature neurons. We propose that beta-catenin contributes to neuronal excitability not only by a local action at the synapse but also by activating gene expression in thalamic neurons.

MODOMICS: a database of RNA modification pathways. 2008 update.

MODOMICS, a database devoted to the systems biology of RNA modification, has been subjected to substantial improvements. It provides comprehensive information on the chemical structure of modified nucleosides, pathways of their biosynthesis, sequences of RNAs containing these modifications and RNA-modifying enzymes. MODOMICS also provides cross-references to other databases and to literature. In addition to the previously available manually curated tRNA sequences from a few model organisms, we have now included additional tRNAs and rRNAs, and all RNAs with 3D structures in the Nucleic Acid Database, in which modified nucleosides are present. In total, 3460 modified bases in RNA sequences of different organisms have been annotated. New RNA-modifying enzymes have been also added. The current collection of enzymes includes mainly proteins for the model organisms Escherichia coli and Saccharomyces cerevisiae, and is currently being expanded to include proteins from other organisms, in particular Archaea and Homo sapiens. For enzymes with known structures, links are provided to the corresponding Protein Data Bank entries, while for many others homology models have been created. Many new options for database searching and querying have been included. MODOMICS can be accessed at http://genesilico.pl/modomics.

YgdE is the 2'-O-ribose methyltransferase RlmM specific for nucleotide C2498 in bacterial 23S rRNA.

The rRNAs of Escherichia coli contain four 2'-O-methylated nucleotides. Similar to other bacterial species and in contrast with Archaea and Eukaryota, the E. coli rRNA modifications are catalysed by specific methyltransferases that find their nucleotide targets without being guided by small complementary RNAs. We show here that the ygdE gene encodes the methyltransferase that catalyses 2'-O-methylation at nucleotide C2498 in the peptidyl transferase loop of E. coli 23S rRNA. Analyses of rRNAs using MALDI mass spectrometry showed that inactivation of the ygdE gene leads to loss of methylation at nucleotide C2498. The loss of ygdE function causes a slight reduction in bacterial fitness. Methylation at C2498 was restored by complementing the knock-out strain with a recombinant copy of ygdE. The recombinant YgdE methyltransferase modifies C2498 in naked 23S rRNA, but not in assembled 50S subunits or ribosomes. Nucleotide C2498 is situated within a highly conserved and heavily modified rRNA sequence, and YgdE's activity is influenced by other modification enzymes that target this region. Phylogenetically, YgdE is placed in the cluster of orthologous groups COG2933 together with S-adenosylmethionine-dependent, Rossmann-fold methyltransferases such as the archaeal and eukaryotic RNA-guided fibrillarins. The ygdE gene has been redesignated rlmM for rRNA large subunit methyltransferase M.

YbeA is the m3Psi methyltransferase RlmH that targets nucleotide 1915 in 23S rRNA.

Pseudouridines in the stable RNAs of Bacteria are seldom subjected to further modification. There are 11 pseudouridine (Psi) sites in Escherichia coli rRNA, and further modification is found only at Psi1915 in 23S rRNA, where the N-3 position of the base becomes methylated. Here, we report the identity of the E. coli methyltransferase that specifically catalyzes methyl group addition to form m(3)Psi1915. Analyses of E. coli rRNAs using MALDI mass spectrometry showed that inactivation of the ybeA gene leads to loss of methylation at nucleotide Psi1915. Methylation is restored by complementing the knockout strain with a plasmid-encoded copy of ybeA. Homologs of the ybeA gene, and thus presumably the ensuing methylation at nucleotide m(3)Psi1915, are present in most bacterial lineages but are essentially absent in the Archaea and Eukaryota. Loss of ybeA function in E. coli causes a slight slowing of the growth rate. Phylogenetically, ybeA and its homologs are grouped with other putative S-adenosylmethionine-dependent, SPOUT methyltransferase genes in the Cluster of Orthologous Genes COG1576; ybeA is the first member to be functionally characterized. The YbeA methyltransferase is active as a homodimer and docks comfortably into the ribosomal A site without encroaching into the P site. YbeA makes extensive interface contacts with both the 30S and 50S subunits to align its active site cofactor adjacent to nucleotide Psi1915. Methylation by YbeA (redesignated RlmH for rRNA large subunit methyltransferase H) possibly functions as a stamp of approval signifying that the 50S subunit has engaged in translational initiation.

Functional specialization of domains tandemly duplicated within 16S rRNA methyltransferase RsmC.

RNA methyltransferases (MTases) are important players in the biogenesis and regulation of the ribosome, the cellular machine for protein synthesis. RsmC is a MTase that catalyzes the transfer of a methyl group from S-adenosyl-l-methionine (SAM) to G1207 of 16S rRNA. Mutations of G1207 have dominant lethal phenotypes in Escherichia coli, underscoring the significance of this modified nucleotide for ribosome function. Here we report the crystal structure of E. coli RsmC refined to 2.1 A resolution, which reveals two homologous domains tandemly duplicated within a single polypeptide. We characterized the function of the individual domains and identified key residues involved in binding of rRNA and SAM, and in catalysis. We also discovered that one of the domains is important for the folding of the other. Domain duplication and subfunctionalization by complementary degeneration of redundant functions (in particular substrate binding versus catalysis) has been reported for many enzymes, including those involved in RNA metabolism. Thus, RsmC can be regarded as a model system for functional streamlining of domains accompanied by the development of dependencies concerning folding and stability.

Sequence-structure-function analysis of the bifunctional enzyme MnmC that catalyses the last two steps in the biosynthesis of hypermodified nucleoside mnm5s2U in tRNA.

MnmC catalyses the last two steps in the biosynthesis of 5-methylaminomethyl-2-thiouridine (mnm(5)s(2)U) in tRNA. Previously, we reported that this bifunctional enzyme is encoded by the yfcK open reading frame in the Escherichia coli K12 genome. However, the mechanism of its activity, in particular the potential structural and functional dependence of the domains responsible for catalyzing the two modification reactions, remains unknown. With the aid of the protein fold-recognition method, we constructed a structural model of MnmC in complex with the ligands and target nucleosides and studied the role of individual amino acids and entire domains by site-directed and deletion mutagenesis, respectively. We found out that the N-terminal domain contains residues responsible for binding of the S-adenosylmethionine cofactor and catalyzing the methylation of nm(5)s(2)U to form mnm(5)s(2)U, while the C-terminal domain contains residues responsible for binding of the FAD cofactor. Further, point mutants with compromised activity of either domain can complement each other to restore a fully functional enzyme. Thus, in the conserved fusion protein MnmC, the individual domains retain independence as enzymes. Interestingly, the N-terminal domain is capable of independent folding, while the isolated C-terminal domain is incapable of folding on its own, a situation similar to the one reported recently for the rRNA modification enzyme RsmC.

Human RNA cap1 methyltransferase CMTr1 cooperates with RNA helicase DHX15 to modify RNAs with highly structured 5' termini.

The 5'-cap structure, characteristic for RNA polymerase II-transcribed RNAs, plays important roles in RNA metabolism. In humans, RNA cap formation includes post-transcriptional modification of the first transcribed nucleotide by RNA cap1 methyltransferase (CMTr1). Here, we report that CMTr1 activity is hindered towards RNA substrates with highly structured 5' termini. We found that CMTr1 binds ATP-dependent RNA DHX15 helicase and that this interaction, mediated by the G-patch domain of CMTr1, has an advantageous effect on CMTr1 activity towards highly structured RNA substrates. The effect of DHX15 helicase activity is consistent with the strength of the secondary structure that has to be removed for CMTr1 to access the 5'-terminal residues in a single-stranded conformation. This is, to our knowledge, the first demonstration of the involvement of DHX15 in post-transcriptional RNA modification, and the first example of a molecular process in which DHX15 directly affects the activity of another enzyme. Our findings suggest a new mechanism underlying the regulatory role of DHX15 in the RNA capping process. RNAs with highly structured 5' termini constitute a significant fraction of the human transcriptome. Hence, CMTr1-DHX15 cooperation is likely to be important for the metabolism of RNA polymerase II-transcribed RNAs.This article is part of the theme issue '5' and 3' modifications controlling RNA degradation'.

Novel inhibitors of the rRNA ErmC' methyltransferase to block resistance to macrolides, lincosamides, streptogramine B antibiotics.

In erythromycin-resistant bacteria, the N6 position of A2058 in 23S rRNA is mono- or dimethylated by Erm family methyltransferases. This modification results in cross-resistance to macrolides, lincosamides and streptogramin B. Most inhibitors of Erm methyltransferases developed up-to-date target the cofactor-binding pocket, resulting in a lack of selectivity whereas inhibitors that bind the substrate-binding pocket demonstrate low in vitro activity. In this study, a molecular docking approach followed by biochemical screening was applied to search for inhibitors targeting both cofactor- and substrate-binding pockets of ErmC' methyltransferase. Based on the results of the molecular docking-based virtual screening of the clean-leads subset of the ZINC database, 29 compounds were chosen for experimental verification. Among them inhibitor 28 (ZINC code 32747906), with an IC50 of 100 µM, decreased the minimal inhibitory concentration of erythromycin in the Escherichia coli strain overexpressing ErmC'. Docking analysis of 28 to the ErmC' structure and the competitive ligand binding assay revealed a non-competitive model of inhibition. Inhibitor 28 served as a template for similarity-based virtual screening, which resulted in the identification of two derivatives 3s (ZINC code 62022572) and 4s (ZINC code 49032257) with an IC50 of 116 µM and 110 µM, respectively. Our results provide a basis for the development of inhibitors against the Erm-family of enzymes.

Structural and evolutionary bioinformatics of the SPOUT superfamily of methyltransferases.

BACKGROUND: SPOUT methyltransferases (MTases) are a large class of S-adenosyl-L-methionine-dependent enzymes that exhibit an unusual alpha/beta fold with a very deep topological knot. In 2001, when no crystal structures were available for any of these proteins, Anantharaman, Koonin, and Aravind identified homology between SpoU and TrmD MTases and defined the SPOUT superfamily. Since then, multiple crystal structures of knotted MTases have been solved and numerous new homologous sequences appeared in the databases. However, no comprehensive comparative analysis of these proteins has been carried out to classify them based on structural and evolutionary criteria and to guide functional predictions. RESULTS: We carried out extensive searches of databases of protein structures and sequences to collect all members of previously identified SPOUT MTases, and to identify previously unknown homologs. Based on sequence clustering, characterization of domain architecture, structure predictions and sequence/structure comparisons, we re-defined families within the SPOUT superfamily and predicted putative active sites and biochemical functions for the so far uncharacterized members. We have also delineated the common core of SPOUT MTases and inferred a multiple sequence alignment for the conserved knot region, from which we calculated the phylogenetic tree of the superfamily. We have also studied phylogenetic distribution of different families, and used this information to infer the evolutionary history of the SPOUT superfamily. CONCLUSION: We present the first phylogenetic tree of the SPOUT superfamily since it was defined, together with a new scheme for its classification, and discussion about conservation of sequence and structure in different families, and their functional implications. We identified four protein families as new members of the SPOUT superfamily. Three of these families are functionally uncharacterized (COG1772, COG1901, and COG4080), and one (COG1756 represented by Nep1p) has been already implicated in RNA metabolism, but its biochemical function has been unknown. Based on the inference of orthologous and paralogous relationships between all SPOUT families we propose that the Last Universal Common Ancestor (LUCA) of all extant organisms contained at least three SPOUT members, ancestors of contemporary RNA MTases that carry out m1G, m3U, and 2'O-ribose methylation, respectively. In this work we also speculate on the origin of the knot and propose possible 'unknotted' ancestors. The results of our analysis provide a comprehensive 'roadmap' for experimental characterization of SPOUT MTases and interpretation of functional studies in the light of sequence-structure relationships.