In silico analysis and characterization of GntR family of regulators from Mycobacterium tuberculosis

Mce1R of Mycobacterium tuberculosis prefers long-chain fatty acids as specific ligands: a computational study

The mce1 operon of Mycobacterium tuberculosis, which codes the Mce1 transporter, facilitates the transport of fatty acids. Fatty acids are one of the major sources for carbon and energy for the pathogen during its intracellular survival and pathogenicity. The mce1 operon is transcriptionally regulated by Mce1R, a VanR-type regulator, which could bind specific ligands and control the expression of the mce1 operon accordingly. This work reports computational identification of Mce1R-specific ligands. Initially by employing cavity similarity search algorithm by the ProBis server, the cavities of the proteins similar to that of Mce1R and the bound ligands were identified from which fatty acids were selected as the potential ligands. From the earlier-generated monomeric structure, the dimeric structure of Mce1R was then modeled by the GalaxyHomomer server and validated computationally to use in molecular docking and molecular dynamics simulation analysis. The fatty acid ligands were found to dock within the cavity of Mce1R and the docked complexes were subjected to molecular dynamics simulation to explore their stabilities and other dynamic properties. The data suggest that Mce1R preferably binds to long-chain fatty acids and undergoes distinct structural changes upon binding.

PcaR, a GntR/FadR Family Transcriptional Repressor Controls the Transcription of Phenazine-1-Carboxylic Acid 1,2-Dioxygenase Gene Cluster in Sphingomonas histidinilytica DS-9

In our previous study, the phenazine-1-carboxylic acid (PCA) 1,2-dioxygenase gene cluster (pcaA1A2A3A4 cluster) in Sphingomonas histidinilytica DS-9 was identified to be responsible for the conversion of PCA to 1,2-dihydroxyphenazine (Ren Y, Zhang M, Gao S, Zhu Q, et al. 2022. Appl Environ Microbiol 88:e00543-22). However, the regulatory mechanism of the pcaA1A2A3A4 cluster has not been elucidated yet. In this study, the pcaA1A2A3A4 cluster was found to be transcribed as two divergent operons: pcaA3-ORF5205 (named A3-5205 operon) and pcaA1A2-ORF5208-pcaA4-ORF5210 (named A1-5210 operon). The promoter regions of the two operons were overlapped. PcaR acts as a transcriptional repressor of the pcaA1A2A3A4 cluster, and it belongs to GntR/FadR family transcriptional regulator. Gene disruption of pcaR can shorten the lag phase of PCA degradation. The results of electrophoretic mobility shift assay and DNase I footprinting showed that PcaR binds to a 25-bp motif in the ORF5205-pcaA1 intergenic promoter region to regulate the expression of two operons. The 25-bp motif covers the -10 region of the promoter of A3-5205 operon and the -35 region and -10 region of the promoter of A1-5210 operon. The TNGT/ANCNA box within the motif was essential for PcaR binding to the two promoters. PCA acted as an effector of PcaR, preventing it from binding to the promoter region and repressing the transcription of the pcaA1A2A3A4 cluster. In addition, PcaR represses its own transcription, and this repression can be relieved by PCA. This study reveals the regulatory mechanism of PCA degradation in strain DS-9, and the identification of PcaR increases the variety of regulatory model of the GntR/FadR-type regulator. IMPORTANCE Sphingomonas histidinilytica DS-9 is a phenazine-1-carboxylic acid (PCA)-degrading strain. The 1,2-dioxygenase gene cluster (pcaA1A2A3A4 cluster, encoding dioxygenase PcaA1A2, reductase PcaA3, and ferredoxin PcaA4) is responsible for the initial degradation step of PCA and widely distributed in Sphingomonads, but its regulatory mechanism has not been investigated yet. In this study, a GntR/FadR-type transcriptional regulator PcaR repressing the transcription of pcaA1A2A3A4 cluster and pcaR gene was identified and characterized. The binding site of PcaR in ORF5205-pcaA1 intergenic promoter region contains a TNGT/ANCNA box, which is important for the binding. These findings enhance our understanding of the molecular mechanism of PCA degradation.

