Unraveling the role of the transcriptional regulator VirS in low pH–induced responses of Mycobacterium tuberculosis and identification of VirS inhibitors

The ability of Mycobacterium tuberculosis to respond and adapt to various stresses such as oxygen/nitrogen radicals and low pH inside macrophages is critical for the persistence of this human pathogen inside its host. We have previously shown that an AraC/XylS-type transcriptional regulator, VirS, which is induced in low pH, is involved in remodeling the architecture of the bacterial cell envelope. However, how VirS influences gene expression to coordinate these pH responses remains unclear. Here, using a genetic biosensor of cytoplasmic pH, we demonstrate that VirS is required for the intracellular pH maintenance in response to acidic stress and inside acidified macrophages. Furthermore, we observed that VirS plays an important role in blocking phagosomal-lysosomal fusions. Transcriptomics experiments revealed that VirS affects the expression of genes encoding metabolic enzymes, cell-wall envelope proteins, efflux pumps, ion transporters, detoxification enzymes, and transcriptional regulators expressed under low-pH stress. Employing electrophoretic mobility-shift assays, DNA footprinting, and in silico analysis, we identified a DNA sequence to which VirS binds and key residues in VirS required for its interaction with DNA. A significant role of VirS in M. tuberculosis survival in adverse conditions suggested it as a potential anti-mycobacterial drug target. To that end, we identified VirS inhibitors in a virtual screen; the top hit compounds inhibited its DNA-binding activity and also M. tuberculosis growth in vitro and inside macrophages. Our findings establish that VirS mediates M. tuberculosis responses to acidic stress and identify VirS-inhibiting compounds that may form the basis for developing more effective anti-mycobacterial agents.

Mycobacterium tuberculosis is known to resist the acidic stress encountered in macrophages and multiply in these hostile conditions; however, the mechanisms for its survival in acidic conditions are poorly understood. There are a few genes that have been implicated in acid resistance in M. tuberculosis, which include serine protease Rv3671c, the OmpATb operon, and a putative magnesium transporter, MgtC (1)(2)(3). A transposon mutant of Rv3671c was shown to be impaired in the maintenance of intrabacterial pH under acidic conditions, suggesting its involvement in the acid resistance (1). OmpATb, an acid-responsive porin, was shown to be involved in the adaptation to acidic environment by mediating secretion of ammonia; however, its role was considered to be redundant, as knockout of OmpATb did not influence the virulence of the pathogen in mice (2,4). M. tuberculosis lacking MgtC, a putative magnesium transporter, was found to be attenuated for growth under mild acidic conditions at low Mg 2ϩ (3). Apart from M. tuberculosis, acid resistance has also been shown to be important for other bacteria, such as Helicobacter pylori, which colonizes in the human stomach, having an acidic environment; Streptococcus pneumoniae; and pathogenic strains of Escherichia coli and Salmonella enterica. Studies have shown the involvement of a few proteins in the acid resistance of these bacteria that include urease and ExbD in H. pylori, Gad proteins, and F 0 F 1 -ATPase in E. coli and Mg 2ϩ transporter in Salmonella (5)(6)(7)(8)(9).
VirS (Rv3082c) of M. tuberculosis belongs to the AraC family of transcriptional regulators (10,11). VirS is present divergently upstream of an acid-inducible operon termed the mymA operon, which comprises of seven genes (Rv3083-Rv3089) (12). The transcription of the mymA operon under acidic stress has been earlier shown to be regulated by VirS, which itself is regulated by acidic pH (12). Studies demonstrated that the virS mutant of M. tuberculosis exhibited altered cell-wall structure, altered mycolic acid content, defective intramacrophage survival, and reduced hematogenous dissemination in vivo (13). Importantly, virS expression was induced during chronic and reactivation phases of murine tuberculosis, implicating VirS in persistence and reactivation of tuberculosis (14). Despite these findings, mechanisms of how VirS exerts its influence on gene expression to elicit the response of M. tuberculosis under acid stress remain uncharacterized.
Here, our study has delineated the contribution of VirS in acid stress and how it mediates its influence on gene expression to coordinate pH responses in M. tuberculosis. Our study shows the requirement of VirS in controlling cytoplasmic pH balance and phagosomal maturation during acidic stress. Further, structure-guided mutational studies revealed key residues of VirS required for its interaction with DNA, and structure-based virtual screening identified some molecules inhibiting the DNA binding activity of VirS that could be a subject of lead optimization for future drug design against M. tuberculosis.

VirS is required to survive under acidic stress and to block phagosomal maturation
VirS has been previously reported to be up-regulated under acidic conditions (12). To understand the involvement of VirS in acidic responses of M. tuberculosis, an in vitro growth study was carried out by growing M. tuberculosis Erdman (WT), virS mutant, and virS complemented strain under varying pH conditions (pH 4.5, 5.0, 5.5, and 6.6) in MB7H9 medium, and the survival of the cells under these conditions was evaluated. The growth of all three strains was comparable at pH 6.6 and 5.5. However, the growth of the virS mutant was significantly reduced at pH 5.0 and 4.5, with a pronounced defect at pH 4.5 as compared with parental M. tuberculosis and virS complemented strain (Fig. 1). We also performed survival studies of WT, virS mutant, and virS complemented strain at acidic conditions of pH 4.5 in 7H9 medium containing nonhydrolyzable tyloxapol (7H9 -4.5-Ty) as the dispersing agent instead of Tween 80 to negate the possibility of a growth defect due to hydrolysis of Tween 80 under acidic conditions to free fatty acids, which can be toxic to the cells (1). We monitored survival of these strains after 6 and 9 days of incubation in 7H9 -4.5-Ty medium by cfu enumeration. It was observed that virS mutant displayed a growth defect in acidified medium after 6 and 9 days of incubation, whereas parental and complemented strain did not exhibit any marked effect on their growth under acidic conditions, suggesting the involvement of VirS in the survival of the bacteria under acidic conditions (Fig. S1). Moreover, it was reported earlier that virS mutant showed survival defects specifically in immune-activated macrophages (13). M. tuberculosis is known to survive under acidic conditions in macrophages by arresting the phagosome-lysosome fusion. Hence, we evaluated the role of VirS in arresting phagosomal maturation in M. tuberculosis. For this, we employed FITC-labeled M. tuberculosis, virS mutant, and virS complemented strain and studied the localization of the pathogen in the lysosomal acidic compartments by using LysoTracker Red dye in THP-1 macrophages. In resting macrophages, all three of the strains exhibited similar colocalization of phagosomes with the LysoTracker-rich compartments (data not shown). However, in the case of activated macrophages, parental M. tuberculosis cells exhibited ϳ9% colocalization with the LysoTracker-stained acidified compartments and largely remained in the nonacidified phagosomes, whereas in contrast to the WT strain, virS mutant displayed significantly greater colocalization (ϳ32%) with the acidified compartments (Fig. 2, A and B). The complemented strain also exhibited ϳ6% colocalization, which is similar to that observed in the case of WT. These observations suggest the involvement of VirS in the arrest of phagosomal maturation in the activated macrophages, providing a survival advantage to the pathogen under harsher conditions.

