Metabolic changes induced by oncogenic drivers of cancer contribute to tumor growth and are attractive targets for cancer treatment. Here, we found that increased growth of PTEN-mutant cells was dependent on glutamine flux through the de novo pyrimidine synthesis pathway, which created sensitivity to the inhibition of dihydroorotate dehydrogenase, a rate-limiting enzyme for pyrimidine ring synthesis. S-phase PTEN-mutant cells showed increased numbers of replication forks, and inhibitors of dihydroorotate dehydrogenase led to chromosome breaks and cell death due to inadequate ATR activation and DNA damage at replication forks. Our findings indicate that enhanced glutamine flux generates vulnerability to dihydroorotate dehydrogenase inhibition, which then causes synthetic lethality in PTEN-deficient cells due to inherent defects in ATR activation. Inhibition of dihydroorotate dehydrogenase could thus be a promising therapy for patients with PTEN-mutant cancers.

Significance: We have found a prospective targeted therapy for PTEN-deficient tumors, with efficacy in vitro and in vivo in tumors derived from different tissues. This is based upon the changes in glutamine metabolism, DNA replication, and DNA damage response which are consequences of inactivation of PTEN. Cancer Discov; 7(4); 380–90. ©2017 AACR.

See related article by Brown et al., p. 391.

This article is highlighted in the In This Issue feature, p. 339

The Warburg effect, describing heightened aerobic glycolysis in tumors, played a key role in launching the field of cancer metabolism. Subsequent studies have found that glutamine is also vital for growth by fueling tricarboxylic acid cycle intermediates, phospholipid and nucleotide synthesis, and NADPH (1). Oncogenic signaling pathways have been shown to play a major role in reprogramming glucose and glutamine metabolism, thus connecting genetic mutations with metabolic alterations (2–5). PTEN is one of the most commonly mutated tumor suppressors and is a fulcrum of multiple cellular functions (6, 7). PTEN's canonical role is as a lipid phosphatase for phosphatidylinositol-3,4,5-trisphosphate, central to the PI3K pathway, limiting AKT, mTOR, and RAC signaling (8–11). Inactivation of PTEN enhances glucose metabolism and diminishes DNA repair and DNA damage checkpoint pathways (12–14). Furthermore, deficient homologous recombination in PTEN-mutant cells leads to sensitivity to γ-irradiation and PARP inhibitors (13, 15). The role of PTEN in metabolism, however, is incompletely understood, and in this study we examine the metabolic consequences of PTEN loss and the resulting vulnerability of PTEN-mutant tumors.

To better understand the relationship between PTEN, cell growth, and cellular metabolism, we generated Pten flox/flox primary mouse embryonic fibroblasts (MEF). Pten−/− MEFs proliferated at a higher rate than wild-type (WT) MEFs but showed no difference in cell death (Fig. 1A; Supplementary Fig. S1A–S1C). This increased proliferation was associated with an increase in the proportion of cells within S-phase and higher numbers of replication forks per S-phase cell (Fig. 1B–D; Supplementary Fig. S1D). There was no difference in mitochondrial function between Pten−/− and WT MEFs, suggesting a different source of altered replication (Supplementary Fig. S1E and S1F).

Figure 1.

A, Growth of Pten WT and knockout (KO) MEFs (one-way ANOVA, *, P < 0.0001, n = 3). B, MEFs labeled with EdU. Representative confocal microscopy images. C, Quantification of B (Student t test, *, P < 0.05, n = 6). D, MEFs labeled with EdU; flow cytometry determined the mean fluorescence intensity (MFI) among cells positively stained (Student t test, *, P < 0.01, n = 3). E,Pten WT and KO MEFs in media containing full glutamine (6 mmol/L) or no added glutamine (one-way ANOVA, *, P < 0.0001, n = 3). F, MEFs treated with 12.5 nmol/L CB-839 or control (one-way ANOVA, *, P < 0.0001, n = 3). G, Relative metabolite concentrations of DNA nucleotide precursors (dGMP was unable to be measured so dGTP was used; Student t test, *, P < 0.05, n = 3). H, Relative metabolite levels of glutamine-labeled de novo pyrimidine synthesis intermediates (Student t test, *, P < 0.05, n = 3). Data were also analyzed with IMPaLA: 13C glutamine-derived pyrimidine metabolism enrichment in PTEN−/− MEFs q value = 3.92 × 10−09. I, Schematic of the de novo pyrimidine synthesis pathway. Not every intermediate was measured in our mass spectrometry panel. Data, means ± SD. TCA, tricarboxylic acid.

Figure 1.

A, Growth of Pten WT and knockout (KO) MEFs (one-way ANOVA, *, P < 0.0001, n = 3). B, MEFs labeled with EdU. Representative confocal microscopy images. C, Quantification of B (Student t test, *, P < 0.05, n = 6). D, MEFs labeled with EdU; flow cytometry determined the mean fluorescence intensity (MFI) among cells positively stained (Student t test, *, P < 0.01, n = 3). E,Pten WT and KO MEFs in media containing full glutamine (6 mmol/L) or no added glutamine (one-way ANOVA, *, P < 0.0001, n = 3). F, MEFs treated with 12.5 nmol/L CB-839 or control (one-way ANOVA, *, P < 0.0001, n = 3). G, Relative metabolite concentrations of DNA nucleotide precursors (dGMP was unable to be measured so dGTP was used; Student t test, *, P < 0.05, n = 3). H, Relative metabolite levels of glutamine-labeled de novo pyrimidine synthesis intermediates (Student t test, *, P < 0.05, n = 3). Data were also analyzed with IMPaLA: 13C glutamine-derived pyrimidine metabolism enrichment in PTEN−/− MEFs q value = 3.92 × 10−09. I, Schematic of the de novo pyrimidine synthesis pathway. Not every intermediate was measured in our mass spectrometry panel. Data, means ± SD. TCA, tricarboxylic acid.

Close modal