Transcriptional regulation and drug resistance in Mycobacterium tuberculosis

Bacterial drug resistance is one of the major challenges to present and future human health, as the continuous selection of multidrug resistant bacteria poses at serious risk the possibility to treat infectious diseases in the near future. One of the infection at higher risk to become incurable is tuberculosis, due to the few drugs available in the market against Mycobacterium tuberculosis. Drug resistance in this species is usually due to point mutations in the drug target or in proteins required to activate prodrugs. However, another interesting and underexplored aspect of bacterial physiology with important impact on drug susceptibility is represented by the changes in transcriptional regulation following drug exposure. The main regulators involved in this phenomenon in M. tuberculosis are the sigma factors, and regulators belonging to the WhiB, GntR, XRE, Mar and TetR families. Better understanding the impact of these regulators in survival to drug treatment might contribute to identify new drug targets and/or to design new strategies of intervention.

In Silico Characterization and Virtual Screening of GntR/HutC Family Transcriptional Regulator MoyR: A Potential Monooxygenase Regulator in Mycobacterium tuberculosis

Simple Summary

In an era where the world faces new diseases and pathogens, another emerging challenge is neglected pathogens becoming more notorious. Transcriptional regulators play a vital role in the pathogenesis and survival of these pathogens. Hence, characterizing transcriptional regulators, either in vitro or in silico, is of great importance. Here, we present the first structural characterization of a GntR/HutC regulator in Mycobacterium tuberculosis via in silico methods. We have suggested its possible role and potential as a drug target as well as identified possible drug leads that can be used for further improvements.

Abstract

Mycobacterium tuberculosis is a well-known pathogen due to the emergence of drug resistance associated with it, where transcriptional regulators play a key role in infection, colonization and persistence. The genome of M. tuberculosis encodes many transcriptional regulators, and here we report an in-depth in silico characterization of a GntR regulator: MoyR, a possible monooxygenase regulator. Homology modelling provided a reliable structure for MoyR, showing homology with a HutC regulator DasR from Streptomyces coelicolor. In silico physicochemical analysis revealed that MoyR is a cytoplasmic protein with higher thermal stability and higher pI. Four highly probable binding pockets were determined in MoyR and the druggability was higher in the orthosteric binding site consisting of three conserved critical residues: TYR179, ARG223 and GLU234. Two highly conserved leucine residues were identified in the effector-binding region of MoyR and other HutC homologues, suggesting that these two residues can be crucial for structure stability and oligomerization. Virtual screening of drug leads resulted in four drug-like compounds with greater affinity to MoyR with potential inhibitory effects for MoyR. Our findings support that this regulator protein can be valuable as a therapeutic target that can be used for developing drug leads.

Molecular Cloning, Purification and Characterization of Mce1R of Mycobacterium tuberculosis

The mce1 operon of Mycobacterium tuberculosis, important for lipid metabolism/transport, host cell invasion, modulation of host immune response and pathogenicity, is under the transcriptional control of Mce1R. Hence characterizing Mce1R is an important step for novel anti-tuberculosis drug discovery. The present study reports functional and in silico characterization of Mce1R. In this work, we have computationally modeled the structure of Mce1R and have validated the structure by computational and experimental methods. Mce1R has been shown to harbor the canonical VanR-like structure with a flexible N-terminal 'arm', carrying conserved positively charged residues, most likely involved in the operator DNA binding. The mce1R gene has been cloned, expressed, purified and its DNA-binding activity has been measured in vitro. The K_d value for Mce1R-operator DNA interaction has been determined to be 0.35 ± 0.02 µM which implies that Mce1R binds to DNA with moderate affinity compared to the other FCD family of regulators. So far, this is the first report for measuring the DNA-binding affinity of any VanR-type protein. Despite significant sequence similarity at the N-terminal domain, the wHTH motif of Mce1R exhibits poor conservancy of amino acid residues, critical for DNA-binding, thus results in moderate DNA-binding affinity. The N-terminal DNA-binding domain is structurally dynamic while the C-terminal domain showed significant stability and such profile of structural dynamics is most likely to be preserved in the structural orthologs of Mce1R. In addition to this, a cavity has been detected in the C-terminal domain of Mce1R which contains a few conserved residues. Comparison with other FCD family of regulators suggests that most of the conserved residues might be critical for binding to specific ligand. The max pKd value and drug score for the cavity are estimated to be 9.04 and 109 respectively suggesting that the cavity represents a suitable target site for novel anti-tuberculosis drug discovery approaches.