Involvement of VirS in maintaining the intrabacterial pH of M. tuberculosis in acidic conditions
Activated macrophages are known to have an acidic environment, and M. tuberculosis robustly persists in this acidic milieu, indicating the role of mycobacterial factors in intramycobacterial pH homeostasis. Hence, we next sought to determine the role of VirS in maintaining the pH homeostasis by measuring the intrabacterial pH (pH IB ) 4 of various strains by using a pH-sensitive ratiometric GFP at various intervals under normal and acidic conditions. The pH IB was calculated by measuring the ratio of the fluorescence intensities at excitation wavelengths of 395 and 475 nm with an emission wavelength of 510 nm (1). A standard curve of the ratio of the intensities (excitation wavelengths, 395 and 475 nm; emission wavelength, 510 nm) at varying pH conditions was plotted to determine the unknown pH IB values (Fig. S2A). It was observed that under normal buffer conditions (pH 7.4), the pH IB was maintained nearly neutral in the parental (pH IB ϳ7.0) as well as in the complemented strain even after 8 and 24 h of incubation (Fig. 3A). However, the virS mutant consistently maintained marginally lower intrabacterial pH (pH IB 6.75) even in normal buffer conditions as compared with the parental and complemented strains (Fig. 3A). At pH 4.5, the initial pH IB values of all three strains were similar to that observed in normal conditions (Fig.  3B). However, after 8 and 24 h of incubation in pH 4.5 buffer, the intrabacterial pH of the mutant dropped to ϳ6.4, whereas the pH IB values of parental and complemented strains were still maintained at nearly neutral conditions of ϳ6.9 (Fig. 3B). This sug-gested that both the parental and complemented strains maintain pH homeostasis upon acidic stress, whereas the virS mutant was relatively defective. Taken together, these observations suggest the importance of VirS in maintaining the intrabacterial pH under acidic conditions. We next sought to determine the pH IB of mycobacteria residing in phagosomes of IFN-␥-activated macrophages. The unknown pH IB values were determined by using a standard curve of the ratio of intensities at excitation wavelengths 405 and 476 nm, emission wavelength 535 nm at varying pH conditions (Fig. S2B). It was observed that WT and virS complemented strains mostly maintained their pH IB at pH 6.5-7.25 after 24 h postinfection. However, the majority of the virS mutant mostly maintained its pH IB at pH 6.26 -7.0 after 24 h postinfection ( Fig. 3C) with a small fraction of mutant cells having a pH IB in the range of Ͻ5.5-6.2, suggesting that at 24 h postinfection, a few mutant cells were unable to maintain their intrabacterial pH near neutral. This effect was more pronounced after 48 h postinfection, when we found that the pH IB of parental and virS complemented strain mostly had a pH IB of 6.51-7.25, whereas virS mutant mostly had a pH IB of 5.5 or lower in activated macrophages, suggesting that after 48 h of postinfection, the virS mutant strain was unable to maintain its intrabacterial pH (Fig. 3D). This inability of virS mutant to maintain its intrabacterial pH in activated macrophages indicates its role in the maintenance of pH homeostasis upon acidic stress and could be a probable reason for its poor survival under acidic conditions.

Identification of genes regulated by VirS (VirS-regulon) under acidic conditions by microarray studies
To understand the role of VirS in regulating acidic pHmediated gene expression of M. tuberculosis, we carried out microarray studies by isolating RNA from M. tuberculosis Erd-

Role of VirS in acidic stress and its inhibitor identification
man, virS mutant, and virS complemented strain grown in MB7H9 medium at pH 6.6 as well as at pH 4.5. We shortlisted the VirS-dependent genes induced or repressed at acidic pH relative to neutral pH with Ͼ2.0-fold difference in the expression level when compared with M. tuberculosis Erdman and virS mutant. We then further shortlisted these genes on the basis of their similar expression profile at acidic pH in M. tuberculosis Erdman and virS complemented strain compared with virS mutant, suggesting that indeed the disruption of virS was involved in the change of expression of such genes at acidic pH. We also carried out the classification of all of these induced and repressed genes into various classes based on their function and homology; however, the hypothetical proteins, unannotated genes, and genes of unknown function were not considered for further analysis. All of these analyses revealed that VirS induces expression of 54 genes and represses 39 genes (Ͼ2.0-fold difference; p Ͻ 0.05) at acidic pH (Fig. 4, A and B). The major classes of genes (after removing genes belonging to hypothetical proteins, genes of unknown function, and unannotated genes) found to be induced and repressed were related to cellular metabolism followed by cell wall-associated genes involved in fatty acid oxidation and biosynthesis, which is in agreement with the role of VirS in the remodeling of cell-enve-lope architecture as reported earlier (13). The other classes of genes induced by VirS at acidic pH were related to the PE/PPE family, efflux pumps and ion transporters, detoxification, and transcriptional regulation. Additionally, genes related to energy pathways and respiration and genes related to information pathways were also found to be regulated by VirS at acidic pH.

Identification of VirS-binding region of DNA
Based on the above observations, VirS appears to be a promising drug target regulating a number of important downstream targets. Hence, we further attempted to characterize VirS by evaluating its DNA-binding activity, its DNA-binding site, and crucial residues of VirS required for its activity and to identify . virS mutant failed to maintain its intrabacterial pH inside activated macrophages, whereas WT and virS complemented strain were able to maintain their pH IB levels. There were ϳ115 bacterial groups examined, and data represent two independent experiments.

Role of VirS in acidic stress and its inhibitor identification
its inhibitors. Toward this, we cloned the virS gene in the pET43.1a vector as described under "Experimental procedures" for the expression of fusion protein NusA-VirS, because our earlier studies showed VirS to be an insoluble protein when expressed in E. coli, 5 and the NusA tag is considered to be an effective solubility tag. NusA-VirS was indeed found to be localized in the soluble fraction; however, the NusA-VirS fusion protein failed to bind to the resin, probably due to the nonexposure of the His tag. Hence, cell lysates were used for further studies (data not shown). The DNA-binding activity of NusA-VirS fusion was determined by employing electrophoretic mobilityshift assay (EMSA) by using a 325 bp DNA fragment, including the intergenic region between the virS and mymA genes, as VirS regulates its own expression and activates mymA operon (12). The fusion protein displayed a shift in the mobility of DNA suggesting the binding of VirS to the DNA fragment employed in the assay (Fig. 5A). To rule out the possibility of involvement of NusA tag in the observed shift in the DNA, the lysate with overexpressed NusA was employed as a control that did not display any shift in the mobility of DNA, indicating that the  Microarray studies revealed VirS-mediated differential gene regulation at acidic pH. Up-regulation of VirS at acidic pH causes change in the expression of many genes that were found to be induced and repressed at acidic pH. These genes were categorized into various classes based on their description given by Tuberculist. A, pie chart represents the categorization of various genes that are induced by VirS at acidic pH. B, pie chart represents the categorization of various genes that are repressed by VirS at acidic pH. The percentage distribution of genes is represented in both of the pie charts. C, heat map showing the expression of VirS regulated genes in M. tuberculosis Erdman, M.tb⌬virS, and M.tb⌬virSComp in response to acidic pH relative to neutral pH (Ͼ2-fold difference, p Ͻ 0.05). The figure here depicts complete linkage and was generated by employing Heatmapper software (57).

