Expanding antitumor therapeutic windows by targeting cancer-specific nicotinamide adenine dinucleotide phosphate-biogenesis pathways

Abstract

Nicotinamide adenine dinucleotide phosphate (NADPH) biogenesis is an essential mechanism by which both normal and cancer cells maintain redox balance. While antitumor approaches to treat cancers through elevated reactive oxygen species (ROS) are not new ideas, depleting specific NADPH-biogenesis pathways that control recovery and repair pathways are novel, viable approaches to enhance cancer therapy. However, to elicit efficacious therapies exploiting NADPH-biogenic pathways, it is crucial to understand and specifically define the roles of NADPH-biogenesis pathways used by cancer cells for survival or recovery from cell stress. It is equally important to select NADPH-biogenic pathways that are expendable or not utilized in normal tissue to avoid unwanted toxicity. Here, we address recent literature that demonstrates specific tumor-selective NADPH-biogenesis pathways that can be exploited using agents that target specific cancer cell pathways normally not utilized in normal cells. Defining NADPH-biogenesis profiles of specific cancer-types should enable novel strategies to exploit these therapeutic windows for increased efficacy against recalcitrant neoplastic disease, such as pancreatic cancers. Accomplishing the goal of using ROS as a weapon against cancer cells will also require agents, such as NQO1 bioactivatable drugs, that selectively induce elevated ROS levels in cancer cells, while normal cells are protected.

Keywords:

• reactive oxygen species (ROS)
• NQO1-bioactivatable drugs
• nicotinamide adenine dinucleotide phosphate (NADPH)
• glutathione (GSH)
• biogenic pathways
• antioxidant

Introduction

Reduced nicotinamide adenine dinucleotide phosphate (NADPH) is a necessary cofactor for anabolic reactions, such as lipid and nucleic acid biosynthesis. Additionally, NADPH provides reducing power to oxidation–reduction reactions necessary for protecting cancer cells against the accumulation of reactive oxygen species (ROS) produced during rapid cellular proliferation.¹

While increased ROS in cancer cells may be an important initiating event in carcinogenesis, excessive levels of ROS can be toxic and lead to cell death by causing irreversible damage to DNA, lipids, and proteins.^1–3 Many chemotherapeutic agents act by inducing excessive ROS damage in cancer cells, but lack the ability to differentiate between normal and tumor tissues, leading to a narrow therapeutic window.^4,5 In addition, some cancers in advanced stages may become resistant to intrinsic oxidative stress and can up-regulate canonical antioxidant defenses to protect against ROS-inducing agents. Reduced glutathione (GSH) and thioredoxin (TRX) are essential ROS scavenging molecules in cancer and in normal cells.⁶ GSH and TRX are necessary for peroxidases, thioreductases, and peroxiredoxins to detoxify ROS. GSH and TRX rely on continuous reduction from NADPH to sustain their function as ROS scavengers.⁶ Therefore, the strategies to inhibit NADPH-biogenesis may dramatically alter the ROS scavenging abilities of cancer cells and sensitize them to oxidative damage. However, to achieve therapeutic selectivity, NADPH must be modulated through tumor-specific NADPH-biogenesis pathways that are necessary for cancer cells, but expendable in normal cells. To this end, this review describes cancer-selective alterations in NADPH biogenesis, defines potential therapies that exploit these pathways to sensitize cancer to ROS damage, and provides a method to predict cancer-specific NADPH-biogenesis profiles. We will not focus on pharmacological modulation of de novo GSH and/or TRX pathways, as these topics have been comprehensively reviewed elsewhere.^7–9

NADPH-biogenesis pathways in normal vs cancer cells

Oxidative pentose phosphate pathway (PPP)

A key mechanism of NADPH generation in normal cells is through the oxidative arm of the PPP. The PPP consists of two phases: the oxidative phase and the non-oxidative phase. The non-oxidative phase produces ribose from glucose, while the oxidative phase generates two NADPH molecules for every glucose entering the pathway (Figure 1).¹⁰ NADPH produced from the oxidative PPP is essential for protection against ROS damage arising from mitochondrial respiration, ionizing radiation, and various xenobiotic agents.¹¹ In this pathway, glucose 6-phosphate dehydrogenase (G6PD) and 6-phosphogluconate dehydrogenase (6PGD) reduce NADP⁺ to NADPH while oxidizing glucose-6-phosphate (G6P) and carboxylating 6-phosphogluconate (6PG), respectively (Figure 1).^12,13