Mce2R/Rv0586 of Mycobacterium tuberculosis is the functional homologue of FadRE. coli

Lipid metabolism is critical to Mycobacterium tuberculosis survival and infection. Unlike Escherichia coli, which has a single FadR, the M. tuberculosis genome encodes five proteins of the FadR sub-family. While the role of E. coli FadR as a regulator of fatty acid metabolism is well known, the definitive functions of M. tuberculosis FadR proteins are still under investigation. An interesting question about the M. tuberculosis FadRs remains open: which one of these proteins is the functional homologue of E. coli FadR? To address this, we have applied two different approaches. The first one was the bioinformatics approach and the second one was the classical molecular genetic approach involving complementation studies. Surprisingly, the results of these two approaches did not agree. Among the five M. tuberculosis FadRs, Rv0494 shared the highest sequence similarity with FadR^{E. coli} and Rv0586 was the second best match. However, only Rv0586, but not Rv0494, could complement E. coli ∆fadR, indicating that Rv0586 is the M. tuberculosis functional homologue of FadR^{E. coli}. Further studies showed that both regulators, Rv0494 and Rv0586, show similar responsiveness to LCFA, and have conserved critical residues for DNA binding. However, analysis of the operator site indicated that the inter-palindromic distance required for DNA binding differs for the two regulators. The differences in the binding site selection helped in the success of Rv0586 binding to fadB upstream over Rv0494 and may have played a critical role in complementing E. coli ∆fadR. Further, for the first time, we report the lipid-responsive nature of Rv0586.

Mycobacterium tuberculosis Rv1152 is a Novel GntR Family Transcriptional Regulator Involved in Intrinsic Vancomycin Resistance and is a Potential Vancomycin Adjuvant Target

Novel factors involved in Mycobacteria antibiotics resistance are crucial for better targets to combat the ever-increasing drug resistant strains. Mycobacterium tuberculosis Rv1152, a novel GntR family transcriptional regulator and a promising vancomycin adjuvant target, was firstly characterized in our study. Overexpression of Rv1152 in Mycobacterium smegmatis decreased bacterial susceptibility to vancomycin. Moreover, a deficiency in MSMEG_5174, an Rv1152 homolog made M. smegmatis more sensitive to vancomycin, which was reverted by complementing the MSMEG_5174 deficiency with Rv1152 of M. tuberculosis. Rv1152 negatively regulated four vancomycin responsive genes, namely genes encoding the ribosome binding protein Hsp, small unit of sulfate adenylyltransferase CysD, L-lysine-epsilon aminotransferase Lat, and protease HtpX. Taken together, Rv1152 controls the expression of genes required for the susceptibility to vancomycin. This is the first report that links the GntR family transcriptional factor with vancomycin susceptibility. Inhibitors of Rv1152 might be ideal vancomycin adjuvants for controlling multi-drug resistant Mycobacterial infections.

A GntR family transcription factor positively regulates mycobacterial isoniazid resistance by controlling the expression of a putative permease

Background

Bacteria use transcriptional regulation to respond to environmental stresses. Specifically, exposure to antibacterial drugs is deemed to be an atypical stress, and altering transcriptional regulation in response to such stress can increase bacterial drug resistance. However, only a few transcription factors that regulate drug resistance have been reported.

Results

