BACKGROUND
tuberculosis  is an infectious disease that mainly afflicts populations in third-world countries.
although most infected people will never develop an active form of the disease, the global tb death rate is substantial, being around two million people per year.
the situation has been further exacerbated by the emergence and spread of extremely resistant mycobacterium tuberculosis strains   <cit>  which, together with the lack of effective antibiotics, urges the development of new alternatives for controlling this worldwide public health threat.
consequently, most research efforts have focused on understanding m. tuberculosis biology with the aim of identifying new therapeutic and vaccine targets.
current management of tb cases consists of a two-month intensive treatment comprising first-line drugs such as rifampicin , isoniazid , pyrazinamide  and ethambutol /streptomycin  which seeks to ensure that mono-resistant strains do not proliferate  <cit> .
when rif and inh are ineffective, standard treatment guidelines recommend a combination of second-line drugs  <cit> .
among these therapeutic choices, aminoglycosides such as capreomycin  and viomycin  are given in combination with other antibiotics to treat multidrug-resistant strains ; they are also effective against non-replicating m. tuberculosis strains, as well as being useful in the treatment of latent tb infections  <cit> .
cmn- and vmn-resistant bacilli are classified as xdr-tb strains and are a cause of major concern since these drugs have more toxic side effects and might result in higher death rates, especially among hiv-infected persons  <cit> .
controversially, m. tuberculosis tlya gene product has been annotated both as a virulence factor, due to its ability to lyse red blood cells, as well as a protein involved in ribosome biogenesis  <cit> .
tlya has been assumed to induce hemolytic activity due to its high amino acid similarity  to the pore-forming hemolysin/cytotoxin virulence determinant from serpulina hyodysenteriae.
notably such activity was only experimentally demonstrated with bacterial lysates  <cit>  and another study has even suggested that tlya could be associated with hemolysin expression in escherichia coli  <cit> .
additionally, recent study showed that a concentration-dependent hemolysis of rabbit and human erythrocytes to be induced on incubation with the recombinant rv <dig> protein  when over-expressed in e. coli  <cit> .
the article reported hemolysis inhibition using specific antiserum against the same protein as well as the presence of the recombinant protein on e. coli cell wall; these data were supportive of tlya's hemolytic property, combined with a preliminary in silico analysis.
rv <dig> oligomers were also observed on lysed erythrocyte membranes, as well as the susceptibility of rv1694-expressing e. coli to cmn  <cit> .
despite such thorough in vitro analysis in e. coli, the functional activities suggested for rv <dig> have still not been clearly demonstrated when using m. tuberculosis mutant strains.
by contrast with the proposed hemolytic function, johansen et al.,  <cit>  have reported that m. tuberculosis h37rv and beijing tlya confer susceptibility to cmn and vmn.
these antibiotics are structurally similar and share the same mechanism of action, both inhibiting bacterial proliferation by blocking the peptidyl transferase reaction.
vmn hinders a-site and p-site trna positioning on the 50s ribosomal subunit, therefore hampering mrna translation  <cit> .
cmn resistance arises due to mutations in the tlya gene which encodes a rna 2'-o-methyltransferase  that methylates ribose moieties from nucleotide c <dig> in 16s rrna and c <dig> in 23s rrna.
the absence of these rrna methylations confers resistance against cmn and vmn given that binding sites for these antibiotics are delineated by the localization of such methylated riboses  <cit> .
additional evidence has been provided by mutant complementation, showing that cmn susceptibility is restored by introducing an active copy of the tlya gene  <cit> .
tlya methylation clearly enhances cmn and vmn susceptibility and is responsible, at least in part, for their efficacy against mycobacterial infections as drug binding is impaired when the enzyme is not produced.
resistance thus emerges when tlya activity is lost.
these studies show the need for a reevaluation of tlya's function in virulence  <cit> .
RESULTS
tlya related sequences
the blast search for tlya-related sequences in the tb and genbank databases showed that although tlya-related sequences are found mainly in prokaryotic organisms, homologous protein sequences are also present in eukaryotic organisms such as plants and algae .
these m. tuberculosis-related tlya protein sequences are mainly classified as hemolysin a inside the ftsj-like class based on domain similarity and have been annotated as such in genbank, embl, uniprot and swiss-prot databases according to sequence similarity with serpulina hyodysenteriae tlya annotated as being a hemolysin.
paradoxically, although hemolysin a contains similar domain architecture to tlya, it is predicted to be a rrna methylase instead of a contact-dependent hemolysin .
tlya amino acid sequences from m. tuberculosis h37rv, h37ra, cdc <dig>  m. bovis and m. bovis-bcg reference strain alignment showed a 100% identity with the harleem, f <dig> and c m. tuberculosis strains .
the same conservation pattern was obtained when the corresponding nucleotide sequences were aligned .
similarity percentages found in the multiple amino acid alignment of mycobacterium species different to m. tuberculosis, such as m. smegmatis , m. abscessus , m. leprae , m. avium  and m. vanbaalenii pyr- <dig> , showed high tlya conservation at genus level.
a glutamine residue insertion was found at position  <dig> in the protein's c-terminal region in all the aforementioned species, while a 12-residue insertion was found only in m. ulcerans  and m. marinum , which are pathogenic mycobacteria .
tlya transcriptional unit revealed some rearrangements in the mycobacterium iprj-recn operon
from the synteny of the iprj-recn operon from m. tuberculosis, m. bovis bcg, m. bovis, m. leprae, m. ulcerans and m. marinum  it was evident that almost all gene products, except for tlya, are involved in intermediary metabolism, cell wall biosynthesis and signal transduction rather than virulence.
moreover, no trafficking proteins were found within the operon or near tlya that could be suggestive of tlya transportation either to the cell membrane or the extracellular milieu.
the iprj-recn operon also had common gene organization in m. tuberculosis and m. bovis, while in m. leprae most genes have lost their function .
such high rate of loss of function in genes involved in intermediary metabolism, cellular respiration and cell-wall biosynthesis may suggest that the loss of such genes confers an adaptive advantage on mycobacterium during host cell infection or transmission  <cit> .
no other genes reported to be involved in bacilli virulence were found within this operon.
tlya three-dimensional structure
the three-dimensional structure of a putative hemolysin from streptococcus thermophilus  was used as template for homology modeling.
the obtained identity  suggested that the selected hemolysin structure was a suitable template for tlya, and that the obtained modeled structure could reflect an experimentally-obtained m. tuberculosis tlya structure.
although this protein seemed to be a hemolysin, its domain architecture was also similar to rna methyltransferase family, fibrillarin homologues and tlya proteins.
a consensus sequence was established which showed the presence of a common core, comprising six β-sheets, the first five of which were found parallel and the sixth anti-parallel between two layers of α-helices and random coiled regions.
such structural organization is commonly referred as to an adomet-dependent methyltransferase  fold  <cit> .
tlya has an overall globular fold and comprises two domains exhibiting different spatial arrangements.
the smaller domain, named s <dig>  is n-terminally located and is overlapped by the catalytic domain which exhibits structural similarity to various unrelated rna-binding proteins, while the large catalytic domain, named ftsj-like methyltransferase domain is common to 2'-o-mtases and exhibits a α/β fold with a deep pocket .
the s <dig> domain consists of  <dig> residues, starting at a- <dig> and ending at h- <dig>  this domain is structurally formed by one β-sheet and three α-helices predominantly composed of positively charged residues  whose probable function is to mediate rna binding and provide stability.
in line with this, the rnabindr server identified three short sequences  inside the s <dig> domain that might be implicated in rna binding .
the s <dig> domain has been specifically identified in some bacterial and eukaryotic ribosomal proteins, pseudouridine synthases, rna methylases, bacterial tyrosyl-trna synthetases and may also be involved in translation regulation  <cit> .
the analysis of the protein's electrostatic surface showed that positively charged residues were localized towards the n-terminal extreme, forming bulges and long finger-like projections that probably extend into the rrna core to stabilize its structure .
in fact, charged residues like r and k exhibited the highest interface propensities, which is consistent with their ability to participate in interactions both with nucleotide bases and the negatively charged rna phosphate backbone  <cit> .
the ftsj-like methyltransferase domain is located between residues  <dig> and  <dig> and includes the catalytic tetrad responsible for tlya's enzymatic activity  <cit> .
moreover, the ftsj-like methyltransferase domain exhibits a rossmann fold composed of five parallel β-strand layers linked by two α-helices following an α/β secondary structure topological order  <cit> .
our comparative analysis of tlya with the fibrillarin family showed the presence of a gastg adomet binding motif comprising residues 90- <dig> and a d located at position  <dig> juxtaposed to the above mentioned motif , which might stabilize tlya coenzyme binding  <cit> .
adomet was located between β-strands  <dig> and  <dig> within the inner depression inside the catalytic pocket formed by gastg residues' main-chain nh groups .
moreover, substrate positioning might be driven by the rossmann configuration fold and the pertinent adomet-binding/catalytic sites, suggesting substrate selectivity in the small ribosomal subunit  <cit> .
likewise, domain structure was in agreement with s-adenosyl-methionine -dependent methyltransferases.
adomet binding and affinity might indicate how an enzyme such as tlya may methylate the two adjacent ribose rings on both ribosomal subunits.
in this reaction s-adenosyl-methionine  acts as methyl donor, as it is converted into s-adenosyl-l-homocysteine ; tlya could thus be also classified as being an rna-adomet-mtase class i  <cit> .
this hypothesis is supported by the high structural similarity  between tlya and ftsj rna-mtase complexed with adomet   .
two expected structural motifs were also identified in this structure which was formed by six buried residues between g- <dig> and f- <dig> which are pivotal for adomet interaction and therefore conserved in mycobacterium tlya proteins .
nest analysis of two structural motifs detected in tlya functionally important regions for adomet binding region.
results were obtained on profunc server  <cit> .
the tlya model was validated by using the rampage server which evaluates a protein's 3d-structure based on the stereochemistry quality of torsion angles and geometry  <cit> .
the ramachandran plot showed  <dig> % feasibility  in favorable positions and  <dig> %  in acceptable regions, with just  <dig> % in an outlier region , supporting the high quality of the 3d-model .
this in silico structure analysis of m. tuberculosis proteins allowed a deeper understanding of tlya's function and challenges the initially described function of tlya-related proteins.
rather than its hypothetical function as a hemolysin, our data suggested functions as an mtase based on in silico prediction.
predicting subcellular localization
subcellular localization was predicted with tbpred  and the results suggested tlya localization in the cytoplasmatic compartment , also supporting the fact that tlya acts as a methyltransferase  <cit> .
a recent mycobacterial protein analysis led to ascertaining several predictors' reliability for subcellular localization, providing highly accurate results for gpos-ploc and pa-sub v. <dig>  and signalp  <dig>  as a confirmatory tool  <cit> .
however, in this study, pa-sub v. <dig>  predicted tlya as being an extracellular protein , which could have been biased as the training strategy used involved results from the swiss-prot  <cit>  database where tlya is mis-annotated.
likewise, in agreement with blast scores for the query sequence ,  <dig> out of  <dig> homologous proteins were designated as putative ribosomal rna methyltransferases and rna binding proteins.
only one submitted sequence seemed to be similar to hemolysin a .
nd: non defined value.
pse-pssm: pseudo position-specific score matrix.
svm: support vector machines.
pssm: position-specific scoring matrix.
meme/mast: multiple em for motif elicitation/motif alignment and search tool.
hm: hybrid module.
pseaac: pseudo amino acid composition.
evolutionary insights into tlya genes
evolutionary analysis of tlya genes from  <dig> mycobacterial species,  <dig> species taxonomically related to mycobacterium,  <dig> adomet-rna mtases and  <dig> pfp showed similar topology by both neighbor-joining  and maximum parsimony  methods, suggesting that the phylogenetic relationship observed was highly reliable .
the consensus tree showed the tlya gene in m. tuberculosis h37rv, cdc <dig>  c, haarlem, beijing, kzn and inh-rif strains as well as in the m. bovis af <dig> and bcg strains to fit in a monophyletic cluster.
the tlya tree exhibited short genetic distances  for m. tuberculosis complex species, including those separated according to clinical considerations.
phylogenetic analysis showed that mycobacterium tlya genes are closely related to those of nocardia and rhodococcus genus; such branching patterns were supported by high bootstrap values and agreed with the results obtained from 16s rdna phylogeny.
no differences were observed between antibiotic susceptible and resistant m. tuberculosis strains.
mycobacterium tlya proteins fell into  <dig> groups; each group formed a well-delineated branch supported by high bootstrap values.
m. abscessus was the most distant and tlya-ancestral species from the mycobacterium genus compared to the closest species to m. tuberculosis .
serpulina hyodysenteriae and lawsonia intracellularis lsaa antigen sequences were included due to their sequence identity to m. tuberculosis tlya , sharing a common ancestor, although having experimentally-attributed different roles  <cit>  .
pattern branching of the phylogenetic tree showed that the tlya protein formed a monophyletic cluster with adomet rna mtases, suggesting a common evolutionary origin; these findings might also imply a similar functional role in translational machinery.
pfp formed a distant cluster where diverse proteins were clustered on a branch with the most divergent proteins .
the synonymous and non-synonymous substitutions were calculated to identify the action of neutral, positive or purifying selection on tlya genes; their synonymous, pi, and non-synonymous diversity, pi, was  <dig>  and  <dig> , respectively, corresponding to a  <dig>  pi/pi rate which can be inferred as purifying selection acting on this gene among mycobacterium species.
this suggests that conservation at sequence level might show mycobacterium tlya gene homogeneity and show low inter-specific variation between the species compared in this study.
discussion
the aim of the present study was to describe structural and phylogenetic features of m. tuberculosis tlya, using bioinformatics tools.
the structural and phylogenetic analysis shown here highlights the need for re-annotating this protein family and redefining it as an rna-binding ftsj-like methyltransferase in bioinformatics databases; this is based on the high degree of amino acid sequence similarity found in the mycobacterium species analyzed here.
it can be surmised that physical organization, domain spatial arrangement and protein folding on tlya is structurally conserved.
previous sequence analysis of 2'-o-mtase families has indicated that they use a similar catalytic mechanism and have inherited a common function from the same ancestral ribose 2'-o-mtase.
feder et al., suggested that tlya is a 2'-o-mtase, identifying a catalytic tetrad consisting of lysine- <dig>  aspartic acid- <dig>  lysine- <dig> and glutamic acid in position  <dig>   <cit> .
based on amino acid sequence alignment, our comparative analysis with the rrmj family and other related rna 2'-o-mtases  <cit>  showed changes in d154-k182-e <dig> residue positions .
however, despite being located in different sequence positions, the scaffold formed by these residues was maintained at the catalytic site in the modeled protein 3d structure without affecting protein function  <cit> .
despite the diversity of existing mtases, the majority maintain the same catalytic mechanism where k- <dig> binds to a hydrogen from the ribose 2-oh' group; this nucleophilic attack generates a transition state .
the lysine side chain shifts e- <dig> so that its negatively charged carboxyl group becomes available to promote the attack of the methyl donor  <cit> ; furthermore, this conserved amino acid residue pattern was also observed in our evolutionary trace analysis.
tlya three-dimensional structure resembled an rna methyltransferase
the s. hyodysenteriae tlya postulated mechanism for hemolysin action is said to be similar to that for pore-forming protein  mechanisms  <cit> ; tlya folding, however, is inconsistent with pfps due its lack of membrane-binding domains, cholesterol-recognition and insertion motifs, and other characteristic features such as a conserved undecapeptide  near the c- or n- terminal ends.
it also lacks repeat domains in the toxin  and calcium-binding glycine-rich motifs necessary for oligomerization  <cit> .
this in silico structure analysis of m. tuberculosis proteins led to a deeper understanding of tlya's role and challenged the initially described function for tlya-related proteins.
a combination of bioinformatics analysis and experimental structure elucidation may suggest alternative functional roles for previously-annotated enzymes, as with -dependent methyltransferase, phosphatase and n-acetyltransferase structures from m. tuberculosis  <cit> .
theoretical evidence has thus been provided for tlya's role in rna modification instead of it just being a hemolytic factor.
subcellular localization suggested a cytoplasmatic localization consistent with ribosome-associated activity
predicting subcellular localization is important since several characteristics can be deduced, such as protein function and genome annotation.
it also aids experimental design for proteomics platforms, particularly for identifying new candidates for vaccine development and drug targets  <cit> .
protein sequence analysis did not reveal a signal peptide, suggesting that tlya is not secreted; however, some authors have described tlya as being a probable pfp  <cit> , supported by the presence of over-expressed tlya on the e. coli wall surface  <cit> .
tlya lacks transmembrane helices and β-barrels , suggesting that this protein might not be a membrane-embedded or contact-dependent hemolysin, as has previously been proposed  <cit> .
tlya also has a hydrophobic region in the c-terminal domain  which may not form a membrane-spanning α-helix.
nonetheless, such a region might allow the adenine present in adomet to fit properly inside the catalytic pocket, thus suggesting a different role for these hydrophobic residues.
virulence-related functional annotation, based on all the above-mentioned data, remains controversial.
conservation of tlya genes and functions was expected according to our phylogenetic and evolutionary trace analyses, the initial hypothesis being that homologous proteins have related functions.
however, the extent to which this is true has not been assessed in detail up to now; annotation based on protein homology might thus lead to finding unexpected variations in function.
CONCLUSIONS
the main concern of our work was to decipher structural information about the tlya family by homology comparison.
our structural analysis and finding residue conservation within the active site supports recent experimental work  <cit>  about tlya's role as rna 2`-o-mtase in m. tuberculosis.
this is further supported by the domain composition also suggesting methyltransferase activity and rna-binding.
the distribution of charged residues also defines a likely rna interaction instead of hemolytic activity, little similarity being exhibited with hemolysin proteins, perhaps representing inaccurate annotations in several databases, as observed with other protein families  <cit> .
since genome annotation is based on sequence comparison, it is the most commonly used approach for determining functional homology.
annotations are assigned if proteins surpass similarity cut-off and the propagation of original errors could be increased by each new entry.
in such cases, protein function needs to be inferred from their common three-dimensional structures and manual curation.

 tuberculosis is an infectious disease that mainly afflicts populations in third-world countries. although most infected people will never develop an active form of the disease, the global tb death rate is substantial, being around two million people per year. the situation has been further exacerbated by the emergence and spread of extremely resistant mycobacterium tuberculosis strains which, together with the lack of effective antibiotics, urges the development of new alternatives for controlling this worldwide public health threat. consequently, most research efforts have focused on understanding m. tuberculosis biology with the aim of identifying new therapeutic and vaccine targets. current management of tb cases consists of a two-month intensive treatment comprising first-line drugs such as rifampicin, isoniazid, pyrazinamide and ethambutol /streptomycin which seeks to ensure that mono-resistant strains do not proliferate. when rif and inh are ineffective, standard treatment guidelines recommend a combination of second-line drugs. among these therapeutic choices, aminoglycosides such as capreomycin and viomycin are given in combination with other antibiotics to treat multidrug-resistant strains; they are also effective against non-replicating m. tuberculosis strains, as well as being useful in the treatment of latent tb infections. cmn- and vmn-resistant bacilli are classified as xdr-tb strains and are a cause of major concern since these drugs have more toxic side effects and might result in higher death rates, especially among hiv-infected persons. controversially, m. tuberculosis tlya gene product has been annotated both as a virulence factor, due to its ability to lyse red blood cells, as well as a protein involved in ribosome biogenesis. tlya has been assumed to induce hemolytic activity due to its high amino acid similarity to the pore-forming hemolysin/cytotoxin virulence determinant from serpulina hyodysenteriae. notably such activity was only experimentally demonstrated with bacterial lysates and another study has even suggested that tlya could be associated with hemolysin expression in escherichia coli. additionally, recent study showed that a concentration-dependent hemolysis of rabbit and human erythrocytes to be induced on incubation with the recombinant rv protein when over-expressed in e. coli. the article reported hemolysis inhibition using specific antiserum against the same protein as well as the presence of the recombinant protein on e. coli cell wall; these data were supportive of tlya's hemolytic property, combined with a preliminary in silico analysis. rv oligomers were also observed on lysed erythrocyte membranes, as well as the susceptibility of rv1694-expressing e. coli to cmn. despite such thorough in vitro analysis in e. coli, the functional activities suggested for rv have still not been clearly demonstrated when using m. tuberculosis mutant strains. by contrast with the proposed hemolytic function, johansen et al., have reported that m. tuberculosis h37rv and beijing tlya confer susceptibility to cmn and vmn. these antibiotics are structurally similar and share the same mechanism of action, both inhibiting bacterial proliferation by blocking the peptidyl transferase reaction. vmn hinders a-site and p-site trna positioning on the 50s ribosomal subunit, therefore hampering mrna translation. cmn resistance arises due to mutations in the tlya gene which encodes a rna 2'-o-methyltransferase that methylates ribose moieties from nucleotide c in 16s rrna and c in 23s rrna. the absence of these rrna methylations confers resistance against cmn and vmn given that binding sites for these antibiotics are delineated by the localization of such methylated riboses. additional evidence has been provided by mutant complementation, showing that cmn susceptibility is restored by introducing an active copy of the tlya gene. tlya methylation clearly enhances cmn and vmn susceptibility and is responsible, at least in part, for their efficacy against mycobacterial infections as drug binding is impaired when the enzyme is not produced. resistance thus emerges when tlya activity is lost. these studies show the need for a reevaluation of tlya's function in virulence.  tlya related sequences the blast search for tlya-related sequences in the tb and genbank databases showed that although tlya-related sequences are found mainly in prokaryotic organisms, homologous protein sequences are also present in eukaryotic organisms such as plants and algae. these m. tuberculosis-related tlya protein sequences are mainly classified as hemolysin a inside the ftsj-like class based on domain similarity and have been annotated as such in genbank, embl, uniprot and swiss-prot databases according to sequence similarity with serpulina hyodysenteriae tlya annotated as being a hemolysin. paradoxically, although hemolysin a contains similar domain architecture to tlya, it is predicted to be a rrna methylase instead of a contact-dependent hemolysin. tlya amino acid sequences from m. tuberculosis h37rv, h37ra, cdc m. bovis and m. bovis-bcg reference strain alignment showed a 100% identity with the harleem, f and c m. tuberculosis strains. the same conservation pattern was obtained when the corresponding nucleotide sequences were aligned. similarity percentages found in the multiple amino acid alignment of mycobacterium species different to m. tuberculosis, such as m. smegmatis, m. abscessus, m. leprae, m. avium and m. vanbaalenii pyr-, showed high tlya conservation at genus level. a glutamine residue insertion was found at position in the protein's c-terminal region in all the aforementioned species, while a 12-residue insertion was found only in m. ulcerans and m. marinum, which are pathogenic mycobacteria. tlya transcriptional unit revealed some rearrangements in the mycobacterium iprj-recn operon from the synteny of the iprj-recn operon from m. tuberculosis, m. bovis bcg, m. bovis, m. leprae, m. ulcerans and m. marinum it was evident that almost all gene products, except for tlya, are involved in intermediary metabolism, cell wall biosynthesis and signal transduction rather than virulence. moreover, no trafficking proteins were found within the operon or near tlya that could be suggestive of tlya transportation either to the cell membrane or the extracellular milieu. the iprj-recn operon also had common gene organization in m. tuberculosis and m. bovis, while in m. leprae most genes have lost their function. such high rate of loss of function in genes involved in intermediary metabolism, cellular respiration and cell-wall biosynthesis may suggest that the loss of such genes confers an adaptive advantage on mycobacterium during host cell infection or transmission. no other genes reported to be involved in bacilli virulence were found within this operon. tlya three-dimensional structure the three-dimensional structure of a putative hemolysin from streptococcus thermophilus was used as template for homology modeling. the obtained identity suggested that the selected hemolysin structure was a suitable template for tlya, and that the obtained modeled structure could reflect an experimentally-obtained m. tuberculosis tlya structure. although this protein seemed to be a hemolysin, its domain architecture was also similar to rna methyltransferase family, fibrillarin homologues and tlya proteins. a consensus sequence was established which showed the presence of a common core, comprising six β-sheets, the first five of which were found parallel and the sixth anti-parallel between two layers of α-helices and random coiled regions. such structural organization is commonly referred as to an adomet-dependent methyltransferase fold. tlya has an overall globular fold and comprises two domains exhibiting different spatial arrangements. the smaller domain, named s is n-terminally located and is overlapped by the catalytic domain which exhibits structural similarity to various unrelated rna-binding proteins, while the large catalytic domain, named ftsj-like methyltransferase domain is common to 2'-o-mtases and exhibits a α/β fold with a deep pocket. the s domain consists of residues, starting at a- and ending at h- this domain is structurally formed by one β-sheet and three α-helices predominantly composed of positively charged residues whose probable function is to mediate rna binding and provide stability. in line with this, the rnabindr server identified three short sequences inside the s domain that might be implicated in rna binding. the s domain has been specifically identified in some bacterial and eukaryotic ribosomal proteins, pseudouridine synthases, rna methylases, bacterial tyrosyl-trna synthetases and may also be involved in translation regulation. the analysis of the protein's electrostatic surface showed that positively charged residues were localized towards the n-terminal extreme, forming bulges and long finger-like projections that probably extend into the rrna core to stabilize its structure. in fact, charged residues like r and k exhibited the highest interface propensities, which is consistent with their ability to participate in interactions both with nucleotide bases and the negatively charged rna phosphate backbone. the ftsj-like methyltransferase domain is located between residues and and includes the catalytic tetrad responsible for tlya's enzymatic activity. moreover, the ftsj-like methyltransferase domain exhibits a rossmann fold composed of five parallel β-strand layers linked by two α-helices following an α/β secondary structure topological order. our comparative analysis of tlya with the fibrillarin family showed the presence of a gastg adomet binding motif comprising residues 90- and a d located at position juxtaposed to the above mentioned motif, which might stabilize tlya coenzyme binding. adomet was located between β-strands and within the inner depression inside the catalytic pocket formed by gastg residues' main-chain nh groups. moreover, substrate positioning might be driven by the rossmann configuration fold and the pertinent adomet-binding/catalytic sites, suggesting substrate selectivity in the small ribosomal subunit. likewise, domain structure was in agreement with s-adenosyl-methionine -dependent methyltransferases. adomet binding and affinity might indicate how an enzyme such as tlya may methylate the two adjacent ribose rings on both ribosomal subunits. in this reaction s-adenosyl-methionine acts as methyl donor, as it is converted into s-adenosyl-l-homocysteine; tlya could thus be also classified as being an rna-adomet-mtase class i. this hypothesis is supported by the high structural similarity between tlya and ftsj rna-mtase complexed with adomet. two expected structural motifs were also identified in this structure which was formed by six buried residues between g- and f- which are pivotal for adomet interaction and therefore conserved in mycobacterium tlya proteins. nest analysis of two structural motifs detected in tlya functionally important regions for adomet binding region. results were obtained on profunc server. the tlya model was validated by using the rampage server which evaluates a protein's 3d-structure based on the stereochemistry quality of torsion angles and geometry. the ramachandran plot showed % feasibility in favorable positions and % in acceptable regions, with just % in an outlier region, supporting the high quality of the 3d-model. this in silico structure analysis of m. tuberculosis proteins allowed a deeper understanding of tlya's function and challenges the initially described function of tlya-related proteins. rather than its hypothetical function as a hemolysin, our data suggested functions as an mtase based on in silico prediction. predicting subcellular localization subcellular localization was predicted with tbpred and the results suggested tlya localization in the cytoplasmatic compartment, also supporting the fact that tlya acts as a methyltransferase. a recent mycobacterial protein analysis led to ascertaining several predictors' reliability for subcellular localization, providing highly accurate results for gpos-ploc and pa-sub v. and signalp as a confirmatory tool. however, in this study, pa-sub v. predicted tlya as being an extracellular protein, which could have been biased as the training strategy used involved results from the swiss-prot database where tlya is mis-annotated. likewise, in agreement with blast scores for the query sequence, out of homologous proteins were designated as putative ribosomal rna methyltransferases and rna binding proteins. only one submitted sequence seemed to be similar to hemolysin a. nd: non defined value. pse-pssm: pseudo position-specific score matrix. svm: support vector machines. pssm: position-specific scoring matrix. meme/mast: multiple em for motif elicitation/motif alignment and search tool. hm: hybrid module. pseaac: pseudo amino acid composition. evolutionary insights into tlya genes evolutionary analysis of tlya genes from mycobacterial species, species taxonomically related to mycobacterium, adomet-rna mtases and pfp showed similar topology by both neighbor-joining and maximum parsimony methods, suggesting that the phylogenetic relationship observed was highly reliable. the consensus tree showed the tlya gene in m. tuberculosis h37rv, cdc c, haarlem, beijing, kzn and inh-rif strains as well as in the m. bovis af and bcg strains to fit in a monophyletic cluster. the tlya tree exhibited short genetic distances for m. tuberculosis complex species, including those separated according to clinical considerations. phylogenetic analysis showed that mycobacterium tlya genes are closely related to those of nocardia and rhodococcus genus; such branching patterns were supported by high bootstrap values and agreed with the results obtained from 16s rdna phylogeny. no differences were observed between antibiotic susceptible and resistant m. tuberculosis strains. mycobacterium tlya proteins fell into groups; each group formed a well-delineated branch supported by high bootstrap values. m. abscessus was the most distant and tlya-ancestral species from the mycobacterium genus compared to the closest species to m. tuberculosis. serpulina hyodysenteriae and lawsonia intracellularis lsaa antigen sequences were included due to their sequence identity to m. tuberculosis tlya, sharing a common ancestor, although having experimentally-attributed different roles. pattern branching of the phylogenetic tree showed that the tlya protein formed a monophyletic cluster with adomet rna mtases, suggesting a common evolutionary origin; these findings might also imply a similar functional role in translational machinery. pfp formed a distant cluster where diverse proteins were clustered on a branch with the most divergent proteins. the synonymous and non-synonymous substitutions were calculated to identify the action of neutral, positive or purifying selection on tlya genes; their synonymous, pi, and non-synonymous diversity, pi, was and, respectively, corresponding to a pi/pi rate which can be inferred as purifying selection acting on this gene among mycobacterium species. this suggests that conservation at sequence level might show mycobacterium tlya gene homogeneity and show low inter-specific variation between the species compared in this study. discussion the aim of the present study was to describe structural and phylogenetic features of m. tuberculosis tlya, using bioinformatics tools. the structural and phylogenetic analysis shown here highlights the need for re-annotating this protein family and redefining it as an rna-binding ftsj-like methyltransferase in bioinformatics databases; this is based on the high degree of amino acid sequence similarity found in the mycobacterium species analyzed here. it can be surmised that physical organization, domain spatial arrangement and protein folding on tlya is structurally conserved. previous sequence analysis of 2'-o-mtase families has indicated that they use a similar catalytic mechanism and have inherited a common function from the same ancestral ribose 2'-o-mtase. feder et al., suggested that tlya is a 2'-o-mtase, identifying a catalytic tetrad consisting of lysine- aspartic acid- lysine- and glutamic acid in position. based on amino acid sequence alignment, our comparative analysis with the rrmj family and other related rna 2'-o-mtases showed changes in d154-k182-e residue positions. however, despite being located in different sequence positions, the scaffold formed by these residues was maintained at the catalytic site in the modeled protein 3d structure without affecting protein function. despite the diversity of existing mtases, the majority maintain the same catalytic mechanism where k- binds to a hydrogen from the ribose 2-oh' group; this nucleophilic attack generates a transition state. the lysine side chain shifts e- so that its negatively charged carboxyl group becomes available to promote the attack of the methyl donor; furthermore, this conserved amino acid residue pattern was also observed in our evolutionary trace analysis. tlya three-dimensional structure resembled an rna methyltransferase the s. hyodysenteriae tlya postulated mechanism for hemolysin action is said to be similar to that for pore-forming protein mechanisms; tlya folding, however, is inconsistent with pfps due its lack of membrane-binding domains, cholesterol-recognition and insertion motifs, and other characteristic features such as a conserved undecapeptide near the c- or n- terminal ends. it also lacks repeat domains in the toxin and calcium-binding glycine-rich motifs necessary for oligomerization. this in silico structure analysis of m. tuberculosis proteins led to a deeper understanding of tlya's role and challenged the initially described function for tlya-related proteins. a combination of bioinformatics analysis and experimental structure elucidation may suggest alternative functional roles for previously-annotated enzymes, as with -dependent methyltransferase, phosphatase and n-acetyltransferase structures from m. tuberculosis. theoretical evidence has thus been provided for tlya's role in rna modification instead of it just being a hemolytic factor. subcellular localization suggested a cytoplasmatic localization consistent with ribosome-associated activity predicting subcellular localization is important since several characteristics can be deduced, such as protein function and genome annotation. it also aids experimental design for proteomics platforms, particularly for identifying new candidates for vaccine development and drug targets. protein sequence analysis did not reveal a signal peptide, suggesting that tlya is not secreted; however, some authors have described tlya as being a probable pfp, supported by the presence of over-expressed tlya on the e. coli wall surface. tlya lacks transmembrane helices and β-barrels, suggesting that this protein might not be a membrane-embedded or contact-dependent hemolysin, as has previously been proposed. tlya also has a hydrophobic region in the c-terminal domain which may not form a membrane-spanning α-helix. nonetheless, such a region might allow the adenine present in adomet to fit properly inside the catalytic pocket, thus suggesting a different role for these hydrophobic residues. virulence-related functional annotation, based on all the above-mentioned data, remains controversial. conservation of tlya genes and functions was expected according to our phylogenetic and evolutionary trace analyses, the initial hypothesis being that homologous proteins have related functions. however, the extent to which this is true has not been assessed in detail up to now; annotation based on protein homology might thus lead to finding unexpected variations in function.  the main concern of our work was to decipher structural information about the tlya family by homology comparison. our structural analysis and finding residue conservation within the active site supports recent experimental work about tlya's role as rna 2`-o-mtase in m. tuberculosis. this is further supported by the domain composition also suggesting methyltransferase activity and rna-binding. the distribution of charged residues also defines a likely rna interaction instead of hemolytic activity, little similarity being exhibited with hemolysin proteins, perhaps representing inaccurate annotations in several databases, as observed with other protein families. since genome annotation is based on sequence comparison, it is the most commonly used approach for determining functional homology. annotations are assigned if proteins surpass similarity cut-off and the propagation of original errors could be increased by each new entry. in such cases, protein function needs to be inferred from their common three-dimensional structures and manual curation. 
