Bacterioferritin from Mycobacterium smegmatis contains zinc in its di-nuclear site

Bacterioferritins, also known as cytochrome b (1), are oligomeric iron-storage proteins consisting of 24 identical amino acid chains, which form spherical particles consisting of 24 subunits and exhibiting 432 point-group symmetry. They contain one haem b molecule at the interface between two subunits and a di-nuclear metal binding center. The X-ray structure of bacterioferritin from Mycobacterium smegmatis (Ms-Bfr) was determined to a resolution of 2.7 A in the monoclinic space group C2. The asymmetric unit of the crystals contains 12 protein molecules: five dimers and two half-dimers located along the crystallographic twofold axis. Unexpectedly, the di-nuclear metal binding center contains zinc ions instead of the typically observed iron ions in other bacterioferritins.     

Identification of small RNAs in Mycobacterium smegmatis using heterologous Hfq

Gene regulation by small RNAs (sRNAs) has been extensively studied in various bacteria. However, the presence and roles of sRNAs in mycobacteria remain largely unclear. Immunoprecipitation of RNA chaperone Hfq to enrich for sRNAs is one of the effective methods to isolate sRNAs. However, the lack of an identified mycobacterial hfq restricts the feasibility of this approach. We developed a novel method that takes advantage of the conserved inherent sRNAs-binding capability of heterologous Hfq from Escherichia coli to enrich sRNAs from Mycobacterium smegmatis, a model organism for studying Mycobacterium tuberculosis. We validated 12 trans-encoded and 12 cis-encoded novel sRNAs in M. smegmatis. Many of these sRNAs are differentially expressed at exponential phase compared with stationary phase, suggesting that sRNAs are involved in the growth of mycobacteria. Intriguingly, five of the cis-encoded novel sRNAs target known transposases. Phylogenetic conservation analysis shows that these sRNAs are pathogenicity dependent. We believe that our findings will serve as an important reference for future analysis of sRNAs regulation in mycobacteria and will contribute significantly to the development of sRNAs prediction programs. Moreover, this novel method of using heterologous Hfq for sRNAs enrichment can be of general use for the discovery of bacterial sRNAs in which no endogenous Hfq is identified.     

Development of a method to investigate the hydrolysis of xenobiotic esters by a Mycobacterium smegmatis homogenate

One of the main problems in combating tuberculosis is caused by a poor penetration of drugs into the mycobacterial cells. A prodrug approach via activation inside mycobacterial cells is a possible strategy to overcome this hurdle and achieve efficient drug uptake. Esters are attractive candidates for such a strategy and we and others communicated previously the activity of esters of weak organic acids against mycobacteria. However very little is known about ester hydrolysis by mycobacteria and no biological model is available to study the activation of prodrugs by these microorganisms. To begin filling this gap, we have embarked in a project to develop an in vitro method to study prodrug activation by mycobacteria using Mycobacterium smegmatis homogenates. Model ester substrates were ethyl nicotinate and ethyl benzoate whose hydrolysis was monitored and characterized kinetically. Our studies showed that in M. smegmatis most esterase activity is associated with the soluble fraction (cytosol) and is preserved by storage at 5 °C or at room temperature for one hour, or by storage at − 80 °C up to one year. In the range of homogenate concentrations studied (5-80% in buffer), k_obs varied linearly with homogenate concentration for both substrates. We also found that the homogenates showed Michaelis-Menten kinetics behavior with both prodrugs. Since ethyl benzoate is a good substrate for the mycobacterial esterases, this compound can be used to standardize the esterasic activity of homogenates, allowing results of incubations of prodrugs with homogenates from different batches to be readily compared.     

Evaluation of the role of sigma B in Mycobacterium smegmatis

The alternate sigma factor, sigB, is known to play a crucial role in maintaining the stationary phase in mycobacteria. In this communication, we have studied the proteomics of Mycobacterium smegmatis mc(2)155 and its two derivatives, one of which has a disrupted sigB gene and the other, PMVSigB, which contains a multicopy plasmid containing sigB. We have identified by two-dimensional gel analyses, several proteins that are over-expressed in PMVSigB compared to mc(2)155. These proteins are either stress proteins or participate actively in different metabolic pathways of the organisms. On the other hand, when sigB deleted mycobacteria were grown until the stationary phase and its two-dimensional protein profile was compared to that of mc(2)155, few DNA binding proteins were found to be up-regulated. We have shown recently that upon over-expressing sigB, the cell surface glycopeptidolipids of M. smegmatis are hyperglycosylated, a situation similar to what was observed for nutritionally starved bacteria. Gene expression profile through quantitative PCR presented here identified a Rhamnosyltransferase responsible for this hyperglycosylation.     