In the present study, a GntR family transcription factor, encoded by the MSMEG_0535 (Ms0535) gene, was shown to be an isoniazid (INH) resistance regulator in Mycobacterium smegmatis. When the Ms0535 gene was overexpressed, cells showed a significant increase in INH resistance. First, the interaction between Ms0535 and its own promoter was determined, and a conserved 26-bp palindromic DNA binding motif was identified using electrophoretic mobility shift and DNaseI footprinting assays. Second, quantitative reverse transcription-PCR assays showed that Ms0535 acted as a transcriptional activator, and positively regulated its own expression, as well as that of a permease encoded by the MSMEG_0534 (Ms0534) gene. Similar to the case for the Ms0535 gene, a recombinant Ms0534-overexpressing strain also exhibited increased INH resistance compared with the wild-type strain. Furthermore, we showed that Ms0535 and Ms0534 deletion strains were more sensitive to INH than the wild-type strain. Interestingly, overexpressing Ms0534 in the Ms0535 deletion strain enhanced its INH resistance. In contrast, the Ms0534 deletion strain was still sensitive to INH even when Ms0535 was overexpressed. These findings suggest that Ms0534 is an effector protein that affects INH resistance in M. smegmatis.

Conclusions

In summary, the GntR transcriptional regulator Ms0535 positively regulates INH resistance by transcriptionally regulating the expression of the Ms0534 permease in M. smegmatis. These results improve our understanding of the role of transcriptional regulation in INH drug resistance in mycobacteria.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0556-8) contains supplementary material, which is available to authorized users.

Transcriptional regulators of GntR family in Streptomyces coelicolor A3(2): analysis in silico and in vivo of YtrA subfamily

Transcriptional factors of the GntR family regulate numerous physiological and morphological processes in response to the nutrient state of bacterial cells. The number of GntR transcriptional factors in genomes of soil-dwelling actinomycetes is one of the highest among bacteria, reflecting both the large size of their chromosomes and the complex ecological niche that they occupy. However, very little is known about the roles of GntRs in actinomycete biology. Here, we analyzed the genome of model actinomycete, Streptomyces coelicolor A3(2), in an attempt to gain new insights into the function of GntR family. All 56 GntR proteins of M145 strain were classified into FadR, HutC, MocR, YtrA, and DevA subfamilies according to their secondary structure. We then checked for the presence of GntR orthologs in six other sequenced Streptomyces and one Kitasatospora genomes, revealing that 12 GntRs were conserved in all analyzed strains. Genomic analysis of the less studied YtrA type regulators revealed 160 sequences present in 88 members of Coriobacteridae, Rubrobacteridae, and Actinobacteridae subclasses. These proteins form seven dense clusters on the consensus phylogenetic tree and their genes are usually co-located with the genes for transport proteins. Probable operator sites were identified for orthologous groups of Sco0823 and Sco3812 proteins. All S. coelicolor YtrA-like regulatory genes (SCO0823, SCO1728, SCO3812) were analyzed at transcriptional level, knocked out, and introduced on moderate copy number plasmid in M145 strain. Also, gene SCO0824, a part of putative SCO0823 operon, was studied. Results of these experiments are discussed here.

Rv0494 is a starvation-inducible, auto-regulatory FadR-like regulator from Mycobacterium tuberculosis

Fatty acid metabolism plays an important role in the survival and pathogenesis of Mycobacterium tuberculosis. Lipids are assumed to be the major source of energy during dormancy. Here, we report the characterization of a starvation-inducible, lipid-responsive transcriptional regulator, Rv0494, divergently transcribed from the Rv0493c probable operon. The striking difference in the transcriptional regulatory apparatus between mycobacteria and other well-studied organisms, such as Escherichia coli, is the organization of mycobacterial promoters. Mycobacterial promoters have diverse architectures and most of these promoters function inefficiently in E. coli. In this study, we characterized the promoter elements of Rv0494 along with the sigma factors required for transcription initiation. Rv0494 promoter activity increased under nutrient starvation conditions and was transcribed via two promoters: the promoter proximal to the translational start site was active under standard growth conditions, whilst both promoters contributed to the increased activity seen during starvation, with the major contribution from the distal promoter. Furthermore, Rv0494 translation initiated at a codon located 9 bp downstream of the annotated start codon. Rv0494 bound to its upstream sequence to auto-regulate its own expression; this binding was responsive to long-chain fatty acyl-CoA molecules. We further report Rv0494-mediated transcriptional regulation of the Rv2326c gene - a probable transmembrane ATP-binding transporter encoding gene.

