Phylogenetic relationships and classification of thiolases and thiolase-like proteins of Mycobacterium tuberculosis and Mycobacterium smegmatis

Thiolases are enzymes involved in lipid metabolism. Thiolases remove the acetyl-CoA moiety from 3-ketoacyl-CoAs in the degradative reaction. They can also catalyze the reverse Claisen condensation reaction, which is the first step of biosynthetic processes such as the biosynthesis of sterols and ketone bodies. In human, six distinct thiolases have been identified. Each of these thiolases is different from the other with respect to sequence, oligomeric state, substrate specificity and subcellular localization. Four sequence fingerprints, identifying catalytic loops of thiolases, have been described. In this study genome searches of two mycobacterial species (Mycobacterium tuberculosis and Mycobacterium smegmatis), were carried out, using the six human thiolase sequences as queries. Eight and thirteen different thiolase sequences were identified in M. tuberculosis and M. smegmatis, respectively. In addition, thiolase-like proteins (one encoded in the Mtb and two in the Msm genome) were found. The purpose of this study is to classify these mostly uncharacterized thiolases and thiolase-like proteins. Several other sequences obtained by searches of genome databases of bacteria, mammals and the parasitic protist family of the Trypanosomatidae were included in the analysis. Thiolase-like proteins were also found in the trypanosomatid genomes, but not in those of mammals. In order to study the phylogenetic relationships at a high confidence level, additional thiolase sequences were included such that a total of 130 thiolases and thiolase-like protein sequences were used for the multiple sequence alignment. The resulting phylogenetic tree identifies 12 classes of sequences, each possessing a characteristic set of sequence fingerprints for the catalytic loops. From this analysis it is now possible to assign the mycobacterial thiolases to corresponding homologues in other kingdoms of life. The results of this bioinformatics analysis also show interesting differences between the distributions of M. tuberculosis and M. smegmatis thiolases over the 12 different classes. (C) 2014 Elsevier Ltd. All rights reserved.

Targeting Mycobacterium tuberculosis nucleoid-associated protein HU with structure-based inhibitors

The nucleoid-associated protein HU plays an important role in maintenance of chromosomal architecture and in global regulation of DNA transactions in bacteria. Although HU is essential for growth in Mycobacterium tuberculosis (Mtb), there have been no reported attempts to perturb HU function with small molecules. Here we report the crystal structure of the N-terminal domain of HU from Mtb. We identify a core region within the HU-DNA interface that can be targeted using stilbene derivatives. These small molecules specifically inhibit HU-DNA binding, disrupt nucleoid architecture and reduce Mtb growth. The stilbene inhibitors induce gene expression changes in Mtb that resemble those induced by HU deficiency. Our results indicate that HU is a potential target for the development of therapies against tuberculosis.

Crystal structure of a putative aspartic proteinase domain of the Mycobacterium tuberculosis cell surface antigen PE_PGRS16

We report the crystal structure of the first prokaryotic aspartic proteinase-like domain identified in the genome of Mycobacterium tuberculosis. A search in the genomes of Mycobacterium species showed that the C-terminal domains of some of the PE family proteins contain two classic DT/SG motifs of aspartic proteinases with a low overall sequence similarity to HIV proteinase. The three-dimensional structure of one of them, Rv0977 (PE_PGRS16) of M. tuberculosis revealed the characteristic pepsinf-old and catalytic site architecture. However, the active site was completely blocked by the N-terminal His-tag. Surprisingly, the enzyme was found to be inactive even after the removal of the N-terminal His-tag. A comparison of the structure with pepsins showed significant differences in the critical substrate binding residues and in the flap tyrosine conformation that could contribute to the lack of proteolytic activity of Rv0977. (C) 2013 The Authors. Published by Elsevier B.V. on behalf of Federation of European Biochemical Societies. All rights reserved.

Mycobacterium bovis BCG promotes tumor cell survival from tumor necrosis factor-alpha-induced apoptosis

