# Protein kinase B controls Mycobacterium tuberculosis growth via phosphorylation of the transcriptional regulator Lsr2 at threonine 112

Kawther Alqaseer, Obolbek Turapov, Philippe Barthe, Heena Jagatia, Angélique de Visch, Christian Roumestand, Malgorzata Wegrzyn, Iona Bartek, Martin Voskuil, Helen O'Hare, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Kawther Alqaseer, Obolbek Turapov, Philippe Barthe, Heena Jagatia, Angélique de Visch, et al. Protein kinase B controls Mycobacterium tuberculosis growth via phosphorylation of the transcriptional regulator Lsr2 at threonine 112. Molecular Microbiology, 2019, 112 (6), pp.1847-1862. 10.1111/mmi.14398. hal-02363543

HAL Id: hal-02363543

https://hal.science/hal-02363543

Submitted on 2 Mar 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Protein kinase B controls *Mycobacterium tuberculosis* growth via phosphorylation of the transcriptional regulator Lsr2 at threonine 112

Kawther Alqaseer,<fnref n="1" />,<fnref n="2" />,† Obolbek Turapov,<fnref n="1" />,†

Philippe Barthe,<fnref n="3" /> Heena Jagatia,<fnref n="4" />

Angélique De Visch,<fnref n="3" /> Christian Roumestand,<fnref n="3" />

Malgorzata Wegrzyn,<fnref n="5" /> Iona L. Bartek,<fnref n="6" />

Martin I. Voskuil,<fnref n="6" /> Helen M. O’Hare,<fnref n="1" />,<fnref n="7" /> Paul Ajuh,<fnref n="8" />

Andrew R. Bottrill,<fnref n="9" /> Adam A. Witney,<fnref n="10" />

Martin Cohen-Gonsaud,<fnref n="3" />,‡ Simon J. Waddell<fnref n="4" />,‡ and
Galina V. Mukamolova<fnref n="1" />,*,‡

## Summary

*Mycobacterium tuberculosis* (*Mtb*) is able to persist in the body through months of multi-drug therapy. Mycobacteria possess a wide range of regulatory proteins, including the protein kinase B (PknB) which controls peptidoglycan biosynthesis during growth. Here, we observed that depletion of PknB resulted in specific transcriptional changes that are likely caused by reduced phosphorylation of the H-NS-like regulator Lsr2 at threonine 112. The activity of PknB towards this phosphosite was confirmed with purified proteins, and this site was required for adaptation of *Mtb* to hypoxic conditions, and growth on solid media. Like H-NS, Lsr2 binds DNA in sequence-dependent and non-specific modes. PknB phosphorylation of Lsr2 reduced DNA binding, measured by fluorescence anisotropy and electrophoretic mobility shift assays, and our NMR structure of phosphomimetic T112D Lsr2 suggests that this may be due to increased dynamics of the DNA-binding domain. Conversely, the phosphoablative T112A Lsr2 had increased binding to certain DNA sites in ChIP-sequencing, and *Mtb* containing this variant showed transcriptional changes that correspond with the change in DNA binding. In summary, PknB controls *Mtb* growth and adaptations to the changing host environment by phosphorylating the global transcriptional regulator Lsr2.

## Introduction

*Mycobacterium tuberculosis* (*Mtb*) is a slow-growing bacterium that can replicate in humans and cause tuberculosis. The pathogen is able to rapidly shut down its growth to persist in non-replicating states in infected individuals, which can be modelled in the laboratory (Wayne and Sohaskey, 2001). *Mtb* adaptation to non-permissive conditions is accompanied by dramatic changes in global protein phosphorylation but the importance of these modifications is poorly defined (Prisic *et al.*, 2010).

<figure><img src="image_3.png" /></figure>

<footnote n="1">¹Leicester Tuberculosis Research Group, Department of
Respiratory Sciences, University of Leicester, Leicester,
LE2 9HN, UK.</footnote>

<footnote n="2">²Department of Basic Science, Faculty of Nursing,
University of Kufa, Najaf Governorate, P.O. Box 21,
Kufa, Najaf, Iraq.</footnote>

<footnote n="3">³Centre de Biochimie Structurale, CNRS, INSERM,
University of Montpellier, 34090, Montpellier, France.</footnote>

<footnote n="4">⁴Wellcome Trust Brighton and Sussex Centre for
Global Health Research, Brighton and Sussex Medical
School, University of Sussex, Brighton, BN1 9PX, UK.</footnote>

<footnote n="5">⁵Core Biotechnology Services, University of Leicester,
University Road, Leicester, LE1 7RH, UK.</footnote>

<footnote n="6">⁶Department of Immunology and Microbiology, University
of Colorado School of Medicine, Aurora, CO 80045,
USA.</footnote>

<footnote n="7">⁷LISCB, Department of Molecular and Cell
Biology, University of Leicester, University Road,
Leicester, LE1 7RH, UK.</footnote>

<footnote n="8">⁸Gemini Biosciences Ltd, Liverpool Science Park,
Liverpool, L7 8TX, UK.</footnote>

<footnote n="9">⁹Protein Nucleic Acid Laboratory, University of Leicester,
Leicester, LE1 7RH, UK.</footnote>

<footnote n="10">¹⁰Institute for Infection and Immunity, St George’s
University of London, London, SW17 0RE, UK.</footnote>

<footnote>Accepted 25 September, 2019. *For correspondence. E-mail gym4@le.ac.uk; Tel. +44(0)116 229 7053; Fax +44(0)1162525030. †Kawther Alqaseer and Obolbek Turapov should be considered joint first author.</footnote>

<footnote>‡Martin Cohen-Gonsaud, Simon J. Waddell and Galina V. Mukamolova should be considered joint senior author.</footnote>
*Mtb* has 11 serine/threonine protein kinases and they play significant roles in growth, virulence and metabolism (Richard-Greenblatt and Av-Gay, 2017). In particular, protein kinase B (PknB) is reported to be essential for *Mtb* growth (Fernandez *et al.*, 2006; Forti *et al.*, 2009) due to its critical function in the regulation of peptidoglycan biosynthesis (Gee *et al.*, 2012; Boutte *et al.*, 2016; Turapov *et al.*, 2018). It is also important for *Mtb* survival in hypoxic conditions and resuscitation during reaeration (Ortega *et al.*, 2014). However, the molecular mechanism for PknB-mediated adaptation to hypoxia is unknown.

We have recently shown that PknB-depleted *Mtb* can grow in osmoprotective sucrose magnesium medium (SMM) (Turapov *et al.*, 2018). Comparative phosphoproteomic analysis of PknB-producing against PknB-depleted mycobacteria revealed substantial changes. Specifically, the transcriptional regulator Lsr2 showed increased phosphorylation in PknB-producing mycobacteria, indicating that this protein may be a PknB substrate.

Lsr2 is a DNA-binding protein that combines the properties of a nucleoid-associated protein (Kriel *et al.*, 2018) and a global transcriptional regulator (Bartek *et al.*, 2014). Lsr2 has over 1000 binding sites in *Mtb* (Gordon *et al.*, 2010; Minch *et al.*, 2015). The precise role of Lsr2 in mycobacterial biology remains unclear, nevertheless parallels may be drawn with similar proteins from other bacteria. Lsr2 represents the first example of an H-NS-like protein identified outside Gram-negative bacteria; moreover, *lsr2* was able to complement an *hns* mutant in *Escherichia coli* (Gordon *et al.*, 2008). Similar to the H-NS proteins, Lsr2 has been proposed to bind to the minor groove of DNA (Gordon *et al.*, 2011) and to possess DNA bridging properties (Chen *et al.*, 2008). Additionally, Lsr2 has been shown to protect DNA from reactive oxygen species, and overexpression of Lsr2 improved survival of mycobacteria treated with hydrogen peroxide (Colangeli *et al.*, 2009). Deletion of *lsr2* in *Mtb* resulted in severe growth impairment on solid media, defects in persistence and adaptation to changing oxygen levels, all of which were accompanied by differential expression of genes involved in cell wall remodelling, respiration and lipid biosynthesis (Bartek *et al.*, 2014).

Here, we profiled the transcriptional changes that resulted from PknB depletion, and investigated the role of Lsr2 in coordinating these changes, as suggested by reduced phosphorylation of Lsr2 at T112 during PknB depletion (Turapov *et al.*, 2018). We probed the role of phosphorylation at this site in regulation of the structure and DNA-binding properties of Lsr2 and in governing growth and survival of *Mtb* in different conditions. Based on our data, we propose that PknB-mediated phosphorylation controls Lsr2 binding to DNA in *Mtb*, providing a functional link between serine/threonine protein kinase signalling in replicating bacilli and regulatory networks

that enable *Mtb* to survive dynamic environments during infection.

## Results

*Transcriptome profiling of PknB-depleted Mtb revealed an Lsr2-regulated gene expression signature*

PknB is essential for growth in standard conditions; however, we have recently developed an osmoprotective medium (SMM) that supported growth of PknB-depleted *Mtb* and allowed us to identify PknB substrates (Turapov *et al.*, 2018). Using the same system, we compared the transcriptional profile of PknB-depleted versus PknB-producing *Mtb* (Fig. 1A, Tables 1 and S1). PknB-depletion led to specific and significant changes in gene expression: 65 genes were induced and 34 repressed (Fig. 1B, Table S1). Two functional classes were overrepresented amongst the induced genes compared to the genome as a whole: regulatory proteins and proteins involved in lipid metabolism. The induced genes annotated as transcriptional regulators were csoR, rv1129c, rv1460, rv2017, rv2250c, rv3334, sigB, whiB3 and whiB6 (Table S1). These transcription factors regulate copper homeostasis (CsoR) (Marcus *et al.*, 2016), iron-sulphur cluster biogenesis (Rv1460) (Willemse *et al.*, 2018), cholesterol catabolism (Rv1129c/PrpR) (Masiewicz *et al.*, 2012), the enduring hypoxic response (Rv3334) (Rustad *et al.*, 2008), multiple stress responses (SigB) (Lee *et al.*, 2008), redox stress and complex lipid biosynthesis (WhiB3) (Mehta and Singh, 2018) and virulence factor expression (WhiB6) (Bosserman *et al.*, 2017).