Biosynthesis of d-arabinose in Mycobacterium smegmatis: specific labeling from d-glucose

d-Arabinose is a major sugar in the cell wall polysaccharides of Mycobacterium tuberculosis and other mycobacterial species. The reactions involved in the biosynthesis and activation of d-arabinose represent excellent potential sites for drug intervention since d-arabinose is not found in mammalian cells, and the cell wall arabinomannan and/or arabinogalactan appear to be essential for cell survival. Since the pathway involved in conversion of d-glucose to d-arabinose is unknown, we incubated cells of Mycobacterium smegmatis individually with [1-(14)C]glucose, [3,4-(14)C]glucose, and [6-(14)C]glucose and compared the specific activities of the cell wall-bound arabinose. Although the specific activity of the arabinose was about 25% lower with [6-(14)C]glucose than with other labels, there did not appear to be selective loss of either carbon 1 or carbon 6, suggesting that arabinose was not formed by loss of carbon 1 of glucose via the oxidative step of the pentose phosphate pathway, or by loss of carbon 6 in the uronic acid pathway. Similar labeling patterns were observed with ribose isolated from the nucleic acid fraction. Since these results suggested an unusual pathway of pentose formation, labeling studies were also done with [1-(13)C]glucose, [2-(13)C]glucose, and [6-(13)C]glucose and the cell wall arabinose was examined by NMR analysis. This method allows one to determine the relative (13)C content in each carbon of the arabinose. The labeling patterns suggested that the most likely pathway was condensation of carbons 1 and 2 of fructose 6-phosphate produced by the transaldolase reaction with carbons 4, 5, and 6 (i.e., glyceraldehyde 3-phosphate) formed by fructose-1,6 bisphosphate aldolase. Cell-free enzyme extracts of M. smegmatis were incubated with ribose 5-phosphate, xylulose 5-phosphate, and d-arabinose 5-phosphate under a variety of experimental conditions. Although the ribose 5-phosphate and xylulose 5-phosphate were converted to other pentoses and hexoses, no arabinose 5-phosphate (or free arabinose) was detected in any of these reactions. In addition, these enzyme extracts did not convert arabinose 5-phosphate to any other pentose or hexose. In addition, incubation of [(14)C]glucose 6-phosphate and various nucleoside triphosphates (ATP, CTP, GTP, TTP, and UTP) with cytosolic or membrane fractions from the mycobacterial cells did not result in formation of a nucleotide form of arabinose, although other radioactive sugars including rhamnose and galactose were found in the nucleotide fraction. Furthermore, no radioactive arabinose was found in the nucleotide fraction isolated from M. smegmatis cells grown in [(3)H]glucose, nor was arabinose detected in a large-scale extraction of the sugar nucleotide fraction from 300 g of cells. The logical conclusion from these studies is that d-arabinose is probably produced from d-ribose by epimerization of carbon 2 of the ribose moiety of polyprenylphosphate-ribose to form polyprenylphosphate-arabinose, which is then used as the precursor for formation of arabinosyl polymers.     

Metabolism of cyclic AMP in non-pathogenic Mycobacterium smegmatis

The intracellular concentration of cyclic AMP reached a maximum in 3.5-day old cultures of Mycobacterium smegmatis grown in the presence of glycerol as the main source of carbon. Glucose-grown cells exhibited decreased cyclic AMP levels at all stages of growth. When M. smegmatis cells were incubated with various metabolites, pyruvate increased whereas glucose, citric acid, succinic acid and lactic acid decreased intracellular cyclic AMP levels. No cyclic AMP was detected in the incubation medium. The presence of a cyclic AMP-binding protein was demonstrated in cellfree extracts of M. smegmatis.     

Mycobacterium tuberculosis ESAT-6 Exhibits a Unique Membrane-interacting Activity That Is Not Found in Its Ortholog from Non-pathogenic Mycobacterium smegmatis