RNA sequencing reveals differences between the global transcriptomes of Salmonella enterica serovar enteritidis strains with high and low pathogenicities

Salmonella enterica serovar Enteritidis is one of the important causes of bacterial food-borne gastroenteritis worldwide. Field strains of S. Enteritidis are relatively genetically homogeneous; however, they show extensive phenotypic diversity and differences in virulence potential. RNA sequencing (RNA-Seq) was used to characterize differences in the global transcriptome between several genetically similar but phenotypically diverse poultry-associated field strains of S. Enteritidis grown in laboratory medium at avian body temperature (42°C). These S. Enteritidis strains were previously characterized as high-pathogenicity (HP; n = 3) and low-pathogenicity (LP; n = 3) strains based on both in vitro and in vivo virulence assays. Using the negative binomial distribution-based statistical tools edgeR and DESeq, 252 genes were identified as differentially expressed in LP strains compared with their expression in the HP strains (P < 0.05). A majority of genes (235, or 93.2%) showed significantly reduced expression, whereas a few genes (17, or 6.8%) showed increased expression in all LP strains compared with HP strains. LP strains showed a unique transcriptional profile that is characterized by significantly reduced expression of several transcriptional regulators and reduced expression of genes involved in virulence (e.g., Salmonella pathogenicity island 1 [SPI-1], SPI-5, and fimbrial and motility genes) and protection against osmotic, oxidative, and other stresses, such as iron-limiting conditions commonly encountered within the host. Several functionally uncharacterized genes also showed reduced expression. This study provides a first concise view of the global transcriptional differences between field strains of S. Enteritidis with various levels of pathogenicity, providing the basis for future functional characterization of several genes with potential roles in virulence or stress regulation of S. Enteritidis.

Identification and characterization of Rv0494: a fatty acid-responsive protein of the GntR/FadR family from Mycobacterium tuberculosis

Escherichia coli FadR, a member of the GntR family of transcription factors, plays dual roles in fatty acid metabolism. FadR-DNA binding is inhibited by fatty acyl-CoAs, and thus FadR acts as a sensor of the fatty acid level in bacteria. We have identified FadR-binding sites in the upstream regions of genes showing altered expression after the disruption of fatty acid biosynthesis in Mycobacterium tuberculosis. A FadR homologue in M. tuberculosis, Rv0494, was identified, which binds to its operator in the upstream region of the kas operon. We have shown that FadRMt (Rv0494) directly binds to long-chain fatty acyl-CoA and that binding quenches the intrinsic fluorescence of the purified protein. FadR-DNA binding can be impaired by long-chain fatty acyl-CoA compounds. Overexpression of Rv0494 in Mycobacterium bovis BCG reduced the basal level expression of kas operon genes, thereby suggesting the repressor nature of this protein in fatty acid synthase II regulation. This is the first report, to the best of our knowledge, of a GntR/FadR family protein acting as a fatty acid-responsive transcriptional regulator in M. tuberculosis, suggesting a possible role for this protein in mycolic acid biosynthesis.

Equilibrium binding and kinetic characterization of putative tetracycline repressor family transcription regulator Fad35R from Mycobacterium tuberculosis

