Chemical Systems Biology Studies of Triple Negative Breast Cancer Cell Lines

Triple negative breast cancer (TNBC) is a highly aggressive type of breast cancer that accounts for 15-20% of breast cancer cases. Targeted therapy remains to be established for TNBC that lacks estrogen receptors, progesterone receptors and human epidermal growth factor receptor HER2. This limits the therapy to traditional chemotherapy, radiation and surgery, which is only beneficial to a fraction of TNBC patients. Transcriptomics-based subtyping of TNBC into six classes, but it is unclear how the transcriptomics-based subtypes link to effective therapeutic strategies, resulting in a poor clinical prognosis in comparison to other breast cancer subtypes. Hence, there is an imminent need for identifying molecular markers and druggable targets against TNBC. This study is focused on the establishment of functional profiling of TNBC cell lines based on their drug vulnerabilities, and to identify novel druggable signaling nodes. We studied a panel of 16 TNBC cell lines using a functional profiling approach in which we measured the responses of TNBC cells to 304 oncology compounds and 355 GSK published kinase inhibitors. The clustering analysis based on overall drug-responses did not match the transcriptomics-based subtypes, suggesting the presence of extensive heterogeneity in TNBC and that the genomic or transcriptomic profiles do not always reflect the functional behavior of these cells. First, to go beyond standard anti-proliferative drug effects, we established a multiplexed readout for both cell viability and cytotoxicity. We identified many drug classes (such as anti-mitotics, anti-metabolites), which generally are assumed to have cytotoxic effects, mostly exhibited strong effects on cell viability but failed to kill the cells. However, in a subset of the cell lines, they induced a selective cell death. In those cases, we identified differential levels of protein markers linked to the cytotoxic responses (e.g. high level PAI-1 linked to anti-mitotics), suggesting their potential use in clinics for therapeutic decision. These results highlighted that simple multiplexed cell viability and cytotoxicity measurements provide more insight in cellular responses towards the treatment and thereby may help in providing better translationally predictive readouts. Second, we devised a novel drug response metric, called normalized drug response (NDR), which accounts for many kinds of screening artifacts such as signal growth rate differences in positive and negative control, as well as in drug-treated conditions. We found that the NDR metric is a time-independent method and it significantly improved the drug response curve fitting. Our NDR will be of great value in cell-based high throughput drug screening approaches as it cuts down the cost and time for the replicate experiments and further validation with cytotoxicity assay. Lastly, we used computational approach to decipher the kinase signal addiction of breast cancer cell lines by integrating vulnerabilities to kinase inhibitors and their polypharmacology data. We developed the kinase inhibition sensitivity score (KISS) to predict single and combinatorial signal addictions. For this study, we used 40 kinase inhibitors with well-defined target selectivities. With this approach, we predicted and validated novel synergistic inhibitor combinations against TNBC cells, such as dasatinib with axitinib, bosutinib with foretinib combinations for HCC1937 cells. This study suggests that drug sensitivity profiling is a powerful strategy for de-convolving cancer cell specific target addictions.N

A normalized drug response metric improves accuracy and consistency of anticancer drug sensitivity quantification in cell-based screening

Accurate quantification of drug effects is crucial for identifying pharmaceutically actionable cancer vulnerabilities. Current cell viability-based measurements often lead to biased response estimates due to varying growth rates and experimental artifacts that explain part of the inconsistency in high-throughput screening results. We developed an improved drug scoring model, normalized drug response (NDR), which makes use of both positive and negative control conditions to account for differences in cell growth rates, and experimental noise to better characterize drug-induced effects. We demonstrate an improved consistency and accuracy of NDR compared to existing metrics in assessing drug responses of cancer cells in various culture models and experimental setups. Notably, NDR reliably captures both toxicity and viability responses, and differentiates a wider spectrum of drug behavior, including lethal, growth-inhibitory and growth-stimulatory modes, based on a single viability readout. The method will therefore substantially reduce the time and resources required in cell-based drug sensitivity screening. Abhishekh Gupta et al. present a normalized drug response (NDR) metric for accurate quantification of drug sensitivity in cell-based high-throughput assays. They show that NDR captures both toxicity and viability responses to improve drug effect classification over existing methods.Peer reviewe