Upon testing the potential role of glutamine for explaining the increased growth of Pten−/− cells, we found that the growth advantage of Pten−/− MEFs was dependent on glutamine: Depletion of glutamine or the addition of the glutaminase inhibitor CB-839 collapsed the growth difference between Pten−/− and WT MEFs (Fig. 1E and F). To better understand the relationship between PTEN and glutamine, we performed steady-state metabolomic profiling of 292 aqueous metabolites to determine whether the loss of PTEN triggers abnormal cellular metabolism to increase growth. Unbiased global metabolic assessment of WT and Pten−/− MEFs revealed that seven of the ten most upregulated pathways in Pten−/− MEFs involved nucleotide synthesis and DNA metabolism, including a higher concentration of pyrimidine 2-deoxyribonucleotides in Pten−/− MEFs (Fig. 1G; Supplementary Fig. S2A and S2B). Because glutamine contributes both nitrogen and carbon to pyrimidines (16), we performed metabolic flux analysis with heavy isotope 15N or 13C-labeled glutamine, which showed increased synthesis of dihydroorotate, orotate, and other components of the de novo pyrimidine synthesis pathway in Pten−/− MEFs relative to WT (Fig. 1H and I; Supplementary Fig. S2C). In addition, the pyrimidine metabolism gene set was upregulated in mRNA from Pten−/− MEFs (Supplementary Fig. S2D). Although Pten−/− fibroblasts had somewhat elevated steady-state glucose metabolism and glycolytic flux relative to WT, depletion of glucose from the medium did not rescue the differences in cell growth, suggesting that glutamine was more critical for the growth advantage of Pten−/− cells (Supplementary Fig. S2A, S2E, and S2F). Nucleotide synthesis is a prerequisite for cellular growth, and Pten−/− MEFs appear to channel glutamine for this purpose.

The fourth step of de novo pyrimidine synthesis in mammals is the conversion of dihydroorotate to orotate, catalyzed by dihydroorotate dehydrogenase (DHODH; ref. 17). To see whether orotate contributes to the growth effects observed, the effect of DHODH inhibitors on cell proliferation was examined. Pten−/− MEFs were about 3-fold more sensitive to leflunomide, a DHODH inhibitor (18), than WT MEFs were (Fig. 2A; Supplementary Fig. S3A and S3B). Pten−/− MEFs were also more sensitive to A771726, the active metabolite of leflunomide (18, 19), as well as another DHODH inhibitor, brequinar (18), indicating that the observed effects were likely through DHODH (Fig. 2A).

Figure 2.

A,Pten WT and knockout (KO) cells treated with dose titrations of leflunomide, A771726, or brequinar to determine GI50s (Student t test, *, P < 0.05, n = 3). B and C, Cells treated with dose titrations of leflunomide to determine GI50s (Student t test, *, P values on figures, n = 3). D and E, Cells treated with 100 μmol/L leflunomide and DRAQ7 to monitor accumulation of cell death, in intervals of 6 hours (one-way ANOVA, *, P values on the figures). F, Human breast cancer cell line growth rates. G, Immunoblots of pAKT in nuclear fractions of Pten−/− and Pik3ca-mutant MEFs. H, Cells treated with 50 μmol/L leflunomide in combination with 0 or 640 μmol/L orotate. Confluence of cells after 5 days of treatment was measured (Student t test, *, P < 0.05, n = 3). I, Cells treated with 50 μmol/L leflunomide in combination with 0, 31.25, 62.5, or 125 μmol/L orotate. Confluence of cells after 5 days was measured (Student t test, *, P < 0.05, n = 3). J, Cells treated with 100 μmol/L leflunomide in combination with 0 or 3.125 mmol/L uridine. Confluence of cells after 5 days of treatment was measured (Student t test, *, P < 0.05, n = 3). K, Cells treated with 100 μmol/L leflunomide in combination with 0, 3.125, or 6.25 mmol/L uridine. Confluence of cells after 5 days was measured (Student t test, *, P < 0.05, n = 3). Data, means ± SD.

Figure 2.

A,Pten WT and knockout (KO) cells treated with dose titrations of leflunomide, A771726, or brequinar to determine GI50s (Student t test, *, P < 0.05, n = 3). B and C, Cells treated with dose titrations of leflunomide to determine GI50s (Student t test, *, P values on figures, n = 3). D and E, Cells treated with 100 μmol/L leflunomide and DRAQ7 to monitor accumulation of cell death, in intervals of 6 hours (one-way ANOVA, *, P values on the figures). F, Human breast cancer cell line growth rates. G, Immunoblots of pAKT in nuclear fractions of Pten−/− and Pik3ca-mutant MEFs. H, Cells treated with 50 μmol/L leflunomide in combination with 0 or 640 μmol/L orotate. Confluence of cells after 5 days of treatment was measured (Student t test, *, P < 0.05, n = 3). I, Cells treated with 50 μmol/L leflunomide in combination with 0, 31.25, 62.5, or 125 μmol/L orotate. Confluence of cells after 5 days was measured (Student t test, *, P < 0.05, n = 3). J, Cells treated with 100 μmol/L leflunomide in combination with 0 or 3.125 mmol/L uridine. Confluence of cells after 5 days of treatment was measured (Student t test, *, P < 0.05, n = 3). K, Cells treated with 100 μmol/L leflunomide in combination with 0, 3.125, or 6.25 mmol/L uridine. Confluence of cells after 5 days was measured (Student t test, *, P < 0.05, n = 3). Data, means ± SD.

Close modal

To determine whether PTEN genotype is predictive of sensitivity to DHODH inhibition in cancer cells, we tested human breast cancer, glioblastoma, and prostate cancer cell lines with DHODH inhibitors. Consistently, the GI50 of the PTEN-mutant lines was lower than that of WT (Fig. 2B; Supplementary Fig. S3C). The mouse cancer lines MCCL-357 (Myc, Pten−/−) and CaP8 (Pten−/−) were also more sensitive than the mouse cancer lines MCCL-278 (Myc, Pik3caH1047R) and Myc-CaP (Myc) were (Fig. 2C; Supplementary Fig. S3D and S3E; refs. 20, 21). Moreover, Pten−/− MEFs, PTEN-mutant human breast cancer cell lines, and Pten−/− mouse breast cancer cell lines displayed increased cell death over time upon treatment with leflunomide (Fig. 2D and E; Supplementary Fig. S3F). It is important to note that sensitivity to leflunomide was not associated with the proliferation rates of human breast, mouse breast, or mouse prostate tumor cell lines (Fig. 2F; Supplementary Fig. S3G and S3H). In addition, consistent with previous reports (22), we found that Pten homozygous deletion caused greater AKT phosphorylation than Pik3ca missense mutation did. This was particularly prominent in the nuclear fractions, where AKT may phosphorylate nuclear substrates (Fig. 2G; Supplementary Fig. S3I and S3J).

