Skip to main content
  • AACR Journals
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

AACR logo

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Focus on Radiation Oncology
      • Novel Combinations
      • Reviews
      • Editors' Picks
      • "Best of" Collection
  • First Disclosures
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

  • AACR Journals
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Molecular Cancer Therapeutics
Molecular Cancer Therapeutics
  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Meeting Abstracts
    • Collections
      • COVID-19 & Cancer Resource Center
      • Focus on Radiation Oncology
      • Novel Combinations
      • Reviews
      • Editors' Picks
      • "Best of" Collection
  • First Disclosures
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • Editors' Picks
    • OnlineFirst
    • Citation
    • Author/Keyword
    • RSS Feeds
    • My Alert Summary & Preferences
  • News
    • Cancer Discovery News
  • COVID-19
  • Webinars
  • Search More

    Advanced Search

Models and Technologies

Cell Panel Profiling Reveals Conserved Therapeutic Clusters and Differentiates the Mechanism of Action of Different PI3K/mTOR, Aurora Kinase and EZH2 Inhibitors

Joost C.M. Uitdehaag, Jeroen A.D.M. de Roos, Martine B.W. Prinsen, Nicole Willemsen-Seegers, Judith R.F. de Vetter, Jelle Dylus, Antoon M. van Doornmalen, Jeffrey Kooijman, Masaaki Sawa, Suzanne J.C. van Gerwen, Jos de Man, Rogier C. Buijsman and Guido J.R. Zaman
Joost C.M. Uitdehaag
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jeroen A.D.M. de Roos
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Martine B.W. Prinsen
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicole Willemsen-Seegers
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Judith R.F. de Vetter
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jelle Dylus
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Antoon M. van Doornmalen
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jeffrey Kooijman
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Masaaki Sawa
2Carna Biosciences, Inc., Kobe, Japan.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Suzanne J.C. van Gerwen
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jos de Man
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rogier C. Buijsman
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Guido J.R. Zaman
1Netherlands Translational Research Center B.V., Kloosterstraat, the Netherlands.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: guido.zaman@ntrc.nl
DOI: 10.1158/1535-7163.MCT-16-0403 Published December 2016
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Cancer cell line panels are important tools to characterize the in vitro activity of new investigational drugs. Here, we present the inhibition profiles of 122 anticancer agents in proliferation assays with 44 or 66 genetically characterized cancer cell lines from diverse tumor tissues (Oncolines). The library includes 29 cytotoxics, 68 kinase inhibitors, and 11 epigenetic modulators. For 38 compounds this is the first comparative profiling in a cell line panel. By strictly maintaining optimized assay protocols, biological variation was kept to a minimum. Replicate profiles of 16 agents over three years show a high average Pearson correlation of 0.8 using IC50 values and 0.9 using GI50 values. Good correlations were observed with other panels. Curve fitting appears a large source of variation. Hierarchical clustering revealed 44 basic clusters, of which 26 contain compounds with common mechanisms of action, of which 9 were not reported before, including TTK, BET and two clusters of EZH2 inhibitors. To investigate unexpected clusterings, sets of BTK, Aurora and PI3K inhibitors were profiled in biochemical enzyme activity assays and surface plasmon resonance binding assays. The BTK inhibitor ibrutinib clusters with EGFR inhibitors, because it cross-reacts with EGFR. Aurora kinase inhibitors separate into two clusters, related to Aurora A or pan-Aurora selectivity. Similarly, 12 inhibitors in the PI3K/AKT/mTOR pathway separated into different clusters, reflecting biochemical selectivity (pan-PI3K, PI3Kβγδ-isoform selective or mTOR-selective). Of these, only allosteric mTOR inhibitors preferentially targeted PTEN-mutated cell lines. This shows that cell line profiling is an excellent tool for the unbiased classification of antiproliferative compounds. Mol Cancer Ther; 15(12); 3097–109. ©2016 AACR.

This article is featured in Highlights of This Issue, p. 2823

Introduction

In cell panel profiling, dose–response curves are determined of a compound in a panel of cell line proliferation assays, and this is an important tool to study the mechanism of action and selectivity of novel cancer therapies, in addition to find drug response biomarkers before initiating expensive and time-consuming experiments in animal models or trials in patients (1, 2).

Cell panel profiling has a long history as a tool in cancer drug discovery. Over the period between 1989 and 2014, tens of thousands of compounds have been profiled in the 60 cell line panel of the National Cancer Institute (NCI-60), which resulted in the discovery of bortezomib and eribulin (3, 4). In parallel, a 39-cell line panel was set up by the Japanese Foundation for Cancer Research (JFCR-39; refs. 5, 6).

With the advent of genomics, attention increasingly focused to coupling cell line biology to compound response. High-throughput cell line profiling showed that cell panels could identify patient subpopulations, for instance BRAF-mutant cell lines are in particular sensitive to BRAF inhibitors and EGFR-mutant cell lines are relatively sensitive to EGFR inhibitors (7, 8). As a result, even larger and more fully genetically characterized cell panels were developed, encompassing up to a thousand cell lines (9–11). The Genomics of Drug Sensitivity in Cancer (GDSC; ref. 9) and the Cancer Cell Line Encyclopedia (CCLE; ref. 10) analyzed 24 drugs in 1,036 cell lines and 138 drugs in 727 cell lines, respectively. Other studies include the Cancer Therapeutics Response Portal study of 481 drugs in 860 cell lines (CTRP; refs. 11, 12) and that of GlaxoSmithKline, of 19 compounds in 311 cell lines (GSK; refs. 8,13). These studies discovered, for instance, that PARP inhibitors were active in cells with the EWS–FLI1 translocation, which has led to clinical evaluation of PARP inhibitors in Ewing sarcoma (9, 14).

Despite the successes, the use of drug screening in large cell line panels has come under debate recently. As the shutdown of the NCI-60 panel was announced, many important drug candidates and selective tool inhibitors remain unscreened, making it harder to find the best tool inhibitors to probe new biological pathways (1, 15, 16). Second, a comparison of the CCLE and GDSC panels prompted a discussion about data reproducibility (17–21). It was pointed out that the rank correlation between IC50s on the same drug in the same cell lines could vary between 0.1 and 0.6, with a median of 0.28, which is poor (17). One cause could be that the GDSC maintains universal dose ranges and extrapolates IC50s if they fall outside the standard range. For example, for the potent compounds bortezomib and docetaxel, about 40% of the profiling IC50s consist of extrapolated values. However, also when the area-under-the-curve (AUC) activity measures are used (as e.g. in refs. 10, 11) which are not extrapolated, rank correlations remained poor (17). A tell-tale sign is that, after data clustering by the GDSC and CTRP, the highly related taxanes docetaxel and paclitaxel fell into two different therapeutic subsets (9, 11). This is improbable and inconsistent with earlier profilings (4, 6).

