- 1Faculty of Health, School of Biomedical Sciences, Institute of Health and Biomedical Innovation, Translational Research Institute, Queensland University of Technology, Woolloongabba, QLD, Australia
- 2Faculty of Science and Engineering, School of Chemistry and Physics, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, QLD, Australia
- 3Cancer Services, Princess Alexandra Hospital, Woolloongabba, QLD, Australia
- 4Department of Clinical Medicine, Trinity Translational Medicine Institute, Trinity College Dublin, Dublin, Ireland
- 5Thoracic Oncology Research Group, Labmed Directorate, St. James’s Hospital, Dublin, Ireland
Platinum-based chemotherapy remains the cornerstone of treatment for most people with non-small cell lung cancer (NSCLC), either as adjuvant therapy in combination with a second cytotoxic agent or in combination with immunotherapy. Resistance to therapy, either in the form of primary refractory disease or evolutionary resistance, remains a significant issue in the treatment of NSCLC. Hence, predictive biomarkers and novel combinational strategies are required to improve the effectiveness and durability of treatment response 6for people with NSCLC. The aim of this study was to identify novel biomarkers and/or druggable proteins from deregulated protein networks within non-oncogene driven disease that are involved in the cellular response to cisplatin. Following exposure of NSCLC cells to cisplatin, in vitro quantitative mass spectrometry was applied to identify altered protein response networks. A total of 65 proteins were significantly deregulated following cisplatin exposure. These proteins were assessed to determine if they are druggable targets using novel machine learning approaches and to identify whether these proteins might serve as prognosticators of platinum therapy. Our data demonstrate novel candidates and drug-like molecules warranting further investigation to improve response to platinum agents in NSCLC.
Introduction
Lung cancer is a significant health care burden, accounting for 18.4% of all cancer-related deaths (1). The most commonly diagnosed form of this disease is non-small cell lung cancer (NSCLC) constituting 85% of all cases. Treatment options for NSCLC include front-line immune checkpoint inhibitors as monotherapy, particularly in patients with programmed death receptor ligand (PD-L1) expression in >50% of tumor cells detected by immunohistochemistry (2), or in combination with platinum-based chemotherapy in the advanced setting (2–4), maintenance immunotherapy for unresectable stage III disease after radical chemoradiotherapy (5) or adjuvant platinum-based chemotherapy for those with high risk resected early stage NSCLC (6). However, while survival rates are improving with advances such as screening, improved imaging technology and surgical techniques along with the use of in targeted therapeutics and immunotherapy agents (7), 5-year survival rates remain poor at 20.5% (8).
Cisplatin, cis-Diammineplatinum (II) dichloride, is the most widely studied platinum agent (9–13), and is currently one of the most effective chemotherapeutic agents available for treating NSCLC. However, objective response rates for patients with advanced NSCLC remain poor at 40%–50% (14, 15). Even in combination with immunotherapy the objective response rates remain between 35% and 70% pending the level of PD-L1 expression (2, 5, 16). The development of refractory disease is primarily due to intrinsic or acquired drug resistance which reduces the effectiveness of platinum agents used alone or in combination during in the management of NSCLC (17, 18). Thus, clinical trials are ongoing to identify novel combination strategies for use with traditional chemotherapy (19). Likewise, patients with resected stage II or III disease receive adjuvant platinum-based chemotherapy with the aim of reducing risk of future recurrence. In this setting, meta-analysis data suggest a benefit in terms of survival of 5% at 5 year with the addition of chemotherapy after a successful complete surgical resection. Many patients will still experience disease recurrence after this therapy suggesting a disease biology that is refractory to chemotherapy (6). Further strategies are warranted to amplify/complement the effects of platinum-based chemotherapy.
The recognized mechanism of action for cisplatin is via the induction of DNA damage. Cisplatin induces anti-tumor effects by predominantly binding to DNA or RNA to form DNA-platinum adducts (Stordal & Davey, 2007). These intra- or inter-strand crosslinks (ICLs) induce DNA damage response pathways by blocking the unwinding of double stranded DNA which is necessary for the cellular processes of transcription and DNA replication. Pathways implicated in adduct removal include the base excision and nucleotide excision repair (BER and NER) pathways and Fanconi anemia (FA) pathway (20, 21). The extent of damage or failure to remove these blockages and repair damaged DNA account for the cytotoxic effects of cisplatin resulting in tumor cell death. The DNA damage response (DDR) is also associated with resistance to cisplatin, with upregulation of these pathways promoting repair of DNA adducts and enabling tumor cell survival (22–24).
Alongside the DDR, over 147 mechanisms of platinum resistance have been suggested (25). These varied and complex mechanisms also include reduced uptake or increased drug efflux via copper transporters to reduce intracellular drug concentrations yielding acquired resistance (26–32). Recent unbiased screening (33–36) has highlighted platinum resistance in lung cancer arises via predominant cellular programs of the aforementioned DNA repair as well as upregulating expression of cell cycle-associated genes, transforming growth factor (TGF)-β signaling and apoptosis avoidance (22). Hence, blocking or exploiting the function of proteins within these pathways and cellular programs might serve as a useful approach to improve the sensitivity of tumors to platinum agents.
Novel strategies to enhance platinum-based therapy response rates will rely upon identifying those people with NSCLC who will best respond to platinum-based therapy, as well as providing alternate treatment options for those less likely to benefit from platinum-based regimens. As such, employing predictive biomarkers and introducing synergistic therapies might improve the effectiveness of therapy by reducing unnecessary treatment or doses of therapy. The aim of this study was to identify protein networks involved in the cellular response to cisplatin by mining quantitative mass spectrometry data. Our study focused on cisplatin response in non-oncogene driven disease where platinum-based chemotherapy is employed in the absence of targeted therapies. We sought to examine whether these deregulated proteins identified in our study were druggable using novel machine learning approaches and to identify whether these proteins might serve as prognosticators of platinum therapy. Based upon our approach, these data demonstrate there are several novel protein candidates and drug-like molecules that warrant further investigation. With additional preclinical evaluation and future clinical investigations, our findings may lead to potentially improved responses to platinum agents.
Materials and Methods
Antibodies and Reagents
Antibodies against aldehyde dehydrogenase 3A1 (ALDH3A1); ab76976) and TP53I3 (ab64798) were purchased from Abcam. The α-tubulin antibody (T9026) was purchased from Sigma Aldrich. Donkey anti-rabbit and anti-mouse Alexa Fluor 680 and 800 antibodies were purchased from Life Technologies. Complete EDTA-free protease inhibitor mixture was from Roche Applied Sciences and phosphatase inhibitor cocktail Cell Signaling Technology. cis-Diammineplatinum (II) dichloride (cisplatin) was purchased from Merck.
Cell Culture and Treatment
All cell lines were purchased from the American Type Culture Collection (ATCC) and maintained at 37°C in a humidified atmosphere containing 5% CO2 in RPMI 1640 medium supplemented with 1% L-Glutamine and 1% non-essential amino acids (NEAA). Cell culture media were supplemented with 10% (v/v) heat-inactivated fetal bovine serum (FBS).
For differential mass spectrometry-based quantitative proteomics, H460 cells were treated cisplatin (7.5 µM) for 24 h. For in vitro drug treatments prior to western blot analysis, cells were treated with cisplatin (5 µM) for 12 h.
Nano Liquid Chromatography Mass Spectrometry Analysis
MS analysis was performed using an AB Sciex 5600+ TripleTOF mass spectrometer interfaced to an Ekspert™ NanoLC system. Sample preparation and library construction for differential quantitative liquid chromatography (LC)-MS/MS was performed and analyzed as previously described (37). Briefly, equal amounts of sample in triplicate were prepared by filter-aided sample preparation (FASP) (38) and tryptic digested peptides were analyzed by LC-MS/MS. For each biological replicate, a spectral library was generated with a traditional data- dependent acquisition (DDA) approach. The DDA data was searched against a human proteome library (UP000005640; UniProt.org) using ProteinPilot 5.0 (SCIEX). Following generation of the spectral library, the identification and quantification of cisplatin-treated and untreated H460 triplicate cell lysate derived peptides were quantified using an LC-MS/MS approach known as variable sequential window acquisition of all theoretical mass spectra (SWATH-MS) with a targeted data extraction strategy to mine the resulting fragment ion data set. SWATH-MS transition, peptide, and protein level summarization was performed using Skyline (39) by normalization at the peptide level using median normalization and Tukey median polish for summarization over all features in a run.
Log2 transformation was performed prior to further statistical analysis. Differentially regulated proteins were identified by applying empirical Bayes moderated t-statistics tests (40, 41) implemented in the limma package using R statistical environment (version 3.5.2) and the Benjamini-Hochberg correction (42) was applied to control the false discovery rate [FDR (or q-value)]. List of all quantified proteins (including differentially regulated proteins) by LC-MS/MS are listed in Supplemental Table 1.
Functional Analysis
Statistically significant proteins from the comparison of cisplatin treated cells versus untreated cells at a q-value threshold of ≤ 0.1 (65 proteins) were used for gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG) and Pathway over representation analysis. For protein network analysis, all proteins at a q-value threshold ≤0.1 were selected for StringDB analysis. Network was generated with default confidence (medium confidence, 0.4) and evidence from active interaction sources included, text mining, experiments, databases, co-expression, neighborhood, gene fusion, and co-occurrence.
Signaling network analysis was conducted using Ingenuity® Pathway Analysis (IPA) (QIAGEN). Right-tailed Fisher’s exact test was used to determine over representation probability of canonical pathways in the protein dataset. Z-score analysis was also performed in IPA to predict activation or inhibition of significant canonical pathways based on Log fold-change values of the input molecules from the quantitative proteomics data set. All IPA canonical pathways are listed in Supplemental Table 2.
GO biological processes and KEGG enrichment analyses were performed using the ClueGo (ver. 2.5.6) (43) plugin in Cytoscape (ver. 3.7.2). Functionally grouped and annotated networks were assessed by hypergeometric testing (enrichment/depletion, two-sided hypergeometric test and Bonferroni step down correction). Networks were visualized using the default ClueGo layout with functionally grouped network, GO/KEGG terms as nodes and linked based on their kappa score level. Percentage of mapped genes belonging to each term were also represented. Node size is representative of the term enrichment significance. A kappa coefficient of 0.4 was used as a threshold values while redundant groups with >50% overlap were merged.
Collection of Lysates and Western Blot Analysis
Prior to western blot analysis, cells were washed once with ice-cold phosphate buffered saline and lysed in lysis buffer (50 mM HEPES (pH 7.5), 150 mM KCl, 5 mM EDTA, 0.05% IGEPAL CA-630 (v/v), 1x protease inhibitor cocktail and 1x phosphatase inhibitor cocktail). Lysates were briefly sonicated before low speed centrifugation (900 g, 10 min at 4°C) with total protein concentration determined by Bicinchoninic Acid (BCA) protein assay. Equal amounts of lysate were denatured in 1x Laemmli buffer containing 8% β-mercaptoethanol by heating for 5 min at 80°C. Samples were subjected to SDS-PAGE analysis using Bolt 4%–12% Bis-Tris Plus pre-cast gels (Life Technologies) before protein transfer to nitrocellulose membrane (GE Healthcare Life Sciences) using the Novex system (Life Technologies). Membranes were blocked in Odyssey blocking buffer (Li-Cor Biosciences) and incubated with primary antibodies overnight at 4°C with gentle agitation. All antibodies were used at a 1:1,000 dilution with the exception of the anti-α-tubulin antibody (1:5,000). Membranes were washed with PBS containing 0.1% Tween-20, incubated with appropriate secondary antibodies and scanned using an Odyssey CLx imaging system (Li-Cor Biosciences). Images were subjected to densitometric analysis using ImageJ software.
Prediction of Druggable Proteins
The druggability and developability space is rapidly evolving in the post-structural genomics era. The development of computational methods using sequence, protein structure and ligand (including Lipinski's rule of five) has helped characterize protein targets and describe strategies for the optimal integration of protein druggability data (44). As an initial step, a recently described machine learning prediction model was utilized for the extrapolation of potentially druggable proteins from the upregulated protein list (45). This machine learning approach was derived from a set of 666 known druggable as well as 219 non-druggable protein sequences. In our study, all sequences for up and downregulated proteins identified by quantitative mass spectrometry were collected in FASTA format. Druggable features of these proteins were evaluated using available python scripts (https://github.com/muntisa/machine-learning-for-druggable-proteins). Briefly, these scripts applied thirteen machine learning classification types which were used to learn and recognize characteristic traits for the prediction and classification of proteins into druggable categories. Prediction performance was evaluated against several breast cancer proteins, cancer-driving proteins and RNA-binding proteins. Druggable and non-druggable sequences were given the labels 1 and 0 respectively; therefore, scores closest to 1 are proteins with greatest druggable potential and are more likely to be a promising drug target candidate. This approach allows unbiased prediction of proteins that are less studied with limited information about their function. Moreover, this approach can be applied to predict druggability of proteins with unresolved or incomplete structure which would otherwise affect other approaches such as the structure-based druggability prediction method.
To complement the machine learning modelling, two additional approaches were used to assess protein structural druggability potential and identify compounds known to target the deregulated proteins. First, proteins with three-dimensional models and ligandable cavities were evaluated using the CanSAR knowledgebase (46) to confirm druggability potential. The previous studies revealed that cancer proteins tend to interact more strongly than other categories of proteins including essential and control proteins in the human interactome (47), therefore CanSAR network-based protein-protein interactions data can serve as potential therapeutic targets in cancer among certain important pro-carcinogenic signaling pathways. The Binding Database was also examined to identify investigational and/or approved chemical entities from the literature for a druggable target (48). All identified compounds are listed in Tables 1 and 2 with additional compounds and structures listed in Supplemental Table 3.
Bioinformatics and Survival Analysis
Survival analysis was performed on Illumina gene expression data and clinical information with the University of Texas Lung Specialized Program of Research Excellence (UT Lung SPORE) dataset and obtained from the GEO database (GSE42127) (49). Within this cohort, 127 patients were observed following curative resection (OBS) whereas 49 patients received adjuvant platinum-based chemotherapy (ACT). Illumina probes corresponding to the identified differentially regulated proteins are listed in Supplemental Table 4. High versus low transcript expression was stratified according to median expression with OBS and ACT overall survival curve comparison performed by the Kaplan-Meier method and log-rank test where P-value > 0.05 was considered significant.
Results
Cisplatin-Induced Regulation of NSCLC Proteins
To identify proteins that might be exploited to improve platinum response, we sought to identify proteins from NSCLC cell lines that were deregulated upon acute cisplatin exposure. We reasoned that inhibiting the function of these identified deregulated proteins which, while not considered classical targets, but may serve to promote cancer cell survival following cisplatin exposure, could improve therapy response. By following our earlier approach, the H460 NSCLC cell line was treated with the appropriate IC50 concentration of cisplatin. This cell line was selected as a commonly utilized non-squamous cell line for in vitro assays. Whole cell lysates were collected from untreated and cisplatin treated cells and subjected to quantitative SWATH-MS to identify differentially expressed proteins. Of the 1081 proteins robustly identified, 430 proteins were upregulated following cisplatin exposure, while 586 proteins were downregulated (Supplemental Table 1). A total of 65 differentially regulated proteins were statistically significant with a q-value of less than 0.1 (Figure 1A). As shown in Figure 1B, cisplatin exposure induced the upregulation of 26 of these proteins and the downregulation of 39 proteins. The top three proteins based on log2 fold change and FDR significance were the ALDH3A1, TP53I3 and ferredoxin reductase (FDXR) upregulated proteins (Figure 1C), whereas the top three downregulated proteins were sulfiredoxin-1 (SRXN1), heat shock protein 90α (HSP90AA1) and phosphoglycerate dehydrogenase (PHGDH) (Figure 1D).
Figure 1 Identification of proteins deregulated by cisplatin exposure. (A) Pie chart showing proportion of identified proteins up or downregulated following cisplatin exposure and those deregulated proteins considered significant (q-value < 0.1). (B) Volcano scatter plot of log2 fold protein changes (cisplatin treatment versus untreated) ranked by significance (-log10 P value) with proteins of q-value < 0.1 highlighted in black. (C) List of top three cisplatin-induced upregulated proteins ranked by Q-value. (D) List of top three cisplatin-induced downregulated proteins ranked by Q-value. (E) Representative endogenous ALDH3A1 and TP53I3 western blot analysis from lysates of NSCLC cell lines treated in the absence or presence of cisplatin. α-Tubulin used as loading control. (F) Densitometry quantification of (E), with dot points representing average log2 of relative ALDH3A1 levels from three independent experiments with lines connecting respective untreated and cisplatin treated cell lines (paired Student’s t Test, two tailed **P = 0.0081). (G) Densitometry quantification of (E), with dot points representing average log2 of relative TP53I3 levels from three independent experiments with lines connecting respective untreated and cisplatin treated cell lines (paired Student’s t Test, two tailed **P = 0.0014).
Based on the SWATH-MS differentially regulated protein candidates, we sought to orthogonally validate our findings from the single H460 cell line by using western blot analysis and a wider panel of NSCLC cell lines that were treated with or without cisplatin. The protein levels of ALDH3A1 and TP53I3 were selected for examination by western blot analysis given our access to commercially available antibodies and that these proteins were identified as significantly upregulated by cisplatin. Consistent with the SWATH-MS, cisplatin significantly induced the upregulation of ALDH3A1 in all cell lines where the protein was detectable (p = 0.0081; Figures 1E, F). Notably, ALDH3A1 expression was not detected in the squamous cell carcinoma cell lines EBC-1, SK-MES-1, and H1703. Consistently, protein levels for TP53I3 were also significantly upregulated following cisplatin exposure in all cell lines evaluated (p = 0.0014; Figures 1E, G). Taken together, these data corroborate with the differential global quantitative MS profile obtained using H460 treated cells and suggest, that at least for ALDH3A1 and TP53I3, modulation of these protein levels are a common in vitro response to cisplatin exposure.
Comparative Pathway Analysis
To examine the functional importance of deregulated proteins following cisplatin exposure, all q-value ≤ 0.1 proteins were subjected to STRINGdb analysis. This database was examined to evaluate potential protein-protein interactions between identified proteins (50). These data indicated there to be clusters of networks existing between 77% of the significantly deregulated proteins following cisplatin treatment (Supplemental Figure 1). Next, we sought to identify the function for these proteins by performing pathway analysis using IPA. A total of 141 canonical pathways were identified (Supplemental Table 2) with 39 of these pathways reaching the significance threshold (p < 0.05; Supplemental Table 2). Of the top 10 canonical pathways ranked by significance, an activation z-score was calculated for three pathways (Figure 2A). EIF2 signaling was predicted to be inhibited following cisplatin exposure (z-score = -2.83) whereas xenobiotic metabolism was predicted to be activated in response to cisplatin (constitutive androstane receptor (CAR) signaling and pregnane X receptor (PXR) signaling pathways with respective z-scores of 0.45 and 2.0). These findings are consistent with the observations that cisplatin inhibits eukaryotic translation (51) and, like other chemotherapeutic agents, induces the upregulation of drug metabolism pathways (52, 53).
Figure 2 Pathway overrepresentation and gene ontology analysis of cisplatin-induced significantly deregulated proteins. (A) Top 10 canonical pathways identified by Ingenuity pathway analysis associated with cisplatin response ranked by significance (-log10 P value) and predicted pathway activation or inhibition indicated by Z-score. (B) Gene ontology enrichment analysis assessed by Gene ontology (GO) biological processes with functional networks visualized using ClueGo. GO biological process terms listed as nodes with genes belonging to each term represented. Node size represents the term enrichment significance. (C) Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses assessing gene ontology with functional networks visualized using ClueGo. KEGG terms listed as nodes with genes belonging to each term represented. Node size represents the term enrichment significance. (D) Top five molecular and cellular functions identified by Ingenuity pathway analysis (IPA) ranked by P value range and listing the number of molecules identified.
The cisplatin-induced differentially regulated proteins were also classified using over representation analysis to enrich for clusters based on protein ontology. As shown in Figure 2B, protein function identified by GO biological process indicated predominant clusters involved in DNA replication, disulphide oxidoreductase activity, tRNA biosynthesis and signal-recognition particle (SRP)-dependent regulation of translation. To further evaluate biological function of the differentially regulated proteins, we also applied KEGG enrichment analysis. As shown in Figure 2C, several pathways were identified consistent with the GO analysis including ribosome regulation, DNA replication and tRNA biosynthesis. In addition, protein clusters involved in ferroptosis, drug metabolism and cysteine and methionine metabolism pathways were identified. Supporting the GO and KEGG enrichment analysis, IPA analysis of cellular and molecular function identified the top pathway functions as RNA damage and repair, protein synthesis and pathways involved in regulating cellular proliferation, survival and development (Figure 2D). Notably, IPA analysis also identified p53 as the top upstream regulator linking this protein with modulation of tumor cell proliferation (Supplemental Figure 2). Taken together, the top cisplatin-induced differentially regulated proteins modulate translation, replication, pathways exposure induced pathways regulating DNA damage, translation and metabolic control.
Identification of Druggable Proteins
Having identified biological pathways that are impacted in response to cisplatin exposure, we postulated that the significantly differentially regulated proteins might prove exploitable as potential therapeutic targets to enhance chemotherapy sensitivity. The differentially regulated proteins were separated into those proteins upregulated (Table 1) or downregulated (Table 2) in response to cisplatin exposure. We applied a recently reported machine learning approach (45) to identify which proteins could be inhibited with drugs/ligands by having possible targetable regions within the protein based upon their protein sequence. This approach predicts druggable proteins by analyzing the primary protein sequence with confidence scores close to 1.0 considered druggable. The top 10 predicted druggable targets identified as significantly upregulated by our proteomics are listed in Table 1. Validating the predictive machine learning approach, structural ligandability was also explored based on CanSAR database. Structural ligandability refers to identification of protein pockets with the 3D structures of target proteins or domains and are summarized in Supplemental Table 3. Of these proteins, structures of full-length or certain domains have been solved, for eight proteins. Besides machine learning and structural ligandability, we also compiled ligand information from CanSAR and BindingDB for these targets wherein four of the top 5 proteins have identified compounds which target these proteins. For example, the anti-seizure drug lacosamide is indicated to specifically bind dihydropyrimidinase-related protein 2 (DPYSL2) (54). Patent literature also exists for the drugs identified that target the proteins DPYSL2 and ALDH3A1. We have provided datasets which includes structure, literature reference and activity data compiled from BindingDB as supporting information for the top ranked up- and downregulated druggable proteins.
As shown in Table 2, the top 10 predicted druggable targets were also identified for the proteins downregulated in response to cisplatin exposure. Like the upregulated proteins, all but four of the proteins predicted to be druggable had a solved protein structure. Moreover, compounds have been identified that target three of the top 10 proteins, alanyl-tRNA synthetase (AARS), ubiquitin carboxy-terminal hydrolase L5 (UCHL5) and tyrosyl-tRNA synthetase (YARS). Unlike the upregulated proteins, these compounds do not have patent coverage. Besides the protein structure and ligand based approaches, protein-protein interactions form signaling nodes and hubs to transmit pathophysiological cues along molecular networks, thereby promoting cancer progression, invasion, and/or metastasis (55). We report the number of signaling proteins in the interactome curated from CanSAR (Tables 1 and 2) for the up and downregulated proteins. Disruption of protein-protein interactions is critical for cancer and offers a novel and effective strategy to curtail the transmission of oncogenic signals. Taken together, by employing these machine learning, structure and ligand-based approaches, we have identified novel therapeutic targets that are regulated by cisplatin. Further investigation is required to identify specific and potent drug-like molecules with therapeutic potential and determine the combinational potential of these compounds with chemotherapy.
Prognostic Benefit of Cisplatin-Induced Deregulated Factors
We next sought to determine whether the proteins identified as differentially regulated had prognostic benefit to identify people with NSCLC who would benefit from platinum-based chemotherapy. To investigate prognostic potential, we examined transcript expression for each significantly differentially regulated protein across microarray profiling of 176 evaluable NSCLC cases from the UT Lung SPORE cohort (GEO database GSE42127). In this cohort, NSCLC cases underwent curative resection and either received adjuvant platinum-based chemotherapy (ACT) or observation (OBS) alone. Cases were stratified based on median expression for each gene and survival analysis performed using the Kaplan-Meier method. The top 10 transcripts of upregulated proteins ranked by prognostic significance are listed in Figure 3A. Elevated levels of four of the top five prognostic transcripts had a hazard ratio (HR) of less than one suggesting that in these cases, patients benefited from adjuvant chemotherapy (TP53I3, ALDH3A1, FK506 binding protein 10 (FKBP10) and integrin beta-1 (ITGB1); Figure 3A). Interestingly, lower levels of five of the top 10 transcripts, including DPYSL2, transmembrane protein 205 (TMEM205), cytokeratin-18 (KRT18), cathepsin L (CTSL) and glutathione S-transferase M3 (GSTM3), were prognostic for patients treated with adjuvant chemotherapy Of the upregulated proteins, elevated TP53I3 levels were the most significant prognostic transcript (HR = 0.41, 95% CI: 0.03–0.59, P = 0.002) with adjuvant chemotherapy having little benefit in those cases with lower levels (HR = 1.31, 95% CI: 0.55–3.11, P = 0.54).
Figure 3 Identification of factors with prognostic benefit for treatment of NSCLC patients with adjuvant chemotherapy. (A) Forest plot showing hazard ratio (HR) and 95% confidence interval (CI) determined from overall survival Kaplan-Meier analysis of transcripts for proteins identified as upregulated following cisplatin exposure. Cases stratified by median transcript expression and compared according to those receiving adjuvant chemotherapy (ACT) or observation alone (OBS) and ranked P value determined by log-rank method. (B) Forest plot determined from Kaplan-Meier analysis as in (A), comparing ACT versus OBS patients stratified by median expression of transcript for proteins identified as significantly downregulated following cisplatin exposure.
The prognostic benefit for transcripts of those proteins downregulated by cisplatin exposure was also investigated. As shown in Figure 3B, lower levels of four of the top 10 prognostic transcripts, including stathmin (STMN1), 40S ribosomal protein S20 (RPS20), caprin-1 (CAPRIN1) and 60S ribosomal protein L35 (RPL35), were associated with patients that benefited from adjuvant chemotherapy with a HR value of less than one. Elevated transcripts of inosine-5′-monophosphate dehydrogenase 2 (IMPDH2), SRXN1, endothelial differentiation-related factor 1 (EDF1), phosphoserine aminotransferase (PSAT1), Clu1/CluA homologue (CLUH) and glycyl-tRNA synthetase 1 (GARS1) were significantly associated with patients that benefited with adjuvant chemotherapy. Of the downregulated proteins, lower transcript levels of STMN1 had the most significant prognostic benefit for patients undergoing adjuvant platinum-based chemotherapy (HR = 0.08, 95% CI: 0.01–0.62, P = 0.002). Taken together, by assessing the prognostic value for those differentially regulated proteins, we have identified transcripts for these proteins that, at least in this patient cohort, have biomarker potential for improving platinum-based chemotherapy response.
Discussion
Novel therapeutic strategies combined with predictive biomarkers of platinum-based chemotherapy response are needed to transform the clinical management of NSCLC. Such strategies are key to selecting the right patients for platinum-based chemotherapy and for overcoming therapy resistance which continues to affect health outcomes for people with this disease. We sought to identify potentially exploitable proteins that have a biological function in response to cisplatin exposure. Although not mutated and considered classically targetable such mutant epidermal growth factor receptor (EGFR), inhibiting the function of deregulated proteins could sensitize to platinum-based therapy. For this reason, we applied a SWATH-MS proteomic approach to identify exploitable proteins and biological pathways that are differentially regulated following cisplatin exposure. Previous screening using genetic-based approaches have confirmed mechanisms of action for cisplatin and its platinum analogues. For example, the anti-tumor effects of cisplatin and carboplatin result from the DNA damage response whereas oxaliplatin kills cells via ribosome biogenesis stress (56). Our quantitative analysis at the protein level in lung cancer revealed that alongside significant downregulation of ribosomal and translation control proteins, cisplatin negatively impacts DNA replication and amino acid synthesis while inducing upregulation of metabolic pathways. By further evaluating available clinical data and applying machine learning approaches, we identified potentially prognostic markers (Figure 3) and novel druggable targets (Tables 1 and 2) from the differentially regulated pathways which might prove actionable to enhance the effectiveness of platinum-based chemotherapy. While at an early exploratory stage, our study provides a platform for future preclinical evaluation. Our work will also complement ongoing and future drug discovery and biomarker development programs to ultimately improve future health outcomes for people with NSCLC.
Our proteomics analysis identified that cisplatin treatment yielded a marked downregulation of translation and protein synthesis. This is consistent with prior work indicating cisplatin, but not carboplatin or oxaliplatin negatively affect translation in in vitro systems (51). Of the proteins associated with translational control, five were found to be potentially druggable, potentially prognostic for chemotherapeutic response, or both: AARS, YARS, 60S ribosomal protein L7 (RPL6), 60S ribosomal protein L7 (RPL7), and 60S ribosomal protein L27a (RPL27a). AARS and YARS are involved in tRNA charging whilst RPL6, RPL7, and RPL27a are involved in EIF2 signaling for regulating both global and specific mRNA translation. AARS and YARS function to catalyze the alanine or tyrosine to tRNA attachment, respectively, in a two-step reaction: first by activating alanine or tyrosine via ATP to form Ala-AMP or Tyr-AMP, then transferring the activated alanine or tyrosine to tRNA. Not much is known about the role of AARS in lung cancer or any other solid malignancies. However, YARS is elevated in malignant gastric tissue compared to adjacent non-malignant tissue. The knockdown of YARS results in repressed proliferation and invasion in gastric cancer whilst enhancing apoptosis. The upregulation of YARS yields the reverse effects (57). RPL6, RPL7, and RPL27a are all components comprising the large ribosomal subunit. Increased expression of RPL6 promotes G1 to S phase transition of gastric cancer cells and leads to accelerated growth (58). RPL27a suppression leads to delayed tumorigenesis in vivo (59). Whilst little is known about the involvement of RPL7 in lung cancer, it is involved in microsatellite instability in colorectal cancer (60). Taken together, although translation control is predicted to be inhibited following cisplatin, the malignant function for these proteins alongside our analysis suggest the identified proteins might serve as therapeutic options.
Further pathways identified by our analysis were proteins associated with cell death and survival. Of these proteins, HSP90AA1, ITGB1, cytokeratin-8 (KRT8) and receptor for activated C-kinase (RACK1) are identified prognostic markers for overall survival in solid malignancies including lung cancer. HSP90AA1 and KRT8 associate with a poorer rate of survival for people with NSCLC whereby higher protein expression is predictive of poorer prognosis and decreased survival (61, 62). Higher levels of KRT8 also significantly associate with tumor progression in NSCLC (62). ITGB1 and RACK1 are significant prognostic markers for survival in NSCLC and early-stage NSCLC, respectively (63, 64). These proteins are also prognostic markers in other solid malignancies. In gastric cancer, higher KRT8 expression also associates with poorer prognosis (65). Similarly, high ITGB1 expression associates with poor patient prognosis and independently correlates with shortened overall survival and shortened disease-free survival in people with colorectal cancer (66). In people with breast cancer, RACK1 is also a significant prognostic marker for survival (67). Indeed, our study identified elevated levels of ITGB1 as a significant prognosticator for overall survival in people treated with platinum-based therapy. Moreover, our findings indicate that ITGB1, alongside KRT8 and RACK1, have druggable potential. As such, further investigation is required to determine the diagnostic or therapeutic potential to improve platinum response.
Another key mechanism identified by our study is the predicted activation of cisplatin-induced drug metabolism. Our findings reveal the drug metabolism enzymes ALDH3A1, GSTM3 and glutathione S-transferase K1 (GSTK1) as significantly upregulated by cisplatin. Expression of drug detoxification and xenobiotic metabolizing enzymes are linked with patient response to chemotherapy and the onset of drug resistance by clearing reactive and cytotoxic compounds (52, 53). For example, ALDH3A1 activity is reported to increase upon cisplatin treatment, suggesting that this enzyme can metabolize cytotoxic aldehydes and confer resistance to cisplatin (68–70). Activity of glutathione transferases (GST), particularly GSTP1, is also associated with cisplatin sensitivity by detoxifying cisplatin to form an inactive cisplatin-glutathione conjugate (71, 72). Of the enzymes identified as deregulated, ALDH3A1 has druggable potential in NSCLC. While this might prove useful to improve platinum-based chemotherapy, targeting isoforms of the ALDH superfamily has been proposed for therapy resistant cancers (73). In addition, our findings also suggest stratifying patients according to ALDH3A1 or GSTM3 levels might prove beneficial. Indeed, in the case of GSTM3, a null genotype yielding low expression, is reportedly associated with reduced lung cancer risk (74). The enzymes identified in our study warrant further investigation in NSCLC as potential biomarkers and as druggable targets to improve patient chemotherapy response.
Our study has also identified the deregulated proteins CLUH and 3′-phosphoadenosine 5′-phosphosulfate (PAPS) synthetase 2 (PAPSS2) with other functional roles in metabolic regulation. CLUH is involved in mitochondrial biogenesis and distribution with genetic depletion of the protein resulting in loss of mitochondrial enzymes yielding oxidative phosphorylation defects and a dysfunctional Krebs cycle (75–77). Radio and chemotherapy resistant cancers, at least in triple negative breast cancer (78) and acute myeloid leukemia (79), are suggested to rely upon oxidative phosphorylation for survival and progression. As such, targeting CLUH or other proteins regulating oxidative phosphorylation might enhance sensitivity to therapy in lung cancer (80). The other protein we identified was PAPSS2. PAPSS2 is an isoform of PAPSS1, which synthesizes PAPS from ATP and inorganic sulphate (81, 82). While further investigation for PAPSS2 is required, depletion of the PAPSS1 isoform sensitized NSCLC to cisplatin suggesting that targeting sulfation reactions might improve activity of cisplatin or other DNA damaging agents (83). In line with these findings, our study indicated that both CLUH and PAPSS2 exhibited high confidence as potential drug targets.
While we have identified deregulated proteins as potential biomarkers or therapeutic targets with roles in drug metabolism, translational control, cell death, and survival, and DNA damage there were deregulated proteins which have no reported association with chemotherapy response or extensively reported link with solid malignancies. Those with no known association with chemotherapy or chemoresistance include heterogenous nuclear ribonucleoprotein A/B (HNRNPAB), nuclear distribution gene C (NUDC), TP5313, transferrin receptor 1 (TFRC), FDXR, and PAPSS2. The only protein found to have no known role in lung cancer or other solid malignancies was AARS, as mentioned previously. Of these proteins, only TP53I3 had prognostic potential whereas each of these proteins were identified as potentially targetable. Although further analysis is required to determine the precise role of these proteins, it remains possible that targeting these proteins alongside chemotherapy might enhance platinum agent activity.
Our exploratory study has applied a SWATH-MS proteomic based approach to identify proteins deregulated following cisplatin exposure which might serve as potential biomarkers or therapeutic targets in NSCLC. It is worth noting that our proteomics approach was not able to detect an entire proteome. Moreover, another limitation of our study was the use of a single cell line for our proteomic analysis. Given that the H460 cell line is not representative of all non-squamous NSCLC, further proteomics analysis of laboratory models of disease or clinically relevant samples would strengthen our findings. One such approach might be to also incorporate the comparison of in vitro models and ex vivo patient tissues pre- and post-therapy, to validate our findings and further identify novel druggable proteins. However, we note that validation of our SWATH-MS by western blot analysis and usage of a wider cell line panel confirmed our single cell line proteomics findings. These data point to the possibility that we have identified a common cisplatin response network, at least for some of the identified significantly deregulated proteins. Nonetheless, our study provides a rationale for further investigation of these deregulated and potentially exploitable proteins. Indeed, while the findings of this study will not immediately impact clinical decision making, the work provides a platform and impetus for further investigation and future drug development. Analysis of further cohorts is also necessary to determine the predictive potential of the biomarkers identified in our study. Moreover, further research is required to confirm the therapeutic potential of the identified proteins and to develop small molecules to inhibit them. Overall, future inhibition of, or testing biopsies for these proteins could be useful to enhance platinum-based chemotherapy response and improve health outcomes for people impacted by NSCLC.
Data Availability Statement
The proteomics data presented in the study are deposited in the Paranorma Public repository, ProteomeXchange ID: PXD024209 (https://urldefense.com/v3/__https://panoramaweb.org/lGlNod.url__;!!NVzLfOphnbDXSw!RlKUaKl0OWpX2gD9cUB6_1rO--VMvbjHMnHX7KSfhgqrwfhpkVfx6TIais4-k4mrV5Y$).
Author Contributions
MA, KO’B, and DR conceived and designed the study. MA supervised the study. S-LR, KD, SB, KS, CM, CO’L, ES, EB, A-MB, and MA performed the experiments and interpreted the data. MG, MM, NG, and EB performed the in silico analyses and bioinformatics. S-LR, KD, and MA analyzed the data. DR, KO’B, and MA wrote the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare competing financial interests; KO’B and DR are founders of CARP Pharmaceuticals. EB, DR, and KO’B are founders of Carpe Vitae Pharmaceuticals. MA, EB, KO’B, and DR are inventors on patent applications filed by Queensland University of Technology.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank members of the Cancer and Ageing Research Program for critical review of the manuscript. The work was supported by a Senior Clinical Research Fellowship (to KO’B), the Princess Alexandra Hospital (to KO’B), project grant funding from the Australian National Health and Medical Research Council and Cure Cancer Australia through the Cancer Australia Priority-driven Cancer Research Scheme (to MA) and an International Association for the Study of Lung Cancer Foundation (IASLC) Foundation Award (to MA and KO’B). NG acknowledges funding support from Advance Queensland Industry Research fellowship. We acknowledge computing time from the QUT high-performance computing facility and the NCI, Canberra.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.615967/full#supplementary-material
Supplementary Figure 1 | StringDB network analysis of proteins identified by quantitative mass spectrometry as significantly deregulated following cisplatin exposure. Line thickness is indicative of confidence in evidence linking proteins within a network.
Supplementary Figure 2 | Identification of TP53 by IPA analysis as the predicted top upstream regulator of differentially regulated protein.
References
1. 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:394–424. doi: 10.3322/caac.21492
2. Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csőszi T, Fülöp A, et al. Pembrolizumab versus Chemotherapy for PD-L1-Positive Non-Small-Cell Lung Cancer. N Engl J Med (2016) 375:1823–33. doi: 10.1056/NEJMoa1606774
3. Paz-Ares L, Luft A, Vicente D, Tafreshi A, Gümüş M, Maziéres J, et al. Pembrolizumab plus Chemotherapy for Squamous Non-Small-Cell Lung Cancer. N Engl J Med (2018) 379:2040–51. doi: 10.1056/NEJMoa1810865
4. Gandhi L, Rodríguez-Abreu D, Gadgeel S, Esteban E, Felip E, De Angelis F, et al. Pembrolizumab plus Chemotherapy in Metastatic Non-Small-Cell Lung Cancer. N Engl J Med (2018) 378:2078–92. doi: 10.1056/NEJMoa1801005
5. Antonia SJ, Villegas A, Daniel D, Vicente D, Murakami S, Hui R, et al. Durvalumab after Chemoradiotherapy in Stage III Non-Small-Cell Lung Cancer. N Engl J Med (2017) 377:1919–29. doi: 10.1056/NEJMoa1709937
6. Pignon J-P, Tribodet H, Scagliotti GV, Douillard J-Y, Shepherd FA, Stephens RJ, et al. Lung Adjuvant Cisplatin Evaluation: A Pooled Analysis by the LACE Collaborative Group. J Clin Oncol (2008) 26:3552–9. doi: 10.1200/jco.2007.13.9030
7. Zhang C, Leighl NB, Wu Y-L, Zhong W-Z. Emerging therapies for non-small cell lung cancer. J Hematol Oncol (2019) 12:45–5. doi: 10.1186/s13045-019-0731-8
8. National Cancer Institute. Cancer Stat Facts: Lung and Bronchus Cancer (2019). Available at: https://seer.cancer.gov/statfacts/html/lungb.html.
9. Fennell DA, Summers Y, Cadranel J, Benepal T, Christoph DC, Lal R, et al. Cisplatin in the modern era: The backbone of first-line chemotherapy for non-small cell lung cancer. Cancer Treat Rev (2016) 44:42–50. doi: 10.1016/j.ctrv.2016.01.003
10. Dasari S, Tchounwou PB. Cisplatin in cancer therapy: molecular mechanisms of action. Eur J Pharmacol (2014) 740:364–78. doi: 10.1016/j.ejphar.2014.07.025
11. Nakashima K, Murakami H, Omori S, Wakuda K, Ono A, Kenmotsu H, et al. Doublet chemotherapy with cisplatin and pemetrexed is associated with a favorable outcome in patients with advanced non-squamous non-small-cell lung cancer who are eligible for bevacizumab and maintenance therapy. Mol Clin Oncol (2016) 5:575–8. doi: 10.3892/mco.2016.1001
12. Mai L, Luo M, Wu JJ, Yang JH, Hong LY. The combination therapy of HIF1alpha inhibitor LW6 and cisplatin plays an effective role on anti-tumor function in A549 cells. Neoplasma (2019) 2019:776–84. doi: 10.4149/neo_2018_180921N708
13. Zhang YQ, Jiang LJ, Jiang SX, Xu YF, Zhou BB, Huang GH, et al. Gefitinib with or without Transarterial Infusion Chemotherapy (Cisplatin) for Large Nonsmall Cell Lung Cancer with Epidermal Growth Factor Receptor Mutations. J Vasc Interventional Radiol JVIR (2019) 30(7):1004–12. doi: 10.1016/j.jvir.2018.12.705
14. Galluzzi L, Senovilla L, Vitale I, Michels J, Martins I, Kepp O, et al. Molecular mechanisms of cisplatin resistance. Oncogene (2012) 31:1869–83. doi: 10.1038/onc.2011.384
15. Lee HY, Mohammed KA, Goldberg EP, Kaye F, Nasreen N. Cisplatin loaded albumin mesospheres for lung cancer treatment. Am J Cancer Res (2015) 5:603–15.
16. Socinski MA, Jotte RM, Cappuzzo F, Orlandi F, Stroyakovskiy D, Nogami N, et al. Atezolizumab for first-line treatment of metastatic nonsquamous NSCLC. New Engl J Med (2018) 378:2288–301. doi: 10.1056/NEJMoa1716948
17. Makarova KS, Wolf YI, Alkhnbashi OS, Costa F, Shah SA, Saunders SJ, et al. An updated evolutionary classification of CRISPR-Cas systems. Nat Rev Micro (2015) 13:722–36. doi: 10.1038/nrmicro3569
18. MacDonagh L, Gallagher MF, Ffrench B, Gasch C, Breen E, Gray SG, et al. Targeting the cancer stem cell marker, aldehyde dehydrogenase 1, to circumvent cisplatin resistance in NSCLC. Oncotarget (2017) 8:72544–63. doi: 10.18632/oncotarget.19881
19. Lazzari C, Karachaliou N, Bulotta A, Viganó M, Mirabile A, Brioschi E, et al. Combination of immunotherapy with chemotherapy and radiotherapy in lung cancer: is this the beginning of the end for cancer? Ther Adv Med Oncol (2018) 10:1758835918762094–1758835918762094. doi: 10.1177/1758835918762094
20. Slyskova J, Sabatella M, Ribeiro-Silva C, Stok C, Theil AF, Vermeulen W, et al. Base and nucleotide excision repair facilitate resolution of platinum drugs-induced transcription blockage. Nucleic Acids Res (2018) 46:9537–49. doi: 10.1093/nar/gky764
21. Smogorzewska A. Fanconi Anemia: A Paradigm for Understanding DNA Repair During Replication. Blood (2019) 134:SCI–32-SCI-32. doi: 10.1182/blood-2019-121229
22. Gonzalez-Rajal A, Hastings JF, Watkins DN, Croucher DR, Burgess A. Breathing New Life into the Mechanisms of Platinum Resistance in Lung Adenocarcinoma. Front Cell Dev Biol (2020) 8:305:305. doi: 10.3389/fcell.2020.00305
23. O’Grady S, Finn SP, Cuffe S, Richard DJ, O'Byrne KJ, Barr MP. The role of DNA repair pathways in cisplatin resistant lung cancer. Cancer Treat Rev (2014) 40:1161–70. doi: 10.1016/j.ctrv.2014.10.003
24. Shukla A, Nguyen THM, Moka SB, Ellis JJ, Grady JP, Oey H, et al. Chromosome arm aneuploidies shape tumour evolution and drug response. Nat Commun (2020) 11:449. doi: 10.1038/s41467-020-14286-0
25. Stewart DJ. Mechanisms of resistance to cisplatin and carboplatin. Crit Rev Oncol Hematol (2007) 63:12–31. doi: 10.1016/j.critrevonc.2007.02.001
26. Ishida S, Lee J, Thiele DJ, Herskowitz I. Uptake of the anticancer drug cisplatin mediated by the copper transporter Ctr1 in yeast and mammals. Proc Natl Acad Sci USA (2002) 99:14298–302. doi: 10.1073/pnas.162491399
27. Kruh GD. Lustrous Insights into Cisplatin Accumulation: Copper Transporters. Clin Cancer Res (2003) 9:5807–9.
28. Lin X, Okuda T, Holzer A, Howell SB. The copper transporter CTR1 regulates cisplatin uptake in Saccharomyces cerevisiae. Mol Pharmacol (2002) 62:1154–9. doi: 10.1124/mol.62.5.1154
29. Safaei R, Holzer AK, Katano K, Samimi G, Howell SB. The role of copper transporters in the development of resistance to Pt drugs. J inorganic Biochem (2004) 98:1607–13. doi: 10.1016/j.jinorgbio.2004.05.006
30. Holzer AK, Samimi G, Katano K, Naerdemann W, Lin X, Safaei R, et al. The copper influx transporter human copper transport protein 1 regulates the uptake of cisplatin in human ovarian carcinoma cells. Mol Pharmacol (2004) 66:817–23. doi: 10.1124/mol.104.001198
31. Katano K, Kondo A, Safaei R, Holzer A, Samimi G, Mishima M, et al. Acquisition of resistance to cisplatin is accompanied by changes in the cellular pharmacology of copper. Cancer Res (2002) 62:6559–65.
32. Komatsu M, Sumizawa T, Mutoh M, Chen ZS, Terada K, Furukawa T, et al. Copper-transporting P-type adenosine triphosphatase (ATP7B) is associated with cisplatin resistance. Cancer Res (2000) 60:1312–6.
33. Cheng H, Zhang Z, Rodriguez-Barrueco R, Borczuk A, Liu H, Yu J, et al. Functional genomics screen identifies YAP1 as a key determinant to enhance treatment sensitivity in lung cancer cells. Oncotarget (2016) 7:28976. doi: 10.18632/oncotarget.6721
34. Jin L, Chun J, Pan C, Li D, Lin R, Alesi GN, et al. MAST1 Drives Cisplatin Resistance in Human Cancers by Rewiring cRaf-Independent MEK Activation. Cancer Cell (2018) 34:315–330.e317. doi: 10.1016/j.ccell.2018.06.012
35. Marini KD, Croucher DR, McCloy RA, Vaghjiani V, Gonzalez-Rajal A, Hastings JF, et al. Inhibition of activin signaling in lung adenocarcinoma increases the therapeutic index of platinum chemotherapy. Sci Transl Med (2018) 10(451):eaat3504. doi: 10.1126/scitranslmed.aat3504
36. Ding H, Zhao J, Zhang Y, Yu J, Liu M, Li X, et al. Systematic Analysis of Drug Vulnerabilities Conferred by Tumor Suppressor Loss. Cell Rep (2019) 27:3331–44.e3336. doi: 10.1016/j.celrep.2019.05.043
37. Ryan SL, Beard S, Barr MP, Umezawa K, Heavey S, Godwin P, et al. Targeting NF-kappaB-mediated inflammatory pathways in cisplatin-resistant NSCLC. Lung Cancer (2019) 135:217–27. doi: 10.1016/j.lungcan.2019.07.006
38. Wisniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods (2009) 6:359–62. doi: 10.1038/nmeth.1322
39. MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, et al. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinf (Oxford England) (2010) 26:966–8. doi: 10.1093/bioinformatics/btq054
40. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43:e47. doi: 10.1093/nar/gkv007
41. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol (2004) 3. doi: 10.2202/1544-6115.1027
42. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Society: Ser B (Methodological) (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
43. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinf (Oxford England) (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101
44. Abi Hussein H, Geneix C, Petitjean M, Borrel A, Flatters D, Camproux A-C. Global vision of druggability issues: applications and perspectives. Drug Discovery Today (2017) 22:404–15. doi: 10.1016/j.drudis.2016.11.021
45. López-Cortés A, Cabrera-Andrade A, Cruz-Segundo CM, Dorado J, Pazos A, Gonzáles-Díaz H, et al. Prediction of druggable proteins using machine learning and functional enrichment analysis: a focus on cancer-related proteins and RNA-binding proteins. bioRxiv (2019) 825513. doi: 10.1101/825513
46. Coker EA, Mitsopoulos C, Tym JE, Komianou A, Kannas C, Di Micco P, et al. canSAR: update to the cancer translational research and drug discovery knowledgebase. Nucleic Acids Res (2019) 47:D917–22. doi: 10.1093/nar/gky1129
47. Xia J, Sun J, Jia P, Zhao Z. Do cancer proteins really interact strongly in the human protein-protein interaction network? Comput Biol Chem (2011) 35:121–5. doi: 10.1016/j.compbiolchem.2011.04.005
48. Liu T, Lin Y, Wen X, Jorissen RN, Gilson MK. BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities. Nucleic Acids Res (2007) 35:D198–201. doi: 10.1093/nar/gkl999
49. Tang H, Xiao G, Behrens C, Schiller J, Allen J, Chow CW, et al. A 12-gene set predicts survival benefits from adjuvant chemotherapy in non-small cell lung cancer patients. Clin Cancer Res (2013) 19:1577–86. doi: 10.1158/1078-0432.CCR-12-2321
50. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res (2019) 47:D607–d613. doi: 10.1093/nar/gky1131
51. Becker JP, Weiss J, Theile D. Cisplatin, oxaliplatin, and carboplatin unequally inhibit in vitro mRNA translation. Toxicol Lett (2014) 225:43–7. doi: 10.1016/j.toxlet.2013.11.015
52. Kaur G, Gupta SK, Singh P, Ali V, Kumar V, Verma M. Drug-metabolizing enzymes: role in drug resistance in cancer. Clin Trans Oncol (2020) 22:1667–80. doi: 10.1007/s12094-020-02325-7
53. Li Y, Steppi A, Zhou Y, Mao F, Miller PC, He MM, et al. Tumoral expression of drug and xenobiotic metabolizing enzymes in breast cancer patients of different ethnicities with implications to personalized medicine. Sci Rep (2017) 7:4747. doi: 10.1038/s41598-017-04250-2
54. Wilson SM, Khanna R. Specific binding of lacosamide to collapsin response mediator protein 2 (CRMP2) and direct impairment of its canonical function: implications for the therapeutic potential of lacosamide. Mol Neurobiol (2015) 51:599–609. doi: 10.1007/s12035-014-8775-9
55. Ivanov AA, Khuri FR, Fu H. Targeting protein-protein interactions as an anticancer strategy. Trends Pharmacol Sci (2013) 34:393–400. doi: 10.1016/j.tips.2013.04.007
56. Bruno PM, Liu Y, Park GY, Murai J, Koch CE, Eisen TJ, et al. A subset of platinum-containing chemotherapeutic agents kills cells by inducing ribosome biogenesis stress. Nat Med (2017) 23:461–71. doi: 10.1038/nm.4291
57. Zhang C, Lin X, Zhao Q, Wang Y, Jiang F, Ji C, et al. YARS as an oncogenic protein that promotes gastric cancer progression through activating PI3K-Akt signaling. J Cancer Res Clin Oncol (2020) 146:329–42. doi: 10.1007/s00432-019-03115-7
58. Gou Y, Shi Y, Zhang Y, Nie Y, Wang J, Song J, et al. Ribosomal protein L6 promotes growth and cell cycle progression through upregulating cyclin E in gastric cancer cells. Biochem Biophys Res Communicationss (2010) 393:788–93. doi: 10.1016/j.bbrc.2010.02.083
59. Terzian T, Dumble M, Arbab F, Thaller C, Donehower LA, Lozano G, et al. Rpl27a mutation in the sooty foot ataxia mouse phenocopies high p53 mouse models. J Pathol (2011) 224:540–52. doi: 10.1002/path.2891
60. Yu C, Hong H, Zhang S, Zong Y, Ma J, Lu A, et al. Identification of key genes and pathways involved in microsatellite instability in colorectal cancer. Mol Med Rep (2019) 19:2065–76. doi: 10.3892/mmr.2019.9849
61. Liu K, Kang M, Li J, Qin W, Wang R. Prognostic value of the mRNA expression of members of the HSP90 family in non−small cell lung cancer. Exp Ther Med (2019) 17:2657–65. doi: 10.3892/etm.2019.7228
62. Fukunaga Y, Bandoh S, Fujita J, Yang Y, Ueda Y, Hojo S, et al. Expression of cytokeratin 8 in lung cancer cell lines and measurement of serum cytokeratin 8 in lung cancer patients. Lung Cancer (2002) 38:31–8. doi: 10.1016/S0169-5002(02)00109-5
63. Choi Y-Y, Lee SY, Lee WK, Jeon H-S, Lee EB, Lee HC, et al. RACK1 is a candidate gene associated with the prognosis of patients with early stage non-small cell lung cancer. Oncotarget (2015) 6:4451. doi: 10.18632/oncotarget.2865
64. Zheng W, Jiang C, Li R. Integrin and gene network analysis reveals that ITGA5 and ITGB1 are prognostic in non-small-cell lung cancer. Onco Targets Ther (2016) 9:2317–27. doi: 10.2147/OTT.S91796
65. Fang J, Wang H, Liu Y, Ding F, Ni Y, Shao S. High KRT8 expression promotes tumor progression and metastasis of gastric cancer. Cancer Sci (2017) 108:178–86. doi: 10.1111/cas.13120
66. Liu Q-Z, Gao X-H, Chang W-J, Gong H-F, Fu C-G, Zhang W, et al. Expression of ITGB1 predicts prognosis in colorectal cancer: a large prospective study based on tissue microarray. Int J Clin Exp Pathol (2015) 8:12802–10.
67. Cao X-X, Xu J-D, Liu X-L, Xu J-W, Wang W-J, Li Q-Q, et al. RACK1: A superior independent predictor for poor clinical outcome in breast cancer. Int J Cancer (2010) 127:1172–9. doi: 10.1002/ijc.25120
68. Moreb JS, Muhoczy D, Ostmark B, Zucali JR. RNAi-mediated knockdown of aldehyde dehydrogenase class-1A1 and class-3A1 is specific and reveals that each contributes equally to the resistance against 4-hydroperoxycyclophosphamide. Cancer Chemother Pharmacol (2007) 59:127–36. doi: 10.1007/s00280-006-0233-6
69. Kim J, Shin JH, Chen C-H, Cruz L, Farnebo L, Yang J, et al. Targeting aldehyde dehydrogenase activity in head and neck squamous cell carcinoma with a novel small molecule inhibitor. Oncotarget (2017) 8:52345–56. doi: 10.18632/oncotarget.17017
70. Zhou L, Sheng D, Wang D, Ma W, Deng Q, Deng L, et al. Identification of cancer-type specific expression patterns for active aldehyde dehydrogenase (ALDH) isoforms in ALDEFLUOR assay. Cell Biol Toxicol (2019) 35:161–77. doi: 10.1007/s10565-018-9444-y
71. Sawers L, Ferguson M, Ihrig B, Young H, Chakravarty P, Wolf C, et al. Glutathione S-transferase P1 (GSTP1) directly influences platinum drug chemosensitivity in ovarian tumour cell lines. Br J Cancer (2014) 111:1150–8. doi: 10.1038/bjc.2014.386
72. De Luca A, Parker LJ, Ang WH, Rodolfo C, Gabbarini V, Hancock NC, et al. A structure-based mechanism of cisplatin resistance mediated by glutathione transferase P1-1. Proc Natl Acad Sci (2019) 116:13943–51. doi: 10.1073/pnas.1903297116
73. Dinavahi SS, Bazewicz CG, Gowda R, Robertson GP. Aldehyde dehydrogenase inhibitors for cancer therapeutics. Trends Pharmacol Sci (2019) 40:774–89. doi: 10.1016/j.tips.2019.08.002
74. Lewis SJ, Cherry NM, Niven RM, Barber PV, Povey AC. GSTM1, GSTT1 and GSTP1 polymorphisms and lung cancer risk. Cancer Lett (2002) 180:165–71. doi: 10.1016/s0304-3835(02)00028-9
75. Gao J, Schatton D, Martinelli P, Hansen H, Pla-Martin D, Barth E, et al. CLUH regulates mitochondrial biogenesis by binding mRNAs of nuclear-encoded mitochondrial proteins. J Cell Biol (2014) 207:213–23. doi: 10.1083/jcb.201403129
76. Schatton D, Pla-Martin D, Marx M-C, Hansen H, Mourier A, Nemazanyy I, et al. CLUH regulates mitochondrial metabolism by controlling translation and decay of target mRNAs. J Cell Biol (2017) 216:675–93. doi: 10.1083/jcb.201607019
77. Wakim J, Goudenege D, Perrot R, Gueguen N, Desquiret-Dumas V, Chao de la Barca JM, et al. CLUH couples mitochondrial distribution to the energetic and metabolic status. J Cell Sci (2017) 130:1940–51. doi: 10.1242/jcs.201616
78. Lee KM, Giltnane JM, Balko JM, Schwarz LJ, Guerrero-Zotano AL, Hutchinson KE, et al. MYC and MCL1 Cooperatively Promote Chemotherapy-Resistant Breast Cancer Stem Cells via Regulation of Mitochondrial Oxidative Phosphorylation. Cell Metab (2017) 26:633–47.e637. doi: 10.1016/j.cmet.2017.09.009
79. Farge T, Saland E, de Toni F, Aroua N, Hosseini M, Perry R, et al. Chemotherapy-Resistant Human Acute Myeloid Leukemia Cells Are Not Enriched for Leukemic Stem Cells but Require Oxidative Metabolism. Cancer Discovery (2017) 7:716–35. doi: 10.1158/2159-8290.Cd-16-0441
80. Kalainayakan SP, FitzGerald KE, Konduri PC, Vidal C, Zhang L. Essential roles of mitochondrial and heme function in lung cancer bioenergetics and tumorigenesis. Cell Biosci (2018) 8:56–6. doi: 10.1186/s13578-018-0257-8
81. Xu Z-H, Freimuth RR, Eckloff B, Wieben E, Weinshilboum RM. Human 3′-phosphoadenosine 5′-phosphosulfate synthetase 2 (PAPSS2) pharmacogenetics: gene resequencing, genetic polymorphisms and functional characterization of variant allozymes. Pharmacogenetics Genomics (2002) 12:11–21. doi: 10.1097/00008571-200201000-00003
82. Xu Z-H, Otterness DM, Freimuth RR, Carlini EJ, Wood TC, Mitchell S, et al. Human 3′-Phosphoadenosine 5′-Phosphosulfate Synthetase 1 (PAPSS1) and PAPSS2: Gene Cloning, Characterization and Chromosomal Localization. Biochem Biophys Res Commun (2000) 268:437–44. doi: 10.1006/bbrc.2000.2123
Keywords: non-small cell lung cancer, platinum-based chemotherapy, cisplatin, quantitative proteomics, biomarkers, therapeutic targets, machine learning
Citation: Ryan S-L, Dave KA, Beard S, Gyimesi M, McTaggart M, Sahin KB, Molloy C, Gandhi NS, Boittier E, O’Leary CG, Shah ET, Bolderson E, Baird A-M, Richard DJ, O’Byrne KJ and Adams MN (2021) Identification of Proteins Deregulated by Platinum-Based Chemotherapy as Novel Biomarkers and Therapeutic Targets in Non-Small Cell Lung Cancer. Front. Oncol. 11:615967. doi: 10.3389/fonc.2021.615967
Received: 10 October 2020; Accepted: 19 January 2021;
Published: 11 March 2021.
Edited by:
Jessica Desiree Menis, Veneto Institute of Oncology (IRCCS), ItalyReviewed by:
Pasquale Pisapia, University of Naples Federico II, ItalyMatthias Scheffler, University Hospital of Cologne, Germany
Copyright © 2021 Ryan, Dave, Beard, Gyimesi, McTaggart, Sahin, Molloy, Gandhi, Boittier, O’Leary, Shah, Bolderson, Baird, Richard, O’Byrne and Adams. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Derek J. Richard, derek.richard@qut.edu.au; Kenneth J. O’Byrne, k.obyrne@qut.edu.au; Mark N. Adams, mn.adams@qut.edu.au