To independently test whether DHODH inhibition is detrimental to PTEN-deficient cells, we performed a rescue experiment with orotate, the metabolite directly downstream of DHODH, as well as with uridine. Increasing concentrations of orotate or uridine rescued growth inhibition by leflunomide (Fig. 2H–K; Supplementary Fig. S4A–S4C). In addition, siRNA against DHODH preferentially killed PTEN-mutant cells, verifying that DHODH was the target of the small-molecule inhibitors (Supplementary Fig. S4D). There was no difference in endogenous DHODH protein level between Pten−/− and WT MEFs, and A771726 did not affect PI3K signaling (Supplementary Fig. S4E and S4F). Consistent with prior reports, CAD phosphorylation downstream of mTORC1 was increased in Pten−/− cells, likely contributing to the push of glutamine flux into the pyrimidine synthesis pathway (Supplementary Fig. S4G and S4H; ref. 5). Pten−/− cells were more sensitive than WT cells to the mTOR inhibitor RAD001, as expected, but RAD001 did not synergize with leflunomide (Supplementary Fig. S4I–S4K; refs. 10, 23). Interestingly, treatment with nucleotide analogue inhibitors, 5-fluorouracil or mercaptopurine, did not show a differential sensitivity, demonstrating that Pten−/− MEFs are selectively vulnerable to inhibition of de novo pyrimidine synthesis (Supplementary Fig. S4L and S4M).

MYC activation is known to cause glutamine addiction (4). CaP8 (Pten−/−) cells were nearly as sensitive to glutamine deprivation as Myc-CaP (Myc oncogene transformed) cells were, substantiating that a notable level of glutamine dependency is also elicited by PTEN loss (Supplementary Fig. S4N). As Myc-CaP cells were resistant to leflunomide, it seems it is not the entry of glutamine alone but its flux into pyrimidines that is important (Supplementary Fig. S3D). Although MYC is known to largely direct glutamine to the tricarboxylic acid cycle and phospholipid synthesis (4), our data suggest that Pten loss in MEFs causes glutamine to cascade through the de novo pyrimidine synthesis pathway, creating the point of vulnerability to DHODH inhibition.

To determine how clinically relevant leflunomide may be as a targeted cancer therapy, we grew patient-derived glioblastomas as three-dimensional neurospheres. Formation of neurospheres was inhibited at lower concentrations of leflunomide in PTEN-deficient samples (Fig. 3A; Supplementary Fig. S5A). In addition, we treated two PTEN-mutant triple-negative breast cancer xenografts with leflunomide, dosing orally as is done clinically. Tumors slowed or regressed upon treatment; remarkably, even very large tumors (4 × 107 photons) regressed after only 1 week of treatment, indicating that leflunomide may have use for neoadjuvant therapy (Fig. 3B and C; Supplementary Fig. S5B). To ensure the effect in vivo is specific to PTEN loss, MCCL-357 and MCCL-278 xenografts were treated with leflunomide; MCCL-357 xenografts had a 4-fold better response than MCCL-278 xenografts did (Supplementary Fig. S5C).

Figure 3.

A, Dispersed (single-cell suspension) glioblastomas were treated with DMSO or 50, 100, or 200 μmol/L leflunomide for 5 days. The number of formed three-dimensional tumor spheres was quantified and normalized to untreated samples. (Student t test, *, P < 0.05, n = 3). B, SUM149 xenografts. Mice were treated with 100 mg/kg leflunomide or vehicle on days indicated with arrows (one-way ANOVA with multiple t tests, corrected for multiple comparisons, *, P < 0.01 for ANOVA and t tests, n = 6). C, MDA-MB-468 xenografts expressing luciferase, normalized to control. Treatment was started on day 7, with 100 mg/kg leflunomide or vehicle for four consecutive days each week (one-way ANOVA with multiple t tests, corrected for multiple comparisons; *, P < 0.05 for ANOVA and t tests, n = 5). Right, luminescence of treated and control mice after 2 weeks of treatment. Data, means ± SD for A and ± SEM for B and C.

Figure 3.

A, Dispersed (single-cell suspension) glioblastomas were treated with DMSO or 50, 100, or 200 μmol/L leflunomide for 5 days. The number of formed three-dimensional tumor spheres was quantified and normalized to untreated samples. (Student t test, *, P < 0.05, n = 3). B, SUM149 xenografts. Mice were treated with 100 mg/kg leflunomide or vehicle on days indicated with arrows (one-way ANOVA with multiple t tests, corrected for multiple comparisons, *, P < 0.01 for ANOVA and t tests, n = 6). C, MDA-MB-468 xenografts expressing luciferase, normalized to control. Treatment was started on day 7, with 100 mg/kg leflunomide or vehicle for four consecutive days each week (one-way ANOVA with multiple t tests, corrected for multiple comparisons; *, P < 0.05 for ANOVA and t tests, n = 5). Right, luminescence of treated and control mice after 2 weeks of treatment. Data, means ± SD for A and ± SEM for B and C.

Close modal

It is logical that a blockade of pyrimidine synthesis would stop cells from dividing, and leflunomide has been previously established as a cytostatic drug (18). What is more enigmatic, however, is why it would cause PTEN−/− cells to die. Consistent with prior reports (24), Pten−/− MEFs had a higher level of γH2AX, an indicator of DNA damage (Fig. 4A). We hypothesized that the dearth of pyrimidine deoxynucleotides caused by DHODH inhibition would exacerbate this defect, and discovered that leflunomide (or A771726) augmented DNA damage to a significantly greater degree in PTEN-deficient cells and that this damage colocalized with replication forks labeled with EdU (Fig. 4B–D; Supplementary Fig. S6A and S6B). Leflunomide-induced DNA damage was rescued by uridine, demonstrating that damage is likely instigated by pyrimidine depletion (Fig. 4E). The greater number of replication forks we described in Pten−/− MEFs remained intact after 24 hours of treatment with leflunomide, showing that the cells continue to replicate despite the presence of DNA damage (Figs. 1B and 4F; Supplementary Fig. S6C and S6D).

