- Open Access
PKM2 inhibition may reverse therapeutic resistance to transarterial chemoembolization in hepatocellular carcinoma
Journal of Experimental & Clinical Cancer Research volume 39, Article number: 99 (2020)
Therapeutic options for patients with hepatocellular carcinoma (HCC) are limited. Transarterial chemoembolization (TACE) is an interventional procedure used to deliver chemotherapy and embolizing agents directly to the tumor and is the procedure of choice for patients with intermediate stage HCC. While effective, more than 40% of patients do not respond to therapy, highlighting the need to investigate possible mechanisms of resistance. We sought to evaluate mechanisms of TACE resistance and evaluate a potential therapeutic target to overcome this resistance.
Using a prognostic gene signature which predicts TACE response (TACE Navigator) in a cohort of HCC patients who received TACE, patients were classified as responders and non-responders. Transcriptomic and gene pathway analysis were used to identify potential drivers of TACE resistance. Knockdown of the gene encoding rate limiting enzyme PKM2 using shRNA in HCC cell lines, as well as pharmacologic inhibition of PKM2 with shikonin using an in vitro TACE model measured response to chemotherapy under hypoxia. Finally, we replicated the TACE model with shikonin using patient derived cell line organoids (PDC). Functional studies were performed in vitro using immunoblotting, quantitative polymerase chain reaction, glycolysis and hypoxia assays.
In patient non-responders, we identified enrichment of the glycolysis pathway, specifically of the gene encoding the rate-limiting enzyme PKM2. We identified four HCC cell lines which recapitulated a TACE responder-like and non-responder-like phenotype. PKM2 knockdown in HCC cell lines demonstrated a less proliferative and aggressive phenotype as well as improved drug sensitivity to both doxorubicin and cisplatin. In vitro TACE model demonstrated that TACE non-responder-like cells overcame therapeutic resistance and rendered them susceptible to therapy through PKM2 knockdown. Lastly, we obtained similar results using a pharmacologic PKM2 inhibitor, shikonin in both cell lines, and PDC organoids.
Elevated PKM2 is associated with treatment resistance and abbreviated survival in patients receiving TACE. Elevated PKM2 in vitro is associated with increased utilization of the glycolysis pathway, resulting in oxygen independent cell metabolism. Through PKM2 knockdown as well as with pharmacologic inhibition with shikonin, non-responder cells can be reprogrammed to act as responders and could improve TACE efficacy in patients.
Hepatocellular carcinoma (HCC) is the most common primary tumor of the liver and with approximately 850,000 new cases per year, represents the fourth leading cause of cancer related death globally [1, 2]. In the United States, liver cancer is the fastest rising malignancy and while the incidence pales in comparison to other malignancies, it is among the top causes of cancer related death . Despite advances in identifying the etiology of disease, advances in durable curative therapy has not been as fruitful. One area of promise is in the treatment of intermediate stage HCC (Barcelona Clinic Liver Cancer Stage B) with transarterial chemoembolization (TACE). TACE is a procedure in which a catheter is directed into the feeding arterial branches of the tumor and chemotherapy, commonly doxorubicin, is directly infused into the tumor, followed by embolization of the feeding artery . With treatment related mortality less than 5%, TACE can be expected to result in a median survival of 11–20 months [5, 6]. While this procedure is life prolonging for many, more than 40% of patient have no objective response to treatment and are left with far less effective systemic therapies .
Due to the large number of patients that do not respond to TACE, it is imperative to identify the drivers that promote resistance. Recently, we developed a gene signature predictive of HCC patient response to TACE treatment . We showed hypoxia signaling and glycolysis/gluconeogenesis-related pathways were activated in HCC patients that did not respond to TACE. Accordingly, we hypothesized that preferential utilization of the glycolysis pathway allows HCC to evade TACE. In the present study, we demonstrate that TACE resistance is associated with alterations in glucose metabolism, specifically through the enrichment of the glycolysis related gene, Pyruvate kinase muscle isozyme M2 (PKM2). PKM2 is a splice variant of the of pyruvate kinase and its role in cancer cell metabolism has been associated with propagation of the Warburg effect, allowing for a selective growth advantage of malignant cells . Using transcriptomic data from HCC patients, we determined that TACE non-responders have elevated PKM2 expression and demonstrated that PKM2 inhibition via gene silencing or a pharmacological agent shikonin sensitized HCC cell lines and patient-derived cell lines (PDCs) to chemotherapeutic treatment in a hypoxic environment mimicking the TACE procedure. We believe that these findings may lead to a novel treatment adjunct to TACE, warranting further clinical testing.
Patient and cell line cohorts
The Liver Cancer Institute (LCI) cohort is a publicly available HCC dataset available at the Gene Expression Omnibus (GEO, GSE14520; http://www.ncbi.nlm.nih.gov/geo). Data from this cohort has been previously described . For this study, only the 105 patients who underwent TACE were selected from the LCI cohort. The Hong Kong cohort consists of 47 patients who underwent TACE, and has previously been described in detail . Patients from these cohorts were assigned as TACE responders based on the previously described TACE Navigator gene signature utilizing NanoString technology (Seattle, Washington) . NanoString was also utilized to evaluate for PKM2 expression in the Hong Kong cohort and PDCs. Gene expression as measured by NanoString counts, were log2 transformed and converted to z-score within each cohort.
HCC cell line transcriptomic data was downloaded from the Cancer Cell Line Encyclopedia (CCLE) and log2 transformed . Expression data was then z-score transformed and the TACE Navigator prognostic signature was then applied. Based on this data, cell lines were classified as responder and non-responder like and used for downstream analysis.
Cell lines and plasmid
HUH7 were cultured in Dulbecco’s modified Eagles Medium (Life Technologies) supplemented with 10% fetal bovine serum (FBS), penicillin, streptomycin and L-glutamine. Hep3B were cultured in Minimal Essential Medium (Life Technologies) supplemented with 10% FBS, penicillin, streptomycin, L-glutamine, non-essential amino acids and sodium pyruvate. SNU-387 were cultured in RPMI-1640 (Life Technologies) supplemented with 10% FBS, penicillin and streptomycin. SNU-475 were cultured in RPMI-1640 supplemented with 10% heat inactivated FBS, penicillin and streptomycin. PDCs (HCC3501, HCC 3796, HCC 4006, and HCC3258) were cultured in Dulbecco’s modified Eagles Medium: Nutrient Mixture F-12 (Life Technologies) supplemented with 10% FBS, non-essential amino acids, penicillin and streptomycin. pGFP-C-shLenti-PKM2 shRNA was purchased from Origene and packaged using FuGENE 6 transfection reagent (Promega). Cells lines were passaged less than 15X after thaw.
Immuno-blotting and quantitative RT-PCR
Protein lysates were separated on Bis-tris 4–12% gels (Life Technologies) and transferred to a nitrocellulose membrane (Life Technologies). Protein detection was performed using anti-PKM2 (Cell Signaling Technology, cat# 4053S) and anti-β-actin (Sigma-Aldrich, cat# A5316). All immunoblots are cropped for viewing. Total RNA was extracted using Trizol (Life Technologies) according to the manufacturer’s protocol. cDNA was synthesized using High Capacity cDNA Reverse Transcription Kit (Life Technologies). Gene expression was quantified by quantitative real-time transcription-polymerase chain reaction (qRT-PCR) using CFX384 Real Time System (BioRad) and SYBR Green PCR Master Mix (GenePharma). Gene expression levels were normalized to β-actin. The relative gene expression levels were detected and calculated using the ΔΔCt comparative method. The RT-PCR primer sequences are as follows: PKM2, 5′- ACTGTCCTCACCAAGTCTGG-3′ (forward) 5′-GAAGATGCCACGGTACAGGT-3′ (reverse); β-Actin, 5′-TTGTTACAGGAAGTCCCTTGCC-3′ (forward) 5′-ATGCTATCACCTCCCCTGTGT-3′ (reverse).
Glycolysis was measured using a commercially available kit (Cayman Chemical, cat#600450). Briefly, cells were seeded at 20,000 cells/well in triplicate in a 96-well plate in FBS free media appropriate for each cell line. Cells were placed in 37 °C incubator for 48 h. For hypoxia induced measurements cells were placed in a 37 °C incubator for 48 h in a closed 1% oxygen environment. After 48 h, plates were centrifuged at 1000 rpm for 5 min and 10 μl of supernatant was removed from each well and added to the reaction solution. The plate was incubated with shaking for 30 min and the absorbance at 490 nm was measured.
Cell Colony formation
Cells were seeded at 500 cells/well in a 6 well plate and cultured for 10 days in both normoxia and 1% oxygen in parallel. Cells were then washed with ice cold 1x PBS (Life Technologies) followed by ice cold fixation with methanol for 30 min. After fixation, cells were washed with dH2O followed by 2 h of 0.05% crystal violet staining and subsequent washings with dH2O. Colonies were quantified via manual counting in triplicates.
Cell proliferation and migration/invasion assays
Xcelligence (ACEA Bioscience) assays were performed for cell proliferation, cell invasion and cell migration and measured by cell index. Cell index is a measure of cell proliferation within each well of the plate as defined by (impedance at time point – impedance without cells) / nominal impedance value. Assays were performed in closed 1% hypoxic environment. Cells were seeded at 1000 cells/well in E-plates (ACEA Bioscience, cat # 5232368001), with appropriate medium in quadruplets. For TACE simulation assay, appropriate wells were treated with doxorubicin (Sigma, cat# D1515) and/or shikonin, (Sigma, cat#S7576) at the IC50 appropriate for each cell line (Supplemental Table 1). Cell specific IC50 was determined based on data from the Genomics of Drug Sensitivity in Cancer database . Proliferation was measured over 72 h. For migration and invasion assay, 30,000 cells/well were plated on a CIM plate (ACEA Bioscience, cat #5665817001) with Matrigel (for invasion only) in a 1:200 dilution with appropriate media. Migration and invasion were measured over 24 h. For cell proliferation, migration and invasion experiments were performed in triplicate.
Relative apoptosis activity was measured via Caspase-3/7 activity by using the Apo-ONE homogeneous Caspase 3/7 assay (Promega, cat#G7792). Briefly cells were seeded in cell specific media at 7500 cells/well in a 96 well plate and placed in a 37 °C incubator overnight. Cells were then placed in a hypoxic incubator for 48 h. The assay was then performed per the manufacturer’s protocol and fluorescence was measured (excitation: 485 nm, emission: 535 nm).
Cell viability and drug sensitivity
Cells were seeded in cell specific media at 7500 cells/well in a white 96 well plate and placed in a 37 °C incubator for 24 h. For doxorubicin-only experiments, two-fold drug dilutions were performed. For doxorubicin and shikonin experiments, a similar two-fold doxorubicin dilution was performed. Following dilutions, shikonin was added at the cell specific IC50. One plate was placed in a normal 37 °C incubator and the other plate was placed in hypoxia. Normoxia/hypoxia experiments were performed concurrently. After 48-h incubation, plates were removed from the incubator and allowed to reach room temperature (RT). At this point Cell-Titer Glo (Promega) was added to each well. Plates were placed on an orbital shaker for 2 min and then incubated for 10 min at RT. Luminescence was then measured. Viability assays are relative to the normoxia untreated control. For combination index experiments, a shikonin two-fold dilution was performed in addition to the previously described. To calculate the combination index, a constant ratio of 1:1 doxorubicin and shikonin was used. Viability and drug sensitivity assays were also performed with cisplatin (Milipore Sigma, cat# 232120) in the same manner. Combination index constant ratio was performed at a ratio of 4:1 cisplatin to shikonin given the relatively high IC50 of cisplatin in HCC cell lines.
Organotypic culture and viability
PDCs were cultured in AlgiMatrix 3D Culture System (Life Technologies) in the method described by Takai et al. . Briefly, one million cells were seeded and cultured in cell specific media for 2 weeks at which point they were treated with doxorubicin or shikonin alone as well as in combination and placed in hypoxia for 48 h. Haematoxylin and eosin staining was performed before and after treatment. For the collection of spheroids, the matrix was dissolved using Algimatrix Dissolving Buffer (Life Technologies) per the manufacturer’s instructions. Organoids were then pelleted by centrifugation at 300 g for 5 min. After collection, organoids were seeded in a white 96 cell plate with the addition of Cell Titer Glo 3D Cell Viability Assay (Promega, cat# G9683). Plates were placed on an orbital shaker and allowed to incubate per the manufacturer’s instructions. Luminescence was then measured, and viability was calculated relative to the untreated control.
For all statistical analyses, a p < 0.05 was considered statistically significant. Statistical tests included 2-tailed Student’s t-test for two groups or one-way ANOVA for multiple groups where appropriate. For clinical data, descriptive statistics were calculated for all variables of interest. Continuous variables were summarized using means, whereas categorical variables were summarized using frequency and percentages. Comparisons of categorical variables were performed using the chi-square or Fisher’s exact test whereas continuous variables were compared with the two-sided Student’s t-test. Overall survival (OS) was calculated based on the provided survival times utilizing the Kaplain-Meier method and log-rank test. Following a univariable analysis, a Cox proportional hazards model was constructed using all variables with p < 0.05. Genomic analyses were performed using BRB-ArrayTools version 4.6.0 (Bethesda, MD). Experimental statistics were performed using GraphPad Prism 7.01 (San Diego, CA). Clinical statistics were performed using STATA 14.0 (College Station, TX). Combination index was calculated using the Chou-Talalay method and Compusyn software [14, 15].
Alterations in glucose metabolism are correlated with TACE response
We first set out to identify signaling pathways associated with TACE resistance by utilizing the TACE Navigator gene signature, which includes 14 genes (ASNS, CDK1, DNASE1L3, FBXL5, GOT2, GRHPR, IARS, LGALS3, LHFPL2, MFGE8, MKI67, PEBP1, TNFSF10 and UBB). The TACE navigator was used to identify two distinct subpopulations, TACE responders and non-responders. Gene expression analyses between the responders and non-responders resulted in 1726 differentially expressed genes, as described previously . To identify genes that are important drivers of resistance and functionally linked to TACE, we correlated the gene expression of the 15 genes from the TACE Navigator to the 1726 differentially expressed genes using Spearman Rank. We reasoned that only genes that are functionally important for TACE resistance will be correlated with all 15 genes. We filtered out genes that were correlated with survival and correlated with all 14 genes from the TACE navigator, resulting in 105 genes (FDR < 0.05) (Fig. 1a). We then performed hierarchical clustering of TACE patients based on our 105 correlated genes and observed two distinct populations (Fig. 1b). We next performed KEGG pathway gene set enrichment analyses to identify unique signaling pathways associated with TACE response. We found correlation with 11 unique KEGG pathways, including glycolysis and gluconeogenesis (Fig. 1c). Given our previous findings that TACE resistance is associated with hypoxic reprograming , we investigated the glycolysis and gluconeogenesis pathway, which includes PKM2, Glucose-6-Phosphatase Catalytic Subunit (G6PC) and Phosphoenolpyruvate Carboxykinase 1 (PCK1) (Fig. 1d). These genes were among the most differentially expressed between responders and non-responders (Fig. 1b) and we chose to further investigate PKM2 due to its known role in hypoxia and its interaction with the hypoxia regulatory gene HIF-1α .
We next set out to examine the functional significance of each of these genes by examining their relative gene expression in responders and non-responders and whether their expression correlated with overall patient survival. We found that in a cohort of patients that received TACE, PKM2 was elevated in non-responders compared to responders (Fig. 1e; log fold change 1.14, p < 0.01). PKM2 expression using a median cutoff was also associated with significantly attenuated overall survival in this patient cohort (Fig. 1f; median survival High: 27.5 months vs Undefined, p < 0.0001). Conversely, we found that gluconeogenesis related genes, G6PC and PCK1, had significantly higher expression amongst responders, which was associated with better overall survival (PCK1: median survival High: Undefined, Low: 47.1 months, p = 0.05; G6PC: High: Undefined, Low: 47.1 months, p = 0.04) (Supplementary Figure 1). These findings indicate that PKM2 is more likely associated with non-responders and the glycolysis pathway. To ensure that these findings were not specific to a single subset of HCC patients, we set out to validate our findings in the Hong Kong cohort (n = 47). From this analysis, we found PKM2 was significantly elevated in non-responders and associated with a worse overall survival (High: 27.8 months, Low: Undefined, p = 0.05) (Supplementary Figure 2), consistent with the LCI cohort.
Finally, to establish the influence of PKM2 on survival amongst other known prognostic clinical variables, we constructed a Cox Proportional-Hazards Model (Supplemental Table 2). In the univariable analysis, elevated PKM2, BCLC stage B and C, local invasion and cirrhosis were all found to impact survival, which were included in the multivariable model. Multivariable analysis revealed high PKM2 expression was associated with a three-fold increased risk of death (HR 3.02, 95% CI: 1.50–6.04, p = 0.002). In addition, BCLC stage C was also found to be significantly associated with death but to a lesser extent (HR 2.62, 95% CI 1.30–5.27, p = 0.007). Given these findings, we hypothesized that PKM2 was associated with TACE resistance and warranted further investigation.
PKM2 and glycolysis in non-responders
To investigate TACE resistance, we first determined HCC cell lines with the similar transcriptome profiles that we observed within our patient cohorts using the 27 HCC cell lines. Using the 14 TACE Navigator genes, we performed hierarchical clustering analyses, which yielded two groups defined as responder and non-responder based on the expression patterns resembling of HCC patients (Fig. 2a). We then selected two responder-like cell lines, Hep3B and HUH7, as well as two non-responder-like cell lines, SNU-387 and SNU-475 for downstream functional studies as these lines were available in our cell repository. To confirm the responder-like and non-responder-like phenotype, we tested the IC50 of the most commonly used TACE chemotherapeutics, doxorubicin and cisplatin, in both normoxia and hypoxia, a condition mimicking the TACE procedure. We observed that in the responder-like cell lines, there was a decrease in the IC50 of doxorubicin (Fold change: Hep3B: -1.7, HUH7: − 3.45) and cisplatin (Fold change: Hep3B -1.14, HUH7–4.89) whereas in the non-responder-like cell lines, the IC50 increased for both doxorubicin (Fold change: SNU-387: 1.58, SNU-475: 1.74) and cisplatin (Fold change: SNU-387: 1.45, SNU-475: 3.27) when subjected to hypoxia (Fig. 2b). We next characterized the cell lines by determining relative PKM2 mRNA and protein expression under both normoxic and hypoxic conditions. Accordingly, the non-responder-like cell lines had higher PKM2 mRNA and protein expression in normoxia compared to responders, which was further enhanced under hypoxic conditions (Fig. 2c). To assess the functional role of PKM2 expression in glucose metabolism, glycolysis activity was assessed. The non-responder like cells have a higher baseline glycolysis activity as well as enhanced activity in hypoxia (Fig. 2d). Together, these findings are indicative of a difference in glycolysis utilization between responder and non-responder like cell lines.
PKM2 knockdown results in a less aggressive phenotype
We next investigated the functional role of PKM2 in both responder- and non-responder-like HCC cell lines. We measured lactate production in HCC cells treated with shPKM2 lentivirus and scramble control. PKM2 knockdown resulted in a significant decrease in lactate production across all cell lines (Fig. 3a). Among responder-like cells, Hep3B demonstrated a 24.2% reduction (p < 0.01) while HUH7 demonstrated a similar 20.3% reduction (p < 0.01). Conversely, non-responder-like cells experienced a greater lactate reduction (42.0% in SNU-387, 40.4% in SNU-475, p < 0.01). Apoptosis activity was also increased in all cell lines with a greater effect in the non-responder-like cells (% increase: SNU-387: 55.6%, p < 0.01; SNU-475: 53.5%, p < 0.01) than in responder-like cells (% increase: Hep3B: 11.7%, p = 0.02; HUH7: 14.0%, p = 0.02) (Fig. 3b). To test the effect of decreased glycolysis activity on cellular viability and metastatic potential, we performed colony formation, migration and invasion assays. Cell viability, as measured by colony formation, demonstrated that all cell lines were affected, with non-responder-like cells illustrating a greater effect (% decrease: SNU-387: 59.1% p = 0.05; SNU-475: 43.9%, p = 0.02) than what was observed in responder-like cells (% decrease: Hep3B: 20.8%, p < 0.01; HUH7: 20.8%, p < 0.01) (Fig. 3c). Cell migration was significantly decreased amongst all cell lines with PKM2 knockdown (Fold change: Hep3B 0.61, p < 0.01; HUH7: 0.66, p = 0.02; SNU-387: 0.73, p = 0.03; SNU-475: 0.45, p < 0.01) (Fig. 3d), whereas a reduction in PKM2 had no significant effect on invasion (Fold change: Hep3B: 0.76, p = 0.20; HUH7: 0.77, p = 0.06; SNU-387: 0.78, p = 0.07; SNU-475: 0.74, p = 0.10) (Fig. 3e). These data indicate that PKM2 knockdown disrupts glycolytic activity, thus affecting cancer-associated phenotypes. In addition, these data suggest that inhibition of PKM2 may be a useful step to overcome TACE resistance.
PKM2 knockdown sensitizes non-responder cells to doxorubicin and improves TACE-like response
We hypothesized that PKM2 knockdown could improve TACE efficacy. We assessed the effect of PKM2 knockdown on the cell’s response to doxorubicin and cisplatin, two commonly used chemotherapeutic drugs in TACE. IC50 values were calculated in hypoxia for each cell line treated with or without shPKM2 knockdown and scramble control. In all cell lines, we noted a reduction in the hypoxia IC50 of these two chemotherapeutics (Fig. 4a-b). Among responder-like cells treated with shPKM2, Hep3B demonstrated a 0.56 (p < 0.01) and 0.41 (p < 0.01) fold change in the IC50 of doxorubicin and cisplatin compared to scramble control, respectively, whereas HUH7 demonstrated a 0.32 (p < 0.01) and 0.31 (p < 0.01) fold change in the IC50 of doxorubicin and cisplatin. In addition, non-responder like cells treated with shPKM2, SNU-387 demonstrated a 0.13 (p < 0.01) and 0.18 (p < 0.01) fold change in the IC50 of doxorubicin and cisplatin while SNU-475 demonstrated a 0.07 (p < 0.01) and 0.08 (p < 0.01) fold change in the IC50 of doxorubicin and cisplatin respectively in hypoxia.
Next, we designed an in vitro TACE simulation assay to capture the effect of both chemotherapy and acute hypoxia in which cells were treated at their IC50 of doxorubicin and placed into a closed 1% oxygen system to induce acute hypoxia. Given that doxorubicin is the most commonly used TACE agent worldwide, and its slightly larger reduction in IC50 in PKM2 knockdown compared to cisplatin, we elected to perform further characterization of in vitro TACE model using doxorubicin. Cellular proliferation was then measured over 72 h. As expected, the responder cell lines had a significant response to doxorubicin alone (Fold change: Hep3B: 0.22, p < 0.01; HUH7: 0.01, p < 0.01), which was slightly improved when cells were treated with shPKM2 (Fold change: Hep3B: 0.09, p < 0.01; HUH7: 0.05, p < 0.01) (Fig. 4c Left panels). Conversely, amongst the non-responder cell lines, SNU387 had no response to doxorubicin (Fold Change:1.24 p = 0.33) whereas SNU475 had a minimal response (Fold Change: 0.87, p < 0.01) (Fig. 4c, right panels). When non-responder cells were treated with shPKM2 and doxorubicin, both non-responder cell lines had a drastic reduction of cell growth (Fold Change: SNU-387: 0.03, p < 0.01); SNU-475: 0.12, p < 0.01) (Fig. 4c, Right panels). These data suggest that PKM2 inhibition was able to inhibit cell growth in all cell lines, resulting in improved TACE-like efficacy in vitro (Fig. 4d).
Pharmacologic inhibition of PKM2 sensitizes non-responder cells to doxorubicin and improves TACE response in vitro
Shikonin, a naphthoquinone, has been shown to inhibit glycolysis through PKM2 specific inhibition . We next sought to evaluate whether shikonin had similar inhibitory effects in our TACE response cell models. We first tested the effect of shikonin on the sensitization of cancer cells to either doxorubicin (Fig. 5a) or cisplatin (Fig. 5c). In the responder-like cell lines, Hep3B demonstrated a reduction in the area under the curve (AUC) in hypoxia of 31% (p < 0.01) with doxorubicin and 31% (p < 0.01) with cisplatin whereas HUH7 demonstrated a 53% reduction (p < 0.01) and a 5% reduction (p = 0.02) with doxorubicin and cisplatin, respectively. The non-responder-like cell lines also experienced a decrease in the AUC with combination chemotherapy and shikonin. SNU-387 demonstrated an 82% (p < 0.01) and 85% (p < 0.01) reduction in AUC with doxorubicin and cisplatin, respectively, whereas SNU-475 demonstrated a 94% (p < 0.01) reduction in AUC with doxorubicin and a 95% (p < 0.01) AUC reduction with cisplatin. To determine the interaction between both doxorubicin and shikonin and cisplatin and shikonin in the non-responder like cell lines, we calculated the combination index via the Chou-Talalay method (Fig. 5b, d). We observed that doxorubicin and shikonin are synergistic at higher level of cytotoxicity (Fig. 5b), whereas the combination of shikonin and cisplatin interact in an additive manner (Fig. 5d). These data indicate that the combination therapy may be a viable option to improve TACE.
Next, we performed in vitro TACE simulation with the addition of shikonin and doxorubicin. Consistent with the above observations, responder-like cells had a significant decrease in cellular proliferation with doxorubicin alone (Fold change: Hep3B: 0.27, p < 0.01; HUH7: 0.31, p < 0.01), which was further enhanced with the combination therapy of shikonin and doxorubicin (Fold change: Hep3B: 0.15, p < 0.01; HUH7:0.10, p < 0.01) (Fig. 5e, Left panels). Conversely, doxorubicin alone had a small reduction in cellular proliferation in both non-responder like cell lines (Fold change: SNU-387: 0.64, p < 0.01; SNU-475: 0.81, p = 0.04) but with the addition of shikonin therapeutic efficacy was significantly increased which resulted in a 0.11 (p < 0.01) fold change in the SNU-387 cells and a 0.13 (p < 0.01) fold change in the SNU-475 cells (Fig. 5e, Right panels). These data suggest that resistance can be overcome with the addition of a PKM2 inhibitor during the TACE procedure, which has the potential to greatly improve treatment response. Additionally, we observed that combination therapy with doxorubicin and shikonin significantly improved in vitro TACE response in all cell lines (Fig. 5f). This data may suggest that combination shikonin and doxorubicin has the potential to improve TACE efficacy even in those patients with expected favorable outcomes.
Patient derived non-responder like organoids can undergo therapeutic reprogramming
To test the hypothesis that combination shikonin and doxorubicin improves TACE response in pre-clinical models, we used a 3D model utilizing PDCs to better approximate a tumor-like environment. Given our findings that PKM2 enrichment both confers a negative prognostic outcome as well as is associated with TACE resistance, we elected to screen 33 HCC PDCs for PKM2 expression. Utilizing the NanoString platform, median relative PKM2 mRNA expression was 13.33 (Fig. 6a). We then selected two representative lines from both the PKM2-high and PKM2-low groups and generated 3D organoids. In a similar manner to our 2D TACE in vitro assay, organoids were treated at cell line specific IC50s of both shikonin and doxorubicin and subjected to hypoxia. Organoid viability was then assessed. PKM2 low cells behaved in a similar manner to the responder-like cells in 2D culture with both cell lines demonstrating a statistically significant decrease in viability upon doxorubicin treatment alone (Viability fold change: HCC 3796: 0.64, p < 0.01; HCC 4006: 0.67, p = 0.03) (Fig. 6b, Left panels). Conversely, PKM2 high cells behaved in a similar manner when compared to the TACE non-responder cells. Both HCC 3258 and HCC 3501 demonstrated resistance to doxorubicin treatment alone (Viability fold change: HCC 3258: 0.74, p = 0.37; HCC 3501: 0.91, p = 0.34) whereas when treated with the combination of doxorubicin and shikonin, resistance is overcome resulting in a 0.27 (p = 0.01) and 0.49 (p < 0.01) fold change in HCC 3258 and HCC 3501 respectively (Fig. 6b, Right panels). In addition to the improved response in non-responders, we also noted a significant difference in the response to doxorubicin and combination treatment regardless of PKM2 expression (Fig. 6c). Finally, we examined the histologic characteristics of the PKM2 High organoids both before and after combination therapy. In a pretreatment condition, HCC spheres displayed a well circumscribed structure whereas following combination treatment they appeared more disorganized and irregular. Similarly, the results seen in our 2D model, we observed increased efficacy of our in vitro TACE assay when non-responder like organoids were treated with combination PKM2 inhibitor and standard chemotherapy.
Unlike many other solid organ tumors where advances in detection and treatment have resulted in improved survival, HCC deaths among both men and woman continues to rise [3, 18]. In patients with advanced disease, the only option was Sorafenib, a multi tyrosine kinase inhibitor, which in the landmark SHARP trial conveyed a less than three-month survival advantage compared to supportive care . The recent introduction of nivolumab and regorafenib in the second line and lenvatinib in the first line offer only modest benefit [20,21,22]. Where systemic therapy has failed, locoregional therapies such as TACE for intermediate HCC have had some success. As previously stated, while benefit can be seen, patient selection is paramount, with more than 40% of patient having no response to therapy. In addition, the definition of intermediate HCC comprises only 10–12% of patients at the time of diagnosis, further limiting treatment . This highlights the need for both improvement in the efficacy of TACE as well as improved methods of patients’ selection to better identify who will benefit from treatment. Our approach to improve TACE efficacy is focused on overcoming the mechanism of resistance. We demonstrate that through manipulation of cellular biology we are able to greatly increase chemotherapeutic sensitivity.
Since first described in the 1920’s by Otto Warburg, the propensity for cancer cells to preferentially utilize glycolysis even in the presence of oxygen has been well documented [24, 25]. While still a topic of much debate, the notion that cancer cells undergo metabolic reprogramming has been dubbed an “emerging hallmark” of cancer . From our data we demonstrate that glucose metabolism dysregulation is heterogenous amongst different patient populations. We also demonstrate that the degree of glycolysis dysregulation is associated with TACE resistance. PKM2 represents the more common of the two splice variants of PKM gene and is abundant in tumors . PKM2 expression has been shown to be elevated in many cancer types and specifically in HCC has be shown to portray a negative prognosis . Our findings show that within the specific subset of patients treated with TACE, elevated PKM2 expression is associated with attenuated survival. In addition, we have demonstrated experimentally that increased expression of PKM2 is associated with increased glycolysis activity. With the knowledge that glycolysis enrichment is related to TACE resistance and PKM2 expression is associated with glycolysis activity, these findings further support the role of PKM2 as a contributor to TACE resistance. In addition, we have observed that all cell lines experience a change in cancer associated phenotypes when PKM2 was knocked down, indicating that this gene may be a suitable candidate to target. In fact, in SNU-387 cells, PKM2 knockdown revealed a relative change in IC50 of both doxorubicin and cisplatin that was more dramatic than either of the responder-like cell lines. A possible explanation for this response may be related to the relative expression of PKM2 in responder-like compared to non-responder-like cell lines, allowing for a differential effect with PKM2 knockdown, however, further mechanisms remain be elucidated.
It has been well established that traditional systemic chemotherapy is ineffective in HCC . Conversely, when administered directly into hepatic arterial circulation, it does appear to convey a survival advantage as compared to embolization alone . Therefore, there is a need to improve the response to TACE chemotherapy, as it appears that hypoxia alone does not induce a sufficient enough tumor response to result in prolonged survival. In the United States, TACE is most commonly performed with a single agent chemotherapy, with doxorubicin and cisplatin representing the two most common agents . Current efforts to improve TACE have focused on the reformulation of doxorubicin as a drug eluting bead (DEB-TACE) with the hypothesis that a more controlled doxorubicin release would improve response rates. While superior to bland embolization, the Precision V trial of DEB-TACE to conventional TACE failed to show improved response [29, 30]. Here, we have demonstrated that PKM2 inhibition resulted in improved doxorubicin and cisplatin sensitivity in HCC cell lines. In addition, with an in vitro TACE assay we demonstrated that previously resistant cell lines become sensitized to therapy. Therefore, we believe that the inhibition of PKM2 inhibition may not only improve response rate but may also improve the magnitude of response in select patients.
As with any proposed targeted therapy, a drug should ideally be specific, minimizing off-target effects while maintaining an acceptable burden of adverse events. While RNA interference has proven to be a successful tool in identification and validation of potential lethal targets in cancer research, its application in clinical practice is limited [31, 32]. We therefore set out to demonstrate that a pharmacologic inhibitor of PKM2 could produce similar results. Currently there are three available PKM2 inhibitors, the small molecule Compound 3, a naphthoquinone, shikonin and its analog alkannin . While shikonin’s anti-cancer properties has been studied for some time, Chen et al. identified its mechanism as an inhibitor of glycolysis by targeting PKM2 . Additionally, shikonin and its analog were observed to have a much greater effect on PKM2 activity than Compound 3. We demonstrate that doxorubicin in combination with shikonin had a greater effect on the proliferation of cell lines and viability of PDC organoids than either treatment alone when used together in a TACE simulation and thus this combination represents a reasonable therapeutic approach to test in clinical trials. Lastly, given that we opted to characterize our in vitro TACE model with doxorubicin, it is possible that using different chemotherapy agents will produce differential effects in combination with PKM2 inhibition. Future work will aim at further characterizing this model using a larger selection of chemotherapy agents, as well as combination therapy to identify optimal agents for combination with PKM2 inhibition.
We demonstrate a potential strategy by which the TACE procedure can be augmented through targeted therapeutic reprogramming through the inhibition of PKM2. We have established that enrichment of PKM2 is a poor prognostic indicator and represents a pivotal role in the altered glucose metabolism observed in non-responders. Furthermore, we have demonstrated that the inhibition of PKM2 results in decreased cell viability and improved response to commonly used TACE chemotherapeutics. More importantly, we demonstrated that when PKM2 is inhibited, treatment resistance can be overcome, and cell death can be induced. Finally, we showed that with a known pharmacologic inhibitor, shikonin, TACE efficacy in vitro is greatly improved. We believe that these findings lay the groundwork for future clinical trials, with the potential to improve patient outcomes.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Barcelona clinic liver cancer stage
Patient derived cell line(s)
Pyruvate kinase muscle isoenzyme M2
Liver cancer institute
Cancer cell line encyclopedia
Fetal bovine serum
quantitative rea-time transcription-polymerase chain reaction
Glucose-6-Phosphatase catalytic subunit
Phosphoenolpyruvate carboxykinase 1
Area under the curve
Drug eluting bead TACE
Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.
Villanueva A. Hepatocellular Carcinoma. N Engl J Med. 2019;380(15):1450–62.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69(1):7–34.
Llovet JM, Bruix J. Systematic review of randomized trials for unresectable hepatocellular carcinoma: chemoembolization improves survival. Hepatology. 2003;37(2):429–42.
Marelli L, Stigliano R, Triantos C, Senzolo M, Cholongitas E, Davies N, et al. Transarterial therapy for hepatocellular carcinoma: which technique is more effective? A systematic review of cohort and randomized studies. Cardiovasc Intervent Radiol. 2007;30(1):6–25.
Llovet JM, Di Bisceglie AM, Bruix J, Kramer BS, Lencioni R, Zhu AX, et al. Design and endpoints of clinical trials in hepatocellular carcinoma. J Natl Cancer Inst. 2008;100(10):698–711.
Gillmore R, Stuart S, Kirkwood A, Hameeduddin A, Woodward N, Burroughs AK, et al. EASL and mRECIST responses are independent prognostic factors for survival in hepatocellular cancer patients treated with transarterial embolization. J Hepatol. 2011;55(6):1309–16.
Fako V, Martin SP, Pomyen Y, Budhu A, Chaisaingmongkol J, Franck S, et al. Gene signature predictive of hepatocellular carcinoma patient response to transarterial chemoembolization. Int J Biol Sci. 2019;15(12):2654–63.
Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, et al. The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth. Nature. 2008;452(7184):230–3.
Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, et al. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Cancer Res. 2010;70(24):10202–12.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483(7391):603–7.
Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of drug sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res. 2013;41(Database issue):D955–61.
Takai A, Fako V, Dang H, Forgues M, Yu Z, Budhu A, et al. Three-dimensional Organotypic culture models of human hepatocellular carcinoma. Sci Rep. 2016;6:21174.
Chou TC. Drug combination studies and their synergy quantification using the Chou-Talalay method. Cancer Res. 2010;70(2):440–6.
Chou TC. Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacol Rev. 2006;58(3):621–81.
Luo W, Semenza GL. Emerging roles of PKM2 in cell metabolism and cancer progression. Trends Endocrinol Metab. 2012;23(11):560–6.
Chen J, Xie J, Jiang Z, Wang B, Wang Y, Hu X. Shikonin and its analogs inhibit cancer cell glycolysis by targeting tumor pyruvate kinase-M2. Oncogene. 2011;30(42):4297–306.
Cronin KA, Lake AJ, Scott S, Sherman RL, Noone AM, Howlader N, et al. Annual report to the nation on the status of Cancer, part I: national cancer statistics. Cancer. 2018;124(13):2785–800.
Llovet JM, Ricci S, Mazzaferro V, Hilgard P, Gane E, Blanc JF, et al. Sorafenib in advanced hepatocellular carcinoma. N Engl J Med. 2008;359(4):378–90.
Bruix J, Qin S, Merle P, Granito A, Huang YH, Bodoky G, et al. Investigators RESOURCE, Regorafenib for patients with hepatocellular carcinoma who progressed on sorafenib treatment (RESORCE): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet. 2017;389(10064):56–66.
El-Khoueiry AB, Sangro B, Yau T, Crocenzi TS, Kudo M, Hsu C, et al. Nivolumab in patients with advanced hepatocellular carcinoma (CheckMate 040): an open-label, non-comparative, phase 1/2 dose escalation and expansion trial. Lancet. 2017;389(10088):2492–502.
Kudo M, Finn RS, Qin S, Han KH, Ikeda K, Piscaglia F, et al. Lenvatinib versus sorafenib in first-line treatment of patients with unresectable hepatocellular carcinoma: a randomised phase 3 non-inferiority trial. Lancet. 2018;391(10126):1163–73.
Llovet JM, Real MI, Montana X, Planas R, Coll S, Aponte J, et al. Arterial embolisation or chemoembolisation versus symptomatic treatment in patients with unresectable hepatocellular carcinoma: a randomised controlled trial. Lancet. 2002;359(9319):1734–9.
Warburg O. The metabolism of carcinoma cells. J Cancer Res. 1925;9(1):148–63.
Denko NC. Hypoxia, HIF1 and glucose metabolism in the solid tumour. Nat Rev Cancer. 2008;8(9):705–13.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Lu DH, Lv WW, Li WX, Gao YD. High PKM2 expression is independently correlated with decreased overall survival in hepatocellular carcinoma. Oncol Lett. 2018;16(3):3603–10.
Zhu AX. Systemic therapy of advanced hepatocellular carcinoma: how hopeful should we be? Oncologist. 2006;11(7):790–800.
Malagari K, Pomoni M, Kelekis A, Pomoni A, Dourakis S, Spyridopoulos T, et al. Prospective randomized comparison of chemoembolization with doxorubicin-eluting beads and bland embolization with BeadBlock for hepatocellular carcinoma. Cardiovasc Intervent Radiol. 2010;33(3):541–51.
Lammer J, Malagari K, Vogl T, Pilleul F, Denys A, Watkinson A, et al. Prospective randomized study of doxorubicin-eluting-bead embolization in the treatment of hepatocellular carcinoma: results of the PRECISION V study. Cardiovasc Intervent Radiol. 2010;33(1):41–52.
Wang Z, Rao DD, Senzer N, Nemunaitis J. RNA interference and Cancer therapy. Pharm Res. 2011;28(12):2983–95.
Takai A, Dang H, Oishi N, Khatib S, Martin SP, Dominguez DA, et al. Genome-wide RNAi Screen identifies PMPCB as a therapeutic vulnerability in EpCAM+ hepatocellular carcinoma. Cancer Res. 2019;79(9):2379–91.
Marshonna Forgues for laboratory support.
This work was supported by grants (ZIA BC 010877, ZIA BC 010876, ZIA BC 010313, and ZIA BC 011870) from the intramural research program of the Center for Cancer Research, National Cancer Institute of the United States.
Ethics approval and consent to participate
This study utilized deidentified publicly available human data and is IRB exempt.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Table S1. IC50 Concentrations Used For Invitro Experiments. Table S2. Univariable and Multivariable Cox Proportional Hazards Model of LCI cohort. Figure S1. Gene expression and survival analysis of gluconeogenesis related genes correlated with TACE navigator. Figure S2. Gene expression and survival analysis of PKM2 in Hong Kong cohort. Figure S3. Dose effect curves for patient derived cell lines.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Martin, S.P., Fako, V., Dang, H. et al. PKM2 inhibition may reverse therapeutic resistance to transarterial chemoembolization in hepatocellular carcinoma. J Exp Clin Cancer Res 39, 99 (2020). https://doi.org/10.1186/s13046-020-01605-y
- Hepatocellular carcinoma
- Therapeutic reprogramming