Fatty acids play critical role in the survival and virulence of Mycobacterium tuberculosis (Mtb). Activation of fatty acids by acyl-CoA synthetases (Fad) into fatty acyl-CoA is the first and one of the crucial steps in fatty acid metabolism. Mtb possesses 36 fatty acyl-CoA synthetases, unlike Escherichia coli, which has single enzyme. However, the mechanisms by which the expression of these multiple Fad genes is regulated remain uncharacterized. We characterized the DNA- and ligand-binding properties of a putative tetracycline repressor family regulator, named Fad35R, located upstream of the Fad35 gene and ScoA-citE operon. We identified a palindromic regulatory motif upstream of Fad35 and characterized the binding of Fad35R to this motif. Equilibrium binding studies show that Fad35R binds to this motif with high affinity (K(d) ∼ 0.033 μm) and the specificity of binding was confirmed by an electromobility gel shift assay. Kinetic studies indicate that faster association (k(a,avg) ∼ 5.4 × 10(4) m(-1) · s(-1)) and slower dissociation rates (k(d,avg) ∼ 5.84 × 10(-4) s(-1)) confer higher affinity. The affinity for the promoter is maximum at 300 mm NaCl but decreases rapidly beyond this range. Ligand-binding studies indicate that Fad35R binds specifically to tetracycline and also binds to fatty acid derivatives. The promoter-binding affinity is decreased significantly in the presence of palmityl-CoA, suggesting that Fad35R can sense the levels of activated fatty acids and alter its DNA-binding activity. Our results suggest that Fad35R may be the functional homologue of FadR and controls the expression of genes in a metabolite-dependent manner.

A copper-responsive global repressor regulates expression of diverse membrane-associated transporters and bacterial drug resistance in mycobacteria

Sequencing of entire bacterial genomes has led to the identification of many membrane-associated transporters, including several multidrug resistance transport proteins, in recent years. However, the regulators and signaling pathways involved in the expression of these genes remain largely unknown. In this study, we have identified Ms2173, a GntR/FadR family transcription factor, as a novel global regulator in Mycobacterium smegmatis. Ms2173 was found to specifically recognize a 15-bp palindromic motif and to broadly regulate expression of 292 genes, including 37 genes that encode membrane-associated transport proteins. Copper ions induced Ms2173 to form inactive proteins lacking DNA-binding activity. Ms2173 was shown to function as a repressor of its target genes. Interestingly, we found that the function of Ms2173 was linked to mycobacterial drug resistance. Compared with the substantially enhanced drug resistance in the Ms2173-deleted mutant strain, the strains overexpressing Ms2173 were more sensitive to anti-tuberculosis drugs than the wild-type strain. Additionally, copper ions could partially counteract the in vivo function of Ms2173. We have thus characterized the first mycobacterial GntR/Fad-like transcription factor that functions as a copper ion-responsive global repressor that we have renamed GfcR. These findings further enhance our understanding of membrane-associated transporter regulation and drug resistance in mycobacteria.

Molecular mechanisms underlying the function diversity of transcriptional factor IclR family

The IclR family transcriptional factor is widespread and involves in diverse bacterial physio-pathological events, such as primary and secondary metabolism, virulence, quorum sensing, sporulation. Unlike other transcriptional factors which function as either activators or repressors, IclR can assume both role simutaneously. Its N-terminal domain possesses a helix-turn-helix DNA binding motif which can dimerize or tetramerize to bind target promoters, while the C-terminal domain is for the effector binding. The function of IclR varies with the effectors bound. Escherichia coli transcription factor IclR is the archetype of this family, which regulates the aceBAK operon responsible for the glyoxylate shunt. The sophisticated regulatory mechanisms underlying iclR was largely based on E. coli iclR. Information concerning the pathogen IclR, especially those of Mycobacterium tuberculosis is poor, and is pivotal to our understanding of its biology and development of new effective TB control measures.

GntR family regulators of the pathogen of fish tuberculosis Mycobacterium marinum

Mycobacterium marinum is a slow-growing pathogenic mycobacterium. It was first isolated by Aronson in 1926 from fish, fish mycobacteriosis or called fish tuberculosis is the common causative agent of bacterial disease in many species of freshwater and marine fish. M. marinum can infect wild fish, aquaculture and ornamental fish, and it has a close relative of the causative agent of human tuberculosis, Mycobacterium tuberculosis. The recently sequenced genome of M. marinum has been shown to contain several putative GntR regulators. This family named after gluconate regulator has a helix-turn-helix structure. Characterization of transcription regulators and their network is an important step towards the complete understanding of cellular physiology. The regulator of this family shares a similar and conserved N-terminal DNA-binding domain, but has a highly diverse C-terminal effector-binding and oligomerization domain. According to the heterogeneity, we classify the M. marinum GntR family to four subfamilies: FadR, HutC, MocR, and YtrA, and these regulators are encoded by 8, 3, 1 and 1 genes, respectively. Thus this study extends the annotation of M. marinum GntR family proteins, and can help to understand the pathogenic role of this family in M. marinum and facilitate future drug design against this pathogen.