Figure 4.

A, Cells were labeled with a γH2AX antibody. Flow cytometry determined the mean fluorescence intensity (MFI; Student t test, *, P < 0.05, n = 3). B and C, Cells treated with 100 μmol/L leflunomide or A771726 were labeled with a γH2AX antibody. Flow cytometry determined the mean fluorescence intensity (Student t test, *, P values in figures, n = 3). D, MEFs treated with 150 μmol/L A771726 for 24 hours, labeled with EdU and γH2AX. Left, representative confocal microscopy images; right, quantified EdU and γH2AX-colocalized foci (Student t test; *, P < 0.05, n = 3). E, Cells treated with 100 μmol/L leflunomide with or without uridine and labeled with a γH2AX antibody. Flow cytometry determined the mean fluorescence intensity (Student t test, *, P values on figures, n = 3). F, MEFs treated with 100 μmol/L leflunomide or control for 48 hours and labeled with EdU. Left, representative confocal microscopy images; right, quantification of the number of foci per cell (Student t test, P > 0.05, n = 6). G, Cells were labeled with a pTOPBP1 S1159 antibody. Flow cytometry determined the mean fluorescence intensity (Student t test, *, P < 0.05, n = 3). H, Cells treated with 150 μmol/L A771726 for times indicated and labeled with antibodies to RPA and γH2AX. Flow cytometry determined the percentage of the cell population positively stained for RPA alone or both RPA and γH2AX (Student t test; *, P < 0.05, n = 4). I, pCHK1 immunoblot after 150 μmol/L A771726 treatment for times indicated. J and K, Quantified chromosomal breaks and multiradial formations per haploid genome (Student t test; *, P values on figure, cells scored/replicate > 100). L,PTEN-mutant cells were transfected with either WT TOPBP1 and CHK1, or mutants incapable of being phosphorylated by AKT, and labeled with a γH2AX antibody after 100 μmol/L leflunomide treatment. Flow cytometry determined the mean fluorescence intensity (Student t test; *, P < 0.05, n = 3). M,PTEN-mutant cells were transfected with either WT TOPBP1 and CHK1, or mutants incapable of being phosphorylated by AKT, and DRAQ7 was used to monitor accumulation of cell death in intervals of 6 hours (one-way ANOVA, *, P < 0.05, n = 3). Data, means ± SD.

Figure 4.

A, Cells were labeled with a γH2AX antibody. Flow cytometry determined the mean fluorescence intensity (MFI; Student t test, *, P < 0.05, n = 3). B and C, Cells treated with 100 μmol/L leflunomide or A771726 were labeled with a γH2AX antibody. Flow cytometry determined the mean fluorescence intensity (Student t test, *, P values in figures, n = 3). D, MEFs treated with 150 μmol/L A771726 for 24 hours, labeled with EdU and γH2AX. Left, representative confocal microscopy images; right, quantified EdU and γH2AX-colocalized foci (Student t test; *, P < 0.05, n = 3). E, Cells treated with 100 μmol/L leflunomide with or without uridine and labeled with a γH2AX antibody. Flow cytometry determined the mean fluorescence intensity (Student t test, *, P values on figures, n = 3). F, MEFs treated with 100 μmol/L leflunomide or control for 48 hours and labeled with EdU. Left, representative confocal microscopy images; right, quantification of the number of foci per cell (Student t test, P > 0.05, n = 6). G, Cells were labeled with a pTOPBP1 S1159 antibody. Flow cytometry determined the mean fluorescence intensity (Student t test, *, P < 0.05, n = 3). H, Cells treated with 150 μmol/L A771726 for times indicated and labeled with antibodies to RPA and γH2AX. Flow cytometry determined the percentage of the cell population positively stained for RPA alone or both RPA and γH2AX (Student t test; *, P < 0.05, n = 4). I, pCHK1 immunoblot after 150 μmol/L A771726 treatment for times indicated. J and K, Quantified chromosomal breaks and multiradial formations per haploid genome (Student t test; *, P values on figure, cells scored/replicate > 100). L,PTEN-mutant cells were transfected with either WT TOPBP1 and CHK1, or mutants incapable of being phosphorylated by AKT, and labeled with a γH2AX antibody after 100 μmol/L leflunomide treatment. Flow cytometry determined the mean fluorescence intensity (Student t test; *, P < 0.05, n = 3). M,PTEN-mutant cells were transfected with either WT TOPBP1 and CHK1, or mutants incapable of being phosphorylated by AKT, and DRAQ7 was used to monitor accumulation of cell death in intervals of 6 hours (one-way ANOVA, *, P < 0.05, n = 3). Data, means ± SD.

Close modal

Depletion of nucleotide pools normally activates the ATR checkpoint at replication forks in S-phase cells (25). ATR checkpoint activation at stalled forks requires two signals, one through single-strand DNA-binding protein [replication protein A (RPA)] interaction with single-strand DNA to recruit the ATRIP–ATR complex, and a second signal through TOPBP1 interaction with the ATR activation domain (26, 27). Deletion of PTEN in cells is known to cause poor ATR checkpoint activation, and AKT phosphorylation of TOPBP1 on serine 1159 and CHK1 on serine 280 inhibits their function (14, 24, 28, 29). Prior work showed increased CHK1 serine 280 phosphorylation in PTEN/− cells that reduced CHK1 activity (24). We found greater TOPBP1 S1159 phosphorylation and concomitantly less TOPBP1 localization to replication forks in Pten−/− cells (Fig. 4G; Supplementary Fig. S6E and S6F). Diminished AKT activity through PI3K inhibition also reduced leflunomide-induced DNA damage (Supplementary Fig. S6G). To further investigate the response to DNA damage occurring at Pten−/− forks, we examined the interaction of RPA and γH2AX by flow cytometry. An increase in RPA signal was first achieved regardless of PTEN genotype in the presence of A771726, followed by a shift toward both RPA and γH2AX-positive cells in Pten−/− MCCL-357, but not in Pten WT MCCL-278 cells (Fig. 4H). Moreover, γH2AX appeared almost exclusively in RPA-positive MCCL-357 cells treated with A771726 (Supplementary Fig. S6H). A771726 also triggered ATR phosphorylation of CHK1 at serine 345 in Pten WT but to a much lesser extent in Pten−/− cells (Fig. 4I; Supplementary Fig. S6I). Thus, Pten−/− cells appear to be incapable of generating an appropriate activation of the ATR–CHK1 checkpoint at replication forks. Activation of CHK1 in MCCL-278 cells declined as RPA declined, suggesting that Pten WT cells eventually recovered from DHODH inhibition, whereas Pten−/− cells instead accumulated damage at 18 hours (Fig. 4I). By 48 hours, this genomic stress manifested in a greater number of chromosome gaps, breaks, and multiradial formations in MCCL-357 cells treated with A771726 compared with MCCL-278 cells (Fig. 4J and K; Supplementary Fig. S6J and S6K). These findings are consistent with the sensitivity to hydroxyurea that occurs in the setting of an ATR inhibitor (30). Furthermore, we were able to rescue DNA damage and cell death in leflunomide-treated PTEN-mutant cells by transfecting cells with TOPBP1 and CHK1 mutants incapable of being phosphorylated by AKT (S1159A and S280A, respectively), demonstrating that the synthetic lethality between pyrimidine depletion and mutation of PTEN is due to the AKT-mediated defects in the ATR pathway (Fig. 4L and M).