The transcriptional signature of PknB depletion resembled features of intracellular growth (Table S1), with a significant overlap with RNA profiles from several studies of *Mtb* in macrophages as reflected by hypergeometric probability values: $6.7 \times 10^{-23}$ (Tailleux *et al.*, 2008), $7.34 \times 10^{-18}$ (Schnappinger *et al.*, 2003) and $3.57 \times 10^{-17}$ (Rohde *et al.*, 2007). For example, there was induction of pathways involved in mycobactin synthesis (*mbtB/C/D*), complex lipid phthiocerol dimycocerosate (PDIM) biosynthesis (*fadD26*, *ppsA/B/C/D*), metabolism of alternative lipid carbon sources, the glyoxylate shunt (*icl*), the methylcitrate cycle (*prpD/C*, *prpR*) and triacylglycerol synthase (*tgs1*). The isoniazid inducible genes (*iniB/A/C*) that respond to cell wall stress (Colangeli *et al.*, 2007), and four of the nine genes coding for alternative ribosomal proteins, *rpmB1*, *rpmB2*, *rpmG1*, *rpsN2* (Prisic *et al.*, 2015) were also induced.

The 34 genes that were significantly repressed in PknB-depleted bacteria included *pknB* itself (sixfold change, whereas no other protein kinases were significantly changed Table 1); *nuoA/B/C*, encoding subunits of NADH dehydrogenase I, which is part of the aerobic
respiratory chain, and several genes involved in inter-
mediary metabolism (Table S1). Comparison of gene
expression and protein abundance of selected targets
showed good agreement (Fig. 1C). Overall the number
of differentially expressed genes was comparable to the

number with differential expression when other regulators
were similarly disrupted, for example DosR (Park *et al.*,
2003). This is in contrast to the large-scale changes in
gene expression after treatment with an inhibitor of PknB
and PknA (Carette *et al.*, 2018), which would likely impact

<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Application of omics tools to characterise the function of *pknB* essential in standard growth media without affecting *Mtb* viability.</figcaption></figure>

A. Experimental set-up for sample preparation and analysis. Conditional PknB mutant (*pknB*-CM) was grown in sucrose-magnesium medium with pristinamycin (SMMpri, activation of PknB expression) or without pristinamycin (SMM, PknB depletion). Phosphoproteomics analysis was previously described (Turapov et al., 2018).

B. PknB depletion results in alteration of *Mtb* transcriptome. Transcriptional impact of PknB depletion in *Mtb*. Volcano plot showing 65 genes significantly induced (red) and 34 genes repressed (blue) by PknB-depletion in replicating *Mtb*. Significantly differentially expressed genes were identified using a moderated *t*-test (*P*-value < 0.05 with Benjamini and Hochberg multiple testing correction) and fold change >1.8 from three biological replicates. Fold change (log2) comparing *Mtb pknB*-CM in SMM with and without pristinamycin is plotted on the *x*-axis; corrected *p*-value (–log10) on the *y*-axis. Application of TFOE tool Lsr2 as potential regulator of the observed gene expression patterns Lsr2.

C. Expression of selected targets was validated using qRT-PCR and proteomics approaches (see for detail Experimental procedures).
Expression of *pknB*, *nuoB*, *hsp*, *csoR* was significantly different in SMM compared with SMMpri (*P* < 0.01, *t*-test).

<table><caption>**Table 1.** The impact of PknB depletion on the gene and protein expression levels of serine/threonine protein kinases and Lsr2 in *Mtb*.</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">Protein</th><th rowspan="2">Description</th><th colspan="2">Transcriptomics</th><th colspan="2">Proteomics</th></tr><tr><th>Fold difference SMMpri vs SMM<sup>a</sup></th><th>Fold difference SMMpri vs SMM</th><th>Fold difference SMMpri vs SMM</th></tr></thead><tbody><tr><td>Rv3597c</td><td>Lsr2</td><td>Protein Lsr2</td><td>1.01</td><td></td><td>1.20</td><td></td></tr><tr><td>Rv0015c</td><td>PknA</td><td>Protein kinase PknA</td><td>1.27</td><td></td><td>0.80</td><td></td></tr><tr><td>Rv0014c</td><td>PknB</td><td>Protein kinase PknB</td><td>5.97</td><td></td><td>6.60</td><td></td></tr><tr><td>Rv0931c</td><td>PknD</td><td>Protein kinase PknD</td><td>1.02</td><td></td><td>0.70</td><td></td></tr><tr><td>Rv1743</td><td>PknE</td><td>Protein kinase PknE</td><td>1.04</td><td></td><td>1.00</td><td></td></tr><tr><td>Rv1746</td><td>PknF</td><td>Protein kinase PknF</td><td>0.77</td><td></td><td>0.40</td><td></td></tr><tr><td>Rv0410c</td><td>PknG</td><td>Protein kinase PknG</td><td>0.93</td><td></td><td>1.30</td><td></td></tr><tr><td>Rv1266c</td><td>PknH</td><td>Protein kinase PknH</td><td>1.02</td><td></td><td>0.70</td><td></td></tr><tr><td>Rv2914c</td><td>PknI</td><td>Protein kinase PknI</td><td>0.84</td><td></td><td>N/D</td><td></td></tr><tr><td>Rv2088</td><td>PknJ</td><td>Protein kinase PknJ</td><td>0.89</td><td></td><td>N/D</td><td></td></tr><tr><td>Rv3080c</td><td>PknK</td><td>Protein kinase PknK</td><td>1.07</td><td></td><td>N/D</td><td></td></tr><tr><td>Rv2176</td><td>PknL</td><td>Protein kinase PknL</td><td>0.90</td><td></td><td>N/D</td><td></td></tr></tbody></table>

<footnote><sup>a</sup>Fold change values, comparing *Mtb pknB*-CM grown in sucrose magnesium medium with (SMMpri) or without pristinamycin (SMM) are derived from microarray and proteomics profiles.</footnote>
Mtb viability. In summary, PknB depletion in replicating bacteria resulted in co-ordinated changes to the transcriptome with similarities to intracellular adaptations, suggesting that PknB may control the induction of alternative gene regulatory pathways.

Application of the Transcription Factor Over-Expression (TFOE) output tool (Rustad *et al.*, 2014) predicted Rv0081 (Galagan *et al.*, 2013), DosR (Park *et al.*, 2003) and Lsr2 (Bartek *et al.*, 2014) as potential regulators of the observed gene expression patterns (Fig. 1, Table S1). We next focussed on the involvement of Lsr2 in PknB-mediated transcriptional adaptation, since Lsr2 was identified as a putative PknB substrate in our earlier phosphoproteomic work (Turapov *et al.*, 2018): PknB depletion decreased Lsr2 phosphorylation 2.54-fold without impacting Lsr2 protein expression levels (Turapov *et al.*, 2018; Table 1). By contrast, the

phosphorylation of DosR was unchanged during PknB depletion, and there are no reports of regulation of Rv0081 by phosphorylation.

### PknB phosphorylated Lsr2 in vitro

Kinase assays of purified PknB kinase domain with Lsr2 used anti-phosphothreonine antibody to detect phosphorylation and demonstrated that Lsr2 was directly phosphorylated by PknB (Fig. 2A). Interestingly, phosphorylation resulted in a marked change in Lsr2 protein mobility in SDS-PAGE (Fig. 2B) and generated several bands, indicative of multiple phosphorylated forms. Mass spectrometry confirmed the previously observed phosphosite on threonine 112 (Turapov *et al.*, 2018) and detected additional phosphorylations at threonine 8, threonine 22 and threonine 31 (Fig. 2C and D).

<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Identification of Lsr2 as a substrate of PknB.</figcaption></figure>

A. Recombinant Lsr2 was phosphorylated by recombinant PknB<sub>KD</sub>. Phosphorylated proteins were detected by western blot using a phosphothreonine antibody. Recombinant CwIM was used as positive control. 1 – protein markers; 2 – Lsr2 incubated with PknB and ATP; 3 – Lsr2 incubated with PknB without ATP; 4 – CwIM incubated with PknB and ATP; 5 – CwIM incubated with PknB without ATP.

B. SDS-PAGE revealed a shift in Lsr2 mobility upon phosphorylation (lanes identical to panel A).

C. Schematic presentation of phosphosites identified in phophoproteomics studies (top) (Turapov *et al.*, 2018), and *in vitro* (bottom).

D. Phosphopeptides detected by mass spectrometry; phosphorylated residues shown in bold font.
Phosphosite threonine 112 was necessary for Lsr2 function in Mtb

The functional importance of the identified phosphorylation sites in Lsr2 was further investigated by constructing a panel of phosphoablative *Lsr2* variants using the pMV306 plasmid that integrates at *attB* site of *Mtb* chromosome (Table S2) and measuring their ability to complement the phenotypic changes caused by *lsr2* deletion in *Mtb*. *Lsr2* deletion mutant containing the empty pMV306 plasmid (Δ*lsr2*<sub>pMV306</sub>) was used as a control. *Lsr2* deletion significantly impaired *Mtb* growth on solid media, similarly to a previous study (Bartek *et al.*, 2014). Expression of wild-type *lsr2* allele at *attB* site (the resultant strain designated as Δ*lsr2*<sub>WT</sub>) fully complemented the defect (Figs 3A and S1A). *Lsr2* deletion mutant expressing phosphoablative T112A Lsr2 variant (Δ*lsr2*<sub>T112A</sub>) had impaired growth on solid media (Figs 3A and S1A),

whereas all other phosphoablative variants (T8A, T22A, T31A) complemented growth fully (Fig. 3A). Growth of *Lsr2* deletion mutant expressing a T112D phosphomimetic variant of Lsr2 (Δ*lsr2*<sub>T112D</sub>) was indistinguishable from growth of wild-type *Mtb* (Fig. S1A). Notably, Δ*lsr2*<sub>pMV306</sub> and Δ*lsr2*<sub>T112A</sub> transformants were recovered from liquid medium, since these strains failed to produce colonies on solid media unlike transformation with plasmids carrying wild-type *Lsr2* and other variants in Fig. 3.

