The National Cancer Institute ALMANAC: A Comprehensive Screening Resource for the Detection of Anticancer Drug Pairs with Enhanced Therapeutic Activity

Association of changes in expression of HDAC and SIRT genes after drug treatment with cancer cell line sensitivity to kinase inhibitors

Histone deacetylases (HDACs) and sirtuins (SIRTs) are important epigenetic regulators of cancer pathways. There is a limited understanding of how transcriptional regulation of their genes is affected by chemotherapeutic agents, and how such transcriptional changes affect tumour sensitivity to drug treatment. We investigated the concerted transcriptional response of HDAC and SIRT genes to 15 approved antitumor agents in the NCI-60 cancer cell line panel. Antitumor agents with diverse mechanisms of action induced upregulation or downregulation of multiple HDAC and SIRT genes. HDAC5 was upregulated by dasatinib and erlotinib in the majority of the cell lines. Tumour cell line sensitivity to kinase inhibitors was associated with upregulation of HDAC5, HDAC1, and several SIRT genes. We confirmed changes in HDAC and SIRT expression in independent datasets. We also experimentally validated the upregulation of HDAC5 mRNA and protein expression by dasatinib in the highly sensitive IGROV1 cell line. HDAC5 was not upregulated in the UACC-257 cell line resistant to dasatinib. The effects of cancer drug treatment on expression of HDAC and SIRT genes may influence chemosensitivity and may need to be considered during chemotherapy.

New methods for drug synergy prediction: A mini-review

In this mini-review, we explore the new prediction methods for drug combination synergy relying on high-throughput combinatorial screens. The fast progress of the field is witnessed in the more than thirty original machine learning methods published since 2021, a clear majority of them based on deep learning techniques. We aim to put these articles under a unifying lens by highlighting the core technologies, the data sources, the input data types and synergy scores used in the methods, as well as the prediction scenarios and evaluation protocols that the articles deal with. Our finding is that the best methods accurately solve the synergy prediction scenarios involving known drugs or cell lines while the scenarios involving new drugs or cell lines still fall short of an accurate prediction level.

Large-scale Pan-cancer Cell Line Screening Identifies Actionable and Effective Drug Combinations

Oncology drug combinations can improve therapeutic responses and increase treatment options for patients. The number of possible combinations is vast and responses can be context-specific. Systematic screens can identify clinically relevant, actionable combinations in defined patient subtypes. We present data for 109 anticancer drug combinations from AstraZeneca's oncology small molecule portfolio screened in 755 pan-cancer cell lines. Combinations were screened in a 7 × 7 concentration matrix, with more than 4 million measurements of sensitivity, producing an exceptionally data-rich resource. We implement a new approach using combination Emax (viability effect) and highest single agent (HSA) to assess combination benefit. We designed a clinical translatability workflow to identify combinations with clearly defined patient populations, rationale for tolerability based on tumor type and combination-specific "emergent" biomarkers, and exposures relevant to clinical doses. We describe three actionable combinations in defined cancer types, confirmed in vitro and in vivo, with a focus on hematologic cancers and apoptotic targets.

SIGNIFICANCE

We present the largest cancer drug combination screen published to date with 7 × 7 concentration response matrices for 109 combinations in more than 750 cell lines, complemented by multi-omics predictors of response and identification of "emergent" combination biomarkers. We prioritize hits to optimize clinical translatability, and experimentally validate novel combination hypotheses. This article is featured in Selected Articles from This Issue, p. 695.

PermuteDDS: a permutable feature fusion network for drug-drug synergy prediction

MOTIVATION

Drug combination therapies have shown promise in clinical cancer treatments. However, it is hard to experimentally identify all drug combinations for synergistic interaction even with high-throughput screening due to the vast space of potential combinations. Although a number of computational methods for drug synergy prediction have proven successful in narrowing down this space, fusing drug pairs and cell line features effectively still lacks study, hindering current algorithms from understanding the complex interaction between drugs and cell lines.

RESULTS

In this paper, we proposed a Permutable feature fusion network for Drug-Drug Synergy prediction, named PermuteDDS. PermuteDDS takes multiple representations of drugs and cell lines as input and employs a permutable fusion mechanism to combine drug and cell line features. In experiments, PermuteDDS exhibits state-of-the-art performance on two benchmark data sets. Additionally, the results on independent test set grouped by different tissues reveal that PermuteDDS has good generalization performance. We believed that PermuteDDS is an effective and valuable tool for identifying synergistic drug combinations. It is publicly available at https://github.com/littlewei-lazy/PermuteDDS .

SCIENTIFIC CONTRIBUTION

First, this paper proposes a permutable feature fusion network for predicting drug synergy termed PermuteDDS, which extract diverse information from multiple drug representations and cell line representations. Second, the permutable fusion mechanism combine the drug and cell line features by integrating information of different channels, enabling the utilization of complex relationships between drugs and cell lines. Third, comparative and ablation experiments provide evidence of the efficacy of PermuteDDS in predicting drug-drug synergy.

SNSynergy: Similarity network-based machine learning framework for synergy prediction towards new cell lines and new anticancer drug combinations

