INTRODUCTION
Tuberculosis (TB), an infectious disease caused by the bacterium
Mycobacterium tuberculosis, is one of the leading causes of death from infectious diseases (
1). The fact that TB treatment requires at least a six-month regimen with four antibiotics is partly due to the intrinsic antibiotic tolerance of
M. tuberculosis (
2,
3). Stressed
M. tuberculosis cells can achieve a dormant or slow-growing state (
4,
5) that exhibits antibiotic tolerance (
6), cell wall thickening (
7), and altered cell wall staining (
4).
The currently accepted cell wall structure of
M. tuberculosis (
8) is composed of three covalently linked layers (
9) as follows: surrounding the plasma membrane, a (i) peptidoglycan (PG) layer is covalently bound to an (ii) arabinogalactan layer; a (iii) lipid layer composed of mycolic acids surrounds the arabinogalactan layer, and the inner leaflet of this layer is covalently linked to the arabinogalactan (
10); and the outer leaflet of the mycolic acid layer contains free mycolic acids, trehalose mycolates and other lipids, glycolipids, glycans and proteins (
11). The mycolic acid layer, or mycomembrane, is the outer membrane of mycobacteria and is the major contributor to impermeability of the cell wall (
12–14).
In addition to serving as a permeability barrier, the regulation of the cell wall likely contributes to antibiotic tolerance either through further changes in permeability (
15) or by changing the activity of antibiotic targets (
16). Several studies have observed changes in the cell wall under stress (
7,
15,
17,
18). These cell wall changes have been shown to correlate with increased antibiotic tolerance (
19–21). This has led the prevalent model that stress-induced regulation of the cell wall contributes to antibiotic tolerance (
22). While most of the extant data to support this model is correlative, we recently identified a mutant in
M. smegmatis that specifically upregulates peptidoglycan metabolism in starvation and also causes decreased antibiotic tolerance in that condition (
23). This shows there is a causal relationship between cell wall regulation and antibiotic tolerance, at least in limited conditions in
M. smegmatis.
Reversible protein phosphorylation is a key regulatory tool used by bacteria for environmental signal transduction to regulate cell growth (
24–27). In
M. tuberculosis, serine/threonine phosphorylation is important in cell wall regulation (
28).
M. tuberculosis has 11 serine/threonine protein kinases (STPKs) (PknA, PknB, and PknD to PknL) and only one serine/threonine protein phosphatase (PstP) (
29,
30).
Among the STPKs, PknA and PknB are essential for growth and phosphorylate many substrates involved in cell growth and division (
23,
31–34). Some of these substrates are enzymes whose activity is directly altered by phosphorylation. For example, all the enzymes in the fatty acid synthase II (FAS-II) system of mycolic acid biosynthesis are inhibited by threonine phosphorylation (
35–40). There are also cell wall regulators that are not enzymes, but whose phosphorylation by STPKs affects cell shape and growth. For example, the regulator CwlM, once it is phosphorylated by PknB, activates MurA (
23), the first enzyme in PG precursor biosynthesis (
41). In the transition to starvation, CwlM is rapidly dephosphorylated in
M. smegmatis (
23). Misregulation of MurA activity increases sensitivity to antibiotics in early starvation (
23), implying that phosphoregulation of CwlM promotes antibiotic tolerance. CwlM may also regulate other steps of PG synthesis (
42). A recent phosphoproteomic study showed that transcriptional repression of the operon that contains both
pstP and
pknB leads to increased phosphorylation of CwlM (
43). While the effects of the individual genes were not separated (
43), this suggests that PstP could dephosphorylate CwlM.
PstP is essential in
M. tuberculosis and
M. smegmatis (
44,
45). It is a member of the protein phosphatase 2C (PP2C) subfamily of metal-dependent protein serine/threonine phosphatases (
46) that strictly require divalent metal ions for activity (
47,
48). PP2C phosphatases are involved in responding to environmental signals, regulating metabolic processes, sporulation, cell growth, division, and stress response in a diverse range of prokaryotes and eukaryotes (
49–54). PstP
Mtb shares structural folds and conserved residues with the human PP2Cα (
55), which serves as the representative of the PP2C family. PstP
Mtb has an N-terminal cytoplasmic enzymatic domain, a transmembrane pass, and a C-terminal extracellular domain (
46).
Many of the proteins known to be dephosphorylated by PstP (
35,
46,
56–59) are involved in cell wall metabolism; however, the effects of this activity seem to differ. For example, dephosphorylation of CwlM should decrease PG metabolism in stasis (
23), but dephosphorylation of the FAS-II enzymes (
35–40) should upregulate lipid metabolism in growth. However, PG and lipid metabolism are expected to be coordinated (
22). Therefore, PstP must be able to alter substrate specificity in growth and in stasis.
PstP
Mtb is itself phosphorylated on threonine residues 137, 141, 174, and 290 (
56). We hypothesized that phosphorylation of the threonine residues of PstP might help coordinate activity against different substrates through changes in access to substrates, or through toggling catalytic activity against substrates.
We report here that phosphoablative and phosphomimetic mutations at the phosphosite T171 of PstPMsmeg (T174 in PstPMtb) alter growth rate, cell length, cell wall metabolism, and antibiotic tolerance in M. smegmatis. Strains of M. smegmatis with pstP T171E alleles grow slowly, are unable to properly downregulate PG metabolism, and upregulate antibiotic tolerance in the transition to starvation. We observed that the same mutation has nearly opposite effects on mycolic acid layer metabolism. We also report that PstPMtb dephosphorylates CwlMMtb.
DISCUSSION
Previous studies on mycobacterial phosphoregulation suggest that PstP could play a critical role in modulating cell wall metabolism in the transition between growth and stasis (
18,
22,
23,
35,
43,
56,
59,
79). In this work, we explored how the phosphorylation of PstP contributes to this regulation. We report here that the phosphosite T171 of PstP
Msmeg impacts growth, cell wall metabolism, and antibiotic tolerance. We found that the PG master regulator CwlM
Mtb is a substrate of PstP
Mtb. Our findings indicate that the phosphorylation on PstP affects PG metabolism in stasis and mycolic acid metabolism during growth.
PG is regulated by phosphorylation factors at several points along the biosynthesis pathway (
23,
42,
80,
81), mostly by PknB. PknB’s kinase activity is responsive to lipid II that it detects in the periplasm (
82). PstP is a global negative regulator of STPK phosphorylation (
43) and has been proposed to be the cognate phosphatase of PknB in regulating cell growth (
22,
43,
59,
83). Our data suggest that mutations at T171 of PstP do not affect PG metabolism in growth (
Fig. 4A; see Fig. S2A and S3A in the supplemental material) but that the PstP
Msmeg T171E strain fails to downregulate PG in starvation (
Fig. 4C and
E; Fig. S3C). We expect that the activity of PstP against the PG regulator CwlM (
Fig. 2A, top panel) should be critical for this downregulation because it should deactivate MurA, the first enzyme in PG precursor synthesis (
23,
41).
The
in vitro biochemistry (
Fig. 2A and
B) predicts the log-phase staining data (
Fig. 4A; Fig. S2A and S3A), where the
pstP T171E variant shows no difference in apparent PG activity. The proximity of a phosphosite to the substrate-binding site of an enzyme may affect the catalytic activity directly (
84), but T174 maps to the β-sheet core (β8) in PstP
Mtb, which is distant from the active site (
Fig. 1A) (
55). The PG staining in starvation (
Fig. 4C and
E; Fig. S3C) suggests that the PstP
Msmeg T171E phosphomimetic variant might dephosphorylate CwlM more slowly
in vivo (
Fig. 5B and
D), but this is not what we see
in vitro (
Fig. 2A and
B). Therefore, it is possible that, in starvation, phosphorylation at this site affects interaction (
85) with other regulatory proteins (
86–88) that could modulate PstP’s activity against PG substrates, or it could affect access to substrates via localization changes.
Synthesis of the various mycobacterial cell wall layers are likely synchronized (
22,
68). PknB almost surely plays a crucial role in connecting PG and mycolic acid metabolism during growth. If PG metabolism is slowed, PknB could sense the accumulation of periplasmic lipid II (
82) and signal to halt mycolic acid biosynthesis by inactivating the FAS-II enzymes and the trehalose monomycolate transporter MmpL3 (
89) via phosphorylation (
35–40,
83). Our data imply that PstP helps balance the effects of these inhibitory phosphorylations to allow coordinated synthesis of mycolic acids in log phase (
Fig. 4B and
5E; Fig. S2B and S3B). Misphosphorylation of PstP likely disrupts this coordination and seems to decrease mycolic acid layer metabolism. This may partly explain the slow growth of the
pstP T171E mutants (
Fig. 1B).
DMN-Tre incorporation into the mycomembrane is directly catalyzed by secreted (
90) Ag85 enzymes (
66). DMN-Tre fluorescence depends on mycomembrane hydrophobicity and is affected by inhibition of cytoplasmic mycolic acid synthesis (
66). Hydrophobicity can be affected by the glycolipid composition of the mycomembrane (
91). The differences in DMN-Tre fluorescence that we see (
Fig. 4B; Fig. S2B and S3B) could be due to changes in mycolic acid synthesis or due to changes in other glycolipids that affect hydrophobicity (
91). Our DMN-Tre data clearly indicate that the trehalose mycolate leaflet of the mycomembrane is affected due to phosphomisregulation of PstP, although it does not yield detailed information about how it is affected.
We propose that PstP’s regulation of mycolic acid layer biosynthesis occurs in the cytoplasm. PstP and all the STPKs work in the cytoplasm, and there are currently no known systems whereby secreted proteins like Ag85 can be regulated by phosphorylation. All the enzymes of the FAS-II complex, which elongates fatty acids into the long lipids used in mycolic acids (
75), are downregulated by phosphorylation (
35–40), and two are biochemically verified substrates of PstP (
35). MmpL3, the mycolic acid flippase (
89), is also inhibited by phosphorylation (
83). It is likely that PstP could affect the activity of the entire FAS-II complex, including the target of isoniazid, InhA, which is inactivated by threonine phosphorylation (
36,
38). Isoniazid is a small hydrophilic drug and undergoes active diffusion via the porins (
78,
92); therefore, alterations in mycomembrane permeability are not likely to contribute substantially to differences in isoniazid sensitivity. Although our data (
Fig. 5E) do not reveal the exact misregulated spot in the mycolic acid synthesis and transport pathway, the higher susceptibility of the phosphomimetic strain (
Fig. 5E) to isoniazid suggests that this metabolic pathway is affected. Our DMN-Tre staining also suggests there should be a balance of the non-phospho- and phosphoforms of PstP
Msmeg T171 (
Fig. 4B) during growth to regulate mycomembrane biosynthesis.
PstP may dephosphorylate the cell wall substrates directly and/or by deactivating their kinases (
93) in both the PG and mycolic acid biosynthesis pathways. All these data combined suggest a complex cross-talk of the STPKs and PstP to regulate diverse cell wall substrates.
MATERIALS AND METHODS
Bacterial strains and culture conditions.
All
Mycobacterium smegmatis mc
2155 ATCC 700084 cultures were started in 7H9 (Becton, Dickinson, Franklin Lakes, NJ) medium containing 5 g/liter bovine serum albumin (BSA), 2 g/liter dextrose, 0.003 g/liter catalase, 0.85 g/liter NaCl, 0.2% glycerol, and 0.05% Tween 80 and incubated at 37°C until log phase. Hartmans-de Bont (HdB) minimal medium made as described previously (
94) without glycerol was used for starvation assays. Serial dilutions of all CFU counts were plated on LB Lennox agar (Fisher BP1427-2).
E. coli Top10, XL1-Blue, and Dh5α strains were used for cloning and E. coli strain BL21 Codon Plus was used for protein expression. Antibiotic concentrations for M. smegmatis were 25 μg/ml kanamycin, 50 μg/ml hygromycin, and 20 μg/ml zeocin. Antibiotic concentrations for E. coli were 50 μg/ml kanamycin, 25 μg/ml zeocin, 20 μg/ml chloramphenicol, and 140 μg/ml ampicillin.
Strain construction.
The PstP
Msmeg knockdown strain was made first by creating a merodiploid strain and then by deleting the native
pstPMsmeg gene from its chromosomal location. The merodiploid strain was generated by introducing a constitutively expressing
pstPMtb gene cloned on an StrR plasmid at the L5 attB integration site. The
pstPMsmeg gene (MSMEG_0033) at the native locus was then deleted by RecET-mediated double-stranded recombineering approach using a 1.53-kb loxP-hyg-loxP fragment carrying a 125-bp region flanking the
pstPMsmeg gene, as described previously (
95). The recombineering substrate was generated by two sequential overlapping PCRs of the loxP-hyg-loxP substrate present in the plasmid pKM342. The downstream flanking primer used in the first PCR also carried an optimized mycobacterial ribosome binding site in front of the start codon of MSMEG_0032 to facilitate the expression of the genes present downstream of
pstPMsmeg in the
M. smegmatis pstP-pknB operon.
Deletion of the
pstPMsmeg gene was confirmed by PCR amplification and sequencing of the 5′ and 3′ recombinant junctions, and the absence of an internal wild-type
pstPMsmeg PCR product. The
pstPMtb allele present at the L5 site was then swapped, as described previously (
96), with a Tet-regulatable
pstPMtb allele (RevTetR-P750-
pstPMtb-DAS tag-L5-Zeo plasmid). The loxP-flanked
hyg marker present in the chromosomal locus was then removed by expressing Cre from pCre-sacB-Kan, and the Cre plasmid was subsequently cured from this strain by plating on sucrose. We named this strain CB1175.
Different alleles of
pstP were attained by swapping the wild-type (WT) allele at the L5 site of CB1175 as described previously (
97). In order to do so, WT and the phosphoablative alleles of
pstPMsmeg were first cloned individually into a kanamycin resistance-marked L5 vector pCT94 under the control of a TetO promoter to generate vectors pCB1206-1208 and pCB1210, which would swap out the zeocin resistance-marked vector at the L5 site in CB1175. The strong TetO promoter in the vectors pCB1206-1208 and pCB1210 was swapped with an intermediate-strength promoter p766TetON6 (cloned from the vector pCB1030 [pGMCgS-TetON-6 sspB]) to generate the L5 vectors pCB1282 to 85. pCB1285 was used as the parent vector later on to clone in the phosphomimetic
pstPMsmeg alleles under the control of p766TetON6.
These kanamycin resistance-marked vector constructs were then used to swap out the zeocin resistance-marked vector at the L5 site of CB1175 to attain different allelic strains of
pstPMsmeg as described previously (
97).
Growth curve assay.
At least three biological replicates of different pstPMsmeg allele variants (T171A, T171E, and WT) were grown in 7H9 medium up to log phase. The growth curves were performed in nontreated 96-well plates using a plate reader (BioTek Synergy neo2 multi mode reader) in 200 μl of 7H9 medium starting at an optical density at 600 nm (OD600) of 0.1. An exponential growth equation was used to calculate the doubling times of each strain using the least squared ordinary fit method in GraphPad Prism (version 7.0d). P values were calculated using two-tailed, unpaired t tests.
Cell staining.
Three biological replicate strains of each pstP allelic variant (T171A, T171E, and WT) were used for this assay. For staining cells in log phase, 100 μl of culture in 7H9 medium was incubated at 37°C with 1 μl of 10 mM DMN-Tre for 30 min and 1 μl of 10 mM HADA for 15 min. Cells were then pelleted and resuspended in 1× phosphate-buffered saline (PBS) supplemented with 0.05% Tween 80 and fixed with 10 μl of 16% paraformaldehyde (PFA) for 10 min at room temperature. Cells were then washed and resuspended in PBS plus Tween 80.
For starvation microscopy, cultures were shaken for 4 h in HdB medium without glycerol at 37°C. Aliquots (500 μl) of each culture were pelleted and concentrated to 100 μl, then incubated at 37°C with 1 μl of 10 mM DMN-Tre for 1 h and 3 μl of 10 mM HADA for 30 min. Cells were then washed and fixed as above. The total time of starvation before fixation was 5.5 h.
Microscopy and image analysis.
Cells were imaged with a Nikon Ti-2 widefield epifluorescence microscope with a Photometrics Prime 95B camera and a Plan Apo 100×, 1.45 -numerical-aperture (NA) lens objective. The green fluorescence images for DMN-Tre staining were taken with a 470/40 nm excitation filter and a 525/50 nm emission filter. Blue fluorescence images for HADA staining were taken using a 350/50 nm excitation filter and a 460/50 nm emission filter. All images were captured using NIS Elements software and analyzed using FIJI and MicrobeJ (
98). For cell detection in MicrobeJ, appropriate parameters for length, width, and area were set. The V-snapping cells were split at the septum so that each daughter cell could be considered a single cell. Any overlapping cells were excluded from analysis.
Length and mean intensities of HADA and DMN-Tre signals of 300 cells from each of pstPMsmeg T171A, pstPMsmeg T171E, and pstPMsmeg WT (100 cells from each of three biological replicate strains of each genotype) were quantified using MicrobeJ. The values of the mean intensities of 300 cells of each pstP allelic mutant and WT are represented in the graph as percentages of the highest mean intensity from all the cells in that experiment. GraphPad Prism (v7.0d) was used to generate the graphs and perform t tests. To compare two groups (pstPMsmeg T171A versus pstPMsmeg WT and pstPMsmeg T171E versus pstPMsmeg WT), the lengths and the percentage-intensity values of 300 cells per pstP allelic variant genotype (100 per biological replicate strain per genotype) were used as data array inputs to perform unpaired, two-tailed t tests in GraphPad Prism.
Medial intensity profiles of DMN-Tre and NADA signals in cells from different pstP allelic strains in log phase and starvation analyzed with MicrobeJ were plotted on the y axis over relative positions of cells using the “XStatProfile” plotting feature in MicrobeJ to show subcellular localization of fluorescent intensities.
Demographs of DMN-Tre and NADA signal intensity across cell lengths in log phase and starvation were built using the “Demograph” feature of MicrobeJ by plotting the medial intensity profiles of DMN-Tre and NADA signals.
Western blotting.
Cultures were grown in 7H9 medium to OD600 = 0.8 in 10 ml of 7H9 medium, pelleted and resuspended in 500 μl PBS with 1 mM phenylmethylsulfonyl fluoride (PMSF), and lysed (MiniBeadBeater-16, model 607, Biospec). Supernatants from the cell lysates were run on 12% resolving Tris-glycine gels and then transferred onto polyvinylidene difluoride (PVDF) membranes (GE Healthcare). Rabbit anti-Strep-tag II antibody (1:1,000, Abcam, ab76949) in Tris-buffered saline with Tween 20 (TBST) buffer with 0.5% milk and goat anti-rabbit IgG (H+L) horseradish peroxidase (HRP)-conjugated secondary antibody (1:10,000, Thermo Fisher Scientific 31460) in TBST were used to detect PstP–Strep-tag II. For starvation experiments, cultures were first grown to log phase, then starved in HdB no-glycerol medium starting at OD600 = 0.5 for 1.5 h.
For Western blotting of in vitro assays, samples were run on 12% SDS gel (Mini-Protean TGX, Bio-Rad, 4561046) and then transferred onto a PVDF membrane (GE Healthcare). Mouse anti-His antibody (1:1,000, Genscript A00186) in TBST buffer with 0.5% BSA and goat anti-mouse IgG (H+L) HRP-conjugated secondary antibody (1:10,000, Invitrogen A28177) were used to detect His-tagged proteins on the blot. The blots were stripped (Thermo Fisher Scientific, 21059) and reprobed with rabbit antiphosphothreonine antibody (1:1,000, Cell Signaling number 9381) and goat anti-rabbit IgG (H+L) HRP-conjugated secondary antibody (1:10,000, Thermo Fisher Scientific 31460) to detect phosphorylation on the blots.
Antibiotic assays.
Biological triplicates of each pstP allelic variant were used for all antibiotic assays. For antibiotic assays in log phase, log phase cultures were diluted in 7H9 medium to OD600 = 0.05 before treatment. For starvation assays, cells were grown to OD600 = 0.5, pelleted, washed, and resuspended in HdB starvation medium (with no glycerol and 0.05% Tween) at OD600 = 0.3 and incubated at 37°C for a total of 5.5 h. The cultures were then diluted to OD600 = 0.05 before antibiotic treatment. Meropenem (8 μg/ml and 45 μg/ml), isoniazid (10 μg/ml and 90 μg/ml), d-cycloserine (100 μg/ml and 900 μg/ml), and trimethoprim (50 μg/ml and 360 μg/ml) were used separately to treat log-phase and starved cultures, respectively. Samples from the cultures were serially diluted and plated on LB agar before and after treatment, and CFU were calculated.
Protein purification.
All the proteins were expressed using E. coli BL21 Codon Plus cells.
N-terminally His-MBP-tagged PknB
Mtb was expressed and purified as described previously (
80). His-PstP
cWT
Mtb (1 to 300 amino acids of the cytosolic domain [
99]) and His-SUMO-CwlM
Mtb were both expressed overnight by IPTG (isopropyl-β-
d-thiogalactopyranoside) induction (1 mM and 1.3 mM, respectively), purified on Ni-nitrilotriacetic acid (Ni-NTA) resin (G-Biosciences, number 786-940 in 5 ml Bio-Scale Mini Cartridges, Bio-Rad number 7324661), then dialyzed, concentrated, and run over size exclusion resin (GE Biosciences Sephacryl S-200 in HiPrep 26/70 column) to obtain soluble proteins. The buffer for His-SUMO-CwlM
Mtb was 50 mM Tris pH 8, 350 mM NaCl, 1 mM dithiothreitol (DTT), and 10% glycerol. The buffer for His-PstP
cWT
Mtb was 50 mM Tris pH 7.5, 350 mM NaCl, 1 mM DTT, and 10% glycerol. Imidazole (20 mM) was added to each buffer for lysis and application to the Ni-NTA column, and 250 mM imidazole was added for elution. His-PstP
cT174E
Mtb was expressed and purified using the same conditions and buffers used for His-PstP
cWT
Mtb.
In vitro dephosphorylation assay.
Purified His-SUMO-CwlMMtb was phosphorylated with the purified kinase His-MBP-PknBMtb for 1 h at room temperature in the presence of 0.5 mM ATP, 1 mM MnCl2, and buffer (50 mM Tris [pH 7.5], 250 mM NaCl, and 1 mM DTT). The amount of kinase was one-tenth of the amount of substrate in the phosphorylation reaction. To stop the kinase reaction by depleting ATP, 0.434 unit of calf intestinal alkaline phosphatase (Quick CIP; New England BioLabs, MO525S) per μg of His-SUMO-CwlMMtb was added to the reaction mixture and incubated for 1 h at 37°C. The reaction mixture was then divided into five parts for the different phosphatase samples and a control with buffer.
Two individually expressed and purified batches of both His-PstPcWTMtb and His-PstPcT1714EMtb were used as biological replicates to perform the dephosphorylation assay. The reaction was carried out at room temperature for up to 90 min in the presence of phosphatase buffer (50 mM Tris [pH 7.5], 10 mM MnCl2, and 1 mM DTT). The amount of phosphatase used was half the amount of His-SUMO-CwlMMtb.
The intensities of the anti-His and the antiphosphothreonine signals on the blots were quantified with FIJI. The intensities of the anti-His and the antiphosphothreonine signals at each time point were normalized against the respective antibodysignal intensity at 0 min. These relative intensities were used to calculate antiphosphothreonine/anti-His for each time point and the values were plotted over time using GraphPad Prism (version 7.0d).