- Research article
- Open Access
Functional characterization of two members of histidine phosphatase superfamily in Mycobacterium tuberculosis
© Coker et al.; licensee BioMed Central Ltd. 2013
- Received: 5 July 2013
- Accepted: 7 December 2013
- Published: 11 December 2013
Functional characterization of genes in important pathogenic bacteria such as Mycobacterium tuberculosis is imperative. Rv2135c, which was originally annotated as conserved hypothetical, has been found to be associated with membrane protein fractions of H37Rv strain. The gene appears to contain histidine phosphatase motif common to both cofactor-dependent phosphoglycerate mutases and acid phosphatases in the histidine phosphatase superfamily. The functions of many of the members of this superfamily are annotated based only on similarity to known proteins using automatic annotation systems, which can be erroneous. In addition, the motif at the N-terminal of Rv2135c is ‘RHA’ unlike ‘RHG’ found in most members of histidine phosphatase superfamily. These necessitate the need for its experimental characterization. The crystal structure of Rv0489, another member of the histidine phosphatase superfamily in M. tuberculosis, has been previously reported. However, its biochemical characteristics remain unknown. In this study, Rv2135c and Rv0489 from M. tuberculosis were cloned and expressed in Escherichia coli with 6 histidine residues tagged at the C terminal.
Characterization of the purified recombinant proteins revealed that Rv0489 possesses phosphoglycerate mutase activity while Rv2135c does not. However Rv2135c has an acid phosphatase activity with optimal pH of 5.8. Kinetic parameters of Rv2135c and Rv0489 are studied, confirming that Rv0489 is a cofactor dependent phosphoglycerate mutase of M. tuberculosis. Additional characterization showed that Rv2135c exists as a tetramer while Rv0489 as a dimer in solution.
Most of the proteins orthologous to Rv2135c in other bacteria are annotated as phosphoglycerate mutases or hypothetical proteins. It is possible that they are actually phosphatases. Experimental characterization of a sufficiently large number of bacterial histidine phosphatases will increase the accuracy of the automatic annotation systems towards a better understanding of this important group of enzymes.
- Acid phosphatase
- Hypothetical protein
- Histidine phosphatase superfamily
- Phosphoglycerate mutase
- Mycobacterium tuberculosis
Mycobacterium tuberculosis remains a threat to global health despite efforts directed towards its eradication. Although several works have been done in recent years towards understanding the genetic repertoire of this organism, many of its strategies involved in virulence, pathogenesis and resistance to both host pressure and antibiotics remain elusive . Mycobacterial genome has been completely sequenced for over a decade . However, the functions of many of its genes are annotated based only on similarity to known proteins using automatic annotation systems. This method of function annotation can be erroneous [3, 4]. Errors in automatic function annotation to genes in bacterial genomes are well documented. They often lead to misinformation that may hamper the understanding of the roles played by many bacterial genes [5–8]. Experimental characterization of additional mycobacterial proteins is needed to aid deeper understanding of the organism.
Histidine phosphatase superfamily is a large family of proteins with diverse functions that are important. This superfamily comprises two branches. The larger branch consists of proteins which function in metabolic regulations, intermediary metabolism and developmental processes. Examples include cofactor dependent phosphoglycerate mutases, alpha-ribazole phosphatase (CobC), mannitol-1-phosphatase, fructose-2,6-bisphosphatase and acid phosphatase (PhoE) [9–11]. The smaller branch consists mainly of phosphatases and phytases with functions ranging from extracellular metabolism to involvement in developmental processes [9, 12]. Examples include human testicular acid phosphatase and lysosomal acid phosphatase [9, 13, 14]. The functions of enzymes in this superfamily are based on a conserved catalytic histidine residue in the motif ‘RHG’ present at the N terminal, which becomes phosphorylated during the reaction [9, 15].
Members of the histidine phosphatase superfamily that have been studied in M. tuberculosis, include Rv0489. The crystal structure of Rv0489 at 1.7 Å resolution reveals the catalytic residues superimposing with those of the cofactor dependent phosphoglycerate mutase of E. coli, with which it shares 42% amino acid identity . However, its biochemical characteristics remain unknown. Other members include Rv3214c, an acid phosphatase with unknown specific substrate  and Rv2419c which was characterized as glucosyl-3-phosphoglycerate phosphatase in lipopolysaccharide biosynthesis with an optimum pH of 7.0 . Rv2135c is a paralog of the aforementioned members of the superfamily, but it is annotated as a hypothetical protein in the genomic database of M. tuberculosis. Bioinformatics similarity searches show that it is a probable cofactor dependent phosphoglycerate mutase. However, there have been reports that proteins annotated as cofactor dependent phosphoglycerate mutases by sequence similarity actually perform the functions of an acid phosphatase when assayed in vitro. Examples in M. tuberculosis are Rv2419c  and Rv3214c . In other organisms, examples include PhoE of Bacillus stearothermophillus, and PfPGM2 of Plasmodium falciparum[4, 19].
Rv2135c was found in Triton X-114 fractions of M. tuberculosis H37Rv strain and reported as one of the cell envelope associated hypothetical proteins . Rv2135c contains a catalytic histidine motif similar to proteins in histidine phosphatase superfamily. Nevertheless, its motif is ‘RHA’ unlike ‘RHG’ commonly found in histidine phosphatase superfamily. These motivate the need to investigate its function in the metabolism of M. tuberculosis.
Phosphoglycerate mutases (EC 220.127.116.11) primarily interconvert 3-phosphoglyceric acid (3-PGA) and 2-phosphoglyceric acid (2-PGA) in both glycolysis and gluconeogenesis [12, 21]. Two different types of phosphoglycerate mutase have been identified. One depends on the cofactor, 2,3-bisphosphoglyceric acid, for activity (dPGMs) while the other does not (iPGMs) [12, 21]. The cofactor-dependent form is found in vertebrates, budding yeast, and bacterial species, while the cofactor-independent form is the only phosphoglycerate mutase present in higher plants. Some bacteria like E. coli, however, possess both forms . There is no amino acid sequence similarity between these two types of PGMs and their structures are also quite different. Deficiencies in dPGM in E. coli and yeast have been linked to severely impaired growth [23, 24].
Acid phosphatases (EC 18.104.22.168) catalyze the hydrolysis of phosphate monoesters or transfer of phosphate groups between phosphoester and alcohols. The enzymes catalyze optimally at acidic conditions and are completely and structurally different from alkaline phosphatases (EC 22.214.171.124), which work optimally at alkaline conditions [25–27]. Unlike the alkaline phosphatases, the acid phosphatases, do not utilize metal ions in their catalysis. They rather utilize histidine residue to form a phospho-histidine-enzyme intermediate which is essential for their catalysis. In contrast, alkaline phosphatases make use of a phospho-serine-enzyme intermediate for their catalysis and have a binuclear Zn (II) active site [26, 28]. Phosphatases are important in the physiology of an organism as they function in many catalytic reactions relating to activation or deactivation of enzymes. Deficiencies in phosphate metabolism have been reported to be related to reduction of virulence in many bacterial species such as Listeria monocytogenes, Streptococcus pneumoniae, Vibrio cholerae, Proteus mirabilis and M. tuberculosis[29–34].
The fact that histidine acid phosphatases and cofactor dependent phosphoglycerate mutases share similar catalytic amino acid residues and mechanism of catalysis warrants their placement in the same superfamily . This often leads to some difficulties in predicting the function of an enzyme that belongs to the superfamily. Thus, biochemical characterization of purified enzymes is necessary before the function of any member of histidine phosphatase superfamily can be ascertained. In this study, we report the first cloning, purification and characterization of M. tuberculosis Rv2135c. In addition, we cloned and characterized Rv0489. Its role as a cofactor dependent phosphoglycerate mutase was confirmed.
The histidine phosphatase motif in Rv2135c
Cloning and expression of C-His-Rv2135c and C-His-Rv0489
The purification table of C-His-Rv0489c from 1 liter culture of E. coli BL21(DE3)-89c
Specific activity (U/mg)
Eluted fractions from IMAC
The purification table of C-His-Rv2135c from 1 liter culture of E. coli BL21(DE3)-35c
Specific activity (U/mg)
Eluted fractions from IMAC
Enzymatic activities of C-His-Rv2135c and C-His-Rv0489
Kinetic parameters for the phosphoglycerate mutase activity of C-His-Rv0489
kcat/Km (mM-1 min-1)
0.40 ± 0.02
250460 ± 8100
626100 ± 20300
Kinetic parameters for the acid phosphatase activity of C-His-Rv2135c at pH 5.8 using pNPP as substrate
kcat/Km (mM-1 min-1)
10.60 ± 0.07
4170 ± 100
392 ± 10
Substrates for C-His-Rv2135c
Kinetic parameters for the activities of C-His-Rv2135c with different substrates at pH 5.8
Specific activity (mol/min/mg)
0.23 ± 0.07
10.60 ± 0.07
0.09 ± 0.002
11.25 ± 0.75
0.05 ± 0.002
14.00 ± 0.00
0.08 ± 0.009
7.75 ± 0.75
Native molecular mass and stability
The size of the native form of C-His-Rv2135c was estimated by gel filtration to be 104.70 kDa. With the amino acid calculated size of 25.95 kDa, this suggests that C-His-Rv2135c forms a tetramer in the native state. This conforms to the results obtained by ND-PAGE, which provided the estimated native size of 103.85 kDa. The molecular mass of the native form of C-His-Rv0489 estimated from the gel filtration is 56.02 kDa. This indicates that C-His-Rv0489 forms a dimer, given both calculated and SDS-PAGE estimated molecular mass of the monomer of 28 kDa.
The acid phosphatase activity of C-His-Rv2135c at pH 5.8 was found to be enhanced by 15% in the presence of 10 mM magnesium ion. The enzyme was found to be stable in 50% glycerol at −20°C for up to 4 months with no significant change in activity.
In addition to Rv2419c  and Rv3214  characterized recently, we have presented the study of a new mycobacterial phosphatase belonging to the histidine phosphatase superfamily. We report the first cloning, expression and characterization of Rv2135c, annotated as hypothetical in the genome database of M. tuberculosis. Simple NCBI BLAST [35, 38] reveals that most of the proteins similar to Rv2135c are annotated as hypothetical proteins or phosphoglycerate mutases. We demonstrated that C-His-Rv2135c possesses neither phosphoglycerate mutase nor phosphoglycerate phosphatase activity. However, it has phosphatase activity in acidic condition. Our findings support the necessity to experimentally characterize enzymes before their biochemical functions can be ascertained. This is important especially for the histidine phosphatase superfamily whose members can perform different metabolic functions [3, 4, 9, 19].
C-His-Rv2135c has 6 more histidine residues at the C- terminal region than the native protein. The method of C-terminal tagging is commonly used for facilitating purification of enzymes and generally does not affect enzyme specificities. The specific acid phosphatase activity of C-His-Rv2135c (0.23 μmol/min/mg) is about 10 times less than that of Rv3214 (2.6 μmol/min/mg). However, some acid phosphatases of other pathogenic microorganisms are known to possess less specific activities than that of C-His-Rv2135c. Examples include the phosphatases of Francisella tularensis with specific activity of 0.002 μmol/min/mg  and Entamoeba histolytica with specific activity of 0.012 μmol/min/mg . It should also be noted that the histidine phosphatase superfamily typically contains the characteristic motif ‘RHG’ at the N-terminal region. However, the motif present in Rv2135c is ‘RHA’ as found in the yet uncharacterized phosphoglycerate domain containing protein of C. parvum (GAN CAD98474). The replacement of glycine with alanine, another non-polar amino acid with a small side chain, may occur without any effect on the specificity of the enzymes in this family. Moreover, Rv2135c contains other residues reported to be important in the phosphatase activities of other members of the superfamily. These include Arg57, Glu82, and a fully conserved His153 at the C-terminal region [3, 9, 36]. Thus, we believe that Rv2135c performs an acid phosphatase function in its native environment.
The substrate specific to Rv2135c is unknown. Its sequence appeared to have little similarity to other previously annotated histidine phosphatases of M. tuberculosis, although the annotations of most of these phosphatases are still computational. Therefore there is no information suggesting the primary substrate of the enzyme. There are few experimentally characterized phosphatases in M. tuberculosis. These include Rv3214 and Rv2419c, which are histidine phosphatases [3, 17], PtpA and PtpB which are tyrosine protein phosphatases [41, 42], and PstP, a serine/threonine protein phosphatases . The specific substrates of these phosphatases have not been identified yet, with the exception of Rv2419c, a glucosyl-3-phosphoglycerate phosphatase .
There are several known functions of histidine acid phosphatases, including extracellular metabolism, scavenging and regulatory functions. Rv2135c was identified as being associated with membrane protein fractions [20, 44]. M. tuberculosis encounters a phosphate deficient acidic environment in an infected macrophage, and has been shown to depend on the acquisition of phosphate groups from the host environment for survival . It is therefore intriguing to further study whether Rv2135c plays some roles in the intramacrophage environment, where it has been shown to be expressed .
Rv2135c and Rv2136c have been predicted to be in the same operon (http://genome.tbdb.org/annotation/genome/tbdb/). Rv2136c is the only mycobacterial gene with the catalytic motif of undecaprenyl pyrophosphate phosphatase. In bacteria, the enzyme hydrolyzes undecaprenyl pyrophosphate to produce undecaprenyl phosphate needed to translocate various cell wall intermediates from the cytosol across the cytoplasmic membrane for polymerization [46, 47]. Despite the apparent essentiality of this function, undecaprenyl pyrophosphatases of many bacteria are known to be non-essential for their growth [48, 49]. Rv2136c has also been shown to be non-essential for the survival of M. tuberculosis. In some bacteria such as E. coli, other membrane associated phosphatases were shown to successfully perform the function of undecaprenyl pyrophosphatase when it was knocked out [49, 51–53]. These include PgpB, YbjG and YeiU of E. coli, which belong to type 2 phosphatidic acid phosphatase family . As Rv2135c and Rv2136c are predicted to be in the same operon, it may be possible that membrane associated Rv2135c performs a role similar to Rv2136c.
According to String Prokaryotic Operon Predictor (http://operons.ibt.unam.mx/OperonPredictor/), homologs of Rv2135c are identified in the same operon as the homologs of Rv2136c (undecaprenyl pyrophosphate phosphatase gene) in some other mycobacteria. These include M. marinum, M. ulcerans, M. smegmatis and M. leprae, but not M. avium. Using tblastx [35, 38], it was found that homologs of Rv2135c and Rv2136c share adjacent positions in the genome of a number of other bacteria belonging to the actinomycetales such as Nocardioides, Micrococcus, Cellulomonas, Geodermatophilus, etc. Additional experiments are needed to investigate the functional relationship between these two genes.
Using Phyre2 , Rv2135c was modeled as a globular protein with a fairly large and hydrophobic pocket on its surface, which might provide a binding space for an undecaprenyl (see Additional file 2). A novel type of phosphoserine phosphatase of Hydrogenobacter thermophiles was also identified as the most similar protein with known crystallographic structural data. However, the possible tetrameric structure of Rv2135c in the native form warrants further biochemical, computational and crystallographic studies in order to ascertain the natural substrate of this enzyme.
The crystal structure of Rv0489 was previously determined at 1.7 Å resolution. The residues at its active site were demonstrated to superimpose with corresponding residues of E. coli cofactor dependent phosphoglycerate mutase . This study presents the first report of its biochemical activity and kinetic parameters, confirming it as a mycobacterial cofactor dependent phosphoglycerate mutase. Rv0489 was earlier found to be essential for the in vitro growth of H37Rv strain of M. tuberculosis by Himar1-based transposon mutagenesis , making it a putative target for drug development. Information about its kinetic parameters may be useful for formulating target-based screening assay for new drug discovery. This study shows that Rv0489 forms a dimer in solution. However, previous crystallization study carried out on Rv0489 showed it as a tetramer and referred to it as a dimer of dimers . Cofactor dependent phosphoglycerate mutases from E. coli and Homo sapiens have been shown to be dimers [57, 58] while those from Saccharomyces cerevisae and Lactococcus lactis are tetramers [59, 60].
Most well-characterized histidine acid phosphatases were reported from eukaryotes . A bacterial histidine phosphatase is usually labeled as a phosphoglycerate mutase by automatic annotation systems. This is partly due to the much more abundance of bacterial phosphoglycerate mutases than bacterial histidine acid phosphatases in the databases. Orthologous proteins to Rv2135c, identified by reciprocal BLAST, are found widely in other mycobacteria as well as various taxa of bacteria, including Staphylococcus aureus and E. coli. Most of them are annotated as phosphoglycerate mutases or hypothetical proteins. It is possible that they are actually phosphatases. Experimental characterization of a sufficiently large number of bacterial histidine phosphatases will increase the accuracy of the automatic annotation systems towards a better understanding of this important group of enzymes.
Bacteria strains and culture conditions
E. coli strain DH5α was used for the maintenance and cloning of plasmids. Plasmid pET23b (Novagen, USA) was used as expression vector. It contains an inbuilt optional C-terminal hexahistidine tag for ease of protein purification. E. coli BL21 (DE3) was used as recipient hosts for recombinant protein expression . E. coli was grown in Luria-Bertani (LB) medium. M. tuberculosis H37Ra (ATCC 25177) was grown on Middlebrook 7H11 agar supplemented with 10% Middlebrook OADC [Oleic acid Albumin Dextrose Catalase] Enrichment (Difco BBL, USA). M. tuberculosis genomic DNA was prepared as previously described .
Identification of histidine phosphatase motif in Rv2135c
Using NCBI BLAST [35, 38], Rv2135c protein was found to be similar to proteins of histidine phosphatase superfamily. Some of the similar proteins were aligned with Rv2135c using ClustalX2 with the default parameters . The similar proteins included in the alignment are some experimentally characterized and predicted members of the superfamily. These are M. tuberculosis probable co-factor dependent phosphoglycerate mutase Rv0489 (GenBank accession number (GAN) CAE55288.1) , E. coli cofactor dependent phosphoglycerate mutase (E.colidpgM, Swissprot P62707), PhoE a broad specificity phosphatase from B. stearothermophilus (Protein data bank (PDB)1H2E_A) , Rv3214, (GAN CAE55568) a M. tuberculosis acid phosphatase , an acid phosphatase from Bacillus licheniformis (Bacillusap, GAN EID46354), newly characterized glucosyl-3-phosphoglycerate phosphatase of M. tuberculosis, Rv2419c  (Swissprot P71724), and Rv3837c (GAN CAB06204) an uncharacterized paralog of Rv2135c. Members of histidine phosphatase superfamily from eukaryotes, the cofactor dependent phosphoglycerate mutase of Saccharomyces arboricola (YDR051pgm) (GAN EJS44264) and phosphoglycerate mutase domain containing protein of Cryptosporidium parvum (Cryparpgm) (GAN CAD98474) were also included.
Cloning of Rv2135c and Rv0489
The open reading frame of Rv2135c and Rv0489 in the virulent strain H37Rv of M. tuberculosis is completely identical to the non-virulent strain H37Ra. The genomic DNA of H37Ra was, therefore, used as the template for the amplification of Rv2135c and Rv0489 by polymerase chain reaction (PCR) using primers 2135EF (AGCCGCCATATGACCGTCATCTTGCTACG) and 2135CER(ACTCGAGGTCGGTGGAACCGCCGATC), and primers Rv0489-F (CTTGCCCATATGGCAAACACTGGCAGCCTGG) and Rv0489-R (CTCAAACTCGAGCCCGCGGCCCTGGCCGGCCA) respectively. Recognition sequences of the restriction enzymes NdeI (CATATG) and XhoI (CTCGAG) are in italics. The PCR products were digested with NdeI and XhoI, ligated to appropriately digested expression plasmid pET23b with C-terminal histidine tag, and transformed into E. coli DH5α. Transformants were selected with ampicillin. Plasmids were purified from the transformants and sequenced to confirm the presence of correct genes tagged at the 3′ end with 6 histidine codons, designated as pET23b-35c and pET23b-89 respectively. E. coli BL21(DE3) strain was transformed with pET23b-35c and pET23b-89 to obtain strains BL21(DE3)-35c and BL21(DE3)-89 used for protein expression and from which the recombinant C-terminal histidine tagged proteins C-His-Rv2135c and C-His-Rv0489 were purified respectively.
Expression and purification
A liter of LB medium with 100 μg/ml ampicillin was inoculated with an overnight culture of BL21(DE3)-35c to a final OD600nm of about 0.03. The culture was incubated at 37°C with shaking speed of 200 rpm until OD600nm reached about 0.6. The expression of the protein was then induced by the addition of IPTG to a final concentration of 0.4 mM. The culture was further incubated at 25°C at the shaking speed of 200 rpm for 8 hours. Cells were harvested by centrifugation at 3500 rpm at 4°C, washed with PBS pH 7.4 and stored at −20°C.
Similar treatment of BL21(DE3)-89 was done and resulted in precipitation of expressed protein after lysis. In order to obtain C-His-Rv0489 in the soluble fraction of cell lysate, BL21(DE3)-89 was cultured in the same media as above with the addition of 10% sucrose to OD600nm of about 0.03. After the OD600nm reached about 0.6, the expression of C-His-Rv0489 was induced with 0.03 mM of IPTG at 18°C overnight. Cells were harvested by centrifugation at 3500 rpm at 4°C, washed with PBS pH 7.4 and stored at −20°C.
Frozen cells were thawed on ice and suspended in the lysis buffer (20 mM Tris–HCl pH 7.0, 100 mM NaCl, 1 mM PMSF, 5 mM imidazole). The suspended cells were lysed by sonication using Misonix Sonicator 3000 (Qsonica LLC, USA) with 30 sec pause intervals until a clear lysate was obtained. The lysate was centrifuged at 11,000 rpm at 4°C for 20 min. The supernates, which contained the expressed histidine tagged protein, C-His-Rv2135c and C-His-Rv0489, were separated from other soluble proteins by immobilized metal affinity chromatography (IMAC). Briefly, the crude extracts were applied to cobalt charged resin column (Talon® Superflow column, GE Healthcare, Sweden) pre-equilibrated with the wash buffer (20 mM Tris–HCl pH 7.0, 100 mM NaCl, 5 mM imidazole). The column was then washed with 4 volumes of the wash buffer. For C-His-Rv2135c, the progress of purification was monitored by fast protein liquid chromatography (FPLC) using AKTA system (GE Healthcare, Sweden). Elution of the bound protein was done using 8 volumes of elution buffer I (20 mM Tris–HCl pH 7.0, 100 mM NaCl, containing a gradient of 0–60 mM imidazole). Eluted fractions were collected and loaded on SDS-PAGE to determine the purity of eluted proteins.
For C-His-Rv0489, after washing with 4 column volumes of lysis buffer, elution was done with elution buffer II (20 mM Tris–HCl pH 7.0, 100 mM NaCl, 150 mM of imidazole). The fractions with highest amount of recombinant C-His-Rv0489, determined by SDS PAGE were pooled and diluted to the imidazole concentration of 15 mM. The pooled fractions were then applied a second time to the cobalt charged resin column pre-equilibrated with wash buffer. The process of purification was repeated as the first column application to obtain pure C-His-Rv0489.
Purified C-His-Rv2135c and C-His-Rv0489 were concentrated using Amicon–Ultra 4 centrifugal filter unit (Merck Millipore USA) and stored in 20 mM Tris–HCl pH 7.0 containing 50% glycerol.
Phosphoglycerate mutase activity: Phosphoglycerate mutase activities of C-His-Rv2135c and C-HisRv0489 in the 3-PGA to 2-PGA (forward) direction were monitored using an assay coupled to the oxidation of NADH as earlier described . The assay was done in 500 μl of reaction mixture, containing 30 mM Tris–HCl pH 7.0, 20 mM KCl, 5 mM MgSO4, 1 mM ADP, 0.15 mM NADH, 0.2 mM 2,3-bisphophoglyceric acid, 2.5 U enolase (Sigma), 2.5 U pyruvate kinase (Sigma), 2.5 U lactate dehydrogenase (Sigma)  with ten concentrations of 3-phosphoglyceric acid (Sigma) (0.019, 0.039, 0.078, 0.156, 0.312, 0.625, 1.25, 2.5, 5 and 10 mM). Changes in absorbance at 340 nm using spectrophotometer (Thermo Electron Corporation, USA) were used in monitoring the oxidation of NADH. The values of absorbance of test solutions were corrected by the absorbance of the solution without enzymes. The assays were carried out in triplicate.
Acid phosphatase assay: The phosphatase activity was measured by monitoring the release of p-nitrophenol from p-nitrophenyl phosphate (pNPP) at a range of pH (3.0-7.5) as earlier described . 25 mM sodium citrate buffer was used at pH 3.0-6.2 while 25 mM Tris–HCl was used at pH 7.0 and 7.5. The reaction, carried out at 37°C was started by the addition of the enzymes to the pre-warmed reaction buffer with eight concentrations of pNPP (New England Biolabs, USA) (0.78, 1.56, 3.125, 6.25, 12.5, 25, 50 and 100 Mm) in a total volume of 200 μl. The mixture was incubated for 60 min, and stopped with the addition of 600 μl of 1 N NaOH. Potato acid phosphatase (Sigma) was used as a positive control at pH 4.8 with 25 mM sodium citrate buffer. The amounts of released p-nitrophenol were estimated from the change in absorbance at 405 nm, corrected by the absorbance of the solution without the enzymes incubated at 37°C for the same period of time. All assays were carried out in triplicate.
Malachite green assay: The activities of C-His-Rv2135c with other substrates were investigated. The amount of phosphate group released from each of the substrates was detected by using the malachite green phosphate detection system (R&D Systems). It is based on quantification of the green complex formed between malachite green, molybdate and free orthophosphate as earlier described . Phosphatase reaction was carried out in 25 mM sodium citrate buffer pH 5.8 at 37°C for 60 min in the presence of eight concentrations (0.78, 1.56, 3.125, 6.25, 12.5, 25, 50 and 100 mM) of glycerol-1-phosphate, glucose-6-phosphate, fructose-6-phosphate, adenosine diphosphate (ADP), phosphoenolpyruvate and 3-phosphoglyceric acid. The detection system was used according to the manufacturer’s instruction to detect the amount of released orthophosphate. The rapid color formation from the reaction was measured by the change in absorbance at 600 nm using a microplate reader (Glomax Multi Detection System, Promega, USA). The amounts of orthophosphate hydrolyzed were estimated in relation to a standard curve constructed with phosphate standard, according to the manufacturer’s instruction. All absorbance results were corrected for enzyme-unrelated absorbance change and all assays were carried out in triplicate.
Estimation of the kinetic parameters: The rate constants (Km) were estimated using Michaelis-Menten kinetics by plotting the values of reaction rates obtained against the concentrations of substrates. The curves were fit non-linearly by generalized reduced gradient (GRG) solving method using the Solver add-in in Microsoft Excel. Km was determined for each experiment and averaged. The specific activities, turnover numbers (kcat) and the catalytic efficiencies (kcat/Km) were estimated using Michaelis-Menten kinetics.
Determination of molecular mass
The native molecular mass of C-His-Rv2135c was determined under non-denaturing condition by gel filtration chromatography and native polyacrylamide gel electrophoresis (ND-PAGE) while gel filtration only was used for the determination of the molecular mass of C-His-Rv0489 in solution.
Pre-packed 10 mm X 30 cm column of Superdex 200 HR 10/30 equilibrated in 20 mM sodium phosphate buffer, pH 7.0, containing 0.1 M NaCl was used with four standard protein markers: catalase (232 kDa), lactate dehydrogenase (140 kDa), bovine serum albumin (66 kDa) from Sigma and MPT83 (50 kDa) , a mycobacterial protein purified in our laboratory. Proteins were eluted at the buffer flow rate of 0.2 ml/min. The void volume of the column was determined by loading blue dextran unto the column. A standard curve was constructed by plotting the molecular masses versus the ratio Ve/Vo for the standard protein markers, while Ve is the volume of elution of each protein and Vo is the void volume of the column. The Ve/Vo for C-His-Rv2135c and C-His-Rv0489 were used in determining their molecular weight from the standard curve.
ND-PAGE was done as previously described . Briefly, 4 μg of the purified protein and standard proteins were loaded on four native gels of different acrylamide concentrations. The concentrations of the acrylamide in resolving gels were 6, 7, 8 and 9%. 3% acrylamide was used for stacking in each resolving gel. The relative migrations of the purified protein and the standard proteins in each gel, designated as Rf, were estimated from each gel  by dividing the migration distance of the protein standards by the migration distance of the dye front. 100log(RfX100) values for each protein standard and C-His-Rv2135c were plotted against the gel concentrations. The negative slope obtained for the standard protein was plotted against their molecular weight values to obtain a standard curve. The molecular weight of C-His-Rv2135c was estimated from the standard curve.
This work was supported by the CPMO (P-10-10647 and P-00-20209), National Science and Technology Development Agency (NSTDA), Thailand and Center for Emerging Bacterial Infections (EBI), Faculty of Science, Mahidol University. We thank Dr. Pimchai Chaiyen, Dr. Danaya Pakotiprapha, Dr. Nat Smittipat and Mr. Tada Juthayothin for their technical assistance. We also thank Dr. Daniel Anderson of UCLA-DOE Institute for Genomics & Proteomics, USA for his support.
- Santos LG, Pires GN, Azeredo Bittencourt LR, Tufik S, Andersen ML: Chronobiology: relevance for tuberculosis. Tuberculosis (Edinb). 2012, 92 (4): 293-300. 10.1016/j.tube.2012.03.006.View ArticleGoogle Scholar
- Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, Gordon SV, Eiglmeier K, Gas S, Barry CE, et al: Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998, 393 (6685): 537-544. 10.1038/31159.PubMedView ArticleGoogle Scholar
- Watkins HA, Baker EN: Structural and functional analysis of Rv3214 from Mycobacterium tuberculosis, a protein with conflicting functional annotations, leads to its characterization as a phosphatase. Journal of bacteriology. 2006, 188 (10): 3589-3599. 10.1128/JB.188.10.3589-3599.2006.PubMedPubMed CentralView ArticleGoogle Scholar
- Hills T, Srivastava A, Ayi K, Wernimont AK, Kain K, Waters AP, Hui R, Pizarro JC: Characterization of a new phosphatase from Plasmodium. Mol Biochem Parasitol. 2011, 179 (2): 69-79. 10.1016/j.molbiopara.2011.06.001.PubMedView ArticleGoogle Scholar
- Richardson EJ, Watson M: The automatic annotation of bacterial genomes. Briefings in bioinformatics. 2013, 14 (1): 1-12. 10.1093/bib/bbs007.PubMedPubMed CentralView ArticleGoogle Scholar
- Gilks WR, Audit B, de Angelis D, Tsoka S, Ouzounis CA: Percolation of annotation errors through hierarchically structured protein sequence databases. Mathematical biosciences. 2005, 193 (2): 223-234. 10.1016/j.mbs.2004.08.001.PubMedView ArticleGoogle Scholar
- Poptsova MS, Gogarten JP: Using comparative genome analysis to identify problems in annotated microbial genomes. Microbiology. 2010, 156 (Pt 7): 1909-1917.PubMedView ArticleGoogle Scholar
- Friedberg I: Automated protein function prediction–the genomic challenge. Briefings in bioinformatics. 2006, 7 (3): 225-242. 10.1093/bib/bbl004.PubMedView ArticleGoogle Scholar
- Rigden DJ: The histidine phosphatase superfamily: Structure and function. Biochem J. 2008, 409 (2): 333-348. 10.1042/BJ20071097.PubMedView ArticleGoogle Scholar
- Pilkis SJ, Lively MO, El-Maghrabi MR: Active site sequence of hepatic fructose-2,6-bisphosphatase. Homology in primary structure with phosphoglycerate mutase. The Journal of biological chemistry. 1987, 262 (26): 12672-12675.PubMedGoogle Scholar
- Fothergill LA, Harkins RN: The amino acid sequence of yeast phosphoglycerate mutase. Proc R Soc Lond B Biol Sci. 1982, 215 (1198): 19-44. 10.1098/rspb.1982.0026.PubMedView ArticleGoogle Scholar
- Fothergill-Gilmore LA, Watson HC: The phosphoglycerate mutases. Adv Enzymol Relat Areas Mol Biol. 1989, 62: 227-313.PubMedGoogle Scholar
- Fleisig H, El-Din El-Husseini A, Vincent SR: Regulation of ErbB4 phosphorylation and cleavage by a novel histidine acid phosphatase. Neuroscience. 2004, 127 (1): 91-100. 10.1016/j.neuroscience.2004.04.060.PubMedView ArticleGoogle Scholar
- Suter A, Everts V, Boyde A, Jones SJ, Lullmann-Rauch R, Hartmann D, Hayman AR, Cox TM, Evans MJ, Meister T, et al: Overlapping functions of lysosomal acid phosphatase (LAP) and tartrate-resistant acid phosphatase (Acp5) revealed by doubly deficient mice. Development. 2001, 128 (23): 4899-4910.PubMedGoogle Scholar
- Bazan JF, Fletterick RJ, Pilkis SJ: Evolution of a bifunctional enzyme: 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase. Proc Natl Acad Sci U S A. 1989, 86 (24): 9642-9646. 10.1073/pnas.86.24.9642.PubMedPubMed CentralView ArticleGoogle Scholar
- Muller P, Sawaya MR, Pashkov I, Chan S, Nguyen C, Wu Y, Perry LJ, Eisenberg D: The 1.70 angstroms X-ray crystal structure of Mycobacterium tuberculosis phosphoglycerate mutase. Acta Crystallogr D Biol Crystallogr. 2005, 61 (Pt 3): 309-315.PubMedView ArticleGoogle Scholar
- Mendes V, Maranha A, Alarico S, da Costa MS, Empadinhas N: Mycobacterium tuberculosis Rv2419c, the missing glucosyl-3-phosphoglycerate phosphatase for the second step in methylglucose lipopolysaccharide biosynthesis. Sci Rep. 2011, 1: 177-PubMedPubMed CentralView ArticleGoogle Scholar
- Lew JM, Kapopoulou A, Jones LM, Cole ST: TubercuList–10 years after. Tuberculosis (Edinb). 2010, 91 (1): 1-7.View ArticleGoogle Scholar
- Rigden DJ, Bagyan I, Lamani E, Setlow P, Jedrzejas MJ: A cofactor-dependent phosphoglycerate mutase homolog from Bacillus stearothermophilus is actually a broad specificity phosphatase. Protein Sci. 2001, 10 (9): 1835-1846. 10.1110/ps.15701.PubMedPubMed CentralView ArticleGoogle Scholar
- Malen H, Pathak S, Softeland T, de Souza GA, Wiker HG: Definition of novel cell envelope associated proteins in Triton X-114 extracts of Mycobacterium tuberculosis H37Rv. BMC Microbiol. 2010, 10: 132-10.1186/1471-2180-10-132.PubMedPubMed CentralView ArticleGoogle Scholar
- Jedrzejas MJ: Structure, function, and evolution of phosphoglycerate mutases: Comparison with fructose-2,6-bisphosphatase, acid phosphatase, and alkaline phosphatase. Prog Biophys Mol Biol. 2000, 73 (2–4): 263-287.PubMedView ArticleGoogle Scholar
- Fraser HI, Kvaratskhelia M, White MF: The two analogous phosphoglycerate mutases of Escherichia coli. FEBS Lett. 1999, 455 (3): 344-348. 10.1016/S0014-5793(99)00910-2.PubMedView ArticleGoogle Scholar
- Gautam N: Mutated forms of phosphoglycerate mutase in yeast affect reversal of metabolic flux. Effect of reversible and irreversible function of an enzyme on pathway reversal. The Journal of biological chemistry. 1988, 263 (30): 15400-15406.PubMedGoogle Scholar
- Foster JM, Davis PJ, Raverdy S, Sibley MH, Raleigh EA, Kumar S, Carlow CK: Evolution of bacterial phosphoglycerate mutases: non-homologous isofunctional enzymes undergoing gene losses, gains and lateral transfers. PloS one. 2010, 5 (10): e13576-10.1371/journal.pone.0013576.PubMedPubMed CentralView ArticleGoogle Scholar
- Vincent JB, Crowder MW, Averill BA: Hydrolysis of phosphate monoesters: A biological problem with multiple chemical solutions. Trends Biochem Sci. 1992, 17 (3): 105-110. 10.1016/0968-0004(92)90246-6.PubMedView ArticleGoogle Scholar
- Bodansky O: Acid phosphatase. Adv Clin Chem. 1972, 15: 43-147.PubMedView ArticleGoogle Scholar
- Vinopal RT: Microbial metabolism: phosphate metabolism and cellular regulation in microorganisms. Science. 1988, 239 (4839): 513-514.PubMedView ArticleGoogle Scholar
- Coleman JE: Structure and mechanism of alkaline phosphatase. Annu Rev Biophys Biomol Struct. 1992, 21: 441-483. 10.1146/annurev.bb.21.060192.002301.PubMedView ArticleGoogle Scholar
- Lamarche MG, Wanner BL, Crepin S, Harel J: The phosphate regulon and bacterial virulence: a regulatory network connecting phosphate homeostasis and pathogenesis. FEMS Microbiol Rev. 2008, 32 (3): 461-473. 10.1111/j.1574-6976.2008.00101.x.PubMedView ArticleGoogle Scholar
- Dubail I, Berche P, Charbit A: Listeriolysin O as a reporter to identify constitutive and in vivo-inducible promoters in the pathogen Listeria monocytogenes. Infect Immun. 2000, 68 (6): 3242-3250. 10.1128/IAI.68.6.3242-3250.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Polissi A, Pontiggia A, Feger G, Altieri M, Mottl H, Ferrari L, Simon D: Large-scale identification of virulence genes from Streptococcus pneumoniae. Infect Immun. 1998, 66 (12): 5620-5629.PubMedPubMed CentralGoogle Scholar
- Talaat AM, Lyons R, Howard ST, Johnston SA: The temporal expression profile of Mycobacterium tuberculosis infection in mice. Proc Natl Acad Sci U S A. 2004, 101 (13): 4602-4607. 10.1073/pnas.0306023101.PubMedPubMed CentralView ArticleGoogle Scholar
- Merrell DS, Hava DL, Camilli A: Identification of novel factors involved in colonization and acid tolerance of Vibrio cholerae. Mol Microbiol. 2002, 43 (6): 1471-1491. 10.1046/j.1365-2958.2002.02857.x.PubMedView ArticleGoogle Scholar
- Burall LS, Harro JM, Li X, Lockatell CV, Himpsl SD, Hebel JR, Johnson DE, Mobley HL: Proteus mirabilis genes that contribute to pathogenesis of urinary tract infection: identification of 25 signature-tagged mutants attenuated at least 100-fold. Infect Immun. 2004, 72 (5): 2922-2938. 10.1128/IAI.72.5.2922-2938.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. Journal of molecular biology. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Kuznetsova E, Xu L, Singer A, Brown G, Dong A, Flick R, Cui H, Cuff M, Joachimiak A, Savchenko A, et al: Structure and activity of the metal-independent fructose-1,6-bisphosphatase YK23 from Saccharomyces cerevisiae. The Journal of biological chemistry. 2010, 285 (27): 21049-21059. 10.1074/jbc.M110.118315.PubMedPubMed CentralView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Gish W, States DJ: Identification of protein coding regions by database similarity search. Nature genetics. 1993, 3 (3): 266-272. 10.1038/ng0393-266.PubMedView ArticleGoogle Scholar
- Reilly TJ, Felts RL, Henzl MT, Calcutt MJ, Tanner JJ: Characterization of recombinant Francisella tularensis acid phosphatase A. Protein expression and purification. 2006, 45 (1): 132-141. 10.1016/j.pep.2005.05.001.PubMedView ArticleGoogle Scholar
- Aguirre-Garcia MM, Cerbon J, Talamas-Rohana P: Purification and properties of an acid phosphatase from Entamoeba histolytica HM-1:IMSS. Int J Parasitol. 2000, 30 (5): 585-591. 10.1016/S0020-7519(00)00032-1.PubMedView ArticleGoogle Scholar
- Grundner C, Ng HL, Alber T: Mycobacterium tuberculosis protein tyrosine phosphatase PtpB structure reveals a diverged fold and a buried active site. Structure. 2005, 13 (11): 1625-1634. 10.1016/j.str.2005.07.017.PubMedView ArticleGoogle Scholar
- Cowley SC, Babakaiff R, Av-Gay Y: Expression and localization of the Mycobacterium tuberculosis protein tyrosine phosphatase PtpA. Res Microbiol. 2002, 153 (4): 233-241. 10.1016/S0923-2508(02)01309-8.PubMedView ArticleGoogle Scholar
- Boitel B, Ortiz-Lombardia M, Duran R, Pompeo F, Cole ST, Cervenansky C, Alzari PM: PknB kinase activity is regulated by phosphorylation in two Thr residues and dephosphorylation by PstP, the cognate phospho-Ser/Thr phosphatase, in Mycobacterium tuberculosis. Mol Microbiol. 2003, 49 (6): 1493-1508. 10.1046/j.1365-2958.2003.03657.x.PubMedView ArticleGoogle Scholar
- de Souza GA, Leversen NA, Malen H, Wiker HG: Bacterial proteins with cleaved or uncleaved signal peptides of the general secretory pathway. J Proteomics. 2011, 75 (2): 502-510. 10.1016/j.jprot.2011.08.016.PubMedView ArticleGoogle Scholar
- Schnappinger D, Ehrt S, Voskuil MI, Liu Y, Mangan JA, Monahan IM, Dolganov G, Efron B, Butcher PD, Nathan C, et al: Transcriptional adaptation of Mycobacterium tuberculosis within macrophages: insights into the phagosomal environment. J Exp Med. 2003, 198 (5): 693-704. 10.1084/jem.20030846.PubMedPubMed CentralView ArticleGoogle Scholar
- Anderson RG, Hussey H, Baddiley J: The mechanism of wall synthesis in bacteria. The organization of enzymes and isoprenoid phosphates in the membrane. Biochem J. 1972, 127 (1): 11-25.PubMedPubMed CentralView ArticleGoogle Scholar
- Swiezewska E, Danikiewicz W: Polyisoprenoids: structure, biosynthesis and function. Prog Lipid Res. 2005, 44 (4): 235-258. 10.1016/j.plipres.2005.05.002.PubMedView ArticleGoogle Scholar
- Chalker AF, Ingraham KA, Lunsford RD, Bryant AP, Bryant J, Wallis NG, Broskey JP, Pearson SC, Holmes DJ: The bacA gene, which determines bacitracin susceptibility in Streptococcus pneumoniae and Staphylococcus aureus, is also required for virulence. Microbiology. 2000, 146 (Pt 7): 1547-1553.PubMedView ArticleGoogle Scholar
- El Ghachi M, Derbise A, Bouhss A, Mengin-Lecreulx D: Identification of multiple genes encoding membrane proteins with undecaprenyl pyrophosphate phosphatase (UppP) activity in Escherichia coli. The Journal of biological chemistry. 2005, 280 (19): 18689-18695. 10.1074/jbc.M412277200.PubMedView ArticleGoogle Scholar
- Darby CM, Venugopal A, Ehrt S, Nathan CF: Mycobacterium tuberculosis gene Rv2136c is dispensable for acid resistance and virulence in mice. Tuberculosis (Edinb). 2011, 91 (5): 343-347. 10.1016/j.tube.2011.06.007.View ArticleGoogle Scholar
- Bernard R, El Ghachi M, Mengin-Lecreulx D, Chippaux M, Denizot F: BcrC from Bacillus subtilis acts as an undecaprenyl pyrophosphate phosphatase in bacitracin resistance. The Journal of biological chemistry. 2005, 280 (32): 28852-28857. 10.1074/jbc.M413750200.PubMedView ArticleGoogle Scholar
- Tatar LD, Marolda CL, Polischuk AN, van Leeuwen D, Valvano MA: An Escherichia coli undecaprenyl-pyrophosphate phosphatase implicated in undecaprenyl phosphate recycling. Microbiology. 2007, 153 (Pt 8): 2518-2529.PubMedView ArticleGoogle Scholar
- Touze T, Blanot D, Mengin-Lecreulx D: Substrate specificity and membrane topology of Escherichia coli PgpB, an undecaprenyl pyrophosphate phosphatase. The Journal of biological chemistry. 2008, 283 (24): 16573-16583. 10.1074/jbc.M800394200.PubMedView ArticleGoogle Scholar
- Kelley LA, Sternberg MJ: Protein structure prediction on the Web: a case study using the Phyre server. Nature protocols. 2009, 4 (3): 363-371. 10.1038/nprot.2009.2.PubMedView ArticleGoogle Scholar
- Chiba Y, Horita S, Ohtsuka J, Arai H, Nagata K, Igarashi Y, Tanokura M, Ishii M: Structural units important for activity of a novel-type phosphoserine phosphatase from Hydrogenobacter thermophilus TK-6 revealed by crystal structure analysis. The Journal of biological chemistry. 2013, 288 (16): 11448-11458. 10.1074/jbc.M112.449561.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffin JE, Gawronski JD, Dejesus MA, Ioerger TR, Akerley BJ, Sassetti CM: High-resolution phenotypic profiling defines genes essential for mycobacterial growth and cholesterol catabolism. PLoS pathogens. 2011, 7 (9): e1002251-10.1371/journal.ppat.1002251.PubMedPubMed CentralView ArticleGoogle Scholar
- Sheibley RH, Hass LF: Isolation and partial characterization of monophosphoglycerate mutase from human erythrocytes. The Journal of biological chemistry. 1976, 251 (21): 6699-6704.PubMedGoogle Scholar
- Bond CS, White MF, Hunter WN: Mechanistic implications for Escherichia coli cofactor-dependent phosphoglycerate mutase based on the high-resolution crystal structure of a vanadate complex. Journal of molecular biology. 2002, 316 (5): 1071-1081. 10.1006/jmbi.2002.5418.PubMedView ArticleGoogle Scholar
- Rigden DJ, Alexeev D, Phillips SE, Fothergill-Gilmore LA: The 2.3 A X-ray crystal structure of S. cerevisiae phosphoglycerate mutase. Journal of molecular biology. 1998, 276 ((2): 449-459.PubMedView ArticleGoogle Scholar
- Solem C, Petranovic D, Koebmann B, Mijakovic I, Jensen PR: Phosphoglycerate mutase is a highly efficient enzyme without flux control in Lactococcus lactis. J Mol Microbiol Biotechnol. 2010, 18 (3): 174-180. 10.1159/000315458.PubMedView ArticleGoogle Scholar
- Studier FW, Rosenberg AH, Dunn JJ, Dubendorff JW: Use of T7 RNA polymerase to direct expression of cloned genes. Methods in enzymology. 1990, 185: 60-89.PubMedView ArticleGoogle Scholar
- van Soolingen D, de Haas PE, Hermans PW, van Embden JD: DNA fingerprinting of Mycobacterium tuberculosis. Methods in enzymology. 1994, 235: 196-205.PubMedView ArticleGoogle Scholar
- Rigden DJ, Mello LV, Setlow P, Jedrzejas MJ: Structure and mechanism of action of a cofactor-dependent phosphoglycerate mutase homolog from Bacillus stearothermophilus with broad specificity phosphatase activity. Journal of molecular biology. 2002, 315 (5): 1129-1143. 10.1006/jmbi.2001.5290.PubMedView ArticleGoogle Scholar
- White MF, Fothergill-Gilmore LA: Development of a mutagenesis, expression and purification system for yeast phosphoglycerate mutase. Investigation of the role of active-site His181. Eur J Biochem. 1992, 207 (2): 709-714. 10.1111/j.1432-1033.1992.tb17099.x.PubMedView ArticleGoogle Scholar
- Geladopoulos TP, Sotiroudis TG, Evangelopoulos AE: A malachite green colorimetric assay for protein phosphatase activity. Anal Biochem. 1991, 192 (1): 112-116. 10.1016/0003-2697(91)90194-X.PubMedView ArticleGoogle Scholar
- Kao FF, Mahmuda S, Pinto R, Triccas JA, West NP, Britton WJ: The secreted lipoprotein, MPT83, of Mycobacterium tuberculosis is recognized during human tuberculosis and stimulates protective immunity in mice. PloS one. 2012, 7 (5): e34991-10.1371/journal.pone.0034991.PubMedPubMed CentralView ArticleGoogle Scholar
- Hedrick JL, Smith AJ: Size and charge isomer separation and estimation of molecular weights of proteins by disc gel electrophoresis. Arch Biochem Biophys. 1968, 126 (1): 155-164. 10.1016/0003-9861(68)90569-9.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.