On the basis of our data, we propose that the inhibition of DHODH in PTEN-deficient cells first causes stalled forks due to inadequate nucleotide pools required to support replication, and that sustained treatment leads to insufficient ATR activation due to AKT phosphorylation of TOPBP1 and CHK1, leading to a buildup of DNA damage and cell death. PTEN WT cells do not exhibit this dependency on pyrimidine synthesis and have fewer forks per cell, perhaps because ATR–CHK1 coordinates origin firing during S-phase (31). In PTEN WT cells, treatment initially increased the RPA signal and triggered transient phosphorylation of CHK1, whereas longer treatment led to abated RPA with little concurrent increase in γH2AX, explaining the largely unaffected WT population upon DHODH inhibition (Supplementary Fig. S7). Although Pik3ca-mutant cells also exhibit AKT signaling, their relative resistance to DHODH inhibitors suggests that a PI3K signaling dosage effect due to their lower level of AKT activation may be important.

PTEN loss leads to chemoresistance in prostate cancer, a poorer response to trastuzumab in triple-negative breast cancer, and a shorter survival time in patients with gliomas (32). Targeting consequences of PTEN loss could be beneficial, particularly as the standard of care for the aforementioned cancers is primarily chemotherapy and radiation. Inhibiting DHODH has the advantage of affecting a specific pathway of glutamine flux downstream of glutaminase, thus preserving glutamine's other important functions in the cell. This increases the specificity of DHODH inhibitors to cells which are dependent on glutamine's role in pyrimidine synthesis per se, and is perhaps why their toxicity is low enough to be taken as a daily medication by patients with rheumatoid arthritis or multiple sclerosis (19). It has been shown that leflunomide inhibits B-cell and T-cell proliferation, contributing to its immunomodulating effects (33, 34). It is possible that preexisting B and T cells can still function in the presence of DHODH inhibitors, thus potentially arguing for the benefit of immunotherapy in combination with a DHODH inhibitor, perhaps in a metronomic therapy pattern. Here, we show that high activation of AKT toward TOPBP1 and CHK1 that downregulates ATR activation at replication forks, compounded with enhanced pyrimidine flux that both occur as a consequence of PTEN inactivation contribute to the observed synthetic lethality between PTEN mutation and DHODH inhibition. We hope that DHODH inhibitors will be a promising therapy for patients with PTEN-deficient cancers.

Immunoblotting

Samples were lysed in 2× Laemelli sample buffer before separation by SDS-PAGE and transferring to PVDF membranes, blocked with 10% milk in TBST for 1 hour, and incubated with primary antibodies overnight. The following antibodies were used: PTEN 6H2.1 (Millipore 04-035; 1:1,000), DHODH (Protein Tech 14877-1-AP; 0.5 ng/uL), vinculin (Sigma; 1:50,000), pCHK1 (Cell Signaling Technology, 2341; 1:500), CHK1 G-4 (Santa Cruz Biotechnology, sc-8408; 0.4 ng/uL), pCAD (Cell Signaling Technology, 12662; 1:500), CAD (Cell Signaling Technology, 11933; 1:500). HRP-conjugated secondary antibodies were used to detect protein signals.

Cell Culture

MEFs and mouse breast tumor lines: DMEM (Corning mt10013cv) supplemented with 10% FBS (Atlanta Biologicals), 1% pencillin/streptomycin (Fisher Scientific, 30002ci), and 2 mmol/L l-glutamine (total 6 mmol/L; Fisher Scientific MT25005CI). MDA-MB-468, MDA-MB-231, Myc-CaP (2015), and U87: DMEM supplemented with 10% FBS and 1% penicillin/streptomycin. HCC1419, HCC1187, HCC 1937, HCC 1806, BT549, ZR75-1, PC3, LNCAP, DBTRG: RPMI (Fisher Scientific, 10040cv) supplemented with 10% FBS and 1% penicillin/streptomycin. CaP8 cells (2015): DMEM with 10% FBS, 1% penicillin/streptomycin, and 5 μg/mL insulin (Sigma I9278). Neurospheres: stem cell media with 10 μg/mL FGF (R&D Systems 233-FB-025), 20 μg/mL EGF (PeproTech AF-100-15), and heparin (from R. Yong, May 2015). All cells were cultured in a 37°C incubator with humidity and 5% CO2. Cell lines were obtained from ATCC (which authenticates cell lines using several methods, including DNA fingerprinting) in 2006, with the exception of MEFs, MCCL-278, and MCCL-357, which were produced in our laboratory from mice (2012–2016). Cell lines were clear of Mycoplasma as determined by the Lonza kit (LT07-418) within 6 months of their use. Cell lines were further authenticated in 2015 by LabCorp using a short tandem repeat method.

Mouse Embryonic Fibroblasts

Embryos were harvested 14 days after conception from Pten flox/flox mice from Jackson Labs. MEFs were treated with adenovirus diluted in growth media supplemented with polybrene with or without Cre recombinase and studied passages 2–5 post infection. Please see Supplementary Data for additional details.