Mycobacterium tuberculosis ESAT-6 (MtbESAT-6) reportedly shows membrane/cell-lysis activity, and recently its biological roles in pathogenesis have been implicated in rupture of the phagosomes for bacterial cytosolic translocation. However, molecular mechanism of MtbESAT-6-mediated membrane interaction, particularly in relation with its biological functions in pathogenesis, is poorly understood. In this study, we investigated the pH-dependent membrane interaction of MtbESAT-6, MtbCFP-10, and the MtbESAT-6/CFP-10 heterodimer, by using liposomal model membranes that mimic phagosomal compartments. MtbESAT-6, but neither MtbCFP-10 nor the heterodimer, interacted with the liposomal membranes at acidic conditions, which was evidenced by release of K⁺ ions from the liposomes. Most importantly, the orthologous ESAT-6 from non-pathogenic Mycobacterium smegmatis (MsESAT-6) was essentially inactive in release of K⁺. The differential membrane interactions between MtbESAT-6 and MsESAT-6 were further confirmed in an independent membrane leakage assay using the dye/quencher pair, 8-aminonapthalene-1,3,6 trisulfonic acid (ANTS)/p-xylene-bis-pyridinium bromide (DPX). Finally, using intrinsic and extrinsic fluorescence approaches, we probed the pH-dependent conformational changes of MtbESAT-6 and MsESAT-6. At acidic pH conditions, MtbESAT-6 underwent a significant conformational change, which was featured by an increased solvent-exposed hydrophobicity, while MsESAT-6 showed little conformational change in response to acidification. In conclusion, we have demonstrated that MtbESAT-6 possesses a unique membrane-interacting activity that is not found in MsESAT-6 and established the utility of rigorous biochemical approaches in dissecting the virulence of M. tuberculosis.     

Effect of morphine on <Emphasis Type="Italic">Mycobacterium smegmatis</Emphasis> infection in mice and macrophages

The immunomodulatory effects of opioids are known in various infections. However, little is known about the effects of opioids in tuberculosis (TB). In the present study, we report the effects of morphine in Mycobacterium smegmatis infection in mice and macrophages. Morphine exerted a dose-dependent suppression of infection in vivo: 50 and 100 mg/kg morphine exerted significant (P<0.05) suppression whereas 5 mg/kg morphine showed no effect. Analogous to the in vivo effects, incubation of M. smegmatis-infected mouse peritoneal macrophages with morphine (100 μM) showed significant reduction in intramacrophage CFU counts. However, morphine did not show any direct antimycobacterial activity in broth dilution assay upto 100 μM concentration. Further, morphine-induced intramacrophage killing of M. smegmatis was abrogated by naloxone and aminoguanidine indicating the involvement of opioid-receptor activation and nitric oxide production in protective effects of morphine. In conclusion, morphine suppressed the progression of experimental TB in both mice and macrophage models.     

Turnover of lipids in Mycobacterium smegmatis CDC 46 and Mycobacterium phlei ATCC 354

The rates of breakdown and renewal of individual lipids in cultures of Mycobacterium smegmatis CDC 46 and Mycobacterium phlei ATCC 354 were investigated by means of a pulse labelling technique using palmitate-1-¹⁴C. The results indicated that in growing cultures of both strains phospholipids were broken down, and cardiolipin had a very rapid turnover. In chase experiments, almost 45% and 40% of the radioactivity of this component were lost respectively from M. smegmatis and M. phlei during one generation time of the cell. The other two major components, phosphatidyl ethanolamine and phosphatidylinositol mannosides showed relatively low turnover. The loss of radioactivity from phosphatidylinositol mannosides was greater in M. phlei than in M. smegmatis but the loss of radioactivity from phosphatidyl ethanolamine was higher in M. smegmatis. The pattern of loss of radioactivity from lipids was almost the same in both strains, the difference being only in the extent of loss. The differences in the cellular localization of the phospholipids indicate their different roles within the cell. Results obtained with the glyceride fraction indicated a very rapid turnover of triglycerides in both strains.Abbreviations CL Cardiolipin - PE Phosphatidyl ethanolamine - PIMx phosphatidylinositol mannosides - PIM₂A phosphatidylinositol dimannoside tetra acylated - PIM₂B phosphatidylinositol dimannoside tri acylated - PIM₅ phosphatidylinositol pentamannoside tetra acylated     

Snapshots of RecA protein involving movement of the C-domain and different conformations of the DNA-binding loops: crystallographic and comparative analysis of 11 structures of Mycobacterium smegmatis RecA