Functional analysis of an intergenic non-coding sequence within mce1 operon of M.tuberculosis

Background

The mce operons play an important role in the entry of M. tuberculosis into macrophages and non-phagocytic cells. Their non-redundant function as well as complex regulation is implied by the phenotype of mce mutants. Recently, mce1 operon was found to extend over 13 genes, fadD5 (Rv0166) being the first gene of the operon. The presence of a non-coding sequence of 200 base pairs between Rv0166 and Rv0167 is peculiar to mce1 among the four mce operons of M.tuberculosis. We have examined the function of this region.

Results

We predicted putative promoter activity of the 200 base pairs of non-coding, intergenic region between Rv0166 and Rv0167 in silico using MEME software and designate it as intergenic promoter, IGPr. We demonstrate both promoter activity and a putative negative regulatory function of this fragment by reporter assays carried out in the surrogate host M.smegmatis. We find that the repressive elements not only control the native promoter but also repress a heterologous promoter of M.smegmatis. The higher activity of the intergenic promoter in a clinical isolate in comparison with the wild type sequence from M.tuberculosis H37Rv could be correlated with a point mutation within the negative element. We have mapped two transcription start sites for mce1 operon both of which are utilized in M.tuberculosis H37Rv as well as the clinical isolate VPCI591. Our studies show that the promoter activity in the non-coding region is relevant not only in reporter gene expression but also in the expression of mce1 operon in M. tuberculosis cells grown in synthetic medium.

Conclusion

The mce operon of M.tuberculosis H37Rv potentially can be transcribed from two promoters P1 and P2, former mapping upstream of Rv0166 and the latter in the non-coding intergenic region between Rv0166 and Rv0167. The transcription initiation from P1 results in a transcript with Rv0166 while that from P2 will be without it. The sequences between the translation start site of Rv0167 and the promoter P2 have a negative regulatory role, as point mutation within the sequence leads to enhanced activity of P2 as well as a heterologous promoter from M.smegmatis. The mutation detected in the clinical isolate VPCI591 therefore behaves like a gain-of-function mutation.

Mce2R from Mycobacterium tuberculosis represses the expression of the mce2 operon

The mce2 operon is one of the four mce operons present in Mycobacterium tuberculosis that encode exported proteins with a probable role in the virulence mechanisms of this bacterium. In the present study we demonstrated that Rv0586, which encodes a putative GntR-like regulator, is part of the mce2 operon. By using a promoter-lacZ fusion approach and bioinformatics tools, we found that Rv0586 represses the expression of Mce2 proteins and of a putative endonuclease IV, encoded by end (Rv0670) gene. For this reason, we have re-named the repressor protein Mce2R. By gel-shift experiments Mce2R binding was determined to be located within the mce2 promoter region. In addition, two FadR-like operator motifs were identified within the promoter regions of both the mce2 operon and the end gene. These motifs overlap putative -10 and -35 promoter boxes. M. tuberculosis carrying mce2 and end promoter-lacZ fusions were used to infect J774 macrophage-like cells. Expression of beta-galactosidase was induced after phagocytocis, suggesting that some cellular factor could be a key component of the molecular switch regulation expression of the mce2 operon. In conclusion, these results add novel evidence of the complex regulation of mce operon expression.

Genome bias influences amino acid choices: analysis of amino acid substitution and re-compilation of substitution matrices exclusive to an AT-biased genome