Proliferation Assay

Mouse cells (1,500 cells per well) or human cells (3,000 cells per well) were plated in 96-well plates (Corning 720089). Growth rates were determined using the phase-confluency readings on an IncuCyte ZOOM (Essen Biosciences) on live cells over time.

Metabolite Labeling

For glutamine flux, media without added glutamine (Corning 17-207-CV) were supplemented with 13C glutamine (fully labeled) or 15N glutamine (amide labeled; Cambridge Isotope Labs). For glucose flux, media without added glucose (Corning 17-207-CV) were supplemented with 13C glucose (fully labeled; Cambridge Isotope Labs). Cells were plated in 10-cm dishes and grown in normal media. One hour prior to metabolite extraction, media were aspirated and replaced with heavy isotope–labeled media.

Metabolic Extraction

Metabolites were extracted in methanol. Please see Supplementary Data for details.

Targeted Mass Spectrometry

Mass spectrometry was performed by the core facility at Beth Israel Deaconess Medical Center (New York, NY). Please see Supplementary Data for details.

Cell-Cycle Analysis

Instructions for the FlowCellect Bivariate Cell Cycle Kit (Millipore FCCH025102) were followed. Fluorescence was measured on a Guava flow cytometer. Bromodeoxyuridine was pulsed for 18 hours.

Cell Death

Instructions for the FlowCellect Annexin Red Kit (Millipore FCCH100108) were followed. Fluorescence was measured on a Guava flow cytometer.

Seahorse Analysis

A Seahorse XF (Agilent) was used to determine the oxygen consumption and extracellular acidification rates. Cells were plated in the Seahorse cartridge and pH-adjusted media were added. Oligomycin, FCCP, and rotenone (Seahorse XF Kit, Agilent) were injected into cartridge wells. Controls for calibration and establishing baseline were used.

Drug Response Assays

Cells were plated in 96-well plates at a density of 1,500 or 3,000 cells per well. Leflunomide (Sigma PHR1378-1G), A771726 (Sigma SML0936), mercaptopurine (Sigma 852678), brequinar (Sigma SML0113), 5-fluorouracil (Millipore 343922), RAD001 and GDC0941 (obtained from Stand Up to Cancer PI3K Dream Team pharmacy), and CB-839 (MedChemexpress HY-12248) were dissolved in DMSO. Sensitivity was determined by a dose–response titration for each cell line, with an equivalent amount of DMSO in each well: 300 μL media with drug was added to one column of wells, and 150 μL media with equivalent DMSO was added to remaining wells. Serial dilutions of 150 μL resulted in a gradient with half the drug concentration as the previous column while maintaining the same amount of DMSO. GI50 values were calculated by linear interpolation: The maximum growth confluence for a cell line prior to growth plateau was divided by 2 to obtain the 50% confluence value. A linear regression curve was calculated using drug concentrations as x-values and confluence as y-values for points surrounding the 50% value. Linear interpolation using the regression line yielded the GI50 concentration. For cell death assays, DRAQ7 (Cell Signaling Technology, 7406S) was added to the media at a 1:200 dilution and red fluorescence was measured in addition to phase in live-cell imaging to measure accumulation of dead cells. An IncuCyte ZOOM was used.

γH2AX Measurement

Instructions for the FlowCellect Cell Cycle Checkpoint H2A.X DNA Damage Kit (Millipore FCCH12542) were followed. Briefly, cells were fixed and permeabilized, followed by staining with an anti-phospho-H2A.X antibody and propidium iodide. For costaining with RPA, an additional step was performed during which the cells were incubated with an RPA antibody (Abcam ab79398) for 1 hour and a secondary antibody for 1 hour (propidium iodide was not used in this setting). Fluorescence was measured on a Guava flow cytometer.

EdU Detection

Instructions for the EdU Cell Proliferation Kit (Millipore 17-10525) were followed. Cells were fixed and permeabilized following a 45- minute EdU pulse and a click chemistry reaction was used to add a fluorescent tag. Fluorescence was measured on a Guava flow cytometer or by immunofluorescence.

Immunofluorescence

Cells were plated on cover slips in media. For detecting replication forks, following a 45-minute EdU pulse, cover slip–attached cells were fixed and permeabilized, and detected after azide conjugation to EdU. For detecting γH2AX or TOPBP1, cells were incubated with primary antibody (BioLegend and Bethyl A300-111A-M, respectively) overnight at 4°C and with secondary antibody for 2 hours at room temperature. Images were taken using a Zeiss LSM880 Airyscan confocal microscope at 63×, and foci number and colocalization was quantified with ImageJ.

pTOPBP1 Measurement

Cells were fixed and permeabilized, followed by incubation with primary antibody (Abgent AP3774a) for 2 hours at room temperature and secondary antibody for 2 hours at room temperature. Fluorescence was measured on a Guava flow cytometer.

Karyotyping

Chromosomal analysis was performed on cultured cells treated with colcemid by Dr. Murty Vundavalli at Columbia University (New York, NY). Please see Supplementary Data for details.

Orotate Rescue

Orotate (Sigma O2750) was dissolved in DMSO. Cells were plated at fixed-concentration leflunomide with increasing concentrations of orotate, keeping DMSO constant in all wells.

Uridine Rescue

Uridine (Sigma U 3750) was dissolved in media. Cells were plated at fixed-concentration leflunomide with increasing concentrations of uridine, keeping DMSO constant in all wells.

RNA Interference

siRNA for DHODH was purchased from Qiagen. Cells were transfected using Lipofectamine (Invitrogen 11668-019) and knockdown was confirmed at 48 hours. Scrambled siRNA was used as a control.

Transfection

Plasmids were electroporated into cells (1 million) using an Amaxa Nucleofector 2b (Lonza) and Cell Line Nucleofector Kit V (Lonza VCA-1003). An mCherry plasmid was cotransfected to determine transfection efficiency and to gate transfected cells for flow cytometry experiments.

Xenografts

Six-week-old female nu/nu mice were engrafted orthotopically with either 5 million SUM149, 5 million MDA-MB-468-luciferase, 1 million MCCL-357, or 0.75 million MCCL-278 cells. Mice were treated by oral gavage with 100 mg/kg leflunomide or vehicle (1% carboxymethylcelluose in water). Animal experiments were approved by the Institutional Animal Care and Use Committee of Icahn School of Medicine at Mount Sinai.