Comparison of the growth rates in liquid 7H9 medium (Fig. 3B) revealed a similar pattern: the phosphomimetic T112D variant grew at a similar rate to the strain with wild-type Lsr2 (Δ*lsr2*<sub>T112D</sub> 0.039 ± 0.009 h⁻¹ and Δ*lsr2*<sub>WT</sub> 0.041 ± 0.004 h⁻¹), whereas the phosphoablative T112A variant had a similar defect to the plasmid control (Δ*lsr2*<sub>T112A</sub> 0.017 ± 0.004 h⁻¹ and Δ*lsr2*<sub>pMV306</sub> 0.016 ± 0.03 h⁻¹). This failure of Lsr2 T112A to complement the growth defect of Δ*lsr2* in liquid or solid medium,

<figure><img src="image_6.png" /><figcaption>**Fig. 3.** Phosphoablative T112A variant does not complement growth and survival defect of *lsr2* deletion mutant.</figcaption></figure>

A. Lsr2 phosphoablative mutants were serially diluted and plated on 7H10 agar. Growth of *lsr2* deletion mutant expressing wild-type, T8A, T22A, T31A and T112A variants was compared with growth of the deletion mutant containing the empty vector on 7H10 agar. Experiment was repeated using two biological replicates.

B. *Mtb Δlsr2<sub>pMV</sub> Δlsr2<sub>WT</sub>, Δlsr2<sub>T112A</sub> and Δlsr2<sub>T112D</sub>* (~6 × 10⁵ cells/ml of each strain) were inoculated in 7H9 liquid medium supplemented with ADC and Tween 80 and incubated at 37°C without shaking.

C. Expression of Lsr2 and Lsr2 variants from pMV306 was verified by western blotting of *Mtb* lysates using an anti-Lsr2 antibody. Expression of GroEL was used as loading control.

D. T112A mutation impairs *Mtb* survival in the Wayne model of non-replicating persistence. *Mtb* Lsr2 mutants were incubated in sealed tubes with gentle mixing for up to 24 weeks. B and D Data presented as mean ± SEM (N = 6, two independent experiments done with biological triplicates).

***Statistically different in Δ*lsr2*<sub>pMV</sub> or Δ*lsr2*<sub>T112A</sub> compared with Δ*lsr2*<sub>WT</sub> and Δ*lsr2*<sub>T112D</sub>.
despite expression at similar levels (Fig. 3C), indicates a
requirement for phosphorylatable threonine at this posi-
tion for the function of Lsr2 in growing *Mtb*, suggesting
a possible mechanism by which PknB might regulate
transcription via phosphorylation of Lsr2. The phospho-
mimetic T112D Lsr2 had reduced mobility in SDS PAGE,
similar to phosphorylated Lsr2 (Figs 3C and 2B).

In separate experiments, we investigated whether expression of T112D Lsr2 from pMV306 would complement the growth defect of *pknB*-CM on solid SMM agar. However, introduction of pMV306::*lsr2*<sub>WT</sub> or pMV306::*lsr2*<sub>T112D</sub> did not improve growth of *pknB*-CM; moreover pMV306::*lsr2*<sub>T112A</sub> failed to produce any transformants, suggesting that T112A Lsr2 was toxic for *pknB*-CM.

*Lsr2 phosphorylation at threonine 112 is crucial for Mtb adaptation to hypoxic conditions but not in prolonged stationary phase*

We hypothesised that regulation of Lsr2 by PknB could
account for the defects in survival of oxygen depletion
reported when *pknB* and *lsr2* function was abrogated
(Bartek *et al.*, 2014; Ortega *et al.*, 2014). Thus, we
assessed survival of *Mtb* carrying Lsr2 variants using the
Wayne model of non-replicating persistence in hypoxia
(Wayne and Sramek, 1994). The viable counts of Δ*lsr2*<sub>WT</sub>
and Δ*lsr2*<sub>T112D</sub> increased after 7 weeks of incubation,
proving that *Mtb* grew during gradual depletion of oxygen.
In opposite, Δ*lsr2*<sub>pMV306</sub> showed a dramatic drop of CFU
counts below the initial inoculum. Further incubation for
24 weeks resulted in near-complete loss of Δ*lsr2*<sub>pMV306</sub>
bacteria but did not significantly alter the survival of
Δ*lsr2*<sub>WT</sub> and Δ*lsr2*<sub>T112D</sub>. Interestingly, T112A Lsr2 was
apparently unable to grow during initial 7 weeks, suggest-
ing a requirement for T112 phosphorylation during adap-
tation for decreasing oxygen levels.

We also assessed the survival of *Mtb* expressing Lsr2 variants in late stationary phase of aerobically grown cultures by MPN and CFU counting. Δ*lsr2*<sub>pMV306</sub> showed impaired survival compared to Δ*lsr2*<sub>WT</sub>, resulting in a 1.5 order of magnitude difference in viable counts. In this model, the survival of T112A or T112D variants (Δ*lsr2*<sub>T112A</sub>, Δ*lsr2*<sub>T112D</sub>) was not significantly different from wild-type Lsr2 (Fig. S1B). Our results suggest that phosphorylation of Lsr2 at T112 may be specifically required during adaptation to hypoxic conditions but not for survival in prolonged stationary phase.

*Increased DNA binding in the T112A Lsr2 variant of Mtb-altered gene expression*

Phosphorylation of nucleoid-associated proteins is known
to influence their interaction with DNA (Dilweg and Dame,
2018), thus providing a mechanism by which PknB might

regulate transcription via Lsr2. The DNA-binding profile of Lsr2, and the influence of phosphorylation upon its DNA-binding profile, were investigated using ChIP-seq and a custom anti-Lsr2 antibody. Lsr2-binding peaks were found at intergenic spaces and running through several open reading frames, and were well conserved across biological replicates (Fig. S2 and Table S3), whereas no DNA was precipitated from strains lacking *lsr2*, confirming antibody specificity. The putative regulon of Lsr2, including genes with a Lsr2-binding site within or immediately upstream of the coding sequence, was defined as 1178 genes (Table S3), and was consistent with previously identified Lsr2-binding patterns (Gordon *et al.*, 2010; Minch *et al.*, 2015) (Table S4). The regulon was significantly enriched in genes that are differentially expressed upon inactivation of Lsr2 (Bartek *et al.*, 2014), macrophage infection (Tailleux *et al.*, 2008), in sputum (Garton *et al.*, 2008) and under acid-nitrosative stress (Cossu *et al.*, 2013; Table S4), overlapping with the transcriptional signature of PknB depletion.

To measure the effect of phosphorylation of Lsr2 on its DNA-binding pattern, we compared the ChIP-seq data above with parallel experiments using phosphomimetic Lsr2 (Δ*lsr2*<sub>T112D</sub>) and phosphoablative Lsr2 (Δ*lsr2*<sub>T112A</sub>). The results for phosphomimetic strain were essentially the same as those with wild-type Lsr2 (Δ*lsr2*<sub>WT</sub> *Mtb*, no significant differences in sequence abundance), whereas DNA precipitated from the phosphoablative strain showed a significant increase in the abundance of 226 putative Lsr2-binding sites (Fig. S2), suggesting that greater DNA binding of this variant. These binding sites may affect the expression of 94 genes (Table S5), which are involved in pathways of cell wall biosynthesis, lipid metabolism, PE/PPE protein synthesis and intermediary metabolism. Application of Motif-based sequence analysis predicted a binding site where T122A variant was preferably binding (Fig. 4A).

To determine whether altered binding of Lsr2 variants resulted in corresponding changes in RNA abundance, we tested two of these genes for differential expression in the presence/absence of Lsr2 and in the presence of wild-type versus phosphoablative Lsr2. *sodA* and *leuS* were selected for investigation on the basis of their divergent patterns of Lsr2 binding (Fig. 4B). *sodA* was upregulated but *leuS* did not change significantly in Δ*lsr2*<sub>WT</sub> compared with Δ*lsr2*<sub>pMV306</sub> (Fig. 4C). Both genes were differentially expressed in Δ*lsr2*<sub>T112A</sub> versus Δ*lsr2*<sub>WT</sub>, but *leuS* expression was lower, while *sodA* expression was higher (Fig. 4C). These results suggest that the influence of phosphorylation on Lsr2 DNA binding can up- or downregulate gene expression, depending on promoter structure and other regulatory factors. Genes known to be differentially expressed upon *lsr2* disruption were included alongside (*rpfA* and *rpfC*) (Bartek *et al.*, 2014), as well as the principle sigma factor gene *sigA*, which has a putative Lsr2-binding site in its promoter (Rustad *et al.*, 2014). However, expression of these
<figure><img src="image_7.png" /><figcaption>**Fig. 4.** T112A mutation alters Lsr2 binding to DNA and gene expression patterns.</figcaption></figure>