In defense of the large screens, it has been argued that using rank correlation as measure for data consistency is not appropriate, and a Pearson correlation on log IC50s should be used (20). In addition, the data serve their purpose of identifying genomic markers (20, 21). Care should be taken therefore not to dismiss these valuable data sets. Unfortunately, and in contrast to the NCI-60 and JFCR-39 screens, none of the large screens provide internal replicate data, which is essential to get insight into the limits of the reproducibility of cell line screening.

Here, we present a data set of inhibition profiles of 122 compounds on a panel of 44 or 66 cell line proliferation assays (Oncolines). Of these, 61 compounds are registered drugs, 42 are investigational drugs, while 19 are in discovery phase or are tool compounds. Compared with other profiles, this data set contains 38 inhibitors for which no data have been reported in the literature, including many inhibitors with novel mechanisms of actions, such as EZH2 inhibitors, BET inhibitors, isoform-selective PI3K inhibitors and Aurora kinase inhibitors (Table 1; Supplementary Tables S1 and S2). Compared with the larger panels, the panel shows higher data quality, as evidenced by the correlations of duplicate profiles of 16 chemically and mechanistically different inhibitors, over a period of three years. Clustering analyses of all library profiles revealed separate clusters of different therapeutic modalities, i.e., taxanes, platins, topoisomerase and EGFR inhibitors, which is highly consistent with earlier results (4, 7, 11) and new clusters comprising compounds sharing novel biochemical targets. This compound response database therefore serves as a valuable reference set that can be used to benchmark other cancer cell line panels and provides an up to date and unbiased view of therapeutic clusters available in antiproliferative therapies.

View this table:
  • View inline
  • View popup
Table 1.

Compounds tested in this study

Materials and Methods

Experimental data

Cell lines were obtained from the American Type Culture Collection (ATCC) from 2011 to 2014 (Supplemental Table S1B) and cultured in ATCC-recommended media. All experiments were carried out within nine passages of the original vials from ATCC who authenticated all cell lines by short tandem repeat analysis. Cell proliferation assays were carried out as described (22) using ATPlite 1step (Perkin Elmer). Exposure time was 72 hours for all compounds, except for epigenetic compounds, for which exposure time was 120 hours (Supplementary Table S2). Percentage growth was calculated, relative to the growth of unexposed cells, and relative IC50s were fitted using a four-parameter logistics curve. From these and seeding cell density data, absolute GI50s were calculated (3, 22). For AUC calculations, all percent-effect data within the test range were summed. If ranges were adapted, AUC data were not calculated.

In the course of three years, 122 inhibitors were tested on 44 cell lines, and after a panel extension, on 66 cell lines (Fig. 1; Supplementary Table S1). Compounds were dissolved in DMSO generally 48 hours before testing. Cisplatin and carboplatin were dissolved on the day of testing. An overview of all IC50 data is presented in Supplementary Table S3A.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Reproducibility of Oncolines cell panel data. A, Monitoring of cell replication rate for all cell lines in the panel. Each data point represents a profiling experiment. In 0.1% of cases, growth speed deviated more than a factor two from the average, in which cases profilings were redone. Twenty-two cell lines were added later to the panel (from experiment no. 122 onward). B, Overlay of doxorubicin data in cell line A375 measured at three different time points. C, Pearson correlations between replicate profiles. Indicated are compound names and days between replicates. Different bars represent correlations based on different cell response measures. AUC were not calculated when dose ranges differed between replicates. Average correlations are 0.80 (IC50s), 0.88 (reinterpreted IC50s), 0.88 (GI50), 0.75 (DSS), and 0.76 (AUC).

When compounds were inactive, their IC50 was set to 31,600 nmol/L. Finally, 10log IC50s (in nmol/L) were used for all subsequent analyses. When compounds had been profiled twice, the data of the largest panel was chosen, and, if equal, the earliest data set.

Bioinformatics analysis

To compare our profiles with each other and with literature, Pearson correlations between 10log IC50s were used (ρ; Figs. 1 and 2). All further calculations were performed in R (23). Drug-sensitivity scores were calculated in the package DSS (24). Hierarchical clustering (Fig. 3) used the Ward method and 1 − ρ as distance measure. The resulting tree was validated with multiscale bootstrap resampling (R package pvclust; ref. 25). A cutoff for the minimally significant Pearson correlation between profiles (ρmin = 0.29; P < 0.05) was obtained by converting the Pearson correlation to a t statistic (26). Then t was translated into P values with the Student t distribution function. Affinity propagation clustering (APC) was performed using the package apcluster, setting a fixed number of 43 clusters (Table 2; Supplementary Table S4A; ref. 27). Network trees were generated using the Fruchterman–Reingold algorithm, in the package igraph (Fig. 4; ref. 28). Principal component analysis was performed using the function princomp (Fig. 4; ref. 23).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Comparison of Oncolines IC50s with IC50s measured in other large scale panels. A, Correlation with overlapping data from 51 compounds and 18 cell lines from the NCI-60 panel (3, 4). B, Correlation with overlapping data from 39 compounds and 62 cell lines from the GDSC data set (9). For other correlation see Supplementary Fig. S1H. In A and B, Axes represent 10logIC50 (nmol/L). Numbers in frames indicate Pearson values. C, Variation of measured IC50s between eight compounds in five cell lines (A549, ACHN, BT-549, NCI-H460, and OVCAR-3). Data from five different panels are compared. The y-axis shows the difference between the pertinent logIC50 and the average logIC50 of that cell line–compound combination across all data sets. MTX: methotrexate. D, Pearson correlation matrix of 51 compounds analyzed in both the Oncolines and NCI-60 panels (blue: high correlation, red: negative correlation). The left triangle shows clusters and correlations using data from the Oncolines panel. The right triangle is identical to the left one, only based on NCI-60 data. Both data sets reveal similar clusters (some classes are indicated).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Hierarchical clustering of 122 compounds profiled on 44 or 66 cell lines. The dotted line represents the Pearson correlation cutoff used for identifying compound clusters (P < 0.05). The red brackets on the right indicate significant clusters according to multiscale bootstrap validation (P < 0.05).

View this table:
  • View inline
  • View popup
Table 2.

Mechanism-based compound clusters observed in this and other cell profiling studies

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Network trees of three important inhibitor classes, showing subclassifications and biochemical origins. A, EZH2 inhibitors. B, Aurora inhibitors. C, PI3K/Akt/mTOR inhibitors. In all trees, connections are drawn between compounds with a Pearson correlation ≥ 0.5. Only compounds within two connections of the investigated compounds (red) are shown. Other circle colors represent clusters as in Fig. 3. D, Biochemical selectivity of Aurora kinase inhibitors. Literature: refs. 15, 44, 49. E, Biochemical selectivity of PI3K inhibitors. Literature: refs. 45, 46, 50. F, Drug sensitivity for PTEN- and G, PIK3CA- and H,PIK3R1-mutant cells. The x-axis shows the geometrically averaged IC50 of mutant cells, divided by the geometrically averaged IC50 of the wild-type cells. The y-axis shows the significance (22). Orange: compounds that are significantly more potent in mutant cells. Blue: other compounds of interest. I, Second and third principal components (PC2 and PC3) of the variation in logIC50s in the profiling data set are related to cancer hallmarks (green).