Figure 1 NADPH production from the oxidative PPP and one-carbon serine catabolism pathway.

Notes: Oxidative PPP uses glucose to generate NADPH via G6PD and 6PGD. G6PD is inhibited by then FDA-approved drug, 6-AN. NADP⁺ is generated through the NAD⁺ salvage pathway, where nicotinamide is converted to NMN via NAMPT. NADP⁺ is then formed by NADK. FK866 and GMX1778 inhibit NAMPT to block the production of NADP⁺, and therefore NADPH. During ROS stress, p53 positively regulates TIGAR to shunt glycolytic flux into the oxidative PPP. PKM2, which is overexpressed in many cancers, is inhibited by ROS, allowing glycolytic flux to be shuttled into the oxidative PPP for NADPH generation. The small-molecule compounds, ML-202/203/265, can positively modulate PKM2, thereby decreasing glycolytic flux into the oxidative PPP and blunting NADPH biogenesis during ROS.
Abbreviations: PPP, pentose phosphate pathway; NADPH, nicotinamide adenine dinucleotide phosphate; G6PD, glucose-6-phosphate dehydrogenase; 6PGD, 6-phosphogluconate dehydrogenase; 6-AN, 6-aminonicotinamide; NMN, nicotinamide mononucleotide; NAMPT, nicotinamide phosphoribosyltransferase; NADK, NAD⁺-kinase; ROS, reactive oxygen species; TIGAR, TP53-induced glycolysis and apoptosis regulator; PKM2, pyruvate kinase 2; G6P, glucose-6-phosphate; 6PG, 6-phosphogluconate; R5P, ribulose-5-phosphate; F16BP, fructose-1,6-bisphosphate; PEP, phosphoenolpyruvate; FDA, food and drug administration; NAD, nicotine adenine dinucleotide.

Pyruvate kinase (PK) is an essential glycolytic enzyme for conversion of phosphoenolpyruvate (PEP) to pyruvate (Figure 1). The M2 isoform of PK (PKM2) is found in many cancer cells and self-renewing cells, but is expressed in an inactive state in normal adult tissues.¹⁴ In many human cancers, PKM2 can be inactivated by ROS, which diverts glycolytic flux back into the oxidative PPP to generate NADPH and detoxify ROS (Figure 1).¹⁵ After ROS stress, PKM2 is an essential in cancer, but not normal, cells to maintain cell viability via redox scavenging. It could provide a potential efficacious antitumor therapeutic window for ROS-inducing agents.¹⁵ PKM2 overexpression ensures that rapidly proliferating cancer cells create enough NADPH to match oxidative metabolism-generating ATP production, protecting the cell from attack by oxidative damage.^16,17,11

The tumor suppressor, p53, can also regulate flux into the oxidative PPP. During genotoxic stress, p53 induces TP53-induced glycolysis and apoptosis regulator (TIGAR), which encodes a protein that degrades fructose-2,6-bisphosphate (Figure 1).¹⁸ Low fructose-2,6-bisphosphate levels inhibit the activity of phosphofructokinase 1 (PFK1), a rate-limiting enzyme in glycolysis that leads to shuttling of earlier glycolytic metabolites into the oxidative PPP to generate NADPH. Overexpression of TIGAR was observed in colon, breast, and glioblastoma cancers.^19–21 Consisten with the enzyme’s role in redox balance, TIGAR knockdown dramatically sensitized glioma cells to ionizing radiation.²²

In cancers that overexpress PKM2, activating PKM2 with ML202, ML203, or other PKM2 activators blocks inhibition of PKM2 from ROS-inducing agents and decreases the flux of glucose through the oxidative PPP. This attenuates production of NADPH during oxidative damage, thereby sensitizing cancer cells to ROS-inducing agents. In contrast, normal cells that have inactive PKM2 are not sensitized to ROS-inducing agents.^15,23,24 Alternatively, the FDA-approved G6PD inhibitor, 6-aminonicotinamide (6-AN), may be utilized in cancers with PKM2 or TIGAR overexpression, thus directly inhibiting NADPH production via the oxidative PPP pathway. This, in turn, also sensitizes cancer cells to ROS-inducing agents (Figures 1 and 2).²⁵ The utility of this latter strategy needs to be empirically determined as G6PD is a major NADPH source in normal cells as well and toxicity concern will be a major factor in its efficacy.