Background: Increased incidence of lung cancer among pulmonary tuberculosis patients suggests mycobacteria-induced tumorigenic response in the host. The alveolar epithelial cells, candidate cells that form lung adenocarcinoma, constitute a niche for mycobacterial replication and infection. We thus explored the possible mechanism of M. bovis Bacillus Calmette-Guerin (BCG)-assisted tumorigenicity in type II epithelial cells, human lung adenocarcinoma A549 and other cancer cells. Methods: Cancer cell lines originating from lung, colon, bladder, liver, breast, skin and cervix were treated with tumor necrosis factor (TNF)-alpha in presence or absence of BCG infection. p53, COP1 and sonic hedgehog (SHH) signaling markers were determined by immunoblotting and luciferase assays, and quantitative real time PCR was done for p53-responsive pro-apoptotic genes and SHH signaling markers. MTT assays and Annexin V staining were utilized to study apoptosis. Gain-and loss-of-function approaches were used to investigate the role for SHH and COP1 signaling during apoptosis. A549 xenografted mice were used to validate the contribution of BCG during TNF-alpha treatment. Results: Here, we show that BCG inhibits TNF-alpha-mediated apoptosis in A549 cells via downregulation of p53 expression. Substantiating this observation, BCG rescued A549 xenografts from TNF-alpha-mediated tumor clearance in nude mice. Furthermore, activation of SHH signaling by BCG induced the expression of an E3 ubiquitin ligase, COP1. SHH-driven COP1 targeted p53, thereby facilitating downregulation of p53-responsive pro-apoptotic genes and inhibition of apoptosis. Similar effects of BCG could be shown for HCT116, T24, MNT-1, HepG2 and HELA cells but not for HCT116 p53(-/-) and MDA-MB-231 cells. Conclusion: Our results not only highlight possible explanations for the coexistence of pulmonary tuberculosis and lung cancer but also address probable reasons for failure of BCG immunotherapy of cancers.

Genomic mapping of cAMP receptor protein (CRPMt) in Mycobacterium tuberculosis: relation to transcriptional start sites and the role of CRPMt as a transcription factor

Chromatin immunoprecipitation identified 191 binding sites of Mycobacterium tuberculosis cAMP receptor protein (CRPMt) at endogenous expression levels using a specific alpha-CRPMt antibody. Under these native conditions an equal distribution between intragenic and intergenic locations was observed. CRPMt binding overlapped a palindromic consensus sequence. Analysis by RNA sequencing revealed widespread changes in transcriptional profile in a mutant strain lacking CRPMt during exponential growth, and in response to nutrient starvation. Differential expression of genes with a CRPMt-binding site represented only a minor portion of this transcriptional reprogramming with similar to 19% of those representing transcriptional regulators potentially controlled by CRPMt. The subset of genes that are differentially expressed in the deletion mutant under both culture conditions conformed to a pattern resembling canonical CRP regulation in Escherichia coli, with binding close to the transcriptional start site associated with repression and upstream binding with activation. CRPMt can function as a classical transcription factor in M. tuberculosis, though this occurs at only a subset of CRPMt-binding sites.

Mycobacterium tuberculosis Rho Is an NTPase with Distinct Kinetic Properties and a Novel RNA-Binding Subdomain

Two mechanisms - factor independent and dependent termination - ensure the completion of RNA synthesis in eubacteria. Factor-dependent mechanism relies on the Rho protein to terminate transcription by interacting with RNA polymerase. Although well studied in Escherichia coli, the properties of the Rho homologs from most bacteria are not known. The rho gene is unusually large in genus Mycobacterium and other members of actinobacteria, having,150 additional residues towards the amino terminal end. We describe the distinct properties of Rho from Mycobacterium tuberculosis. It is an NTPase with a preference for purine nucleoside triphosphates with kinetic properties different from E. coli homolog and an ability to use various RNA substrates. The N-terminal subdomain of MtbRho can bind to RNA by itself, and appears to contribute to the interaction of the termination factor with RNAs. Furthermore, the interaction with RNA induces changes in conformation and oligomerization of MtbRho.

