Abstract
Oxidative phosphorylation (OXPHOS) is an active metabolic pathway in many cancers. RNA from pretreatment biopsies from patients with triple-negative breast cancer (TNBC) who received neoadjuvant chemotherapy demonstrated that the top canonical pathway associated with worse outcome was higher expression of OXPHOS signature. IACS-10759, a novel inhibitor of OXPHOS, stabilized growth in multiple TNBC patient-derived xenografts (PDX). On gene expression profiling, all of the sensitive models displayed a basal-like 1 TNBC subtype. Expression of mitochondrial genes was significantly higher in sensitive PDXs. An in vivo functional genomics screen to identify synthetic lethal targets in tumors treated with IACS-10759 found several potential targets, including CDK4. We validated the antitumor efficacy of the combination of palbociclib, a CDK4/6 inhibitor, and IACS-10759 in vitro and in vivo. In addition, the combination of IACS-10759 and multikinase inhibitor cabozantinib had improved antitumor efficacy. Taken together, our data suggest that OXPHOS is a metabolic vulnerability in TNBC that may be leveraged with novel therapeutics in combination regimens.
These findings suggest that triple-negative breast cancer is highly reliant on OXPHOS and that inhibiting OXPHOS may be a novel approach to enhance efficacy of several targeted therapies.
Introduction
Approximately 10%–15% of women diagnosed with breast cancer present with triple-negative breast cancer (TNBC), which lacks expression of estrogen and progesterone receptors or amplification/overexpression of HER2 (1). There are limited targeted therapy options currently for patients with TNBC. For early-stage patients with TNBC, neoadjuvant (preoperative) cytotoxic chemotherapy can achieve complete pathologic responses (pCR) in 30%–40% of patients with TNBC, resulting in improved long-term outcomes (2). However, patients with TNBC who have significant residual disease following neoadjuvant chemotherapy have substantially worse oncologic outcomes (3). There is thus a clinical need to find factors associated with residual TNBC to identify novel therapeutic strategies to improve TNBC outcomes.
Recent studies have suggested that some tumor cells have an enhanced dependence on oxidative phosphorylation (OXPHOS; refs. 4–6). The increase in OXPHOS/mitochondria genes can provide/restore energy and metabolites (aspartate) needed for biosynthesis and proliferation (7). This increased reliance on OXPHOS can result from passenger genomic alterations associated with genomic deletions of neighboring tumor suppressors (7, 8). These coalterations result in decreased glycolytic potential alongside the tumorigenic events (4). However, other studies have shown at least a partial reliance on OXPHOS even without these genomic alterations (9, 10).
Here, we show that OXPHOS is associated with a higher risk of recurrence and death in TNBC using gene expression profiling of treatment-naïve tumors. We found that TNBC can be functionally reliant on OXPHOS using a panel of molecularly diverse patient-derived xenografts (PDX) generated from chemotherapy-treated, residual TNBC. IACS-10759 is a novel inhibitor of OXPHOS that is currently in clinical trials (11). IACS-10759 is clinical grade inhibitor of complex I of the mitochondrial electron transport chain. Molecular studies have shown that IACS-10759 inhibits the OXPHOS by binding mitochondrial respiratory complex 1 near or in the entrance of the ubiquinone channel thereby inhibiting ubiquinone function (7). We found that targeting OXPHOS with IACS-10759 inhibits growth of TNBC PDXs across a range of molecular subtypes with greatest antitumor efficacy seen in the basal-like TNBC molecular subtype and in tumors with high mitochondrial RNA expression. We systematically assessed potential predictive and pharmacodynamic markers of response to IACS-10759, using a panel of PDXs, to identify potential avenues to monitor ongoing OXPHOS inhibition. We also explored potential combination therapies by targeting putative mechanisms of intrinsic resistance and by performing an in vivo synthetic lethality screen for combination partners for IACS-10759 to improve antitumor efficacy.
Materials and Methods
Patient cohort for RNA sequencing
We performed RNA sequencing (RNA-seq) on pretreatment biopsies from 43 patients with operable TNBC who received sequential taxane- and anthracycline-based neoadjuvant chemotherapy. All patients gave informed consent for treatment on an open-label randomized clinical trial of standard neoadjuvant chemotherapy with paclitaxel followed by FEC versus the combination of paclitaxel and everolimus followed by FEC in women with TNBC (NCT00499603; ref. 12). The study did not show a difference in pCR rate between the arms, so pretreatment biopsies from both arms were analyzed together for analysis of transcriptomic features associated with relapse and overall survival.
PDXs
PDXs were generated from surgical samples obtained from patients after prospective informed consent as described previously (13). Investigators obtained informed written consent from all subjects. For this study, female athymic nude mice purchased from Harlan/Envigo were used for treatment testing. All animal experiments were approved by the University of Texas MD Anderson Cancer Center (Houston, TX) Institutional Animal Care and Use Committee. For AXL suppression in BCX.010, we tested three short hairpin RNAs (shRNA) in pGIPz lentiviral vector and chose one to proceed to in vivo testing. We used a firefly luciferase targeted control shRNA. We choose the cell lines with greatest knockdown of AXL by qRT-PCR and injected 2.5 × 106 cells on to the mammary fat pad of each mouse.
Animal treatments
IACS-10759 was suspended in 0.5% methycellulose and given orally for 5 days followed by 2 days off every week. Cabozantinib was dissolved in DMSO and then added to 0.5% methylcellulose to a final DMSO concentration of 10%. Cabozantinib was given orally daily. BGB-324 (50 mg/kg) was dissolved in 0.5% methylcellulose and given orally 7 days a week. Palbociclib (50 mg/kg) was dissolved in 0.5% methycellulose and given 7 days a week. Talazoparib (0.3 mg/kg) was dissolved in DMAC and then added to 20% Solutol, and given orally 7 days a week. Entinostat (20 mg/kg) was suspended in PBS and given orally 7 days a week. BAY87-2243 (5 mg/kg) was dissolved in 0.5% methycellulose and given orally for 5 days followed by 2 days off every week. Tumor volume was calculated by the formula: TV (mm3) = ((width)2 × length)/2.
Molecular analysis
For RNA analysis, small fragments of frozen tumor were placed in lysis buffer and homogenized manually. RNA was isolated from the frozen tumors using Norgen BIOTEK Total RNA Purification Plus Kit. RNA was quantified using Qubit system. Methodology and gene selection for Nanostring OXPHOS panel has been described previously (7).
Reverse phase protein arrays (RPPA) were performed at the MD Anderson Cancer Center Functional Proteomics Core Facility. Frozen tumor pieces were cut into approximately 3 × 3 mm size fragments and placed in bead lysis tubes for protein extraction. We chose the RPPA results originated from rabbit and goat antibodies (230 antibodies) for PDXs because validated mouse antibodies have significant background in PDX samples due to the mouse stromal tissue.
For IHC, we used monoclonal rabbit antibodies for phosphoHistone H3 (Bethyl, catalog no. IHC-00061), cleaved caspase 3 (CC3; Cell Signaling Technology #9664), and PCNA (KLJ). Sample were collected 4 hours after IACS-10759 treatment and 12 days after start of treatments. The assessment was performed by a clinically trained pathologist.
In vivo synthetic lethality screen
We utilized a barcoded FDAome shRNA library correspond to targets linked to drugs approved by the FDAto determine potential IACS-10759 synthetic lethal partners. Breast cancer cell line BCX.010-CL was derived from BCX.010 as described previously (14). BCX.010-CL cells were infected with our FDAome shRNA library, reinjected into mice (2 × 106 cells) and treated with either vehicle or IACS-10759. After 10 days, xenografts were harvested and DNA was extracted. Deep sequencing was performed to determine shRNA abundance (barcode abundance) on the basis of shRNA performance supervised analysis. Cumulative Wilcoxon test was done to score statistical gene level significance.
In vitro testing
Breast cancer cell lines were cultured in DMEM/F-12 supplemented with 10% FBS at 37°C and in a humidified incubator containing 5% CO2. BCX.010-CL cells were isolated as described previously (14); origin validation was done by short tandem repeat fingerprinting. MDA-MB-468 cells were obtained from ATCC. Laboratory cell cultures were regularly tested for mycoplasma and passaged a maximum of ten times after thawing early passage stocks. For in vitro testing, all agents were dissolved in DMSO. For cell viability assays, cells were plated into 96-well plates at 5,000 cells/well. The cells were allowed to attach overnight, and then treated with the selected drugs for 72 hours. Cell viability was evaluated by sulforhodamine B assay. IC50s from dose–response curves were determined using GraphPad Prism v.6.05 software. Combination index (CI) was calculated using the IC50-based Chou-Talalay model.
Results
Oxidative phosphorylation is associated with recurrence in patients with TNBC
To determine molecular features associated with higher risk of relapse in operable TNBC, we performed RNA sequencing on pretreatment biopsies from 43 patients with TNBC who received sequential taxane-and anthracycline-based neoadjuvant chemotherapy. At >5 years median follow-up, 14 patients had recurred and of those all but 2 had died. At a FDR of 0.05 (q value <0.05), 33 genes were differentially expressed between patients who did and did not have a subsequent recurrence (Fig. 1). Twenty seven genes were differentially expressed between patients based on overall survival (Supplementary Fig. S1).
One of the top canonical pathways that was associated with higher risk of relapse and lower risk of survival was higher expression of OXPHOS signature (Padjusted < 0.001) as determined by Ingenuity Pathway Analysis (IPA). The patients that recurred had significantly higher levels of expression of mitochondrial genes: MT-ND1 (adjusted P = 0.007); MT-ND5 (P = 0.03) and MT-ND4 (P = 0.04). Similarly, the patients that had died had higher expression of several mitochondrial genes: MT-ND5 (P = 0.001); MT-ND4 (P = 0.005), MT-ND4 L (P = 0.015), MT-ND6 (P = 0.018), and MT-ATP6 (P = 0.03). All of these mitochondrial genes (excluding MT-ATP6) encode members of Complex 1. For these reasons, we hypothesized that OXPHOS pathway, and specifically mitochondrial Complex 1, may represent a potential therapeutic target in TNBC.
Inhibiting oxidative phosphorylation has broad antitumor effect in TNBC PDXs
To determine the functional necessity of OXPHOS in TNBC, we treated a panel of 10 TNBC PDXs with OXPHOS inhibitor IACS-10759. Eight of these 10 PDXs were generated from residual tumors collected after neoadjuvant chemotherapy; these patients by definition did not have a complete pathological response and thus these tumors represent relatively chemoresistant disease. Although the models differ in their chemotherapy sensitivity, most have been previously demonstrated to be resistant to several standard chemotherapy regimens (13, 15). BCX.092 and BCX.024 were generated from patients that did not receive neoadjuvant therapy (13). Known clinically relevant genomic alterations in these models are listed in Supplementary Table S1. We treated the PDXs with IACS-10759 using a small cohort (n = 2–4) screening approach (Fig. 2). Five of 10 PDX models had either stable disease (did not increase ≥20% in size) or regressed at 28 days of treatment. One PDX (BCX.070), which is known to be resistant to several chemotherapeutic agents (13), drastically regressed with IACS-10759 treatment.
To elucidate the specificity of these effects in vivo, we treated a PDX that was less sensitive to IACS-10759 (BCX.011) and another that is more sensitive to IACS-10759 (BCX.017) with BAY87-2243, an alternative complex 1 inhibitor (16). Similar to IACS-10759, BAY87-2243 did not prevent progression of BCX.011 but did cause disease stabilization in BCX.017 (Supplementary Fig. S2). This suggests that indeed the differences we observed with IACS-10759 sensitivity may be representative of intrinsic differences in metabolism and reliance on OXPHOS for tumor growth.
Baseline biomarkers associated with intrinsic or acquired resistance and sensitivity to IACS-10759
We next performed analyses to identify potential predictive biomarkers of sensitivity using baseline RNA-seq data for the PDXs. We first determined the association between IACS-10759 sensitivity and clinically relevant RNA expression-based molecular subtypes of TNBC based on the Lehman classification as described previously (TNBCtype; refs. 13, 17). All of the IACS-10759 sensitive PDXs were TNBCtype basal-like 1 (BL1).
The baseline expression of the oxidative phosphorylation signature by IPA was not significantly different between IACS-10759 sensitive versus resistant PDX models. However, PDXs sensitive to IACS-10759 had higher expression of protein-encoding mitochondrial genes, including the specific mitochondrial genes that we found to be also differentially expressed higher in the patients more likely to develop recurrent disease (Fig. 3A; Supplementary Fig. S3).
As loss of alpha enolase (ENO1) expression has been linked to reliance on OXPHOS (7), we assessed ENO1 expression by RNA-seq and IHC. The more sensitive models trended towards lower ENO1 RNA expression, with the regressive model, BCX.070, having lowest expression (Supplementary Fig. S4A). The IACS-10759 sensitive models had lower ENO1 staining intensity by IHC (P < 0.001; Supplementary Fig. S4B). Using RNA-seq data, we assessed a published proliferation index (18); we found no significant difference in baseline proliferation between more and less sensitive PDXs (Supplementary Fig. S5).
Comparing more versus less IACS-10759–sensitive models, we identified 202 genes that are higher in the less sensitive models at FDR = 0.01 (Supplementary Table S2). By IPA, epithelial–mesenchymal transition (EMT) signature as well as hypoxia signatures were differentially expressed between sensitive and less sensitive PDXs (Supplementary Table S3). PDXs that were less sensitive to IACS-10759 expressed higher levels of TWIST1, a key transcriptional factor associated with EMT, and AXL, a tyrosine kinase that has been associated with EMT in breast cancer.
AXL expression has been associated with more aggressive tumor behavior, as well as therapeutic resistance to chemotherapy and targeted therapy, and AXL inhibition has been shown to enhance sensitivity to targeted therapy (19, 20). Therefore we also assessed AXL protein expression in IACS-10759 sensitive versus resistant PDXs. We found that AXL protein expression was higher in the less sensitive models by RPPA (Fig. 3B) and by Western blotting (Supplementary Fig. S6). We also assessed AXL phosphorylation at Tyrosine 702 but were unable to detect baseline phosphorylation of AXL in these samples.
To determine whether AXL plays a role in acquired resistance to IACS-10759, we treated an IACS-10759 sensitive PDX (BCX.087) with IACS-10759 until tumors progressed to greater than 100% starting tumor volume (for up to 92 days; Fig. 3C). We analyzed AXL expression in regrown BCX.087 by qRTPCR and found significant higher levels compared with the control (Fig. 3D).
These data taken together suggest that baseline molecular profiles of tumors may help separate tumors with greater reliance on OXPHOS signaling. Potential biomarkers of OXPHOS dependence include basal TNBC type and higher baseline mitochondrial gene expression while EMT signature and higher AXL expression was associated with less sensitivity to OXPHOS inhibition.
Pharmacodynamic markers of response to OXPHOS inhibition as determined by IHC
We next used a multipronged approach to determine pharmacodynamic changes induced by OXPHOS inhibition to find clinically relevant, early response markers for IACS-10759. We chose three models that are more sensitive to IACS-10759 and three models that are less sensitive to IACS-10759 for early pharmacodynamic marker assessment. We treated these PDXs with IACS-10759 for 12 days (10 treatment days) and collected the tumors four hours following final treatment. We analyzed these samples for traditional markers of apoptosis (nuclear CC3; Fig. 3E) and proliferation (PCNA, pHistone H3; Fig. 3F and G) by IHC. The most significant increase in CC3 was seen in BCX.070, which regressed on therapy. Phosphorylated Histone H3 (pHistone mitosis) decreased with IACS-10759 treatment in two of three sensitive models. For PCNA, BCX.017 (more sensitive) and BCX.084 (less sensitive) had significantly lower expression with IACS-10759 treatment.
Pharmacodynamic markers of response to OXPHOS inhibition as determined by RNA analysis
We next compared the gene expression profiles of seven PDXs, three less sensitive models and four more sensitive models, treated with IACS-10759 for 12 days by RNA-seq. On the basis of global mRNA expression profiles, sensitive models clustered separately from the less sensitive models (Supplementary Fig. S7). In addition, treated and untreated models clustered as separate groups for each sensitive model. We performed gene set enrichment analysis with these data using the Reactome and Hallmark databases. We found several profiles downregulated by IACS-10759 treatment, including a lipid metabolism profile and hypoxia profile (Supplementary Fig. S8).
We assessed the role that proliferation plays in response to OXPHOS inhibition by calculating a proliferation index for each model using RNA-seq data (18). We found that the proliferation index was significantly decreased upon IACS-10759 treatment in more sensitive models but not in the resistant models (Supplementary Fig. S9). We next removed the genes comprising the proliferation index and performed cluster analysis. Even without these proliferation specific genes, the more sensitive models and less sensitive models clustered separately (Supplementary Fig. S10).
We then analyzed the tumors using a Nanostring gene set that was specifically developed to assess genes found to be consistently altered by inhibition of OXPHOS in acute myeloid leukemia and pancreatic cancer (7). When this dataset was subjected to unsupervised clustering, the IACS-10759–treated samples for the most part clustered together, away from their controls similar to RNA-seq (Supplementary Fig. S11). The four samples that clustered away from the respective treatment groups were from the less sensitive PDXs.
When the seven PDX models were analyzed together, there were 146 genes differentially regulated with IACS-10759 treatment (Supplementary Fig. S12). Only three of four sensitive models had differentially regulated genes at FDR = 0.05 (Supplementary Fig. S13–S15). Ten genes were differentially expressed with IACS-10759 treatment in three of the four sensitive models (Supplementary Table S4). The differentially expressed genes included metabolic regulator ASNS that encodes asparagine synthetase, and TRIB3 (Tribbles pseudokinase 3) and DDIT3 (DNA damage-inducible transcript 3), genes known to play a role in cellular stress response.
Pharmacodynamic markers of response to OXPHOS inhibition as determined by RPPA
We analyzed the samples based on protein expression and phosphorylation as determined by RPPA. On unsupervised clustering, the IACS-10759–treated samples clustered together with matched untreated controls. Notably RPPA clustered IACS-10759 sensitive versus IACS-10759 PDX models separately (Supplementary. Fig. S16).
We went on to identify a set of 31 total or phosphoproteins that are differentially expressed between treated and untreated samples regardless of sensitivity (Supplementary Fig. S 17). The expression of AXL protein was not significantly changed with IACS-10759 treatment versus matched controls. Phosphorylated 5′ AMP-activated protein kinase (pAMPK), well-established marker of metabolic stress, was increased by OXPHOS inhibition. SRC phosphorylation on Tyr527, a phosphorylation site that renders the SRC enzyme less active, was decreased (21). MAPK and PI3K downstream signaling was repressed as demonstrated by the decrease in phospho-p44/42 MAPK (ERK1/2) Thr202/Tyr204) and phospho-S6 Ser235/236 and phospho-S6 Ser240/244, and phosphorylation of tumor suppressor Rb was decreased. In addition, we identified proteomic changes unique to each sensitive model (Supplementary Figs. S18–S20). Three of four sensitive models had differentially expressed genes at FDR = 0.05. However, the only proteomic change occurring in all three of these models was decreased phosphorylation of AKT Thr308, a key mediator of PI3K signaling pathway.
Identifying rational combinations with IACS-10759
Building on the finding that AXL is highly expressed in PDXs that are less reliant on OXPHOS, we hypothesized that targeting AXL may enhance sensitivity to IACS-10759. We thus tested IACS-10759 in combination with cabozantinib, a multi-kinase inhibitor that targets AXL and is approved by the FDA for the treatment of renal cell carcinoma, medullary thyroid cancer, and hepatocellular carcinoma. We tested the combination in two less sensitive models with high expression of AXL (BCX.010 and BCX.084) and one more sensitive model (BCX.087). In both less OXPHOS-dependent PDXs, we found that the combination of cabozantinib and IACS-10759 significantly prolonged tumor inhibition compared with either single agent (Fig. 4A and B). We also observed tumor regression in the more OXPHOS-sensitive model (BCX.087) compared with stable disease with single-agent treatment (Fig. 4C). We assessed signaling in the BCX.010-treated samples by RPPA and found that S6 phosphorylation was inhibited significantly more in the combination group (Fig. 4D).
We also tested IACS-10759 plus an AXL-specific inhibitor (BGB-324) in AXL high BCX.010. BGB-324 added no benefit over IACS-10759 alone (Supplementary Fig. S21A). We also assessed effect of AXL knocked down by shRNA. AXL played an important growth regulatory role, as demonstrated by the significant growth inhibition seen by AXL shRNA in BCX.010 (Supplementary Fig. S21B–S21C). However, IACS-10759 did not further enhance growth inhibitory effect.
In vivo synthetic lethality to identify therapies that sensitize to OXPHOS inhibition
To identify additional combination partners for OXPHOS inhibition, we also performed an in vivo synthetic lethality screen using a set of shRNA linked to FDA-approved agents (Fig. 5A). To do so, we transduced BCX.010 cells with molecular-barcoded shRNA at a multiplicity of infection of 0.2. We injected the cells into mice and treated the mice with vehicle or IACS-10759 for 12 days. We collected the tumors and compared remaining cells with shRNA integrations by sequencing of the barcodes. We identified several potential targets with linked FDA-approved agents with high breast cancer relevance (Fig. 5B), including CDK4, PARP1 and 2, and HDAC3.
We tested the in vitro growth inhibitory efficacy of IACS-10759 in combination with palbociclib (targeting CDK4/6), talazoparib (targeting PARP), and entinostat (targeting HDAC3) in two breast cancer cell lines, BCX.010-CL, a cell line derived from BCX-10 (Fig. 5C) and MDA-MB-468 (Fig. 5D). All three combinations were synergistic, with combination indices of 0.16, 0.05, and 0.27 for palbociclib, talazoparib, and entinostat respectively in BCX.010-CL. Using extracellular O2 consumption assay, we noted that while IACS-10759 inhibited oxygen consumption rate (OCR) in MDA-MB-468 cells, none of the other agents directly inhibited OCR, but the combinations of agent led to a greater inhibition of OCR (Supplementary Fig. S22).
Next, we tested these three combinations on BCX.010 in vivo. We found that the combination of IACS-10759 with palbociclib significantly enhanced growth inhibition at doses tested compared with palbociclib or IACS-10759 alone in BCX.010 (Fig. 6A) and an additional Rb wild-type model, BCX.080 with low (10%) ER expression (Fig. 6B). We also developed a PDX from a patient with metastatic hormone receptor–positive (HR+) breast cancer who had received palbociclib and progressed (i.e., had intrinsic resistance). We treated this PDX with palbociclib and IACS-10759 and found enhanced response to the combination (Fig 6C).
Although talazoparib plus IACS-10759 had greater growth inhibition in BCX.010 PDXs in vivo than with either agent alone, this difference was not statistically significant (Fig. 6D). To assess the role of IACS-10759 in modifying sensitivity to PARP inhibitors, we also tested the combination of talazoparib and IACS-10759 using two talazoparib acquired resistance models developed in our lab. These models were developed by treating two TNBC BL1 PDXs that regressed with talazoparib and viably collecting tumors that subsequently progressed (13). The combination with talazoparib had no improved efficacy in vivo using these models (Fig. 6E and F), but IACS-10759 alone did show statistically significant growth inhibition compared with PARP inhibition in these acquired resistance models. The combinations of entinostat and IACS-10759 did not increase efficacy over the single agent on BCX.010 in vivo (Supplementary Fig. S23).
Discussion
We show here that chemoresistant TNBC can be reliant on OXPHOS to varying extents; importantly, complete regression was observed in a PDX with IACS-01759 treatment. Thus several metabolic pathways may play a role in chemoresistant breast cancer but OXPHOS may be the absolute necessary determinant in some resistant specific TNBC. Notably, progression-free survival and overall survival in TNBC may be influenced by many clinical parameters such as tumor size and nodal status as well as molecular features. Further study is needed to see how OXPHOS activity varies with other parameters. The association of OXPHOS and outcomes in this small study is hypothesis generating and would need much larger studies for validation, also taking into account to rapidly evolving treatment algorithms. However, our finding is consistent with a recent report by Echeverria and colleagues that showed that OXPHOS was the most significantly upregulated pathway in residual tumors after TNBC PDXs were treated with chemotherapy compared with vehicle-treated tumors (10).
We show that several basal-like PDXs comprised of cancers resistant to standard neoadjuvant therapy were sensitive to IACS-10759. Moreover, we found that basal-like TNBC, which are most often defective in DNA damage repair response pathways and often sensitive to FDA-approved PARP inhibitors, are sensitive to OXPHOS inhibition. While we found in vitro synergy with talazoparib in TNBC models, we did not observe enhanced efficacy with the combination results in vivo. However, the sensitivity of models with acquired talazoparib resistance to IACS-10759 highlights that there is no cross-resistance to PARP inhibitors.
Previous studies have shown that genomic passenger alterations can result in reliance on OXPHOS. While we found molecular features associated with sensitivity to IACS-10759 in TNBC, including high mitochondrial RNA and low ENO1 protein expression among others, further studies are required to determine whether any biomarkers can reproducibly assess reliance on OXPHOS in TNBC. We show evidence that higher expression of mitochondrial protein coding RNA is correlated with functional reliance on OXPHOS in TNBC PDXs, and we believe that this demonstrates higher mitochondrial numbers in sensitive models but could also represent overcompensation due to lower mitochondrial function. We also observed an increase in phosphorylated AMPK in all models with IACS-10759 treatment. ATP, which is produced via oxidative phosphorylation, inhibits the phosphorylation of AMPK. Thus, the increased level of phosphorylated AMPK may suggests that in these models ATP deficiencies induced by OXPHOS are not fully compensated for by alternative mechanisms.
Changes in traditional IHC pharmacodynamics (CC3, pHistone H3, PCNA) did not perfectly correlate with increased sensitivity to IACS-10759. However, two of three more sensitive PDXs had significantly decreased pHistone H3, but no less sensitive PDXs had changes in pHistone H3. Interesting, the PDX that regressed with IACS-10759 treatment (BCX.070) had the largest increase in CC3. These findings suggested that both CC3 and pHistone H3, along with gene expression changes observed, have potential as pharmacodynamic markers in OXPHOS inhibitor trials as potential early molecular predictors of clinical response/benefit.
RPPA analysis revealed that IACS-10759 reduces activity of several well-characterized signaling pathways. In the set of PDXs tested, which includes a PIK3CA-mutant PDX and a PDX with PTEN loss, the PI3K/AKT pathway was inhibited by IACS-10759 treatment as evidenced by decreased phosphorylation of downstream targets such S6 ribosomal subunit, PRAS40 (pT246), and RICTOR (pT1135). Studies have suggested that reducing energy expenditure, via reduced mTOR activity, can allow cancer cells to escape energy stress caused by OXPHOS inhibition (8). Therefore, a decrease in this important protein translation pathway may be a mechanism of adaptive response/resistance to IACS-10759.
SRC has been established as a major regulator of mitochondrial energy metabolism and OXPHOS (22). On RPPA, we also observed a decrease in SRC (pT527), suggesting increased SRC activity. Several oxidative phosphorylation components have been shown to be direct targets of SRC (23). SRC has been shown to be a preeminent inhibitor of pyruvate entering oxidative phosphorylation via inhibition of PDH (24). SRC also has been shown to interact with and phosphorylate hexokinases HK1 and HK2, the rate-limiting enzymes in glycolysis (25). Thus, increased SRC activity with IACS-10759 treatment may demonstrate the cells are attempting to shift reliance on OXPHOS. Further studies are needed to determine effects of SRC modulation in combination with IACS-10759.
Our PDX testing demonstrated the potential effectiveness of inhibiting OXPHOS in TNBC with IACS-10759. However, several TNBC PDXs were less sensitive to OXPHOS inhibition alone, suggesting the need for future combination strategies with IACS-10759. We thus identified two potential combination partners that are FDA approved for other indications: cabozantinib and palbociclib.
We hypothesized that cabozantinib would synergize with IACS-10759 due to high expression of AXL in the resistant PDXs. AXL is among proteins known to increase glycolytic function in cancer and be expressed in some TNBC (26, 27). AXL is also known to promote breast cancer metastasis and chemoresistance (28, 29). AXL has been shown to upregulate two key enzymes of aerobic glycolysis (glucose transporter type 4 and pyruvate Dehydrogenase Kinase 1) via Tensin 2 (TNS2) phosphorylation (30). In melanoma, AXL expression is inversely correlated with MITF expression, and low MITF expression is linked to low OXPHOS (31). Low MITF1/AXL expression ratio is linked to resistance to multiple drug targets, strengthening the hypothesis that AXL may be associated IACS-10759 resistance. AXL suppression caused strong antitumor efficacy alone and may not properly represent what could be achievable by pharmacologic inhibition. However, AXL-specific inhibitor and AXL suppression via shRNA were not synergistic with IACS-10759. These data suggest AXL upregulation alone may not be sufficient to drive OXPHOS resistance, and may be part of a larger reprogramming such as EMT in association with drug resistance. This is supported by upregulation of other EMT markers such as TWIST in the resistant models (Supplementary Table S3). Further work is needed to assess whether AXL expression or other EMT-associated markers could potentially serve as biomarkers of resistance to OXPHOS inhibition. Notably, cabozantinib is a multi-kinase inhibitor, and its other kinase targets such as VEGFR2 and MET and/or the resultant anti-angiogenesis effect may have mediated the enhanced efficacy of the combination of cabozantinib and IACS-10759. Cabozantinib is approved for hepatocellular cancer, medullary thyroid cancer, and renal cell carcinoma and is being investigated in TNBC and multiple other tumor types. Thus, these findings may have implications for many tumor types.
Palbociclib is approved for HR+ breast cancer. In our study, we first found efficacy of the palbociclib and IACS-10759 combination in hormone receptor–negative Rb1-positive breast cancer. Other studies have shown potential utility for CDK inhibition in combination with other targets in TNBC (32–34), and this could be an avenue to explore in combination with OXPHOS inhibition. However, we also showed that PDXs generated from a patient with HR+ breast cancer were sensitive to this combination and thus this combination could be explored in the HR+ setting as well.
In summary, our study provided evidence for OXPHOS dependence in chemotherapy-resistant TNBC. Future metabolic studies are needed to better elucidate the fuel for this increased metabolic activity. We have identified and validated biomarkers that can be analyzed as a measure of OXPHOS dependence as well as correlative studies in clinical trials of OXPHOS inhibitors for assessing target engagement as well as retrospectively accessing predictive reliance. Importantly, we illuminated potential targeted therapies where OXPHOS inhibition may represent a novel combination strategy.
Authors' Disclosures
K.W. Evans reports grants from NCI during the conduct of the study. C. Tapia is an employee of Epizyme Inc. and has shares with Epizyme Inc. C. Tapia performed contract work for Armo Bioscience on another study. C. Tapia has received salary support from Merck on another study. B. Lim reports contracted research with Puma, Pfizer, Amgen, Merck, Takeda, Genentech, Calithera. D. Tripathy reports grants from Novartis and Polyphor, and personal fees from Pfizer, AstraZeneca, GlaxoSmithKline, Gilead, OncoPep, Exact Sciences, and Novartis outside the submitted work. T.A. Yap reports grants from Artios, Constellation, Cyteir, Eli Lilly, Forbius, GlaxoSmithKline, Genentech, ImmuneSensor, Ipsen, Jounce, Karyopharm, Kyowa, Novartis, Ribon Therapeutics, Regeneron, Sanofi, Scholar Rock, Tesaro, Vertex Pharmaceuticals; grants and personal fees from AstraZeneca, Bayer, Clovis, EMD Serono, F-Star, Merck, Pfizer, Repare, Seattle Genetics; personal fees from Almac, Aduro, Atrin, Axiom, Bristol Myers Squibb, Calithera, Cybrexa, Guidepoint, Ignyta, I-Mab, Jansen, Roche, Rubius, Schrodinger, Varian, and Zai Labs outside the submitted work. M.E. Di Francesco reports a patent for WO2014/031936 issued. G.F. Draetta reports personal fees from Alligator Biosciences, Biovelocita, Blueprint Medicines, Forma Therapeutics, Frontier Medicines, Helsinn, Karyopharm, Metabomed, Nurix, Orionis, Symphogen, Taiho Pharm Co, and Tessa Therapeutics outside the submitted work. P. Jones reports a patent for US9663504 issued and is an inventor on above patents, and potentially would benefit from any upside as part of UT revenue share policy. T.P. Heffernan reports other support from Boehringer Ingelheim and Taiho Pharmaceuticals; personal fees from Cullgen Inc. and Silicon Therapeutics outside the submitted work. J.R. Marszalek reports a patent for IACS-010759 issued. F. Meric-Bernstam reports personal fees from AbbVie, Aduro BioTech Inc., Alkermes, AstraZeneca, DebioPharm, eFFECTOR Therapeutics, F. Hoffman-La Roche Ltd., Genentech Inc., IBM Watson, Infinity Pharmaceuticals, Jackson Laboratory, Kolon Life Science, OrigiMed, PACT Pharma, Parexel International, Pfizer Inc., Samsung Bioepis, Seattle Genetics Inc., Tyra Biosciences, Xencor, Zymeworks, Black Diamond, Eisai, Immunomedics, Inflection Biosciences, Karyopharm Therapeutics, Mersana Therapeutics, OnCusp Therapeutics, Puma Biotechnology Inc., Seattle Genetics, Silverback Therapeutics, Spectrum Pharmaceuticals, Zentalis, grants from Aileron Therapeutics, Inc., AstraZeneca, Bayer Healthcare Pharmaceutical, Curis Inc., CytomX Therapeutics Inc., Daiichi Sankyo Co. Ltd., Debiopharm International, eFFECTOR Therapeutics, Genentech Inc., Guardant Health Inc., Klus Pharma, Takeda Pharmaceutical (formerly Millennium Pharmaceutical), Novartis, Puma Biotechnology Inc., Taiho Pharmaceutical Co., and Chugai Biopharmaceuticals, Mayo Clinic, Rutgers Cancer Institute of New Jersey outside the submitted work. No disclosures were reported by the other authors.
Authors' Contributions
K.W. Evans: Investigation, writing-original draft. E. Yuca: Investigation. S.S. Scott: Investigation. M. Zhao: Investigation. N. Paez Arango: Investigation. C.X. Cruz Pico: Investigation. T. Saridogan: Investigation. M. Shariati: Investigation. C.A. Class: Investigation. C.A. Bristow: Investigation. C.P. Vellano: Resources. X. Zheng: Data curation. A.M. Gonzalez-Angulo: Conceptualization. X. Su: Data curation. C. Tapia: Investigation. K. Chen: Data curation. A. Akcakanat: Investigation. B. Lim: Investigation, writing–review and editing. D. Tripathy: Investigation. T.A. Yap: Conceptualization, writing–review and editing. M.E. Di Francesco: Resources. G.F. Draetta: Resources. P. Jones: Resources. T.P. Heffernan: Resources. J.R. Marszalek: Conceptualization, resources. F. Meric-Bernstam: Conceptualization, resources, supervision, writing–original draft.
Acknowledgments
The authors appreciate Susanna Brisendine for helping with article preparation. They appreciate insights from our collaboration with Charles Perou. For funding, the authors acknowledge the following: NCI PDX Development and Trial Center U54 #CA224065, NIH T32 CA009599, The MD Anderson Cancer Center Sheikh Khalifa Bin Zayed Al Nahyan Institute for Personalized Cancer, MD Anderson Breast Cancer Moonshot Program, NIH Clinical Translational Science Award 1UL1TR003167, the Nellie B. Connally Breast Cancer Research Endowment, the Barr funds, Cancer Prevention & Research Institute of Texas grant RP180248, and the MD Anderson Cancer Center Support Grant (P30 CA016672).
The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).