Figure 2 Agents targeting specific NADPH-biogenesis pathways.

Abbreviations: NADPH, nicotinamide adenine dinucleotide phosphate; BPTES, bis-2-(5-phenylacetamido-1,2,4-thiadiazol-2-yl)ethyl sulfide 3; PPP, pentose phosphate pathway; PKM2, pyruvate kinase 2; ME, malic enzyme; GLS1, glutaminase 1; NAMPT, nicotinamide phosphoribosyltransferase; TS, thymidylate synthase; DHFR, dihydrofolate reductase; TIGAR, TP53-induced glycolysis and apoptosis regulator; PDAC, pancreatic ductal adenocarcinomas; SHMT2, serine hydroxymethyltransferase 2; MTHFD2, methylenetetrahydrofolate dehydrogenase 2.

Serine catabolism

Serine-driven, one-carbon metabolism was recently shown to be a major source of NADPH in dividing cells.²⁶ Serine is metabolized in the cytoplasm or mitochondria to methylene-tetrahydrofolate (methylene-THF) by serine hydroxymethyltransferase (SHMT) 1 or 2 (cytoplasmic and mitochondrial, respectively), which then forms 10- formyl-THF via methylenetetrahydrofolate dehydrogenase (MTHFD) 1 or 2 (cytoplasmic and mitochondrial, respectively). The flux through MTHFD generates NADPH in the cytoplasm or mitochondria. 10-formyl-THF is an essential for purine biogenesis, and MTHFD’s most important function was thought to facilitate purine biosynthesis. However, the NADPH generated from this reaction is also an integral source of cellular reducing power in dividing cells (Figure 3), including normal tissues that turn over quickly, like the colon.²⁶

Figure 3 One-carbon serine catabolism pathway.

Notes: Mitochondrial one-carbon serine metabolism pathway highlighting the production of NADPH from serine and folic acid through mitochondrial MTHFD2. THF is produced from DHF via folic acid by DHFR, which is inhibited by PEM or MTX. Additionally, DHF can be generated from dUMP via TS, which is also inhibited by PEM.
Abbreviations: NADPH, nicotinamide adenine dinucleotide phosphate; MTHFD2, methylenetetrahydrofolate dehydrogenase; THF, tetrahydrofolate; DHF, dihydrofolate; DHFR, dihydrofolate reductase; PEM, pemetrexed; MTX, methotrexate; dUMP, deoxyuridine monophosphate; TS, thymidylate synthase; SHMT2, serine hydroxymethyltransferase 2.

Recently, it was reported that this serine catabolism pathway can regulate mitochondrial redox control during hypoxia in Myc-driven cancers.²⁷ Specifically, SHMT2 was essential in maintaining mitochondrial NADPH and reduced GSH levels during hypoxia. SHMT2 expression was transcriptionally regulated by the coordinated activities of Myc and HIF-1α. Indeed, silencing SHMT2 in neuroblastoma cell lines significantly decreased growth in vitro under hypoxic conditions, and in a xenograft model of neuroblastoma (Figure 3).²⁷ Additionally, the authors demonstrated that high SHMT2 levels correlated with a poorer prognostic outcome in neuroblastoma patients, providing a clinical context for targeting this pathway. While SHMT2 is a source of NADPH in normal dividing cells, inhibiting SHMT2 in normal cells should not significantly alter NADPH-biogenesis, since normal cells have robust compensatory mechanisms for redox balance, unlike cancer cells.²⁶