Characterizing the pocketome of Mycobacterium tuberculosis and application in rationalizing polypharmacological target selection

Polypharmacology is beginning to emerge as an important concept in the field of drug discovery. However, there are no established approaches to either select appropriate target sets or design polypharmacological drugs. Here, we propose a structural-proteomics approach that utilizes the structural information of the binding sites at a genome-scale obtained through in-house algorithms to characterize the pocketome, yielding a list of ligands that can participate in various biochemical events in the mycobacterial cell. The pocket-type space is seen to be much larger than the sequence or fold-space, suggesting that variations at the site-level contribute significantly to functional repertoire of the organism. All-pair comparisons of binding sites within Mycobacterium tuberculosis (Mtb), pocket-similarity network construction and clustering result in identification of binding-site sets, each containing a group of similar binding sites, theoretically having a potential to interact with a common set of compounds. A polypharmacology index is formulated to rank targets by incorporating a measure of druggability and similarity to other pockets within the proteome. This study presents a rational approach to identify targets with polypharmacological potential along with possible drugs for repurposing, while simultaneously, obtaining clues on lead compounds for use in new drug-discovery pipelines.

Mycobacterium tuberculosis DinG Is a Structure-specific Helicase That Unwinds G4 DNA IMPLICATIONS FOR TARGETING G4 DNA AS A NOVEL THERAPEUTIC APPROACH

The significance of G-quadruplexes and the helicases that resolve G4 structures in prokaryotes is poorly understood. The Mycobacterium tuberculosis genome is GC-rich and contains >10,000 sequences that have the potential to form G4 structures. In Escherichia coli, RecQ helicase unwinds G4 structures. However, RecQ is absent in M. tuberculosis, and the helicase that participates in G4 resolution in M. tuberculosis is obscure. Here, we show that M. tuberculosis DinG (MtDinG) exhibits high affinity for ssDNA and ssDNA translocation with a 5' -> 3' polarity. Interestingly, MtDinG unwinds overhangs, flap structures, and forked duplexes but fails to unwind linear duplex DNA. Our data with DNase I footprinting provide mechanistic insights and suggest that MtDinG is a 5' -> 3' polarity helicase. Notably, in contrast to E. coli DinG, MtDinG catalyzes unwinding of replication fork and Holliday junction structures. Strikingly, we find that MtDinG resolves intermolecular G4 structures. These data suggest that MtDinG is a multifunctional structure-specific helicase that unwinds model structures of DNA replication, repair, and recombination as well as G4 structures. We finally demonstrate that promoter sequences of M. tuberculosis PE_PGRS2, mce1R, and moeB1 genes contain G4 structures, implying that G4 structures may regulate gene expression in M. tuberculosis. We discuss these data and implicate targeting G4 structures and DinG helicase in M. tuberculosis could be a novel therapeutic strategy for culminating the infection with this pathogen.

Structure, interactions and evolutionary implications of a domain-swapped lectin dimer from Mycobacterium smegmatis

Crystal structure determination of the lectin domain of MSMEG_3662 from Mycobacterium smegmatis and its complexes with mannose and methyl-alpha-mannose, the first effort of its kind on a mycobacterial lectin, reveals a structure very similar to beta-prism II fold lectins from plant sources, but with extensive unprecedented domain swapping in dimer formation. The two subunits in a dimer often show small differences in structure, but the two domains, not always related by 2-fold symmetry, have the same structure. Each domain carries three sugar-binding sites, similar to those in plant lectins, one on each Greek key motif. The occurrence of beta-prism II fold lectins in bacteria, with characteristics similar to those from plants, indicates that this family of lectins is of ancient origin and had evolved into a mature system before bacteria and plants diverged. In plants, the number of binding sites per domain varies between one and three, whereas the number is two in the recently reported lectin domains from Pseudomonas putida and Pseudomonas aeruginosa. An analysis of the sequences of the lectins and the lectin domains shows that the level of sequence similarity among the three Greek keys in each domain has a correlation with the number of binding sites in it. Furthermore, sequence conservation among the lectins from different species is the highest for that Greek key which carries a binding site in all of them. Thus, it would appear that carbohydrate binding influences the course of the evolution of the lectin.