Kinase activity assays

The inhibitory activity of compounds on Aurora A and C (Carna Biosciences, Inc.) was determined with LANCE Ultra TR-FRET assays (PerkinElmer) at KM,ATP and using ULight-labeled PLK substrate peptide (PerkinElmer, cat. no. TRF-0110; ref. 29). IC50s from 9-point dose–response curves were converted to KD using the Cheng–Prusoff equation (Fig. 4D).

Kinase binding assays

For surface plasmon resonance experiments, biotin-labelled kinases (Carna) were immobilized on streptavidin-coated chips (GE Healthcare; Cat. no. BR100531) on a Biacore T200 (GE Healthcare) to level of 4000 response units (RU) using Biacore buffer (50 mmol/L Tris pH 7.5, 0.05% (v/v) Tween-20, 150 mmol/L NaCl and 5 mmol/L MgCl2), except PI3Kδ, which was immobilized to a level of 8,000 RU. Remaining streptavidin was blocked with biocytin. The kinetic constants of the compounds were determined by single-cycle kinetic experiments as described previously (29). Data reported are the geometric averages of two independent experiments (Fig. 4D and E).

Cell line genetics

Mutation and gene copy number data of our cell lines were based on COSMIC v.75 of the GDSC (9) and were additionally validated by DNA sequence analysis of a number of important cancer genes (Supplementary Table S1C). For the analysis in Fig. 4F–H, cell lines were classified as having a “wild-type” or a “mutated” genotype, and the logIC50s in both sets were submitted to a t test, followed by a Benjamini–Hochberg correction for multiple testing (22, 30).

Processing of literature data

NCI-60 data were downloaded from https://dtp.cancer.gov/databases_tools/. All other data came from supplementary data from indicated references. IC50s were transformed to units of 10log (nmol/L), in parallel with our data. In case of NCI-60 and JFCR-39, data of multiple replicates were averaged (Fig. 2).

Results

Precision of cell line profiling

The high-throughput Oncolines cell panel has been running over the past three years with 44, and later 66 cell lines according to well-defined protocols (22). Biological variation was monitored by the average proliferation rates of each cell line, showing overall stability of growth rates over the course of three years (Fig. 1A).

Variation in IC50s was probed with the reference doxorubicin, which was chosen because it gives a potent, full sigmoidal dose–response curve in nearly all cell lines (Fig. 1B). The variation in 10logIC50s of six replicate profilings of doxorubicin was comparable with that of the NCI-60 panel, and similar to reported values in standardized high-throughput assays (Supplementary Fig. S1A; ref. 31).

Next, Pearson correlations (ρ) between replicate compound signatures were investigated (20). Any Pearson correlation above 0.6 is considered “fair” and above 0.8 is considered “near-perfect” (17). The doxorubicin replicates show ρ = 0.99 after 6 days and 0.81 after 643 days (Fig. 1C; Supplementary Table S3B). The replicate profiling of a total of 15 diverse (candidate) drugs, separated by 392 to 854 days, shows correlations between 0.58 and 0.93 (average 0.8; Fig. 1C).

Because it has been claimed that IC50 is not an optimal response metric (24, 32), we investigated if using other metrics would improve correlations. We calculated, from the same raw data, three other response measures: AUC (11, 22), DSS (24), and GI50 (3). AUC is based on percent-growth signals and does not use data-fitting. DSS is an area-based method that uses fitting. GI50 is the concentration of 50% cell growth inhibition and accounts for the starting cell density. For AUC and DSS, lower correlations are found (averages 0.76 and 0.75, Fig. 1C; Supplementary Table S3C). For AUC, this is consistent with literature (20), probably because an AUC value is more sensitive to variations in maximum and minimum signals in the assay. The fitting applied in DSS does not ameliorate this, maybe because DSS was implemented without manual curve curation. However, GI50-based correlations were higher than their IC50 equivalents (average 0.88 compared with 0.80), demonstrating that correcting for cell line growth rates improves reproducibility (32).

Next, we investigated if other factors also contribute to IC50 reproducibility. Plotting the standard deviations of replicate log IC50s shows that they depend on cell line growth rate, and on time between experimental replicates, but not on the standard deviation of cell growth rate in each cell line (Supplementary Fig. S1B–S1D). As this suggested a non-biological source of error, we studied the influence of data interpretation, by repeating curve fitting of all replicate profiles, as data interpretation rules tended to shift over the years of the experiments (Fig. 1C; see Supplementary Table S3B for data and Supplementary Fig. S1E–S1G for final rule set and some illustrative case studies). This resulted in significant improvements. After refitting, only 2 out of 15 replicates have a correlation below 0.8, and none have a correlation below 0.7. The average is 0.88, similar to the GI50-based correlations. This shows that it is feasible to have “near-perfect” (17, 20) correlations between profiles in a cell line panel, and that to increase reproducibility, more effort should be put into uniform and unambiguous computational determination of response parameters (refs. 24, 32; Fig. 1C).

Accuracy of cell line profiling

As a next step, we compared the IC50s of our cell panel profiling to those of other studies. The overlap of our 122 compounds with other cell line profiling studies is given in Supplementary Table S1A. Analysis of 10logIC50s measured on identical compounds in identical cell lines shows the best Pearson correlations with data from NCI-60 (ρ = 0.82) and JFCR-39 (ρ = 0.85), followed by data from the GDSC (ρ = 0.73), CCLE (ρ = 0.75), and GSK studies (ρ = 0.78; Fig. 2A and B; Supplementary Fig. S1H). Data from CTRP, which consist only of AUC data, were not analyzed because comparison of AUC metrics requires use of identical dose ranges, which was not the case. Given the fact that all platforms differ in read-out technology, incubation time, cell line passage, cell densities and plate formats, this shows that cell line profiling can be reproducible across platforms (20, 21).

As five cell lines (A549, ACHN, BT-549, NCI-H460, and OVCAR-3) were profiled in all these studies, we next compared the potency across platforms for eight highly active compounds that were profiled in multiple panels (Fig. 2C). For the Oncolines data, mitomycin-C and methotrexate deviate from values measured in at least three other panels (Fig. 2C). Data refitting improved concordance for mitomycin-C, but not for methotrexate (Fig. 2C). Therefore, methotrexate was reprofiled (Supplementary Table S3B). The correlation of the new data with the earlier replicate is high (ρ = 0.87). Thus, it appears that we measure a significantly more potent activity (factor of 15) of methotrexate than that presented in the literature (Fig. 2C). One reason might be that our solutions were always freshly made and that methotrexate can easily degrade when light-exposed (33, 34).