In Myc-driven neuroblastoma, inhibiting SHMT2 or MTHFD2 would decrease NADPH biogenesis derived from one-carbon serine catabolism. While there are currently no known specific inhibitors of SHMT2 or MTHFD2, targeting production of serine’s obligate reaction partner, THF, may offer a strategy to decrease NADPH production from serine catabolism in a tumor-selective manner. For example, inhibiting dihydrofolate reductase (DHFR) with the anti-folate methotrexate (MTX) will decrease THF production, thereby decreasing the flux through SHMT2 and MTHFD2. This would attenuate NADPH levels specifically in cancer cells with elevated SHMT2 expression (Figures 2 and 3). Indeed, MTX exposure leads to cytostasis in overactive inflammatory cells seen in autoimmune diseases by decreasing GSH production, presumably due to decreased NADPH-biogenesis.^28,29 Alternatively, the new-generation anti-folate, pemetrexed (PEM), can also attenuate NADPH production from THF/serine catabolism by inhibiting both thymidylate synthase (TS) and DHFR, enzymes essential in THF synthesis (Figures 2 and 3).³⁰ Utilizing anti-folate drugs in SHMT2 overexpressing cancers in combination with ROS-inducing agents may provide a robust antitumor therapeutic window to exploit using these agents.

Malic enzymes

Another source of cellular NADPH is the NADP-dependent family of malic enzymes. This family of enzymes catalyzes the oxidative decarboxylation of malate to generate CO₂ and pyruvate, while reducing NAD⁺ or NADP⁺ to NADH or NADPH in the process (Figure 4).³¹ Three isoforms were identified in mammalian systems: cytosolic NADP⁺- dependent (ME1), mitochondrial NAD(P)⁺-dependent (ME2), and mitochondrial NADP⁺-dependent malic enzyme (ME3).³²

Figure 4 KRAS-reprogrammed glutamine metabolism in pancreatic cancer.

Notes: In KRAS-mutant pancreatic cancer, mitochondrial glutamine flux is reprogrammed to predominantly flux through GOT2, instead of the canonical GLUD1 pathway. The Asp produced from this reaction is shuttled to the cytoplasm to produce pyruvate and NADPH from malate via cytosolic ME1. TCA cycle-derived malate can also produce NADPH from mitochondrial ME2. Wild-type p53 inhibits both ME1 and ME2. BPTES and CB-839 are non-competitive small-molecule inhibitors of mitochondrial GLS1.
Abbreviations: GOT2, mitochondrial aspartate transaminase; GLUD1, glutamate dehydrogenase 1; Asp, aspartate; NADPH, nicotinamide adenine dinucleotide phosphate; ME1, malic enzyme 1; ME2, malic enzyme 2; BPTES, bis-2-(5-phenylacetamido-1,2,4-thiadiazol-2-yl)ethyl sulfide 3; GLS1, glutaminase 1; MDH1, malate dehydrogenase 1; GOT1, cytosolic aspartate transaminase; OAA, oxaloacetate; αKG, α-ketoglutarate; TCA, tricarboxylic acid cycle.

A recent report demonstrated the requirement of the cytosolic malic enzyme (ME1) in utilizing glutamine as an upstream metabolite to generate NADPH and to maintain redox balance in KRAS-mutated pancreatic ductal adenocarcinomas (PDAC), but not in normal pancreatic tissue (Figure 4).³³ The canonical metabolism of glutamine generates α-ketoglutarate (αKG) via the upstream activity of glutamate dehydrogenase 1 (GLUD1) to drive anaplerosis to replenish the tricarboxylic acid cycle.³⁴ However, in KRAS-mutated PDAC, glutamine flux is primarily driven through mitochondrial aspartate transaminase (GOT2) to generate mitochondrial αKG and aspartate from glutamate and oxaloacetate (OAA). Aspartate is then shuttled to the cytoplasm and then acted on by cytosolic aspartate transaminase (GOT1), which is converted back to OAA (Figure 4). OAA is then converted to malate by malate-dehydrogenase 1 (MDH1) and then to pyruvate and NADPH by ME1 (Figure 4). The depletion of ME1 in these PDAC cancer cells suppressed cell line growth and tumor growth in vivo by ROS accumulation from loss of NADPH.³³ Furthermore, the inhibition of these enzymes in normal pancreatic cells did not significantly alter NADPH concentrations. Intriguingly, KRAS-mutated PDACs have dramatically decreased glucose flux into the oxidative PPP, suggesting that this non-canonical glutamine pathway is compensated for decreased NADPH production from the oxidative PPP.³⁵ In PDAC with activating KRAS mutations (which is ~90% of all PDACs), inhibiting ME1 decreases the utilization of glutamine for NADPH production and sensitizes cells to oxidative damage. While there are currently no known ME1 inhibitors, inhibiting the upstream utilization of glutamine via glutaminase 1 (GLS1) with bis-2-(5-phenylacetamido-1,2,4-thiadiazol-2-yl)ethyl sulfide 3 (BPTES), Compound 968, CB-839, or other GLS1 inhibitors would sensitize KRAS-mutated PDAC to ROS-inducing agents in a tumor-specific manner (Figures 2 and 4).^36–38