Characterization of a dual-active enzyme, DcpA, involved in cyclic diguanosine monophosphate turnover in Mycobacterium smegmatis

We have reported previously that the long-term survival of Mycobacterium smegmatis is facilitated by a dual-active enzyme MSDGC-1 (renamed DcpA), which controls the cellular turnover of cyclic diguanosine monophosphate (c-di-GMP). Most mycobacterial species possess at least a single copy of a DcpA orthologue that is highly conserved in terms of sequence similarity and domain architecture. Here, we show that DcpA exists in monomeric and dimeric forms. The dimerization of DcpA is due to non-covalent interactions between two protomers that are arranged in a parallel orientation. The dimer shows both synthesis and hydrolysis activities, whereas the monomer shows only hydrolysis activity. In addition, we have shown that DcpA is associated with the cytoplasmic membrane and exhibits heterogeneous cellular localization with a predominance at the cell poles. Finally, we have also shown that DcpA is involved in the change in cell length and colony morphology of M. smegmatis. Taken together, our study provides additional evidence about the role of the bifunctional protein involved in c-di-GMP signalling in M. smegmatis.

Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

Molecular Dissection of Mycobacterium tuberculosis Integration Host Factor Reveals Novel Insights into the Mode of DNA Binding and Nucleoid Compaction

Background: mIHF belongs to a subfamily of proteins, distinct from E. coli IHF. Results: Functionally important amino acids of mIHF and the mechanism(s) underlying DNA binding, DNA bending, and site-specific recombination are distinct from that of E. coli IHF. Conclusion: mIHF functions could contribute beyond nucleoid compaction. Significance: Because mIHF is essential for growth, the molecular mechanisms identified here can be exploited in drug screening efforts. The annotated whole-genome sequence of Mycobacterium tuberculosis revealed that Rv1388 (Mtihf) is likely to encode for a putative 20-kDa integration host factor (mIHF). However, very little is known about the functional properties of mIHF or the organization of the mycobacterial nucleoid. Molecular modeling of the mIHF three-dimensional structure, based on the cocrystal structure of Streptomyces coelicolor IHF duplex DNA, a bona fide relative of mIHF, revealed the presence of Arg-170, Arg-171, and Arg-173, which might be involved in DNA binding, and a conserved proline (Pro-150) in the tight turn. The phenotypic sensitivity of Escherichia coli ihfA and ihfB strains to UV and methyl methanesulfonate could be complemented with the wild-type Mtihf but not its alleles bearing mutations in the DNA-binding residues. Protein-DNA interaction assays revealed that wild-type mIHF, but not its DNA-binding variants, binds with high affinity to fragments containing attB and attP sites and curved DNA. Strikingly, the functionally important amino acid residues of mIHF and the mechanism(s) underlying its binding to DNA, DNA bending, and site-specific recombination are fundamentally different from that of E. coli IHF. Furthermore, we reveal novel insights into IHF-mediated DNA compaction depending on the placement of its preferred binding sites; mIHF promotes DNA compaction into nucleoid-like or higher order filamentous structures. We therefore propose that mIHF is a distinct member of a subfamily of proteins that serve as essential cofactors in site-specific recombination and nucleoid organization and that these findings represent a significant advance in our understanding of the role(s) of nucleoid-associated proteins.

The adenylyl cyclase Rv2212 modifies the proteome and infectivity of Mycobacterium bovis BCG