Role of VirS in acidic stress and its inhibitor identification
observed shift in the mobility was due to the binding of VirS to the DNA (Fig. 5A). Further, a DNase I protection assay was performed to determine the VirS-binding region of DNA. The lysate of overexpressed fusion protein (i.e. NusA-VirS) was used for determining the binding region of VirS that displayed a protected region (Fig. 5B). The lysate of only overexpressed NusA protein was used as a negative control, which did not display any protection. It was found that VirS protected a 46-bp sequence that was determined to be 5Ј-CTGTCGGAAAAT-GATAAAAGCGTGTCGCAAAGTGTCAATACGTGGC-3Ј (Fig. 5B). Based on this VirS-binding DNA sequence, we determined a VirS-binding consensus logo by using the software Weblogo (15) taking the upstream regions of a few identified downstream targets of VirS for alignment (16) (Fig. 5C). For this, we performed multiple-sequence alignment of this DNA sequence with the DNA sequence of micF and mar promoters used in the crystal structure determination of Rob and MarA, respectively, which provided us with a 28-bp sequence exhibiting alignment and conservation in DNA-binding regions of these three proteins. This 28-bp region was employed to search for homology in the upstream regions of a few of the VirS-regulated genes, namely icl, ctpG, ahpC, infC, sirA, pdc, Rv1405c, adh, fadD9, and fadE35, identified by microarray. These genes were randomly selected to represent a variety of different functional categories.
We observed that all of these genes carried the 28-bp sequence required for binding of VirS although with a few variations in the sequences (Fig. 5C). We further explored the regulation of a few of the genes identified by microarray studies by VirS. The genes selected were Rv2428 (ahpD), Rv1706 (ppe23), Rv1994c (cmtR), Rv1992c (ctpG), Rv2352c (ppe38), and Rv2590 (fadD9). For determining whether these genes are directly or indirectly regulated by VirS, we selected a ϳ400 -600-bp upstream region of these genes that comprised their respective promoters as predicted by the online software NNPP (http://www.fruitfly.org/seq_tools/ promoter.html) 6 (29). These upstream regions were separately amplified by employing a Cy5 PCR Labeling Kit (Jena Biosci-ences). The labeled product was then evaluated for its ability to bind VirS by employing EMSA as described under "Experimental procedures." It was observed that the upstream regions of all of the selected genes displayed a shift in mobility, suggesting their direct interaction with VirS ( Fig. 5D). As a negative control to negate a general binding of VirS to DNA, we used a Cy5-labeled 32-bp intergenic region between MbtA and MbtB that contains the IdeR (iron-dependent transcription regulator)-binding region, which did not display any shift in its mobility in the presence of VirS, indicating the specificity of VirS to its own downstream targets (Fig. 5E).

Homology modeling of VirS
We generated a three-dimensional comparative homology model of VirS for performing virtual screening studies and also for carrying out identification of the crucial residues involved in its DNA-binding activity. Because the crystal structure of M. tuberculosis VirS is not available, we generated the three-dimensional structure of VirS by employing comparative homology modeling, as described under "Experimental procedures." It was found by bioinformatic analysis that the C terminus of VirS (amino acids 250 -340) contained the DNA-binding region and that it was homologous to the C terminus of the transcriptional regulators of the AraC/XylS family. Hence, this C-terminal region was employed for multiple-sequence alignment by ClustalW (16) and for model generation by Modeler 9.0 (17) (Fig. 6A). The generated 3D homology model of VirS contained two helix-turn-helix motifs (Fig. 6B). The coordinates of DNA bound to the MarA structure were used to predict the DNA contact sites of VirS (Fig. 6C). The predicted structure of VirS was also validated by using Rampage (18), Errat (19), Verify 3D (20, 21), Pro-Q (22), and Prosa (23,24) software (Table S1). The results obtained from these programs implied that the generated structure of VirS was in a good configuration.

Identification of the critical residues of VirS involved in its DNA-binding activity
On the basis of multiple-sequence alignment and the predicted structure of VirS, 15 residues were selected that might affect the DNA-binding activity of VirS, and structure-guided 6 Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party hosted site.

Role of VirS in acidic stress and its inhibitor identification
mutational studies were performed ( Fig. 7A and Table 3). These residues were substituted with alanine for mutational studies because alanine is a nonpolar and small amino acid that takes care of charge as well as the size of the side chain. For carrying out the site-directed mutagenesis study, we subcloned the virS gene from the pET43.1a-virS construct to another vector, pBEn-SET3a, as described under "Experimental procedures," because the size of the pBEn-3a-virS construct (6.8 kb) is smaller than the pET43.1a-virS (8.3 kb) and thus more appropriate for site-directed mutagenesis studies. The 3a-VirS fusion protein was purified and assessed for its DNAbinding activity with the VirS-binding region of DNA identified by a DNA footprint assay containing 46 bp flanked by 10 bp at each side (66 bp). For ease of work, the ends of this 66-bp region were end-labeled with fluorophore Cy5. The Cy5-labeled DNA was then employed for the determination of DNA-binding activity of VirS by EMSA studies. The purified protein was found to be functionally active and displayed the gel shift that also verified the identified binding region of DNA (data not shown). 0.5 g of protein was selected for further work.
Fifteen VirS mutants were generated by using site-directed mutagenesis as described under "Experimental procedures." The sequences of the primers used for mutagenesis are given in Table S2. All of the mutations were confirmed by DNA C, a consensus DNA sequence for binding of VirS was identified. The protected region identified by footprint analysis was used to define the VirS-binding sites in its downstream targets. The consensus sequence logo was generated by the software WebLogo, in which the height of the stack indicates the degree of conservation at that position in the consensus and the height of an individual letter within a stack represents the relative frequency of that nucleotide at that position. D, the upstream regions of a few of the downstream targets of VirS (ahpD, ppe23, cmtR, ctpG, ppe38, and fadD9) were selected and evaluated for their binding to 1 g of purified VirS by employing EMSA. All of the genes showed a positive binding with VirS, suggesting their direct interaction with VirS. E, as a negative control, an intergenic region between MbtA and MbtB containing the IdeR-binding region was employed, which did not display any shift in the mobility of the DNA in the presence of VirS. As a positive control, VirS and its binding sequence were used, which showed shift in the mobility of its DNA.

Role of VirS in acidic stress and its inhibitor identification
sequencing, and the mutant proteins were separately expressed and purified successfully by using nickel-nitrilotriacetic acid (Ni-NTA) affinity chromatography (data not shown). Mutants Y61A, R73A, and R80A showed poor binding to Ni-NTA resin due to probable changes in the protein folding as a result of mutation. Thus, their lysates were prepared and employed for EMSA studies, and lysate of WT 3a-VirS was used as a control for comparison. The mutants were evaluated for their DNAbinding activity by using EMSA (Fig. 7B).
Mutants R16A, Y61A, and R69A displayed more than 90% reduction in the DNA-binding activity of VirS, suggesting that the corresponding residues are critical for the DNA binding (Fig. 7B). Mutant R20A displayed 84% reduction in the activity, indicating the importance of Arg-20 in DNA binding. Mutants R28A and R73A showed moderate effect on the activity, as they displayed reductions of 34 and 30%, respectively (Fig. 7B). Fig. 7C depicts the bar graph representing the percentage of DNA-binding activity exhibited by the mutants when compared with the WT VirS protein. Mutants H30A, E34A, R38A, and S70A did not display any influence on the DNA-binding activity, suggesting that these residues are not important for the DNA-binding activity of VirS. In conclusion, four residues (Arg-16, Arg-20, Tyr-61, and Arg-69), Figure 6. Prediction of a 3D model of VirS. A, multiple-sequence alignment of amino acid sequences of VirS and its homologous proteins was generated by employing ClustalW. This alignment file was then used to generate the 3D structure of VirS by using the software Modeler 9.0. B, the 3D structure of VirS was generated by comparative homology modeling by using the software Modeler 9.0. The predicted structure of VirS contains two HTH motifs. C, structure of VirS with the predicted DNA-binding sites. B and C were generated by PyMOL (58).

Role of VirS in acidic stress and its inhibitor identification
when substituted with alanine, resulted in the loss of VirS DNA binding activity by more than 84% compared with the WT VirS. These residues of VirS might directly make contacts with the cognate DNA and can be used as specific target sites for rational design of inhibitors against VirS for anti-tubercular drug development.

Identification of the inhibitory molecules against VirS by structure-based virtual screening and their evaluation against DNA-binding activity of VirS and M. tuberculosis growth
Among the homologous proteins of VirS, only three structures (i.e. MarA (PDB entry 1BL0), Rob (PDB entry 1D5Y), and AdpA (PDB entry 3W6V)) were available in ligand-bound form (25)(26)(27). The DNA contact sites of VirS were predicted by superimposition with DNA-bound MarA structure (Fig. 8A). On superimposing the 3D structure of VirS with the structure of MarA, two DNA contact sites in VirS were observed (Fig.  8B). Hence, virtual screening was carried out at both of the DNA contact sites of VirS by using Autodock4.2 software (http://autodock.scripps.edu/) 6 (28). The grid covering DNA contact site 1 of VirS is represented in Fig. 8, C and D. The residues (i.e. Glu-5, Gln-19, Arg-28, and His-30), were surrounding the grid and found to be covered in the generated grid (Fig. 8E). Similarly, another grid was generated against DNA contact site 2 of VirS (Fig. 8F), and its enlarged representation is shown in Fig. 8G. The residues found to surround and be covered in the grid were Ala-56, Ser-53, Tyr-61, Ser-62, Glu-63, and Gln-64 (Fig. 8H). A filtered library of the NCI (National Institutes of Health) Open Database was employed for virtual screening against the DNA contact sites of VirS (30). A total of 114 compounds, 57 compounds each from docking at site 1 and site 2, were procured from NCI-DTP, numbered as V1A-V57A (site 1) and V1B-V57B (site 2), and evaluated for their ability to inhibit the DNA-binding activity of VirS at a concentration of 100 g/ml (Fig. 9). 26 compounds (i.e. 10 from site 1 and 16 from site 2) displayed more than 60% inhibition with 15 com-pounds (4 from site 1 and 11 from site 2) exhibiting greater than 80% inhibition at 100 g/ml (Fig. 9, A and B and Table S3). These 26 compounds were further subjected to dose-response studies for the determination of their IC 50 values (Table S3). Compound V7B with an IC 50 value of 1.8 g/ml was found to be the most potent one, whereas compound V29B displayed an IC 50 value of 16.38 g/ml. The rest of the compounds showed IC 50 values in the range of 17-75 g/ml (Table S3).
These 26 compounds were also evaluated for their ability to inhibit the growth of M. tuberculosis in broth culture by subjecting them to a resazurin microtiter assay (31)(32)(33). Rifampicin was used as a positive control in the assay, which exhibited an MIC 90 of 0.25 g/ml, and DMSO was used as a vehicle control. Fourteen compounds displayed growth inhibition of M. tuberculosis with varying MIC 90 values (Table 4 and Fig. S3). The most potent inhibition against M. tuberculosis was exhibited by compound V15B with an MIC 90 value of 5 g/ml. Compound V7B, which was the most potent compound in the in vitro assay, displayed an MIC 90 of 40 g/ml. To summarize, we identified 14 compounds (listed in Table 4) that exhibited growth inhibition against M. tuberculosis as well as inhibited DNA-binding activity of VirS in biochemical assays.
We then further evaluated the cytotoxicity of eight compounds (V4A, V6A, V8A, V44A, V55A, V7B, V15B, and V56B) that exhibited an MIC 90 Յ 40 g/ml and were also active against the DNAbinding activity of VirS. Various mammalian cell lines were employed, and IC 50 values of the compounds were determined as described under "Experimental procedures." The most potent inhibitor of VirS i.e. V7B (IC 50 ϭ 1.8 g/ml, MIC 90 ϭ 40 g/ml) was found to be noncytotoxic until a concentration of 100 g/ml in all of the cell lines (Table S4). Another potent inhibitor of M. tuberculosis growth (i.e. V15B (MIC 90 ϭ 5 g/ml and IC 50 ϭ 31.21 g/ml)) was also found to be noncytotoxic until a concentration of 150 g/ml in all of the cell lines (Table S4). Compounds V4A and V8A also displayed negligible cytotoxicity in all of the employed cell lines until 200 and 150 g/ml, respectively. The rest of the compounds displayed moderate cytotoxicity in at least one or two of the three employed cell lines (Table S4).
Four compounds (V4A, V8A, V7B, and V15B) that displayed inhibition of VirS DNA-binding activity as well as growth of the pathogen and did not display any significant cytotoxicity were further evaluated for their potency to inhibit the growth of the pathogen inside macrophages as described under "Experimental procedures." It was found that compounds V7B and V15B inhibited the intracellular growth of the pathogen at Ն80 and Ն40 g/ml, respectively (Fig. 10). To rule out the possibility of this inhibition due to cytotoxic effects of the compounds, uninfected THP-1 cells were also subjected to the same concentrations of the compounds. V7B and V15B were found to be noncytotoxic to THP-1 cells until 320 and 100 g/ml, respectively. Compounds V4A and V8A failed to inhibit the intracellular growth of the pathogen even at a con-

Role of VirS in acidic stress and its inhibitor identification
centration of 100 g/ml (data not shown). Thus, these results suggested that the compounds V7B and V15B also possess the ability to inhibit the intracellular growth of the pathogen; however, further optimization of these compounds is required to improve their potency.

Structure-activity relationship studies for compounds V7B and V15B
On the basis of IC 50 , MIC 90 , cytotoxicity, and macrophage studies, two hit molecules (i.e. compounds V7B and V15B) were shortlisted for further SAR studies. For this, they were subjected to a similarity search for identifying their analogs by using the NCI database employing the software Open Babel. The molecules were selected on the basis of a Tanimoto coefficient of Ն0.7. Seventeen molecules similar to compound V7B and 10 molecules similar to compound V15B were procured from NCI-DTP and evaluated for their ability to inhibit VirS DNA-binding activity to identify the structural features respon-  Figure 10. Evaluation of the inhibitory potential of compounds against intracellular survival of the pathogen. Compounds V7B and V15B were evaluated for their inhibitory potential against intracellular survival of M. tuberculosis inside macrophages. A, compound V7B displayed inhibition of growth of the pathogen at 80 g/ml and above. Cytotoxicity of compound V7B was determined by resazurin dye. The bottom panel depicts the cytotoxicity exhibited by compound V7B to the macrophages. Except at the highest concentration employed, compound V7B did not display any cytotoxicity to the macrophages. B, compound V15B displayed inhibition of growth of the pathogen at 40 g/ml and above. The bottom panel depicts the cytotoxicity exhibited by compound V15B to the macrophages. Compound V15B did not display any cytotoxicity to the macrophages.

Role of VirS in acidic stress and its inhibitor identification
sible for inhibiting VirS activity. These molecules were named V7B-1 to V7B-17 and V15B-1 to V15B-10. The procured molecules were evaluated for their ability to inhibit DNA-binding activity of VirS at a concentration of 100 g/ml, and it was found that four analogs of V7B (V7B-4, V7B-5, V7B-8, and V7B-9) and one analog of V15B (V15B-2) displayed Ͼ60% inhibition of VirS DNA-binding activity at 100 g/ml, whereas the rest of them exhibited lesser inhibition (Tables S5 and S6). These five analogs were further evaluated for the determination of their IC 50 values. They were found to display a reduced inhibitory potential against VirS DNA-binding activity compared with their parent molecules, which is tabulated in Table S5 and S6. It was found that none of these molecules displayed an enhanced potential to inhibit VirS; however, their structures were carefully analyzed to derive the structure-activity relationship.
A detailed comparison of the structures of V7B and its analogs revealed the presence of a benzene sulfonamide moiety in all of the structures except in V7B-5 and V7B-17. The structure comparison of a few of the analogs of V7B with its parent compound is shown in Fig. 12. The absence of acetamide moiety in compound V7B-9 (IC 50 ϭ 2.9 g/ml) led to reduction in the potency of V7B-9 compared with its parent compound V7B (IC 50 ϭ 1.8 g/ml). Compound V7B-4 contains 2-amino-3-chloro-1,4-napthoquinone moiety attached to benzene sulfonamide, which resulted in a poor IC 50 . Substitution of methane diamine by benzene sulfonamide attached to 2-amino-3chloro-1,4-napthoquinone moiety in compound V7B-8 also resulted in poor IC 50 (12.4 g/ml) compared with its parent molecule. Replacement of benzene sulfonamide with sulfonyl benzene and modification of 1,4-dihydronapthalene-1-one to 1,4-dihydronapthalene-1,4-dione in parent compound resulted in reduced potency to inhibit VirS DNA-binding activity, as observed in compound V7B-5. Thus, the structural comparisons of all of these molecules with their activity emphasize the importance of the acetamide ring attached to benzene sulfonamide, which might be essential for the potent inhibition of DNA-binding activity of VirS. Second, it was also found that bulkier modifications to the acetamide group attached to benzene sulfonamide may lead to poor inhibition of VirS DNAbinding activity, as observed in the cases of V7B-10, V7B-13, V7B-15, and V7B-17.
Among the procured analogs of compound V15B, only one compound was found to be functionally active against the VirS DNA-binding activity whose structure comparison with its parent compound is shown in Fig. 13. Replacement of the isoquinoline ring from the parent molecule (V15B) with a benzene ring, pyridine ring, and methyl pyridine displayed poor inhibi-

Role of VirS in acidic stress and its inhibitor identification
tion of VirS DNA-binding activity compared with its parent molecule, as observed in compounds V15B-1, V15B-3, V15B-5, and V15B-6. Further modification of the isoquinoline ring to 1-methyl-3,4-dihydroisoquinoline-4-one and its attachment to the 2-(cyclohexa-1,3-dien-1-yl) pyridine moiety in compound V15B-2 showed a potential to inhibit the DNA-binding activity of VirS; however, it displayed a slightly lower IC 50 than the parent molecule. It was also observed that modifications to the 2-(4,5-dihydro-1H-pyrrol-2-y) pyridine moiety failed to show any inhibition against DNA-binding activity of VirS, as observed in the cases of V15B-8 and V15B-10. Thus, the structural comparisons of all of these molecules with their activities emphasize the importance of the isoquinoline ring attached to 3-nitroso-2-phenyl-pyrrole as present in the parent molecule,

Role of VirS in acidic stress and its inhibitor identification
which is essential for inhibition of the DNA-binding activity of VirS.

Discussion
In this study, we attempted to dissect the role of VirS in the survival of M. tuberculosis under acidic stress and in the maintenance of pH homeostasis of the pathogen. We found that the loss of virS reduced the ability of M. tuberculosis to block phagosomallysosomal fusion in the activated macrophages as well as its ability to survive in acidic conditions, and its absence also led to a disturbance in the maintenance of the intrabacterial pH under acidic stress and inside activated macrophages, thereby suggesting the role of VirS in evading stressful acidic conditions.
To understand how VirS mediates the regulation of gene expression in acid stress, microarray studies were carried out. Transcriptome profiling has provided insights into the role of VirS in the regulation of various pathways, such as those involved in cellular metabolism; cell-wall and lipid metabolism; PE/PPE family proteins; energy pathways and respiration; efflux pumps and ion transporters; detoxification; information pathways; and transcription regulators in acidic stress.
First, we found a striking influence of VirS on the regulation of a number of genes (induced or repressed) involved in fatty acid oxidation and biosynthesis that can contribute to remodeling of the M. tuberculosis cell wall. It has been a general belief that the cell wall of mycobacteria remodels itself when exposed to hostile environments to become impervious to outside harsher conditions (13,34). We observed VirS-dependent up-regulation of genes involved in fatty acid ␤-oxidation (fadB2 and fadE35), lipase and hydrolyzing enzymes (Rv2672 and Rv0867c), and membrane-associated proteins (Rv2911, Rv0314c, Rv0531, and Rv2180c) and down-regulation of genes encoding enzymes involved in polyketide synthase (pks2, pks3, pks6, and papA1), membrane proteins (mmpS3, Rv2869c, Rv0412c, Rv2091c, and Rv0677c), phenolic glycolipids (fadD22), biosynthesis and export of siderophores (mmpS5 and mmpL5), and lipoprotein (Rv2290). VirS-regulated expression of genes involved in lipid metabolism suggests its pivotal role in remodeling of the cell envelope that can facilitate intracellular survival under acidic stress.
Our expression profile also revealed the regulation of metabolic enzymes by VirS at acidic pH. It has been reported that there is an induction of genes involved in metabolism via an anaplerotic node at acidic pH in M. tuberculosis. We observed VirS-dependent induction of a glyoxylate shunt comprising isocitrate lyase (icl) at acidic pH, which replenishes intermediates of the TCA cycle and gluconeogenesis and decreases carbon flux through the TCA cycle (35). The gene most up-regulated by VirS in acid stress is Rv1405c, which has also been reported previously to be induced under acid stress and is implicated in adaptation of M. tuberculosis to the intramacrophage environment; however, the substrate for Rv1405c is still not known (36 -38). It has been reported that fermentative pathways become up-regulated under acidic stress in B. subtilis and Bacillus cereus (39,40). The genes for alcohol dehydrogenase were found to be induced that catalyze the conversion of pyruvate to ethanol, generating CO 2 and dissipating H ϩ to deal with low intracellular pH or restoration of NAD ϩ /NADH bal-ance (41). Consistent with this fact, we found VirS-dependent induction of pyruvate decarboxylase (pdc) and alcohol dehydrogenases (adhE1 and Rv1530) at acidic pH, which may reduce acidity through NAD(P)H that transfers electrons to electron transport chain and pumps protons out of the cell (40). We also found a strong influence of VirS on the expression of genes related to cellular metabolism. Genes encoding enzymes for the synthesis of amino acids feeding into the TCA cycle (argJ), nitrate reductase (narX), for the synthesis of heme and porphyrin (hemB) were found to be down-regulated by VirS at acidic pH. Genes for sulfur reduction are essential for the synthesis of cysteine, methionine, mycothiol, and sulfolipids, and these reduced sulfur products contribute to M. tuberculosis pathogenesis. We found that genes for sulfur reduction (sirA and cysH) were also down-regulated in a VirS-dependent manner. These data suggest that changes in the central metabolism by VirS enable M. tuberculosis to adapt to the acidic intraphagosomal environment.
In the context of maintenance of intrabacterial pH, we observed VirS-dependent induction of efflux pumps and ion transporters (kdpC and ctpG). It is reported that a decrease in pH leads to reduction in membrane potential and K ϩ uptake, which induces the genes encoding K ϩ transporters (42). Induction of potassium transporters is necessary for the maintenance of internal pH homeostasis and promotes survival in the acidic environment (42). Kdp is a K ϩ uptake system known to be induced in an acidic environment, which may be required for bacterial adaptation at low pH levels (42). As Kdp is a K ϩ /H ϩ antiporter, a maintenance of cytoplasmic pH is expected by pumping protons out of the cell. The metal transporter CtpG that transports Cd 2ϩ across the mycobacterial cell membrane (43), which is also regulated by VirS. Such transporters have been previously reported to be critical for ion homeostasis and bacterial survival (43,44). Thus, regulation of these cell-wall transporters suggests the role of VirS in response to phagosomal acidification and in the maintenance of intrabacterial pH.
Transcriptome profiling also revealed that VirS further activates other transcriptional regulators (Rv1994c and Rv3765c); thus, it might be acting as a master regulator under these conditions inducing the acidic responses of the cell. This is in agreement with the earlier study by Talaat et al. (14), which implicated that VirS could be the master regulator of reactivation of tuberculosis. CmtR (Rv1994c) is an AsrR-SmtB repressor that senses cadmium and modulates the expression of genes related to the efflux and detoxification of this metal ion (45). The regulation of CmtR and CtpG by VirS indicates its importance in the efflux of the cadmium ions at acidic pH to regulate ion homeostasis.
A down-regulation in the VirS-dependent gene expression was observed in ribosome biogenesis pathway (rpsk) and protein translation (infC). Because the genes involved in ribosome biogenesis and protein translation are the major consumer of cellular energy, their down-regulation is also expected (46,47). Repression of these genes indicates a slowing down of the bacterial growth at acidic stress, which is mediated by VirS. It is well-understood that there is a strong connection between acid stress and oxidative stress (40,48), and our transcriptome pro-

Role of VirS in acidic stress and its inhibitor identification
filing revealed VirS-dependent up-regulation of genes involved in protection from oxidative stress (i.e. ahpC, ahpD, and Rv3177), which may aid bacterial survival at acidic pH. We also observed a decreased expression of aerobic respiration (nuok), and such repression at acidic pH has also been documented earlier (49). We also found VirS-dependent down-regulation of cytochrome assembly ATP-binding protein (cydD) at acidic pH, which promotes proton export for the maintenance of intrabacterial pH (40).
Based on the analysis of our transcriptional profile, we propose a model for the survival of M. tuberculosis under acidic conditions mediated by VirS, although further investigations are needed to confirm these proposed mechanisms and the involvement of the downstream genes of VirS in pH homeostasis (Fig. 14).
The importance of VirS for the survival of bacterial cells under acidic stress and its role in regulating multiple downstream targets signify that it is an attractive drug target for the development of new anti-tuberculosis drugs. Toward this, we characterized it further by evaluating the DNA-binding activity of VirS and by determining the DNA sequence to which VirS binds by employing DNA footprinting. Further, we generated a consensus logo by WebLogo to map the conservation of nucleotides in the sequences employed (Fig. 5C) (15). Additionally, we performed EMSA with the upstream regions of a few downstream targets of VirS, which were found to exhibit direct interaction with VirS, as evident by a shift in the EMSA studies (Fig.  5D).
Site-directed mutational studies revealed four residues (i.e. Arg-16, Arg-20, Tyr-61, and Arg-69) as the critical residues for binding of VirS to its cognate DNA. Residue Tyr-61 was found to be a conserved residue in other homologous proteins. Arg-20 of VirS corresponds to Arg-46 in MarA, which penetrates into the major groove of DNA and makes hydrogen bonds with bases. Similarly, Arg-69 of VirS is homologous to Arg-96 of MarA, which hydrogen-bonds with the base of DNA. Similarly, Arg-16 of VirS corresponds to Trp-42 in MarA; side chains of Trp-42 are involved in the van der Waals contacts with the Exposure to acidic pH causes up-regulation of VirS, which in turn regulates a number of genes probably involved in combating acidic stress faced by the pathogen. First, VirS regulates genes related to cellular metabolism that may help bacteria to adapt to acidic stress. Second, it regulates the genes related to various cell-wall components that might help in remodeling of the cell envelope to make it impervious to acidic stress. It also up-regulates the expression of efflux pumps and ion transporters that may help in the maintenance of intrabacterial pH by exporting an excess of protons outside the cell. VirS also regulates genes related to detoxification, the information pathway, and other transcriptional regulators that may help bacteria to resist and adapt to acidic stress.

Role of VirS in acidic stress and its inhibitor identification
DNA. The essentiality of the identified critical residues for DNA-binding activity suggests their importance in the interaction of VirS with the DNA. Hence, these residues could be utilized for rational design of inhibitors to target the DNA-binding activity of VirS.
Because VirS appears as an attractive target for the discovery of new inhibitors against M. tuberculosis, we performed structure-based virtual screening for the identification of small molecules against the DNA binding activity of VirS, which may serve as hit molecules for novel anti-tuberculosis drug development. We identified two hit molecules (i.e. V7B and V15B) that exhibited potent inhibition of VirS activity and inhibited the growth of M. tuberculosis in vitro as well as inside macrophages. However, IC 50 and MIC values of both of these compounds did not correlate well, and the most probable reason for such disparity could be attributed to the chemical nature of the potent enzyme inhibitors leading to their poor penetration into M. tuberculosis cells, resulting in weak MIC values. Also, the target specificity of these compounds inside M. tuberculosis still needs to be validated. Hence, these two compounds can be further explored to increase their potency and could be a subject of lead optimization.
Thus, we have explored the promising drug target VirS by identifying a few crucial residues for its activity that could help in designing rational inhibitors as well as by identifying few moieties that inhibit VirS activity as well as M. tuberculosis growth. This is the first report of exploring VirS as a target for drug development against M. tuberculosis.

Experimental procedures
Oligonucleotides for PCR and for site-directed mutagenesis, FITC, IFN-␥, poly(dI-dC), tyloxapol, and complete mini EDTA-free protease inhibitor mixture tablets were purchased from Sigma-Aldrich. DpnI and Pfu polymerase were from New England Biolabs (Ipswich, MA). Cy5-labeled dsDNA was customized from IDT (Coralville, IA). Ni-NTA was obtained from IBA Life Sciences (Goettingen, Germany). RPMI 1640, Dulbecco's modified Eagle's medium, LysoTracker Red DND-99, Prolong Gold antifade reagent with DAPI, antibiotic-antimycotic, and fetal bovine serum (FBS) were obtained from Thermo Fisher Scientific (Waltham, MA). Difco Middlebrook 7H9 (MB7H9) medium, 7H11 agar, and ADC were obtained from BD Biosciences (Franklin Lakes, NJ). The compounds were obtained from the Drug Synthesis and Chemistry Branch, Developmental Therapeutics Program, Division of Cancer Treatment and Diagnosis, NCI, National Institutes of Health (Bethesda, MD).

Bacterial strains and growth conditions
M. tuberculosis strains were grown in MB7H9 medium supplemented with 1ϫ ADC, 0.2% glycerol and 0.2% Tween 80 with constant shaking at 200 rpm or on 7H11 agar plates containing 10% oleic acid/albumin/dextrose/catalase (OADC) and 0.2% glycerol at 37°C. The pH of MB7H9 medium was adjusted by using HCl, and instead of Tween 80, 0.02% tyloxapol was used, wherever required. Hygromycin B and kanamycin were used at a concentration of 50 and 25 g/ml, respectively, for mycobacteria. Phosphate-citrate buffers were prepared from 200 mM sodium phosphate and 100 mM citric acid. E. coli strains XL-1 Blue and BL21(DE3) were grown in Luria-Bertani (LB) broth or on LB agar. Ampicillin was used at a concentration of 50 g/ml for E. coli.

Growth kinetics and colocalization studies by confocal microscopy
The virS mutant of M. tuberculosis and virS complemented strain were earlier generated in our laboratory (12,13). The M. tuberculosis Erdman, virS mutant, and virS complemented strain were grown in MB7H9 medium adjusted to a pH of 6.6, 5.5, 5.0, and 4.5 at 37°C with constant shaking at 200 rpm. The growth was monitored daily by measuring the A 600 nm . In case of cfu studies, the cells were appropriately diluted and plated at various time points on MB7H11 agar plates. The plates were kept at 37°C for 3-4 weeks for the appearance of colonies. The phagosomal maturation arrest of M. tuberculosis was studied as described previously (50). For colocalization studies in activated macrophages, PMA-activated macrophages were washed once with plain RPMI 1640 medium, and cells were further incubated for 24 h in RPMI 1640 medium containing 10% FBS and 10 ng/ml IFN-␥ before infection. After infection and Lyso-Tracker Red treatment, coverslips were mounted by using Prolong Gold antifade reagent and visualized by a Leica TCS Sp5 confocal laser-scanning microscope (Leica Microsystems, Mannheim, Germany). The percentage colocalization of FITClabeled mycobacteria with LysoTracker Red was calculated by analyzing ϳ100 phagosomes.

Intrabacterial pH measurement assays
The plasmid encoding ratiometric pH-GFP (a kind gift from Dr. Sabine Ehrt, Department of Microbiology and Immunology, Weill Cornell Medical College, New York; received from Dr. Amit Singh, Department of Microbiology and Cell Biology, Indian Institute of Science, Bengaluru, India) was electroporated into M. tuberculosis Erdman, virS mutant, and complemented strain, and the recombinants were selected on the basis of fluorescence of GFP on 7H11-OADC-hygromycin B plates. The selected recombinant clones were grown in MB7H9 medium in the presence of hygromycin until A 600 nm of 0.3-0.5. The cultures were harvested, and pellets were resuspended in phosphate-citrate buffer, pH 7.4 and 4.5, supplemented with ADC, 0.2% Tween 80, and hygromycin B. The intrabacterial pH measurements were carried out after incubation of the bacteria in the respective buffers for different time intervals. The bacterial cultures at different time points were harvested and concentrated to 20 times by centrifugation to increase GFP signal. Fluorescence was analyzed in a Cary Varian spectrofluorometer with excitation wavelengths of 395 and 475 nm at an emission wavelength of 510 nm. The intrabacterial pH was derived by interpolating the ratio of 395/475-nm fluorescence intensity on a standard curve. The standard curve for intrabacterial pH measurement was generated by incubating the lysate containing GFP with the phosphate-citrate buffers of pH ranging from 5.5 to 8.5 in a 96-well plate (1). The plate was read at excitation wavelengths of 395 and 475 nm and an emission wavelength of 510 nm.

Role of VirS in acidic stress and its inhibitor identification
For pH IB measurements of bacteria residing inside macrophages, 5 ϫ 10 5 PMA-differentiated THP-1 cells were seeded in poly-L-lysinecoated glass bottom culture dishes of 35-mm thickness. PMA-activated macrophages were washed once with plain RPMI 1640 medium, and cells were further incubated for 24 h in RPMI 1640 medium containing 10% FBS and 10 ng/ml IFN-␥ before infection. Macrophages were infected with M. tuberculosis at an MOI of 1:5 (macrophage/bacteria) for 4 h followed by amikacin treatment of 1 h. After different time points postinfection, macrophages were analyzed under confocal microscopy to determine the intensities of GFP (exciters D405 and D476, emitter 535). All of the images were acquired within an experiment under identical conditions. A bacterial group was defined as at least one bacterium but may consist of 2-5 bacilli. The bacterial groups were analyzed for their emission intensities when excited at 405 and 476 nm, and the ratios of intensities at 405:476 nm were calculated. pH IB were derived by interpolating the 405:476 nm intensity ratios on a standard curve.

Microarray studies
M. tuberculosis Erdman, virS mutant, and virS complemented strain were grown until an A 600 nm of 0.3-0.5, and cells were harvested by centrifugation at 4500 rpm for 10 min. Harvested cells were exposed to MB7H9 medium adjusted to pH 6.6 and 4.5 and grown at 37°C with constant shaking at 200 rpm for 1 h. Total RNA was isolated from all of the mycobacterial strains (grown at pH 6.6 and 4.5) and was processed and hybridized to a M. tuberculosis whole-genome Gene Expression Profiling microarray G2509F (AMADID: G2509F_43893, Agilent Technologies). DNA microarrays and data analysis were performed at the University of Delhi South Campus Microarray Centre (UDSCMAC). Slides were scanned on a microarray scanner (Agilent Technologies) and analyzed using the Gene-Spring software.

Cloning and expression of virS
The virS gene was PCR-amplified by employing M. tuberculosis H37Rv genomic DNA as the template and by using the primers F-virS containing SpeI, thrombin site, and EcoRI followed by a His 6 tag and R-virS-containing HindIII site. The sequence of the primers is given in Table S2. The amplified PCR product after digestion with SpeI and HindIII was cloned into pET43.1a vector digested with the same enzymes to create pET43.1a-virS (NusA-VirS) for the synthesis of N-terminal His-tagged VirS. To subclone virS from pET43.1a-virS construct into vector pBEn-SET3a, the positive clone of pET43.1a-virS was digested with EcoRI and HindIII to obtain a fallout of ϳ1.0 kb, corresponding to the size of virS. The pBEn-SET3a vector was also digested with the same enzymes. The digested vector pBEn-SET3a and digested virS were ligated to obtain the recombinant clone pBEn-virS (3a-VirS). Positive clones were confirmed by restriction digestion. For expression, E. coli BL21(DE3) cells transformed with pET43.1a-virS and pBEn-virS, were grown at 37°C in Luria-Bertani medium containing 50 g/ml ampicillin until A 600 nm of 0.8. The cultures (i.e. NusA-VirS and 3a-VirS) were then induced with 1 mM isopropyl-1-thio-␤-D-galactopyranoside and were allowed to grow for 16 h at 16 and 18°C, respectively. The cells were harvested by centrifugation at 4°C and 6000 ϫ g for 10 min.

Purification of 3a-VirS
E. coli BL21(DE3) cells were transformed with pBEn-virS, and the cells were grown and induced as described above. For purification, the cells from the induced culture were harvested and resuspended in lysis buffer containing 1ϫ PBS, pH 7.4, 10 mM imidazole, 5 mM ␤-mercaptoethanol, 1 mM phenylmethylsulfonyl fluoride, and complete mini EDTA-free protease inhibitor mixture tablet and lysed by sonication followed by centrifugation to remove cell debris (15, 000 ϫ g, 45 min, 4°C). Clarified lysate was loaded onto a Ni-NTA resin pre-equilibrated with lysis buffer. Washing of the column was performed initially with buffer containing 1ϫ PBS, 10 mM imidazole followed by one wash each with the same buffer containing 20 mM and 50 mM imidazole. The protein was eluted with buffer containing 250 mM imidazole, and the purity of the protein was analyzed by SDS-PAGE by using a 12.5% polyacrylamide gel. The fractions containing VirS were pooled, and 10% glycerol was added. The pooled protein was dialyzed overnight against 20 mM Tris, pH 8.0, and stored at Ϫ80°C.

Gel-shift assay for NusA-VirS fusion protein
DNA probes for gel mobility-shift assays were amplified by PCR from vector pSG10 comprising the intergenic region between virS and mymA, which is hypothesized to be bound by VirS (12). Two primers (UP6 and DW2) were utilized for the amplification of a 325-bp intergenic region. (Table S2). Primer DW2 was end-labeled with [␥-32 P]ATP by using T4 polynucleotide kinase. The PCR product was run on 1.5% agarose gel, and the amplified band was cut and gel-eluted by using GFX columns and quantitated by a ␤-counter (TRI-CARB 2900TR Liquid Scintillation Counter, PerkinElmer Life Sciences). Lysates of E. coli cells overexpressing either the NusA-VirS fusion protein or only NusA protein were prepared. The reaction mixture consisted of 120,000 cpm of labeled DNA, binding mix, and different amounts (i.e. 13, 26, and 52 g) of the lysate. The binding mix consisted of 1ϫ binding buffer, 1 mg/ml BSA, 1 mM DTT, 4% PEG-8000, and 1 g of poly(dI-dC). The samples were incubated for 30 min at 37°C. Samples were then electrophoresed on a 3.5% polyacrylamide gel, dried, and autoradiographed.

DNA footprinting assay
Binding reaction mixture was prepared as described above in a total volume of 30 l that consisted of 126,000 cpm of labeled DNA probe and 52 g of lysate and incubated at 37°C for 30 min. After incubation, the binding reaction mix was subjected to DNase I digestion (0.2 units) for 1 min, and digestion was terminated by the addition of 180 l of stop solution (300 mM sodium acetate, pH 7.0, 10 mM EDTA, pH 8.0, 0.3 mg/ml yeast t-RNA, 0.01% SDS). Samples were extracted with phenol-chloroform, and precipitation of nucleic acid was performed with ethanol. DNA was resuspended in water, and formamide dye was added. A GϩA ladder was also prepared by using the same probe, and formamide dye was mixed. The samples and GϩA ladder were preheated at 90°C for 10 min before loading. These samples were separated on 6% polyacrylamide sequencing gel

Role of VirS in acidic stress and its inhibitor identification
and electrophoresed (in 0.5ϫ TBE buffer and 0.5 M sodium acetate in the lower chamber and 0.5ϫ TBE in the upper chamber) at 1050 V. Gel was dried and exposed for 14 h, followed by autoradiography.

Homology modeling of VirS
The amino acid sequence similarity search was performed by using the NCBI-BLAST server in the PDB structure database. Based on the BLAST results, the first few homologs of VirS were identified, namely ToxT of Vibrio cholera (PDB entry 3GBG) (51), gadX of E. coli (PDB entry 3MKL), 7 AdpA of Streptomyces griseus (PDB entry 3W6V) (27), unknown protein of Listeria innocula (PDB entry 3OOU), 8 AraC of E. coli (PDB entry 2K9S) (52), MarA of E. coli (PDB entry 1BL0) (25), yesn of Fusobacterium nucleatum (PDB entry 3LSG), 9 ternary complex of MarA, DNA, and C terminus of RNA polymerase of E. coli (PDB entry 1XS9) (53), Rob of E. coli (PDB entry 1D5Y) (26), transcription regulator of Chromobacterium violaceum (PDB entry 3OIO), 10 and XylR of E. coli (PDB entry 4FE4) (54). The sequences of all of these homologs were aligned with the C-terminal region (amino acids 250 -340) of the M. tuberculosis VirS sequence by using the software ClustalW. The result of sequence alignments was then submitted to the software Modeler 9.0 to generate a 3D structure of VirS, and the generated homology model was validated by employing Rampage, Errat, Verify 3D, Pro-Q, and Prosa software.

Site-directed mutagenesis
Amino acid substitutions were generated by site-directed mutagenesis kit provided by Stratagene (La Jolla, CA) according to the manufacturer's instructions. The E. coli XL1 Blue cells were transformed to get the colonies on ampicillin plates. These colonies were grown in LB medium, and DNA was isolated. The mutations in the DNA were confirmed by DNA sequencing. The primers used in mutagenesis are listed in Table S2.

EMSA for purified 3a-VirS
0.5 g of purified 3a-VirS protein was taken in a reaction mixture that contained 1ϫ binding buffer (10 mM Tris acetate, pH 8.0, 50 mM KCl, 5% glycerol), 4% PEG, 1 g of poly(dI-dC), 100 g/ml BSA, 1 mM DTT, and 0.8 pmol of 5Ј Cy5 fluorophore-labeled DNA and incubated for 30 min at 37°C. This labeled dsDNA represents the sequence to which VirS binds as identified by a DNA footprinting study with 10-bp flanking regions on both sides (66-bp dsDNA), which were labeled with Cy5 fluorophore at the 5Ј ends. After incubation, samples were loaded onto a 5% polyacrylamide gel and electrophoresed in 1ϫ TAE buffer overnight at 32 V. The next day, fluorescence was detected by a phosphor imager (model FLA-9000, FUJIFILM Corp., Minato-ku, Tokyo, Japan).

Regulation of downstream genes by VirS using EMSA
Upstream regions of the genes were PCR-amplified by using the Cy5 PCR Labeling Kit (Jena Biosciences). The amplified products were gel-eluted and quantified. Forty fmol of labeled product was incubated with 1 g of purified VirS in the binding buffer (composition same as mentioned above) for 30 min at 37°C. Samples were electrophoresed onto a 5% polyacrylamide gel in 1ϫ TAE buffer overnight. The bands were visualized using a phosphor imager.

Virtual screening
The 3D structure generated by Modeler 9.0 was superimposed with the crystal structure of MarA (PDB entry 1BL0), and the coordinates of DNA bound to MarA structure were used to predict the DNA contact sites of VirS. Virtual screening was carried out against both of the DNA contact sites by using a small-molecule library of the NCI Open Database consisting of 260,071 compounds. These compounds were shortlisted by using the FAF server, which resulted in 95,748 compounds that were used for screening against the DNA contact sites of VirS by using Autodock4.2 (30). The top compounds were procured from NCI-DTP, depending upon the availability for conducting inhibition studies. Fifty-seven compounds each on the basis of docking at site 1 and site 2 were procured from NCI-DTP. Compounds obtained from docking at site 1 were denoted as V1A-V57A. Compound V7A was not processed further due to its toxicity. Compounds obtained from docking at site 2 were denoted as V1B-V57B.

Assay for inhibition of DNA-binding activity of VirS
All the compounds were dissolved in DMSO to prepare a stock solution of 2 mg/ml and were screened for their ability to inhibit the DNA-binding activity of VirS at a concentration of 100 g/ml. The reaction mixture was same as described above. Binding buffer, purified protein (3a-VirS), and the individual test compounds were incubated at 37°C for 30 min followed by the addition of the rest of the components and further incubated for 30 min at 37°C. The rest of the steps were the same as mentioned above. Sample containing DMSO was used as a negative control for the assay. IC 50 is defined as the concentration of the compound that exhibits 50% inhibition of the activity.

Assay for inhibition of growth of M. tuberculosis
The inhibition of M. tuberculosis growth was evaluated by employing a resazurin microtiter assay as described previously (30,55). Blue color depicts no growth, whereas conversion to pink color depicts viable cells. Fluorescence was measured at an emission wavelength of 590 nm with an excitation wavelength of 530 nm, and percentage inhibition was calculated. Rifampicin was employed as a positive control, whereas wells with DMSO were used as a negative control. MIC 90 is defined as the concentration of compound at which growth of the pathogen was inhibited by 90%.

Evaluation of cytotoxicity of the compounds
Compounds were evaluated for their cytotoxicity against three mammalian cell lines, namely THP-1 (human monocytic cell line), HEK (human embryonic kidney cell line), and HeLa (cervical cell line), by employing the protocol described previously (30,55). IC 50 is defined as the concentration of the compound at which 50% growth of the cell line is inhibited.

Evaluation of inhibitory potential of the compounds against intracellular survival of M. tuberculosis inside macrophages
The ability of the compounds to inhibit the intracellular growth of the pathogen was evaluated as described previously (56).
Author contributions-S. S., G. K., and A. K. T. designed the experiments. S. S. and N. G. performed the experiments. S. S. wrote the manuscript; G. K. and A. K. T. proofread, corrected, and edited the manuscript. G. K. and A. K. T. provided overall supervision of the work.