A. Lsr2 T112A binding site was predicted using the MEME Suite (http://meme-suite.org) and sequence of 25 DNA fragments which were >6-fold enriched in Lsr2 T112A compared with Lsr2 WT.

B. Representative plots describe Lsr2 binding upstream of *rpfA*, intragenic binding in *leuS*, or both intergenic and intragenic binding in *sodA*, showing greater Lsr2 binding in three biological replicates of phosphoablative Δ*lsr2*<sub>T112A</sub> (green) compared to Δ*lsr2*<sub>WT</sub> (yellow). A black asterisk marks the position of a Lsr2-binding site in *leuS*. Plots adapted from Integrative Genomics Viewer IGV, (Robinson et al., 2011). Binding patterns in Δ*lsr2*<sub>T112D</sub> were identical to those in Δ*lsr2*<sub>WT</sub> and not shown for clarity.

C. Expression of *leuS*, *rpfA*, *rpfC*, *sigA* and *sodA* relative to Δ*lsr2*<sub>pMV306</sub> measured by quantitative RT-PCR in Δ*lsr2*<sub>WT</sub> and Δ*lsr2*<sub>T112A</sub> and normalised to 16s rRNA and Δ*lsr2*<sub>pMV306</sub>. Data presented as mean ± SEM (N = 6). **Statistically different in Δ*lsr2*<sub>T112A</sub> compared with Δ*lsr2*<sub>WT</sub> (P < 0.01).

genes did not significantly change in Δ*lsr2*<sub>T112A</sub> versus Δ*lsr2*<sub>WT</sub> in accordance with our ChIP-sequencing data.

*PknB-mediated phosphorylation of Lsr2 or phosphomimetic mutation of Lsr2 reduced its DNA binding in vitro*

To determine how phosphorylation of Lsr2 affects the affinity and sequence specificity of DNA binding, we performed electrophoretic mobility shift assays (EMSA) with a range

of DNA fragments containing putative Lsr2-binding sites. These fragments were obtained by PCR or by annealing pairs of oligonucleotides to produce double-stranded DNA (Table S2). Lsr2 reduced the mobility of all fragments when added at concentrations above 1.9 µM and there was no apparent difference in affinity between the DNA fragments (Fig. S3). The same pattern was observed when a shorter double-stranded DNA containing a putative Lsr2-binding site within the *leuS* gene (Reddy *et al.*, 2009) was used (Fig. 5A). We tested various DNA fragments, including a mutated *leuS*
<figure><img src="image_8.png" /><figcaption>**Fig. 5.** Phosphorylated Lsr2 or its phosphomimetic T112D variant do not bind DNA.</figcaption></figure>

A. Lsr2 was mixed with *leuS* fragment containing a putative binding site (AATTCGGCAAAATCGGTAAG), position of which is marked with an asterisk in Fig. 4A.

B. Lsr2 was mixed with the mutated *leuS*<sub>MUT</sub> fragment (AACTCGGCGAGGTCGGTCAG).

C. Lsr2 was phosphorylated by PknB and mixed with *leuS*-binding site.

D. Lsr2 T112D variant was mixed with *leuS* fragment. Lsr2 was added to DNA at a range of concentrations (0.95–7.6 µM). Representative results from three independent experiments shown.

E. Quantification of Lsr2<sub>BD</sub> (DNA-binding domain) interaction with DNA by fluorescence anisotropy. Titration of 5' Alexa Fluor 488 double-stranded DNA (CGCGCATATATGCG) (4 nM) by Lsr2<sub>BD</sub> WT (open circles), Lsr2<sub>BD</sub> T112A (black circles) and Lsr2<sub>BD</sub> T112D (black diamonds).

F. Control experiments with GC-rich double-stranded DNA fragment (AACTCGGCGAGGTCGGTCAG).

site (Fig. 5B), and all were shifted by Lsr2 regardless of their sequence, demonstrating that Lsr2 binds non-specifically to DNA, as previously suggested (Colangeli *et al.*, 2007).

We next assessed the effect of Lsr2 phosphorylation on binding to these DNA fragments. As shown in Fig. 5C, phosphorylation of Lsr2 by PknB completely abolished binding (Figs 5C and S3), whereas Lsr2 T112D variant showed reduced DNA binding (Fig. 5D). This recalls the ChIP-sequencing results where Lsr2 T112D and wild-type Lsr2 (that is phosphorylated in cells), precipitated less DNA than phosphoablative T112A Lsr2.

*In vitro* phosphorylation of Lsr2 resulted in phosphorylation of both domains. Like Lsr2, HN-S proteins also consist of an oligomerisation domain and DNA-binding domain.

In the HN-S family, the two domains have distinct roles: shaping the nucleoid and regulating gene expression respectively (Winardhi *et al.*, 2015). In the case of Lsr2, threonine 112 in the DNA-binding domain was required for function in *Mtb*, whereas other putative phosphorylation sites in the oligomerisation domain were mutated without loss of function (Fig. 3). To discriminate between the effects of phosphorylation on these two domains, we purified the DNA-binding domain (Lsr2<sub>BD</sub>) to compare DNA binding of phosphorylated and unphosphorylated protein. A fluorescence anisotropy approach was used, since truncated Lsr2<sub>BD</sub> was not suitable for EMSA.

Previous studies demonstrated that Lsr2<sub>BD</sub> mainly recognised AT-rich DNA sequences that formed a
hook-like structure (Gordon et al., 2011). We measured
to ability of Lsr2<sub>BD</sub> and T112 mutants to bind an Alexa
488N labelled AT-rich double-stranded DNA fragment
(5'-CGCGCATATATGCG-3') (Fig. 5E). At pH 7.5 the Kd
value for Lsr2<sub>BD</sub> was 0.21 ± 0.06 µM, that for Lsr2<sub>BD</sub>T112A
was 0.02 ± 0.008 µM and Lsr2<sub>BD</sub> T112D showed no sig-
nificant DNA binding (Kd could not be determined). In
control experiments, Lsr2<sub>BD</sub> or Lsr2<sub>BD</sub>T112D showed no
binding to a GC-rich DNA sequence (Fig. 5F), confirming
Lsr2<sub>BD</sub> preferential binding to AT-rich DNA. These results
show that the T112D mutation, mimicking the Lsr2 phos-
phorylated state, reduced Lsr2<sub>BD</sub>T112D binding to DNA *in vitro*.
In summary, phosphorylation controls both Lsr2 and
Lsr2<sub>BD</sub> binding to DNA.

*Phosphomimetic T112D variant changed the*
*conformation entropy of the Lsr2 DNA-binding domain*

To elucidate a molecular basis for the reduced binding
of phosphomimetic Lsr2 to DNA we compared nuclear
magnetic resonance (NMR) structures of Lsr2<sub>BD</sub>WT and
Lsr2<sub>BD</sub>T112D. In accordance with the previously pub-
lished structure (Gordon et al., 2010), our data confirmed
that Lsr2<sub>BD</sub> consists of two perpendicular α-helices
(α1, residues 78–89; α2, residues 102–112) linked by a
long loop (residues 90–101). The two major components
involved in Lsr2 to DNA binding are residues Arg97–
Gly98–Arg99 that are inserted into the minor grove of
DNA and Arg77, Ser 80, Arg84 and Ser95 that interact
with the phosphate-sugar backbone of the minor grove
(Gordon et al., 2010). This organisation was preserved
in Lsr2<sub>BD</sub>T112D (Table S6). However, T112D mutation
resulted in a shorter α2 helix (Fig. S4), which ended with
alanine 111 in Lsr2<sub>BD</sub>T112D compared with threonine
112 in Lsr2<sub>BD</sub>WT (Fig. 6A and B). Both variants were
monomeric and the N-terminal segment (residues 66-75)

upstream of the DNA-binding domain was disordered. We
also found that in Lsr2<sub>BD</sub> the methyl group of threonine
112 interacted with tyrosine 108, while the hydroxyl group
of threonine 112 interacted with tryptophan 86 (Fig. 6C).
These interactions did not form in Lsr2<sub>BD</sub>T112D, likely
accounting for the shorter helix. Lsr2<sub>BD</sub> T112A, which
must lack both interactions, bound DNA in our anisotropy
experiments (Fig. 5E), demonstrating that the interactions
themselves are not required for DNA binding, while sug-
gesting that changes in conformation or dynamics related
to the shorter helix could account for changes in DNA
binding.

<sup>15</sup>N heteronuclear NMR relaxation analysis was per-
formed to assess the dynamic behaviour of the two pro-
teins (Fig. S5). In both isoforms, two α helices showed
similar amplitudes for internal motion. However, in
Lsr2<sub>BD</sub> mainly the C-terminal helix was affected, while
in Lsr2<sub>BD</sub>T112D the internal motion was extended to
the N-terminal helix. A previous study on the catabolite
activator protein (CAP) demonstrated that different pro-
tein mutants with the same structure of interaction inter-
face displayed very different affinity for their target DNA
(Tzeng and Kalodimos, 2012). The authors showed
that changes of the binding affinity were linked to fast
internal dynamics (conformation entropy). Similarly, the
T112D mutation and, presumably, phosphorylation of
threonine 112, resulted in a shorter α2 helix and a more
mobile loop that increased the Lsr2<sub>BD</sub> dynamics and
impaired DNA binding.

Discussion

*Mtb* can subvert the immune system to survive in the host for many years. This remarkable ability is determined by mechanisms that allow *Mtb* to respond to multiple environments and adjust metabolic activity and cell division.

<figure><img src="image_9.png" /><figcaption>**Fig. 6.** Solution structure of Lsr2<sub>BD</sub> (A) and Lsr2<sub>BD</sub>T112D (B). Superimposition of the 20 best calculated structures in cartoon representation with the last residues represented as sticks (A and B). α1 helix and linker region are involved in DNA binding; the structure of these regions is not affected by T112D mutation. In Lsr2<sub>BD</sub>, the threonine T112 side chain interacts with both tyrosine Y108 (3.7 Å) and tryptophan W86 (3.0 Å), all represented in sticks (C).</figcaption></figure>
One of these regulatory systems is protein phosphorylation involving 11 serine/threonine kinases, including protein kinase B. PknB is expressed during active replication (Kang et al., 2005) and is essential for growth in standard media (Fernandez et al., 2006; Forti et al., 2009). PknB has been implicated in the regulation of peptidoglycan biosynthesis (Bellinzoni et al., 2019); its main substrate CwIM stimulates biosynthesis of peptidoglycan precursors (Boutte et al., 2016) and may facilitate their transport to the cell surface (Turapov et al., 2018). This regulation is crucial for adjusting bacterial growth and synthesis of the cell wall. The external PASTA domain of PknB (Barthe et al., 2010; Prigozhin et al., 2016) is essential for PknB-mediated signalling and its disruption results in bacterial death and alteration of antimicrobial susceptibility (Chawla et al., 2014; Turapov et al., 2015; Turapov et al., 2018). Furthermore, PknB can phosphorylate sigma factor SigH and its cognate anti-sigma factor RshA; phosphorylation disrupts interaction of these two proteins and result in increased expression of the SigH regulon (Park et al., 2008).

Here, we present data demonstrating that PknB also controls global gene expression via another substrate, the DNA-binding protein Lsr2. Transcriptomic analysis of PknB-depleted *Mtb* suggest that PknB phosphorylation may silence alternative pathways, which are not important for logarithmic growth but which may be critical for mediating stress responses and virulence. These include enzymes involved in alternative metabolic pathways, synthesis of complex lipids, regulators of stress responses and antimicrobial tolerance (Fig. 1, Table S1). PknB depletion resulted in the altered expression of several transcriptional regulators and the regulons of Rv0081, DosR and Lsr2. However only Lsr2 was more phosphorylated in PknB-producing *Mtb* compared with PknB-depleted *Mtb*. According to our transcriptomic and proteomic signatures, none of the other serine/threonine protein kinases were upregulated in these conditions (Table 1), suggesting that PknB was responsible for Lsr2 phosphorylation in growing *Mtb*.

We demonstrate that Lsr2 could be phosphorylated by PknB at several threonines but only the T112 was essential for growth and survival of *Mtb*. The phosphoablative T112A mutant could not complement the growth defect of the *Mtb* *lsr2* deletion mutant and was impaired in the Wayne hypoxia model. To investigate the molecular mechanism responsible for the growth defect in Δ*lsr2*<sub>T112A</sub> we conducted ChIP-Seq analysis and compared DNA-binding patterns in Δ*lsr2*<sub>T112A</sub>, Δ*lsr2*<sub>T112D</sub> and Δ*lsr2*<sub>WT</sub>. In agreement with previously published data (Gordon et al., 2010; Minch et al., 2015) we detected multiple occurrences of DNA binding in all three Lsr2 backgrounds supporting both nucleoid shaping and gene regulatory functions of Lsr2. However, phosphoablative Lsr2 T112A variant

had increased DNA binding, potentially directly affecting expression of 94 genes. Interestingly, most of these genes were not essential for growth *in vitro* (DeJesus et al., 2017). Moreover, deletion of some of these genes has been previously shown to be advantageous for growth (DeJesus et al., 2017), including genes of unknown function (Rv0888, Rv1958, Rv1957), PE/PPE genes (Rv0878, Rv1983), *aprA* and genes controlling transport of PDIM (*drrB* and *drrC*). While products of these genes might be disadvantageous for growth *in vitro*, they likely play a critical role in adaptations to stress and virulence (Camacho et al., 2001; Sassetti and Rubin, 2003; Rohde et al., 2007).

Our data demonstrate that PknB phosphorylation of Lsr2 *in vitro* completely abolished DNA binding, while the phosphomimetic mutation reduced Lsr2–DNA interactions (Fig. 5). We have not investigated the effect of T8, T22 and T31 phosphorylations on DNA binding. Based on previously published data we hypothesise that phosphorylation of these sites in the oligomerisation domain might be important for controlling nucleoid shape and DNA bridging properties. The results of our study suggest that phosphorylation of T112 in the DNA-binding domain controls interaction of Lsr2 with AT-rich DNA sequences which modifies gene expression (Fig. 4C). Our structural studies further confirm that phosphomimetic T112D Lsr2 variant had a shorter C-terminal helix and increased dynamics of the DNA-binding domain, leading to impaired Lsr2–DNA binding (Figs 6 and S5).

Post-translational modifications are common mechanisms for the regulation of DNA binding both in eukaryotes (Bannister and Kouzarides, 2011) and prokaryotes (Dilweg and Dame, 2018). Phosphorylation or nitrosylation of transcriptional regulators abolish DNA binding (Leiba et al., 2014; Smith et al., 2017). H-NS protein, a homologue of Lsr2 in *E. coli*, has been shown to be acetylated and phosphorylated; however, the precise function of these modifications remains to be characterised (Dilweg and Dame, 2018). Our data show that phosphorylation of Lsr2 is important for *Mtb* growth and that this may be a key mechanism for controlling mycobacterial adaptations to permissive and non-permissive environments. Thus, PknB mediates two critical components of mycobacterial growth, peptidoglycan biosynthesis and gene expression of alternative pathways.

Based on our data, we propose that PknB controls *Mtb* growth by phosphorylating Lsr2. Like other H-NS-like proteins Lsr2 plays a dual role in mycobacterial biology, it shapes and protects the nucleoid and it controls gene expression (Bartek et al., 2014; Kriel et al., 2018). However, unlike H-NS proteins in Gram-negative bacteria that mainly silence the expression of foreign DNA (Lucchini et al., 2006), Lsr2 regulates expression of genes that are essential for growth, virulence and adaptation (Bartek et al., 2014). Our study suggests that
phosphorylation of T112 might be important for tuning gene expression during growth, and the dynamic change between phosphorylated and non-phosphorylated Lsr2 may help to adjust transcriptional patterns according to growth conditions. Reduced T112 phosphorylation, for example during starvation, may increase Lsr2 binding and upregulate pathways that are critical for Mtb survival under these conditions. Our data suggest that PknB is the main serine/threonine kinase responsible for phosphorylation of Lsr2 at T112 during growth, however we cannot exclude that other kinases can phosphorylate Lsr2 at this or other sites under different conditions as it has previously been shown for other substrates (Baer *et al.*, 2014).

While there are many outstanding questions on the precise mechanisms of PknB-mediated regulation of gene expression and Lsr2 binding to DNA, our findings provide a functional link between serine/threonine protein kinase signalling and transcriptional regulatory pathways that enable Mtb to survive the varied environments encountered during infection.

## Experimental procedures

### Strains and media

Mtb H37Rv was grown in Middlebrook 7H9 (Becton, Dickinson and Company) liquid medium supplemented with 10% (v/v) Albumin-Dextrose Complex (ADC), 0.2% (v/v) glycerol and 0.1% (w/v) Tween 80 at 37°C with shaking or in SMM containing hygromycin with or without pristinamycin. SMM comprised of 0.3 M sucrose, 20 mM MgSO₄·4H₂O, 0.1% Tween 80 (w/v), 10% (v/v) ADC in standard 7H9 broth. Bacterial growth was measured by absorbance at 580nm, or by colony-forming unit (CFU) counting on 7H10 agar (Becton, Dickinson and Company), or by most probable number (MPN) counting using established protocols (Loraine *et al.*, 2016) and the MPN calculator program (Jarvis *et al.*, 2010). *Escherichia coli* OverExpress™ C41(DE3) and DH5α were grown in Lysogeny broth. For protein expression, *E. coli* was grown to mid-log phase (OD₆₀₀ 0.6-0.8) at 37°C with shaking at 200 rpm before adding 0.5 mM isopropyl β-D-1-thiogalactopyranoside followed by incubation at 18°C overnight. Antibiotics were used at the following concentrations (µg/ml): pristinamycin, 0.5; kanamycin, 50; hygromycin, 50; ampicillin, 50. Wayne model of non-replicating persistence was set up as previously described (Wayne and Sramek, 1994). CFU and MPN counts were determined at 0-, 7- and 12-week time points.

### Genetic manipulations

Previously described *pknB*-CM (Forti *et al.*, 2009) and *Δlsr2* mutants (Bartek *et al.*, 2014) were used in this study. In *Δlsr2* mutant, a gene fragment encoding the C-terminal DNA-binding domain (corresponding to 174–268 bp region) was deleted. The *Rv3597* (*lsr2*) coding sequence with additional

200 bp upstream region containing the putative promoter was amplified from the *Mtb* H37Rv genome using Platinum Taq-HF polymerase (ThermoFisher) and cloned into the pMV306 plasmid that integrates at *attB* site of *Mtb* chromosome (primers in Table S2). Lsr2 *Mtb* variants were obtained using the GeneArt™ Site-Directed Mutagenesis System (Thermo Fisher Scientific) according to the manufacturer's instructions. All constructs were sequenced before transformation into an *Mtb* *lsr2* deletion mutant. Transformants were selected on 7H10 agar or in supplemented 7H9 containing kanamycin. For expression of Lsr2 proteins *Mtb* H37Rv *lsr2*-coding sequence or shorter fragments were PCR amplified and cloned into pET15-TEV and the resultant constructs were transformed into *E. coli* OverExpress™ C41 (DE3). Bacterial strains generated in this study are shown in Table S2.

## Preparation of samples for proteomics and transcriptomics

Mtb *pknB*-CM was grown to OD₅₀₀nm ~0.7 in SMM with or without pristinamycin. For proteomics studies, washed *pknB*-CM bacteria were resuspended in buffer containing 20 mM TrisCl, pH 7.5, 1 M NaCl, 8 M urea and proteinase/ phosphatase inhibitors. After bead beating, lysates were cleared by centrifugation and filtration and treated using the FASP protocol, as described previously (Turapov *et al.*, 2018; Iswahyudi *et al.*, 2019). For transcriptomics analysis, bacterial cultures were incubated with four volumes of guanidine thiocyanate (GTC) for 30 minutes prior centrifugation.

## Quantitative label-free proteomics analysis

Analysis was performed as previously described (Turapov *et al.*, 2018; Iswahyudi *et al.*, 2019). For quantification, all peptides of an identified protein were included and the total cumulative abundance was calculated by summing the abundance of all peptides allocated to the respective protein. Additionally, Scaffold Q+ (version Scaffold_4.3.4, Proteome Software Inc) was used for peptide and protein identifications as previously described (Turapov *et al.*, 2015). Peptide probabilities from Mascot were assigned by the Scaffold Local FDR algorithm.

## Transcriptomic analyses

RNA was extracted from three biological replicates using the GTC/Trizol method (Waddell and Butcher, 2010). *Mtb* RNA (2 µg) was enzymatically labelled with Cy3 fluorophore and hybridised to a *Mtb* complex microarray (ArrayExpress accession number A-BUGS-41) as previously described (Salina *et al.*, 2014). For quantitative RT-PCR, total RNA was isolated from triplicate *Mtb* cultures; cDNA was generated using Superscript Reverse Transcriptase II and mycobacterial genome-directed primers (Rachman *et al.*, 2006). qPCR was performed in a Corbett Rotor Gene 6000 real-time thermocycler using Absolute qPCR SYBR Green mix and gene expression values were normalised to 16S rRNA expression. For qPCR three biological and two replicates were assessed.
## ChIP-Seq analysis

DNA-Lsr2 interactions in Δ*lsr2*<sub>T112A</sub>, Δ*lsr2*<sub>T112D</sub>, Δ*lsr2*<sub>WT</sub> were assayed using ChIP-seq methods as previously described (Minch *et al.*, 2015). Briefly, mid log-phase *Mtb* H37Rv cultures (OD<sub>600</sub> 0.4–0.6) were crosslinked with 1% formaldehyde, followed by incubation with 125 mM glycine for 5 minutes at 37°C. The cells were mechanically lysed and then sonicated to produce 200–500 bp fragments. Input control samples were taken for each genotype before antibody was added to assess antibody specificity. Samples were immunoprecipitated using a polyclonal anti-rabbit anti-Lsr2 antibody and protein-G agarose beads. The Lsr2 complexes were de-crosslinked by heating at 65°C overnight and proteins removed by treatment with proteinase K (10 mg/ml) for 2 h at 55°C. The DNA samples were column-purified (Qiagen) and the quality of purified IP-Lsr2 DNA verified using the Qubit DNA HS quantification assay and Nanodrop spectrophotometer. Libraries were prepared and sequenced using Illumina HiSeq SE50, 20 million reads (Novogene, Hong Kong). Raw fastq files were aligned to the *Mtb* H37Rv (NC_000962.3) reference genome using bwa samse (Li and Durbin, 2009). MACS2 version 2.1.1.20160309 (Zhang *et al.*, 2008) was used to compare each of the input controls to the immunoprecipitated samples, identifying Lsr2-binding sites (callpeak) using default parameters but including -g 4.41e+06 -nomodel -extsize 147' (Table S3). Differential peaks comparing Δ*lsr2*<sub>T112A</sub> to Δ*lsr2*<sub>WT</sub> were then identified using 'macs2 bdgdiff'.

## Purification of recombinant proteins and phosphorylation in vitro

Recombinant Lsr2 proteins were purified using immobilised metal affinity chromatography (Ni-NTA agarose, Qiagen) and size exclusion chromatography. The recombinant catalytic domain of PknB was purified using Glutathione Sepharose 4B GST-tagged protein purification resin (GE Healthcare). For phosphorylation recombinant Lsr2 (10 µM) was mixed with the recombinant catalytic domain of PknB (5 µM) in a kinase buffer (20 mM Tris-HCl, pH 8.0; 0.5 mM DTT; 10 mM MgCl₂; 0.2 mM ATP) and incubated at 37°C for 1h. Phosphorylation was confirmed by western blot analysis. Phosphorylated residues were identified in trypsin-digested proteins using LTQ-Orbitrap-Velos mass spectrometer.

## Protein electrophoresis and western blotting

Proteins were separated on 4–20% gradient SERVA gels and transferred onto a nitrocellulose membrane using a Trans-Blot® Turbo™ Transfer System (Bio-Rad). SIGMAFAST™ BCIP®/NBT or SignalFire™ Elite ECL Reagent were used to visualise proteins on C-DiGit Chemiluminescent Blot Scanner (LI-COR Biosciences). The following antibodies were used: custom polyclonal antibody raised against Lsr2 in rabbit (Gemini Biosciences); monoclonal murine anti-polyhistidine antibody (Sigma-Aldrich); phospho-threonine antibody (Cell Signaling Technology); monoclonal anti-*MtbGroEL2* (*Rv0440*), clone IT-70 (BEIResources); mouse anti-rabbit IgG antibody:alkaline phosphatase (Sigma-Aldrich; anti-mouse IgG (whole

molecule:alkaline phosphatase antibody produced in rabbit (Sigma-Aldrich), and anti-rabbit IgG, HRP-linked antibody (Cell Signaling Technology).

## Electrophoretic mobility shift assay

Electrophoretic mobility shift assays were carried out with DNA fragments amplified by PCR or annealed oligonucleotides (see Table S2 for primer details). DNA (1.2 nM) was mixed with indicated amounts of Lsr2 in a total volume of 20µl reaction buffer containing (10 mM Tris-HCl, pH 7.5, 50 mM KCl, 1 mM DTT, 5 mM MgCl₂, and 2.5% glycerol). The mixture was incubated for 30 min at room temperature followed by native polyacrylamide gel electrophoresis using 8% gels in 0.5 x Tris-Borate-EDTA buffer, pH 7.5 for 24 min at 120 V. The gels were stained with SYBR Safe DNA stain (Thermo Fisher Scientific) and visualised using a ChemiDoc system (Bio-Rad).

## Fluorescence anisotropy

Custom made 5' Alexa Fluor 488 succinimidyl ester labelled oligonucleotide probe (sequence 5'-CGCATATATGCGCG-3') was purchased from Integrated DNA Technologies. Steady-state fluorescence anisotropy-binding titrations were performed on a Tecan Saphire II microplate reader, using a 470 nm LED for excitation and a monochromator set at 530 nm (bandwidth 20 nm) for emission.

## Determination of solution structures and dynamics of Lsr2<sub>BD</sub> and Lsr2<sub>BD</sub>T112D

All ¹H-¹⁵N double-resonance NMR experiments were performed at 20°C on Bruker Avance III spectrometers (700 or 800 MHz) using previously described methods (Barthe *et al.*, 1999; Gordon *et al.*, 2010). NMR samples of 0.5mM ¹⁵N-labelled protein dissolved in 25mM sodium phosphate buffer (pH 6.8), 150mM NaCl with 10% D₂O for the lock. ¹H chemical shifts were directly referenced to the methyl resonance of DSS, while ¹⁵N chemical shifts were referenced indirectly to the absolute ¹⁵N/¹H frequency ratio. All NMR spectra were processed with GIFA (Pons *et al*, 1996). Chemical shift assignments were made using standard NOESY, TOCSY experiments performed on the ¹⁵N-labelled protein sample. NOE cross-peaks identified on 3D [¹H, ¹⁵N] NOESY-HSQC (mixing time 160 ms) were assigned through automated NMR structure calculations with CYANA 2.1 (Güntert, 2004). Backbone φ and ϕ torsion angle constraints were obtained from a database search procedure on the basis of backbone (¹⁵N, HN, Hα) chemical shifts using TALOS+ (Shen *et al.*, 2009). For each protein, a total of 200 three-dimensional structures were generated using the torsion angle dynamics protocol of CYANA 2.1. The 20 best structures of each protein (based on the final target penalty function values) were minimised with CNS 1.2. All statistical parameters are summarised in (Table S6). Relaxation rate constant measurements were performed on a 0.5mM protein sample, at 18.8 T (800 MHz). The pulse sequences used to determine ¹⁵N R<sub>N</sub>(N<sub>z</sub>) (R<sub>1</sub>), R<sub>N</sub>(N<sub>xy</sub>) (R<sub>2</sub>), and ¹⁵N{¹H} NOE values were similar to
those described (Barthe *et al.*, 1999). The ¹⁵N longitudinal relaxation rates ($R_N(N_Z)$) were obtained from 10 standard inversion recovery experiments, with relaxation delays ranging from 18 to 1206 ms. The ¹⁵N transverse relaxation experiments ($R_N(N_{xy})$) were obtained from 10 standard CPMG experiments, with relaxation delays ranging from 16 to 160 ms. Both series of experiments were acquired in two single interleaved matrices to ensure uniformity of the experimental conditions. Heteronuclear ¹⁵N{¹H} NOE were determined from the ratio of two experiments, with and without saturation.

### Statistical analysis

Calculation of the protein p-values was performed on the sum of the normalised abundance across all runs using one-way ANOVA. Significantly differentially expressed genes (Table S1) were identified using a moderated *t*-test (*P*-value < 0.05 with Benjamini and Hochberg multiple testing correction), and fold change> 1.8 in GeneSpring 14.5 (Agilent Technologies). Hypergeometric probability and TFOE analysis (Rustad *et al.*, 2014) were used to identify significantly enriched signatures. An unpaired *t*-test was performed to compare gene expression in $\Delta lsr2_{WT}$ and $\Delta lsr2_{T112A}$. For evaluation of growth parameters and survival in the Wayne model, one-way ANOVA (GraphPad Prism) was used, comparing $\Delta lsr2_{pMV}$ or $\Delta lsr2_{T112A}$ to $\Delta lsr2_{WT}$ and $\Delta lsr2_{T112D}$.

### Acknowledgement

The following reagent was obtained through the NIH Biodefense and Emerging Infections Research Resources Repository, NIAID, NIH: Monoclonal Anti-*Mycobacterium tuberculosis* GroEL2 (Gene Rv0440), Clone IT70 (DCA4) (produced *in vitro*), NR-13657; Genomic DNA from *Mtb*, Strain H37Rv, NR-48669. We acknowledge the Centre for Core Biotechnology Services at the University of Leicester for support with containment level 3 experiments and analysis of mycobacterial proteins. The project was supported by the UK Biotechnology and Biological Sciences Research Council grants BB/K000330/1 and BB/P001513/1 (GVM); the French Infrastructure for Integrated Structural Biology (FRISBI) ANR-10-INBS-05 grant (MCG), the Ministry of Higher Education and Scientific Research Iraq (KA), and the Wellcome Trust 204538/Z/16/Z (SJW).

### Author contributions

Conceptualisation, OT, MCG, SJW and GVM; Methodology, ARB, PA, HJ, CR, AAW, MW; Investigation KA, OT, PB, HJ, ADV, GVM; Analysis, OT, AAW, KA, HJ, MCG, SJW, GVM; Resources, MW, ILB and MIV; Writing – Original Draft, GVM, MCG, SWJ; Writing – Review and Editing, GVM, SJW, MIV, HMO. Funding Acquisition, KA, MCG, SJW, GVM. Supervision, GVM, OT, HMO, MCG and SJW.

### Conflict of interest

Paul Ajuh is the director and shareholder in Gemini Biosciences Limited, Liverpool, UK. Other authors declare no competing interests.

### Data availability statement

Database: Mycob browser release 3 2018-06-05, 8 https://mycob browser.epfl.ch (Kapopoulou *et al.*, 2011). The accession numbers for the mass spectrometry proteomics data reported in this paper are ProteomeXchange Consortium via the PRIDE partner repository PXD009239 and 10.6019/PXD009239 (http://www.proteomexchange.org). The accession numbers for microarray data – ArrayExpress, E-MTAB-7627, https://www.ebi.ac.uk/arrayexpress; ChIP-Seq datasets – the European Nucleotide Archive (ENA), PRJEB31102, http://www.ebi.ac.uk/ena/data/view/PRJEB3110; for Lsr2 protein structures – the Protein Data Bank https://www.wwpdb.org, PDB6QKP and PDB6QKQ; the chemical shifts – the BMRB, BMRB ID 34358 and BMRB ID 34358.

### References

- Baer, C.E., Iavarone, A.T., Alber, T. and Sassetti, C.M. (2014) Biochemical and spatial coincidence in the provisional Ser/Thr protein kinase interaction network of *Mycobacterium tuberculosis*. *Journal of Biological Chemistry*, **289**, 20422–20433.

- Bannister, A.J. and Kouzarides, T. (2011) Regulation of chromatin by histone modifications. *Cell Research*, **21**, 381–395.

- Bartek, I.L., Woolhiser, L.K., Baughn, A.D., Basaraba, R.J., Jacobs, W.R. Jr., Lenaerts, A.J., *et al.* (2014) *Mycobacterium tuberculosis* Lsr2 is a global transcriptional regulator required for adaptation to changing oxygen levels and virulence. *mBio*, **5**, e01106–e01114.

- Barthe, P., Chiche, L., Declerck, N., Delsuc, M.A., Lefevre, J.F., Malliavin, T., *et al.* (1999) Refined solution structure and backbone dynamics of 15N-labeled C12A–p8MTCP1 studied by NMR relaxation. *Journal of Biomolecular NMR*, **15**, 271–288.

- Barthe, P., Mukamolova, G.V., Roumestand, C. and Cohen-Gonsaud, M. (2010) The structure of PknB extracellular PASTA domain from *Mycobacterium tuberculosis* suggests a ligand-dependent kinase activation. *Structure*, **18**, 606–615.

- Bellinzoni, M., Wehenkel, A.M., Duran, R. and Alzari, P.M. (2019) Novel mechanistic insights into physiological signaling pathways mediated by mycobacterial Ser/Thr protein kinases. *Genes and Immunity*, **20**, 383–393.

- Bosserman, R.E., Nguyen, T.T., Sanchez, K.G., Chirakos, A.E., Ferrell, M.J., Thompson, C.R., *et al.* (2017) WhiB6 regulation of ESX-1 gene expression is controlled by a negative feedback loop in *Mycobacterium marinum*. *Proceedings of the National Academy of Sciences*, **114**, 10772–10781.
Boutte, C.C., Baer, C.E., Papavinasasundaram, K., Liu, W., Chase, M.R., Meniche, X., et al. (2016). A cytoplasmic peptidoglycan amidase homologue controls mycobacterial cell wall synthesis. *eLife*, **5**, e14590. https://doi.org/10.7554/eLife.14590.

Camacho, L.R., Constant, P., Raynaud, C., Laneelle, M.A., Triccas, J.A., Gicquel, B., et al. (2001) Analysis of the phthiocerol dimycocerosate locus of *Mycobacterium tuberculosis*. Evidence that this lipid is involved in the cell wall permeability barrier. *Journal of Biological Chemistry*, **276**, 19845-19854.

Carette, X., Platig, J., Young, D.C., Helmel, M., Young, A.T., Wang, Z., et al. (2018) Multisystem analysis of *Mycobacterium tuberculosis* reveals kinase-dependent re-modeling of the pathogen-environment interface, *mBio*, **9**, e02333-17. https://doi.org/10.1128/mBio.02333-17.

Chawla, Y., Upadhyay, S., Khan, S., Nagarajan, S.N., Forti, F. and Nandicoori, V.K. (2014) Protein kinase B (PknB) of *Mycobacterium tuberculosis* is essential for growth of the pathogen in vitro as well as for survival within the host. *Journal of Biological Chemistry*, **289**, 13858-13875.

Chen, J.M., Ren, H., Shaw, J.E., Wang, Y.J., Li, M., Leung, A.S., et al. (2008) Lsr2 of *Mycobacterium tuberculosis* is a DNA-bridging protein. *Nucleic Acids Research*, **36**, 2123-2135.

Colangeli, R., Helb, D., Vilcheze, C., Hazbon, M.H., Lee, C.G., Safi, H., et al. (2007) Transcriptional regulation of multi-drug tolerance and antibiotic-induced responses by the histone-like protein Lsr2 in *M. tuberculosis*. *PLoS Path*, **3**, e87. https://doi.org/10.1371/journal.ppat.0030087.

Colangeli, R., Haq, A., Arcus, V.L., Summers, E., Magliozzo, R.S., McBride, A., et al. (2009) The multifunctional histone-like protein Lsr2 protects mycobacteria against reactive oxygen intermediates. *Proceedings of the National Academy of Sciences*, **106**, 4414-4418.

Cossu, A., Sechi, L.A., Bandino, E., Zanetti, S. and Rosu, V. (2013) Expression profiling of *Mycobacterium tuberculosis* H37Rv and *Mycobacterium smegmatis* in acid-nitrosative multi-stress displays defined regulatory networks. *Microbial Pathogenesis*, **65**, 89-96.

DeJesus, M.A., Gerrick, E.R., Xu, W., Park, S.W., Long, J.E., Boutte, C.C., et al. (2017) Comprehensive essentiality analysis of the *Mycobacterium tuberculosis* genome via saturating transposon mutagenesis, *mBio*, **8**, e02133-16. https://doi.org/10.1128/mBio.02133-16.

Dilweg, I.W. and Dame, R.T. (2018) Post-translational modification of nucleoid-associated proteins: an extra layer of functional modulation in bacteria? *Biochemical Society Transactions*, **46**, 1381-1392.

Fernandez, P., Saint-Joanis, B., Barilone, N., Jackson, M., Gicquel, B., Cole, S.T. et al. (2006) The Ser/Thr protein kinase PknB is essential for sustaining mycobacterial growth. *Journal of Bacteriology*, **188**, 7778-7784.

Forti, F., Crosta, A. and Ghisotti, D. (2009) Pristinamycinducible gene regulation in mycobacteria. *Journal of Biotechnology*, **140**, 270-277.

Galagan, J.E., Minch, K., Peterson, M., Lyubetskaya, A., Azizi, E., Sweet, L., et al. (2013) The *Mycobacterium tuberculosis* regulatory network and hypoxia. *Nature*, **499**, 178-183.

Garton, N.J., Waddell, S.J., Sherratt, A.L., Lee, S.M., Smith, R.J., Senner, C., et al. (2008) Cytological and transcript analyses reveal fat and lazy persister-like bacilli in tuberculous sputum. *PLoS Medicine*, **5**, e75. https://doi.org/10.1371/journal.pmed.0050075.

Gee, C.L., Papavinasasundaram, K.G., Blair, S.R., Baer, C.E., Falick, A.M., King, D.S., et al. (2012) A phosphorylated pseudokinase complex controls cell wall synthesis in mycobacteria. *Science Signaling*, **5**, ra7. https://doi.org/10.1126/scisignal.2002525.

Gordon, B.R., Imperial, R., Wang, L., Navarre, W.W. and Liu, J. (2008) Lsr2 of *Mycobacterium* represents a novel class of H-NS-like proteins. *Journal of Bacteriology*, **190**, 7052-7059.

Gordon, B.R., Li, Y., Wang, L., Sintsova, A., van Bakel, H., Tian, S., et al. (2010) Lsr2 is a nucleoid-associated protein that targets AT-rich sequences and virulence genes in *Mycobacterium tuberculosis*. *Proceedings of the National Academy of Sciences*, **107**, 5154-5159.

Gordon, B.R., Li, Y., Cote, A., Weirauch, M.T., Ding, P., Hughes, T.R., et al. (2011) Structural basis for recognition of AT-rich DNA by unrelated xenogeneic silencing proteins. *Proceedings of the National Academy of Sciences*, **108**, 10690-10695.

Güntert, P. (2004) Automated NMR structure calculation with CYANA. *Methods in Molecular Biology*, **278**, 353-378.

Iswahyudi, G.V.M., Straatman-Iwanowska, A.A., Allcock, N., Ajuh, P., Turapov, O., and O'Hare, H.M. (2019) Mycobacterial phosphatase PstP regulates global serine threonine phosphorylation and cell division. *Scientific Reports*, **9**, 8337. https://doi.org/10.1038/s41598-019-44841-9.

Jarvis, B., Wilrich, C. and Wilrich, P.T. (2010) Reconsideration of the derivation of most probable numbers, their standard deviations, confidence bounds and rarity values. *Journal of Applied Microbiology*, **109**, 1660-1667.

Kang, C.M., Abbott, D.W., Park, S.T., Dascher, C.C., Cantley, L.C. and Husson, R.N. (2005) The *Mycobacterium tuberculosis* serine/threonine kinases PknA and PknB: substrate identification and regulation of cell shape. *Genes and Development*, **19**, 1692-1704.

Kapopoulou, A., Lew, J.M. and Cole, S.T. (2011) The MycoBrowser portal: a comprehensive and manually annotated resource for mycobacterial genomes. *Tuberculosis* (*Edinburgh*), **91**, 8-13.

Kriel, N.L., Gallant, J., van Wyk, N., van Helden, P., Sampson, S.L., Warren, R.M., et al. (2018) Mycobacterial nucleoid associated proteins: an added dimension in gene regulation. *Tuberculosis* (*Edinburgh*), **108**, 169-177.

Lee, J.H., Karakousis, P.C. and Bishai, W.R. (2008) Roles of SigB and SigF in the *Mycobacterium tuberculosis* sigma factor network. *Journal of Bacteriology*, **190**, 699-707.

Leiba, J., Carrere-Kremer, S., Blondiaux, N., Dimala, M.M., Wohlkonig, A., Baulard, A., et al. (2014) The *Mycobacterium tuberculosis* transcriptional repressor EthR is negatively regulated by Serine/Threonine phosphorylation. *Biochemical and Biophysical Research Communications*, **446**, 1132-1138.

Li, H. and Durbin, R. (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics*, **25**, 1754-1760.
Loraine, J., Pu, F., Turapov, O. and Mukamolova, G.V. (2016) Development of an *in vitro* assay for detection of drug-induced Resuscitation-promoting-factor-dependent mycobacteria. *Antimicrobial Agents and Chemotherapy, **60**, 6227–6233.

Lucchini, S., Rowley, G., Goldberg, M.D., Hurd, D., Harrison, M. and Hinton, J.C. (2006) H-NS mediates the silencing of laterally acquired genes in bacteria. *PLoS Path*, **2**, e81. https://doi.org/10.1371/journal.ppat.0030038.

Marcus, S.A., Sidiropoulos, S.W., Steinberg, H. and Talaat, A.M. (2016) CsoR is essential for maintaining copper homeostasis in *Mycobacterium tuberculosis*. *PLoS ONE*, **11**, e0151816.

Masiewicz, P., Brzostek, A., Wolanski, M., Dziadek, J. and Zakrzewska-Czerwinska, J. (2012) A novel role of the PrpR as a transcription factor involved in the regulation of methylcitrate pathway in *Mycobacterium tuberculosis*. *PLoS ONE*, **7**, e43651.

Mehta, M. and Singh, A. (2018) *Mycobacterium tuberculosis* WhiB3 maintains redox homeostasis and survival in response to reactive oxygen and nitrogen species. *Free Radical Biology and Medicine*, **131**, 50–58.

Minch, K.J., Rustad, T.R., Peterson, E.J., Winkler, J., Reiss, D.J., Ma, S., et al. (2015) The DNA-binding network of *Mycobacterium tuberculosis*. *Nature Communications*, **6**, 5829. https://doi.org/10.1038/ncomms6829.

Ortega, C., Liao, R., Anderson, L.N., Rustad, T., Ollodart, A.R., Wright, A.T., et al. (2014) *Mycobacterium tuberculosis* Ser/Thr protein kinase B mediates an oxygen-dependent replication switch. *PLoS Biology*, **12**, e1001746.

Park, H.D., Guinn, K.M., Harrell, M.I., Liao, R., Voskuil, M.I., Tompa, M., et al. (2003) Rv3133c/dosR is a transcription factor that mediates the hypoxic response of *Mycobacterium tuberculosis*. *Molecular Microbiology*, **48**, 833–843.

Park, S.T., Kang, C.M. and Husson, R.N. (2008) Regulation of the SigH stress response regulon by an essential protein kinase in *Mycobacterium tuberculosis*. *Proceedings of the National Academy of Sciences*, **105**, 13105–13110.

Pons, J.L., Malliavin, T.E. and Delsuc, M.A. (1996) Gifa V. 4: A complete package for NMR data set processing. *Journal of Biomolecular NMR*, **8**, 445–452.

Prigozhin, D.M., Papavinasasundaram, K.G., Baer, C.E., Murphy, K.C., Moskaleva, A., Chen, T.Y., et al. (2016) Structural and genetic analyses of the *Mycobacterium tuberculosis* protein kinase B sensor domain identify a potential ligand-binding site. *Journal of Biological Chemistry*, **291**, 22961–22969.

Prisic, S., Dankwa, S., Schwartz, D., Chou, M.F., Locasale, J.W., Kang, C.M., et al. (2010) Extensive phosphorylation with overlapping specificity by *Mycobacterium tuberculosis* serine/threonine protein kinases. *Proceedings of the National Academy of Sciences*, **107**, 7521–7526.

Prisic, S., Hwang, H., Dow, A., Barnaby, O., Pan, T.S., Lonzanida, J.A., et al. (2015) Zinc regulates a switch between primary and alternative S18 ribosomal proteins in *Mycobacterium tuberculosis*. *Molecular Microbiology*, **97**, 263–280.

Rachman, H., Lee, J.S., Angermann, J., Kowall, J. and Kaufmann, S.H. (2006) Reliable amplification method for bacterial RNA. *Journal of Biotechnology*, **126**, 61–68.

Reddy, T.B., Riley, R., Wymore, F., Montgomery, P., DeCaprio, D., Engels, R., et al. (2009) TB database: an integrated platform for tuberculosis research. *Nucleic Acids Research*, **37**, D499–D508.

Richard-Greenblatt, M. and Av-Gay, Y. (2017) Epigenetic phosphorylation control of *Mycobacterium tuberculosis* infection and persistence. *Microbiology Spectrum*, **5**, https://doi.org/10.1128/microbiolspec.TBTB2-0005-2015.

Robinson, J.T., Thorvaldsdottir, H., Winckler, W., Guttman, M., Lander, E.S., Getz, G., et al. (2011) Integrative genomics viewer. *Nature Biotechnology*, **29**, 24–26.

Rohde, K.H., Abramovitch, R.B. and Russell, D.G. (2007) *Mycobacterium tuberculosis* invasion of macrophages: linking bacterial gene expression to environmental cues. *Cell Host and Microbe*, **2**, 352–364.

Rustad, T.R., Harrell, M.I., Liao, R. and Sherman, D.R. (2008) The enduring hypoxic response of *Mycobacterium tuberculosis*. *PLoS ONE*, **3**, e1502.

Rustad, T.R., Minch, K.J., Ma, S., Winkler, J.K., Hobbs, S., Hickey, M., et al. (2014) Mapping and manipulating the *Mycobacterium tuberculosis* transcriptome using a transcription factor overexpression-derived regulatory network. *Genome Biology*, **15**, 502. http://genomebiology.com/2014/15/11/502.

Salina, E.G., Waddell, S.J., Hoffmann, N., Rosenkrands, I., Butcher, P.D. and Kaprelyants, A.S. (2014) Potassium availability triggers *Mycobacterium tuberculosis* transition to, and resuscitation from, non-culturable (dormant) states. *Open Biology*, **4**, 140106. http://dx.doi.org/10.1098/rsob.140106.

Sassetti, C.M. and Rubin, E.J. (2003) Genetic requirements for mycobacterial survival during infection. *Proceedings of the National Academy of Sciences*, **100**, 12989–12994.

Schnappinger, D., Ehrt, S., Voskuil, M.I., Liu, Y., Mangan, J.A., Monahan, I.M., et al. (2003) Transcriptional adaptation of *Mycobacterium tuberculosis* within macrophages: insights into the phagosomal environment. *Journal of Experimental Medicine*, **198**, 693–704.

Shen, Y., Delaglio, F., Cornilescu, G. and Bax, A. (2009) TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts. *Journal of Biomolecular NMR*, **44**, 213–223.

Smith, L.J., Bochkareva, A., Rolfe, M.D., Hunt, D.M., Kahramanoglou, C., Braun, Y., et al. (2017) Cmr is a re-dox-responsive regulator of DosR that contributes to *M. tuberculosis* virulence. *Nucleic Acids Research*, **45**, 6600–6612.

Tailleux, L., Waddell, S.J., Pelizzola, M., Mortellaro, A., Withers, M., Tanne, A., et al. (2008) Probing host pathogen cross-talk by transcriptional profiling of both *Mycobacterium tuberculosis* and infected human dendritic cells and macrophages. *PLoS ONE*, **3**, e1403.

Turapov, O., Loraine, J., Jenkins, C., Barthe, P., McFeely, D., Forti, F., et al. (2015) The external PASTA domain of the essential serine/threonine protein kinase PknB regulates mycobacterial growth. *Open Biology*, **5**, 150025.

Turapov, O., Forti, F., Kadhim, B., Ghisotti, D., Sassine, J., Straatman-Iwanowska, A., et al. (2018) Two faces of CwIM, an essential PknB substrate, in *Mycobacterium tuberculosis*. *Cell Reports*, **25**, 57–67.
Tzeng, S.R. and Kalodimos, C.G. (2012) Protein activity regulation by conformational entropy. *Nature*, **488**, 236–240.

Waddell, S.J. and Butcher, P.D. (2010) Use of DNA arrays to study transcriptional responses to antimycobacterial compounds. I. In Antibiotic resistance protocols (S.H. Gillespie and T.D. McHugh Eds.). *Methods Molecular Biology*, **642**, 75–91.

Wayne, L.G. and Sohaskey, C.D. (2001) Nonreplicating persistence of *Mycobacterium tuberculosis*. *Annual Review of Microbiology*, **55**, 139–163.

Wayne, L.G. and Sramek, H.A. (1994) Metronidazole is bactericidal to dormant cells of *Mycobacterium tuberculosis*. *Antimicrobial Agents and Chemotherapy*, **38**, 2054–2058.

Willemse, D., Weber, B., Masino, L., Warren, R.M., Adinolfi, S., Pastore, A., et al. (2018) Rv1460, a SufR homologue, is

a repressor of the Suf operon in *Mycobacterium tuberculosis*. *PLoS ONE*, **13**, e0200145.

Winardhi, R.S., Yan, J. and Kenney, L.J. (2015) H-NS regulates gene expression and compacts the nucleoid: insights from single-molecule experiments. *Biophysical Journal*, **109**, 1321–1329.

Zhang, Y., Liu, T., Meyer, C.A., Eeckhout, J., Johnson, D.S., Bernstein, B.E., et al. (2008) Model-based analysis of ChIP-Seq (MACS). *Genome Biology*, **9**, R137. https://doi.org/10.1186/gb-2008-9-9-r137.

## Supporting Information

Additional supporting information may be found online in the Supporting Information section at the end of the article