Multiobjective optimization identifies cancer-selective combination therapies

Author summary Cancer is diagnosed in nearly 40% of people in the U.S at some point during their lifetimes. Despite decades of research to lower cancer incidence and mortality, cancer remains a leading cause of deaths worldwide. Therefore, new targeted therapies are required to further reduce the death rates and toxic effects of treatments. Here we developed a mathematical optimization framework for finding cancer-selective treatments that optimally use drugs and their combinations. The method uses multiobjective optimization to identify drug combinations that simultaneously show maximal therapeutic potential and minimal non-selectivity, to avoid severe side effects. Our systematic search approach is applicable to various cancer types and it enables optimization of combinations involving both targeted therapies as well as standard chemotherapies. Combinatorial therapies are required to treat patients with advanced cancers that have become resistant to monotherapies through rewiring of redundant pathways. Due to a massive number of potential drug combinations, there is a need for systematic approaches to identify safe and effective combinations for each patient, using cost-effective methods. Here, we developed an exact multiobjective optimization method for identifying pairwise or higher-order combinations that show maximal cancer-selectivity. The prioritization of patient-specific combinations is based on Pareto-optimization in the search space spanned by the therapeutic and nonselective effects of combinations. We demonstrate the performance of the method in the context of BRAF-V600E melanoma treatment, where the optimal solutions predicted a number of co-inhibition partners for vemurafenib, a selective BRAF-V600E inhibitor, approved for advanced melanoma. We experimentally validated many of the predictions in BRAF-V600E melanoma cell line, and the results suggest that one can improve selective inhibition of BRAF-V600E melanoma cells by combinatorial targeting of MAPK/ERK and other compensatory pathways using pairwise and third-order drug combinations. Our mechanism-agnostic optimization method is widely applicable to various cancer types, and it takes as input only measurements of a subset of pairwise drug combinations, without requiring target information or genomic profiles. Such data-driven approaches may become useful for functional precision oncology applications that go beyond the cancer genetic dependency paradigm to optimize cancer-selective combinatorial treatments.Peer reviewe

Phenotypic Screening Combined with Machine Learning for Efficient Identification of Breast Cancer-Selective Therapeutic Targets

The lack of functional understanding of most mutations in cancer, combined with the non-druggability of most proteins, challenge genomics-based identification of oncology drug targets. We implemented a machine-learning-based approach (idTRAX), which relates cell-based screening of small-molecule compounds to their kinase inhibition data, to directly identify effective and readily druggable targets. We applied idTRAX to triple-negative breast cancer cell lines and efficiently identified cancer-selective targets. For example, we found that inhibiting AKT selectively kills MFM-223 and CAL148 cells, while inhibiting FGFR2 only kills MFM-223. Since the effects of catalytically inhibiting a protein can diverge from those of reducing its levels, targets identified by idTRAX frequently differ from those identified through gene knockout/knockdown methods. This is critical if the purpose is to identify targets specifically for small-molecule drug development, whereby idTRAX may produce fewer false-positives. The rapid nature of the approach suggests that it may be applicable in personalizing therapy.Peer reviewe

Multiobjective optimization identifies cancer-selective combination therapies