A recent report demonstrated that a subset of lung tumors overexpress ME2 relative to normal lung tissue. A similar overexpression of ME2 was observed in melanoma vs normal skin, suggesting an important role for ME2 in these cancer-types.³⁹ Indeed, when ME2 was knocked down in the A549 lung cancer cell line, the cellular NADPH/NADP⁺ ratio decreased three-fold compared to non-targeting control cells, indicative of a pro-oxidant state in the absence of ME2.

It was recently shown that ME1 and ME2 are negatively regulated by wild-type p53, and that the absence of a functional p53 led to a dramatic up-regulation of ME1/2 expression.⁴⁰ Consistent with this finding, the authors demonstrated that ME1/2 enzymes were essential for NADPH maintenance in the absence of functional p53.⁴⁰ In the context of cancer, this is an important observation as p53 is a commonly mutated tumor suppressor and the loss of its function may lead to a cancer cell-specific mechanism of NADPH biogenesis via ME1/2 de-repression.

Isocitrate dehydrogenases

NADPH production can also be driven by the conversion of isocitrate to αKG by NADP⁺-dependent cytosolic isocitrate dehydrogenase 1 (IDH1) and mitochondrial isocitrate dehydrogenase 2 (IDH2) (Figure 5).^41,42 While NADPH generation has well known roles in the reduction of ROS, αKG can also detoxify ROS by scavenging hydrogen peroxide (H₂O₂) through non-enzymatic decarboxylation to form water and succinate.^43,44

Figure 5 NADPH biogenesis via IDH1.

Note: Cytoplasmic NADPH generated by IDH1 and mitochondrial NADPH generated by IDH2.
Abbreviations: IDH1, isocitrate dehydrogenase 1; NADPH, nicotinamide adenine dinucleotide phosphate; IDH2, isocitrate dehydrogenase 2; αKG, α-ketoglutarate.

IDH1 and IDH2 are mutated in ~80% of cases of adult glioma and secondary glioblastoma, and in 30% of cases of acute myeloid leukemia (AML).^45,46 It was originally believed that these mutants led to enzymatic loss of function through dominant-negative inhibition of wild-type IDH1 and IDH2.^46–48 However, it is now believed that IDH1 and IDH2 mutants confer to these enzymes the ability to convert αKG to the novel oncometabolite, 2-hydroxyglutarate (2-HG).⁴⁹ This change causes mutated IDH1 and IDH2 enzymes to consume rather than produce NADPH, altering the cellular redox balance and leading to a pro-oxidant state in the cancer cell.^50,51 Additionally, overexpression of the IDH1 mutant protein in glioma cell lines sensitizes these cells to the ROS-inducing effects of ionizing radiation.⁵²

Glioma patients with IDH1^R132 mutations have prolonged survival compared to patients with wild-type IDH1.^47,50,53 A hypothesis for this observation could be that IDH1 mutants are defective in generating protective concentrations of NADPH to maintain reduced GSH and thus are more sensitive to oxidative damage. Thus, glioma or AML patients with a R132 mutation in IDH1 might benefit from ROS-inducing agents early during the course of treatment.

Nicotinamide phosphoribosyltransferase (NAMPT)

NADPH generation can be driven through the NAD⁺ salvage pathway via NAMPT, which catalyzes the transfer of the phosphoribosyl group from 5-phosphoribosyl-1-pyrophosphate to nicotinamide, forming nicotinamide mononucleotide (NMN), and pyrophosphate (Figure 1).⁵⁴ NAD⁺ generation can then be coupled with NAD⁺ kinase (NADK) activity to generate NADP⁺ that can then be reduced to NADPH through the enzymes discussed above (Figure 1).⁵⁵