Cluster analysis of all 122 inhibitor profiles provides 26 validated therapeutic clusters

Next, we analyzed the data of 122 anticancer agents targeting all important oncogenic signaling pathways. A total of 38 of these compounds were not part of any profiling study before (Table 1 and Supplementary Table S1A). The 10log IC50 values were submitted to unsupervised hierarchical clustering (Fig. 3). When duplicate profiles are included, nearly all replicates are neighbors, showing that the data reproducibility is of such quality that we can pick up “identical” compounds (Supplementary Fig. S2A). From the tree without replicates (Fig. 3A), we wanted to isolate the minimum amount of relevant clusters. Therefore, we applied a cutoff at the correlation level that is minimally significant (ρ = 0.29; P < 0.05; see Methods and Supplementary Table S3D).

The cutoff leads to a total of 44 clusters, of which 35 contain at least two compounds, and maximally eight (Table 2). In total, 26 clusters contain two or more compounds with known and similar biochemical targeting, such as EGFR, topoisomerase, or MEK inhibitors (Fig. 3; Table 2). That compounds are grouped according to mechanism, not potency, is illustrated by the MEK inhibitors selumetinib, PD-0325901, and trametinib that have an average IC50 on the full panel of, respectively, 5.3 μmol/L, 0.9 μmol/L, and 0.49 μmol/L, but which still belong to the same cluster (Fig. 3; Table 2). We label the 26 mechanistically defined clusters as ‘highly validated clusters’. This approach to validation conceptually resembles the use of biochemical target profiles to pinpoint highly validated clusters (11). Our highly validated clusters include for instance a joint group of paclitaxel and docetaxel (Table 2), in contrast with GDSC and CTRP data, and in accordance with NCI-60 and JFCR-39 data (4, 6, 911).

Because clusters might depend on the clustering method used, we next evaluated affinity propagation clustering (APC; refs. 9, 27). This results in 46 clusters, of which 42 are highly similar to the hierarchical clustering (Table 2; Supplementary Table S4A). Also, with this method docetaxel and paclitaxel cluster together. However, cisplatin and carboplatin do not, which is in contrast to the hierarchical clustering and the literature (Fig. 3; refs. 4, 6).

To further compare the robustness of clustering across platforms, we calculated hierarchical clustering trees and correlation matrices for the 51 compounds that our data set has in common with the NCI-60 panel. Despite the differences in the platforms, clear parallels can be seen between groupings of compounds (Fig. 2D). In the 51-compound NCI-60 data, 11 “highly validated” clusters appear, i.e., that contain at least two compounds with similar biochemical mechanism. These are all represented in the clustering tree of Fig. 3 (see also Supplementary Table S4B). These clusters comprise EGFR, BRAF/MEK, mTOR, and multikinase inhibitors, as well as topoisomerase inhibitors, platins, taxanes, purine analogues and histone deacetylase (HDAC) inhibitors (Table 2). The fact that these clusters can be found in two different platforms with different read-out technology and comprising different cell lines provides good validation that they are truly distinct mechanistic approaches of cancer therapy.

Unexpected clusterings give insight into compound mechanism

The hierarchical clustering groups compounds with similar biochemical targets, but also reveals some surprises. For instance ibrutinib, an FDA-registered irreversible BTK inhibitor, clusters with the irreversibe EGFR inhibitor afatinib, and other reversible EGFR and HER2 inhibitors (Fig. 3). To further investigate this, we profiled ibrutinib on a panel of biochemical kinase assays (ref. 35; Supplementary Table S3E), which showed that aside from BTK, ibrutinib is also a potent EGFR and HER2 inhibitor, as described before (36). Because of these activities, ibrutinib behaves as an EGFR inhibitor in the cell panel.

The registered ALK/MET kinase inhibitor crizotinib clusters with danusertib, a pan-Aurora inhibitor (Fig. 3). This may be related to the high Aurora A activity of crizotinib (37) and the tyrosine kinase activity of danusertib (38), making a cluster with dual Aurora/Tyrosine kinase inhibitors.

Another paradoxical cluster is that of the ALK inhibitor ceritinib and the IGF1R inhibitors GSK-1838705A and NVP-ADW-742. Ceritinib could be expected to cluster with other ALK inhibitors such as crizotinib. However, recently it was stressed that ceritinib also has IGF1R activity (39). Therefore, this cluster seems to be based on common IGF1R inhibition. In the CTRP panel, the ALK inhibitor NVP-TAE684 was also observed to cluster with IGF1R inhibitors (11). As combined ALK/IGF1R inhibition showed synergistic effects, it was suggested that cluster proximity can predict synergy (11).

Some compounds with common biochemical mechanism are dispersed over various clusters. For instance, all alkylating agents tested, such as, e.g., dacarbazine, melphalan, and temozolomide cluster in different groups. Another group that is spread across clusters are the CDK inhibitors (PHA-793887, dinaciclib, and palbociclib).

Novel clusters distinguish epigenetic approaches such as EZH2 inhibitors

In the hierarchical clustering tree, many clusterings are seen of compounds that were profiled for the first time. Examples are the clusters of the two TTK inhibitors NMS-P715 and MPI-0479605, and the cluster of the two CHK1 inhibitors CHIR-124 and SCH-900776 (Fig. 3). This confirms that inhibition of these cell-cycle targets generates distinct cellular profiles (29).

A substantial set of epigenetic modulators was tested, including HDAC, BET, DOT1L, and EZH2 inhibitors. The HDAC inhibitors (entinostat, vorinostat, and panobinostat) form one cluster with the PARP inhibitor olaparib. This is based on their inhibition of among others hematopoietic cell lines, consistent with HDAC inhibitors being most frequently approved for use in hematological cancers (40). Furthermore, as HDAC and PARP inhibitors work synergistically in many settings (41), their proximity supports that clustering can be used to identify synergistic pairs (11).

The BET inhibitors JQ1 and I-BET-762 (also known as BRD4 inhibitors, which is one of the BET-family members) form a separate cluster close to the BCL-2 inhibitors (ABT-737 and navitoclax). All are very potent inhibitors of leukemia and lymphoma lines such as MOLT-4, Jurkat E6-1, CCRF-CEM, and SR, which is in line with the clinical trials for BET and BCL-2 inhibitors, that focus on leukemias and lymphomas (42).

