Abstract
Purpose: The goal of this study is to enhance the efficacy of imipridones, a novel class of AKT/ERK inhibitors that displayed limited therapeutic efficacy against glioblastoma (GBM).
Experimental Design: Gene set enrichment, LC/MS, and extracellular flux analyses were used to determine the mechanism of action of novel imipridone compounds, ONC206 and ONC212. Orthotopic patient-derived xenografts were utilized to evaluate therapeutic potency.
Results: Imipridones reduce the proliferation of patient-derived xenograft and stem-like glioblastoma cell cultures in vitro and in multiple xenograft models in vivo. ONC212 displayed the highest potency. High levels of c-myc predict susceptibility to growth inhibition and apoptosis induction by imipridones and increased host survival in orthotopic patient-derived xenografts. As early as 1 hour, imipridones elicit on-target inhibition, followed by dephosphorylation of GSK3β at serine 9. GSK3β promotes phosphorylation of c-myc at threonine 58 and enhances its proteasomal degradation. Moreover, inhibition of c-myc by BRD4 antagonists sensitizes for imipridone-induced apoptosis in stem-like GBM cells in vitro and in vivo. Imipridones affect energy metabolism by suppressing both glycolysis and oxidative phosphorylation, which is accompanied by a compensatory activation of the serine-one carbon-glycine (SOG) pathway, involving the transcription factor ATF4. Interference with the SOG pathway through novel inhibitors of PHGDH results in synergistic cell death induction in vitro and in vivo.
Conclusions: These results suggest that c-myc expression predicts therapeutic responses to imipridones and that imipridones lead to suppression of tumor cell energy metabolism, eliciting unique metabolic vulnerabilities that can be exploited for clinical relevant drug combination therapies. Clin Cancer Res; 24(21); 5392–406. ©2018 AACR.
Imipridones have recently entered clinical testing for hematologic and solid malignancies, including glioblastoma, the most common primary brain tumor that requires novel treatments. In the current study, we provide evidence that chemically modified imipridones, ONC206 and ONC212, exert more potent antiglioma activity than the original lead compound, ONC201. Moreover, we demonstrate that high levels of c-myc predict therapeutic responses in preclinical model systems of glioblastoma. Notably, single treatment of ONC212 significantly extends host survival in an orthotopic patient-derived glioblastoma xenograft model. Combined treatment with a clinically validated BRD4 antagonist, OTX015, and imipridones induces enhanced reduction of glioma growth in vitro and in vivo. These observations support the notion that chemically modified imipridones display preclinical activity singly or as part of combination therapies in model systems of glioblastoma. Therefore, clinical testing of these novel compounds either singly or in combination therapies is warranted.
Introduction
Imipridone derivatives are a new class of molecules that have entered clinical testing for solid and nonsolid malignancies. Representative examples are the lead compound, called ONC201 (1), and the chemically modified and presumably enhanced derivatives, ONC206 and ONC212 (2). The original compound, ONC201, was identified in a screen for molecules that induce TRAIL, a cytokine that in the field of oncology was considered as the holy grail for cancer therapy because it specifically kills tumor cells in a rapid manner, but not nonneoplastic cells (1). On the molecular level, ONC201 elicited an increase of TRAIL transcription through a FOXO3A-dependent mechanism, involving the inhibition of the kinases ERK and AKT (1). Early after the discovery of ONC201, it became evident that this compound exerts additional anticancer properties, e.g., the inhibition of proliferation in cells that do not undergo apoptosis upon ONC201 administration (3–5). These additional features are in part attributed to the fact that imipridones inhibit dopamine receptor signaling, such as DRD2 (6).
Tumor cell metabolism has become a central focus in cancer research because it bears the potential to target metabolic aberrancies in tumors, thus providing a therapeutic window. Although cancer cells often depend on glycolysis in the presence of abundant oxygen (called Warburg effect), it has become evident that other metabolic pathways are targetable in malignant cells, even oxidative phosphorylation (OXPHOS), the most efficient form of energy generation in cells. Quite prominent is the metabolic dependency on glutamine by many malignant tumors (7–9), which requires the presence of a functional OXPHOS. Being a master regulator of tumor stem cells, the transcription factor c-myc (10–12) is at the core of regulation of tumor cell metabolism and affects the expression of enzymes related to glycolysis, glutamine, and OXPHOS.
In this report, we provide evidence that novel imipridones with different potency are highly effective against state-of-the-art model systems of glioblastoma, involving in particular glioma stem-like cells. We have unraveled a novel mechanism by which imipridones elicit their antiglioma activity, involving suppression of the master regulator, c-myc, as early as 1 hour, and shut down of tumor energy cell metabolism, leading to a state of energy deprivation and tumor cell cytostasis and apoptosis. Moreover, high levels of c-myc predict both apoptotic response and inhibition of proliferation of GBM cells to imipridone derivatives and extended host survival in orthotopic patient-derived xenografts (PDX) of glioblastoma, positioning c-myc as a potential therapeutic predictive marker, and that further inhibition of c-myc along with imipridones induces synergistic growth reduction in vitro and in vivo. In addition, we have unraveled that upon imipridone-induced energy deprivation, the serine-one carbon-glycine pathway is increased, eliciting a novel unique metabolic vulnerability.
Materials and Methods
Reagents
Imipridones, ONC201, ONC206, and ONC212, were received from Oncoceutics. OTX015 was purchased from Selleckchem. A 10 mmol/L working solution in dimethylsulfoxide (DMSO) was prepared for all reagents prior to storage at −20°C. Final concentrations of DMSO were below 0.1% (v/v). The plasmids for wild-type c-myc and mutant c-myc were obtained from Addgene (ID: 45597 and 45598).
Cell cultures and growth conditions
All cells were cultured as described (13–17). The identities of the cell cultures were confirmed by the respective source of purchase. LN229 and U87 cells are established glioblastoma cells. SF188 is an MYC-amplified pediatric glioblastoma cell culture. NCH644, NCH690, and NCH421K glioma stem-like cells were cultured in MG-43 medium (CLS) for both maintenance and experiments (13–15, 18). GBM12 and GBM14 are PDX tumors as described elsewhere (13–17). Human astrocytes were obtained from ScienCell Research Laboratories, Inc. and cultured as recommended by the provider.
Cell viability assays
In order to examine cellular proliferation, CellTiter-Glo assays were performed as previously described. ATP levels were determined as performed in refs. 13–17.
Measurement of apoptosis and mitochondrial membrane potential
Annexin V/propidium iodide, propidium iodide, and tetramethylrhodamine, ethyl ester (TMRE) staining (for mitochondrial membrane potential) were performed as previously described (13–17) or in accordance with the manufacturer instructions for TMRE staining (Cell Signaling Technology). The data were analyzed with the FlowJo software (version 8.7.1; Tree Star).
Extracellular flux analysis
Extracellular flux analysis was performed on the Seahorse XFe24 Analyzer. The “XF cell mito stress test kit” (Agilent Technologies) was utilized to determine parameters relevant to oxidative phosphorylation and determined as described earlier in ref. 19. GBM cells were incubated with Seahorse XF base medium supplemented with 5 mmol/L glucose, 1 mmol/L pyruvate, and 2 mmol/L l-glutamine in a CO2-free incubator for 1 hour before the assay. During the course of the assay, 2 μmol/L oligomycin, 2 μmol/L FCCP, and 0.5 μmol/L rotenone/antimycin were added sequentially. Regarding glycolysis, the “XF cell glycolysis stress test kit” (Agilent Technologies) was used in accordance with manufacturer's instructions. GBM cells were incubated with Seahorse XF base medium supplemented with 1 mmol/L l-glutamine in a CO2-free incubator for 1 hour before the assay. During the course of the assay, 10 mmol/L glucose, 1 μmol/L oligomycin, and 50 mmol/L 2-DG were added sequentially.
LC/MS analysis of metabolites
Metabolite analysis was carried out on a Thermo Scientific QExactive Orbitrap in a manner as described earlier by others (20).
Western blot analysis and capillary electrophoresis on Wes instrument (Proteinsimple)
Specific protein expression in cell lines was determined by Western blot analysis or capillary electrophoresis as described before. Capillary electrophoresis was run on the Wes instrument (Proteinsimple). The following antibodies were used on the Wes instrument: p-Akt (serine 473; 1:25; CST, Cell Signaling Technology), Akt (1:50; CST), Mcl-1 (1:50; CST:), Bcl-2 (1:25; R&D Systems), BIM (1:25; CST), Bcl-xL (1:25; CST), c-myc (1:25; CST), Usp9X (1:25; CST), Noxa (1:25, clone 114C307; Calbiochem), p-Akt (1:25; CST), Akt (1:25; CST), p-AMPK (1:25; CST), AMPK (1:25; CST), PHGDH antibody (Novus; #NBP1-87311), PSAT1 Polyclonal Antibody (Invitrogen; #PA5-22124), β-actin (1:250, clone AC15; Sigma Aldrich), and secondary horseradish peroxidase–linked antibodies were purchased from Santa Cruz Biotechnology Inc. For the expression levels of respiratory complexes, the Total OXPHOS Human WB Antibody Cocktail was used (Abcam). Western blots were acquired, using the Azure (C300) imaging system (CCD – camera based).
Real-time PCR analysis
RNA was isolated and reverse-transcribed as previously described (21). For cDNA amplification, c-myc primers were used: Forward: CCT GGT GCT CCA TGA GGA GAC; Reverse: CAG ACT CTG ACC TTT TGC GAG G. Amplification of 18S served as normalization control. For the determination of mtDNA, the following primers were used: Forward cga aag gac aag aga aat aag g; Reverse: ctg taa agt ttt aag ttt tat gcg. The other primers were designed by Origene Technologies.
Microarray and gene set enrichment analysis
Transcriptome analysis and gene set enrichment analysis (GSEA), involving microarrays, were performed as previously described in ref. 21. The related data and cel files are archived through GEO under the following accession numbers: GSE104273 and GSE103963.
Transfections of siRNAs or transductions of shRNAs
Transfections were performed as previously described (22), using either Oligofectamine or Lipofectamine 2000. CMYC siRNA 1 and 2 were purchased from Cell Signaling Technology. DRD2-specific siRNAs were obtained from Dharmacon. Nontargeting siRNA-pool (ON-TARGETplus Non-targeting Pool, # D-001810-10-20) was purchased from Dharmacon. Lentiviral shRNA particles targeting c-myc were purchased from Santa Cruz Biotechnology.
Orthotopic xenograft models of glioblastoma
The individual cells were injected as previously described in refs. 16, 21. Briefly, 50,000 U87 cells or 3 × 105 GBM123 cells in 2 μL of PBS were intracranially injected into nu/nu mouse. After 3 days of engraftment, the drug treatments were initiated. All treatments were administered by i.p. injection.
Subcutaneous xenograft model
NCH644 (GBM stem-like cells), U87-EGFRvIII, or HCT116 (colon carcinoma) cells were suspended 1:1 in Matrigel Matrix (Corning Inc.) and implanted subcutaneously into the flanks of 6- to 8-week-old Nu/Nu mice. Tumors were measured with a caliper and sizes calculated according to the standard formula: (length × width2) × 0.5. All treatments were administered by i.p. injection. Drug compounds were dissolved in 10% DMSO, 32% Cremophor EL (SIGMA), 8% Ethanol (Pharmco-Aaper), and 50% PBS.
Immunohistochemistry
Statistical analysis
Statistical significance was assessed by two-tailed Student t test or the Mann–Whitney test using Prism version 5.04 (GraphPad). A P ≤ 0.05 was considered statistically significant.
Study approval
All procedures were in accordance with Animal Welfare Regulations and approved by the Institutional Animal Care and Use Committee at the Columbia University Medical Center.
Results
High levels of c-myc protein predict apoptosis induction and inhibition of proliferation by imipridones in model systems of GBM
The chemical structures of the different imipridones are shown in Fig. 1A. Therein, ONC201 represents the lead compound, whereas ONC206 and ONC212 are chemically modified derivatives with anticipated enhanced efficacy (Fig. 1A). To get an understanding of the underlying cause of the observed reduction in cellular viability, we conducted cell-cycle analysis of GBM cultures after treatment with various imipridones for 72 hours. Our observations suggested that GBM cells with high c-myc levels were more likely to respond with apoptosis induction upon imipridone (ONC201, ONC206, and ONC212) treatment (Fig. 1B–G; Supplementary Fig. S1F). While in PDX GBM14 cells, established LN229, T98G and U87 GBM (lower levels) cells we only observed a settle increase in sub-G1 (apoptotic) cells (Fig. 1B and C), this effect was more pronounced in stem-like NCH644, NCH690, and NCH421k cells as well as in a pediatric GBM cell line, SF188 (all very high c-myc expression; Supplementary Fig. S1; Fig. 1B and E–G). Consistent with the highest expression levels of c-myc, NCH421k glioma stem-like cells underwent apoptosis most strikingly upon administration of imipridones (Fig. 1B, E, and G; Supplementary Fig. S1C). Similarly, NCH690 glioma stem-like cells treated with various imipridones displayed a high fraction of Annexin V–positive cells (Supplementary Fig. S1B and S1D), which occurred in a time-dependent manner (Supplementary Fig. S1A–S1C). Notably, c-myc suppression precedes apoptosis induction in NHC421k cells (Supplementary Fig. S1A–S1C). Among imipridones, ONC212 (nanomolar range) was most potent to elucidate these changes when compared with ONC201 (micromolar range; Fig. 1B and G; Supplementary Fig. S1). However, despite the fact that intermediate c-myc–expressing cells (U87) or low c-myc–expressing cells (GBM14) did not undergo apoptosis, they still display inhibition of proliferation (IC80 values) upon imipridone administration in accordance with their expression level of c-myc (Fig. 1E and G). Notably, nonneoplastic human astrocytes displayed the highest resistance against imipridones (Fig. 1G).
We modulated the expression levels of c-myc in an intermediate expressing c-myc line, U87, as well as in MYC-amplified SF188 cells (very high levels). We confirmed that SF188 cells are acutely dependent on c-myc for their survival. To test this hypothesis, we acutely silenced the expression of c-myc by two c-myc targeting siRNAs and subsequently analyzed them for apoptosis induction (Supplementary Fig. S1E). Our findings indicate that silencing of c-myc elicited significant apoptosis induction in these cells, reinforcing their dependence on the transcription factor, c-myc. Next, we created a stable cell c-myc–overexpressing U87 line and evaluated its susceptibility to ONC201 in the context of apoptosis induction (Fig. 1H and I). In agreement with our hypothesis, we found that c-myc–overexpressing cells displayed a higher rate of DNA fragmentation upon imipridone treatment as compared with control-transduced cells. In a reciprocal approach, we created a stable SF188 line with silenced c-myc expression (lentiviral shRNA) that in turn becomes less addicted to c-myc. As anticipated, SF188 GBM cells with silenced c-myc expression were less sensitive to apoptosis induction by imipridones as compared with cells transduced with a nontargeting shRNA (Fig. 1J and K). These findings suggest that imipridones mediate their effects in part by inhibition of proliferation and apoptosis induction in a cell-type and c-myc–dependent manner.
Imipridones regulate phosphorylation of AKT, ERK, and FOXO3A, mTORC1 signaling, expression of anti- and proapoptotic Bcl-2 family members, and c-myc protein levels in GBM model systems
To assess the pharmacodynamics properties of imipridones in model systems of GBM, we performed protein capillary electrophoresis for the expression of phosphorylated Akt, total Akt, phosphorylated ERK, total ERK, phosphorylated FOXO3A (downstream targets of AKT and ERK), and total FOXO3A in U87 and stem-like GBM cell, NCH644, after treatment with ONC201, ONC206, and ONC212 at 200 nmol/L for 72 hours (Fig. 2A; Supplementary Fig. S2A). Although, as anticipated, ONC201 did not show activity, ONC212 was most potent to affect the phosphorylation status of FOXO3A, Akt, and ERK, in keeping with the cellular viability assays. However, higher concentrations of ONC201 resembled the above findings related to ONC212. Next, we assessed the activation status of mTORC1 signaling after imipridone exposure to U87 and stem-like GBM cell, NCH644, because earlier findings suggested that ONC201/TIC10 interfered with this signaling cascade. We found that mTORC1 targets, mTOR, 4EBP1, S6K, and S6, were most potently and consistently inhibited by ONC212 in U87 and NCH644 stem-like GBM cells (Fig. 2A; Supplementary Fig. S2A).
To determine if the potency of imipridones correlates with the regulation of apoptosis mediators (Bcl-2 family proteins; ref. 13), including their chaperones, Usp9X and Bag3, we assessed as to whether or not ONC201 derivatives, ONC206 and ONC212, would regulate the protein levels of those molecules in a more potent manner. We found that in U87 and stem-like GBM cells, NCH644, ONC206 and ONC212 at 200 nmol/L suppressed the expression levels of antiapoptotic, Bcl-2, Bcl-xL and Mcl-1, and the deubiquitinase Usp9X, whereas in contrast at the same concentration, ONC201 was less effective (Fig. 2B; Supplementary Fig. S2B). Under these conditions, proapoptotic Bcl-2 family members, BIM and Noxa, were downregulated after ONC206 and ONC212 treatment (Fig. 2B; Supplementary Fig. S2B). However, ONC201 (10 μmol/L) upregulated BIM in U87 cells, priming cells for potential proapoptotic stimuli (Fig. 2B; Supplementary Fig. S2B). In essence, these findings support the idea that imipridones convert the cells into a more “apoptosis favorable” state.
Imipridones downregulate protein levels of the transcription factor, c-myc, at the posttranslational level
We found substantial evidence for suppression of c-myc downstream targets by imipridones in accordance with their potency (Fig. 2C). We conducted a time course experiment to assess c-myc protein levels and detected a downregulation of c-myc as early as 2 hours and extending to 72 hours after imipridone treatment in established and stem-like GBM cells (Fig. 2D and H). Of note, ONC212 had the strongest impact on c-myc levels, consistent with its strongest potency among imipridones. We hypothesized that based on our transcriptome analysis, most likely imipridones modulate the stability of c-myc protein because this transcription factor is known to have a very short half-life and is highly susceptible to proteasomal degradation upon phosphorylation at threonine 58. Therefore, we conducted a time course analysis of phosphorylated and total c-myc protein levels upon stimulation with imipridones (Fig. 2H). Although ONC201 led to a rapid depletion of total c-myc protein (within 2 hours) and decrease of phosphorylated c-myc at serine 62, it increased phosphorylation of c-myc at threonine 58, which was preceded and coincided with a reduction in phosphorylation of GSK3β at serine 9, a kinase that is known to phosphorylate c-myc at threonine 58 (Fig. 2H; refs. 23, 24). Given that GSk3β is a substrate of AKT, we assessed AKT phosphorylation status and found that ONC201 decreased phosphorylation of AKT at serine 473, coinciding with GSK3β dephosphorylation (Fig. 2H). Consistently, two pharmacologic inhibitors of GSK3β restored c-myc protein levels upon imipridone treatment (Fig. 2I), confirming the pivotal role of GSK3β in imipridone-mediated reduction of c-myc protein levels. Moreover, imipridones did not suppress protein levels of a T58A c-myc mutant (a hotspot region for mutations), further supporting the role of GSK3β in the mechanism (Fig. 2J). To confirm that these phosphorylation events of c-myc result in a decrease of c-myc stability, U87 GBM cells were treated with ONC212 or DMSO and exposed to protein synthesis inhibitor, cycloheximide. ONC212 reduced the stability of c-myc protein (Fig. 2E and 2F). Next, we treated U87 GBM cells with ONC212 in the presence or absence of the bona fide proteasomal inhibitor, MG132. In keeping with the initial hypothesis, MG132 attenuated ONC212-mediated suppression of c-myc protein (Fig. 2G). Similar findings were made with the lead compound ONC201 (Supplementary Fig. S2D and S2E). To determine the relative contribution of c-myc mRNA suppression by imipridones, we treated U87GBM cells with various concentration of ONC201, ONC206, and ONC212 and assessed c-myc mRNA levels. We found only slight suppression of c-myc mRNA levels, suggesting that the major effect of imipridones on c-myc is nontranscriptional (Supplementary Fig. S3A and S3B). Overall, these findings suggest that imipridone enhances c-myc degradation, most likely through regulation of the proteasomal pathway, involving GSK3β-mediated phosphorylation of threonine 58.
Next, we posed the question as to whether or not simultaneous targeting of c-myc along with imipridones provides a potential therapeutic benefit. To this end, we utilized a bromodomain protein inhibitor (BRD2-4), OTX015, a compound that has been shown to transcriptionally antagonize the expression of c-myc (25, 26), and tested the efficacy of imipridones in the presence or absence of this inhibitor in stem-like GBM cells, NCH644 and NCH421k. NCH644 and NCH421k cells were treated with OTX015, imipridones, or the combination and analyzed for apoptosis by Annexin V/PI staining with subsequent flow cytometry. In the presence of OTX015, imipridone-induced apoptosis was significantly enhanced (Supplementary Fig. S2C), suggesting that dual-targeting c-myc by BRD4 inhibitors and imipridones might be viable treatment strategy that is effective against tumor-initiating cells.
Imipridones reprogram the transcriptome of GBM cells and suppress glycolysis and oxidative phosphorylation
In order to identify the underlying mechanisms by which imipridones elicit their antiglioma effects, we conducted a transcriptome analysis and subsequent GSEA of GBM cells treated with ONC201, ONC206, and ONC212 (at 200 nmol/L) as well as ONC201 at 10 μmol/L. In keeping with the Western blot findings, we only identified transcriptomic changes in cells treated with ONC206 and ONC212 at 200 nmol/L, whereas ONC201 appeared to elicit no effects at that concentration (not shown). However, higher concentrations of ONC201 (10 μmol/L) resembled the effects of ONC206 and ONC212. Our GSEA uncovered additional mechanistic insights into the biology of imipridones. Most notably, we found changes related to cellular metabolism that suggest that imipridones suppress glycolysis as well as oxidative phosphorylation, resulting in a “transcriptional state of energy starvation” (Fig. 3A). These findings were linked to changes in the expression of relevant associated transcription factors. As mentioned earlier, we identified several gene sets to be downregulated related to c-myc targets (Supplementary Fig. S3A), suggesting that first c-myc is likely a key player in imipridone-mediated antiglioma activity and that second c-myc is likely responsible for the transcriptional metabolic reprogramming of GBM cells by imipridones because many of the c-myc targets are genes implicated in metabolism.
In light of these transcriptional findings, we hypothesized that imipridones affect ATP levels of tumor cells. Our results suggest that ONC201 and ONC206 reduce ATP levels in a dose-dependent manner (Fig. 3B). Notably, ONC212 was most efficacious to reduce ATP levels, followed by ONC206 and ONC201 (Fig. 3B). Accompanied by the decline of ATP levels, we found an increase of phosphorylation of the AMP-sensitive kinase, AMPK (threonine 172), starting as early as 24 hours (Supplementary Fig. S3C). In keeping with the ATP levels, ONC212 was most potent in enhancing the phosphorylation of AMPK. It is generally accepted that low energy levels are causal for the accumulation of unfolded proteins in the endoplasmic reticulum. In keeping with that notion, we found evidence of activation of several ER stress–related pathways as manifested by elevation of IRE1α, ATF4, ATF3, and phosphorylated eif2α proteins (Supplementary Fig. S3D).
Because the global transcriptome analysis suggests an inhibition of glycolysis and oxidative energy metabolism (Supplementary Fig. S4A–S4E), we conducted extracellular flux analysis. To this purpose, SF188 (high c-myc) and U87 (intermediate c-myc levels) GBM cells were treated with vehicle, ONC201, ONC206, or ONC212 and assayed for glycolytic activity by extracellular flux analysis (Fig. 3C–E; Supplementary Figs. S4E–S4L and S5A–S5E). We found that imipridones, ONC206 and ONC212, were most potent to suppress baseline extracellular acidification rate (ECAR; Supplementary Fig. S5A–S5E), which was more pronounced in high c-myc–expressing SF188 cells as compared with U87 (Fig. 3E). However, utilizing 10 μmol/L of the lead compound, ONC201, recapitulated the effects of ONC206 and ONC212 on glycolysis (Supplementary Fig. S4F–S4I). In addition, we tested whether or not c-myc–overexpressing U87 cells are more susceptible to glycolysis inhibition by imipridones. As shown in Fig. 3F, c-myc–overexpressing cells are more susceptible to ECAR reduction as compared with vector-transduced cells. We assessed the protein levels of glycolysis-related enzymes and transporters, HK2, LDHA, and GLUT1. To this purpose, U87 cells were treated with various imipridones as above. We found that ONC206 and ONC212 reduced the levels of HK2, GLUT1, and LDH (Fig. 3G). In order to confirm that imipridone-mediated suppression of c-myc is sufficient to affect glycolysis, we performed knockdown experiments, involving two c-myc siRNA oligonucleotides. There was an siRNA-dependent reduction in c-myc protein levels and an associated depletion of HK2, LDHA, and GLUT1, in keeping with the hypothesis that c-myc governs the expression of these proteins and that this is also true in the setting of GBM model systems (Fig. 3). To further delineate the role of c-myc in glycolytic metabolism, we pharmacologically inhibited c-myc by using OTX015. We have previously shown that the concentrations of OTX015 applied in the current studies were sufficient to suppress c-myc protein levels (10, 27). Using this dosage of OTX015, we performed extracellular flux analysis. We found that c-myc inhibition recapitulates the effects on ECAR-mediated imipridones, furthermore suggesting the pivotal role of c-myc in imipridone-mediated glycolytic inhibition (Fig. 3I). To further assess the impact of imipridones on glycolysis, we measured glycolytic metabolites by LC/MS following ONC212 treatment. These results show that ONC212 leads to an accumulation of glycolytic metabolites, suggesting a glycolytic block with an overall reduction of glycolysis (as supported by the extracellular flux analysis with reduction in ECAR; Fig. 3J).
Next, we determined the effects of imipridones on OXPHOS-related metabolism. To this end, U87 GBM cells were treated with vehicle, ONC201, ONC206, and ONC212 for 24 hours and analyzed for basal oxygen consumption rate (OCR), ATP production, and maximal respiration (Fig. 4A–C; Supplementary Fig. S5F). We detected a significant reduction in all these parameters in cells treated with 200 nmol/L ONC206 and ONC212 (Fig. 4A–C; Supplementary Fig. S5F). In contrast, 200 nmol/L of ONC201 had little effect on mitochondrial respiration (Fig. 4A and B). However, 10 μmol/L of ONC201 recapitulated the effects of ONC206 on ONC212 (Fig. 4A–C; Supplementary Fig. S5H), further corroborating the earlier notion that ONC206 and ONC212 are more potent than ONC201. We assessed the relative potency of imipridone OCR reduction in the context of high and intermediate c-myc expression. As anticipated, high c-myc–expressing SF188 cells were more prone to OCR inhibition by imipridones as compared with U87 cells (Fig. 4C).
Given the significant suppression of OXPHOS, we evaluated the underlying mechanisms governing these changes. We conducted Western blot analysis to detect the expression levels for respiratory complexes, pyruvate dehydrogenase (PDH). To this purpose, U87, LN229, and stem-like GBM cells, NCH644, were treated with vehicle, ONC201, ONC206, and ONC212 for 72 hours (Fig. 4D; Supplementary Fig. S5H and S5I). Although no effects on the expression of respiratory complexes were detected with 200 nmol/L of ONC201, the imipridone derivatives, ONC206 and ONC212, suppressed the levels of OXPHOS complexes as well as of PDH at that concentration (Fig. 4D). However, 10 μmol/L of ONC201 suppressed OXPHOS complexes as well (Fig. 4D; Supplementary Fig. S5H and S5I). Similarly, c-myc knockdown also affected the protein expression of respiratory complexes (Fig. 4E and F), in keeping with the effects of imipridones. To assess the functional implication of c-myc on OXPHOS in glioblastoma cells, U87 GBM cells were treated with OTX015 for 24 hours and subsequently analyzed by extracellular flux analysis. OTX-015 reduced the OCR (Fig. 4G). Given the suppression of OXPHOS, we analyzed TCA metabolite levels upon treatment with imipridone derivatives by LC/MS. In keeping with the findings obtained with the extracellular flux analysis, we encountered a suppression of TCA cycle intermediates (Fig. 4H), further corroborating the earlier notion that imipridones significantly impair tumor energy metabolism. Consistent with an impact on mitochondrial metabolism, we found a significant suppression of mtDNA levels in imipridone-treated U87 GBM cells, which is most pronounced in ONC212-treated cells (Fig. 4I). Similarly, mitochondrial membrane potential was lowered by ONC212 (Supplementary Fig. S5G).
Imipridones enhance serine-one carbon-glycine metabolism
Although imipridones suppress the main energetic metabolic pathways of GBM cells, we were hypothesizing that compensatory pathways will be activated. In this vein, GSEA provided information that imipridone-treated GBM cells activate the serine-one carbon-glycine metabolism with an increase of PHGDH, PSAT1, and PSPH and other enzyme-related genes to the folate cycle (Fig. 5A; Supplementary Fig. S6A). These findings were confirmed by real-time PCR and protein expression analysis in both glioblastoma and colonic carcinoma cells, suggesting that these observations are not limited to just GBM (Fig. 5B and D; Supplementary Fig. S7A). Consistently, we found that depriving GBM cells of energy (ATP) by the ATP synthase inhibitor, oligomycin, upregulates PHGDH and PSAT1, in keeping with our hypothesis that energy deprivation leads to the activation of this pathway (Supplementary Fig. S6D). To corroborate these findings, LC/MS analysis was performed, which showed significant increases of metabolites related to the serine-one carbon-glycine metabolism in imipridone-treated cells (Fig. 5C). To address the question by which mechanism imipridones induce this pathway, we consulted our GSEA data and the literature and concluded that likely ATF4 is the pivotal mediator of these transcriptional changes. To test this hypothesis, we silenced ATF4 levels in the presence or absence of imipridones and found that ATF4 silencing (involving one siRNA pool, four distinct siRNAs, and a shRNA construct) attenuates imipridone-mediated increase of PHGDH and PSAT1 (Fig. 5E; Supplementary Figs. S6B and S6C, and S7B and S7C). Next, we determined the translational impact of pharmacologic targeting of the SOG pathway along with imipridones. To this end, we used two recently identified pharmacologic inhibitors of PHGDH, NCT-503, and CBR-5884 (28), and combined them with ONC212. We found that the combination treatment of NCT-503 and ONC212 enhanced induction of apoptosis more potent than each compound alone in stem-cell like GBM cells, NCH421k, NCH644, and NCH690, as well as in U87 GBM cells (Fig. 5F). Similar results were obtained when we administered CBR-5884 in lieu of NCT-503 (Fig. 5G). These synergistic apoptotic effects were not limited to glioblastoma model systems since HCT116 colon carcinoma cells showed a comparable response rate to the combination treatment (Supplementary Fig. S7D). The findings of the activation of the SOG pathway are summarized in Supplementary Fig. S7G.
The novel imipridone derivatives are active either singly or in combination therapies in multiple xenograft models and extend host survival in orthotopic PDXs of GBMs
PDX model systems remain the most relevant preclinical model systems to predict antitumor responses in patients, which is accepted in the studies of GBM as well. Following this strategy, we compared the efficacy of the lead compound, ONC201, with the chemical derivative, ONC212, which was the most potent in all assays related to in vitro testing. To this purpose, we used two PDX models, GBM12 and GBM43, that have different mutational backgrounds and susceptibilities to the standard–of-care treatment (temozolomide; ref. 29). These lines were implanted in immunocompromised mice, and after tumors were established, three treatment groups were formed, consisting of vehicle, ONC201, and ONC212. Treatment was performed 2 times a week. Unanimously, we found that ONC212 was most potent in reducing tumor growth (Supplementary Fig. S8A–S8F), intimating that imipridones are efficacious in relevant preclinical animal models of GBM and that chemical modification results in an improvement of efficacy of imipridones. Finally, to account for all here studied imipridones, we conducted a third animal model, using the highly challenging to treat U87-EGFRvIII model. Tumor cells were implanted subcutaneously, and upon establishment of the tumors, three groups were formed for subsequent treatments with vehicle, ONC206, and ONC212. Upon conclusion of the study, host tumor sizes in animals treated with ONC212 were the smallest as compared with vehicle or ONC206 (Supplementary Fig. S8G–S8I), in keeping with the results in the PDX models. To account for toxicity, we measured the weights of the animals and found no significant weight loss by imipridone derivatives (Supplementary Fig. S8I). With respect to morphologic appearance, we noted that tumors treated with imipridones were smaller, and on the histopathologic level, ONC212 showed the most prominent reduction in cellular density and mitotic rate and an increase in cell death (Supplementary Fig. S9C and S9E). We assessed the effects on imipridones on cellular proliferation in vivo by staining for Ki67. In comparison with vehicle, both ONC201 and ONC212 reduced the proliferation index and number of positive cells (Supplementary Fig. S9C). However, ONC212 was most potent. To assess the amount of apoptosis induced by imipridones in vivo, we stained tumor section for TUNEL. Although vehicle and ONC201 revealed minimal staining, ONC212 displayed significantly more positive tumor cells (Supplementary Fig. S9D). These observations are in keeping with the significantly smaller tumor size in animals that were treated with ONC212. Taken together, these results intimate that imipridone derivatives are potential drug compounds for the treatment of malignant glial brain tumors.
Based on these results, we extended our studies to orthotopic PDXs of glioblastoma with high or low c-myc expression (Supplementary Fig. S8J) in order to test the hypothesis that high levels of c-myc protein predict a better response to imipridone treatment in the most preclinical relevant models. GBM123 cells (high c-myc protein) were more susceptible to imipridones as compared with U87 GBM cells (low-to-intermediate c-myc protein levels; Supplementary Fig. S8J–S8L). For these experiments, we used the most potent imipridone, ONC212, given the well-known challenges of treating orthotopic glioblastoma models. U87 or GBM123 cells were injected into the brain of nude mice. After establishment of tumors, animals were treated with vehicle or ONC212. Although in the U87 model we only detected a small survival benefit, we found a substantial increase of overall survival in the ONC212 treatment group in the GBM123 model (Fig. 6A and B). Notably, when comparing the survival curves of ONC212 in both the U87 and GBM123 model, we found that ONC212 prolonged host survival significantly longer as compared with the U87 model, suggesting that high c-myc protein levels predict a longer host survival upon imipridone treatment (P < 0.05).
Given the role of c-myc in imipridone-mediated cell death, the in vivo implication of high c-myc levels, and our earlier findings that interference with c-myc by BRD4 inhibitors enhances the efficacy of imipridones, we assessed the combination treatment strategy of imipridones and BRD4 inhibitors in a xenograft model, utilizing the NCH644, stem-like GBM cells. Although single treated host animals revealed a slight reduction in tumor growth, the combined treatment of ONC201 along with OTX015 led to a significant reduction in tumor growth with significant cell death/necrosis (Fig. 6C–F). Despite this efficacy, the combination treatment had no effect on the viability of organs, such as brain, heart, lung, liver, intestine, spleen, and kidney, and did not induce significant weight loss (Supplementary Fig. S9A and S9B). Finally, we confirmed that PHGDH inhibitors and imipridones reduce tumor growth synergistically in vivo. To this purpose, we tested the glioblastoma model system, U87-EGFRvIII. Host animals treated with the combination treatment of ONC212 and NCT-503 had significantly smaller tumors than single or vehicle treatments (Fig. 6G). As an additional model system, the HCT116 colon carcinoma xenograft model was tested. In like manner, host animals that received the combination treatment of ONC212 and NCT-503 had significantly smaller tumors than single or vehicle treatments (Fig. 6H–I; Supplementary Fig. S7E and S7F).
Discussion
Glioblastoma remains one of the biggest challenges in medicine, and thus far, researchers have not been successful to provide a substantial improvement in prognosis for this detrimental disease (30). Despite the fact that alike other tumors GBMs have known recurrent genetic alterations and certain point mutations for which an arsenal of inhibitors already exists, none of these compounds has shown a durable response clinically (31, 32). The reason for these failures might be multiple, but certainly related to the fact that GBMs have a substantial, intrinsic heterogeneity that in itself requires either pleiotropic drug compounds that interfere with multiple pathways or a rational drug combination therapy at the onset (27, 33–38).
Imipridone derivatives elegantly fit into this category because these compounds act on several pivotal pathways for gliomagenesis, maintenance, and progression of these neoplasms, e.g., AKT and ERK signaling and the antiapoptotic Bcl-2 family members of proteins (13). The first compound out of this category is TIC10/ONC201 that was discovered in an attempt to identify compounds that increase expression for the death ligand TRAIL (3). Shortly thereafter, it was suggested that ONC201 elicits the induction of an integrated stress response, involving the transcription factor, ATF4 (4), leading to an inhibition of proliferation of cancer cells. However, the underlying mechanism as to how this happens remained elusive, and our current study provides significant additional novel insight about the underlying events.
Having discovered a potent class of molecules with anticancer properties, it was tempting to develop derivatives that are more potent than the lead compound. Examples of this strategy are the imipridones, ONC206 and ONC212 (2). In our studies, we found that ONC212 was the most potent inhibitor with respect to antiglioma activity, followed by ONC206 and ONC201. These results are in keeping with one earlier study that focused on pancreatic adenocarcinoma model systems (2).
In seeking more upstream mechanisms as to how imipridones elicit their antiglioma activity, we conducted transcriptome analysis with subsequent GSEA and metabolite analysis by LC/MS. We made the novel observation that imipridones affected energy metabolism by suppression of genes encoding for enzymes and metabolites related to glycolysis and oxidative phosphorylation. Consistently, the proteins related to those genes were downregulated as well. To the best of our knowledge, these findings have not been reported thus far and mechanistically significantly extend the understanding of how imipridones elicit their antitumor effects. Given that many tumors are “glycolytic” in nature, it is a notable finding that imipridones interfere with both OXPHOS and glycolysis because many compounds that target only one of these pathways, such as 2-DG, will result in the compensatory activation of other major metabolic energy pathways. For instance, 2-DG inhibits glycolysis, and in turn, tumor cells activate oxidative phosphorylation to balance their energy needs. Mechanistically, this is owed in part due to mitochondrial plasticity as 2-DG inhibits mitochondrial fission, leading to enhanced fusion by modulation of the phosphorylation status of Drp1 (39). Our findings also show that imipridones regulate mitochondrial plasticity because they lead to enhanced mitochondrial fission, which in turn dampens OXPHOS activity. Another example is gemcitabine resistance in model systems of pancreatic cancer because these tumors enhance glucose metabolism to escape from therapy (40). In this context, it may be tempting to speculate as to whether or not imipridones might reverse gemcitabine resistance or act synergistically to reduce tumor cell proliferation with this chemotherapeutic drug.
To compensate for this massive impact on tumor cell energy metabolism by imipridones, we found that these compounds activate the serine-one carbon-glycine pathway that is utilized by several types of tumors to support growth. Interestingly, the literature also suggests a role for this pathway to synthesize ATP independent of OXPHOS and glycolysis (41). Therefore, imipridone might utilize this pathway in part to produce energy. In keeping with its prosurvival function (42–45), interference with this pathway through PHGDH inhibition (28), using two novel pharmaceutical inhibitors, leads to synergistic cell death in vitro and in vivo. Although we have not evaluated methotrexate in combination with imipridone, it is likely that this compound might synergize with imipridones as well because methotrexate interferes with the folate cycle and thus with the serine-glycine pathway.
We determined the mechanism how imipridones suspend energy metabolism in tumor cells and found that the transcription factor c-myc (46–48) is involved in this process. Although c-myc was not downregulated at the level of transcription, we found a marked decrease of c-myc on the protein level as early as 1 hour after treatment. Consistently, cells with high levels of c-myc are more susceptible to imipridone treatment, an important finding that has not been shown previously before. Given that GBM stem-like cells are known to be highly dependent on c-myc, our observation that imipridones induce massive apoptosis in these cells is in keeping with this earlier notion. Therefore, our findings suggest that imipridones have the potential to eradicate the stem-cell fraction within a glial neoplasm, thus counteracting the development of treatment resistance and recurrence. Similarly, the pediatric GBM cell line SF188 (c-myc amplified) displays high sensitivity to imipridone derivative, suggesting that c-myc status may predict apoptotic susceptibility to imipridone derivatives. However, it should be noted that even tumor cells with lower expression levels of c-myc display an antiproliferative response to imipridones, albeit without significant apoptosis induction. Our findings also demonstrate that inhibition of c-myc superenhancers by BRD4 inhibitors along with imipridones synergistically reduces the viability of glioma cells in vitro and in vivo, which has not been show in other tumor entities before. Therefore, it might be tempting to speculate whether or not nonsolid malignancies might respond to this novel combination therapy as well.
Finally, we confirmed whether or not imipridones prolong host survival in orthotopic xenograft models. Our findings showed that ONC212 extended host survival in a high c-myc–expressing PDX model, whereas in a low c-myc–expressing model system, these antiglioma effects were significantly less pronounced. These findings support the overall notion that high c-myc levels predict therapeutic responses to imipridones and that chemical modification of imipridones potently enhances antitumor activity.
Taken together, these findings provide a foundation for clinical testing of these compounds either singly or in combination regimes in patients. In addition, utilization of c-myc as a biomarker might be a strategy for prediction of therapeutic responses.
Disclosure of Potential Conflicts of Interest
V.V. Prabhu is an employee of and holds ownership interest (including patents) in Oncoceutics. J.E. Allen holds ownership interest (including patents) in Oncoceutics. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: C.T. Ishida, Y. Zhang, V.V. Prabhu, J.E. Allen, M.D. Siegelin
Development of methodology: C.T. Ishida, G. Karpel-Massler, V.V. Prabhu, M.D. Siegelin
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): C.T. Ishida, Y. Zhang, C. Shu, G. Kleiner, M.J. Sanchez-Quintero, M.D. Siegelin
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): C.T. Ishida, Y. Zhang, M.D. Siegelin
Writing, review, and/or revision of the manuscript: C.T. Ishida, Y. Zhang, E. Bianchetti, C.M. Quinzii, M.-A. Westhoff, G. Karpel-Massler, V.V. Prabhu, J.E. Allen, M.D. Siegelin
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): E. Bianchetti, C. Shu, T.T.T. Nguyen, V.V. Prabhu, M.D. Siegelin
Study supervision: C.M. Quinzii, M.D. Siegelin
Acknowledgments
This work was supported by the NIH NINDS K08NS083732, R01NS095848, R01NS102366, the 2017 American Brain Tumor Association Discovery Grant (DG1700013), and Louis V. Gerstner, Jr. Scholars Program (2017-2020). The flow cytometry experiments were performed in the CCTI Flow Cytometry Core, supported in part by the Office of the Director, NIH under award S10RR027050. Transcriptome analysis was supported by the CTSA grant UL1-TR001430 to the Boston University Microarray and Sequencing Resource Core Facility.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.