Increased NAMPT expression was reported in colorectal, non-small cell lung (NSCL), prostate, and pancreatic cancers.^56–60 In these contexts, NAMPT has been shown to be an important source of reducing equivalents for redox balance within the cancer cell.⁵⁶ In fact, knockdown of NAMPT sensitized prostate and head and neck cancer cell lines to ROS induction from ionizing radiation.^56,61–63

NAMPT inhibitors are undergoing clinical trials as single-agent therapies, but recent results have, unfortunately, not been promising (Figure 2).⁶⁴ The NAMPT inhibitors, FK866, and GMX1778, may have the greatest efficacy when combined with ROS-inducing agents that take advantage of the pro-oxidant state of NAMPT-inhibited tumors. Indeed, pre-clinical studies have validated this strategy utilizing GMX1778 to sensitize breast cancer cells against the ROS production from ionizing radiation therapy, FK866 to sensitize prostate cancer cell against H₂O₂, and FK866 to sensitize neuroblastoma cells against cisplatin.^56,61,65 To enhance the selectivities of NAMPT inhibitors, NADPH:quinone oxidoreductase 1 (NQO1) bioactivatable drugs can be used in combination treatments. This results in cancer-specific lethality of cells that overexpress NQO1,⁶⁶ such as in pancreatic, NSCL, breast, prostate, and head and neck cancers.

Predicting tumor-specific NADPH-biogenesis profiles from publicly available datasets

Known NADPH-biogenesis pathways can be useful when combined with publically available cancer gene expression and patient outcomes data to generate hypotheses for tumor-specific NADPH-biogenesis profiles.

Here, we present a conservative set of criteria for determining candidate genes (Figure 6) using the following approach:

1. Determine if the genes of interest are significantly up-regulated in patient tumor tissue vs associated normal tissue in two or more independent datasets for the cancer-type in question. Only genes with P-values of <1×10⁻⁴ in each dataset will be considered for further analyses.

2. Of the up-regulated genes found in #1, determine which genes significantly predict poor outcomes in patients after radiation treatment.

3. Of the candidates from #2, determine if genes in the same NADPH biogenesis pathway as the candidate gene are coordinately up-regulated. Only positive correlation values of 0.5 and above will be considered for further analysis.

4. Determine if common oncogenic drivers of the tumor type in question drive genes from #3.

Figure 6 Work flow to identify NSCLC-specific NADPH-biogenesis genes in a cancer-specific NADPH-biogenesis screen.

Abbreviations: NSCLC, non-small cell lung cancer; NADPH, nicotinamide adenine dinucleotide phosphate.