Author summaryCancer is diagnosed in nearly 40% of people in the U.S at some point during their lifetimes. Despite decades of research to lower cancer incidence and mortality, cancer remains a leading cause of deaths worldwide. Therefore, new targeted therapies are required to further reduce the death rates and toxic effects of treatments. Here we developed a mathematical optimization framework for finding cancer-selective treatments that optimally use drugs and their combinations. The method uses multiobjective optimization to identify drug combinations that simultaneously show maximal therapeutic potential and minimal non-selectivity, to avoid severe side effects. Our systematic search approach is applicable to various cancer types and it enables optimization of combinations involving both targeted therapies as well as standard chemotherapies.Combinatorial therapies are required to treat patients with advanced cancers that have become resistant to monotherapies through rewiring of redundant pathways. Due to a massive number of potential drug combinations, there is a need for systematic approaches to identify safe and effective combinations for each patient, using cost-effective methods. Here, we developed an exact multiobjective optimization method for identifying pairwise or higher-order combinations that show maximal cancer-selectivity. The prioritization of patient-specific combinations is based on Pareto-optimization in the search space spanned by the therapeutic and nonselective effects of combinations. We demonstrate the performance of the method in the context of BRAF-V600E melanoma treatment, where the optimal solutions predicted a number of co-inhibition partners for vemurafenib, a selective BRAF-V600E inhibitor, approved for advanced melanoma. We experimentally validated many of the predictions in BRAF-V600E melanoma cell line, and the results suggest that one can improve selective inhibition of BRAF-V600E melanoma cells by combinatorial targeting of MAPK/ERK and other compensatory pathways using pairwise and third-order drug combinations. Our mechanism-agnostic optimization method is widely applicable to various cancer types, and it takes as input only measurements of a subset of pairwise drug combinations, without requiring target information or genomic profiles. Such data-driven approaches may become useful for functional precision oncology applications that go beyond the cancer genetic dependency paradigm to optimize cancer-selective combinatorial treatments

Identification of selective cytotoxic and synthetic lethal drug responses in triple negative breast cancer cells

Background: Triple negative breast cancer (TNBC) is a highly heterogeneous and aggressive type of cancer that lacks effective targeted therapy. Despite detailed molecular profiling, no targeted therapy has been established. Hence, with the aim of gaining deeper understanding of the functional differences of TNBC subtypes and how that may relate to potential novel therapeutic strategies, we studied comprehensive anticancer-agent responses among a panel of TNBC cell lines. Method: The responses of 301 approved and investigational oncology compounds were measured in 16 TNBC cell lines applying a functional profiling approach. To go beyond the standard drug viability effect profiling, which has been used in most chemosensitivity studies, we utilized a multiplexed readout for both cell viability and cytotoxicity, allowing us to differentiate between cytostatic and cytotoxic responses. Results: Our approach revealed that most single-agent anti-cancer compounds that showed activity for the viability readout had no or little cytotoxic effects. Major compound classes that exhibited this type of response included anti-mitotics, mTOR, CDK, and metabolic inhibitors, as well as many agents selectively inhibiting oncogene-activated pathways. However, within the broad viability-acting classes of compounds, there were often subsets of cell lines that responded by cell death, suggesting that these cells are particularly vulnerable to the tested substance. In those cases we could identify differential levels of protein markers associated with cytotoxic responses. For example, PAI-1, MAPK phosphatase and Notch-3 levels associated with cytotoxic responses to mitotic and proteasome inhibitors, suggesting that these might serve as markers of response also in clinical settings. Furthermore, the cytotoxicity readout highlighted selective synergistic and synthetic lethal drug combinations that were missed by the cell viability readouts. For instance, the MEK inhibitor trametinib synergized with PARP inhibitors. Similarly, combination of two non-cytotoxic compounds, the rapamycin analog everolimus and an ATP-competitive mTOR inhibitor dactolisib, showed synthetic lethality in several mTOR-addicted cell lines. Conclusions: Taken together, by studying the combination of cytotoxic and cytostatic drug responses, we identified a deeper spectrum of cellular responses both to single agents and combinations that may be highly relevant for identifying precision medicine approaches in TNBC as well as in other types of cancers.Peer reviewe

Endothelial destabilization by angiopoietin-2 via integrin beta 1 activation

Peer reviewe

Multi-modal meta-analysis of cancer cell line omics profiles identifies ECHDC1 as a novel breast tumor suppressor