Mycobacterium smegmatis RecA and its nucleotide complexes crystallize in three different, but closely related, forms characterized by specific ranges of unit cell dimensions. The six crystals reported here and five reported earlier, all grown under the same or very similar conditions, belong to these three forms, all in space group P6(1). They include one obtained by reducing relative humidity around the crystal. In all crystals, RecA monomers form filaments around a 6(1) screw axis. Thus, the c-dimension of the crystal corresponds to the pitch of the RecA filament. As reported for Escherichia coli RecA, the variation in the pitch among the three forms correlates well with the motion of the C-terminal domain of the RecA monomers with respect to the main domain. The domain motion is compatible with formation of inactive as well as active RecA filaments involving monomers with a fully ordered C domain. It does not appear to influence the movement upon nucleotide-binding of the switch residue, which is believed to provide the trigger for transmitting the effect of nucleotide binding to the DNA-binding region. Interestingly, partial dehydration of the crystal results in the movement of the residue similar to that caused by nucleotide binding. The ordering of the DNA-binding loops, which present ensembles of conformations, is also unaffected by domain motion. The conformation of loop L2 appears to depend upon nucleotide binding, presumably on account of the movement of the switch residue that forms part of the loop. The conformations of loops L1 and L2 are correlated and have implications for intermolecular communications within the RecA filament. The structures resulting from different orientations of the C domain and different conformations of the DNA-binding loops appear to represent snapshots of the RecA at different phases of activity, and provide insights into the mechanism of action of RecA.     

Oxygen binding and NO scavenging properties of truncated hemoglobin, HbN, of Mycobacterium smegmatis

Unraveling of microbial genome data has indicated that two distantly related truncated hemoglobins (trHbs), HbN and HbO, might occur in many species of slow-growing pathogenic mycobacteria. Involvement of HbN in bacterial defense against NO toxicity and nitrosative stress has been proposed. A gene, encoding a putative HbN homolog with conserved features of typical trHbs, has been identified within the genome sequence of fast-growing mycobacterium, Mycobacterium smegmatis. Sequence analysis of M. smegmatis HbN indicated that it is relatively smaller in size and lacks N-terminal pre-A region, carrying 12-residue polar sequence motif that is present in HbN of M. tuberculosis. HbN encoding gene of M. smegmatis was expressed in E. coli as a 12.8kD homodimeric heme protein that binds oxygen reversibly with high affinity (P50 approximately 0.081 mm Hg) and autooxidizes faster than M. tuberculosis HbN. The circular dichroism spectra indicate that HbN of M. smegmatis and M. tuberculosis are structurally similar. Interestingly, an hmp mutant of E. coli, unable to metabolize nitric oxide, exhibited very low NO uptake activity in the presence of M. smegmatis HbN as compared to HbN of M. tuberculosis. On the basis of cellular heme content, specific nitric oxide dioxygenase (NOD) activity of M. smegmatis HbN was nearly one-third of that from M. tuberculosis. Additionally, the hmp mutant of E. coli, carrying M. smegmatis HbN, exhibited nearly 10-fold lower cell survival under nitrosative stress and nitrite derived reactive nitrogen species as compared to the isogenic strain harboring HbN of M. tuberculosis. Taken together, these results suggest that NO metabolizing activity and protection provided by M. smegmatis HbN against toxicity of NO and reactive nitrogen is significantly lower than HbN of M. tuberculosis. The lower efficiency of M. smegmatis HbN for NO detoxification as compared to M. tuberculosis HbN might be related to different level of NO exposure and nitrosative stress faced by these mycobacteria during their cellular metabolism.     

Phospholipids of ethambutol-susceptible and resistant strains ofMycobacterium smegmatis

The composition, subcellular distribution and rate of synthesis of phospholipids were compared in ethambutol susceptible and resistant strains ofMycobacterium smegmatis. Significant quantitative alterations in phospholipids accompanied the acquisition of resistance, whereas fatty acyl group composition of total phospholipid remained the same in ethambutol resistant and susceptible strains. Cell wall of resistant strain exhibited an accumulation of phospholipids and a decrease in the degree of unsaturation of phospholipid fatty acyl groups. Changes in the cell wall phospholipid composition may contribute to resistance ofMycobacterium smegmatis to ethambutol.     

An orphan gyrB in the Mycobacterium smegmatis genome uncovered by comparative genomics

DNA gyrase is an essential topoisomerase found in all bacteria. It is encoded bygyrB andgyrA genes. These genes are organized differently in different bacteria. Direct comparison ofMycobacterium tuberculosis andMycobacterium smegmatis genomes reveals presence of an additionalgyrB inM. smegmatis flanked by novel genes. Analysis of the amino acid sequence of GyrB from different organisms suggests that the orphan GyrB inM. smegmatis may have an important cellular role.     