The genomic era has seen a remarkable increase in the number of genomes being sequenced and annotated. Nonetheless, annotation remains a serious challenge for compositionally biased genomes. For the preliminary annotation, popular nucleotide and protein comparison methods such as BLAST are widely employed. These methods make use of matrices to score alignments such as the amino acid substitution matrices. Since a nucleotide bias leads to an overall bias in the amino acid composition of proteins, it is possible that a genome with nucleotide bias may have introduced atypical amino acid substitutions in its proteome. Consequently, standard matrices fail to perform well in sequence analysis of these genomes. To address this issue, we examined the amino acid substitution in the AT-rich genome of Plasmodium falciparum, chosen as a reference and reconstituted a substitution matrix in the genome's context. The matrix was used to generate protein sequence alignments for the parasite proteins that improved across the functional regions. We attribute this to the consistency that may have been achieved amid the target and background frequencies calculated exclusively in our study. This study has important implications on annotation of proteins that are of experimental interest but give poor sequence alignments with standard conventional matrices.

Comparative genomic analysis of Mycobacterium avium subspecies obtained from multiple host species

Background

Mycobacterium avium (M. avium) subspecies vary widely in both pathogenicity and host specificity, but the genetic features contributing to this diversity remain unclear.

Results

A comparative genomic approach was used to identify large sequence polymorphisms among M. avium subspecies obtained from a variety of host animals. DNA microarrays were used as a platform for comparing mycobacterial isolates with the sequenced bovine isolate M. avium subsp. paratuberculosis (MAP) K-10. Open reading frames (ORFs) were classified as present or divergent based on the relative fluorescent intensities of the experimental samples compared to MAP K-10 DNA. Multiple large polymorphic regions were found in the genomes of MAP isolates obtained from sheep. One of these clusters encodes glycopeptidolipid biosynthesis enzymes which have not previously been identified in MAP. M. avium subsp. silvaticum isolates were observed to have a hybridization profile very similar to yet distinguishable from M. avium subsp. avium. Isolates obtained from cattle (n = 5), birds (n = 4), goats (n = 3), bison (n = 3), and humans (n = 9) were indistinguishable from cattle isolate MAP K-10.

Conclusion

Genome diversity in M. avium subspecies appears to be mediated by large sequence polymorphisms that are commonly associated with mobile genetic elements. Subspecies and host adapted isolates of M. avium were distinguishable by the presence or absence of specific polymorphisms.

GntR family of regulators in Mycobacterium smegmatis: a sequence and structure based characterization

Background

Mycobacterium smegmatis is fast growing non-pathogenic mycobacteria. This organism has been widely used as a model organism to study the biology of other virulent and extremely slow growing species like Mycobacterium tuberculosis. Based on the homology of the N-terminal DNA binding domain, the recently sequenced genome of M. smegmatis has been shown to possess several putative GntR regulators. A striking characteristic feature of this family of regulators is that they possess a conserved N-terminal DNA binding domain and a diverse C-terminal domain involved in the effector binding and/or oligomerization. Since the physiological role of these regulators is critically dependent upon effector binding and operator sites, we have analysed and classified these regulators into their specific subfamilies and identified their potential binding sites.

Results

The sequence analysis of M. smegmatis putative GntRs has revealed that FadR, HutC, MocR and the YtrA-like regulators are encoded by 45, 8, 8 and 1 genes respectively. Further out of 45 FadR-like regulators, 19 were classified into the FadR group and 26 into the VanR group. All these proteins showed similar secondary structural elements specific to their respective subfamilies except MSMEG_3959, which showed additional secondary structural elements. Using the reciprocal BLAST searches, we further identified the orthologs of these regulators in Bacillus subtilis and other mycobacteria. Since the expression of many regulators is auto-regulatory, we have identified potential operator sites for a number of these GntR regulators by analyzing the upstream sequences.

Conclusion

This study helps in extending the annotation of M. smegmatis GntR proteins. It identifies the GntR regulators of M. smegmatis that could serve as a model for studying orthologous regulators from virulent as well as other saprophytic mycobacteria. This study also sheds some light on the nucleotide preferences in the target-motifs of GntRs thus providing important leads for initiating the experimental characterization of these proteins, construction of the gene regulatory network for these regulators and an understanding of the influence of these proteins on the physiology of the mycobacteria.