Metallophosphoesterases In Mycobacteria Enigmatic Roles In Regulating Mycobacterial Physiology
Pathogenic bacteria such as M.tuberculosis have evolved several mechanisms to aid their intracellular survival and subvert host defenses. One of the contributing factors is thought to be the production and secretion of large amount of cAMP, Mycobacterial genomes encode a large number of adenylyl cyclases distinct in their structure and regulatory mechanisms. The roles of these enzymes in the physiology and pathogenesis of virulent mycobacteria are only now being elucidated. The roles of phosphodiesterases (PDEs), which serve to lower cAMP levels through degradation are, however, relatively unexplored. The Rv0805 gene was previously shown to code for an active phosphodiesterase from Mycobacterium tuberculosis. Bioinformatics analysis revealed that orthologs of Rv0805 were found even in eukaryotes. Biochemical and structural characterization of Rv0805 revealed that it was a class III cAMP phosphodiesterase. Comparative genomics identified a close ortholog of Rv0805 in M. leprae (ML2210). The genome of M. leprae Encodes only 1,604 predicted proteins and possesses the highest number of pseudogenes, 1,116. The retention of a functional PDE, the ortholog of Rv0805, in the minimal genome of M. leprae is indicative of its importance in cellular physiology. Biochemical characterization of proteins from M. leprae and use of heterologous hosts will help understand this human pathogen better, since there are no tools currently available to genetically manipulate this bacterium. Sequence analysis of ML2210 revealed the presence of conserved motifs and residues known to be critical for catalysis and unique to class III phosphodiesterases. ML2210 shares 83% sequence identity with Rv0805 and 24% sequence identity with the phosphodiesterase from E. coli (cpdA). In vitro biochemical characterization of ML2210 using non-nucleotide colorigenic and cyclic nucleotide substrates revealed that it was an enzymatically active phosphodiesterase. Kinetic parameters of ML2210 with respect ot colorigenic substrates revealed that its catalytic properties were similar to that of Rv0805. However, with respect to hydrolysis of 3’, 5’-cAMP, ML2210 was catalytically more efficient than Rv0805, suggesting that in spite of being orthologs, these enzymes have evolved distinct specificities at their active site. A parallel of monoclonal antibodies raised to Rv0805 was also used understand the differences in the biochemical properties of Rv0805 and ML2210 better. It was observed that only one monoclonal antibody was able to recognize ML2210 by ELISA and not by Western blot analysis. These results revealed that conformational differences between ML2210 and Rv0805 exist. Over-expression of ML2210 in M. smegmatis resulted in a modest decrease in intracellular cAMP levels. Despite the absence of a predicted transmembrane region or a membrane-targeting signal, ML2210 localized to the cell envelop fraction upon over expression in M. smegmatis. Moreover, like Rv0805, over-expression of ML2210 also resulted in perturbation of the cell wall of M. smegmatis, arguing for additional cellular roles of this protein. Orthologs of Rv0805 or ML2210 are found only in slow growing mycobacteria suggesting that other cyclic nucleotide phosphodiesterases could regulate cAMP levels in fast growing mycobacteria like M. smegmatis. Since BLAST results did not retrieve an ortholog of Rv0805 or ML2210, COG1409 (COG database) containing Rv0805 was examined for the presence of other mycobacterial phosphodiesterases. Bioinformatics analysis identified Rv2795c as another PDE from M. tuberculosis. Sequence analysis of Rv2795c revealed the presence of all the motifs conserved in the class III PDEs but Rv2795c shared only 22% sequence identity with Rv0805 and 19% sequence identity with CpdA. Importantly, an ortholog of Rv2795c was identified in M. leprae. Interestingly. Rv2795c and its orthologs branched away from Rv0805, making it phylogenetically distinct and hence warranting further characterization. Recombinant, purified MSMEG_2647 (the Rv2795c ortholog from M. smegmatis) was able to hydrolyze cyclic nucleotides and other phosphodiester substrates in vitro. The Km for colorigenic substrates was higher when compared to the Km of ML2210 or Rv0805 for these substrates. However, the kinetic parameters of MSMEG_2647 for cyclic nucleotides were comparable to those of ML2210 or Rv0805. MSMEG_2647 was a metal dependent enzyme and among the panel of metals tested, Mn2+ supported the highest in vitro catalytic activity of MSMEG_2647. Zn2+ inhibited the catalytic activity of MSMEG_2647. In order to gain insight into the catalysis of MSMEG_2647, the end products of cAMP hydrolysis by MSMEG_2647 were analysed using reverse phase HPLC. The assay revealed that the end products of cyclic nucleotide hydrolysis by MSMEG_2647 were different when compared to the end products of hydrolysis of the same substrates by Rv0805 or ML2210. This suggests differences in the architecture of the active site residues of the mycobacterial MPEs. A mutational anlaysis of the active site residues in MSMEG_2647 was carried out to identify residues involved in substrate recognition and metal coordination. Although Rv0805 and MSMEG_2647 shared only a 22% sequence identity, MSMEG_2647 displayed strict conservation in the core MPE motifs. Mutation of the active residues N97 and H98 in Rv0805 had led to an abrogation of its catalytic activity. However, corresponding mutations of N76A and H77A in MSMEG_2647, did not lead to a loss in its catalytic activity. A third mutation known to be important for the catalytic activity of Rv0805 (D19) was incorporated. The corresponding residue at D19 position was mutated to an alanine. The catalytic activity of MSMEG_2647D19AN76AH77A mutant was abrogated, suggesting that while the core MPE motifs are conserved between mycobacterial PDEs, differences in the ensemble of the active site residues contributing to their catalytic activity exist. Thus, at least two biochemically diverse PDE clades are found in mycobacterial species. In order to decipher the function of MSMEG_2647, its expression was monitored during the growth of M. Smegmatis. The promoter of MSMEG_2647 displayed maximum activity during the logarithmic phase of M. smegmatis growth after which its activity declined as M. smegmatis entered the stationary phase. However in contrast to this, the transcript corresponding to msmeg_2647 mRNA was found at both logarithmic and stationary phases. The MSMEG_2647 protein was also detected at both logarithmic and stationary phases of M. smegmatis. These results suggest that additional factors may contribute to the stability of msmeg_2647 mRNA and protein levels. Localization studies of MSMEG_2647 revealed that MSMEG_2647 was present in the cytosol as well as in the cell envelope fractions. Interestingly, over-expression of MSMEG_2647 did not result in a significant increase in PDE activity in various subcellular fractions, suggesting tight regulation on the in vivo activity in various subcellular fractions, suggesting tight regulation on the in vivo activity of MSMEG_2647. In addition, over-expression of MSMEG_2647 in M. smegmatis led to only a modest decrease in cAMP levels in M. smegmatis. These results suggested additional roles of MSMEG_2647 in the biology of mycobacteria. Overexpression of MSMEG_2647 peturbed the integrity of cell wall as assessed by the use of lipophillic indicators of cell growth, crystal violet and malachite green, and a cell wall targeting antibiotic, isoniazid. Analyzing the gene neighborhood of MSMEG_2647 provided an insight into its putative function. It was observed that the stop codon of msmeg_2647 overlapped with the start codon of msmeg_2648 and stop codon of msmeg-2648 overlapped with the start codon of msmeg_2649. RT PCR was carried out at logarhtimic and stationary phases of M. smegmatis growth, which revealed that a polycistronic mRNA was being transcribed. These results confirmed that msmeg_2647, msmeg_2648 and msmeg_2649 were a part of an operon. Interestingly, these three genes as a gene cluster were confined to only those actinobacteria that produced mycolic acids. An operon often encodes products that form multiprotein complexes and operate in a common pathway. Since there were a part of an operon, a GST pull-down approach was employed to test if MSMEG_2647, MSMEG_2648 and MSMEG_2649 could physically interact. It was observed that MSMEG_2647 interacted with MSMEG_2648 and MSMEG_2649. MSMEG_2648 in turn interacted with MSMEG_2649. A role for MSMEG_2647 as a scaffold recruiting MSMEG_2648 and MSMEG_2649 is therefore proposed. In turn, a complex formation with these proteins may regulate the activity of MSMEG_2647. Attempts to generate a knock out of msmeg_2647 in M. smegmatis by homologous recombination were not successful suggesting either the gene was essential or a polar effect on msmeg_2648(an essential gene for the viability of M. smegmatis) may not allow msmeg_2647 to be deleted from the genome of M. smegmatis. In summary, this study has identified and characterized two new phosphodiesterases from mycobacteria, one from the pathogenic mycobacterium, M. leprae and the other, a PDE from M. smegmatis that is conserved in all species of mycobacteria. Several, key biochemical differences were observed using biochemical and biological approaches. It appears that the cellular roles of mycobacterial phsophodiesterases may extend beyond cAMP hydrolysis, with these proteins not only regulating cell wall properties but also acting as scaffolding proteins in the cell.