key: cord-0003236-n057xgzc authors: Grover, Sonam; Sharma, Tarina; Singh, Yadvir; Kohli, Sakshi; P., Manjunath; Singh, Aditi; Semmler, Torsten; Wieler, Lothar H.; Tedin, Karsten; Ehtesham, Nasreen Z.; Hasnain, Seyed E. title: The PGRS Domain of Mycobacterium tuberculosis PE_PGRS Protein Rv0297 Is Involved in Endoplasmic Reticulum Stress-Mediated Apoptosis through Toll-Like Receptor 4 date: 2018-06-19 journal: mBio DOI: 10.1128/mbio.01017-18 sha: d921cd2d71ffa5577a0d2c534da20f6780b3a564 doc_id: 3236 cord_uid: n057xgzc The genome of Mycobacterium tuberculosis, the causal organism of tuberculosis (TB), encodes a unique protein family known as the PE/PPE/PGRS family, present exclusively in the genus Mycobacterium and nowhere else in the living kingdom, with largely unexplored functions. We describe the functional significance of the PGRS domain of Rv0297, a member of this family. In silico analyses revealed the presence of intrinsically disordered stretches and putative endoplasmic reticulum (ER) localization signals in the PGRS domain of Rv0297 (Rv0297PGRS). The PGRS domain aids in ER localization, which was shown by infecting macrophage cells with M. tuberculosis and by overexpressing the protein by transfection in macrophage cells followed by activation of the unfolded protein response, as evident from increased expression of GRP78/GRP94 and CHOP/ATF4, leading to disruption of intracellular Ca(2+) homeostasis and increased nitric oxide (NO) and reactive oxygen species (ROS) production. The consequent activation of the effector caspase-8 resulted in apoptosis of macrophages, which was Toll-like receptor 4 (TLR4) dependent. Administration of recombinant Rv0297PGRS (rRv0297PGRS) also exhibited similar effects. These results implicate a hitherto-unknown role of the PGRS domain of the PE_PGRS protein family in ER stress-mediated cell death through TLR4. Since this protein is already known to be present at later stages of infection in human granulomas it points to the possibility of it being employed by M. tuberculosis for its dissemination via an apoptotic mechanism. T uberculosis (TB) remains a major public health problem and is caused by infections with the pathogen Mycobacterium tuberculosis. According to the WHO 2015 TB report, 9.6 million new cases of TB infection were reported globally, with approximately 1.5 million deaths (1) . One-third of the world's population is estimated to be latently infected with M. tuberculosis. The problem is further compounded by the emergence of drug-resistant forms, HIV coinfections, and diabetes (2) and the lack of an effective vaccine. Approximately 10% of the coding capacity of the M. tuberculosis genome is dedicated to the PE and PPE gene family members, so termed due to the occurrence of PE and PPE domains close to the N-terminal region (3) (4) (5) . This family is present exclusively in the genus Mycobacterium and nowhere else in the living kingdom (4) . Various PE/PPE proteins of M. tuberculosis have been reported to be expressed upon infection of macrophages and play crucial roles in virulence, antigenic diversity, and modulation of the host immune response (6) (7) (8) . Numerous members of the PE gene family display several copies of polymorphic guanine-cytosine-rich sequences (PGRSs) at their C-terminal ends in the so-called PE_PGRS subfamily (9) . In the past decade, there has been a growing interest in determining the role of PE_PGRS proteins in the pathophysiology of TB due to their limited presence in nonpathogenic mycobacteria (10) . Various members of the PE_PGRS family stimulate strong T-cell responses and immune quorum sensing (3, 11) . The PGRS domain of M. tuberculosis PE_PGRS33 (Rv1818c) is responsible for inducing humoral as well as cellular immune responses in humans, and there is also evidence for the presence of major histocompatibility complex class I (MHC-I)-restricted CD8 ϩ T cells in mice, suggesting their highly immunogenic nature (12, 13) . Mutations in the homologs of PE_PGRS62 and PE_PGRS30 in Mycobacterium marinum resulted in reduced persistence of these bacteria in granulomas (14) . In addition, the Wag22 PE_PGRS antigen-mediated immune response was found to be involved in sustaining latent infection in a mouse model of chronic M. tuberculosis infection (15) . PE_PGRS proteins have also been found at the cell surface of M. tuberculosis and found to affect its cell surface interactions (16) . Their association with the cell wall and surface exposure also leads to trafficking of these proteins out of the mycobacterial phagosome into endocytic compartments. Few of the PE_PGRS proteins are exocytosed into the extracellular environment (17) . Localization of Rv1818c to mitochondria resulted in induction of cell death (18) and was also shown to be involved in enhanced survival of Mycobacterium smegmatis in macrophages (11) . The stress pathway mediated through the endoplasmic reticulum (ER), also known as the unfolded protein response (UPR), is an alternative survival pathway which guards cells from the effects of the accumulation of misfolded/unfolded proteins (19) . However, under certain situations the UPR can cause cell death via apoptosis, which has a significant role in the pathogenicity and survival of intracellular M. tuberculosis (20) . Interestingly, ER stress is induced in specific areas of TB granulomas where apoptotic macrophages accumulate (21) . ER stress-mediated apoptosis of host cells caused by mycobacterial proteins could be a potential mechanism, which has also been reported for two important TB vaccine candidate proteins, ESAT-6 and HBHA (22, 23) . However, little is known about the effect of PE_PGRS proteins on granuloma formation and the cell death mechanism. Proteomic analysis revealed the presence of a hypothetical protein, Rv0297, in lung granulomas 90 days postinfection (24) . It is one of the T-cell antigens for which memory T cells are present in latently M. tuberculosis-infected individuals (25) . Rv0297 is present in the region of the M. tuberculosis-complex-specific genomic island as revealed by horizontal gene transfer studies (26) . It is also a part of the clusters of M. tuberculosis genes which are downregulated 4 h postinfection in human macrophages (THP-1) (27) . However, the biological role of Rv0297 remains unknown. This study was designed to delineate the mechanism of action of Rv0297 on host cell death, particularly the role of its PGRS domain alone or in fusion with the PE domain. In contrast to previous reports on the involvement of the PGRS domain in the mitochondrion-mediated cell death pathway, here we show that the Rv0297PGRS domain is associated with ER stressmediated apoptosis. The Rv0297PGRS domain is associated with upregulation of the UPR pathway within the ER. Having shown the ability of the PGRS domain to mediate ER localization, we then considered the functional consequences of such targeting. To begin with, we determined the effects of the accumulation of the Rv0297PGRS domain alone or in combination with other PE domains on host cell functions. Given the association of the ER with the UPR, the levels of UPR pathway proteins, such as GRP78 and GRP94, were investigated. We simultaneously explored the levels of ATF4 (28) and CHOP (29), well-known ER stress markers. These transcription regulators target genes that are involved in host cell apoptosis (30) . Immunoblotting studies clearly showed the overexpression of GRP78/GRP94 and ATF4/CHOP, compared to the PE control, in cells transfected with the Rv0297PGRS domain alone or in combination with other PE domains ( Fig. 1B and C) . Transfection with the Rv0297PGRS domain alone or the full-length protein in murine macrophages, as opposed to the PE domain, induced overexpression of GRP78 and GRP94 ( Fig. 1B and C, compare lanes 5 and 6 with lane 4) as evident from immunoblotting. A construct carrying the Rv0297PGRS domain fused with Rv1788PE protein ( Fig. 1B and C, lane 2) and the Rv0297PGRS domain fused with the Rv1818cPE domain ( Fig. 1B and C, lane 8) was able to induce the overexpression of these ER-specific chaperones/stress markers compared to their respective native proteins. These results were also verified by densitometric analysis (Fig. S6 ). Having shown that the effect observed was not due to the generic PGRS domain but the specific PGRS Cells with DsRed1 fluorescence (center column) were stained with ER-specific anticalnexin antibodies (left column), and merged images (right column) were created. Magnification, ϫ2,000. (B and C) PGRS domain of Rv0297 induces UPR pathway proteins and ER stress markers in transfected macrophage cells. Immunoblot analysis of murine macrophage RAW 264.7 cells with various constructs after 30 h of transfection. Total protein concentrations were determined, and equal amounts of proteins (25 g/lane) were separated by 10% SDS-PAGE and blotted onto a polyvinylidene difluoride membrane followed by immunoblotting with anti-GRP78 and anti-GRP94 antibodies (B) or anti-ATF4 and anti-CHOP antibodies (C) as markers of the UPR pathway and ER stress, respectively. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) or ␤-actin was taken as a loading control. It can be seen that GRP78/GRP94 and ATF4/CHOP were overexpressed in cells transfected with the Rv0297 PGRS domain alone or in combination with other PE domains, compared to the PE control. domain sequence present in Rv0297, which activated the ER stress pathway, we then investigated the downstream consequences of ER stress. Rv0297PGRS aids in intracellular Ca 2؉ release. The ER lumen is the key repository site of intracellular Ca 2ϩ and Ca 2ϩ binding chaperones which facilitate correct folding of native proteins (31) . Changes in the levels of Ca 2ϩ inside the ER lumen can substantially influence protein folding, leading to activation of cell death pathways. Extending the above results, we measured the levels of Ca 2ϩ ion release from the ER, in cells transfected with various constructs. After 30 h of transfection of RAW 264.7 cells, there was a noticeable increase in Ca 2ϩ after transfection with constructs carrying Rv0297PGRS ( Fig. 2A , columns 2, 3, 5, and 11). Ca 2ϩ efflux levels were not statistically significant ( Fig. 2A, columns 6 , 8, and 10) when constructs having the PGRS domain of Rv1818c were used. As variations in Ca 2ϩ dynamics perform a significant role in the ER stress-mediated cell death mechanism, these results were supportive of the idea that Rv0297PGRS might be involved in host cell apoptosis. Rv0297PGRS domain-mediated Ca 2؉ release leads to ROS generation. Perturbations in the intra-ER calcium species affect the production of ER stress-induced reactive oxygen species (ROS) (32) (33) (34) . Numerous pieces of evidence support a linkage between generation of ROS and the UPR. Oxidative stress and ROS production are vital constituents of ER stress and are not just consequences of ER stress induction. We therefore determined the cellular ROS levels using CellROX green reagent. The results show that Rv0297PGRS alone (Fig. 2B , column 2) or in combination with other PE domains (Fig. 2B , columns 3, 5, and 11) greatly enhanced the cellular ROS levels in transfected macrophages. Compared to the PE domain of Rv0297 (column 1), PGRS domain-harboring cells (columns 2, 3, 5, and 11) showed a 5-to 6-fold increase in the levels of ROS. These results therefore established a link between Rv0297PGRS expression in macrophage cells and a role in endogenous ROS regulation. The Rv0297PGRS domain also causes generation of NO. Increased levels of nitric oxide (NO) disturb ER calcium pump activity and also induce ER stress (35) . We accordingly investigated the role of the PGRS domain present within Rv0297 in NO synthesis. As could be seen, RAW 264.7 cells transfected with a plasmid carrying the Rv0297PE domain alone secreted basal levels of NO (Fig. 2C , column 1), while production of NO was significantly increased~15-fold in those expressing Rv0297PGRS alone or as a full-length protein or fused with Rv1788 PE protein or the Rv1818c PE domain (columns 2, 3, 5, and 11, respectively). These results further indicated an involvement of the Rv0297PGRS domain in activation of cellular NO species. Nuclear blebbing was observed in cells transfected with the Rv0297PGRS domain. Macrophage cells were transfected with various constructs, including the PE domain of Rv0297 (Fig. 3A, row b) , Rv0297PGRS (row c), or full-length Rv0297 protein (row d), along with vector control (row a) for 30 h. DsRed1 signals when merged with Hoechst nuclear stain highlighted the relative intensity of blebbing (right column). Analyses of fluorescence microscopic images revealed that wherever Rv0297PGRS was present, a clear nuclear blebbing was observed (Fig. 3A) . These results clearly pointed to the role of the PGRS domain of Rv0297 in the induction of apoptosis in macrophage cell lines. Caspase-8 is activated by the Rv0297PGRS domain. To confirm the results described above, activation of caspase-8, the effector caspase of the cellular apoptotic pathway, was assayed. Caspase-8 activity was measured in RAW 264.7 cells transfected with the Rv0297PGRS construct. A significantly elevated (5-to 6-fold) increase in the levels of caspase-8 could be seen in RAW 264.7 cells transfected with constructs carrying Rv0297PGRS alone (Fig. 3B , column 3) or full-length Rv0297 (column 4) or Rv0297 fused with the PE domain of Rv1818c (column 5) or Rv1788 PE protein (column 11). The activation of caspase-8 noticed for constructs carrying the PGRS domain of Rv1818c (column 6) or full-length Rv1818c (column 7) or Rv1818c fused to the PE domain of Rv0297 (column 8) or Rv1788 PE protein (column 12) was expected since the Rv1818c PGRS domain is known to cause apoptosis involving the mitochondrial apoptotic pathway. Interestingly, the PE domain alone of Rv0297 (column 2) or Rv1788 PE protein (column 10) was unable to cause caspase-8 activation. These results demonstrated that induction of apoptosis by Rv0297PGRS is facilitated by classical caspase-8-dependent pathways, similarly to Rv1818c-mediated apoptosis (36) . In vivo localization of Rv0297 from M. tuberculosis H 37 Rv to the ER of RAW 264.7 cells. Colocalization of PE_PGRS5 to the ER in pDsRed1C1-Rv0297-transfected cells was observed in previous results. To validate this, in vivo colocalization of Rv0297 to the ER of infected macrophages was performed. Uninfected RAW 264.7 cells were used as a negative control. Confocal microscopy analysis revealed the colocalization of M. tuberculosis Rv0297 into the ER of infected macrophages (Fig. 4) . Merged images clearly indicated the presence of Rv0297 in the ER of infected macrophages (Fig. 4 , upper panels) but not in uninfected ones (Fig. 4 , lower panels). These observations corroborate our previous results where overexpression of Rv0297 was done in HEK293T cells. Consistent with transfection results, the recombinant protein also efficiently induced cell death of RAW 264.7 macrophages. The viability of cells was significantly decreased in the presence of rRv0297PGRS as evident from alamarBlue assay (Fig. 5A ). rRv0297PGRS protein (2.5 g/ml) induced 12.42% Ϯ 1.26% mortality of macrophage cells, and this increased as a direct function of the protein concentration: 16.24% Ϯ 4.29%, 22.79% Ϯ 2.90%, 35.48% Ϯ 0.98%, and 36.07% Ϯ 0.61% with 5, 10, 20, and 40 g/ml rRv0297PGRS, respectively. These results demonstrate that administration of the rRv0297PGRS protein to macrophages leads to cell death in a dose-dependent manner. In order to determine the mode of cell death, annexin V/propidium iodide (PI) staining of rRv0297PGRS-treated RAW 264.7 cells was used. Both early and late apoptosis of cells were observed upon stimulation with 20 and 40 g/ml of rRv0297PGRS protein (Fig. 5B ). The total percentage of apoptotic cells induced by 20 and 40 g/ml of rRv0297PGRS was found to be 43.83% Ϯ 2.54% and 58.37% Ϯ 7.92%, respectively (Fig. 5C ), indicating that Rv0297PGRS is inducing apoptosis rather than necrosis. The effect of rRv0297PGRS protein on generation of ROS and NO and their involvement in ER stress was also assessed. Increased production of ROS and NO macrophages was observed upon stimulation with rRv0297PGRS protein. Fold increments of 1.53 Ϯ 0.22 and 1.80 Ϯ 0.13 in the levels of ROS were observed compared to untreated cells (Fig. 5D ). Almost a 2-fold increase in the levels of NO was observed in cells stimulated with 20 g/ml of rRv0297PGRS protein in contrast to unstimulated cells (Fig. 5E ). Bovine serum albumin (BSA) (80 g/ml) induced negligible amounts of ROS and NO. These results support the role of Rv0297PGRS in ER stress-mediated apoptosis. In silico analysis revealed preferential binding of Rv0297 with TLR4 compared to TLR2. Host-pathogen cross talks mostly involve pathogen recognition receptors (PRRs) or Toll-like receptors (TLRs) of host and pathogen-associated molecular patterns (PAMPs) of microbes. Few of the PE_PGRS proteins of M. tuberculosis were known to interact with TLRs, mainly the TLR2 type. We therefore carried out in silico simulations to predict the identity of the TLR involved in interaction. The I-TASSER-predicted model was chosen on the basis of the high confidence score (C-score), 0.26, and the PSVS results. Verify 3D (reference range, Ϫ1 to ϩ1) also confirmed the good model quality with a score of 0.56. The modeled three-dimensional (3D) structure was refined through molecular dynamics (MD) simulations and showed a stable trajectory from 6.5 to 20 ns with small root mean square deviations (RMSDs) between 0.8 nm and 1 nm (Fig. 6A) . The model quality significantly improved with a Verify 3D score of 0.58 and overall 98.8% residues in the favored or allowed region in comparison to 97.7% before MD. Only 1.2% amino acids were in the disallowed region (earlier value, 2.4%). MD refinement of the TLR4 structure also shows a stable trajectory from 30 ns to 50 ns (Fig. 6B) . Protein-protein docking was carried out using HADDOCK 2.2 for comparative analysis of Rv0297 interaction with TLR4 versus TLR2. We obtained a Ϫ13.9 Ϯ 14.5 HADDOCK score for TLR4 in comparison to Ϫ5.3 Ϯ 17.8 for TLR2. HADDOCK clustered 106 structures into 11 clusters for the TLR4-Rv0297 complex, while it clustered 41 structures into 7 clusters for the TLR2-Rv0297 complex. The TLR4-Rv0297 complex structure was taken from cluster 2 on the basis of the best Z score, Ϫ2.2, while the TLR2-Rv0297 complex from the best cluster, 1, had a Z score of Ϫ1.6. These in silico studies suggested that TLR4 may have a better interaction with Rv0297 than does TLR2. The TLR4-Rv0297 complex was then analyzed for stable interaction after long-term MD simulations. The backbone RMSD of the TLR4-Rv0297 complex shows a stable trajectory of 50 ns (Fig. 6D) . Comparison of pre-MD and post-MD complexes by LigPlotϩ (37) revealed conserved interactions in the form of hydrogen or hydrophobic bonds (Table S4). The surface interaction diagram of the TLR4-Rv0297 complex also disclosed a perfectly fitting interaction between TLR4 and Rv0297 (Fig. 6C) . Having predicted TLR4 interaction with Rv0279PGRS, we demonstrated the same by in vitro interaction assay. RAW 264.7 and RAW-ΔTLR4 cells were incubated with rRv0297PGRS protein and immunostained with anti-Rv0297PGRS sera to assess the interaction of TLR4 localized on the surface of macrophages with Rv0297PGRS protein. Confocal imaging revealed the presence of the interacting complex on the surface of RAW 264.7 cells incubated with rRv0297PGRS protein (Fig. 6E, row 2) . No such interaction was observed without any protein (Fig. 6E, row 1) . In contrast to RAW 264.7 cells, RAW-ΔTLR4 cells did not display such complexes on their surface even in the presence of rRv0297PGRS protein ( Fig. 6E, row 4) . The presence of this interaction on the surface of RAW 264. 7 The downstream effect of action of the PGRS domain is TLR4 dependent. RAW-ΔTLR4 cells were stimulated with 0 to 40 g/ml of rRv0297PGRS protein for 48 h, and the cell survival was assessed by alamarBlue assay. No significant change in cell survival was observed upon rRv0297PGRS protein treatment in comparison to untreated cells and 50-g/ml BSA control (Fig. 7A) . All cells were healthy and alive even upon 40-g/ml rRv0297PGRS protein stimulation. On the other hand, decreased viability was observed in the case of 0.1% SDS treatment (positive control). The results clearly indicate the role of TLR4 in the effect of Rv0297PGRS protein. The levels of ROS and NO were then assessed once again to further validate the role of TLR4 in rRv0297PGRS protein-stimulated ER stress-mediated apoptosis. To investigate this, RAW-ΔTLR4 cells were treated with rRv0297PGRS (0 to 40 g/ml) and the levels of ROS and NO were assayed and found to be exactly similar in both treated and untreated cells ( Fig. 7B and C, respectively) . In contrast to this, 1 M thapsigargin (positive control) efficiently induced almost 50% increase in both ROS and NO levels from RAW-ΔTLR4 cells. These results demonstrate that the downstream effects of Rv0297PGRS protein, such as cell death and ROS and NO production, are TLR4 dependent. Intracellular pathogens often act to prevent host cell apoptosis to reduce the immune response. Previous studies revealed that apoptosis is a macrophage defense process to counter M. tuberculosis infection (38) . Despite the fact that M. tuberculosis is able to inhibit apoptosis in macrophages, it is still capable of inducing apoptosis (39) . Apoptotic macrophages have been found within M. tuberculosis-infected lung granulomas (40, 41) . These findings generate considerable interest in the potential role of apoptosis in M. tuberculosis infection and subsequent pathogenesis. Previous studies also emphasized that ER stress-mediated apoptosis plays an important role in TB pathogenesis (22, 23) . ER stress has been found to induce apoptosis in TB granulomas, in zones where apoptotic macrophages accumulate in mice and humans (21) . The appearance of Rv0297 protein in TB granuloma 90 days postinfection and its absence during initial phases of infection (30 days postinfection) point to its role in granuloma maintenance (24) . However, the exact mechanism by which it helps in M. tuberculosis survival in the human host and maintenance of disease is not known. Using various fusion constructs and colocalization studies, we showed that Rv0297 protein localizes to the ER, independent of its PE domain. In order to decipher the biological effects of localization within the ER of the host cells, we investigated the ER stress pathways in cells harboring various domains of Rv0297 as well as another PE protein (Rv1788) and PE_PGRS protein (Rv1818c). Our results show that the Rv0297 PGRS domain is involved in ER stress. Notably, it has been shown that when stress is excessively high, the UPR pathway shifts from a prosurvival to a prodeath mode (42) . GRP78 and GRP94 are well-established ER chaperone molecules which have a crucial role in maintaining cell viability under numerous cell stresses (43) . Under pathological conditions, altered levels of ATF4 and CHOP-crucial proteins expressed during ER stress-induced apoptosis-have been shown (22) . Using GRP94, GRP78, ATF4, and CHOP as standard UPR markers for ER stress, we showed upregulation of these markers as a consequence of Rv0297PGRS domain expression in macrophages. Together, these data suggest that, in later stages of M. tuberculosis infection, the ER stress sensors might be upregulated for induction of apoptosis as a reaction to M. tuberculosis infection. Disruption of Ca 2ϩ homeostasis is a contributing element of ER stress. We therefore also investigated the possible involvement of the Rv0297PGRS-mediated Ca 2ϩ ion release in ER stress in RAW 264.7 cells. The levels of intracellular Ca 2ϩ ions in RAW 264.7 cells were shown to be elevated using a fluorescent Ca 2ϩ indicator. ROS and NO generation represent essential mechanisms invoked by macrophages to control M. tuberculosis (44) . These antimicrobial species play a critical role in controlling intracel-lular bacteria. Interestingly, ER stress pathways are also induced by ROS or NO (45) . Our results indicate that Rv0297PGRS induces both NO and ROS production in RAW 264.7 cells. We also showed that ectopic expression of the PGRS domain affects the levels of caspase-8, an effector caspase usually activated by the extrinsic apoptotic pathway (46) . Caspase-8 and caspase-3 are the key caspases involved in ER stress-mediated apoptosis (47) . Pathogens use intrinsically disordered proteins to perturb and hijack host cell networks for a productive infection (48, 49) . An increased disordered region present in Rv0297PGRS might be responsible for differences in organelle localization and function from Rv1818cPGRS. Fusion of Rv0297PGRS with other PE domains allowed them to be targeted to the ER. Protein sequence alignment revealed some extra sequences in the Rv0297PGRS domain, possibly containing motifs responsible for ER localization. Rv1818c, a well-characterized protein known to be localized to mitochondria, while sharing some common motifs apparently lacks a few of them (see Table S3 in the supplemental material). While some ER signal sequences are present in Rv0297PGRS, we were unable to identify the precise signature motifs which enable exclusive targeting to the ER. It is presumable that these sequences might not be responsible for organelle-specific targeting but could be responsible for organelle targeting generically. Furthermore, some of these signals might be responsible for retention of Rv0297PGRS in the host ER, leading to generation of the stress response and induction of caspase-8-mediated apoptosis. Very little is known about the functional role of PE_PGRS proteins, as only a few PE_PGRS proteins have been assigned a physiological function so far (3) . For example, the PE domain Rv1818c, a well-studied protein of this family, induces primary necrosis, whereas the PGRS domain is involved in induction of mitochondrion-mediated apoptosis. Recently, it was demonstrated that PE_PGRS30 is required for virulence, as its deletion from M. tuberculosis compromised its capacity to colonize in lung tissue and hampered tissue damage. The inactivation of the gene also leads to inhibition of phagolysosome fusion in macrophages (50) . Similarly, M. tuberculosis PE_PGRS17 promotes the death of host macrophages as well as increases in the secretion of the proinflammatory cytokine tumor necrosis factor alpha (TNF-␣) (51) . On the other hand, PE_PGRS62 reduces phagolysosome maturation and induces secretion of gamma interferon (IFN-␥) (52) . Induction of ER stress-mediated cell death via apoptosis has been described in several diseases, including vascular diseases (53) . The ER stress pathway was initially recognized as a cellular mechanism activated by the accumulation of unfolded proteins in the ER to maintain proper ER functions. In addition, ER stress pathways are also involved in protection of cells by several other cellular stresses. However, when stress reaches a critical level, apoptosis is induced to eliminate injured cells (54) . In later stages of M. tuberculosis infection, the consequence of cell death favors the pathogen. A mature TB granuloma in human patients has slight vascularization and a limited access to immune cells such as macrophages and lymphocytes (25) . The stimulation of cellular apoptosis in the infected and foamy macrophages, which are profusely present, will cause accumulation of caseum and the development of disease pathology. M. tuberculosis takes benefit of this pathology by liquefaction and cavitation of the granuloma, in order to transmit infectious bacilli into the airways. As a result, whereas apoptosis could be advantageous to the host at low bacterial numbers at early stages of infection, it is harmful after the disease has advanced and fibrocaseous granulomas have been generated. Mycobacterial phenolic glycolipid-dependent production of CCL2 in human alveolar macrophages serves as a source for recruitment of permissive macrophages to sustain and disseminate the mycobacterial infection (55) . The earliest step in an infection is the interaction of host cell and microbial components. Macrophages sense or recognize specific conserved moieties of microbes called PAMPs with the help of PRRs expressed on their surface. Interaction of TLRs of macrophages with different M. tuberculosis ligands is required for the modulation of various cellular events such as apoptosis of macrophages (56) , production of ROS and NO intermediates (57) , antigen presentation (58) , and phagolysosomal fusion (59) . Pathogenesis and virulence of M. tuberculosis depend on TLR1, TLR2, TLR4, and TLR9 and their signaling cascades (60) . Induction of host cell apoptosis through TLR2 signaling has been shown by stimulation with mycobacterial ligands and live M. tuberculosis. Few of the PE/PPE proteins such as Rv1818c of M. tuberculosis have been shown to stimulate host cell apoptosis in a TLR2-dependent manner (36) . Maturation and activation of dendritic cells by PE_PGRS11 and PE_PGRS17 were shown to be dependent on TLR2 signaling (61) . However, it is yet to be established whether TLR2 interaction is the common property of PE/PPE proteins. The complex interactive network of PE/PPE proteins with host cell receptors modulates the pathogenesis of disease. In contrast, an important role of both TLR2 and TLR4 in induction of macrophage apoptosis upon M. tuberculosis infection has also been shown. TLR4-dependent signaling is important in maintaining the balance between necrosis and apoptosis of macrophages (62) . We report for the first time the interaction of TLR4 with Rv0297PGRS of M. tuberculosis in execution of host cell apoptosis. In silico analysis revealed that Rv0297 interacts with TLR4, which was also corroborated by our in vitro TLR4 and Rv0297PGRS interaction assay. Confocal images of this interaction clearly indicated the presence of interacting complexes on the surface of RAW 264.7 cells, but they were totally absent on the cells lacking TLR4. Also, the induction of ROS and NO from macrophages was found to be dependent on TLR4 interaction. RAW cells lacking TLR4 were unable to undergo apoptosis upon stimulation with rRv0297PGRS. All these results supported the functioning of Rv0297 in a TLR4-dependent manner. Another mycobacterial PE protein complex, PE9/PE10, has also been shown to induce apoptosis of macrophages via its interaction with TLR4 (63) . The involvement of the PE_PGRS protein Rv0297 in ER stress-mediated induction of apoptosis in macrophages points to a novel role for the PGRS domain of this family of proteins. We propose that Rv0297 expression in human cells leads to its ER localization, followed by disruption of intracellular Ca 2ϩ levels and induction of ROS and NO production (Fig. 8) . The ensuing ER stress response culminates in caspase-8 activation and induces apoptosis in macrophage cells. Apoptosis of infected macrophages in advanced stages of granulomas would support dissemination of the bacteria through apoptotic bodies, thereby playing an important role in spreading the disease. The molecular dissection of the different domains within the PE_PGRS protein of the PE/PPE family, present exclusively in the genus Mycobacterium, and their possible functional significance will aid in a better understanding of the virulence and pathogenesis of M. tuberculosis. Future experiments with lung granulomas will likely provide a better insight into the mode of action. In silico analyses of Rv0297. The GlobPlot 2.3 (http://globplot.embl.de/cgiDict.py) tool was used to predict disordered regions and globularity in the PE_PGRS proteins used in this study. The presence of putative ER localization signals in Rv0297 protein was searched from LocSigDB (64) using the FIMO tool (65) from MEME Suite version 4.11.2 (http://meme-suite.org/tools/fimo) with a P value of Ͻ0.001. The sequence file of the amino acid sequence of Rv0297 was uploaded under "Input the sequences," and ER signal sequences were typed under "Input the motifs" in the data submission form. Protein structure modeling, validation, and protein-protein docking. The Rv0297 protein sequence was obtained from the UniProtKB database (identifier [ID] Q6MX50) (66) , and the 3D structure was modeled using I-TASSER (67) . The human TLR4 protein structure (27 to 627 amino acids [aa]) was obtained from the RCSB protein data bank (PDB ID 3FXI). Rv0297 interactions with TLR4 and TLR2 were compared through the HADDOCK 2.2 web server (68) . Structure validation was done through PSVS (69) , and molecular dynamics (MD) simulations were performed with GROMACS v5.0.7 (70) (71) (72) (73) (74) . Generation of constructs. Genomic DNA of H 37 Ra was provided by Astrid Lewin, RKI, Berlin, Germany. The PE18 gene (Rv1788), Rv0297, Rv1818c full-length gene, and their PE and PGRS domains were cloned as C-terminal fusions into pDsRed1C1 vector. The primer sequences and PCR conditions are listed in Tables S1 and S2 in the supplemental material, respectively. A complete list and description of the constructs prepared are shown in Fig. S1 . For generating rRv0297PGRS protein, the Rv0297PGRS gene was cloned in the pET28a expression vector and expressed in BL21(DE3)pLysS cells. Recombinant protein was purified from inclusion bodies by solubilization in 8 M urea in phosphate-buffered saline (PBS) (pH 7.5) and on-column renaturation using a urea gradient followed by Ni 2ϩ -nitrilotriacetic acid (NTA) chromatography. The protein was treated with polymyxin B at 4°C for 2 h. Immunization and antibody generation. A rabbit was immunized with 500 g/ml of rRv0297PGRS, 3 booster doses of immunization were given at an interval of 15 days, and sera were collected 15 days after the last immunization. Cell culture and transfection. The murine macrophage cell lines RAW 264.7 and RAW-ΔTLR4 and HEK293T (human embryonic kidney cells) were maintained in Dulbecco's modified Eagle's medium (Invitrogen) supplemented with 10% fetal bovine serum (FBS) (Invitrogen) and penicillin-streptomycin (Pen-Strep) (100 g/ml). Cells (5,000/well) were seeded in a 96-well plate. The constructs were transfected into HEK293T cells for localization studies and in the RAW 264.7 cell line for ER stress-related experiments using Lipofectamine 3000 (Invitrogen) for 30 h and 48 h. Immunofluorescence staining. HEK293T, RAW 264.7, and RAW-ΔTLR4 cells were seeded on coverslips to adhere for 16 to 18 h and transfected with constructs shown in Fig. S1 or stimulated with different concentrations of rRv0297PGRS. After the required period of incubation, cells were fixed with 3.7% paraformaldehyde for 15 min, permeabilized, and blocked with 0.1% Triton X-100 -1% BSA-PBS for 10 min. The cells were incubated with the respective primary and Alexa Fluor-conjugated secondary antibodies (75) . Images were visualized using an Olympus FluoView FV1200 laser scanning confocal microscope. In Cell death assay. RAW 264.7 and RAW-ΔTLR4 cells were treated with rRv0297PGRS (2.5 to 40 g/ml) for 48 h. Ten percent resazurin sodium salt (100 mg/ml) was added. Absorbance was monitored at 570 nm with reference at 600 nm after 24 h of reaction. For annexin V/PI staining, treated cells (24 h) were stained with the Alexa Fluor 488 annexin V and PI flow cytometry kit (Invitrogen). Analysis of 10,000 stained cells was performed with a FACSAria 3 cytometer (BD Biosciences, USA). Calcium influx assay. Transfected cells were stained with Fluo-4 NW dye solution (Molecular Probes) at 37°C for 30 min and further incubated at room temperature (RT) for 30 min. Calcium release was detected by measuring the fluorescence intensity using a 494-nm excitation and a 516-nm emission wavelength. Detection of ROS. Transfected or treated cells were stained with 5 mM CellROX green reagent for 30 min at 37°C. Plates were read using a 485-nm excitation and a 520-nm emission wavelength. Results were represented as the mean fluorescent intensity of cells after subtraction of the blank (CellROX green reagent-unstained cells). NO quantitation in macrophages. Cells were either transfected with various combination of PE and PGRS gene constructs or treated with 20 to 80 g/ml of rRv0297PGRS protein for 30 h. Cell-free supernatant (150 l) was mixed with 50 l of Griess reagent for 30 min. The nitrite concentration was measured using sodium nitrite as a standard. Plates were read at 540 nm. Western blot analysis. Western blotting was performed with anti-CHOP, anti-ATF4, anti-GRP94, and anti-GRP78/BiP antibodies (Abcam, USA) and anti-␤-actin (Santa Cruz Biotechnology, USA). Membranes were developed with a chemiluminescent reagent (Millipore). Caspase assay. The activity level of initiator caspase-8 was determined using the Caspase-Glo assay (Promega). Statistical analysis. All data were expressed in the form of means Ϯ standard deviations (SDs) derived from 3 different groups of independent experiments using GraphPad Prism5 software. A one-way analysis of variation (ANOVA) was performed followed by Dunnett's post hoc test in order to calculate the statistical significance at a P value of Ͻ0.05. Supplemental material for this article may be found at https://doi.org/10.1128/mBio .01017-18. We declare no conflict of interest. World Health Organization The co-management of tuberculosis and diabetes: challenges and opportunities in the developing world The PE/PPE multigene family codes for virulence factors and is a possible source of mycobacterial antigenic variation: perhaps more? Transcriptional regulation of Mycobacterium tuberculosis PE/PPE genes: a molecular switch to virulence? Clusters of PE and PPE genes of Mycobacterium tuberculosis are organized in operons: evidence that PE Rv2431c is co-transcribed with PPE Rv2430c and their gene products interact with each other Contrasting function of structured N-terminal and unstructured C-terminal segments of Mycobacterium tuberculosis PPE37 protein Mycobacterium tuberculosis co-operonic PE32/PPE65 proteins alter host immune responses by hampering Th1 response Comparative genomic and proteomic analyses of PE/PPE multigene family of Mycobacterium tuberculosis The PE multigene family: a "molecular mantra" for mycobacteria Roles of PE_PGRS family in Mycobacterium tuberculosis pathogenesis and novel measures against tuberculosis Expression of the PE_PGRS 33 protein in Mycobacterium smegmatis triggers necrosis in macrophages and enhanced mycobacterial survival Evaluation of T-cell responses to peptides with MHC class I-binding motifs derived from PE_PGRS 33 protein of Mycobacterium tuberculosis Detection of interferon gammasecreting CD8(ϩ) T lymphocytes in humans specific for three PE/PPE proteins of Mycobacterium tuberculosis Granuloma-specific expression of Mycobacterium virulence proteins from the glycine-rich PE-PGRS family The PGRS domain of Mycobacterium tuberculosis PE_PGRS Rv1759c antigen is an efficient subunit vaccine to prevent reactivation in a murine model of chronic tuberculosis Evidence that mycobacterial PE_PGRS proteins are cell surface constituents that influence interactions with other cells Mycobacterial surface moieties are released from infected macrophages by a constitutive exocytic event Induction of cell death after localization to the host cell mitochondria by the Mycobacterium tuberculosis PE_PGRS33 protein Human resistin, a proinflammatory cytokine, shows chaperone-like activity Endoplasmic reticulum stress pathway-mediated apoptosis in macrophages contributes to the survival of Mycobacterium tuberculosis Induction of ER stress in macrophages of tuberculosis granulomas Endoplasmic reticulum stress response is involved in Mycobacterium tuberculosis protein ESAT-6-mediated apoptosis Mycobacterial HBHA induces endoplasmic reticulum stress-mediated apoptosis through the generation of reactive oxygen species and cytosolic Ca2ϩ in murine macrophage RAW 264.7 cells Portrait of a pathogen: the Mycobacterium tuberculosis proteome in vivo Memory T cells in latent Mycobacterium tuberculosis infection are directed against three antigenic islands and largely contained in a CXCR3(ϩ)CCR6(ϩ) Th1 subset Contribution of horizontally acquired genomic islands to the evolution of the tubercle bacilli Global transcriptional profile of Mycobacterium tuberculosis during THP-1 human macrophage infection Transcription factor ATF cDNA clones: an extensive family of leucine zipper proteins able to selectively form DNA-binding heterodimers Stressinduced binding of the transcriptional factor CHOP to a novel DNA control element Upregulation of CHOP/GADD153 during coronavirus infectious bronchitis virus infection modulates apoptosis by restricting activation of the extracellular signal-regulated kinase pathway ER stress-induced cell death mechanisms An involvement of oxidative stress in endoplasmic reticulum stress and its associated diseases In vitro cultured Spodoptera frugiperda insect cells: model for oxidative stress-induced apoptosis The baculovirus antiapoptotic p35 gene also functions via an oxidantdependent pathway Nitric oxide and endoplasmic reticulum stress Execution of macrophage apoptosis by PE_PGRS33 of Mycobacterium tuberculosis is mediated by Toll-like receptor 2-dependent release of tumor necrosis factor-alpha LigPlotϩ: multiple ligand-protein interaction diagrams for drug discovery Apoptosis is an innate defense function of macrophages against Mycobacterium tuberculosis Mycobacterium tuberculosis PE25/PPE41 protein complex induces necrosis in macrophages: role in virulence and disease reactivation? Apoptosis of macrophages and T cells in tuberculosis associated caseous necrosis Ipr1 gene mediates innate immunity to tuberculosis Cellular stress responses: cell survival and cell death From endoplasmic-reticulum stress to the inflammatory response Mycobacterial survival strategies in the phagosome: defence against host stresses Endoplasmic reticulum stress induced by oxidative stress in retinal pigment epithelial cells RIPK1 promotes death receptor-independent caspase-8-mediated apoptosis under unresolved ER stress conditions ER stress induces caspase-8 activation, stimulating cytochrome c release and caspase-9 activation Intrinsically unstructured proteins How pathogens use linear motifs to perturb host cell networks PE_PGRS30 is required for the full virulence of Mycobacterium tuberculosis Mycobacterium tuberculosis PE_P-GRS17 promotes the death of host cell and cytokines secretion via Erk kinase accompanying with enhanced survival of recombinant Mycobacterium smegmatis Phagolysosome maturation of macrophages was reduced by PE_PGRS 62 protein expressing in Mycobacterium smegmatis and induced in IFN-␥ priming Mechanisms of ER stress-induced apoptosis in atherosclerosis Apoptosis and aging: increased resistance to apoptosis enhances the aging process Phenolic glycolipid facilitates mycobacterial escape from microbicidal tissue-resident macrophages Regulation of antigen presentation by Mycobacterium tuberculosis: a role for Toll-like receptors Sensing of Mycobacterium tuberculosis and consequences to both host and bacillus Mycobacterium tuberculosis modulates macrophage lipid-sensing nuclear receptors PPAR␥ and TR4 for survival Innate immune recognition of Mycobacterium tuberculosis Interaction of pattern recognition receptors with Mycobacterium tuberculosis PE_PGRS antigens of Mycobacterium tuberculosis induce maturation and activation of human dendritic cells Role of TLR2-and TLR4-mediated signaling in Mycobacterium tuberculosis-induced macrophage death The Mycobacterium tuberculosis protein pair PE9 (Rv1088)-PE10 (Rv1089) forms heterodimers and induces macrophage apoptosis through Toll-like receptor 4 LocSigDB: a database of protein localization signals FIMO: scanning for occurrences of a given motif UniProt: the universal protein knowledge base The I-TASSER suite: protein structure and function prediction The HADDOCK2.2 web server: user-friendly integrative modeling of biomolecular complexes Evaluating protein structures determined by structural genomics consortia The GROMACS development team. GROMACS user manual version 5 Alanine mutation of the catalytic sites of pantothenate synthetase causes distinct conformational changes in the ATP binding region Insight into the inhibitor discrimination by FLT3 F691L Mechanistic principles behind molecular mechanism of rifampicin resistance in mutant RNA polymerase beta subunit of Mycobacterium tuberculosis GROMOS force field Role of UDP-N-acetylglucosamine-2-epimerase/ N-acetylmannosamine kinase (GNE) in beta1-integrin-mediated cell adhesion