- Open Access
Metabolic impairment of non-small cell lung cancers by mitochondrial HSPD1 targeting
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 248 (2021)
The identification of novel targets is of paramount importance to develop more effective drugs and improve the treatment of non-small cell lung cancer (NSCLC), the leading cause of cancer-related deaths worldwide. Since cells alter their metabolic rewiring during tumorigenesis and along cancer progression, targeting key metabolic players and metabolism-associated proteins represents a valuable approach with a high therapeutic potential. Metabolic fitness relies on the functionality of heat shock proteins (HSPs), molecular chaperones that facilitate the correct folding of metabolism enzymes and their assembly in macromolecular structures.
Gene fitness was determined by bioinformatics analysis from available datasets from genetic screenings. HSPD1 expression was evaluated by immunohistochemistry from formalin-fixed paraffin-embedded tissues from NSCLC patients. Real-time proliferation assays with and without cytotoxicity reagents, colony formation assays and cell cycle analyses were used to monitor growth and drug sensitivity of different NSCLC cells in vitro. In vivo growth was monitored with subcutaneous injections in immune-deficient mice. Cell metabolic activity was analyzed through extracellular metabolic flux analysis. Specific knockouts were introduced by CRISPR/Cas9.
We show heat shock protein family D member 1 (HSPD1 or HSP60) as a survival gene ubiquitously expressed in NSCLC and associated with poor patients’ prognosis. HSPD1 knockdown or its chemical disruption by the small molecule KHS101 induces a drastic breakdown of oxidative phosphorylation, and suppresses cell proliferation both in vitro and in vivo. By combining drug profiling with transcriptomics and through a whole-genome CRISPR/Cas9 screen, we demonstrate that HSPD1-targeted anti-cancer effects are dependent on oxidative phosphorylation and validated molecular determinants of KHS101 sensitivity, in particular, the creatine-transporter SLC6A8 and the subunit of the cytochrome c oxidase complex COX5B.
These results highlight mitochondrial metabolism as an attractive target and HSPD1 as a potential theranostic marker for developing therapies to combat NSCLC.
Lung cancer is the most commonly diagnosed cancer worldwide and the main cause of cancer-related death in both men and women . Non-small cell lung cancer (NSCLC) is the most frequent type occurring in about 85% of cases , mainly presenting with an adenocarcinoma (LUAD) or squamous cell lung cancer (SqCC) histology . The standard of care for NSCLC has considerably changed in the last decade with the introduction of target therapies and immunotherapies. These, in synergy with chemo/radiotherapy, allowed to improve the clinical outcome of a fraction of NSCLC patients, in some cases providing a durable disease stabilization . Nevertheless, relapse and occurrence of treatment-resistant phenotypes are still very common, contributing to the current suboptimal scenario of 60 to 70% overall survival for early-stage and between 0 to 10% for patients carrying the advanced disease . This unsatisfactory picture strongly reinforces the need for developing new drugs targeting other fundamental lethal properties of lung cancer cells.
It has long been established that cancer cells alter their metabolism to sustain their increased energetic demands and fuel the malignant phenotype . In contrast, or rather as an expansion of the so-called “Warburg effect”, it is now clear that cancers not only depend on glycolysis for their growth, but also rely on mitochondrial respiration , enabling a hybrid or “plastic” behavior in which an opportunistic energetic adaptability determines or facilitate survival and aggressiveness . However, the metabolic alterations frequently found in tumor cells could be turned into therapeutic opportunities, as the newly acquired cancer-specific vulnerabilities or dependencies have been found to represent, in some cases, excellent drug targets [9, 10]. Recent pivotal studies have shown the power of this approach to induce a significant tumor growth arrest [11,12,13].
Many key cellular metabolic processes are dependent on the functionality of heat shock proteins (HSPs) . HSPs are highly conserved ATP-dependent molecular chaperons mainly involved in maintaining the correct protein structure, for instance avoiding mis-folding or aggregation and protecting cell integrity under stressful conditions . High levels of HSPs have been found in different cancer cells, where they sustain higher metabolic demand and correct misfolded oncoproteins . The mitochondrial heat shock protein family D member 1 (HSPD1), also called HSP60, belongs to this family  and together with the co-chaperonin HSP10 plays an essential role in mitochondrial-imported proteins folding or refolding under mitochondrial stress . HSPD1 is involved in several diseases such as neurodegenerative disorders or cardiovascular diseases . In addition, HSPD1 carries a central role in cancer development, with either pro-survival or pro-apoptotic functions reported in different tumor types . In NSCLC, HSPD1 has been proposed as a biomarker, but little is known about its role or the effects of its suppression [21, 22].
Here, we investigated the tumor-promoting role of HSPD1 and explored if its targeting could interfere with the malignant metabolic mechanisms underlying NSCLC. Our results revealed HSPD1 as a ubiquitous lung cancer dependency gene, which maintains the metabolic fitness and promotes growth/survival, thus highlighting it as an attractive therapeutic target for NSCLC.
Cell culture and chemicals
All cell lines were cultured in media supplemented with 10% FBS, 1% Pen/Strep and 1% L-Glutamine (all from Sigma) at 37 °C and 5% CO2 in a humidified incubator. H460, H1299, H520, SK-MES-1 and Calu-1 (from ATCC) were cultured in RMPI-1640 (Sigma). H23 and H838 (both from ATCC) cells were cultured in RMPI-1640, supplemented with 1 mM sodium pyruvate (Gibco). A549, 293T, BEAS-2B (ATCC), BEN (DSMZ), murine Lewis lung carcinoma cell line LL2, Ladi2.1 and Ladi3.1 cells (derived from p53fl/fl-LSL KRASG12D/ + mouse NSCLC model) were cultured in DMEM (Sigma). Human cells were STR-profiled, used between passages 3 and 15, examined for mycoplasma regularly (detection kit from Invivogen). The A549-Cas9 cell line was generated by lentiviral transduction with lentiCas9-Blast (Addgene #52962) and selection with 100 μg/mL of blasticidin (Sigma) for 3 days. KHS101 hydrochloride (4888) was purchased from Tocris; cisplatin (CAS 15663–27-1), Z-VAD-FKM (CAS 187389–52-2) and 2-Deoxy-D-Glucose (2-DG, sc-202010) were purchased from Santa Cruz Biotechnology; Necrostatin-1 (N9037-10NMG), Ferrostatin-1 (SML0583-5MG) and 1-Fluoro-2,4-dinitrobenzene (DNFB, D1529) were purchased from Sigma. HB072 was obtained as previously shown .
Proliferation assay and in vitro drug treatment
For in vitro drug treatment, cells were plated in a 96-well plate in low density (5–10% initial confluence) and incubated overnight. On the next day, cells were treated either with vehicle or with the drugs. Plates were loaded into the IncuCyte Zoom (Essen Bioscience) and scanned every 2–4 h. For each scan, phase contrast images were acquired from every well and analyzed by IncuCyte Zoom software. IC50 values were calculated after 72 h of treatment as expression of confluency percentage of treated cells normalized to vehicle control (100% confluency) using Graphpad (nonlinear regression curve fitting model) to obtain the values. Proliferation assay was performed as described in . Growth curves were analyzed with IncuCyte Zoom software.
Cytotoxicity death assay
For cytotoxicity death assay, cells were seeded in 96-well plates in low density (5,000 cells/well) and incubated overnight. 1000X Cytotox Green Reagent or Caspase 3/7 Green Apoptosis Assay Reagent (Essen Bioscience) was diluted in medium and working dilutions of the drug were prepared in Cytotox Green or Caspase3/7 Green supplemented media. After treatment, plates were loaded in IncuCyte Zoom and images were acquired in real-time for phase to quantify growth. Activity of green reagent was simultaneously acquired at the green channel to quantify death. IncuCyte Zoom software was used for the analysis and data export.
Extracellular flux assays
Oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) measurements were determined using the XFe96 Extracellular Flux Analyzer (Seahorse Bioscience/Agilent Technologies). Cells were seeded in specialized culture microplates at a density of 10,000–20,000 cells/well and prepared as describe in . For Mito Stress Test (Agilent kit 103015–100) 1 μM oligomycin, 1 μM FCCP, 1 μM rotenone and 1 μM antimycin were sequentially injected at regular intervals. OCR, indicator for mitochondrial respiration, and ECAR, indicative of glycolysis, were measured.
Plasmids for HSPD1 knock down (TRCN0000029444, TRCN0000029445, TRCN0000029446, TRCN0000029447 and TRCN0000029448 for human cell lines) were purchased from Sigma. Scrambled pLKO.1 (referred to as pLKO) was used as non-targeting control. NLRC5 expression vector (EX-A3335-Lv105) and control vector (Ex-Neg-LV105b) were from GeneCopoeia. Plasmids containing RNAs for COX5B (gRNA sequence #1: ACTTCGCGGAGCTGGAACGC; gRNA sequence #2: CAGCCAAAACCAGACGACGC) and SLC6A8 (gRNA sequence #2: ACGGGGCCGTCGCCCTTGGC; gRNA sequence #3: GCCCTTACCATGCAGACCAG) knock-out (pLENTI-CRISPR-V2) were from Genscript. For production of lentiviral particles, 293 T cells were transfected with 8 μg knock-down/expression vectors and 2 μg of packaging vectors (pMDL, pVsVg and pRevRes) in complex with 24 μg PEI (Polysciences) in 0.9% NaCl. After 48 h, supernatant was collected, centrifuged and filtered. For transduction, 150,000 cells were seeded in a 6-well plate and infected in presence of 8 μg/mL polybrene (Sigma). Selection was done with 3 μg/mL puromycin (Sigma) and cells were maintained in 1 μg/mL puromycin.
Western blot analysis
For proteins isolation, cells were lysed in RIPA buffer and quantified using Pierce BCA kit (Thermo-Fisher). For cytosolic/mitochondrial fractionation Cytochrome c Release Assay Kit (abcam, ab65311) was used. Proteins lysates (10–20 μg) were resolved on 5%-12% SDS–PAGE gels and transferred to PVDF membrane (Thermo-Fisher). Membranes were blocked in 5% Milk (BioRad) or 5% BSA (Sigma) in 1X TBST and incubated overnight at 4 °C in primary antibodies. Membranes were then washed with 1X TBST and incubated with secondary antibodies (Southern Biotech) for 1 h. The membranes were developed with ECL reagent (Thermo Fisher) on to X-ray films (Thermo-Fisher) using the chemiluminescence imager, AGFA CP100. Rabbit anti-HSPD1 (ab46798, 1:20,000) and rabbit anti-NLRC5 (ab117624, 1:1000) antibodies were purchased from Abcam; mouse anti-TOMM20 (H00009804-M01, 1:1000) was purchased from Abnova; anti-β-Actin (8H10D10) HRP conjugate (1:10,000) was purchased from Cell Signaling. Protein band density was quantified using ImageJ.
Colony formation assay
To assess clonogenic ability of knockdown cells or KHS101 treated cells, a colony formation assay was performed. Control (pLKO) and knockdown cells were plated in triplicates at a low density (1,500 cells/well) in a 6-well plate and they were allowed to grow until they formed visible colonies. To assess the clonogenic ability of cells upon KHS101 treatment, 200,000 cells were seeded in triplicates in 6-well plates, subsequently treated with 10 μM KHS101 for 5 days, after which media was replaced with normal media, whereas the DMSO cells were collected, counted and seeded in triplicates at low density (1,500 cells/well). Cells were allowed to grow until they formed visible colonies in the controls. Cells were washed with 1X PBS and fixed in 10% (v/v) formalin (Sigma) for 5 min. After washing, colonies were stained with 0.05% crystal violet solution (Biomatik, CAS 548629) for 30 min and washed twice with deionized water. Colonies were photographed and counted.
Cell cycle analysis was performed using the Propidium iodide (PI; Sigma) staining. 50,000 cells were seeded in triplicates in 12-well plates and incubated overnight. Samples were prepared as described in . For ROS staining, 50,000 cells were seeded in 12-well plates and the day after they were treated with KHS101 in combination with cisplatin. After 24 h, samples were washed and stained for 30 min at 37 °C with CM-H2DCFDA 5 μM (Thermo Fisher, C6827). Then they were collected in FACS tubes, washed and re-suspended in FACS buffer (2% FBS, 5 mM EDTA in 1X PBS). Samples were run on Cytoflex FACS machine (Beckman) and data were analyzed using FlowJo software v10.6.
HSP60 (HSPD1) immunohistochemistry was performed on n = 30 NSCLC lung surgical specimens using an automated platform (BenchMark, Ventana Medical Systems, Roche, Basel, Switzerland). Briefly, samples were pre-treated for 36 min with antigen retrieval ULTRA CC1 then they were incubated for 32 min at 36 °C with HSP60 (HSPD1) (AbTA800758, Clone OTI3A2, 1:100 dilution, ORIGENE, Rockville, US) primary antibody.
For mitochondrial staining, 100,000 cells were seeded in 12-well plates. The day after, MitoTracker Green FM (Cell Signaling Technologies, #9074S) was diluted to 400 nM directly into fresh media and added to the cells, which were incubated for 30 min at 37 °C. After incubation, media was changed, and live cells were imaged using IncuCyte-Zoom.
Total RNA was extracted using miRNeasy Mini Kit (Qiagen) according to the manufacturer's instruction. 500 ng of RNA were used to perform reverse transcription with the Tetro cDNA synthesis Kit (Bioline). cDNA was amplified with gene-specific (HK2, FBP1 and PKM) TaqMan probes (Applied Biosystem) using TaqMan Universal Master Mix II (Thermo Fisher Scientific). GAPDH was used as internal control. Quantifications were done by Applied Biosystems 7300 Real Time PCR system and relative mRNA expression level was calculated using the ΔΔCt method.
CL-100 ProLiFiler screening
26 different lung cancer cell lines (Additional file 1: Table S1) were screened to determine their viability upon treatment with KHS101. The drug was added to the cells one day after seeding and the treatment was performed by nanodrop-dispensing using a Tecan Dispenser. 0.1% DMSO (solvent) and Staurosporine (1,0E-05 M) served as high control (100% viability) and low control (0% viability), respectively. After 72 h of incubation with the compound, cell plates were equilibrated to room temperature for one hour and CellTiter-Glo™ Luminescence Cell Viability Reagent (Promega) was added to the cell suspension. The luminescence was measured 1 h later using a luminometer. IC50 values were expressed as percentage of proliferation in presence of solvent alone (100% = high control) as compared to cells treated with 1E-5 M Staurosporine (0% = low control). IC50 calculation was performed using GraphPad Prism software with a variable slope sigmoidal response fitting model.
Generation of CRISPR/Cas9 lentiviral library and CRISPR screen
The CRISPR/Cas9 lentival library was generated as described in . A549 Cas9 cells were transduced with serial dilutions of a virus to find the MOI of ~ 0.3. A549 Cas9 cells were transduced with lentiviral Human GeCKO v2 knockout pooled library part A and part B (a gift from Feng Zhang , Addgene # 1000000049) at MOI of 0.3 in the presence of 10 μg/mL of polybrene for 24 h, then replaced the virus medium with fresh growth medium and continued to culture the cells for 48 h. The cells were selected with 50 μg/mL of puromycin for 3 days. Subsequently, 4 million cells were seeded in 10-cm dishes and treated the day after either with vehicle or with 15/10 μM KHS101. A549 Cas9 pLKO were used as positive control. After 2 weeks, the media containing the drug was replaced with normal media to allow the cells to grow again. 40 million cells were then collected for genomic DNA isolation. Next generation sequencing was performed on the Illumina HiSeq 2500 platform in Deep Sequencing Facility of TU Dresden. The raw FASTQ files were analyzed with MAGeCK-VISPR .
Genomic DNA isolation and PCR amplification
Genomic DNA was extracted with NucleoSpin® Blood XL (Machery Nagel # 740,950.50) according to the manufacturer’s protocol. The first round PCR of Next Generation Sequence (NGS) was performed with 26 separate 100 μL redundant reactions, each containing 5 μg of DNA, 50 μL Q5® Hot Start High-Fidelity 2X Master Mix (NEB # M0494L), and 3 μL of a 10 μM solution of each primer (P5 and P7). The PCR was performed as described in .
Transmission electron microscopy (TEM)
A549 cells were treated either with vehicle (DMSO) or KHS101 10 μM for 48 and 96 h. After treatment the cells were fixed with 2% (v/v) glutaraldehyde in 0.04 M phosphate buffer pH 7.4 for 60 min, scraped off, collected in a tube and post-fixed with the same fixative for 24 h. The cells were then washed in PBS and resuspended in 15% BSA for 10 min. The cells were centrifuged, and the cell pellet fixed by slowly adding the fixative for 48 h at 4 °C. The cell pellets were cut in smaller pieces and prepared for EM like following. Pellet fragments were washed in phosphate buffer, stained with 1% osmium tetraoxide in phosphate buffer for 90 min, dehydrated in ethanol series and acetone and embedded in TAAB 812 Embedding Resin (T030, TAAB). Semithin (1 µm) sections were cut on Leica Ultracut UCT microtome and stained with toluidine blue. Ultrathin (70 nm) sections were collected on uncoated nickel grids (M200-NI, Electron Microscope Sciences). The grids were stained with 3% uranyl acetate for 15 min at 60 °C and 3% lead citrate (Leica Ultrostain 2) for 6 min at room temperature. The sections were analyzed with JEM-1400 Plus electron microscope, equipped with Quemsa TEM CCD camera and images obtained using Radius software (TEM Imaging Platform software).
Normalized gene expression values for HSPD1 were obtained for lung cancer patient samples from GEO (GSE30219) and mRNA z-score values for the TCGA profile (TCGA LUAD, Cell 2018) from the cbioportal platform . Survival curves were generated using Kaplan–Meier estimate with the samples categorized into ‘low’ and ‘high’ HSPD1 expression groups based on the median of the HSPD1 mRNA expression value. Log-rank test was conducted to obtain the significance between the two groups in R software.
Differential gene expression analysis and correlation
Differential gene expression (DGE) analysis was performed between the KHS101-sensitive and KHS101-resistant cell lines identified from the drug sensitivity analysis. Briefly, the top four KHS101-sensitive cell lines (NCI-H460, NCI-H1581, LOU-NH91 and A549) and top four KHS101-resistant cell line (SK-LU-1, NCI-H1563, BEN and NCI-H838) were chosen for the analysis. Differential expression between the sensitive and resistant groups was performed using cancer cell line encyclopedia (CCLE) gene expression profile (GSE36133) with RMA log2 signal intensity using GEO2R package in GEO with a moderated t-test parameter . 33 up-regulated and 26 down-regulated genes were filtered with a p-value of < 0.05 and log2 fold change cutoff of 2. In parallel, Pearson’s correlation analysis was performed between IC50 values of KHS101 from 21 cell lines and 18,569 genes expression from CCLE. Based on a correlation value of 0.5 and p-value < 0.05, 298 positively and 97 negatively correlated genes were selected. An overlap analysis was performed with the DGEs list and correlation list to identify a stringent gene list of 12 genes. Correlation analysis between HSPD1 gene expression value obtained from GEO (GSE47206) and KHS101 dose response IC50 values was determined using GraphPad Prism 8.
Differential metabolite analysis
Similar to the differential gene expression analysis, a differential metabolite analysis was performed for the top four KHS101-sensitive and KHS101-resistant group cell lines with the information of 225 log10 normalized quantitative metabolites levels across 177 CCLE cell lines . 7 differential metabolites were identified with a log10 fold change of 2 and p-value < 0.05.
Gene ontology analysis
Gene ontology analysis was performed for the top 20 gene candidates identified from the CRISPR knockout screen of A549 cells treated with KHS101. Genes were assessed for the enrichment using GO cellular component 2017b in Enrichr .
Fitness gene and cancer gene dependency
HSPD1 fitness was analyzed as described in . HSPD1 gene’s cancer dependency was analyzed in a panel of NSCLC cell lines (n = 52) from PROJECTDRIVE RNAi screen . Using redundant siRNA activity (RSA) gene level metric of threshold -3, the essentiality of the gene in cell lines was determined.
In vivo experiments
NSG strain (JAX) was used as an experimental model to study the tumor growth of HSPD1 knockdown cells (A549 and H1299). For subcutaneous injections, 0.5 × 106 cells re-suspended in 50 μL 0.9% NaCl were mixed with Matrigel (Corning, 734–0268) in a ratio 1:1 (v:v). Cells were injected in flanks of 9–15-weeks-old NSG with 8 mice per group. Caliper measurements were taken twice a week and tumor volume was calculated using the formula (Length × Width2 × π)/6. After 4/5 weeks, mice were euthanized by cervical dislocation and tumors were isolated and weighted.
C57BL/6 strain was used as experimental model to evaluate KHS101 effect on lung metastasis tumor formation. For the tail vein metastasis assay, 0.5 × 106 LL2 cells were re-suspended in 100 μl PBS and injected in the tail vein of female C57BL/6, with 10 mice per group. After 2 weeks either vehicle (5% (v/v) Ethanol-15% (w/v) (2-Hydroxypropyl)-β-cyclo-dextrin (Captisol Technology)) or 6 mg/kg KHS101 (Cellagen Technology) were injected twice a day subcutaneously for 10 days. Lung metastases were monitored by bioluminescence imaging (BLI). Anesthetized mice were intraperitoneally injected with 50 mg/mL D-luciferin (Kayman Chemicals). Bioluminescence images were acquired with Lumina III in vivo Imaging System (IVIS, Perkin Elmer). Time of surrogate survival with disease was calculated as difference between the experimental endpoint and the day of appearance of the bioluminescence signal.
Statistical tests were performed with the GraphPad software (Prism) v.8 comparing groups of different conditions. In all tests at least three biological replicates were used and the statistical significance was considered at p ≤ 0.05.
HSPD1 is a fitness gene and a potential target for NSCLC
To understand whether HSPD1 could be an attractive target for NSCLC, we performed a fitness analysis from a previously published dataset of a pan-cancer and genome-wide CRISPR/Cas9 screening of 18,009 genes, which included n = 21 adenocarcinoma (ADC) and n = 11 squamous cell carcinoma (SqCC) cell lines . As a result, the HSPD1 gene scored high, ranking 25th in ADC and 71st in SqCC cell lines, out of 805 and 559 total ‘fitness’ genes, respectively (Fig. 1A). This initial observation was supported by the results from an independent dataset based on a large-scale RNAi screen of 7,837 genes performed in 52 NSCLC cell lines . HSPD1 was a significant cancer-dependent gene in over 70% of the cells (Fig. 1B, Additional file 1: Table S2), confirming the importance of this gene for NSCLC survival. In addition, survival analysis on HSPD1 mRNA expression showed that the patients with ‘high’ levels of HSPD1 carry a worse prognosis compared to those with ‘low’ HSPD1, in terms of overall (Fig. 1C), relapse-free (Fig. 1D) and disease-specific survival (Fig. 1E). We then investigated the prevalence of HSPD1 expression in tissue samples from a small cohort of NSCLC patients (n = 30) by immunohistochemistry (IHC), and observed positive HSPD1 protein expression in 100% of the cases (Fig. 1F-I). The expression pattern was cytoplasmic, granular or diffuse, with some degree of heterogeneity (semi quantitatively scored in a three-level scale), with very high protein levels identified in 64% of specimens (Fig. 1J). No difference was observed between ADC and SqCC samples in terms of HSPD1 expression levels. To confirm HSPD1 expression in NSCLC cell lines grown in vitro, a panel of 10 different cell lines, including three mouse-derived ones (Ladi3.1, LL2 and Ladi2.1) , was subjected to western blot quantification. The results indicated detectable HSPD1 protein in all samples (Fig. 1K). Separated protein isolation from mitochondria and cytosolic fractions confirmed that HSPD1 is predominantly localized in the mitochondria in NSCLC cells (Fig. 1L). Therefore, HSPD1 is a mitochondrial protein essential for NSCLC with a strong prognostic value and high prevalence of expression in tissue samples, making it a very valuable candidate for targeted inhibition.
HSPD1 promotes NSCLC growth in vitro and in vivo
To study the function of HSPD1 in NSCLC cells in vitro, we analyzed the effects of its knockdown by transducing 5 different NSCLC cell lines with up to 5 independent shRNA sequences (called #44, #45, #46, #47 and #48) (Fig. 2A, Additional file 1: Fig. S1A). Proliferation was evaluated through a real-time cell confluency assay (using the Incucyte system). Shortly after shRNA-based inhibition of HSPD1 expression, the cells showed a decrease in proliferation compared to the controls transduced with a scrambled/non-targeting shRNA, referred to as pLKO (Fig. 2B, Additional file 1: Fig. S1B). In particular, the phenotype was stronger for shRNA sequences providing higher levels of HSPD1 suppression (i.e., #48). HSPD1 knockdown cells were not able to form colonies compared to the control cells when plated at lower density (1,500 cells/well) as observed over two weeks (Fig. 2C-D, Additional file 1: Fig. S1C-D). In line with a strong and durable effect on cell proliferation, a cell cycle analysis (by FACS) showed a significant increase in the percentage of cells in sub-G1 and G1 phases in HSPD1 knockdown cells, in contrast to a reduction in cells residing in S or G2-M phases (Fig. 2E-H, Additional file 1: Fig. S1E-F). To evaluate whether HSPD1 loss-of-function could also affect the growth of normal lung cells, we induced the knockdown in an immortalized non-tumorigenic bronchial epithelial cell line (BEAS-2B). As a result, BEAS-2B cells showed a dramatically reduced growth in vitro (Additional file 1: Fig. S1G-H); however, they were able to completely recover their proliferative capacity as early as one week post infection in contrast to the cancer cells showing a sustained reduction in proliferation capacity (Additional file 1: Fig. S1I).
HSPD1 has a critical function in mitochondrial metabolism . Therefore, to investigate whether the growth suppression observed with HSPD1 knockdown was linked to metabolic alterations, an extracellular metabolic flux analysis was performed using two different cell lines. Throughout the experiments, knockdown cells showed a significant reduction in their basal respiration and in their ATP-linked respiration, suggesting that the alterations in the mitochondrial metabolism may be responsible for the reduction in cell proliferation (Fig. 3A-B). To investigate whether this reduction could be due to a decrease in mitochondrial number, a specific-mitochondrial staining was performed in HSPD1 knockdown cells, which, did not show any difference in the fluorescent intensity of the mitochondrial dye (Additional file 1: Fig. S2A). Finally, we tested whether HSPD1 knockdown inhibited cell growth in vivo. To this end, we injected shHSPD1-expressing A549 and H1299 cells subcutaneously in NSG mice to evaluate the tumor growth rate over time. The results markedly indicated that HSPD1 knockdown cells formed very small and slower-growing tumors compared to control pLKO cells, in terms of both tumor volume and tumor weight (Fig. 3C-H), an effect which was maintained for the entire duration of the experiment, hence confirming that HSPD1 is required for NSCLC growth.
The HSPD1- targeting small molecule KHS101 arrests NSCLC growth
KHS101 is a synthetic small molecule reported to impair the growth of glioblastoma cells in vitro and in vivo via inhibition of HSPD1 chaperone activity . As KHS101 leads to metabolic exhaustion and apoptosis in glioblastoma cells, we first tested as to whether the compound elicits cytotoxicity in NSCLC cell lines. After treatment with KHS101, 6 human NSCLC cell lines showed a significant reduction in cell proliferation that was dependent on the compound concentration (Fig. 4A, Additional file 1: Fig. S2B-C). To verify the specificity of HSPD1 targeting, we exposed the cells to HB072, an HSPD1 non-targeting chemical analogue of KHS101  (Fig. 4B), which did not significantly affect cell growth (Fig. 4C, Additional file 1: Fig. S3A). After 5 days of continuous drug treatment, NSCLC cells lost their clonogenic growth ability, as shown in a colony formation assay (Fig. 4D-E, Additional file 1: Fig. S3B). In line with HSPD1 knockdown experiments, cell cycle analysis showed a significant reduction in the percentage of cells in S and G2-M phases concomitant with an increase of cells in G1 phase (Fig. 4F-I), as assessed 24 h after KHS101 treatment. Finally, since many conventional anti-cancer therapeutic strategies fail to target non-proliferating cancer cells, we sought to test whether the KHS101-induced growth inhibitory effect was maintained altering the proliferative state of the cells in a controlled fashion. Therefore, after being synchronized overnight using a low FBS concentration (0.5%), the cells were treated with KHS101 in presence of media with increasing FBS, ranging from 0 to 5%. As a result, no differences were measured in KHS101 efficacy in the different conditions by a real-time proliferation assay (Additional file 1: Fig. S3C), indicating that KHS101 is effective independently of the cells’ proliferative status. In conclusion, these data indicate that chemical targeting of HSPD1 markedly reduces the growth of NSCLC cell lines.
KHS101 induces metabolic disruption in NSCLC cells leading to cell death
To test whether KHS101 affects energy metabolism in NSCLC cell lines, we performed an extracellular flux analysis 24 h after KHS101 treatment. The results indicated a strong reduction in cellular basal respiration and in ATP-coupled respiration capacity. Concomitantly, a significant increase was observed for basal glycolysis (Fig. 5A-B). We next asked as to how fast KHS101 would induce such alterations and treated the cells for shorter periods (ranging from 0.5 h to 6 h), observing that glycolytic activity increases simultaneously with oxidative phosphorylation (OXPHOS) reduction from less than one hour after treatment (Additional file 1: Fig. S4A-B), and suggesting that the metabolic rewiring induced by KHS101 is an early event, as previously described . We further investigated metabolic alterations evaluating the expression of some key metabolic enzymes in A549 HSPD1 knockdown cells and A549 cells treated with KHS101 for 72 h. We observed a significant increase in hexokinase 2 (HK2) expression (Additional file 1: Fig. S4C), which was consistent with the detected increase in glycolysis. In addition to metabolic flux phenograms, cell morphology was found altered in cells exposed to 10 µM KHS101. Treated cells presented an altered ultrastructure and, after 96 h, severe swelling of mitochondria and sparse presence of cristae was observed. Moreover, an increased number of secondary lysosomes was found indicating a marked cellular degeneration associated with a loss of microvilli (Fig. 5C). These alterations were already detectable as early as 48 h after treatment (Additional file 1: Fig. S4D). As a result of these alterations, and differently to HPSD1 knockdown, KHS101 treatment induced cell death, as indicated by the incorporation of fluorescent Cytotox green dye (Fig. 5D, Additional file 1: Fig. S4E). To test whether the affected cells were undergoing apoptosis, we performed a caspase 3/7 activity assay, which showed specific activation of apoptosis 2 days after the treatment (Fig. 5E, Additional file 1: Fig. S4F-G). Co-administration of the pan-caspase inhibitor Z-VAD-FMK could efficiently suppress this caspases activation (Additional file 1: Fig. S4I) but did not rescue the cells from KHS101 cytotoxicity (Fig. 5F, Additional file 1: Fig. S4H), thus indicating that cell death is not occurring solely as a result of apoptotic induction. Interestingly, no protection was observed also in presence of necroptosis or ferroptosis inhibitors (Additional file 1: Fig. S4J-K), suggesting that the metabolic alterations caused by KHS101 could not be rescued by blocking cell death pathways. Furthermore, to investigate the impact of KHS101 in vivo, we tested the sensitivity of the tumor-forming mouse-derived NSCLC cell line LL2 (Additional file 1: Fig. S5A-B), and subsequently injected them in the tail vein of BL6 mice to evaluate the growth of lung metastasis upon KHS101 treatment. Lung tumors-bearing mice that were treated subcutaneously for 2 weeks with 6 mg/kg KHS101 (twice a day) did not show a significant difference in tumor area, as evaluated by IVIS measurement. However, KHS101-treated mice showed a prolonged disease-specific survival time (Additional file 1: Fig. S5C), in line with a non-significant trend for a reduced number of lung lesions, as macroscopically detected (Additional file 1: Fig. S5D).
Finally, we investigated the in vitro effects of KHS101 in combination with cisplatin, one of the most commonly used chemotherapeutic drug used in relapsing NSCLC . Strikingly, a low cisplatin dose (1 or 3 µM) that alone was not able to strongly impair cell growth or survival, showed to induce a marked growth suppression (Fig. 5G, Additional file 1: Fig. S6A) and increase in cell death (Fig. 5H, Additional file 1: Fig. S6B) when combined with KHS101 (also at sub-lethal doses), resulting in a 2 to 4-fold reduced relative growth compared to cisplatin alone in three different cell lines (Additional file 1: Fig. S6C). Cisplatin exerts its activity via upregulation of reactive oxygen species (ROS)  and by fluorescent ROS staining of treated cells we found that the combination of cisplatin and KHS101 synergistically increased ROS activation (Fig. 5I-J). These data suggest that HSPD1 targeting can also induce cell death and increase the efficacy of chemotherapeutic treatments.
Sensitivity to KHS101 is related to the metabolic state of the cells
We further profiled in vitro sensitivity to KHS101 by performing a screening on 26 different lung cancer cell lines, measuring cell viability 72 h after treatment. Interestingly, none of the cell lines in the panel appeared to be completely resistant (Fig. 6A, Additional file 1: Fig. S7A). KHS101-induced cytotoxicity was independent of the lung cancer type (NSCLC or SCLC) or NSCLC histological subtype (Fig. 6B). However, cells showed different grades of sensitivity towards the compound with IC50 values ranging from 2 to 14 μM (Fig. 6C, Additional file 1: Fig. S7A). These results were independently validated comparing KHS101 sensitivity in the most sensitive and resistant cell lines based on percentage of inhibition of cell proliferation (Additional file 1: Fig. S8A). To find possible determinants of drug resistance, cell lines with different genetic alterations were compared, but no particular genetic background was found as significantly protective (Additional file 1: Table S3). To identify possible molecular determinants of drug sensitivity, we performed a transcriptomic analysis on cell lines based on their IC50 values. HSPD1 mRNA expression correlated with IC50 values obtained from the screening (Fig. 6D). In order to have a stringent and genome-wide comparison, we used a 2-sided statistical approach. First, the 4 most resistant (NCI-H838, BEN, NCI-H1563 and SK-LU-1) and the 4 most sensitive NSCLC cell lines (H460, NCI-H1581, LOU-NH91 and A549) were grouped and the differentially expressed genes between the two groups were identified and ranked according to their fold differences (Additional file 1: Table S4-5). Moreover, using available data from all the cells in the screening, we correlated the expression of each gene with the KHS101 IC50 values and ranked them based on Pearson metrics (Additional file 1: Table S6-7). The results from these two analyses were overlapped (Additional file 1: Fig. S8B), obtaining a list of potential up- and down- regulated genes in cells resistant to KHS101 (Fig. 6E). The top up-regulated resistance-conferring gene was the transcription factor NLRC5, which acts as a transactivator for the MHC class I genes  (some of which were also up-regulated in the analysis, like HLA-J or HLA-F). To validate this finding, we lentivirally overexpressed NLRC5 in sensitive NSCLC cells (Fig. 6F), and detected a significantly decreased cytotoxic effect of KHS101, indicating cellular resistance (Fig. 6G, Additional file 1: Fig. S8C). Among the genes down-regulated in the resistant cell lines, the top hit was SLC6A8 (Fig. 6E), a sodium chloride-dependent creatine transporter 1 . This evidence was validated by introducing a specific CRISPR/Cas9-based knock-out of SLC6A8, which increased cells viability of sensitive cell lines treated with KHS101 (Fig. 6H). Since SLC6A8 is a creatine transporter we speculated that it could alter KHS101 sensitivity by changing the metabolic features of the cells. Analysis of a previously published metabolome profiling of NSCLC cell lines  revealed that cells that we determined to be the most KHS101 sensitive had a higher creatine level (Fig. 6I), which was consistent with the results from the transcriptomic analysis (i.e. higher expression of the importer) suggesting a possible effect of creatine metabolism in the response to KHS101. Creatine can be phosphorylated by the cytosolic creatine kinase (C-CK) or by the mitochondrial isoenzymes (MtCK) to provide a rapid source of ATP  thereby fueling ATP-dependent cellular processes. To test whether differences in cellular creatine levels could be correlated with the cells’ metabolic state, we performed extracellular flux analysis of cells with higher (A549 and H460) and lower creatine (H838 and BEN) baseline levels  and found that cells with higher creatine had higher basal respiration (Fig. 6J). To prove this correlation further functionally, we inhibited creatine kinase activity with the specific inhibitor DNFB which decreased the basal respiration level of A549 cells (Fig. 6K), concomitantly increasing the cells’ resistance to KHS101 as assessed by cell viability assays (Fig. 6L). Taken together these findings suggest that a higher dependency on OXPHOS due to the alteration of creatine metabolism activity could determine the sensitivity to KHS101 of NSCLC cells.
KHS101 sensitivity depends on OXPHOS activity
To further investigate KHS101 resistance at a broader level, we performed a genome-wide CRISPR/Cas9-based dropout screening. A genetically barcoded whole-genome library (GeCKO v2) was transduced in A549 cells overexpressing the Cas9 nuclease. These cells were treated with two different KHS101 concentrations (10 and 15 μM) or DMSO as control (Fig. 7A). None of the cells were able to survive the treatment with the 15 μM KHS101 dose, suggesting that is not possible to develop a complete resistance to this compound by the elimination of a single gene. By contrast, we could observe non-growing and morphologically altered cells surviving the 10 μM dose. After 2 weeks of survival (after all control pLKO-infected Cas9-A549 cells treated with the same KHS101 dose underwent complete cell death), the drug treatment was stopped to allow the selected cells to grow and be processed for the DNA preparation and the next-generation sequencing steps. The sequencing results on untreated cells indicated that it was possible to detect a total of 116,700 sgRNAs (all gRNAs excluding those with no counts) present in the library, with a 97% representation, a level similar to previous screenings , while in the treated group this percentage decreased to 56%, as a result of the drug selection. ‘Hits’ were ranked and gene ontology analysis on the 20 highest ranking ‘hits’ (Additional file 1: Table S8) indicated that the genes whose knockout increased cellular resistance to KHS101 were mainly involved in regulating mitochondrial structure and activity (Fig. 7B). The COX5B gene, a subunit of the cytochrome c oxidase, was the ‘top hit’. In order to independently validate this finding, we knocked out COX5B by two independent gRNAs in A549 cells, and in both cases, we found the emergence of a more KHS101-resistant phenotype (Fig. 7C). COX5B knock-out cells presented a reduction in their basal respiration and ATP-linked respiration (Fig. 7D), strengthening the notion that KHS101 sensitivity is linked with cellular oxidative activity, in line with the results obtained from creatine metabolism. Conversely, blocking glycolysis by 2-deoxy-D-glucose (2-DG) treatment significantly enhanced KHS101 sensitivity in different cell lines (Fig. 7E), further confirming that cancer cells respond to KHS101 proportionally to their OXPHOS-dependence. Therefore, HSPD1 targeting induces a mitochondrial metabolic breakdown in NSCLC, which results in loss of growth potential and cell death, and the small molecule KHS101 is particularly active in OXPHOS-dependent cells.
The identification of novel targets and the development of more effective drugs is of paramount importance for improving the clinical management of NSCLC. Metabolism has recently gained much attention in the oncological field and anti-metabolism strategies have been proposed to reduce the lethal effects of NSCLC [41,42,43]. We report here heat shock protein HSPD1 as a theranostic metabolic marker for NSCLC, whose elimination can have profound effect on outcome of the patients.
Heat shock proteins (HSPs) are molecular chaperones categorized according to their molecular sizes into small HSPs and large HSPs . HSPs can be highly expressed in cancer cells, where they promote cell proliferation, metastasis and drug resistance. HSPD1 is an ATP-dependent HSP mainly localized in the mitochondria, playing an essential role in guaranteeing the correct folding of the mitochondrial-imported proteins  and some data are available on its role as a stemness/metastasis regulator in other cancers [22, 45, 46]. In lung cancer, HSPD1 has been previously found highly expressed in NSCLC tissues  and identified as a predictive marker for survival in both smokers and non-smokers patients . Our findings confirm that HSPD1 is ubiquitously expressed at a high level (as detectable by IHC in all samples) and is associated with adverse prognosis, adding the notion that its expression is essential for NSCLC survival (fitness gene), making it a very attractive target for future therapy. In order to more functionally explore this possibility, we evaluated the effect of HSPD1 targeting achieved by either shRNA-mediated knockdown or by a specific small molecule. In shRNA experiments, the level of growth suppression measured for each shRNA sequence was generally proportional to the degree of HSPD1 knockdown. Importantly, even if a complete protein elimination was not achieved, the cell growth blockage obtained by the best performing shRNA (#48) was durable and very effective both in vitro and in vivo in all the cell lines investigated, suggesting that this behavior is not correlated to NSCLC histotypes or HSPD1 expression level. When cellular metabolic changes were measured, cells with HSPD1 knockdown presented a strong reduction in basal respiration and impaired capacity in synthesizing ATP though OXPHOS, in line with the role of HSPD1 in the (re)folding of mitochondrial-imported oxidative respiration proteins . Noteworthy, while this investigation mainly focused on the effects on metabolism, other functions of HSPD1 might have contributed to the loss of viability observed upon its knockdown and could be further investigated. For example, HSPD1 can act as immunomodulant , can stimulate phagocytosis of senescent cells , and regulate apoptosis and intracellular trafficking [52, 53]. In any case, our data clearly indicate the activity of HSPD1 as essential for maintaining NSCLC metabolic fitness and its loss causes a profound energetic breakdown affecting the ability of the cancer cells to divide and expand, making it an attractive therapeutic target.
An emerging way of chemically disrupting the functionality of HSPD1 is the treatment with the synthetic small compound KHS101. This compound was originally identified through a phenotypical screen as an inducer of neuronal differentiation both in vitro and in vivo . KHS101 promotes HSPD1 aggregation thereby altering the metabolic activity of glioblastoma cells . In NSCLC, we observed that KHS101-treated cells were very severely growth arrested, in analogy to what was observed with the HSPD1 knockdown approach. OXPHOS activity was also efficiently repressed, even with low micromolar KHS101 concentrations. In contrast to HSPD1 knockdown, KHS101-treated cells increased glycolysis, suggesting that they might try to compensate for the loss of ATP production. This increase is, however, not sufficient for the cells to gain resistance and survive the treatment. Indeed, KHS101-treated NSCLC cells underwent significant cell death, probably due to the formation of irreversible HSPD1 multi-aggregates coagulated by KHS101, which can disrupt the energetic capacity of the cancer cells, as previously shown . Notably, co-administration of specific inhibitors of death pathways did not rescue KHS101 treated cells, indicating that the occurring death is strong, multifactorial, and irreversible, all features of crucial importance for a novel anti-cancer therapeutic strategy.
In order to better understand the biology of HSPD1 targeting, two independent unbiased screenings were carried out: a differential gene expression analysis on KHS101 sensitive and resistant cells and a genome-wide CRISPR/Cas9-based dropout screening on drug treated cells. The first approach identified a biomarker of increased sensitivity, the creatine importer SLC6A8, which led us to discover OXPHOS dependency as condition for better KHS101 efficacy. Moreover, a potential role in determining KHS101 resistance of the MHC class I genes transactivator NLRC5 has been here identified. Based on its importance on inflammatory processes in tumors , the molecular details of this link should be in the future determined. The second approach came independently to a similar and converging conclusion, by identifying mitochondrial metabolism genes, and in particular the subunit of the cytochrome c oxidase complex (COX5B), as enriched in cells surviving KHS101. This evidence highlights the targeting of mitochondrial metabolism as an effective therapeutic strategy to achieve a complete, irreversible elimination of NSCLC.
The occurrence of drug-overcoming mutations is a typical feature of NSCLC, which has unfortunately so far limited the success of targeted therapies . It is noteworthy that not a single cell from the whole-genome CRISPR/Cas9 library was able to survive the KHS101 selection when administered at the highest dose (15 μM). At the lower dose (10 μM), a small population of cells was able to persist, but could not grow (and be processed for sequencing) until the drug was entirely removed from the experiment. This is a very important observation in line with the observed inability of death pathways inhibitors to protect in vitro from KHS101, and in line with the data from the two independent fitness analyses on HSPD1 gene and with the drug screening from large collections of NSCLC cell lines of various genetic backgrounds, which confirmed that no single genetic alteration can guarantee a complete protection from the treatment with KHS101. The sensitivity towards KHS101 is also independent from the lung cancer type or NSCLC histotype, suggesting that a HSPD1 targeting approach might be valuable also for other cancer types. Another aspect of fundamental importance to consider for NSCLC drug development is the potential successful combination with chemotherapy, as chemotherapy still represents the backbone of therapeutic regimens, especially in late-stage lung cancer patients [56,57,58]. Even if the data here presented will require future dedicated experimentation, the results are encouraging, as KHS101 acted synergistically with cisplatin in vitro, enhancing the production of reactive oxygen species and reducing the dose of cisplatin needed to slow down or kill the cancer cells. These observations reinforce the potential relevance of designing and developing HSPD1-targeting drugs for NSCLC therapy. This goal, however, is not likely to be achieved by KHS101 itself, which this study found active at relatively high concentrations. KHS101 therefore stands out as a rather excellent tool compound for accelerating discoveries for HSPD1 targeting, but novel drugs (or modifications of the same compound) with improved affinity will need to be designed and tested, for instance by taking advantage of emerging targeting techniques based on induction of selective intracellular proteolysis [59, 60]. The activity of these drugs could be in the future tested beyond NSCLC, as HSPD1 appears to be a fundamental oncoprotein also in other cellular contexts [45, 46].
In conclusion, we showed that HSPD1 elimination interferes with NSCLC metabolic activity causing a strong OXPHOS-dependent energetic impairment, which the cancer cells fail to overcome, highlighting HSPD1 as a powerful theranostic marker for improving lung cancer therapy.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Non-small cell lung cancer
Heat shock proteins
Squamous cell lung cancer
Oxygen consumption rate
Extracellular acidification rate
Differential gene expression
Mitochondrial heat shock protein family D member 1
NOD-Like Receptor C5
Solute Carrier Family 6 Member 8
Cytochrome C Oxidase Subunit 5B
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature. 2018;553(7689):446–54.
Kleczko EK, Kwak JW, Schenk EL, Nemenoff RA. Targeting the Complement Pathway as a Therapeutic Strategy in Lung Cancer. Front Immunol. 2019;10:954.
Osmani L, Askin F, Gabrielson E, Li QK. Current WHO guidelines and the critical role of immunohistochemical markers in the subclassification of non-small cell lung carcinoma (NSCLC): Moving from targeted therapy to immunotherapy. Semin Cancer Biol. 2018;52(Pt 1):103–9.
Duma N, Santana-Davila R, Molina JR. Non-Small Cell Lung Cancer: Epidemiology, Screening, Diagnosis, and Treatment. Mayo Clin Proc. 2019;94(8):1623–40.
Luengo A, Abbott KL, Davidson SM, Hosios AM, Faubert B, Chan SH, et al. Reactive metabolite production is a targetable liability of glycolytic metabolism in lung cancer. Nat Commun. 2019;10(1):5604.
Vander Heiden MG, DeBerardinis RJ. Understanding the Intersections between Metabolism and Cancer Biology. Cell. 2017;168(4):657–69.
Yu L, Lu M, Jia D, Ma J, Ben-Jacob E, Levine H, et al. Modeling the Genetic Regulation of Cancer Metabolism: Interplay between Glycolysis and Oxidative Phosphorylation. Cancer Res. 2017;77(7):1564–74.
Faubert B, Solmonson A, DeBerardinis RJ. Metabolic reprogramming and cancer progression. Science. 2020;368(6487):eaaw5473.
Kim J, DeBerardinis RJ. Mechanisms and Implications of Metabolic Heterogeneity in Cancer. Cell Metab. 2019;30(3):434–46.
Jin N, Bi A, Lan X, Xu J, Wang X, Liu Y, et al. Identification of metabolic vulnerabilities of receptor tyrosine kinases-driven cancer. Nat Commun. 2019;10(1):2701.
Fenouille N, Bassil CF, Ben-Sahra I, Benajiba L, Alexe G, Ramos A, et al. The creatine kinase pathway is a metabolic vulnerability in EVI1-positive acute myeloid leukemia. Nat Med. 2017;23(3):301–13.
Bader DA, Hartig SM, Putluri V, Foley C, Hamilton MP, Smith EA, et al. Mitochondrial pyruvate import is a metabolic vulnerability in androgen receptor-driven prostate cancer. Nat Metab. 2019;1(1):70–85.
Jacob P, Hirt H, Bendahmane A. The heat-shock protein/chaperone network and multiple stress resistance. Plant Biotechnol J. 2017;15(4):405–14.
Dubrez L, Causse S, Borges Bonan N, Dumetier B, Garrido C. Heat-shock proteins: chaperoning DNA repair. Oncogene. 2020;39(3):516–29.
Chatterjee S, Burns TF. Targeting Heat Shock Proteins in Cancer: A Promising Therapeutic Approach. Int J Mol Sci. 2017;18(9):1978.
Wu J, Liu T, Rios Z, Mei Q, Lin X, Cao S. Heat Shock Proteins and Cancer. Trends Pharmacol Sci. 2017;38(3):226–56.
Gomez-Llorente Y, Jebara F, Patra M, Malik R, Nisemblat S, Chomsky-Hecht O, et al. Structural basis for active single and double ring complexes in human mitochondrial Hsp60-Hsp10 chaperonin. Nat Commun. 2020;11(1):1916.
Teng R, Liu Z, Tang H, Zhang W, Chen Y, Xu R, et al. HSP60 silencing promotes Warburg-like phenotypes and switches the mitochondrial function from ATP production to biosynthesis in ccRCC cells. Redox Biol. 2019;24:101218.
Yun CW, Kim HJ, Lim JH, Lee SH. Heat Shock Proteins: Agents of Cancer Development and Therapeutic Targets in Anti-Cancer Therapy. Cells. 2019;9(1):60.
Xu X, Wang W, Shao W, Yin W, Chen H, Qiu Y, et al. Heat shock protein-60 expression was significantly correlated with the prognosis of lung adenocarcinoma. J Surg Oncol. 2011;104(6):598–603.
Fucarino A, Pitruzzella A. Role of HSP60/HSP10 in Lung Cancer: Simple Biomarkers or Leading Actors? J Oncol. 2020;2020:4701868.
Polson ES, Kuchler VB, Abbosh C, Ross EM, Mathew RK, Beard HA, et al. KHS101 disrupts energy metabolism in human glioblastoma cells and reduces tumor growth in mice. Sci Transl Med. 2018;10(454):eaar2718.
Siddiqui A, Gollavilli PN, Schwab A, Vazakidou ME, Ersan PG, Ramakrishnan M, et al. Thymidylate synthase maintains the de-differentiated state of triple negative breast cancers. Cell Death Differ. 2019;26(11):2223–36.
Gollavilli PN, Parma B, Siddiqui A, Yang H, Ramesh V, Napoli F, et al. The role of miR-200b/c in balancing EMT and proliferation revealed by an activity reporter. Oncogene. 2021;40(12):2309–22.
Siddiqui A, Vazakidou ME, Schwab A, Napoli F, Fernandez-Molina C, Rapa I, et al. Thymidylate synthase is functionally associated with ZEB1 and contributes to the epithelial-to-mesenchymal transition of cancer cells. J Pathol. 2017;242(2):221–33.
Sanjana NE, Shalem O, Zhang F. Improved vectors and genome-wide libraries for CRISPR screening. Nat Methods. 2014;11(8):783–4.
Li W, Xu H, Xiao T, Cong L, Love MI, Zhang F, et al. MAGeCK enables robust identification of essential genes from genome-scale CRISPR/Cas9 knockout screens. Genome Biol. 2014;15(12):554.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.
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.
Li H, Ning S, Ghandi M, Kryukov GV, Gopal S, Deik A, et al. The landscape of cancer cell line metabolism. Nat Med. 2019;25(5):850–60.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.
Siddiqui MA, Gollavilli PN, Ramesh V, Parma B, Schwab A, Vazakidou ME, et al. Thymidylate synthase drives the phenotypes of epithelial-to-mesenchymal transition in non-small cell lung cancer. Br J Cancer. 2021;124(1):281–9.
McDonald ER, 3rd, de Weck A, Schlabach MR, Billy E, Mavrakis KJ, Hoffman GR, et al. Project DRIVE: A Compendium of Cancer Dependencies and Synthetic Lethal Relationships Uncovered by Large-Scale, Deep RNAi Screening. Cell. 2017;170(3):577–92.
Behan FM, Iorio F, Picco G, Goncalves E, Beaver CM, Migliardi G, et al. Prioritization of cancer therapeutic targets using CRISPR-Cas9 screens. Nature. 2019;568(7753):511–6.
Senan S, Brade A, Wang LH, Vansteenkiste J, Dakhil S, Biesma B, et al. PROCLAIM: Randomized Phase III Trial of Pemetrexed-Cisplatin or Etoposide-Cisplatin Plus Thoracic Radiation Therapy Followed by Consolidation Chemotherapy in Locally Advanced Nonsquamous Non-Small-Cell Lung Cancer. J Clin Oncol. 2016;34(9):953–62.
Kim SJ, Kim HS, Seo YR. Understanding of ROS-Inducing Strategy in Anticancer Therapy. Oxid Med Cell Longev. 2019;2019:5381692.
Benko S, Kovacs EG, Hezel F, Kufer TA. NLRC5 Functions beyond MHC I Regulation-What Do We Know So Far? Front Immunol. 2017;8:150.
Ji L, Zhao X, Zhang B, Kang L, Song W, Zhao B, et al. Slc6a8-Mediated Creatine Uptake and Accumulation Reprogram Macrophage Polarization via Regulating Cytokine Responses. Immunity. 2019;51(2):272–84.
Snow RJ, Murphy RM. Creatine and the creatine transporter: a review. Mol Cell Biochem. 2001;224(1–2):169–81.
Ramesh V, Brabletz T, Ceppi P. Targeting EMT in Cancer with Repurposed Metabolic Inhibitors. Trends in cancer. 2020;6(11):942–50.
Kawada K, Toda K, Sakai Y. Targeting metabolic reprogramming in KRAS-driven cancers. Int J Clin Oncol. 2017;22(4):651–9.
Obrist F, Michels J, Durand S, Chery A, Pol J, Levesque S, et al. Metabolic vulnerability of cisplatin-resistant cancers. EMBO J. 2018;37(14):e98597.
Marino Gammazza A, Macaluso F, Di Felice V, Cappello F, Barone R. Hsp60 in Skeletal Muscle Fiber Biogenesis and Homeostasis: From Physical Exercise to Skeletal Muscle Pathology. Cells. 2018;7(12):224.
Tsai YP, Yang MH, Huang CH, Chang SY, Chen PM, Liu CJ, et al. Interaction between HSP60 and beta-catenin promotes metastasis. Carcinogenesis. 2009;30(6):1049–57.
Berger E, Rath E, Yuan D, Waldschmitt N, Khaloian S, Allgäuer M, et al. Mitochondrial function controls intestinal epithelial stemness and proliferation. Nat Commun. 2016;7:13171.
Michils A, Redivo M, Zegers de Beyl V, de Maertelaer V, Jacobovitz D, Rocmans P, et al. Increased expression of high but not low molecular weight heat shock proteins in resectable lung carcinoma. Lung Cancer. 2001;33(1):59–67.
Sotgia F, Lisanti MP. Mitochondrial markers predict survival and progression in non-small cell lung cancer (NSCLC) patients: Use as companion diagnostics. Oncotarget. 2017;8(40):68095–107.
Zhou C, Sun H, Zheng C, Gao J, Fu Q, Hu N, et al. Oncogenic HSP60 regulates mitochondrial oxidative phosphorylation to support Erk1/2 activation during pancreatic cancer cell growth. Cell Death Dis. 2018;9(2):161.
Zininga T, Ramatsui L, Shonhai A. Heat Shock Proteins as Immunomodulants. Molecules (Basel, Switzerland). 2018;23(11):2846.
Goh YC, Yap CT, Huang BH, Cronshaw AD, Leung BP, Lai PB, et al. Heat-shock protein 60 translocates to the surface of apoptotic cells and differentiated megakaryocytes and stimulates phagocytosis. Cell Mole Life Sci. 2011;68(9):1581–92.
Saini J, Sharma PK. Clinical, Prognostic and Therapeutic Significance of Heat Shock Proteins in Cancer. Curr Drug Targets. 2018;19(13):1478–90.
Zhang L, Fok JH, Davies FE. Heat shock proteins in multiple myeloma. Oncotarget. 2014;5(5):1132–48.
Wurdak H, Zhu S, Min KH, Aimone L, Lairson LL, Watson J, et al. A small molecule accelerates neuronal differentiation in the adult rat. Proc Natl Acad Sci U S A. 2010;107(38):16542–7.
Lim ZF, Ma PC. Emerging insights of tumor heterogeneity and drug resistance mechanisms in lung cancer targeted therapy. J Hematol Oncol. 2019;12(1):134.
Pirker R. Chemotherapy remains a cornerstone in the treatment of nonsmall cell lung cancer. Curr Opin Oncol. 2020;32(1):63–7.
Visconti R, Morra F, Guggino G, Celetti A. The between Now and Then of Lung Cancer Chemotherapy and Immunotherapy. Int J Mol Sci. 2017;18(7):1374.
Scagliotti GV, Ceppi P, Capelletto E, Novello S. Updated clinical information on multitargeted antifolates in lung cancer. Clin Lung Cancer. 2009;10(Suppl 1):S35-40.
Khan S, He Y, Zhang X, Yuan Y, Pu S, Kong Q, et al. PROteolysis TArgeting Chimeras (PROTACs) as emerging anticancer therapeutics. Oncogene. 2020;39(26):4909–24.
Sun X, Gao H, Yang Y, He M, Wu Y, Song Y, et al. PROTACs: great opportunities for academia and industry. Signal Transduct Target Ther. 2019;4:64.
This work was supported by the Interdisciplinary Center for Clinical Research of the University of Erlangen-Nuremberg (PC) and the German Research Foundation grant DFG-CE-281/5–1 (PC; BP). A. Siddiqui is supported by the Independent Research Fund Denmark grant number 0134-00012B to PC. CP is supported by the European Research Council (MSCA-ITN-ETN PRECODE 861196). Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
The use of retrospective LC tissues for immunohistochemical study was approved by the Research Ethics Committee of the San Luigi Hospital/University of Turin (approval n.26 dated January 18th 2019). Informed consent was obtained from all patients. Samples have been anonymized by a staff member not involved in the study and all patient information are made unavailable to the investigators. In vivo experiments were performed by skilled experimenters trained according to FELASA guidelines. Animal protocols were approved by the Institutional Animal Care and Use Committee of the Regierung von Unterfranken.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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
Parma, B., Ramesh, V., Gollavilli, P.N. et al. Metabolic impairment of non-small cell lung cancers by mitochondrial HSPD1 targeting. J Exp Clin Cancer Res 40, 248 (2021). https://doi.org/10.1186/s13046-021-02049-8
- Non-small cell lung cancer