- 1Department of Cancer Immunology, Institute for Cancer Research, Oslo University Hospital, Oslo, Norway
- 2K.G. Jebsen Centre for B Cell Malignancies, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
- 3Oslo Centre for Biostatistics and Epidemiology, University of Oslo, Oslo, Norway
- 4Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology, Trondheim, Norway
- 5Department of Hematology, St. Olav’s University Hospital, Trondheim, Norway
- 6Department of Immunology and Transfusion Medicine, St. Olav’s University Hospital, Trondheim, Norway
- 7Norwegian Cancer Genomics Consortium, Oslo University Hospital, Oslo, Norway
- 8Department of Tumor Biology, Institute for Cancer Research, Oslo University Hospital, Oslo, Norway
- 9Centre for Cancer Cell Reprogramming, Institute of Clinical Medicine, Faculty of Medicine, University of Oslo, Oslo, Norway
- 10Genomics Core Facility, Department of Core Facilities, Institute for Cancer Research, Oslo University Hospital, Oslo, Norway
- 11Center for Bioinformatics, Department of Informatics, University of Oslo, Oslo, Norway
- 12Bioinformatics Core Facility, Norwegian University of Science and Technology, Trondheim, Norway
Introduction: Multiple myeloma (MM) is a heterogeneous disease where cancer-driver mutations and aberrant signaling may lead to disease progression and drug resistance. Drug responses vary greatly, and there is an unmet need for biomarkers that can guide precision cancer medicine in this disease.
Methods: To identify potential predictors of drug sensitivity, we applied integrated data from drug sensitivity screening, mutational analysis and functional signaling pathway profiling in 9 cell line models of MM. We studied the sensitivity to 33 targeted drugs and their association with the mutational status of cancer-driver genes and activity level of signaling proteins.
Results: We found that sensitivity to mitogen-activated protein kinase kinase 1 (MEK1) and phosphatidylinositol-3 kinase (PI3K) inhibitors correlated with mutations in NRAS/KRAS, and PI3K family genes, respectively. Phosphorylation status of MEK1 and protein kinase B (AKT) correlated with sensitivity to MEK and PI3K inhibition, respectively. In addition, we found that enhanced phosphorylation of proteins, including Tank-binding kinase 1 (TBK1), as well as high expression of B cell lymphoma 2 (Bcl-2), correlated with low sensitivity to MEK inhibitors.
Discussion: Taken together, this study shows that mutational status and signaling protein profiling might be used in further studies to predict drug sensitivities and identify resistance markers in MM.
Introduction
MM is considered a treatable, but generally incurable disease, with a heterogeneous clinical course as one of its hallmarks. Several types of molecular alterations, such as single-nucleotide variants and cytogenetic abnormalities, are responsible for disease initiation, maintenance, and progression (1). Advances in molecular profiling technologies have enabled us to characterize the molecular landscape of aberrations, which in turn improves our understanding of the underlying cancer biology, and indicates potential molecular treatment targets.
Many MM driver mutations occur in known cancer-signaling pathways, such as the Rat sarcoma virus (RAS), mitogen-activated protein kinase (MAPK), and PI3K/AKT pathways, which are known to control proliferation and survival of MM cells (2, 3). Sequencing analysis of MM patient samples has shown that the Kirsten RAS (KRAS) oncogene is the most commonly mutated gene (36%) in the disease, followed by the neuroblastoma RAS (NRAS) (20%), with frequent co-existence of one or more variants in both KRAS and NRAS cancer-driver genes (4–6). Interestingly, refractory MM patients with multi-drug resistance to standard myeloma therapies (e.g. proteasome inhibitors and immunomodulatory drugs) showed a significant increase (72%) in the mutation rate of RAS pathway genes, as compared to newly diagnosed MM cases (7). This suggests that RAS mutations play a role in acquired drug resistance in refractory MM. Targeting the RAS signaling pathway is therefore of potentially high therapeutic interest.
Studies in MM have shown that molecular features, including mutations, translocations and copy-number abnormalities, have both prognostic and predictive value and may enable further improvement in patient outcomes if employed to define personalized treatment strategies (8). Inhibition of RAS effectors, such as MEK, has emerged as a viable strategy for the treatment of KRAS/NRAS mutant MM clones (9, 10). Specifically, treatment of RAS/RAF mutant MM with the MEK inhibitor trametinib has shown good tolerability and durable remission in some patients (9, 11).
Besides the RAS signaling pathway, the PI3K/AKT/mTOR pathway is activated in a significant proportion of MM patients (3). Increased levels of MM-promoting cytokines, including interleukin-6 (IL6) and insulin-like growth factor-1 (IGF1), have been reported to be involved in the activation of this pathway as well as the aberrant up-regulation of other pathways that feed into the PI3K/AKT activation complex (12). Inhibition of PI3K/AKT/mTOR pathway induces apoptosis in MM (13), and targeting this pathway may therefore show therapeutic benefit.
MM mutational status is currently used to guide the use of targeted agents in precision medicine trials, including the MyDrug study (NCT03732703), the CAPTUR study (NCT03297606), and the MATCH screening trial (NCT02465060). Protein expression and phosphorylation levels also impact drug responsiveness, arguing for integrating both genetic and functional analyses in treatment selection strategies (14–16).
Here, we used high-throughput DNA sequencing to characterize the genetic make-up of nine MM cell lines. Based on the spectrum of mutated targets and their druggability, we designed a library of 33 targeted drugs, which was used on the same MM cell lines. Associations between mutational status and drug sensitivity were investigated. Next, we used a dataset of signaling readouts (n=31) available on the same cell lines (16) to study how drug sensitivity correlated with baseline protein expression and phosphorylation levels of cancer-driving signaling proteins.
Our findings showed that mutant NRAS/KRAS MM cells were highly sensitive to MEK inhibition, as compared to wild-type (WT) MM cells. Furthermore, a systematic correlation analysis of drug sensitivities and signaling protein readouts revealed several drug dependencies. Overall, both mutational status and protein phosphorylation/expression status may help elucidate drug-specific sensitivities in a small panel of MM cell lines, which could next guide future studies to predict biomarkers of drug sensitivity in the context of MM.
Materials and methods
Cell lines
Myeloma cell lines used in this study were U266 (17), JJN3 (18), CAG (19), INA6 (20), OH2 (21, 22), IH1 (23), KJON, VOLIN (24), URVIN, and FOLE. The FOLE MM cell line (Misund et al., unpublished), has biallelically-inactive TP53, and was used as a positive control for testing of the MDM2-P53 interaction inhibitor nutlin-3A. The MM cell lines were cultured in RPMI 1640 medium (ThermoFisher Scientific, Waltham, MA, USA) supplemented with 2mM L-glutamine, 1% Penicillin-Streptomycin, 1x Sodium Pyruvate (NaPur), and fetal bovine serum (FBS) at 10% (JJN3, CAG, INA6, VOLIN, URVIN) or 15% FBS (U266), or 10% human serum (Sigma-Aldrich, Saint-Louis, MO, USA) (IH1, OH2, KJON). The culture medium was supplemented with 2 ng/mL IL6 (ThermoFisher Scientific) for culturing of INA6, KJON, IH1, OH2, URVIN, and VOLIN. The cells were expanded, aliquoted and cryopreserved until experimental assays were performed. See Table 1 for MM cell line characteristics.
Cell viability assay
Drugs (n=33) were selected based on mutated targets in the cell lines (see Tables 2, 3 for drugs and their corresponding mutated targets), and were added to 384-well TC-microplates (Greiner #781080) using an acoustic dispenser (Echo 550, Labcyte Inc., CA, USA). Each drug was tested at five concentrations ranging from 1 to 10,000 nM. Experiments on cell lines were done on freshly thawed cells. MM cells (5000 cells per well in 25 ul volume) were transferred into plates using an automatic dispenser (Certus Flex, Fitz Gyger, Thun, Switzerland), and incubated at 37°C for 72h. Cell viability was assessed by the CellTiter-Glo luminescence assay (Promega, WI, USA) according to the manufacturer’s instructions. Luminescence was recorded with an Envision Xcite plate reader (Perkin Elmer, MA, USA). The raw concentration-response data were analyzed using the KNIME software (AG, Zurich, Switzerland) and Rstudio Team (Boston, MA) (25). Normalization of the response readout was done to the negative (0.1% DMSO) and positive (100 µM benzethonium chloride) controls.
Phospho flow
Phospho flow assays were performed on freshly thawed myeloma cells as described previously (26). Antibody-stained samples were run on a BD LSR Fortessa and output data were analyzed in Cytobank (https://cellmass.cytobank.org/cytobank/). Raw data were transformed to an arcsinh ratio relative to the signal of an isotype control, which was set to zero.
Data analysis and statistics
Curve-fits of normalized concentration-response data used the function drm from the R package drc (https://www.r-project.org/) with the four parameter log-logistic model, LL.4, or the logistic model, L.4, where LL.4 failed to converge. To quantify drug responses, a modified drug sensitivity score (DSS) was calculated for each drug (27). In this modified function, area under the curve was calculated using a response-window from 100% to 10%, and a concentration-window from the minimum concentration tested to the concentration where the viability reached 10% (threshold %). DSS type 1 was used, without the term for division by the logarithm of the upper limit. The DSS scores were calculated on a scale of 0-100. A high DSS therefore indicates that MM cells are drug-sensitive, while a low DSS indicates drug-resistant MM cells.
Statistical analyses of output data from drug sensitivity screens and phospho flow assays were performed in GraphPad Prism v8 (San Diego, CA, USA). Welch’s t test [(*p < 0.05, ***p < 0.001, ns (not significant)] was used to compare two means, as indicated in the respective figure legends. ClustVis (28) was used for unsupervised clustering of the DSS values and column annotations indicating genes with the presence of point mutations. No scaling is applied to rows. Columns were clustered using Euclidean distance and Ward linkage. The significance of the correlation coefficients between DSS values and phospho flow readouts was assessed by Pearson’s r test for Figures 1B, C and Supplementary Figure 2, where the p-values prior to multiple-testing correction are shown.
Predictive models for drug response (as defined by the DSS) were constructed using either mutation or phospho flow data as covariates. In each case, a type of multi-response penalized linear regression model known as tree-lasso (29, 30) was used for feature selection and to assess the predictive power of the selected features for drug response. Regression coefficients for all 33 drugs were penalized jointly in a hierarchical fashion, with penalties weighted by the strengths of the correlations among drug responses, as defined by the height parameters in the hierarchical tree structure of the DSS data (Supplementary Figure 1). For optimal results, the tree height parameters were cut-off at 0.7, on a scale where the root of the tree is at a height of 1 and the leaves are at a height of zero.
Both the mutation-dependent and protein-dependent models were trained using leave-one-out cross validation to compute the Mean Squared Error (MSE) as a function of the tuning parameter. The tree-lasso cost function was optimized using the sub-function tree.lasso from the R package IPFStructPenalty, available at https://github.com/zhizuio/IPFStructPenalty. In each model, the tuning parameter Λ was set to a value at which the cross-validated MSE becomes essentially flat and the smoothing parameter for Smoothing Proximal Gradient descent (SPG) optimization was set to μ = 10-4. For the mutation-dependent model, the tuning parameter was set to Λ = 25 (with a cross-validated MSE of 175), and for the protein-dependent model it was set to Λ = 24 (with a cross-validated MSE of 530). The tree-lasso regression coefficients of each model were then computed using the optimal values of the tuning parameter Λ (Figures 2A, 1A, B).
The amount of variance explained by each model was also assessed by the R squared value, which was computed for each tuned model as R2 = 1 - RSS/TSS, where the Residual Sum of Squares (RSS) and Total Sum of Squares (TSS) were both defined using all 9 cell lines at once. The R squared value was R2 = 0.39 for the mutation-dependent model, and R2 = 0.58 for the protein-dependent model.
The covariates for the protein-dependent model were defined as the standardized phospho flow data, i.e., for each protein the phospho flow arcsinh ratios were centralized at their mean value across samples and normalized by their standard deviation prior to the tree-lasso model calculations. As a consistency check, an overall matching between the signs of the regression coefficients and those of the protein-drug Pearson correlation coefficients was observed (Supplementary Figure 1B), since the standardization of the phospho flow data allows analogous interpretations for these two types of coefficients.
Targeted high-throughput sequencing analysis – variant calling
Targeted DNA sequencing was performed using the SureSelect Human Kinome kit (Agilent Technologies), with capture probes targeting 3.2 Mb of the human genome, including exons and untranslated regions (UTRs) of all known kinases and selected cancer‐related genes (n = 612). Paired-end sequencing reads of 100-bp length were aligned to the human reference genome (hg19) with Novoalign (version 2.08.3), followed by filtering and realignment with GATK tools and Picard. Single point mutations were identified using MuTect (version 1.1.4). Variant consequence annotation was performed with ANNOVAR, using RefSeq as the underlying transcript model. Details of the complete variant calling pipeline that was applied on the MM cell lines have been described previously (31). The mean sequencing coverage of the kinome was 477x, and we obtained a minimum coverage of 150x for 87.6% (VOLIN) and 90.0% (KJON-1) of the kinome targeted regions. Variant calling performed on cancer cell lines without utilizing matched normal samples is bound to generate a mix of germline variants and somatic variants. We therefore set up a set of filtering procedures to both i) exclude known germline variants, and ii) enrich for coding, cancer-associated variants. Specifically, we excluded all variants that overlapped with germline variants found in the 1000 Genomes Project phase 3 (minor allele frequency > 1% in any population), and NHLBI Exome Sequencing Project (minor allele frequency > 0.1% in any population) (32, 33). In addition, we excluded variants present in the single nucleotide polymorphism database (dbSNP) (build 138) (34) that had no clinical associations (as given from ClinVar cross-references). Finally, we restricted the variant set to coding variants (missense, stop-gain/stop-loss, frameshift/non-frameshift, splice site donor/acceptor) in known cancer census genes (COSMIC version 68) (35). All variants were subjected to a functional annotation workflow that included UniProt KB (functional protein properties) and Polymorphism Phenotyping v2 (PolyPhen-2) web server (computational predictions of effect of coding variants) (http://genetics.bwh.harvard.edu/pph2/) (36).
Results
Identification of mutations in MM cell lines
In order to identify functionally relevant and actionable mutations in MM, nine MM cell lines (Table 1) were subjected to targeted high-throughput DNA sequencing. Short sequencing reads were processed with a variant calling pipeline and subsequent variant filtering procedure, from which we detected a total of 136 mutated genes in the cell lines (Supplementary Table 1). Mutated genes included NRAS/KRAS, BRAF, RAF-1, TP53, PIK3CA, PIK3R3, MTOR, FLT1, FLT3, FLT4, EGFR, and SYK, which are known to be frequently mutated in MM and other related blood malignancies (37–39), as well as to have a therapeutic potential in MM (39) (Tables 2, 3).
To prioritize the identified gene variants for functional impact, all protein-coding alterations were subjected to analysis with PolyPhen2 (see Materials and Methods), a tool that provides computational predictions for the functional impact of amino acid changes. We identified 31 amino acid changes with predicted damaging effects which were druggable (39), including the most oncogenic RAS mutated isoform, G12 and Q61, as well as PIK3CA catalytic subunit mutations, p.H1047R and p.W590C (40, 41) (Tables 2, 3 and Supplementary Table 1).
Targeting the mutational status in MM cell lines
Based on the mutation analysis, we designed a drug library consisting of 33 clinically relevant drugs targeting identified druggable gene products (Table 3). Drug sensitivity screens were then performed on the 9 MM cell lines.
An unsupervised clustering of the cell lines by their DSS profiles showed that MM cell lines with PI3K and NRAS/KRAS mutations displayed high sensitivity to MEK and PI3K inhibitors, respectively (Figure 1A). Concentration-response effects on the MM cell viability to individual inhibitors with and without associated targets are shown in Figures 1B, C. These findings support our hypothesis that the presence of RAS and PI3K mutations may increase MM cell sensitivity to MEK and PI3K inhibitors.
Figure 1 Targeting point mutations in cancer-driving genes in MM cell line models. A–C, Freshly thawed cells from the indicated MM cell lines were dispensed in 384-well plates pre-coated with a customized drug library of 33 single drugs. At 72h, cell viability was assessed by CellTiter-Glo. A drug sensitivity score (DSS) was calculated ranging from 0 to 100 for the entire drug library (see Materials and Methods). High DSS indicates high drug sensitivity. (A), MM cell lines (n=9) were exposed to single drugs as described above. DSS (rows) was calculated for each drug and cell line and plotted as a heatmap. Columns are clustered using Euclidean distance and the Ward linkage method. Column annotations (top) indicate mutated genes, where blue means mutation present, while light blue absent. (B, C), Concentration-response curves for the effects of taselisib and trametinib, on the viability of indicated MM cell lines (72h) treated as above. Graphs show mean viability ± standard deviation (SD, n=3). MM cell lines with mutations in PI3K/mTOR and RAS target genes are indicated in red, whereas the MM cell lines in different colored lines are wild-type and each line indicates cell viability effects in each individual cell line.
We observed variability in efficacy among drugs targeting the PI3K signaling pathway. The most effective drugs were temsirolimus (mean DSS ± SD = 51.96 ± 13.66), taselisib (mean DSS ± SD = 41.82 ± 17.10), pictilisib (mean DSS ± SD = 32.50 ± 13.66), and copanlisib (mean DSS = 19.83 ± 18.65), while treatment with idelalisib and GSK2636771 resulted in markedly low DSS (mean DSS = 2.29 ± 2.5, mean DSS = 0.26 ± 0.46, respectively) (Figure 2A). Among the MEK inhibitors, trametinib induced the highest response (mean DSS ± SD = 47.20 ± 27.46), followed by binimetinib (mean DSS ± SD = 21.88 ± 21.54), cobimetinib (mean DSS ± SD = 21.11 ± 23.80), and U0126 (mean DSS ± SD = 11.83 ± 13.85) (Figure 2A).
Figure 2 Drug sensitivity profile in MM cell lines. (A), Drug sensitivity to the indicated drugs across the 9 MM cell lines shown in Figure 1. The graph shows DSS for the indicated drugs sorted from most to least effective as mean ± SD. Green and pink bars indicate PI3K/mTOR inhibitors (PI3K/mTORi) and MEK inhibitors (MEKi), respectively. (B), DSS distribution for the three most effective drugs (target indicated): temsirolimus (mTOR, PI3K, AKT), trametinib (MEK1/2), and taselisib (PI3K) inhibitors in the MM cell lines. MM cell lines with targeted mutations are indicated in different colors.
Interestingly, the MM cell lines CAG and URVIN, harboring PIK3CA and PIK3R3 mutations, respectively, exhibited the highest sensitivity to taselisib and temsirolimus treatment, while JJN3, harboring mutations in PIK3CA and mTOR genes, showed less sensitivity to taselisib compared to CAG and URVIN cell lines. On the other hand, JJN3 gained sensitivity to temsirolimus that also target mTOR (Figure 2B). The MM cell lines IH1, OH2, and KJON, with NRAS/KRAS mutations, were found highly sensitive to trametinib (Figure 2B). These results indicate that sensitivity to kinase inhibitors may be associated with corresponding pathway mutations.
Mutational status as a predictor for in vitro drug responsiveness
To study whether the mutational status could predict drug sensitivity, we compared treatment responses to PI3K and MEK inhibitors between cell lines with and without selected mutations. We found that cell lines with PI3K pathway mutations (i.e. PIK3CA, PIK3R3, and mTOR) or RAS gene mutations (i.e. p.G12V, p.G12D, p.G13D, p.Q61K, and p.Q61H) (Table 2) exhibited significantly higher DSS relative to cell lines with WT forms of the listed genes when comparing aggregated data for each class of drugs tested (Figures 3A, B, left panels). The three PI3K inhibitors, taselisib, copanlisib, and pictilisib, each showed a trend towards association with PI3K mutational status, but this was not statistically significant (Figure 3A, right panels). Of the four MEK inhibitors tested, a statistically significant association was found between the response to binimetinib and RAS mutational status whereas trends towards statistical significance were observed for trametinib, cobimetinib, and U0126 (Figure 3B, right panels).
Figure 3 Associations between drug sensitivity and mutational profiles in MM cell lines. (A, B), Pharmacogenomic comparison between PI3K and RAS mutation status and DSS to PI3K and MEK inhibitors. Aggregated (left) and individual drug (middle-left) effects are shown. The graph shows mean DSS ± SD. (C), Heatmap of the tree-lasso regression coefficients for DSS as a function of mutations, where the rows are all the possible selected features (i.e., genes whose regression coefficients are exactly zero for all drugs were excluded from the plot, as well as those coefficients that, for every drug, are lower in absolute value than 5% of the sum of all coefficients for the given drug). Framed rectangles highlight the predictive power of NRAS and PIK3CA for the responses to MEKi and PI3Ki, respectively. (D), Concentration-response curves of cell viability for nutlin-3A in the IH1 and FOLE MM cell lines that differ in TP53 mutational status. The graph shows mean viability ± SD (n=3). (E), Association of the TP53 mutation status with nutlin-3A sensitivity for all 10 MM cell lines. The graph shows mean DSS ± SD. Statistics were performed using Welch’s test. Statistics were performed using Welch’s test. *p < 0.05, **p < 0.01, ***p < 0.001, ns (not significant).
Interestingly, by applying a tree-lasso regression model as described in Materials and Methods, a few mutations were found to be associated in varying degrees with high sensitivity to several drugs, as indicated by the positive values in the sparse matrix of regression coefficients (Figure 3C). In particular, the mutational status of the PIK3CA and NRAS genes were selected as predictors for drug sensitivity to PI3K and MEK inhibitors (Figure 3C).
We also noticed that, while the prediction analyses demonstrated that the mutational status may be useful to predict drug sensitivity, the predictive impact of NRAS mutations on the drug sensitivity to MEK inhibitors was higher than the one observed for PIK3CA mutations on PI3K inhibitors, as indicated by the color-intensity of the heatmap in Figure 3C.
The inhibition of the MDM2-TP53 interaction by nutlin-3A is an attractive strategy to stabilize the P53 mediated apoptosis in various WT TP53 tumors, including MM (42–44), and therefore worth investigating. Interestingly, we observed differing sensitivities to the MDM2-TP53 inhibitor nutlin-3A between the TP53 WT IH1 and the biallelically TP53 mutant MM cell line FOLE. (Figure 3D). When stratifying nutlin-3A responses on TP53 mutational status in all MM cell lines tested, we observed a significantly higher drug sensitivity in the TP53 WT cell lines, in agreement with earlier reports (14, 15, 45) (Figure 3E).
Protein expression level affects drug sensitivities
Having demonstrated that gene mutations in MM cell lines are linked to drug sensitivities towards PI3K and MEK inhibitors, we were then interested in studying how expression levels of the pathway (phospho) proteins correlated with their drug sensitivity.
In order to identify relevant signaling proteins whose observed basal expression or phosphorylation levels across different cell lines can explain corresponding DSS values for different drugs, our first approach was a correlation analysis using Pearson correlation coefficients. Correlation coefficients for all possible protein-drug pairs, including the 33 drugs in our library, and the 31 relevant signaling proteins (16) selected for our study, were computed using the phospho flow readouts and DSS values from all the MM cell lines (Supplementary Figure 1B). Multiple significance tests of the Pearson correlation coefficients followed by a false discovery rate controlling procedure revealed that this method is not sufficiently robust to capture significant correlations. Our solution was to model all the drugs jointly using the tree-lasso regression approach described in Materials and Methods. By penalizing regression coefficients jointly according to the hierarchical clustering tree of the correlations among DSS values of different drugs, we were able to achieve enough sensitivity to capture protein-drug correlations that had been lost in our previous, naïve approach (Figure 4A). We found that PI3K and MEK inhibitors emerged as separate clusters, each with a common set of predictive variables, as indicated by a few non-zero tree-lasso regression coefficients (Figure 4A). We note that since the covariates in our model were defined as the standardized phospho flow readouts, the tree-lasso regression coefficients have a similar interpretation to correlation coefficients (i.e., positive or negative values for the coefficients have similar meanings in both cases), and therefore can be regarded as a special type of correlation.
Figure 4 Association between drug sensitivity and expression or activation profile of selected intracellular cancer-driven proteins. A–C, Freshly thawed MM cell lines (n=9) were fixed, permeabilized and stained with antibodies against the indicated proteins or phosphoprotein epitopes (rows). Signals were analyzed by flow cytometry (see Materials and Methods). (A), Heatmap of the tree-lasso regression coefficients for DSS as a linear function of the standardized phospho flow data, with proteins as rows and drugs as columns. Positive and negative regression coefficients are represented in shades of red and blue, respectively. B,C, Pearson correlation plots for DSS to the indicated MEK inhibitors trametinib (B) and binimetinib (C) versus three signaling readouts inversely associated with drug response.
Notably, using the tree-lasso regression model we found that high levels of phospho-TBK1 inversely correlated with low drug sensitivity to MEK inhibitors, including trametinib, binimetinib, and cobimetinib (Figure 4A). Although these correlations seem to be in agreement with individual Pearson correlation plots of the TBK1 levels against DSS for trametinib and binimetinib (upper plots in Figures 4B, C), we emphasize that their true significance is only revealed by an integrated model that explores the data structure of all the drugs to define a new, more robust, type of correlation.
As also observed in our previous study (16), we found that high levels of Bcl-2 and phospho-Bcl-2 predicted low sensitivity to trametinib, binimetinib, and cobimetinib (Figures 4A-C middle and lower plots). Once again, individual Pearson correlation plots (Figures 4B, C middle and lower plots) are not enough to identify these relationships, since their statistical significance were only seen after using the more robust tree-lasso regression approach for variable selection. Also, worth mentioning is our finding that phospho-MEK1 (pS298) is positively correlated (in the sense of positive tree-lasso regression coefficients) with sensitivity to the same MEK inhibitors (trametinib, binimetinib, and cobimetinib). Moreover, we observed that high phosphorylation levels of AKT (pS473) correlated with high responses to PI3K inhibitors taselisib, pictilisib and copanlisib (Figure 4A). In this case, Pearson correlation p-values are too high even before multiple-testing correction (Supplementary Figures 2A-D), however, tree-lasso regression was still sensitive enough to capture a relatively weak correlation.
Taken together, our results indicate that the expression and phosphorylation levels of signaling proteins can inform on drug sensitivity, and therefore, can provide relevant information as part of functional drug testing studies.
Discussion
Currently, the choice of therapy for MM patients is mostly based on physician’s decision and patient’s clinical status, including age and comorbidities (46, 47). Identification of molecular drivers and biomarkers is therefore needed to improve the outcome of individual patients. Mutated cancer-driver genes are now employed as drug sensitivity biomarkers in clinical trials for MM, leukemia and lymphomas (48–50) (NCT04470947). Moving forward in this context, the EXALT study (NCT03096821) has demonstrated that drug-response testing in cancer cells combined with molecular profiling data is feasible and can improve treatment outcome (51), providing proof-of-concept data for this approach.
Here, we aimed to test 33 targeted drugs and to compare their sensitivity and relationships to mutational status of selected genes in a panel of MM cell lines. We also performed a comprehensive analysis correlating drug sensitivities to basal expression and phosphorylation levels of signaling proteins to explore the impact of activated mitogenic pathways on differential drug responses.
We found that RAS- and PI3K- related mutations supported MM cell sensitivity to MEK and PI3K inhibitors, respectively. We observed drug-class response differences in drugs with the same target. The PI3K pathway inhibitors temsirolimus, taselisib, pictilisib, and copanlisib were very active, whereas idelalisib and GSK2636771 showed little or no activity in this study. This might be due to isoform-specific effects of the PI3K inhibitors tested, as previously reviewed (52, 53). In addition, the impact of PIK3CA mutations on the responses to PI3K inhibitors was moderate compared to that observed for NRAS mutations on the effect of MEK inhibitors. However, the observations and prediction model might be impacted by the lower frequency of MM cell lines harboring PIK3CA mutations (2/9) compared to the corresponding frequency for NRAS mutations (4/9). Nevertheless, our data are in agreement with results from a meta-analysis case in breast cancer, where the predictive role of PIK3CA mutation status on the clinical efficacy of PI3K inhibitors remains controversial (54).
Trametinib induced the highest sensitivity across the MM cell lines, and sensitivity to both binimetinib and trametinib correlated with the presence of RAS mutations, in particular in the NRAS isoform. In contrast, we observed that expression or activation of survival members, such as Bcl-2 or phospho-Bcl-2 as well as phospho-TBK1, correlated with low responses to MEK inhibitors (i.e. trametinib, bimetinib).
With respect to RAS mutations and prediction of MEK inhibitor sensitivity, the cell line IH1 harboring a p.G12V mutation in NRAS had the highest sensitivity to trametinib, followed by the OH2 and KJON cell lines harboring NRAS p.Q61K, and KRAS p.Q61H, respectively. The presence of RAS mutations and RAS pathway activation in lung cancer was previously shown to confer the highest sensitivity to MEK inhibitors among a panel of 500 oncology drugs (55), demonstrating that RAS mutations may be indicators of sensitivity to MEK inhibitors.
Of the RAS WT cell lines, CAG and U266 show the lowest sensitivity to MEK inhibition by trametinib, whereas URVIN shows reduced sensitivity compared to RAS mutant cell lines. In contrast, VOLIN has a sensitivity to trametinib in range with mutant cell lines. This demonstrates that mutational analysis alone is insufficient to predict drug sensitivity, which is in agreement with previous reports (56, 57). Next, we found that the MM cell lines CAG and URVIN with PIK3CA and PIK3R3 mutations were the most sensitive to taselisib-induced inhibition. Taselisib is a selective PI3Kα isoform-specific PI3K inhibitor. Similarly, a study in PIK3CA positive primary uterine carcinomas cells, demonstrated a stronger growth inhibition by taselisib, when compared with WT cells, as well as in vivo tumor growth inhibition in mice with PIK3CA mutation (58). A phase II study on taselisib is ongoing for patients with PIK3CA mutation and advanced refractory solid tumors, lymphomas or MM (NCT02465060). This may indicate that patients with PIK3CA mutations may be more addicted to taselisib, and further studies on this are warranted.
Since molecularly targeted therapies act on proteins, measuring the expression and activation at this level is critical to optimize the selection of targeted therapies. We therefore looked at several relevant signaling proteins and their relationships with drug sensitivity. We found that low sensitivity to the MEK inhibitors trametinib and binimetinib observed in some MM cell lines, including CAG and U266 with RAF-1 (CRAF) and BRAF mutations (see Figure 1A), could possibly be explained by high basal expression levels of phospho-TBK1 and Bcl-2 proteins. It has been suggested that another effective strategy to enhance the efficacy of MEK inhibitors involves simultaneous targeting of proteins that are outside the RAS pathway (40). An example is TBK1, an atypical I-κB kinase family member that acts through the Ral guanine exchange factor (RalGEF) cascade to promote tumor signaling, including activation of AKT (59, 60) and NF-κB (61, 62). Targeting TBK1 pharmacologically or by mRNA knockdown induces apoptosis and reduces cell viability in a subset of acute myeloid leukemia (AML) cells with an activated MYC signaling pathway necessary for survival (63). It has been shown that combined TBK1/MEK inhibitors synergistically enhance apoptosis in several model systems with mutated RAS, including lung cancer, melanoma cells, as well as in BRAF mutant melanoma cells, resistant to MEK (62, 64–66), warranting further investigation of combined use of TBK1 and MEK inhibitors also in MM.
We found that high basal expression levels of Bcl-2 predict low drug sensitivity to the MEK inhibitors trametinib, binimetinib and cobimetinib. Several pre-clinical and clinical studies have indicated cooperative activity between MEK and Bcl-2 antagonists in solid tumors, MM, and leukemias (67, 68) (NCT03312530), (NCT02670044), (NCT04487106). We have previously also reported synergistic effects of trametinib and venetoclax in vitro in chronic lymphocytic leukemia and MM (16). This indicates that our approach can identify markers that affect drug response phenotypes, which may have clinical utility.
In summary, we provide an integrated approach that uses mutational status and profiling of intracellular (phospho) proteins to test how these markers inform on drug sensitivity of targeted treatments. However, the small sample size of the cell lines in this study limited our ability to draw clear conclusions on the predictive value of this approach. The next step will be to improve the prediction significance by including a larger cohort of samples, which will also accommodate for more extensive multiple correlations analyses. While our results suggest that MM cells are sensitive to MEK and PI3K single agents, MM patients are often treated with more than one drug at the same time. Hence it would be of interest to test the efficacy of combined MEK/PI3K inhibition in the future, and also to add conventional agents, including proteasome inhibitors.
Data availability statement
The datasets presented in this article are not readily available because the data are human identifiable and cannot be openly shared according to Norwegian legislation. Data can only be shared on a collaborative basis based on approval from the Regional Ethics Committee. Requests to access the datasets should be directed to the corresponding author.
Author contributions
SS and MG designed the research together with KT. MG and SS designed methodology, performed the experiments, analyzed and interpreted data together with KT and LS. LS performed modeling of prediction analysis. TH, SN, EC, KM, AW, EH, DV, and LM-Z provided the MM cell lines with DNA sequencing data, and performed analyses of raw sequencing data. MG wrote the paper with input from KT and SS. All authors read and commented on draft versions of the manuscript and approved the final version.
Funding
This work was supported by the Regional Health Authority for South-Eastern Norway (2015031), the Research Council of Norway Centre for Digital Life PINpOINT project (294916) and the Stiftelsen Kristian Gerhard Jebsen (grant 19). This project has also received funding from the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie grant agreement (801133), and the Research Council of Norway under the frame of ERA PerMed (322898). MG was a PhD fellow funded by the UiO: Life Science convergence environment PerCaThe when this work was performed. The exome data were produced by the Genomics Core Facility (oslo.genomics.no) as part of the Norwegian Cancer Genomics Consortium cohort.
Acknowledgments
We thank Martine Schrøder for excellent technical assistance, Thea Kristin Våtsveen for helping with cell culturing and for providing the FOLE cell line, and Alexandra Gade at the Chemical Biology Screening Platform at Centre for Molecular Medicine Norway (NCMM), University of Oslo, for assistance with drug sensitivity screens. We also thank Zhi Zhao at the Oslo Center for Biostatistics and Epidemiology, University of Oslo, for useful discussions and for suggesting his own R implementation of tree-lasso regression (available online at https://github.com/zhizuio/IPFStructPenalty). We would also like to thank Arnoldo Frigessi, from the Oslo Center for Biostatistics and Epidemiology, University of Oslo, for fruitful suggestions that helped to significantly improve the presentation of our results.
Conflict of interest
The 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.1040730/full#supplementary-material
References
1. Kumar SK, Rajkumar V, Kyle RA, van Duin M, Sonneveld P, Mateos MV, et al. Multiple myeloma. Nat Rev Dis Primers (2017) 3:17046. doi: 10.1038/nrdp.2017.46
2. Brigle K, Rogers B. Pathobiology and diagnosis of multiple myeloma. Semin Oncol Nurs (2017) 33(3):225–36. doi: 10.1016/j.soncn.2017.05.012
3. Ramakrishnan V, Kumar S. PI3K/AKT/mTOR pathway in multiple myeloma: from basic biology to clinical promise. Leuk Lymphoma (2018) 59(11):2524–34. doi: 10.1080/10428194.2017.1421760
4. Hu Y, Chen W, Wang J. Progress in the identification of gene mutations involved in multiple myeloma. Onco Targets Ther (2019) 12:4075–80. doi: 10.2147/OTT.S205922
5. Kortüm KM, Langer C, Monge J, Bruins L, Zhu YX, Shi CX, et al. Longitudinal analysis of 25 sequential sample-pairs using a custom multiple myeloma mutation sequencing panel (M3P). Ann Hematol (2015) 94(7):1205–11. doi: 10.1007/s00277-015-2344-9
6. Vo JN, Wu YM, Mishler J, Hall S, Mannan R, Wang L, et al. The genetic heterogeneity and drug resistance mechanisms of relapsed refractory multiple myeloma. Nat Commun (2022) 13(1):3750. doi: 10.1038/s41467-022-31430-0
7. Kortüm KM, Mai EK, Hanafiah NH, Shi CX, Zhu YX, Bruins L, et al. Targeted sequencing of refractory myeloma reveals a high incidence of mutations in CRBN and ras pathway genes. Blood (2016) 128(9):1226–33. doi: 10.1182/blood-2016-02-698092
8. Pawlyn C, Davies FE. Toward personalized treatment in multiple myeloma based on molecular characteristics. Blood (2019) 133(7):660–75. doi: 10.1182/blood-2018-09-825331
9. Heuck CJ, Jethava Y, Khan R, van Rhee F, Zangari M, Chavan S, et al. Inhibiting MEK in MAPK pathway-activated myeloma. Leukemia (2016) 30(4):976–80. doi: 10.1038/leu.2015.208
10. Manier S, Salem KZ, Park J, Landau DA, Getz G, Ghobrial IM. Genomic complexity of multiple myeloma and its clinical implications. Nat Rev Clin Oncol (2017) 14(2):100–13. doi: 10.1038/nrclinonc.2016.122
11. Xu J, Pfarr N, Endris V, Mai EK, Md Hanafiah NH, Lehners N, et al. Molecular signaling in multiple myeloma: association of RAS/RAF mutations and MEK/ERK pathway activation. Oncogenesis (2017) 6(5):e337. doi: 10.1038/oncsis.2017.36
12. Eichner R, Fernández-Sáiz V, Targosz BS, Bassermann F. Cross talk networks of mammalian target of rapamycin signaling with the ubiquitin proteasome system and their clinical implications in multiple myeloma. Int Rev Cell Mol Biol (2019) 343:219–97. doi: 10.1016/bs.ircmb.2018.06.001
13. Anwer F, Gee KM, Iftikhar A, Baig M, Russ AD, Saeed S, et al. Future of personalized therapy targeting aberrant signaling pathways in multiple myeloma. Clin Lymphoma Myeloma Leukemia (2019) 19(7):397–405. doi: 10.1016/j.clml.2019.03.017
14. Skånland SS, Cremaschi A, Bendiksen H, Hermansen JU, Thimiri Govinda Raj DB, Munthe LA, et al. An in vitro assay for biomarker discovery and dose prediction applied to ibrutinib plus venetoclax treatment of CLL. Leukemia (2020) 34(2):478–87. doi: 10.1038/s41375-019-0569-7
15. Lukas M, Velten B, Sellner L, Tomska K, Hüellein J, Walther T, et al. Survey of ex vivo drug combination effects in chronic lymphocytic leukemia reveals synergistic drug effects and genetic dependencies. Leukemia (2020) 34(11):2934–50. doi: 10.1038/s41375-020-0846-5
16. Melvold K, Giliberto M, Karlsen L, Ayuda-Durán P, Hanes R, Holien T, et al. Mcl-1 and bcl-xL levels predict responsiveness to dual MEK/Bcl-2 inhibition in b-cell malignancies. Mol Oncol (2021) 16(5):1153–70. doi: 10.1002/1878-0261.13153
17. Nilsson K, Bennich H, Johansson SG, Pontén J. Established immunoglobulin producing myeloma (IgE) and lymphoblastoid (IgG) cell lines from an IgE myeloma patient. Clin Exp Immunol (1970) 7(4):477–89.
18. Jackson N, Lowe J, Ball J, Bromidge E, Ling NR, Larkins S, et al. Two new IgA1-kappa plasma cell leukaemia cell lines (JJN-1 & JJN-2) which proliferate in response to b cell stimulatory factor 2. Clin Exp Immunol (1989) 75(1):93–9.
19. Børset M, Hjertner O, Yaccoby S, Epstein J, Sanderson RD. Syndecan-1 is targeted to the uropods of polarized myeloma cells where it promotes adhesion and sequesters heparin-binding proteins. Blood (2000) 96(7):2528–36. doi: 10.1182/blood.V96.7.2528
20. Burger R, Guenther A, Bakker F, Schmalzing M, Bernand S, Baum W, et al. Gp130 and ras mediated signaling in human plasma cell line INA-6: a cytokine-regulated tumor model for plasmacytoma. Hematol J (2001) 2(1):42–53. doi: 10.1038/sj.thj.6200075
21. Börset M, Waage A, Brekke OL, Helseth E. TNF and IL-6 are potent growth factors for OH-2, a novel human myeloma cell line. Eur J Haematol (1994) 53(1):31–7. doi: 10.1111/j.1600-0609.1994.tb00176.x
22. Våtsveen TK, Tian E, Kresse SH, Meza-Zepeda LA, Gabrea A, Glebov O, et al. OH-2, a hyperdiploid myeloma cell line without an IGH translocation, has a complex translocation juxtaposing MYC near MAFB and the IGK locus. Leuk Res (2009) 33(12):1670–7. doi: 10.1016/j.leukres.2009.03.001
23. Hjertner O, Hjorth-Hansen H, Börset M, Seidel C, Waage A, Sundan A. Bone morphogenetic protein-4 inhibits proliferation and induces apoptosis of multiple myeloma cells. Blood (2001) 97(2):516–22. doi: 10.1182/blood.V97.2.516
24. Våtsveen TK, Børset M, Dikic A, Tian E, Micci F, Lid AHB, et al. VOLIN and KJON–two novel hyperdiploid myeloma cell lines. Genes Chromosomes Cancer (2016) 55(11):890–901. doi: 10.1002/gcc.22388
25. Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing (2018).
26. Skånland SS. Phospho flow cytometry with fluorescent cell barcoding for single cell signaling analysis and biomarker discovery. J Vis Exp (2018) 140:58386. doi: 10.3791/58386
27. Yadav B, Pemovska T, Szwajda A, Kulesskiy E, Kontro M, Karjalainen R, et al. Quantitative scoring of differential drug sensitivity for individually optimized anticancer therapies. Sci Rep (2014) 4:5193. doi: 10.1038/srep05193
28. Metsalu T, Vilo J. ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res (2015) 43(W1):W566–70. doi: 10.1093/nar/gkv468
29. Kim S, Xing EP. Tree-guided group lasso for multi-response regression with structured sparsity, with an application to eQTL mapping. Ann Appl Stat (2012) 6(3):1095–117, 23. doi: 10.1214/12-AOAS549
30. Zhao Z, Zucknick M. Structured penalized regression for drug sensitivity prediction. J R Stat Society: Ser C (Applied Statistics) (2020) 69(3):525–45. doi: 10.1111/rssc.12400
31. Rustad EH, Dai HY, Hov H, Coward E, Beisvag V, Myklebost O, et al. BRAF V600E mutation in early-stage multiple myeloma: good response to broad acting drugs and no relation to prognosis. Blood Cancer J (2015) 5(3):e299–e. doi: 10.1038/bcj.2015.24
32. Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, et al. A global reference for human genetic variation. Nature (2015) 526(7571):68–74. doi: 10.1038/nature15393
33. Fu W, O’Connor TD, Jun G, Kang HM, Abecasis G, Leal SM, et al. Analysis of 6,515 exomes reveals the recent origin of most human protein-coding variants. Nature (2013) 493(7431):216–20. doi: 10.1038/nature11690
34. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res (2001) 29(1):308–11. doi: 10.1093/nar/29.1.308
35. Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, et al. COSMIC: mining complete cancer genomes in the catalogue of somatic mutations in cancer. Nucleic Acids Res (2011) 39(Database issue):D945–50. doi: 10.1093/nar/gkq929
36. Adzhubei I, Jordan DM, Sunyaev SR. Predicting functional effect of human missense mutations using PolyPhen-2. Curr Protoc Hum Genet (2013) Chapter 7:7.20. doi: 10.1002/0471142905.hg0720s76
37. Daver N, Schlenk RF, Russell NH, Levis MJ. Targeting FLT3 mutations in AML: review of current knowledge and evidence. Leukemia (2019) 33(2):299–312. doi: 10.1038/s41375-018-0357-9
38. Liu D, Mamorska-Dyga A. Syk inhibitors in clinical development for hematological malignancies. J Hematol Oncol (2017) 10(1):145. doi: 10.1186/s13045-016-0386-7
39. Lind J, Czernilofsky F, Vallet S, Podar K. Emerging protein kinase inhibitors for the treatment of multiple myeloma. Expert Opin Emerg Drugs (2019) 24(3):133–52. doi: 10.1080/14728214.2019.1647165
40. Molina-Arcas M, Samani A, Downward J. Drugging the undruggable: Advances on RAS targeting in cancer. Genes (2021) 12(6):899. doi: 10.3390/genes12060899
41. Le X, Antony R, Razavi P, Treacy DJ, Luo F, Ghandi M, et al. Systematic functional characterization of resistance to PI3K inhibition in breast cancer. Cancer Discovery (2016) 6(10):1134–47. doi: 10.1158/2159-8290.CD-16-0305
42. Saha MN, Jiang H, Jayakar J, Reece D, Branch DR, Chang H. MDM2 antagonist nutlin plus proteasome inhibitor velcade combination displays a synergistic anti-myeloma activity. Cancer Biol Ther (2010) 9(11):936–44. doi: 10.4161/cbt.9.11.11882
43. Tisato V, Voltan R, Gonelli A, Secchiero P, Zauli G. MDM2/X inhibitors under clinical evaluation: perspectives for the management of hematological malignancies and pediatric cancer. J Hematol Oncol (2017) 10(1):133. doi: 10.1186/s13045-017-0500-5
44. Gluck WL, Gounder MM, Frank R, Eskens F, Blay JY, Cassier PA, et al. Phase 1 study of the MDM2 inhibitor AMG 232 in patients with advanced P53 wild-type solid tumors or multiple myeloma. Invest New Drugs (2020) 38(3):831–43. doi: 10.1007/s10637-019-00840-1
45. Teoh PJ, Chung TH, Sebastian S, Choo SN, Yan J, Ng SB, et al. p53 haploinsufficiency and functional abnormalities in multiple myeloma. Leukemia (2014) 28(10):2066–74. doi: 10.1038/leu.2014.102
46. Moreau P, Kumar SK, San Miguel J, Davies F, Zamagni E, Bahlis N, et al. Treatment of relapsed and refractory multiple myeloma: recommendations from the international myeloma working group. Lancet Oncol (2021) 22(3):e105–e18. doi: 10.1016/S1470-2045(20)30756-7
47. Dimopoulos MA, Moreau P, Terpos E, Mateos M-V, Zweegman S, Cook G, et al. Multiple myeloma: EHA-ESMO clinical practice guidelines for diagnosis, treatment and follow-up. HemaSphere (2021) 5(2):e528–e. doi: 10.1097/HS9.0000000000000567
48. DiNardo CD, Cortes JE. Mutations in AML: prognostic and therapeutic implications. Hematol Am Soc Hematol Educ Program (2016) 2016(1):348–55. doi: 10.1182/asheducation-2016.1.348
49. Auclair D, Anderson KC, Avigan D, Bianchi G, Biran N, Chaudhry M, et al. The myeloma-developing regimens using genomics (MyDRUG) master protocol. J Clin Oncol (2019) 37(15_suppl):TPS8057–TPS. doi: 10.1200/JCO.2019.37.15_suppl.TPS8057
50. Mateo J, Steuten L, Aftimos P, André F, Davies M, Garralda E, et al. Delivering precision oncology to patients with cancer. Nat Med (2022) 28(4):658–65. doi: 10.1038/s41591-022-01717-2
51. Kornauth C, Pemovska T, Vladimer GI, Bayer G, Bergmann M, Eder S, et al. Functional precision medicine provides clinical benefit in advanced aggressive hematological cancers and identifies exceptional responders. Cancer Discov (2021) 12(2):372–87. doi: 10.1158/2159-8290.CD-21-0538
52. Yap TA, Bjerke L, Clarke PA, Workman P. Drugging PI3K in cancer: refining targets and therapeutic strategies. Curr Opin Pharmacol (2015) 23:98–107. doi: 10.1016/j.coph.2015.05.016
53. Yin Y, Athanasiadis P, Karlsen L, Urban A, Xu H, Murali I, et al. Functional testing to characterize and stratify PI3K inhibitor responses in chronic lymphocytic leukemia. Clin Cancer Res (2022) 28(20):4444–55. doi: 10.1158/1078-0432.CCR-22-1221
54. Wang M, Li J, Huang J, Luo M. The predictive role of PIK3CA mutation status on PI3K inhibitors in HR+ breast cancer therapy: A systematic review and meta-analysis. BioMed Res Int (2020) 2020:1598037. doi: 10.1155/2020/7451576
55. East P, Kelly GP, Biswas D, Marani M, Hancock DC, Creasy T, et al. Oncogenic RAS activity predicts response to chemotherapy and outcome in lung adenocarcinoma. bioRxiv (2021) 2021:2021.04.02.437896. doi: 10.1101/2021.04.02.437896
56. Coker EA, Stewart A, Ozer B, Minchom A, Pickard L, Ruddle R, et al. Individualized prediction of drug response and rational combination therapy in NSCLC using artificial intelligence–enabled studies of acute phosphoproteomic changes. Mol Cancer Ther (2022) 21(6):1020–9. doi: 10.1158/1535-7163.MCT-21-0442
57. Ali M, Khan SA, Wennerberg K, Aittokallio T. Global proteomics profiling improves drug sensitivity prediction: results from a multi-omics, pan-cancer modeling approach. Bioinformatics (2018) 34(8):1353–62. doi: 10.1093/bioinformatics/btx766
58. Lopez S, Schwab CL, Cocco E, Bellone S, Bonazzoli E, English DP, et al. Taselisib, a selective inhibitor of PIK3CA, is highly effective on PIK3CA-mutated and HER2/neu amplified uterine serous carcinoma in vitro and in vivo. Gynecol Oncol (2014) 135(2):312–7. doi: 10.1016/j.ygyno.2014.08.024
59. Ou YH, Torres M, Ram R, Formstecher E, Roland C, Cheng T, et al. TBK1 directly engages Akt/PKB survival signaling to support oncogenic transformation. Mol Cell (2011) 41(4):458–70. doi: 10.1016/j.molcel.2011.01.019
60. Cooper JM, Ou YH, McMillan EA, Vaden RM, Zaman A, Bodemann BO, et al. TBK1 provides context-selective support of the activated AKT/mTOR pathway in lung cancer. Cancer Res (2017) 77(18):5077–94. doi: 10.1158/0008-5472.CAN-17-0829
61. Rajurkar M, De Jesus-Monge WE, Driscoll DR, Appleman VA, Huang H, Cotton JL, et al. The activity of gli transcription factors is essential for kras-induced pancreatic tumorigenesis. Proc Natl Acad Sci (2012) 109(17):E1038. doi: 10.1073/pnas.1114168109
62. Alam M, Ansari MM, Noor S, Mohammad T, Hasan GM, Kazim SN, et al. Therapeutic targeting of TANK-binding kinase signaling towards anticancer drug development: Challenges and opportunities. Int J Biol Macromol (2022) 207:1022–37. doi: 10.1016/j.ijbiomac.2022.03.157
63. Liu S, Marneth AE, Alexe G, Walker SR, Gandler HI, Ye DQ, et al. The kinases IKBKE and TBK1 regulate MYC-dependent survival pathways through YB-1 in AML and are targets for therapy. Blood Adv (2018) 2(23):3428–42. doi: 10.1182/bloodadvances.2018016733
64. Vu HL, Aplin AE. Targeting TBK1 inhibits migration and resistance to MEK inhibitors in mutant NRAS melanoma. Mol Cancer Res (2014) 12(10):1509–19. doi: 10.1158/1541-7786.MCR-14-0204
65. Zhu Z, Aref AR, Cohoon TJ, Barbie TU, Imamura Y, Yang S, et al. Inhibition of KRAS-driven tumorigenicity by interruption of an autocrine cytokine circuit. Cancer Discovery (2014) 4(4):452–65. doi: 10.1158/2159-8290.CD-13-0646
66. Eskiocak B, McMillan EA, Mendiratta S, Kollipara RK, Zhang H, Humphries CG, et al. Biomarker accessible and chemically addressable mechanistic subtypes of BRAF melanoma. Cancer Discovery (2017) 7(8):832–51. doi: 10.1158/2159-8290.CD-16-0955
67. Iavarone C, Zervantonakis IK, Selfors LM, Palakurthi S, Liu JF, Drapkin R, et al. Combined MEK and BCL-2/X(L) inhibition is effective in high-grade serous ovarian cancer patient-derived xenograft models and BIM levels are predictive of responsiveness. Mol Cancer Ther (2019) 18(3):642–55. doi: 10.1158/1535-7163.MCT-18-0413
Keywords: drug sensitivity screening, multiple myeloma, mutations, targeted therapy, MEK, PI3K, drug response biomarkers, precision medicine
Citation: Giliberto M, Santana LM, Holien T, Misund K, Nakken S, Vodak D, Hovig E, Meza-Zepeda LA, Coward E, Waage A, Taskén K and Skånland SS (2022) Mutational analysis and protein profiling predict drug sensitivity in multiple myeloma cell lines. Front. Oncol. 12:1040730. doi: 10.3389/fonc.2022.1040730
Received: 09 September 2022; Accepted: 31 October 2022;
Published: 29 November 2022.
Edited by:
Francesco Paolo Tambaro, AORN Santobono-Pausilipon, ItalyReviewed by:
Jia-nan Gong, Chinese Academy of Medical Sciences and Peking Union Medical College, ChinaMarc Raab, German Cancer Research Center (DKFZ), Germany
Copyright © 2022 Giliberto, Santana, Holien, Misund, Nakken, Vodak, Hovig, Meza-Zepeda, Coward, Waage, Taskén and Skånland. 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: Kjetil Taskén, kjetil.tasken@medisin.uio.no