Molecular and functional profiling of cancer cell lines is subject to laboratory‐specific experimental practices and data analysis protocols. The current challenge therefore is how to make an integrated use of the omics profiles of cancer cell lines for reliable biological discoveries. Here, we carried out a systematic analysis of nine types of data modalities using meta‐analysis of 53 omics studies across 12 research laboratories for 2,018 cell lines. To account for a relatively low consistency observed for certain data modalities, we developed a robust data integration approach that identifies reproducible signals shared among multiple data modalities and studies. We demonstrated the power of the integrative analyses by identifying a novel driver gene, ECHDC1, with tumor suppressive role validated both in breast cancer cells and patient tumors. The multi‐modal meta‐analysis approach also identified synthetic lethal partners of cancer drivers, including a co‐dependency of PTEN deficient endometrial cancer cells on RNA helicases.</p

Concurrent Inhibition of IGF1R and ERK Increases Pancreatic Cancer Sensitivity to Autophagy Inhibitors

The aggressive nature of pancreatic ductal adenocarcinoma (PDAC) mandates the development of improved therapies. As KRAS mutations are found in 95% of PDAC and are critical for tumor maintenance, one promising strategy involves exploiting KRAS-dependent metabolic perturbations. The macrometabolic process of autophagy is upregulated in KRAS-mutant PDAC, and PDAC growth is reliant on autophagy. However, inhibition of autophagy as monotherapy using the lysosomal inhibitor hydroxychloroquine (HCQ) has shown limited clinical efficacy. To identify strategies that can improve PDAC sensitivity to HCQ, we applied a CRISPR-Cas9 loss-of-function screen and found that a top sensitizer was the receptor tyrosine kinase (RTK) insulin-like growth factor 1 receptor (IGF1R). Additionally, reverse phase protein array pathway activation mapping profiled the signaling pathways altered by chloroquine (CQ) treatment. Activating phosphorylation of RTKs, including IGF1R, was a common compensatory increase in response to CQ. Inhibition of IGF1R increased autophagic flux and sensitivity to CQ-mediated growth suppression both in vitro and in vivo. Cotargeting both IGF1R and pathways that antagonize autophagy, such as ERK-MAPK axis, was strongly synergistic. IGF1R and ERK inhibition converged on suppression of glycolysis, leading to enhanced dependence on autophagy. Accordingly, concurrent inhibition of IGF1R, ERK, and autophagy induced cytotoxicity in PDAC cell lines and decreased viability in human PDAC organoids. In conclusion, targeting IGF1R together with ERK enhances the effectiveness of autophagy inhibitors in PDAC. Significance: Compensatory upregulation of IGF1R and ERK- MAPK signaling limits the efficacy of autophagy inhibitors chloroquine and hydroxychloroquine, and their concurrent inhibition synergistically increases autophagy dependence and chloroquine sensitivity in pancreatic ductal adenocarcinoma.Peer reviewe

Network pharmacology modeling identifies synergistic Aurora B and ZAK interaction in triple-negative breast cancer

Cancer cells with heterogeneous mutation landscapes and extensive functional redundancy easily develop resistance to monotherapies by emerging activation of compensating or bypassing pathways. To achieve more effective and sustained clinical responses, synergistic interactions of multiple druggable targets that inhibit redundant cancer survival pathways are often required. Here, we report a systematic polypharmacology strategy to predict, test, and understand the selective drug combinations for MDA-MB-231 triple-negative breast cancer cells. We started by applying our network pharmacology model to predict synergistic drug combinations. Next, by utilizing kinome-wide drug-target profiles and gene expression data, we pinpointed a synergistic target interaction between Aurora B and ZAK kinase inhibition that led to enhanced growth inhibition and cytotoxicity, as validated by combinatorial siRNA, CRISPR/Cas9, and drug combination experiments. The mechanism of such a context-specific target interaction was elucidated using a dynamic simulation of MDA-MB-231 signaling network, suggesting a cross-talk between p53 and p38 pathways. Our results demonstrate the potential of polypharmacological modeling to systematically interrogate target interactions that may lead to clinically actionable and personalized treatment options.</p