All organisms have the capacity to sense and respond to environmental changes. These signals often involve the use of second messengers such as cyclic adenosine monophosphate (cAMP). This second messenger is widely distributed among organisms and coordinates gene expression related with pathogenesis, virulence, and environmental adaptation. Genomic analysis in Mycobacterium tuberculosis has identified 16 adenylyl cyclases (AC) and one phosphodiesterase, which produce and degrade cAMP, respectively. To date, ten AC have been biochemically characterized and only one (Rv0386) has been found to be important during murine infection with M. tuberculosis. Here, we investigated the impact of hsp60-driven Rv2212 gene expression in Mycobacterium bovis Bacillus Calmette-Guerin (BCG) during growth in vitro, and during macrophage and mice infection. We found that hsp60-driven expression of Rv2212 resulted in an increased capacity of replication in murine macrophages but an attenuated phenotype in lungs and spleen when administered intravenously in mice. Furthermore, this strain displayed an altered proteome mainly affecting proteins associated with stress conditions (bfrB, groEL-2, DnaK) that could contribute to the attenuated phenotype observed in mice.

Enriching the annotation of Mycobacterium tuberculosis H37Rv proteome using remote homology detection approaches: Insights into structure and function

The availability of the genome sequence of Mycobacterium tuberculosis H37Rv has encouraged determination of large numbers of protein structures and detailed definition of the biological information encoded therein; yet, the functions of many proteins in M. tuberculosis remain unknown. The emergence of multidrug resistant strains makes it a priority to exploit recent advances in homology recognition and structure prediction to re-analyse its gene products. Here we report the structural and functional characterization of gene products encoded in the M. tuberculosis genome, with the help of sensitive profile-based remote homology search and fold recognition algorithms resulting in an enhanced annotation of the proteome where 95% of the M. tuberculosis proteins were identified wholly or partly with information on structure or function. New information includes association of 244 proteins with 205 domain families and a separate set of new association of folds to 64 proteins. Extending structural information across uncharacterized protein families represented in the M. tuberculosis proteome, by determining superfamily relationships between families of known and unknown structures, has contributed to an enhancement in the knowledge of structural content. In retrospect, such superfamily relationships have facilitated recognition of probable structure and/or function for several uncharacterized protein families, eventually aiding recognition of probable functions for homologous proteins corresponding to such families. Gene products unique to mycobacteria for which no functions could be identified are 183. Of these 18 were determined to be M. tuberculosis specific. Such pathogen-specific proteins are speculated to harbour virulence factors required for pathogenesis. A re-annotated proteome of M. tuberculosis, with greater completeness of annotated proteins and domain assigned regions, provides a valuable basis for experimental endeavours designed to obtain a better understanding of pathogenesis and to accelerate the process of drug target discovery. (C) 2014 Elsevier Ltd. All rights reserved.

Paralogous cAMP Receptor Proteins in Mycobacterium smegmatis Show Biochemical and Functional Divergence

The cyclic AMP receptor protein (CRP) family of transcription factors consists of global regulators of bacterial gene expression. Here, we identify two paralogous CRPs in the genome of Mycobacterium smegmatis that have 78% identical sequences and characterize them biochemically and functionally. The two proteins (MSMEG_0539 and MSMEG_6189) show differences in cAMP binding affinity, trypsin sensitivity, and binding to a CRP site that we have identified upstream of the msmeg_3781 gene. MSMEG_6189 binds to the CRP site readily in the absence of cAMP, while MSMEG_0539 binds in the presence of cAMP, albeit weakly. msmeg_6189 appears to be an essential gene, while the ?msmeg_0539 strain was readily obtained. Using promoter-reporter constructs, we show that msmeg_3781 is regulated by CRP binding, and its transcription is repressed by MSMEG_6189. Our results are the first to characterize two paralogous and functional CRPs in a single bacterial genome. This gene duplication event has subsequently led to the evolution of two proteins whose biochemical differences translate to differential gene regulation, thus catering to the specific needs of the organism.