The computational method has been proven to be a promising means for pre-screening large-scale anticancer drug combinations to support precision oncology applications. Pioneering efforts have been made to develop machine learning technology for predicting drug synergy, but high computational cost for training models as well as great diversity and limited size in screening data escalate the difficulty of prediction. To address this challenge, we propose a simple machine learning framework, namely Similarity Network-based Synergy prediction (SNSynergy), for predicting synergistic effects towards new cell lines and new drug combinations by two locally weighted models CLSN and DCSN. This framework only requires a small amount of auxiliary data, like genomics information of cell lines and the molecular fingerprints or targets of drugs. Based on the assumption that similar cell lines and similar drug combinations have similar synergistic effects, CLSN and DCSN predict synergy scores through capturing individual synergy contributions of nearest cell line and drug combination neighbors, respectively. High correlations between predicted and measured synergy scores on two leading cancer cell line pharmacogenomic screening datasets (the O'Neil dataset and the NCI-ALMANAC dataset) demonstrate the effectiveness and robustness of SNSynergy. Many of the identified drug combinations are consistent with previous studies, or have been explored in clinical settings against the specific cancer type, showing that SNSynergy has the potential to supply cost-saving and effective high-throughput screening for prioritizing the most applicable cell lines and the most promising drug combinations.

A Deep Neural Network for Predicting Synergistic Drug Combinations on Cancer

The exploration of drug combinations presents an opportunity to amplify therapeutic effectiveness while alleviating undesirable side effects. Nevertheless, the extensive array of potential combinations poses challenges in terms of cost and time constraints for experimental screening. Thus, it is crucial to narrow down the search space. Deep learning approaches have gained widespread popularity in predicting synergistic drug combinations tailored for specific cell lines in vitro settings. In the present study, we introduce a novel method termed GTextSyn, which utilizes the integration of gene expression data and chemical structure information for the prediction of synergistic effects in drug combinations. GTextSyn employs a sentence classification model within the domain of Natural Language Processing (NLP), wherein drugs and cell lines are regarded as entities possessing biochemical relevance. Meanwhile, combinations of drug pairs and cell lines are construed as sentences with biochemical relational significance. To assess the efficacy of GTextSyn, we conduct a comparative analysis with alternative deep learning approaches using a standard benchmark dataset. The results from a five-fold cross-validation demonstrate a 49.5% reduction in Mean Square Error (MSE) achieved by GTextSyn, surpassing the performance of the next best method in the regression task. Furthermore, we conduct a comprehensive literature survey on the predicted novel drug combinations and find substantial support from prior experimental studies for many of the combinations identified by GTextSyn.

Synergistic Drug Combinations Promote the Development of Resistance in Acute Myeloid Leukemia

Combination therapy is an important part of cancer treatment and is often employed to overcome or prevent drug resistance. Preclinical screening strategies often prioritize synergistic drug combinations; however, studies of antibiotic combinations show that synergistic drug interactions can accelerate the emergence of resistance because resistance to one drug depletes the effect of both. In this study, we aimed to determine whether synergy drives the development of resistance in cancer cell lines using live-cell imaging. Consistent with prior models of tumor evolution, we found that when controlling for activity, drug synergy is associated with increased probability of developing drug resistance. We demonstrate that these observations are an expected consequence of synergy: the fitness benefit of resisting a drug in a combination is greater in synergistic combinations than in nonsynergistic combinations. These data have important implications for preclinical strategies aiming to develop novel combinations of cancer therapies with robust and durable efficacy.

SIGNIFICANCE

Preclinical strategies to identify combinations for cancer treatment often focus on identifying synergistic combinations. This study shows that in AML cells combinations that rely on synergy can increase the likelihood of developing resistance, suggesting that combination screening strategies may benefit from a more holistic approach rather than focusing on drug synergy. See related commentary by Bhola and Letai, p. 81. This article is featured in Selected Articles from This Issue, p. 80.

Ex vivo discovery of synergistic drug combinations for hematologic malignancies

Combination therapies have improved outcomes for patients with acute myeloid leukemia (AML). However, these patients still have poor overall survival. Although many combination therapies are identified with high-throughput screening (HTS), these approaches are constrained to disease models that can be grown in large volumes (e.g., immortalized cell lines), which have limited translational utility. To identify more effective and personalized treatments, we need better strategies for screening and exploring potential combination therapies. Our objective was to develop an HTS platform for identifying effective combination therapies with highly translatable ex vivo disease models that use size-limited, primary samples from patients with leukemia (AML and myelodysplastic syndrome). We developed a system, ComboFlow, that comprises three main components: MiniFlow, ComboPooler, and AutoGater. MiniFlow conducts ex vivo drug screening with a miniaturized flow-cytometry assay that uses minimal amounts of patient sample to maximize throughput. ComboPooler incorporates computational methods to design efficient screens of pooled drug combinations. AutoGater is an automated gating classifier for flow cytometry that uses machine learning to rapidly analyze the large datasets generated by the assay. We used ComboFlow to efficiently screen more than 3000 drug combinations across 20 patient samples using only 6 million cells per patient sample. In this screen, ComboFlow identified the known synergistic combination of bortezomib and panobinostat. ComboFlow also identified a novel drug combination, dactinomycin and fludarabine, that synergistically killed leukemic cells in 35 % of AML samples. This combination also had limited effects in normal, hematopoietic progenitors. In conclusion, ComboFlow enables exploration of massive landscapes of drug combinations that were previously inaccessible in ex vivo models. We envision that ComboFlow can be used to discover more effective and personalized combination therapies for cancers amenable to ex vivo models.

Harnessing machine learning to find synergistic combinations for FDA-approved cancer drugs

Combination therapy is a fundamental strategy in cancer chemotherapy. It involves administering two or more anti-cancer agents to increase efficacy and overcome multidrug resistance compared to monotherapy. However, drug combinations can exhibit synergy, additivity, or antagonism. This study presents a machine learning framework to classify and predict cancer drug combinations. The framework utilizes several key steps including data collection and annotation from the O'Neil drug interaction dataset, data preprocessing, stratified splitting into training and test sets, construction and evaluation of classification models to categorize combinations as synergistic, additive, or antagonistic, application of regression models to predict combination sensitivity scores for enhanced predictions compared to prior work, and the last step is examination of drug features and mechanisms of action to understand synergy behaviors for optimal combinations. The models identified combination pairs most likely to synergize against different cancers. Kinase inhibitors combined with mTOR inhibitors, DNA damage-inducing drugs or HDAC inhibitors showed benefit, particularly for ovarian, melanoma, prostate, lung and colorectal carcinomas. Analysis highlighted Gemcitabine, MK-8776 and AZD1775 as frequently synergizing across cancer types. This machine learning framework provides a valuable approach to uncover more effective multi-drug regimens.

SynergyX: a multi-modality mutual attention network for interpretable drug synergy prediction

Discovering effective anti-tumor drug combinations is crucial for advancing cancer therapy. Taking full account of intricate biological interactions is highly important in accurately predicting drug synergy. However, the extremely limited prior knowledge poses great challenges in developing current computational methods. To address this, we introduce SynergyX, a multi-modality mutual attention network to improve anti-tumor drug synergy prediction. It dynamically captures cross-modal interactions, allowing for the modeling of complex biological networks and drug interactions. A convolution-augmented attention structure is adopted to integrate multi-omic data in this framework effectively. Compared with other state-of-the-art models, SynergyX demonstrates superior predictive accuracy in both the General Test and Blind Test and cross-dataset validation. By exhaustively screening combinations of approved drugs, SynergyX reveals its ability to identify promising drug combination candidates for potential lung cancer treatment. Another notable advantage lies in its multidimensional interpretability. Taking Sorafenib and Vorinostat as an example, SynergyX serves as a powerful tool for uncovering drug-gene interactions and deciphering cell selectivity mechanisms. In summary, SynergyX provides an illuminating and interpretable framework, poised to catalyze the expedition of drug synergy discovery and deepen our comprehension of rational combination therapy.

Augmented drug combination dataset to improve the performance of machine learning models predicting synergistic anticancer effects

Combination therapy has gained popularity in cancer treatment as it enhances the treatment efficacy and overcomes drug resistance. Although machine learning (ML) techniques have become an indispensable tool for discovering new drug combinations, the data on drug combination therapy currently available may be insufficient to build high-precision models. We developed a data augmentation protocol to unbiasedly scale up the existing anti-cancer drug synergy dataset. Using a new drug similarity metric, we augmented the synergy data by substituting a compound in a drug combination instance with another molecule that exhibits highly similar pharmacological effects. Using this protocol, we were able to upscale the AZ-DREAM Challenges dataset from 8798 to 6,016,697 drug combinations. Comprehensive performance evaluations show that ML models trained on the augmented data consistently achieve higher accuracy than those trained solely on the original dataset. Our data augmentation protocol provides a systematic and unbiased approach to generating more diverse and larger-scale drug combination datasets, enabling the development of more precise and effective ML models. The protocol presented in this study could serve as a foundation for future research aimed at discovering novel and effective drug combinations for cancer treatment.

Mapping combinatorial drug effects to DNA damage response kinase inhibitors

One fundamental principle that underlies various cancer treatments, such as traditional chemotherapy and radiotherapy, involves the induction of catastrophic DNA damage, leading to the apoptosis of cancer cells. In our study, we conduct a comprehensive dose-response combination screening focused on inhibitors that target key kinases involved in the DNA damage response (DDR): ATR, ATM, and DNA-PK. This screening involves 87 anti-cancer agents, including six DDR inhibitors, and encompasses 62 different cell lines spanning 12 types of tumors, resulting in a total of 17,912 combination treatment experiments. Within these combinations, we analyze the most effective and synergistic drug pairs across all tested cell lines, considering the variations among cancers originating from different tissues. Our analysis reveals inhibitors of five DDR-related pathways (DNA topoisomerase, PLK1 kinase, p53-inducible ribonucleotide reductase, PARP, and cell cycle checkpoint proteins) that exhibit strong combinatorial efficacy and synergy when used alongside ATM/ATR/DNA-PK inhibitors.

Logic-based modeling and drug repurposing for the prediction of novel therapeutic targets and combination regimens against E2F1-driven melanoma progression

Melanoma presents increasing prevalence and poor outcomes. Progression to aggressive stages is characterized by overexpression of the transcription factor E2F1 and activation of downstream prometastatic gene regulatory networks (GRNs). Appropriate therapeutic manipulation of the E2F1-governed GRNs holds the potential to prevent metastasis however, these networks entail complex feedback and feedforward regulatory motifs among various regulatory layers, which make it difficult to identify druggable components. To this end, computational approaches such as mathematical modeling and virtual screening are important tools to unveil the dynamics of these signaling networks and identify critical components that could be further explored as therapeutic targets. Herein, we integrated a well-established E2F1-mediated epithelial-mesenchymal transition (EMT) map with transcriptomics data from E2F1-expressing melanoma cells to reconstruct a core regulatory network underlying aggressive melanoma. Using logic-based in silico perturbation experiments of a core regulatory network, we identified that simultaneous perturbation of Protein kinase B (AKT1) and oncoprotein murine double minute 2 (MDM2) drastically reduces EMT in melanoma. Using the structures of the two protein signatures, virtual screening strategies were performed with the FDA-approved drug library. Furthermore, by combining drug repurposing and computer-aided drug design techniques, followed by molecular dynamics simulation analysis, we identified two potent drugs (Tadalafil and Finasteride) that can efficiently inhibit AKT1 and MDM2 proteins. We propose that these two drugs could be considered for the development of therapeutic strategies for the management of aggressive melanoma.

A granularity-level information fusion strategy on hypergraph transformer for predicting synergistic effects of anticancer drugs

Combination therapy has exhibited substantial potential compared to monotherapy. However, due to the explosive growth in the number of cancer drugs, the screening of synergistic drug combinations has become both expensive and time-consuming. Synergistic drug combinations refer to the concurrent use of two or more drugs to enhance treatment efficacy. Currently, numerous computational methods have been developed to predict the synergistic effects of anticancer drugs. However, there has been insufficient exploration of how to mine drug and cell line data at different granularity levels for predicting synergistic anticancer drug combinations. Therefore, this study proposes a granularity-level information fusion strategy based on the hypergraph transformer, named HypertranSynergy, to predict synergistic effects of anticancer drugs. HypertranSynergy introduces synergistic connections between cancer cell lines and drug combinations using hypergraph. Then, the Coarse-grained Information Extraction (CIE) module merges the hypergraph with a transformer for node embeddings. In the CIE module, Contranorm is a normalization layer that mitigates over-smoothing, while Gaussian noise addresses local information gaps. Additionally, the Fine-grained Information Extraction (FIE) module assesses fine-grained information's impact on predictions by employing similarity-aware matrices from drug/cell line features. Both CIE and FIE modules are integrated into HypertranSynergy. In addition, HypertranSynergy achieved the AUC of 0.93${\pm }$0.01 and the AUPR of 0.69${\pm }$0.02 in 5-fold cross-validation of classification task, and the RMSE of 13.77${\pm }$0.07 and the PCC of 0.81${\pm }$0.02 in 5-fold cross-validation of regression task. These results are better than most of the state-of-the-art models.

A CRISPR-drug perturbational map for identifying compounds to combine with commonly used chemotherapeutics

Combination chemotherapy is crucial for successfully treating cancer. However, the enormous number of possible drug combinations means discovering safe and effective combinations remains a significant challenge. To improve this process, we conduct large-scale targeted CRISPR knockout screens in drug-treated cells, creating a genetic map of druggable genes that sensitize cells to commonly used chemotherapeutics. We prioritize neuroblastoma, the most common extracranial pediatric solid tumor, where ~50% of high-risk patients do not survive. Our screen examines all druggable gene knockouts in 18 cell lines (10 neuroblastoma, 8 others) treated with 8 widely used drugs, resulting in 94,320 unique combination-cell line perturbations, which is comparable to the largest existing drug combination screens. Using dense drug-drug rescreening, we find that the top CRISPR-nominated drug combinations are more synergistic than standard-of-care combinations, suggesting existing combinations could be improved. As proof of principle, we discover that inhibition of PRKDC, a component of the non-homologous end-joining pathway, sensitizes high-risk neuroblastoma cells to the standard-of-care drug doxorubicin in vitro and in vivo using patient-derived xenograft (PDX) models. Our findings provide a valuable resource and demonstrate the feasibility of using targeted CRISPR knockout to discover combinations with common chemotherapeutics, a methodology with application across all cancers.

SynAI: an AI-driven cancer drugs synergism prediction platform

Summary

The SynAI solution is a flexible AI-driven drug synergism prediction solution aiming to discover potential therapeutic value of compounds in early stage. Rather than providing a finite choice of drug combination or cell lines, SynAI is capable of predicting potential drug synergism/antagonism using in silico compound SMILE (Simplified Molecular Input Line Entry System) sequences. The AI core of SynAI platform has been trained against cell lines and compound pairs listed by NCI (National Cancer Institute)-Almanac and DurgCombDB datasets. In total, the training data consists of over 1 200 000 in vitro synergism tests on 150 cancer cell lines of different organ origins. Each cell line is tested against over 6000 pairs of FDA (Food and Drug Administration) approved compound combinations. Given one or both candidate compound in SMILE sequence, SynAI is able to predict the potential Bliss score of the combined compound test with the designated cell line without the needs of compound synthetization or structural analysis; thus can significantly reduce the candidate screening costs during the compound development. SynAI platform demonstrates a comparable performance to existing methods but offers more flexibilities for data input.

Availability and implementation

The evaluation version of SynAI is freely accessible online at https://synai.crownbio.com.

Targeted Investigational Oncology Agents in the NCI-60: A Phenotypic Systems-based Resource

The NCI-60 human tumor cell line panel has proved to be a useful tool for the global cancer research community in the search for novel chemotherapeutics. The publicly available cell line characterization and compound screening data from the NCI-60 assay have significantly contributed to the understanding of cellular mechanisms targeted by new oncology agents. Signature sensitivity/resistance patterns generated for a given chemotherapeutic agent against the NCI-60 panel have long served as fingerprint presentations that encompass target information and the mechanism of action associated with the tested agent. We report the establishment of a new public NCI-60 resource based on the cell line screening of a large and growing set of 175 FDA-approved oncology drugs (AOD) plus >825 clinical and investigational oncology agents (IOA), representing a diverse set (>250) of therapeutic targets and mechanisms. This data resource is available to the public (https://ioa.cancer.gov) and includes the raw data from the screening of the IOA and AOD collection along with an extensive set of visualization and analysis tools to allow for comparative study of individual test compounds and multiple compound sets.

Potent Personalized Venetoclax Partners for Acute Myeloid Leukemia Identified by Ex Vivo Drug Screening

SUMMARY

High-throughput screens (HTS) have been utilized to assess the efficacy of single drugs against patient tumor samples with the purpose of optimizing precision therapy, but testing the synergy of drug combinations can identify the ideal second drug to add. With novel sophisticated HTS, effective venetoclax combinations can be revealed that provide the cell state, phenotype, and molecular features of the susceptible and resistant cell populations. See related article by Eide, Kurtz et al., p. 452 (14) .

Augmented drug combination dataset to improve the performance of machine learning models predicting synergistic anticancer effects

Combination therapy has gained popularity in cancer treatment as it enhances the treatment efficacy and overcomes drug resistance. Although machine learning (ML) techniques have become an indispensable tool for discovering new drug combinations, the data on drug combination therapy currently available may be insufficient to build high-precision models. We developed a data augmentation protocol to unbiasedly scale up the existing anti-cancer drug synergy dataset. Using a new drug similarity metric, we augmented the synergy data by substituting a compound in a drug combination instance with another molecule that exhibits highly similar pharmacological effects. Using this protocol, we were able to upscale the AZ-DREAM Challenges dataset from 8,798 to 6,016,697 drug combinations. Comprehensive performance evaluations show that Random Forest and Gradient Boosting Trees models trained on the augmented data achieve higher accuracy than those trained solely on the original dataset. Our data augmentation protocol provides a systematic and unbiased approach to generating more diverse and larger-scale drug combination datasets, enabling the development of more precise and effective ML models. The protocol presented in this study could serve as a foundation for future research aimed at discovering novel and effective drug combinations for cancer treatment.

RECOVER identifies synergistic drug combinations in vitro through sequential model optimization

For large libraries of small molecules, exhaustive combinatorial chemical screens become infeasible to perform when considering a range of disease models, assay conditions, and dose ranges. Deep learning models have achieved state-of-the-art results in silico for the prediction of synergy scores. However, databases of drug combinations are biased toward synergistic agents and results do not generalize out of distribution. During 5 rounds of experimentation, we employ sequential model optimization with a deep learning model to select drug combinations increasingly enriched for synergism and active against a cancer cell line-evaluating only ∼5% of the total search space. Moreover, we find that learned drug embeddings (using structural information) begin to reflect biological mechanisms. In silico benchmarking suggests search queries are ∼5-10× enriched for highly synergistic drug combinations by using sequential rounds of evaluation when compared with random selection or ∼3× when using a pretrained model.

Statistical Assessment of Drug Synergy from In Vivo Combination Studies Using Mouse Tumor Models

Drug combination therapy is a promising strategy for treating cancer; however, its efficacy and synergy require rigorous evaluation in preclinical studies before going to clinical trials. Existing methods have limited power to detect synergy in animal studies. Here, we introduce a novel approach to assess in vivo drug synergy with high sensitivity and low false discovery rate. It can accurately estimate combination index and synergy score under the Bliss independence model and the highest single agent (HSA) model without any assumption on tumor growth kinetics, study duration, data completeness and balance for tumor volume measurement. We show that our method can effectively validate in vitro drug synergy discovered from cell line assays in in vivo xenograft experiments, and can help to elucidate the mechanism of action for immune checkpoint inhibitors in syngeneic mouse models by combining an anti-PD-1 antibody and several tumor-infiltrating leukocytes depletion treatments. We provide a unified view of in vitro and in vivo synergy by presenting a parallelism between the fixed-dose in vitro and the 4-group in vivo combination studies, so they can be better designed, analyzed, and compared. We emphasize that combination index, when defined here via relative survival of tumor cells, is both dose and time dependent, and give guidelines on designing informative in vivo combination studies. We explain how to interpret and apply Bliss and HSA synergies. Finally, we provide an open-source software package named invivoSyn that enables automated analysis of in vivo synergy using our method and several other existing methods.

SIGNIFICANCE

This work presents a general solution to reliably determine in vivo drug synergy in single-dose 4-group animal combination studies.

Inhibitory Mechanisms of Lekethromycin in Dog Liver Cytochrome P450 Enzymes Based on UPLC-MS/MS Cocktail Method

Lekethromycin (LKMS) is a synthetic macrolide compound derivative intended for use as a veterinary medicine. Since there have been no in vitro studies evaluating its potential for drug-drug interactions related to cytochrome P450 (CYP450) enzymes, the effect of the inhibitory mechanisms of LKMS on CYP450 enzymes is still unclear. Thus, this study aimed to evaluate the inhibitory effects of LKMS on dog CYP450 enzymes. A cocktail approach using ultra-performance liquid chromatography-tandem mass spectrometry was conducted to investigate the inhibitory effect of LKMS on canine CYP450 enzymes. Typical probe substrates of phenacetin, coumarin, bupropion, tolbutamide, dextromethorphan, chlorzoxazone, and testosterone were used for CYP1A2, CYP2A6, CYP2B6, CYP2C9, CYP2D6, CYP2E1, and CYP3A4, respectively. This study showed that LKMS might not be a time-dependent inhibitor. LKMS inhibited CYP2A6, CYP2B6, and CYP2D6 via mixed inhibition. LKMS exhibited mixed-type inhibition against the activity of CYP2A6 with an inhibition constant (Ki) value of 135.6 μΜ. LKMS inhibited CYP2B6 in a mixed way, with Ki values of 59.44 μM. A phenotyping study based on an inhibition assay indicated that CYP2D6 contributes to the biotransformation of LKMS. A mixed inhibition of CYP2D6 with Ki values of 64.87 μM was also observed. Given that this study was performed in vitro, further in vivo studies should be conducted to identify the interaction between LKMS and canine CYP450 enzymes to provide data support for the clinical application of LKMS and the avoidance of adverse interactions between other drugs.

Inhibition of GCN2 Reveals Synergy with Cell-Cycle Regulation and Proteostasis

The integrated stress response is a signaling network comprising four branches, each sensing different cellular stressors, converging on the phosphorylation of eIF2α to downregulate global translation and initiate recovery. One of these branches includes GCN2, which senses cellular amino acid insufficiency and participates in maintaining amino acid homeostasis. Previous studies have shown that GCN2 is a viable cancer target when amino acid stress is induced by inhibiting an additional target. In this light, we screened numerous drugs for their potential to synergize with the GCN2 inhibitor TAP20. The drug sensitivity of six cancer cell lines to a panel of 25 compounds was assessed. Each compound was then combined with TAP20 at concentrations below their IC50, and the impact on cell growth was evaluated. The strongly synergistic combinations were further characterized using synergy analyses and matrix-dependent invasion assays. Inhibitors of proteostasis and the MEK-ERK pathway, as well as the pan-CDK inhibitors, flavopiridol, and seliciclib, were potently synergistic with TAP20 in two cell lines. Among their common CDK targets was CDK7, which was more selectively targeted by THZ-1 and synergized with TAP20. Moreover, these combinations were partially synergistic when assessed using matrix-dependent invasion assays. However, TAP20 alone was sufficient to restrict invasion at concentrations well below its growth-inhibitory IC50. We conclude that GCN2 inhibition can be further explored in vivo as a cancer target.

A model-agnostic framework to enhance knowledge graph-based drug combination prediction with drug-drug interaction data and supervised contrastive learning

Combination therapies have brought significant advancements to the treatment of various diseases in the medical field. However, searching for effective drug combinations remains a major challenge due to the vast number of possible combinations. Biomedical knowledge graph (KG)-based methods have shown potential in predicting effective combinations for wide spectrum of diseases, but the lack of credible negative samples has limited the prediction performance of machine learning models. To address this issue, we propose a novel model-agnostic framework that leverages existing drug-drug interaction (DDI) data as a reliable negative dataset and employs supervised contrastive learning (SCL) to transform drug embedding vectors to be more suitable for drug combination prediction. We conducted extensive experiments using various network embedding algorithms, including random walk and graph neural networks, on a biomedical KG. Our framework significantly improved performance metrics compared to the baseline framework. We also provide embedding space visualizations and case studies that demonstrate the effectiveness of our approach. This work highlights the potential of using DDI data and SCL in finding tighter decision boundaries for predicting effective drug combinations.

Review of Predicting Synergistic Drug Combinations

The prediction of drug combinations is of great clinical significance. In many diseases, such as high blood pressure, diabetes, and stomach ulcers, the simultaneous use of two or more drugs has shown clear efficacy. It has greatly reduced the progression of drug resistance. This review presents the latest applications of methods for predicting the effects of drug combinations and the bioactivity databases commonly used in drug combination prediction. These studies have played a significant role in developing precision therapy. We first describe the concept of synergy. we study various publicly available databases for drug combination prediction tasks. Next, we introduce five algorithms applied to drug combinatorial prediction, which include traditional machine learning methods, deep learning methods, mathematical methods, systems biology methods and search algorithms. In the end, we sum up the difficulties encountered in prediction models.

Interpreting the Mechanism of Synergism for Drug Combinations Using Attention-Based Hierarchical Graph Pooling

Synergistic drug combinations provide huge potentials to enhance therapeutic efficacy and to reduce adverse reactions. However, effective and synergistic drug combination prediction remains an open question because of the unknown causal disease signaling pathways. Though various deep learning (AI) models have been proposed to quantitatively predict the synergism of drug combinations, the major limitation of existing deep learning methods is that they are inherently not interpretable, which makes the conclusions of AI models untransparent to human experts, henceforth limiting the robustness of the model conclusion and the implementation ability of these models in real-world human-AI healthcare. In this paper, we develop an interpretable graph neural network (GNN) that reveals the underlying essential therapeutic targets and the mechanism of the synergy (MoS) by mining the sub-molecular network of great importance. The key point of the interpretable GNN prediction model is a novel graph pooling layer, a self-attention-based node and edge pool (henceforth SANEpool), that can compute the attention score (importance) of genes and connections based on the genomic features and topology. As such, the proposed GNN model provides a systematic way to predict and interpret the drug combination synergism based on the detected crucial sub-molecular network. Experiments on various well-adopted drug-synergy-prediction datasets demonstrate that (1) the SANEpool model has superior predictive ability to generate accurate synergy score prediction, and (2) the sub-molecular networks detected by the SANEpool are self-explainable and salient for identifying synergistic drug combinations.

Transcriptional Response to Standard AML Drugs Identifies Synergistic Combinations

Unlike genomic alterations, gene expression profiles have not been widely used to refine cancer therapies. We analyzed transcriptional changes in acute myeloid leukemia (AML) cell lines in response to standard first-line AML drugs cytarabine and daunorubicin by means of RNA sequencing. Those changes were highly cell- and treatment-specific. By comparing the changes unique to treatment-sensitive and treatment-resistant AML cells, we enriched for treatment-relevant genes. Those genes were associated with drug response-specific pathways, including calcium ion-dependent exocytosis and chromatin remodeling. Pharmacological mimicking of those changes using EGFR and MEK inhibitors enhanced the response to daunorubicin with minimum standalone cytotoxicity. The synergistic response was observed even in the cell lines beyond those used for the discovery, including a primary AML sample. Additionally, publicly available cytotoxicity data confirmed the synergistic effect of EGFR inhibitors in combination with daunorubicin in all 60 investigated cancer cell lines. In conclusion, we demonstrate the utility of treatment-evoked gene expression changes to formulate rational drug combinations. This approach could improve the standard AML therapy, especially in older patients.

Challenges and Opportunities in Developing Targeted Therapies for Triple Negative Breast Cancer

Triple negative breast cancer (TNBC) is a heterogeneous group of breast cancers characterized by their lack of estrogen receptors, progesterone receptors, and the HER2 receptor. They are more aggressive than other breast cancer subtypes, with a higher mean tumor size, higher tumor grade, the worst five-year overall survival, and the highest rates of recurrence and metastasis. Developing targeted therapies for TNBC has been a major challenge due to its heterogeneity, and its treatment still largely relies on surgery, radiation therapy, and chemotherapy. In this review article, we review the efforts in developing targeted therapies for TNBC, discuss insights gained from these efforts, and highlight potential opportunities going forward. Accumulating evidence supports TNBCs as multi-driver cancers, in which multiple oncogenic drivers promote cell proliferation and survival. In such multi-driver cancers, targeted therapies would require drug combinations that simultaneously block multiple oncogenic drivers. A strategy designed to generate mechanism-based combination targeted therapies for TNBC is discussed.

Pharmacokinetic and Pharmacodynamic Drug-Drug Interactions: Research Methods and Applications

Because of the high research and development cost of new drugs, the long development process of new drugs, and the high failure rate at later stages, combining past drugs has gradually become a more economical and attractive alternative. However, the ensuing problem of drug-drug interactions (DDIs) urgently need to be solved, and combination has attracted a lot of attention from pharmaceutical researchers. At present, DDI is often evaluated and investigated from two perspectives: pharmacodynamics and pharmacokinetics. However, in some special cases, DDI cannot be accurately evaluated from a single perspective. Therefore, this review describes and compares the current DDI evaluation methods based on two aspects: pharmacokinetic interaction and pharmacodynamic interaction. The methods summarized in this paper mainly include probe drug cocktail methods, liver microsome and hepatocyte models, static models, physiologically based pharmacokinetic models, machine learning models, in vivo comparative efficacy studies, and in vitro static and dynamic tests. This review aims to serve as a useful guide for interested researchers to promote more scientific accuracy and clinical practical use of DDI studies.

Heterocyclic Iminoquinones and Quinones from the National Cancer Institute (NCI, USA) COMPARE Analysis

This review uses the National Cancer Institute (NCI) COMPARE program to establish an extensive list of heterocyclic iminoquinones and quinones with similarities in differential growth inhibition patterns across the 60-cell line panel of the NCI Developmental Therapeutics Program (DTP). Many natural products and synthetic analogues are revealed as potential NAD(P)H:quinone oxidoreductase 1 (NQO1) substrates, through correlations to dipyridoimidazo[5,4-f]benzimidazoleiminoquinone (DPIQ), and as potential thioredoxin reductase (TrxR) inhibitors, through correlations to benzo[1,2,4]triazin-7-ones and pleurotin. The strong correlation to NQO1 infers the enzyme has a major influence on the amount of the active compound with benzo[e]perimidines, phenoxazinones, benz[f]pyrido[1,2-a]indole-6,11-quinones, seriniquinones, kalasinamide, indolequinones, and furano[2,3-b]naphthoquinones, hypothesised as prodrugs. Compounds with very strong correlations to known TrxR inhibitors had inverse correlations to the expression of both reductase enzymes, NQO1 and TrxR, including naphtho[2,3-b][1,4]oxazepane-6,11-diones, benzo[a]carbazole-1,4-diones, pyranonaphthoquinones (including kalafungin, nanaomycin A, and analogues of griseusin A), and discorhabdin C. Quinoline-5,8-dione scaffolds based on streptonigrin and lavendamycin can correlate to either reductase. Inhibitors of TrxR are not necessarily (imino)quinones, e.g., parthenolides, while oxidising moieties are essential for correlations to NQO1, as with the mitosenes. Herein, an overview of synthetic methods and biological activity of each family of heterocyclic imino(quinone) is provided.

SynBa: improved estimation of drug combination synergies with uncertainty quantification

MOTIVATION

There exists a range of different quantification frameworks to estimate the synergistic effect of drug combinations. The diversity and disagreement in estimates make it challenging to determine which combinations from a large drug screening should be proceeded with. Furthermore, the lack of accurate uncertainty quantification for those estimates precludes the choice of optimal drug combinations based on the most favourable synergistic effect.

RESULTS

In this work, we propose SynBa, a flexible Bayesian approach to estimate the uncertainty of the synergistic efficacy and potency of drug combinations, so that actionable decisions can be derived from the model outputs. The actionability is enabled by incorporating the Hill equation into SynBa, so that the parameters representing the potency and the efficacy can be preserved. Existing knowledge may be conveniently inserted due to the flexibility of the prior, as shown by the empirical Beta prior defined for the normalized maximal inhibition. Through experiments on large combination screenings and comparison against benchmark methods, we show that SynBa provides improved accuracy of dose-response predictions and better-calibrated uncertainty estimation for the parameters and the predictions.

AVAILABILITY AND IMPLEMENTATION

The code for SynBa is available at https://github.com/HaotingZhang1/SynBa. The datasets are publicly available (DOI of DREAM: 10.7303/syn4231880; DOI of the NCI-ALMANAC subset: 10.5281/zenodo.4135059).

A Comprehensive Review of Interventional Clinical Trials in Patients with Bone Metastases

Bone metastasis is one of the most important factors associated with poor prognosis for patients with prostate, breast, thyroid, and lung cancer. In the past two decades, 651 clinical trials, including 554 interventional trials, were being registered in ClinicalTrials.gov and pharma.id.informa.com to combat bone metastases from different perspectives. In this review, we comprehensively analyzed, regrouped, and discussed all the interventional trials on bone metastases. Clinical trials were re-grouped into bone-targeting agents, radiotherapy, small molecule targeted therapy, combination therapy, and others, based on the different mechanisms of action including modifying the bone microenvironment and preventing the growth of cancer cells. We also discussed the potential strategies that might improve overall survival and progression-free survival of patients with bone metastases in the future.

A landscape of response to drug combinations in non-small cell lung cancer

Combination of anti-cancer drugs is broadly seen as way to overcome the often-limited efficacy of single agents. The design and testing of combinations are however very challenging. Here we present a uniquely large dataset screening over 5000 targeted agent combinations across 81 non-small cell lung cancer cell lines. Our analysis reveals a profound heterogeneity of response across the tumor models. Notably, combinations very rarely result in a strong gain in efficacy over the range of response observable with single agents. Importantly, gain of activity over single agents is more often seen when co-targeting functionally proximal genes, offering a strategy for designing more efficient combinations. Because combinatorial effect is strongly context specific, tumor specificity should be achievable. The resource provided, together with an additional validation screen sheds light on major challenges and opportunities in building efficacious combinations against cancer and provides an opportunity for training computational models for synergy prediction.

[SMILESynergy: Anticancer drug synergy prediction based on Transformer pre-trained model]

The synergistic effect of drug combinations can solve the problem of acquired resistance to single drug therapy and has great potential for the treatment of complex diseases such as cancer. In this study, to explore the impact of interactions between different drug molecules on the effect of anticancer drugs, we proposed a Transformer-based deep learning prediction model-SMILESynergy. First, the drug text data-simplified molecular input line entry system (SMILES) were used to represent the drug molecules, and drug molecule isomers were generated through SMILES Enumeration for data augmentation. Then, the attention mechanism in the Transformer was used to encode and decode the drug molecules after data augmentation, and finally, a multi-layer perceptron (MLP) was connected to obtain the synergy value of the drugs. Experimental results showed that our model had a mean squared error of 51.34 in regression analysis, an accuracy of 0.97 in classification analysis, and better predictive performance than the DeepSynergy and MulinputSynergy models. SMILESynergy offers improved predictive performance to assist researchers in rapidly screening optimal drug combinations to improve cancer treatment outcomes.

CellMinerCDB: NCATS Is a Web-Based Portal Integrating Public Cancer Cell Line Databases for Pharmacogenomic Explorations

Major advances have been made in the field of precision medicine for treating cancer. However, many open questions remain that need to be answered to realize the goal of matching every patient with cancer to the most efficacious therapy. To facilitate these efforts, we have developed CellMinerCDB: National Center for Advancing Translational Sciences (NCATS; https://discover.nci.nih.gov/rsconnect/cellminercdb_ncats/), which makes available activity information for 2,675 drugs and compounds, including multiple nononcology drugs and 1,866 drugs and compounds unique to the NCATS. CellMinerCDB: NCATS comprises 183 cancer cell lines, with 72 unique to NCATS, including some from previously understudied tissues of origin. Multiple forms of data from different institutes are integrated, including single and combination drug activity, DNA copy number, methylation and mutation, transcriptome, protein levels, histone acetylation and methylation, metabolites, CRISPR, and miscellaneous signatures. Curation of cell lines and drug names enables cross-database (CDB) analyses. Comparison of the datasets is made possible by the overlap between cell lines and drugs across databases. Multiple univariate and multivariate analysis tools are built-in, including linear regression and LASSO. Examples have been presented here for the clinical topoisomerase I (TOP1) inhibitors topotecan and irinotecan/SN-38. This web application provides both substantial new data and significant pharmacogenomic integration, allowing exploration of interrelationships.

SIGNIFICANCE

CellMinerCDB: NCATS provides activity information for 2,675 drugs in 183 cancer cell lines and analysis tools to facilitate pharmacogenomic research and to identify determinants of response.

Analysis and modeling of cancer drug responses using cell cycle phase-specific rate effects

Identifying effective therapeutic treatment strategies is a major challenge to improving outcomes for patients with breast cancer. To gain a comprehensive understanding of how clinically relevant anti-cancer agents modulate cell cycle progression, here we use genetically engineered breast cancer cell lines to track drug-induced changes in cell number and cell cycle phase to reveal drug-specific cell cycle effects that vary across time. We use a linear chain trick (LCT) computational model, which faithfully captures drug-induced dynamic responses, correctly infers drug effects, and reproduces influences on specific cell cycle phases. We use the LCT model to predict the effects of unseen drug combinations and confirm these in independent validation experiments. Our integrated experimental and modeling approach opens avenues to assess drug responses, predict effective drug combinations, and identify optimal drug sequencing strategies.

The recent progress of deep-learning-based in silico prediction of drug combination

Drug combination therapy has become a common strategy for the treatment of complex diseases. There is an urgent need for computational methods to efficiently identify appropriate drug combinations owing to the high cost of experimental screening. In recent years, deep learning has been widely used in the field of drug discovery. Here, we provide a comprehensive review on deep-learning-based drug combination prediction algorithms from multiple aspects. Current studies highlight the flexibility of this technology in integrating multimodal data and the ability to achieve state-of-art performance; it is expected that deep-learning-based prediction of drug combinations should play an important part in future drug discovery.

Understanding of Ovarian Cancer Cell-Derived Exosome Tropism for Future Therapeutic Applications

Exosomes, a subtype of extracellular vesicles, ranging from 50 to 200 nm in diameter, and mediate cell-to-cell communication in normal biological and pathological processes. Exosomes derived from tumors have multiple functions in cancer progression, resistance, and metastasis through cancer exosome-derived tropism. However, there is no quantitative information on cancer exosome-derived tropism. Such data would be highly beneficial to guide cancer therapy by inhibiting exosome release and/or uptake. Using two fluorescent protein (mKate2) transfected ovarian cancer cell lines (OVCA4 and OVCA8), cancer exosome tropism was quantified by measuring the released exosome from ovarian cancer cells and determining the uptake of exosomes into parental ovarian cancer cells, 3D spheroids, and tumors in tumor-bearing mice. The OVCA4 cells release 50 to 200 exosomes per cell, and the OVCA8 cells do 300 to 560 per cell. The uptake of exosomes by parental ovarian cancer cells is many-fold higher than by non-parental cells. In tumor-bearing mice, most exosomes are homing to the parent cancer rather than other tissues. We successfully quantified exosome release and uptake by the parent cancer cells, further proving the tropism of cancer cell-derived exosomes. The results implied that cancer exosome tropism could provide useful information for future cancer therapeutic applications.

Predicting anti-cancer drug combination responses with a temporal cell state network model

Cancer chemotherapy combines multiple drugs, but predicting the effects of drug combinations on cancer cell proliferation remains challenging, even for simple in vitro systems. We hypothesized that by combining knowledge of single drug dose responses and cell state transition network dynamics, we could predict how a population of cancer cells will respond to drug combinations. We tested this hypothesis here using three targeted inhibitors of different cell cycle states in two different cell lines in vitro. We formulated a Markov model to capture temporal cell state transitions between different cell cycle phases, with single drug data constraining how drug doses affect transition rates. This model was able to predict the landscape of all three different pairwise drug combinations across all dose ranges for both cell lines with no additional data. While further application to different cell lines, more drugs, additional cell state networks, and more complex co-culture or in vivo systems remain, this work demonstrates how currently available or attainable information could be sufficient for prediction of drug combination response for single cell lines in vitro.

Dose-response prediction for in-vitro drug combination datasets: a probabilistic approach

In this paper we propose PIICM, a probabilistic framework for dose-response prediction in high-throughput drug combination datasets. PIICM utilizes a permutation invariant version of the intrinsic co-regionalization model for multi-output Gaussian process regression, to predict dose-response surfaces in untested drug combination experiments. Coupled with an observation model that incorporates experimental uncertainty, PIICM is able to learn from noisily observed cell-viability measurements in settings where the underlying dose-response experiments are of varying quality, utilize different experimental designs, and the resulting training dataset is sparsely observed. We show that the model can accurately predict dose-response in held out experiments, and the resulting function captures relevant features indicating synergistic interaction between drugs.

SYNDEEP: a deep learning approach for the prediction of cancer drugs synergy

Drug combinations can be the prime strategy for increasing the initial treatment options in cancer therapy. However, identifying the combinations through experimental approaches is very laborious and costly. Notably, in vitro and/or in vivo examination of all the possible combinations might not be plausible. This study presented a novel computational approach to predicting synergistic drug combinations. Specifically, the deep neural network-based binary classification was utilized to develop the model. Various physicochemical, genomic, protein-protein interaction and protein-metabolite interaction information were used to predict the synergy effects of the combinations of different drugs. The performance of the constructed model was compared with shallow neural network (SNN), k-nearest neighbors (KNN), random forest (RF), support vector machines (SVMs), and gradient boosting classifiers (GBC). Based on our findings, the proposed deep neural network model was found to be capable of predicting synergistic drug combinations with high accuracy. The prediction accuracy and AUC metrics for this model were 92.21% and 97.32% in tenfold cross-validation. According to the results, the integration of different types of physicochemical and genomics features leads to more accurate prediction of synergy in cancer drugs.

Harmonizing across datasets to improve the transferability of drug combination prediction

Combination treatment has multiple advantages over traditional monotherapy in clinics, thus becoming a target of interest for many high-throughput screening (HTS) studies, which enables the development of machine learning models predicting the response of new drug combinations. However, most existing models have been tested only within a single study, and these models cannot generalize across different datasets due to significantly variable experimental settings. Here, we thoroughly assessed the transferability issue of single-study-derived models on new datasets. More importantly, we propose a method to overcome the experimental variability by harmonizing dose-response curves of different studies. Our method improves the prediction performance of machine learning models by 184% and 1367% compared to the baseline models in intra-study and inter-study predictions, respectively, and shows consistent improvement in multiple cross-validation settings. Our study addresses the crucial question of the transferability in drug combination predictions, which is fundamental for such models to be extrapolated to new drug combination discovery and clinical applications that are de facto different datasets.

A systematic evaluation of deep learning methods for the prediction of drug synergy in cancer

One of the main obstacles to the successful treatment of cancer is the phenomenon of drug resistance. A common strategy to overcome resistance is the use of combination therapies. However, the space of possibilities is huge and efficient search strategies are required. Machine Learning (ML) can be a useful tool for the discovery of novel, clinically relevant anti-cancer drug combinations. In particular, deep learning (DL) has become a popular choice for modeling drug combination effects. Here, we set out to examine the impact of different methodological choices on the performance of multimodal DL-based drug synergy prediction methods, including the use of different input data types, preprocessing steps and model architectures. Focusing on the NCI ALMANAC dataset, we found that feature selection based on prior biological knowledge has a positive impact-limiting gene expression data to cancer or drug response-specific genes improved performance. Drug features appeared to be more predictive of drug response, with a 41% increase in coefficient of determination (R2) and 26% increase in Spearman correlation relative to a baseline model that used only cell line and drug identifiers. Molecular fingerprint-based drug representations performed slightly better than learned representations-ECFP4 fingerprints increased R2 by 5.3% and Spearman correlation by 2.8% w.r.t the best learned representations. In general, fully connected feature-encoding subnetworks outperformed other architectures. DL outperformed other ML methods by more than 35% (R2) and 14% (Spearman). Additionally, an ensemble combining the top DL and ML models improved performance by about 6.5% (R2) and 4% (Spearman). Using a state-of-the-art interpretability method, we showed that DL models can learn to associate drug and cell line features with drug response in a biologically meaningful way. The strategies explored in this study will help to improve the development of computational methods for the rational design of effective drug combinations for cancer therapy.

Stabilization of E-cadherin adhesions by COX-2/GSK3β signaling is a targetable pathway in metastatic breast cancer

Metastatic progression of epithelial cancers can be associated with epithelial-mesenchymal transition (EMT) including transcriptional inhibition of E-cadherin (CDH1) expression. Recently, EM plasticity (EMP) and E-cadherin-mediated, cluster-based metastasis and treatment resistance have become more appreciated. However, the mechanisms that maintain E-cadherin expression in this context are less understood. Through studies of inflammatory breast cancer (IBC) and a 3D tumor cell "emboli" culture paradigm, we discovered that cyclooxygenase 2 (COX-2; PTGS2), a target gene of C/EBPδ (CEBPD), or its metabolite prostaglandin E2 (PGE2) promotes protein stability of E-cadherin, β-catenin, and p120 catenin through inhibition of GSK3β. The COX-2 inhibitor celecoxib downregulated E-cadherin complex proteins and caused cell death. Coexpression of E-cadherin and COX-2 was seen in breast cancer tissues from patients with poor outcome and, along with inhibitory GSK3β phosphorylation, in patient-derived xenografts (PDX) including triple negative breast cancer (TNBC).Celecoxib alone decreased E-cadherin protein expression within xenograft tumors, though CDH1 mRNA levels increased, and reduced circulating tumor cell (CTC) clusters. In combination with paclitaxel, celecoxib attenuated or regressed lung metastases. This study has uncovered a mechanism by which metastatic breast cancer cells can maintain E-cadherin-mediated cell-to-cell adhesions and cell survival, suggesting that some patients with COX-2+/E-cadherin+ breast cancer may benefit from targeting of the PGE2 signaling pathway.

Kinome inhibition states and multiomics data enable prediction of cell viability in diverse cancer types

Protein kinases play a vital role in a wide range of cellular processes, and compounds that inhibit kinase activity emerging as a primary focus for targeted therapy development, especially in cancer. Consequently, efforts to characterize the behavior of kinases in response to inhibitor treatment, as well as downstream cellular responses, have been performed at increasingly large scales. Previous work with smaller datasets have used baseline profiling of cell lines and limited kinome profiling data to attempt to predict small molecule effects on cell viability, but these efforts did not use multi-dose kinase profiles and achieved low accuracy with very limited external validation. This work focuses on two large-scale primary data types, kinase inhibitor profiles and gene expression, to predict the results of cell viability screening. We describe the process by which we combined these data sets, examined their properties in relation to cell viability and finally developed a set of computational models that achieve a reasonably high prediction accuracy (R2 of 0.78 and RMSE of 0.154). Using these models, we identified a set of kinases, several of which are understudied, that are strongly influential in the cell viability prediction models. In addition, we also tested to see if a wider range of multiomics data sets could improve the model results and found that proteomic kinase inhibitor profiles were the single most informative data type. Finally, we validated a small subset of the model predictions in several triple-negative and HER2 positive breast cancer cell lines demonstrating that the model performs well with compounds and cell lines that were not included in the training data set. Overall, this result demonstrates that generic knowledge of the kinome is predictive of very specific cell phenotypes, and has the potential to be integrated into targeted therapy development pipelines.

CCSynergy: an integrative deep-learning framework enabling context-aware prediction of anti-cancer drug synergy

Combination therapy is a promising strategy for confronting the complexity of cancer. However, experimental exploration of the vast space of potential drug combinations is costly and unfeasible. Therefore, computational methods for predicting drug synergy are much needed for narrowing down this space, especially when examining new cellular contexts. Here, we thus introduce CCSynergy, a flexible, context aware and integrative deep-learning framework that we have established to unleash the potential of the Chemical Checker extended drug bioactivity profiles for the purpose of drug synergy prediction. We have shown that CCSynergy enables predictions of superior accuracy, remarkable robustness and improved context generalizability as compared to the state-of-the-art methods in the field. Having established the potential of CCSynergy for generating experimentally validated predictions, we next exhaustively explored the untested drug combination space. This resulted in a compendium of potentially synergistic drug combinations on hundreds of cancer cell lines, which can guide future experimental screens.

Predicting cell line-specific synergistic drug combinations through a relational graph convolutional network with attention mechanism

Identifying synergistic drug combinations (SDCs) is a great challenge due to the combinatorial complexity and the fact that SDC is cell line specific. The existing computational methods either did not consider the cell line specificity of SDC, or did not perform well by building model for each cell line independently. In this paper, we present a novel encoder-decoder network named SDCNet for predicting cell line-specific SDCs. SDCNet learns common patterns across different cell lines as well as cell line-specific features in one model for drug combinations. This is realized by considering the SDC graphs of different cell lines as a relational graph, and constructing a relational graph convolutional network (R-GCN) as the encoder to learn and fuse the deep representations of drugs for different cell lines. An attention mechanism is devised to integrate the drug features from different layers of the R-GCN according to their relative importance so that representation learning is further enhanced. The common patterns are exploited through partial parameter sharing in cell line-specific decoders, which not only reconstruct the known SDCs but also predict new ones for each cell line. Experiments on various datasets demonstrate that SDCNet is superior to state-of-the-art methods and is also robust when generalized to new cell lines that are different from the training ones. Finally, the case study again confirms the effectiveness of our method in predicting novel reliable cell line-specific SDCs.

Accelerating space radiation countermeasure development through drug repurposing

The discovery of safe and effective radiation countermeasures (MCM) for long-duration spaceflight is challenging due to the complexity of the space radiation biology and high safety requirements. There are few if any clinically-validated molecular targets for this use case, and preclinical models have several known limitations. These challenges make the evaluation of existing FDA-approved drugs for this indication, or drug repurposing, an attractive strategy to accelerate space radiation countermeasure development. Drug repurposing offers several advantages over de novo drug discovery including established manufacturing methods, human clinical safety data, and well-understood dosing and pharmacokinetic considerations. There are limitations working with a fixed set of possible candidate compounds, but some properties of repurposed drugs can be tailored for well-defined new indications through reformulation and development of drug combinations. Drug repurposing is thus an attractive strategy for mitigating the high risks and costs of drug development and delivering new countermeasures to protect human from space radiation in long-term missions.

SYNPRED: prediction of drug combination effects in cancer using different synergy metrics and ensemble learning

Background

In cancer research, high-throughput screening technologies produce large amounts of multiomics data from different populations and cell types. However, analysis of such data encounters difficulties due to disease heterogeneity, further exacerbated by human biological complexity and genomic variability. The specific profile of cancer as a disease (or, more realistically, a set of diseases) urges the development of approaches that maximize the effect while minimizing the dosage of drugs. Now is the time to redefine the approach to drug discovery, bringing an artificial intelligence (AI)–powered informational view that integrates the relevant scientific fields and explores new territories.

Results

Here, we show SYNPRED, an interdisciplinary approach that leverages specifically designed ensembles of AI algorithms, as well as links omics and biophysical traits to predict anticancer drug synergy. It uses 5 reference models (Bliss, Highest Single Agent, Loewe, Zero Interaction Potency, and Combination Sensitivity Score), which, coupled with AI algorithms, allowed us to attain the ones with the best predictive performance and pinpoint the most appropriate reference model for synergy prediction, often overlooked in similar studies. By using an independent test set, SYNPRED exhibits state-of-the-art performance metrics either in the classification (accuracy, 0.85; precision, 0.91; recall, 0.90; area under the receiver operating characteristic, 0.80; and F1-score, 0.91) or in the regression models, mainly when using the Combination Sensitivity Score synergy reference model (root mean square error, 11.07; mean squared error, 122.61; Pearson, 0.86; mean absolute error, 7.43; Spearman, 0.87). Moreover, data interpretability was achieved by deploying the most current and robust feature importance approaches. A simple web-based application was constructed, allowing easy access by nonexpert researchers.

Conclusions

The performance of SYNPRED rivals that of the existing methods that tackle the same problem, yielding unbiased results trained with one of the most comprehensive datasets available (NCI ALMANAC). The leveraging of different reference models allowed deeper insights into which of them can be more appropriately used for synergy prediction. The Combination Sensitivity Score clearly stood out with improved performance among the full scope of surveyed approaches and synergy reference models. Furthermore, SYNPRED takes a particular focus on data interpretability, which has been in the spotlight lately when using the most advanced AI techniques.

Control of cell state transitions

Understanding cell state transitions and purposefully controlling them is a longstanding challenge in biology. Here we present cell state transition assessment and regulation (cSTAR), an approach for mapping cell states, modelling transitions between them and predicting targeted interventions to convert cell fate decisions. cSTAR uses omics data as input, classifies cell states, and develops a workflow that transforms the input data into mechanistic models that identify a core signalling network, which controls cell fate transitions by influencing whole-cell networks. By integrating signalling and phenotypic data, cSTAR models how cells manoeuvre in Waddington's landscape1 and make decisions about which cell fate to adopt. Notably, cSTAR devises interventions to control the movement of cells in Waddington's landscape. Testing cSTAR in a cellular model of differentiation and proliferation shows a high correlation between quantitative predictions and experimental data. Applying cSTAR to different types of perturbation and omics datasets, including single-cell data, demonstrates its flexibility and scalability and provides new biological insights. The ability of cSTAR to identify targeted perturbations that interconvert cell fates will enable designer approaches for manipulating cellular development pathways and mechanistically underpinned therapeutic interventions.