Neurosphere Sensitivity Assay

Neurospheres were disrupted by manual pipetting until single-cell suspension was achieved, and 10,000 cells/well were plated in low-attachment 6-well plates (Fisher Scientific 3471). After 5 days, neurospheres were counted; sphere-forming ability is an indicator of tumorigenicity. Dense clusters >0.05 mm in diameter were counted as true tumor spheres. Spheres were measured on ImageJ and quantified.

Statistical Analysis

ANOVA or Student t tests were used to test means between groups. Correction for multiple comparisons was added where needed. Analysis was done using GraphPad Prism 6 or Microsoft Excel.

L.C. Cantley is a member of the board of directors of Agios and the scientific advisory board at Petra, reports receiving commercial research grants from Astellas and Takeda, has ownership interest (including patents) in Agios and Petra, and is a consultant/advisory board member for Enilibrium. No potential conflicts of interest were disclosed by the other authors.

One of the Editors-in-Chief is an author on this article. In keeping with the AACR's editorial policy, the peer review of this submission was managed by a senior member of Cancer Discovery's editorial team; a member of the AACR Publications Committee rendered the final decision concerning acceptability.

Conception and design: D. Mathur, S. Schoenfeld, R. Parsons

Development of methodology: D. Mathur, S. Ozturk, S. Schoenfeld, J.M. Asara, R. Parsons

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): D. Mathur, E. Stratikopoulos, S. Ozturk, N. Steinbach, S. Pegno, S. Schoenfeld, R. Yong, V.V. Murty, J.M. Asara, R. Parsons

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): D. Mathur, E. Stratikopoulos, S. Ozturk, S. Schoenfeld, J.M. Asara, R. Parsons

Writing, review, and/or revision of the manuscript: D. Mathur, E. Stratikopoulos, N. Steinbach, S. Pegno, R. Yong, J.M. Asara, L.C. Cantley, R. Parsons

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): D. Mathur

Study supervision: D. Mathur, R. Parsons

We would like to thank the members of the Parsons laboratory for providing feedback on this manuscript, and D. Mathur's thesis committee for helping shape the project. Confocal experiments were performed at the Microscopy core facility at Icahn School of Medicine at Mount Sinai. We thank Min Yuan and Susanne Breitkopf for help with mass spectrometry work. We also thank Dr. Gerard Karsenty, Columbia University Medical Center, for use of his Seahorse XF and Dr. Grzegorz Sumara for training us to use it. TOPBP1 WT and mutant plasmids were graciously supplied by Dr. Weei-Chin Lin, Baylor College of Medicine.

This work was supported by R01CA082783, R01CA155117, and P01CA97403 (to R. Parsons), and partially supported by NIH grants 5P01CA120964 and 5P30CA006516 (to J.M. Asara), and R01 GM041890 and the Breast Cancer Research Foundation (to L.C. Cantley).

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.