at present, its dual activity as hemolysin and rna methyltransferase in m. tuberculosis has been indirectly proposed based on in vitro results. three-dimensional modeling of tlya showed a consensus structure, having a common core formed by a six-stranded β-sheet between two α-helix layers, which is consistent with an rna methyltransferase structure. the tlya protein has a controversial function as a virulence factor in mycobacterium tuberculosis. altogether, our manual in silico curation suggested that tlya is involved in ribosomal biogenesis and that there is a functional annotation error regarding this protein family in several microbial and plant genomes, including the m. tuberculosis genome. additionally, the nucleotide substitution rates suggested purifying selection during tlya gene evolution and the absence of a common ancestor between tlya proteins and bacterial pore-forming proteins. phylogenetic analyses showed high conservation of the tlya gene among mycobacterium species. subcellular localization prediction showed that tlya lacks a signal peptide and its hydrophobicity profile showed no evidence of transmembrane helices. bioinformatics analysis of tlya identified a ribosomal protein binding domain, located between residues and as well as an ftsj-like methyltranferase domain encompassing residues and all of which have been previously described in translation machinery-associated proteins.
 the tlya protein has a controversial function as a virulence factor in mycobacterium tuberculosis. at present, its dual activity as hemolysin and rna methyltransferase in m. tuberculosis has been indirectly proposed based on in vitro results. there is no evidence however for tlya relevance in the survival of tubercle bacilli inside host cells or whether both activities are functionally linked. a thorough analysis of structure prediction for this mycobacterial protein in this study shows the need for reevaluating tlya's function in virulence.  bioinformatics analysis of tlya identified a ribosomal protein binding domain, located between residues and as well as an ftsj-like methyltranferase domain encompassing residues and all of which have been previously described in translation machinery-associated proteins. subcellular localization prediction showed that tlya lacks a signal peptide and its hydrophobicity profile showed no evidence of transmembrane helices. these findings suggested that it may not be attached to the membrane, which is consistent with a cytoplasmic localization. three-dimensional modeling of tlya showed a consensus structure, having a common core formed by a six-stranded β-sheet between two α-helix layers, which is consistent with an rna methyltransferase structure. phylogenetic analyses showed high conservation of the tlya gene among mycobacterium species. additionally, the nucleotide substitution rates suggested purifying selection during tlya gene evolution and the absence of a common ancestor between tlya proteins and bacterial pore-forming proteins.  altogether, our manual in silico curation suggested that tlya is involved in ribosomal biogenesis and that there is a functional annotation error regarding this protein family in several microbial and plant genomes, including the m. tuberculosis genome. 