Of the class of EZH2 inhibitors, we profiled five representatives, which form two clusters, the first consisting of GSK-343, EPZ-005687, and UNC-1999. The second consisting of GSK-126 and EPZ-6438 (Fig. 3). Both groups are active in many leukemic cell lines; the second EZH2 cluster also strongly inhibits the osteosarcoma line MG-63. It must be noted that the first group shares a common indazole scaffold, which is not seen in compounds from the second group (Supplementary Table S2), so the clustering might be related to a scaffold-specific off-target activity.

To further verify the EZH2 inhibitor subgroups, we constructed a network tree (Fig. 4A). This tree is an independent analysis from hierarchical clustering because every inhibitor can have more than two neighbours. For clarity, the network contains only inhibitors of which the profiles highly correlate with the EZH2 inhibitors, and correlates thereof (ρ ≥ 0.5, Fig. 4A). This tree confirms that there are two kinds of EZH2 inhibitors.

Aurora inhibitors come in two classes

A total of six Aurora-kinase inhibitors were profiled (tozasertib, AMG-900, GSK-1070916, danusertib, MK-5108, and MLN-8054), which all have been tested in patients (Table 1). The first three form a separate cluster that borders CDK and TTK inhibitors. Danusertib clusters with crizotinib (see above). The other two, MK-5108 and MLN-8054, border BET and BCL2 inhibitors. Follow-up analysis with a network tree shows two groups, one consisting of MK-5108 and MLN-8054, and one of the remaining four inhibitors (Fig. 4B).

To study if subgroups are related to isoform selectivity, we measured Aurora kinase A, B, and C selectivity in enzyme activity assays and binding assays based on surface plasmon resonance (Fig. 4D). Consistent with literature, MK-5108 and MLN-8054 are Aurora A selective inhibitors, whereas tozasertib, AMG-900, GSK-1070916, and danusertib are pan-Aurora inhibitors (Fig. 4D). Thus the two biochemical groups overlap with the groups in the network tree (Fig. 4B) and correlate with Aurora A and pan-Aurora selectivity, respectively.

Studying cell line profiles, it appears that the four inhibitors in the pan-Aurora cluster (Fig. 4D) are more active than the Aurora A inhibitors in diverse lines such as SK-N-AS, LS174T, BxPC-3 and AN3 CA, but also acute lymphoblastic leukemia (ALL) cell lines such as MOLT-4 and CCRF-CEM. This matches the observation that inhibitors of Aurora B (such as pan-Aurora inhibitors) are more effective than selective Aurora A inhibitors in ALL patients (43).

PI3K inhibitors

The PI3K/AKT/mTOR pathway contains many important drug targets. Particularly, mTOR inhibitors and the PI3Kδ inhibitor idelalisib have been approved for use as anticancer drugs, while many other compounds are in clinical trials (44, 45). We profiled a total of 12 inhibitors, of which 3 were not characterized in a larger panel before (Supplementary Table S1A). The hierarchical clustering assigns all inhibitors to the upper part of the tree, as expected for compounds active in the same pathway (Fig. 3). The only exception is the allosteric AKT inhibitor MK-2206, which profile seems most related to the platins (Fig. 3).

With regard to mTOR inhibitors, the hierarchical clustering clearly distinguishes the allosteric agents everolimus and temsirolimus (also known as rapamycin analogues or rapalogs) from two ATP-competitive inhibitors of mTOR (dactolisib and AZD-8055).

PI3K inhibitors are separated into two clusters, the first containing alpelisib, buparlisib, pictilisib, and apitolisib, the second containing duvelisib, idelalisib, and TGX-221. These groups are confirmed in a network analysis (Fig. 4C). Interestingly, buparlisib, the most advanced pan-PI3K inhibitor in the clinic (Table 1; ref. 44) shows correlation with a substantial number of profiles of non-PI3K inhibitors (Fig. 4C).

Because the catalytic subunit of PI3K exists in several isoforms, we compared the affinity for PI3Kα to δ by surface plasmon resonance. This shows that alpelisib, buparlisib, pictilisib, and apitolisib all have pan-PI3K activity and potently inhibit PI3Kα (Fig. 4E). In contrast, TGX-221, duvelisib, and idelalisib all have relatively high inhibitory activity on PI3Kβ, γ, or δ, and low inhibitory activity on PI3Kα (Fig. 4E). Thus, the two PI3K clusters in Fig. 4C correlate with PI3K isoform selectivity.

To study the biological relevance of the PI3K clusters, we characterized our cell lines for mutations in PTEN, PIK3CA, and PIK3R1, which are three of the most frequently mutated genes in cancer, and components of the PI3K pathway (46). Then we studied which of the inhibitors in our library most selectively inhibited mutant cell lines. For PTEN mutants, these are the rapalogs, followed by duvelisib (Fig. 4F). For PIK3CA-mutant cell lines, these are MK-2206, temsirolimus, and the PI3Kα-selective inhibitor alpelisib (Fig. 4G). For PIK3R1-mutant cell lines, this is the PI3Kγ/δ-selective inhibitor duvelisib (Fig. 4H). Thus different inhibitors of the PI3K pathway selectively inhibit subsets of cell lines dependent on PI3K signalling. Notably, only the association between rapalog sensitivity and PTEN-mutation achieves high significance (Fig. 4F), and indeed this is the only biomarker that has achieved success in the clinic so far (44).

Principal components analysis reveals four basic therapeutic vectors

Next, we investigated the mechanistic commonalities that connect the hierarchical clusters. Because the correlations at the root of the hierarchical clustering tree (Fig. 3) are not significant, we studied the variation in the cellular panel profiles by principal component analysis, which summarizes the variation in IC50 response profiles in a few essential vectors. Principal component 1 (PC1) captures mostly differences in potency (Supplementary Fig. S2B). In contrast, PC2 distinguishes between inhibitors of PI3K/AKT/mTOR signalling and RAS/RAF/MEK signalling (Fig. 4I). PC3 distinguishes between EGFR inhibitors and Aurora inhibitors. When the compound contributions to PC2 and PC3 are plotted, it is clear that the four quadrants are each occupied by PI3K/mTOR, EGFR, RAF/MEK, Aurora kinase, and Bcl-2 inhibitors (Fig. 4I). This is consistent with the four major branches in the clustering tree (Fig. 3) and indicates that these four are the major therapeutic vectors in our data.

Discussion

Here, we present the profiling of 122 compounds, of which many are drugs or drug candidates that were never profiled before (Table 1). Because of the controversies surrounding reproducibility of cell panel profiling, we carefully characterized variation in our assays. The metrics generated can serve as a useful benchmark for future studies (Fig. 2; Supplementary Fig. S1). We showed that, using optimized experimental and data interpretation protocols, an average correlation of 0.8 between internal control profiles is attainable (Fig. 1C), leading to data correlations across platforms of >0.7 (Fig. 2A and B), which is a clear improvement over the literature (17, 20) and which is sufficient to produce useful clusterings.