1.
Dang
CV
. 
Links between metabolism and cancer
.
Genes Dev
2012
;
26
:
877
90
.
2.
Shim
H
,
Chun
YS
,
Lewis
BC
,
Dang
CV
. 
A unique glucose-dependent apoptotic pathway induced by c-Myc
.
Proc Natl Acad Sci U S A
1998
;
95
:
1511
6
.
3.
Elstrom
RL
,
Bauer
DE
,
Buzzai
M
,
Karnauskas
R
,
Harris
MH
,
Plas
DR
, et al
Akt stimulates aerobic glycolysis in cancer cells
.
Cancer Res
2004
;
64
:
3892
9
.
4.
Wise
DR
,
DeBerardinis
RJ
,
Mancuso
A
,
Sayed
N
,
Zhang
X-Y
,
Pfeiffer
HK
, et al
Myc regulates a transcriptional program that stimulates mitochondrial glutaminolysis and leads to glutamine addiction
.
Proc Natl Acad Sci U S A
2008
;
105
:
18782
7
.
5.
Ben-Sahra
I
,
Howell
JJ
,
Asara
JM
,
Manning
BD
. 
Stimulation of denovo pyrimidine synthesis by growth signaling through mTOR and S6K1
.
Science
2013
;
339
:
1323
8
.
6.
Li
J
,
Yen
C
,
Liaw
D
,
Podsypanina
K
,
Bose
S
,
Wang
SI
, et al
PTEN, a putative protein tyrosine phosphatase gene mutated in human brain, breast, and prostate cancer
.
Science
1997
;
275
:
1943
7
.
7.
Steck
PA
,
Pershouse
MA
,
Jasser
SA
,
Yung
WKA
,
Lin
H
,
Ligon
AH
, et al
Identification of a candidate tumour suppressor gene, MMAC1, at chromosome 10q23.3 that is mutated in multiple advanced cancers
.
Nat Genet
1997
;
15
:
356
62
.
8.
Maehama
T
,
Dixon
JE
. 
The tumor suppressor, PTEN/MMAC1, dephosphorylates the lipid second messenger, phosphatidylinositol 3, 4, 5-trisphosphate
.
J Biol Chem
1998
;
273
:
13375
8
.
9.
Stambolic
V
,
Suzuki
A
,
De La Pompa
JL
,
Brothers
GM
,
Mirtsos
C
,
Sasaki
T
, et al
Negative regulation of PKB/Akt-dependent cell survival by the tumor suppressor PTEN
.
Cell
1998
;
95
:
29
39
.
10.
Podsypanina
K
,
Lee
RT
,
Politis
C
,
Hennessy
I
,
Crane
A
,
Puc
J
, et al
An inhibitor of mTOR reduces neoplasia and normalizes p70/S6 kinase activity in Pten+/− mice
.
Proc Natl Acad Sci U S A
2001
;
98
:
10320
5
.
11.
Mense
SM
,
Barrows
D
,
Hodakoski
C
,
Steinbach
N
,
Schoenfeld
D
,
Su
W
, et al
PTEN inhibits PREX2-catalyzed activation of RAC1 to restrain tumor cell invasion
.
Sci Signal
2015
;
8
:
ra32
ra
.
12.
Wong
J
,
Kim
P
,
Peacock
J
,
Yau
T
,
Mui
A-F
,
Chung
S
, et al
Pten (phosphatase and tensin homologue gene) haploinsufficiency promotes insulin hypersensitivity
.
Diabetologia
2007
;
50
:
395
403
.
13.
Mendes-Pereira
AM
,
Martin
SA
,
Brough
R
,
McCarthy
A
,
Taylor
JR
,
Kim
JS
, et al
Synthetic lethal targeting of PTEN mutant cells with PARP inhibitors
.
EMBO Mol Med
2009
;
1
:
315
22
.
14.
Shtivelman
E
,
Sussman
J
,
Stokoe
D
. 
A role for PI 3-kinase and PKB activity in the G2/M phase of the cell cycle
.
Curr Biol
2002
;
12
:
919
24
.
15.
Bassi
C
,
Ho
J
,
Srikumar
T
,
Dowling
R
,
Gorrini
C
,
Miller
S
, et al
Nuclear PTEN controls DNA repair and sensitivity to genotoxic stress
.
Science
2013
;
341
:
395
9
.
16.
Cory
JG
,
Cory
AH
. 
Critical roles of glutamine as nitrogen donors in purine and pyrimidine nucleotide synthesis: asparaginase treatment in childhood acute lymphoblastic leukemia
.
In vivo
2006
;
20
:
587
9
.
17.
Stamato
TD
,
Patterson
D
. 
Biochemical genetic analysis of pyrimidine biosynthesis in mammalian cells. II. Isolation and characterization of a mutant of Chinese hamster ovary cells with defective dihydroorotate dehydrogenase (EC 1.3. 3.1) activity
.
J Cell Physiol
1979
;
98
:
459
68
.
18.
Greene
S
,
Watanabe
K
,
Braatz-Trulson
J
,
Lou
L
. 
Inhibition of dihydroorotate dehydrogenase by the immunosuppressive agent leflunomide
.
Biochem Pharmacol
1995
;
50
:
861
7
.
19.
Munier-Lehmann
Hln
,
Vidalain
P-O
,
Tangy
Fdr
,
Janin
YL
. 
On dihydroorotate dehydrogenases and their inhibitors and uses
.
J Med Chem
2013
;
56
:
3148
67
.
20.
Jiao
J
,
Wang
S
,
Qiao
R
,
Vivanco
I
,
Watson
PA
,
Sawyers
CL
, et al
Murine cell lines derived from Pten null prostate cancer show the critical role of PTEN in hormone refractory prostate cancer development
.
Cancer Res
2007
;
67
:
6083
91
.
21.
Stratikopoulos
EE
,
Dendy
M
,
Szabolcs
M
,
Khaykin
AJ
,
Lefebvre
C
,
Zhou
M-M
, et al
Kinase and BET inhibitors together clamp inhibition of PI3K signaling and overcome resistance to therapy
.
Cancer Cell
2015
;
27
:
837
51
.
22.
Stemke-Hale
K
,
Gonzalez-Angulo
AM
,
Lluch
A
,
Neve
RM
,
Kuo
W-L
,
Davies
M
, et al
An integrative genomic and proteomic analysis of PIK3CA, PTEN, and AKT mutations in breast cancer
.
Cancer Res
2008
;
68
:
6084
91
.
23.
Neshat
MS
,
Mellinghoff
IK
,
Tran
C
,
Stiles
B
,
Thomas
G
,
Petersen
R
, et al
Enhanced sensitivity of PTEN-deficient tumors to inhibition of FRAP/mTOR
.
Proc Natl Acad Sci U S A
2001
;
98
:
10314
9
.
24.
Puc
J
,
Parsons
R
. 
PTEN loss inhibits CHK1 to cause double stranded-DNA breaks in cells
.
Cell Cycle
2005
;
4
:
927
9
.
25.
Tibbetts
RS
,
Cortez
D
,
Brumbaugh
KM
,
Scully
R
,
Livingston
D
,
Elledge
SJ
, et al
Functional interactions between BRCA1 and the checkpoint kinase ATR during genotoxic stress
.
Genes Dev
2000
;
14
:
2989
3002
.
26.
Zou
L
,
Elledge
SJ
. 
Sensing DNA damage through ATRIP recognition of RPA-ssDNA complexes
.
Science
2003
;
300
:
1542
8
.
27.
Kumagai
A
,
Lee
J
,
Yoo
HY
,
Dunphy
WG
. 
TopBP1 activates the ATR-ATRIP complex
.
Cell
2006
;
124
:
943
55
.
28.
Liu
K
,
Graves
JD
,
Scott
JD
,
Li
R
,
Lin
W-C
. 
Akt switches TopBP1 function from checkpoint activation to transcriptional regulation through phosphoserine binding-mediated oligomerization
.
Mol Cell Biol
2013
;
33
:
4685
700
.
29.
Liu
K
,
Paik
JC
,
Wang
B
,
Lin
FT
,
Lin
WC
. 
Regulation of TopBP1 oligomerization by Akt/PKB for cell survival
.
EMBO J
2006
;
25
:
4795
807
.
30.
Toledo Luis
I
,
Altmeyer
M
,
Rask
M-B
,
Lukas
C
,
Larsen Dorthe
H
,
Povlsen Lou
K
, et al
ATR prohibits replication catastrophe by preventing global exhaustion of RPA
.
Cell
2013
;
155
:
1088
103
.
31.
Shechter
D
,
Costanzo
V
,
Gautier
J
. 
ATR and ATM regulate the timing of DNA replication origin firing
.
Nat Cell Biol
2004
;
6
:
648
55
.
32.
Keniry
M
,
Parsons
R
. 
The role of PTEN signaling perturbations in cancer and in targeted therapy
.
Oncogene
2008
;
27
:
5477
85
.
33.
Siemasko
KF
,
Chong
AS
,
Williams
JW
,
Bremer
EG
,
Finnegan
A
. 
Regulation of B cell function by the immunosuppressive agent leflunomide
.
Transplantation
1996
;
61
:
635
42
.
34.
Chong
AS
,
Finnegan
A
,
Jiang
X
,
Gebel
H
,
Sankary
HN
,
Foster
P
, et al
Leflunomide, a novel immunosuppressive agent: the mechanism of inhibition of T cell proliferation
.
Transplantation
1993
;
55
:
1361
6
.