Presence of mycobacteriophage I3-like DNA sequences in the genome of its host Mycobacterium smegmatis

The genomes of phage I3 and its host Mycobacterium smegmatis have been compared. From thermal melting studies the GC contents of DNA from mycobacteriophage I3 and its host M. smegmatis were found to be 66%. A new method, based only on the initial rates of reassociation, has been developed for calculating the DNA homology. Analysis of DNA reassociation kinetics suggested the presence of one equivalent of the phage I3 genome within the M. smegmatis genome. Southern analysis revealed the presence of almost all of the phage I3 specific sequences within the host genome.     

Origin of asymmetry in adenylyl cyclases: structures of Mycobacterium tuberculosis Rv1900c

Rv1900c, a Mycobacterium tuberculosis adenylyl cyclase, is composed of an N-terminal alpha/beta-hydrolase domain and a C-terminal cyclase homology domain. It has an unusual 7% guanylyl cyclase side-activity. A canonical substrate-defining lysine and a catalytic asparagine indispensable for mammalian adenylyl cyclase activity correspond to N342 and H402 in Rv1900c. Mutagenic analysis indicates that these residues are dispensable for activity of Rv1900c. Structures of the cyclase homology domain, solved to 2.4 A both with and without an ATP analog, form isologous, but asymmetric homodimers. The noncanonical N342 and H402 do not interact with the substrate. Subunits of the unliganded open dimer move substantially upon binding substrate, forming a closed dimer similar to the mammalian cyclase heterodimers, in which one interfacial active site is occupied and the quasi-dyad-related active site is occluded. This asymmetry indicates that both active sites cannot simultaneously be catalytically active. Such a mechanism of half-of-sites-reactivity suggests that mammalian heterodimeric adenylyl cyclases may have evolved from gene duplication of a primitive prokaryote-type cyclase, followed by loss of function in one active site.     

Functional insights by comparison of modeled structures of 18kDa small heat shock protein and its mutant in Mycobacterium leprae

In this work we are proposing Homology modeled structures of Mycobacterium leprae 18kDa heat shock protein and its mutant. The more closely related structure of the small heat shock protein (sHSP) belonging to the eukaryotic species from wheat sHSP16.9 and 16.3kDa ACR1 protein from Mycobacterium tuberculosis were used as template structures. Each model contains an N-terminal domain, alpha-crystalline domain and a C-terminal tail. The models showed that a single point mutation from serine to proline at 52^nd position causes structural changes. The structural changes are observed in N-terminal region and alpha-crystalline domains. Serine in 52^nd position is observed in β4 strand and Proline in 52^nd position is observed in loop. The number of residues contributing α helix at N-terminal region varies in both models. In 18S more number of residues is present in α helix when compared to 18P. The loop regions between β3 and β4 strands of both models vary in number of residues present in it. Number of residues contributing β4 strand in both models vary. β6 strand is absent in both models. Major functional peptide region of alpha crystalline domains of both models varies. These differences observed in secondary structures support their distinct functional roles. It also emphasizes that a point mutation can cause structural variation.     

Superoxide dismutase from Mycobacterium species,strain Takeo

Superoxide dismutase from Mycobacterium species, strain Takeo, has been purified to homogeneity as judged by disc gel electrophoresis and ultracentrifugation. The enzyme was found to have a molecular weight of approximately 61 500 by sedimentation equilibrium and to contain manganese by atomic absorption and electron spin resonance spectra. The amino acid composition was also determined. The enzyme was considerably stable to the treatment with sodium dodecyl sulfate; unless incubating at 80°C for 2 min, it was not completely dissociated into the subunits. The molecular weight of the subunit was found to be approximately 21 000. Antibodies against the superoxide dismutase were produced by immunization of rabbits with the enzyme, and the -globulin fraction was purified. Superoxide dismutase preparations obtained from various species of mycobacteria and nocardia cross-reacted to different degrees with these antibodies on the Ouchterlony double diffusion plates. Comparative immunological studies indicated that strain Takeo might be most closely related to Myobacterium smegmatis among species of mycobacteria and nocardia tested. The antibodies against superoxide dismutase may be used as a valuable tool for the classification of mycobacteria.     

Binding of the anti-tubercular drug isoniazid to the arylamine N-acetyltransferase protein from Mycobacterium smegmatis