Correction for cell growth rate significantly improves correlations. In addition to that, data interpretation appears a great source of error, something which is supported by the finding that our data correlate best with NCI-60 and JFCR-39 panels, which also comprise carefully manually fitted data despite having technologically the greatest differences with our panel.

Within 122 Oncolines profiles, we identified 44 clusters with correlations >0.3. Although 44 might seem a high number of clusters, another profiling study, also using a correlation cutoff, found 49 profiles from 130 agents tested (7).

Of our 44 clusters, 26 are “highly validated” i.e., contain at least two compounds with similar biochemical mechanisms. This again shows that the molecular target class of inhibitors can predict the in vitro response profile (11, 13). In total, 17 of the highly validated clusters have been observed in earlier profilings (refs. 4, 6, 9, 11, 13; Table 2), demonstrating that cluster identities are independent of cell panel composition and represent a relevant and systematic classification of therapeutic modalities. The other 9 highly validated clusters were never described before, including TTK, multikinase ABL, WEE1/CHK1, and the subclasses of Aurora, PI3K, and EZH2 inhibitors (Table 2). This suggests that many of the novel antiproliferative approaches not only address new biochemical targets but also lead to distinctly new cell line response profiles. It is expected that if additional therapies are profiled, new clusters will appear. To facilitate their identification, we generated an overview of the 72 validated clusters identified so far in the literature (Supplementary Table S4C). It is clear that cancer cell line profiling can be used to test how distinct a new therapy is, and to get an overview of the diversity within cancer therapy.

The finding that biochemical targeting is predictive of cluster identity is supported by our analysis of the cellular profiles and biochemical selectivities of panels of Aurora and PI3K inhibitors, which shows that subclusters can relate to isoform selectivity of these inhibitor classes. Thus, there are cell line populations that are more sensitive to Aurora A or pan-Aurora inhibitors. Others are more sensitive to pan-PI3K, or PI3Kβ-, γ-, or δ- selective inhibitors. It is possible that more classes of PI3K inhibitors will appear if the same inhibitor set is tested on a larger panel, incorporating cell lines representing chronic lymphocytic leukemia, which is one of the target diseases for PI3Kγ/δ inhibitors. It must be noted that isoform selectivity is not relevant for all compounds, as for instance all non-covalent EGFR inhibitors cluster together (Fig. 3), whereas some of them also have additional ERBB2 activity (35–37) and our panel comprises a.o. the ERBB2-overexpressing line AU-565 (10).

Principal component analysis shows that our data set differentiates four major different therapy areas: PI3K/mTOR (related to metabolism), EGFR and MEK (related to growth factor signalling), Aurora (related to cell-cycle signalling), and Bcl2 (related to apoptosis). Intriguingly, this coincides with the antiproliferative subsets of hallmarks of cancer, as defined by Hanahan and Weinberg (Fig. 4I; ref. 47). Thus, these theoretical signalling classes empirically emerge from our unbiased screen, which supports the classification of cancer therapy with the hallmarks concept.

In conclusion, we have addressed controversies concerning the reproducibility of cancer cell line profiling and shown that data can contain meaningful information even if panels run over several years. The most direct proof is that compound profiles cluster according to biochemical mechanism. Clustering not only reveals side activities for known compounds, such as the EGFR activity of ibrutinib, but also shows that isoform-selective compounds are therapeutically distinguishable subclasses, as for PI3K and Aurora inhibitors. Moreover, cell panel profiling can identify compounds with synergistic effects, by combining compounds with similar cellular selectivities (48), or by combining cluster neighbours (11). Thorough cell line profiling is therefore a powerful in vitro tool for categorizing and understanding novel cancer therapeutics.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

Authors' Contributions

Conception and design: J.C.M. Uitdehaag, G.J.R. Zaman

Development of methodology: J.C.M. Uitdehaag

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J.A.D.M. de Roos, M.B.W. Prinsen, N. Willemsen-Seegers, J.R.F. de Vetter, A.M. van Doornmalen, M. Sawa

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J.C.M. Uitdehaag, J.A.D.M. de Roos, M.B.W. Prinsen, N. Willemsen-Seegers, J.R.F. de Vetter, J. Dylus, A.M. van Doornmalen, J. Kooijman, G.J.R. Zaman

Writing, review, and/or revision of the manuscript: J.C.M. Uitdehaag, J.A.D.M. de Roos, J. Dylus, S.J.C. van Gerwen, J. de Man, R.C. Buijsman, G.J.R. Zaman

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J. Dylus, A.M. van Doornmalen, S.J.C. van Gerwen, J. de Man