Utilizing these criteria, we attempted to assess the NADPH-biogenesis profile of non-small cell lung cancer (NSCLC). Using the Oncomine webtool (http://www.oncomine.com), we first determined which NADPH-biogenesis related genes were overexpressed in NSCLC tissue vs associated normal tissue from patient samples. To ensure the robustness of our results, a gene was considered a potential hit only if it was found to be up-regulated in two or more datasets with a P-value <1×10⁻⁴ in each dataset. This narrowed our initial list of ten genes down to four potential genes (Figures 7 and 8A). Next, to determine if our potential targets were clinically relevant in the context of ROS-inducing therapy, we used the KMPLOT software and stratified NSCLC patients into high- and low-expressers of the genes of interest, and compared overall survival outcomes between these groups after radiation therapy (Figures 8A and 9).⁶⁷

Figure 7 A cancer-specific NADPH-biogenesis screen.

Notes: Overexpression status of NADPH-biogenesis genes of interest assessed using Oncomine in NSCLC patients. Datasets used are described.^82–84
Abbreviations: NADPH, nicotinamide adenine dinucleotide phosphate; NSCLC, non-small cell lung cancer; G6PD, glucose-6-phosphate dehydrogenase; PKM2, pyruvate kinase 2; TIGAR, TP53-induced glycolysis and apoptosis regulator; SHMT2, serine hydroxymethyltransferase 2; MTHFD2, methylenetetrahydrofolate dehydrogenase; ME1, malic enzyme 1; ME2, malic enzyme 2; IDH1, isocitrate dehydrogenase 1; IDH2, isocitrate dehydrogenase 2; NAMPT, nicotinamide phosphoribosyltransferase.

Figure 8 Cancer-specific NADPH-biogenesis screen continued.

Notes: (A) Prognosis after radiation therapy in NSCLC patients. Dataset used is described;⁶⁷ (B) Correlation of genes of interest with upstream NADPH-biogenesis pathway members; (C) Mutant KRAS-dependence of MTHFD2 in a mouse model of NSCLC.
Abbreviations: NADPH, nicotinamide adenine dinucleotide phosphate; NSCLC, non-small cell lung cancer; MTHFD2, methylenetetrahydrofolate dehydrogenase; PKM2, pyruvate kinase 2; SHMT2, serine hydroxymethyltransferase 2; IDH2, isocitrate dehydrogenase 2; G6PD, glucose-6-phosphate dehydrogenase; TS, thymidylate synthase.

Figure 9 Screen results for MTHFD2 in NSCLC patients.

Notes: (A and B) mRNA expression of MTHFD2 in normal lung tissue vs NSCLC tissue in two different datasets from Oncomine; (C) Overall survival after radiation therapy in NSCLC patients segregated in high and low MTHFD2 expressers. Dataset used is described;⁶⁷ (D) Correlation of MTHFD2 with upstream NADPH-biogenesis pathway members; (E) Mutant KRAS-dependence of MTHFD2 in a mouse model of NSCLC (GSE40606).
Abbreviations: MTHFD2, methylenetetrahydrofolate dehydrogenase; NSCLC, non-small cell lung cancer; NADPH, nicotinamide adenine dinucleotide phosphate; SHMT2, serine hydroxymethyltransferase 2; TS, thymidylate synthase; HR, hazard ratio.

From this analysis, we selected the genes whose high expression in NSCLC patients led to significantly decreased survival after radiation therapy, suggesting that these genes may confer tumor protection from radiation-induced ROS, presumably through enhanced NADPH biogenesis.⁶⁸ Using this cutoff, we were able to narrow down the gene list to PKM2, which regulates NADPH biogenesis via the oxidative PPP, and MTHFD2, which generates NADPH from serine catabolism (Figure 8B). We then determined if these genes were coordinately up-regulated with other enzymes in their respective NADPH-biogenesis pathways. Utilizing the r2 Genomics Analysis platform’s (r2.amc.nl) co-expression analysis feature in NSCLC patient samples, we found that G6PD was coordinately up-regulated with PKM2 in patient samples with a correlation value of 0.26, and that TS and SHMT2 were co-expressed with MTHFD2 with correlation values of 0.63 and 0.64, respectively (Figures 8B and 9D). TS and SHMT2 also appear to be co-expressed with a correlation value of 0.65. Given that we defined our correlation cutoff at 0.5 and above, we considered MTHFD2 as the top candidate for NSCLC (Figures 8B and 9D).

Next, we investigated whether mutant-KRAS, the most common oncogenic driver in NSCLC,⁶⁹ might regulate MTHFD2. Such an association would provide insight into regulatory mechanisms of MTHFD2 and additional therapeutic targets. For this, we utilized publically available gene expression datasets from the National Center for Biotechnology Information Gene Expression Omnibus (GEO). To determine if mutant-KRAS regulates MTHFD2, we analyzed the mRNA expression-profiling data from a transgenic mouse model of NSCLC expressing doxycycline-inducible KRAS^G12D in the respiratory epithelium (GSE40606).⁷⁰ When administered doxycycline, the mice develop lung tumors that are dependent on constitutive KRAS^G12D expression. Within 48 hours of doxycycline withdrawal, KRAS^G12D expression was extinguished and whole-genome gene expression analyses of lung tumors were performed. The expression levels of MTHFD2 were significantly up-regulated when KRAS^G12D was induced vs 48 hours after KRAS extinction with doxycycline withdrawal, indicating a positive regulatory role for mutant-KRAS in MTHDF2 expression (Figures 8C and 9E).

If this hypothesis is validated through RNAi and redox balance studies, it may suggest that the utilization of serine catabolism inhibitors, such as MTX or PEM, may provide an effective therapeutic strategy to target NADPH biogenesis, specifically in KRAS-mutated NSCLC. This analysis also reveals that more than one NADPH-biogenesis pathway may also be regulating NADPH as observed from the PKM2 data. However, validation experiments will need to be conducted to determine which pathway predominates in NSCLC.

Modulating NADPH biogenesis as a mechanism to potentiate NQO1-bioactivatable drugs

To exploit the metabolic vulnerabilities within a cancer as mentioned above, ROS-inducing agents with the capacity to specifically target tumors would be ideal. We believe that NQO1-bioactivatable drugs represent an ideal class of agents to combine with metabolic inhibitors. NQO1 is an inducible phase II detoxifying enzyme overexpressed in breast, lung, pancreatic, and colon cancers. The two-electron oxidoreductase is capable of reducing quinones by forming stable hydroquinones of the parent quinone.^71,72 β-Lapachone (β-lap, in clinical trials as ARQ761), IB-DNQ, and other NQO1-bioactivatable drugs are unique quinones that are metabolized by NQO1 into an unstable hydroquinone that spontaneously oxidizes back to the parental compound, generating a futile redox cycle in which 1 mole of β-lap generates ~120 moles of superoxide within 2 minutes, consuming >60 moles of NADH and/or NADPH.^73,74 The superoxide (O₂⁻) radicals formed are quickly metabolized by superoxide dismutase (SOD) into H₂O₂.^73,75 The massive levels of H₂O₂ formed causes extensive oxidative DNA damage that hyperactivates poly(ADP-ribosyl) polymerase 1 (PARP1), resulting in a dramatic loss of the intracellular NAD⁺/ATP pools and an inability to repair DNA damage.^76–78 Cell death is stimulated by caspase-independent NAD⁺-Keresis, a form of programmed necrosis.^66,73 Cancer cells with >100 units of NQO1 enzyme activity are sensitive to β-lap lethality, while normal tissues that lack, or express low levels of, NQO1 are spared.⁷⁹ While this class of drugs represents an attractive antitumor strategy, dose-limiting methemoglobinemia caused by non-specific ROS generation at high doses may limit its efficacy in monotherapy.⁸⁰ Strategies for increasing cancer cell- cytotoxicity while maintaining NQO1 specificity could greatly enhance the efficacy of β-lap for use in solid cancers that overexpress NQO1. Thus, combining β-lap with metabolic inhibitors that target cancer-specific NADPH-biogenesis pathways, such as FK866, CB-839, 6-AN, PEM or MTX, may synergistically expand the antitumor therapeutic window for NQO1 bioactivatable drugs, while increasing the tumor selectivity of metabolic inhibitors.

Conclusion

Cancer cells need to tightly regulate NADPH biogenesis to protect themselves against oxidative damage.⁸¹ To sustain protective levels of NADPH, cancer cells rely on various NADPH-biogenesis pathways, including oxidative PPP, serine catabolism, glutamine metabolism, and NAD⁺ salvage pathways (Figures 1, 3, and 4). Strategies to inhibit NADPH biogenesis may dramatically alter the ROS scavenging abilities of cancer cells and sensitize them to oxidative damage. However, to achieve therapeutic selectivity, NADPH must be modulated through tumor-specific NADPH-biogenesis pathways that are necessary for cancer cells but are expendable in normal cells. Thus, by rigorously studying these unique pathways in the context of a specific cancer, we will be able to create novel patient-specific antitumor therapeutic strategies that exploit the ROS balances of tumor tissue while sparing normal tissue in the process. Here, we summarized our current understanding of known cancer-specific NADPH-biogenesis pathways, drugs to specifically target these pathways, and an example of using publically available databases to predict cancer-type specific NADPH-biogenesis genes. It is our belief that studying these pathways and comprehensively profiling tumors based on this understanding will be an essential step forward in designing cancer-specific ROS combination therapies, such as the use of NQO1-bioactivatable drugs.

Acknowledgments

The authors thank Dr Ralph J DeBerardinis for his assistance in reviewing this work. Support for use of Shared Resources within the Simmons Comprehensive Cancer Center was funded by NIH/NCI grant 5P30CA142543. Support for these studies was provided by grant NIH/NCI R01 CA102972, and a grant from the Pancreatic Cancer Action Network (PanCan) to DAB. This work is dedicated to the memories of Mrs Sarah Hildebrand and Rosemary Bouley for their strong support of our research.

Disclosure

The authors declare no conflicts of interest in this work.