Isoniazid is a frontline drug used in the treatment of tuberculosis (TB). Isoniazid is a prodrug, requiring activation in the mycobacterial cell by the catalase/peroxidase activity of the katG gene product. TB kills two million people every year and the situation is getting worse due to the increase in prevalence of HIV/AIDS and emergence of multidrug-resistant strains of TB. Arylamine N-acetyltransferase (NAT) is a drug-metabolizing enzyme (E.C. 2.1.3.5). NAT can acetylate isoniazid, transferring an acetyl group from acetyl coenzyme A onto the terminal nitrogen of the drug, which in its N-acetylated form is therapeutically inactive. The bacterium responsible for TB, Mycobacterium tuberculosis, contains and expresses the gene encoding the NAT protein. Isoniazid binds to the NAT protein from Salmonella typhimurium and we report here the mode of binding of isoniazid in the NAT enzyme from Mycobacterium smegmatis, closely related to the M. tuberculosis and S. typhimurium NAT enzymes. The mode of binding of isoniazid to M. smegmatis NAT has been determined using data collected from two distinct crystal forms. We can say with confidence that the observed mode of binding of isoniazid is not an artifact of the crystallization conditions used. The NAT enzyme is active in mycobacterial cells and we propose that isoniazid binds to the NAT enzyme in these cells. NAT activity in M. tuberculosis is likely therefore to modulate the degree of activation of isoniazid by other enzymes within the mycobacterial cell. The structure of NAT with isoniazid bound will facilitate rational drug design for anti-tubercular therapy.     

Genome analysis reveals three genomospecies in Mycobacterium abscessus

Background

Mycobacterium abscessus complex, the third most frequent mycobacterial complex responsible for community- and health care-associated infections in developed countries, comprises of M. abscessus subsp. abscessus and M. abscessus subsp. bolletii reviously referred as Mycobacterium bolletii and Mycobacterium massiliense. The diversity of this group of opportunistic pathogens is poorly described.

Results

In-depth analysis of 14 published M. abscessus complex genomes found a pan-genome of 6,153 proteins and core-genome of 3,947 (64.1%) proteins, indicating a non-conservative genome. Analysing the average percentage of amino-acid sequence identity (from 94.19% to 98.58%) discriminates three main clusters C1, C2 and C3: C1 comprises strains belonging to M. abscessus, C2 comprises strains belonging to M. massiliense and C3 comprises strains belonging to M. bolletii; and two sub-clusters in clusters C2 and C3. The phylogenomic network confirms these three clusters. The genome length (from 4.8 to 5.51-Mb) varies from 5.07-Mb in C1, 4.89-Mb in C2A, 5.01-Mb in C2B and 5.28-Mb in C3. The mean number of prophage regions (from 0 to 7) is 2 in C1; 1.33 in C2A; 3.5 in C2B and five in C3. A total of 36 genes are uniquely present in C1, 15 in C2 and 15 in C3. These genes could be used for the detection and identification of organisms in each cluster. Further, the mean number of host-interaction factors (including PE, PPE, LpqH, MCE, Yrbe and type VII secretion system ESX3 and ESX4) varies from 70 in cluster C1, 80 in cluster C2A, 74 in cluster C2B and 93 in clusters C3A and C3B. No significant differences in antibiotic resistance genes were observed between clusters, in contrast to previously reported in-vitro patterns of drug resistance. They encode both penicillin-binding proteins targeted by β-lactam antibiotics and an Ambler class A β-lactamase for which inhibitors exist.

Conclusions

Our comparative analysis indicates that M. abscessus complex comprises three genomospecies, corresponding to M. abscessus, M. bolletii, and M. massiliense. The genomics data here reported indicate differences in virulence of medical interest; and suggest targets for the refined detection and identification of M. abscessus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-359) contains supplementary material, which is available to authorized users.     

Analysis of DevR regulated genes in Mycobacterium tuberculosis

The DevRS two component system of Mycobacterium tuberculosis is responsible for its dormancy in host and becomes operative under hypoxic condition. It is experimentally known that phosphorylated DevR controls the expression of several downstream genes in a complex manner. In the present work we propose a theoretical model to show role of binding sites in DevR mediated gene expression. Individual and collective role of binding sites in regulating DevR mediated gene expression has been shown via modeling. Objective of the present work is twofold. First, to describe qualitatively the temporal dynamics of wild type genes and their known mutants. Based on these results we propose that DevR controlled gene expression follows a specific pattern which is efficient in describing other DevR mediated gene expression. Second, to analyze behavior of the system from information theoretical point of view. Using the tools of information theory we have calculated molecular efficiency of the system and have shown that it is close to the maximum limit of isothermal efficiency.