Study supervision: G.J.R. Zaman

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Footnotes

  • Note: Supplementary data for this article are available at Molecular Cancer Therapeutics Online (http://mct.aacrjournals.org/).

  • Received June 22, 2016.
  • Revision received August 8, 2016.
  • Accepted August 21, 2016.
  • ©2016 American Association for Cancer Research.

References

  1. 1.↵
    1. Sharma SV,
    2. Haber DA,
    3. Settleman J
    . Cell line-based platforms to evaluate the therapeutic efficacy of candidate anticancer agents. Nat Rev Cancer 2010;10:241–53.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Weinstein JN
    . Cell lines battle cancer. Nature 2012;483:544–5.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Shoemaker RH
    . The NCI human tumour cell line anticancer screen. Nat Rev Cancer 2006;6:813–23.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Holbeck SL,
    2. Collins JM,
    3. Doroshow JH
    . Analysis of Food and Drug Administration - approved anticancer agents in the NCI-60 panel of human tumor cell lines. Mol Cancer Ther 2010;9:1451–60.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Dan S,
    2. Tsunoda T,
    3. Kitahara O,
    4. Yanagawa R,
    5. Zembutsu H,
    6. Katagiri T,
    7. et al.
    An integrated database of chemosensitivity to 55 anticancer drugs and gene expression profiles of 39 human cancer cell lines. Cancer Res 2002;62:1139–47.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Nakatsu N,
    2. Nakamura T,
    3. Yamazaki K,
    4. Sadahiro S,
    5. Makuuchi H,
    6. Kanno J,
    7. et al.
    Evaluation of action mechanisms of toxic chemicals using JFCR-39, a panel of human cancer cell lines. Mol Pharm 2007;72:1171–80.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Neve RM,
    2. Chin K,
    3. Fridlyand J,
    4. Yeh J,
    5. Baehner FL,
    6. Fevr T,
    7. et al.
    A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell 2006;10:515–27.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. McDermott U,
    2. Sharma SV,
    3. Dowell L,
    4. Greninger P,
    5. Montagut C,
    6. Lamb J,
    7. et al.
    Identification of genotype-correlated sensitivity to selective kinase inhibitors by using high-throughput tumor cell line profiling. Proc Natl Acad Sci U S A 2007;104:19936–41.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Garnett MJ,
    2. Edelman EJ,
    3. Heidorn SJ,
    4. Greenman CD,
    5. Dastur A,
    6. Lau KW,
    7. et al.
    Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature 2012;483:570–5.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Barretina J,
    2. Caponigro G,
    3. Stransky N,
    4. Venkatesan K,
    5. Margolin AA,
    6. Kim S,
    7. et al.
    The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 2013;483:603–7.
    OpenUrl
  11. 11.↵
    1. Seashore-Ludlow B,
    2. Rees MG,
    3. Cheah JH,
    4. Cokol M,
    5. Price EV,
    6. Coletti ME,
    7. et al.
    Harnessing connectivity in a large scale small-molecule sensitivity dataset. Cancer Discov 2015;5:1210–23.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Rees MG,
    2. Seashore-Ludlow B,
    3. Cheah JH,
    4. Adams DJ,
    5. Price EV,
    6. Gill S,
    7. et al.
    Correlating chemical sensitivity and basal gene expression reveals mechanism of action. Nat Chem Biol 2016;12:109–16.
    OpenUrl
  13. 13.↵
    1. Greshock J,
    2. Bachman KE,
    3. Degenhardt YY,
    4. Jing J,
    5. Wen YH,
    6. Eastman S,
    7. et al.
    Molecular target class is predictive of in vitro response profile. Cancer Res 2010;70:3677–86.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Choy E,
    2. Butrynski JE,
    3. Harmon DC,
    4. Morgan JA,
    5. George S,
    6. Wagner AJ,
    7. et al.
    Phase II study of olaparib in patients with refractory Ewing sarcoma following failure of standard chemotherapy. BMC Cancer 2014;14:813.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Uitdehaag JCM,
    2. Verkaar F,
    3. Alwan H,
    4. de Man J,
    5. Buijsman RC,
    6. Zaman GJ,
    7. et al.
    A guide to picking the most selective kinase inhibitor tool compounds for pharmacological validation of drug targets. Br J Pharmacol 2011;166:858–76.
    OpenUrl
  16. 16.↵
    1. Elkins JM,
    2. Fedele V,
    3. Szklarz M,
    4. Abdul-Azeez KR,
    5. Salah E,
    6. Mikolajczyk J,
    7. et al.
    Comprehensive characterization of the published kinase inhibitor set. Nat Biotech 2016;34:95–103.
    OpenUrlPubMed
  17. 17.↵
    1. Haibe-Kains B,
    2. El-Hachem N,
    3. Birkbak NJ,
    4. Jin AC,
    5. Beck AH,
    6. Quackenbush J,
    7. et al.
    Inconsistency in large pharmacogenomic studies. Nature 2013;504:389–93.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Fallahi-Sichani M,
    2. Honarnejad S,
    3. Heiser LM,
    4. Gray JW,
    5. Sorger PK
    . Metrics other than potency reveal systematic variation in responses to cancer drugs. Nat Chem Biol 2013;9:708–14.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Fourches D,
    2. Muratov E,
    3. Tropsha A
    . Curation of chemogenomics data. Nat Chem Biol 2015;11:535.
    OpenUrl
  20. 20.↵
    Cancer Cell Line Encyclopedia Consortium & Genomics of Drug Sensitivity in Cancer Consortium. Pharmacogenomic agreement between two cancer cell line data sets. Nature 2015;528:84–7.
    OpenUrlPubMed
  21. 21.↵
    1. Haverty PM,
    2. Lin E,
    3. Tan J,
    4. Yu Y,
    5. Lam B,
    6. Lianoglou S,
    7. et al.
    Reproducible pharmacogenomic profiling of cancer cell line panels. Nature 2016;533:333–7.
    OpenUrl
  22. 22.↵
    1. Uitdehaag JCM,
    2. de Roos JADM,
    3. van Doornmalen AM,
    4. Prinsen MBW,
    5. de Man J,
    6. Tanizawa Y,
    7. et al.
    Comparison of the cancer gene targeting and biochemical selectivities of all targeted kinase inhibitors approved for clinical use. PLoS One 2014;9:e92146.
    OpenUrlCrossRefPubMed
  23. 23.↵
    R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. 2008.
  24. 24.↵
    1. Yadav B,
    2. Pemovska T,
    3. Szwajda A,
    4. Kulessiy E,
    5. Kontro M,
    6. Karjalainen R,
    7. et al.
    Quantitative scoring of differential drug sensitivity for individually optimized anticancer therapies. Sci Rep 2014;4:5193.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Suzuki R,
    2. Shimodaira H
    . Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 2006;12:1540–2.
    OpenUrl
  26. 26.↵
    1. Fisher RA
    . On the probable error of a coefficient of correlation deduced from a small sample. Metron 1921;4:3–32.
    OpenUrl
  27. 27.↵
    1. Frey BJ,
    2. Dueck D
    . Clustering by passing messages between data points. Science 2007;315:972–6.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Csardi G,
    2. Nepusz T
    . The igraph software package for complex network research. Interjournal 2006;1695.
  29. 29.↵
    1. Maia ARR,
    2. de Man J,
    3. Boon U,
    4. Janssen A,
    5. Song JY,
    6. Omerzu M,
    7. et al.
    Inhibition of the spindle assembly checkpoint kinase TTK enhances the efficacy of docetaxel in a triple negative breast cancer model. Ann Oncol 2015;26:2180–92.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Benjamini Y,
    2. Hochberg Y
    . Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statistic Soc B 1995;57:289–300.
    OpenUrl
  31. 31.↵
    1. Kalliokoski T,
    2. Kramer C,
    3. Gedeck P
    . Comparability of mixed IC50 data - a statistical analysis. PLoS One 2013;8:e61007.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Hafner M,
    2. Niepel M,
    3. Chung M,
    4. Sorger PK
    . Growth rate inhibition metrics correct for confounders in measuring sensitivity to cancer drugs. Nat Methods 2016;13:521–7.
    OpenUrl
  33. 33.↵
    1. Chatterji DC,
    2. Gallelli JF
    . Thermal and photolytic decomposition of methotrexate in aqueous solutions. J Pharm Sci 1978;67:526–31.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. D'Hondt M,
    2. Vangheluwe E,
    3. van Dorpe S,
    4. Boonen J,
    5. Bauters T,
    6. vandenBroucke J,
    7. et al.
    Stability of extemporaneously prepared cytarabine, methotrexate sodium, and methylprednisolone sodium succinate. Am J Health Syst Pharm 2012;69:232–40.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Kitagawa D,
    2. Yokota K,
    3. Gouda M,
    4. Narumi Y,
    5. Ohmoto H,
    6. Nishiwaki E,
    7. et al.
    Activity-based kinase profiling of approved tyrosine kinase inhibitors. Genes Cells 2013;18:110–22.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Liu Q,
    2. Sabnis Y,
    3. Zhao Z,
    4. Zhang T,
    5. Buhrlage SJ,
    6. Jones LH,
    7. et al.
    Developing irreversible inhibitors of the protein kinase cysteinome. Chem Biol 2013;20:146–59.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Davis MI,
    2. Hunt JP,
    3. Herrgard S,
    4. Ciceri P,
    5. Wodicak LM,
    6. Pallares G,
    7. et al.
    Comprehensive analysis of kinase inhibitor selectivity. Nat Biotech 2011;11:1046–51
    OpenUrl
  38. 38.↵
    1. Carpinelli P,
    2. Ceruti R,
    3. Giorgini ML,
    4. Cappella P,
    5. Gianellini L,
    6. Croci V,
    7. et al.
    PHA-739358, a potent inhibitor of Aurora kinases with a selective target inhibition profile relevant to cancer. Mol Cancer Ther 2007;6:3158–68.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Lovly CM,
    2. McDonald NT,
    3. Chen H,
    4. Ortiz-Cuaran S,
    5. Heukamp LC,
    6. Yan Y,
    7. et al.
    Rationale for co-targeting IGF-1R and ALK in ALK fusion positive lung cancer. Nat Med 2014;20:1027–34.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. West AC,
    2. Johnstone RW
    . New and emerging HDAC inhibitors for cancer treatment. J Clin Invest 2014;124:30–9.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Adimoolan S,
    2. Sirisawad M,
    3. Chen J,
    4. Thiemann P,
    5. Ford JM,
    6. Buggy JJ
    . HDAC inhibitor PCI-24781 decreases RAD51 expression and inhibits homologous recombination. Proc Natl Acad Sci U S A 2007;104:19482–7.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Chaidos A,
    2. Caputo V,
    3. Karadimitris A
    . Inhibition of bromodomain and extra-terminal proteins (BET) as a potential therapeutic approach in hematological malignancies: emerging preclinical and clinical evidence. Ther Adv Hematol 2015;6:128–41.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Hartsink-Segers SA,
    2. Zwaan CM,
    3. Exalto C,
    4. Luijendijk MW,
    5. Calvert VS,
    6. Petricoin EF,
    7. et al.
    Aurora kinases in childhood acute leukemia: the promise of aurora B as therapeutic target. Leukemia 2013;27:560–8.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Dienstmann R,
    2. Rodon J,
    3. Serra V,
    4. Tabernero J
    . Picking the point of inhibition: a comparative review of PI3K/AKT/mTOR pathway inhibitors. Mol Cancer Ther 2014;13:1021–31.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Yap TA,
    2. Bjerke L,
    3. Clarke PA,
    4. Workman P
    . Drugging PI3K in cancer: refining targets and therapeutic strategies. Curr Opin Pharmacol 2015;23:98–107.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Lawrence MS,
    2. Stojanov P,
    3. Mermel CH,
    4. Robinson JT,
    5. Garraway LA,
    6. Golub TR,
    7. et al.
    Discovery and saturation analysis of cancer genes across 21 tumour types. Nature 2014;505:495–501.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Hanahan D,
    2. Weinberg RA
    . Hallmarks of cancer: the next generation. Cell 2011;144:646–74.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Uitdehaag JCM,
    2. de Roos JA,
    3. van Doornmalen AM,
    4. Prinsen MBW,
    5. Spijkers-Hagelstein JA,
    6. de Vetter JR
    . Selective targeting of CTNNB1-, KRAS- or MYC-driven cell growth by combinations of existing drugs. PLoS One 2015;10:e0125021.
    OpenUrl
  49. 49.↵
    1. Carpinelli P,
    2. Moll J
    . Aurora kinases and their inhibitors: more than one target and one drug. Adv Exp Med Biol 2008;610:54–73.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Marone R,
    2. Cmiljanovic V,
    3. Giese B,
    4. Wymann MP
    . Targeting phosphoinositide 3-kinase – Moving towards therapy. Biochim Biophys Acta 2008;1784:159–85.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Molecular Cancer Therapeutics: 15 (12)
December 2016
Volume 15, Issue 12
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Editorial Board (PDF)

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Molecular Cancer Therapeutics article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Cell Panel Profiling Reveals Conserved Therapeutic Clusters and Differentiates the Mechanism of Action of Different PI3K/mTOR, Aurora Kinase and EZH2 Inhibitors
(Your Name) has forwarded a page to you from Molecular Cancer Therapeutics
(Your Name) thought you would be interested in this article in Molecular Cancer Therapeutics.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Cell Panel Profiling Reveals Conserved Therapeutic Clusters and Differentiates the Mechanism of Action of Different PI3K/mTOR, Aurora Kinase and EZH2 Inhibitors
Joost C.M. Uitdehaag, Jeroen A.D.M. de Roos, Martine B.W. Prinsen, Nicole Willemsen-Seegers, Judith R.F. de Vetter, Jelle Dylus, Antoon M. van Doornmalen, Jeffrey Kooijman, Masaaki Sawa, Suzanne J.C. van Gerwen, Jos de Man, Rogier C. Buijsman and Guido J.R. Zaman
Mol Cancer Ther December 1 2016 (15) (12) 3097-3109; DOI: 10.1158/1535-7163.MCT-16-0403

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Cell Panel Profiling Reveals Conserved Therapeutic Clusters and Differentiates the Mechanism of Action of Different PI3K/mTOR, Aurora Kinase and EZH2 Inhibitors
Joost C.M. Uitdehaag, Jeroen A.D.M. de Roos, Martine B.W. Prinsen, Nicole Willemsen-Seegers, Judith R.F. de Vetter, Jelle Dylus, Antoon M. van Doornmalen, Jeffrey Kooijman, Masaaki Sawa, Suzanne J.C. van Gerwen, Jos de Man, Rogier C. Buijsman and Guido J.R. Zaman
Mol Cancer Ther December 1 2016 (15) (12) 3097-3109; DOI: 10.1158/1535-7163.MCT-16-0403
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Disclosure of Potential Conflicts of Interest
    • Authors' Contributions
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

  • Pharmacodynamics of PI3K and MAPK Pathway Inhibitors
  • Neutrophils in HPV-Positive Penile Cancer
  • Cabozantinib Enhances Nanoparticle Delivery to Tumor Bed
Show more Models and Technologies
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook  Twitter  LinkedIn  YouTube  RSS

Articles

  • Online First
  • Current Issue
  • Past Issues
  • Meeting Abstracts

Info for

  • Authors
  • Subscribers
  • Advertisers
  • Librarians

About MCT

  • About the Journal
  • Editorial Board
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2021 by the American Association for Cancer Research.

Molecular Cancer Therapeutics
eISSN: 1538-8514
ISSN: 1535-7163

Advertisement