oncoPredict: an R package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data.

Cell line drug screening datasets can be utilized for a range of different drug discovery applications from drug biomarker discovery to building translational models of drug response. Previously, we described three separate methodologies to (1) correct for general levels of drug sensitivity to enable drug-specific biomarker discovery, (2) predict clinical drug response in patients and (3) associate these predictions with clinical features to perform in vivo drug biomarker discovery. Here, we unite and update these methodologies into one R package (oncoPredict) to facilitate the development and adoption of these tools. This new OncoPredict R package can be applied to various in vitro and in vivo contexts for drug and biomarker discovery.

Discovering novel pharmacogenomic biomarkers by imputing drug response in cancer patients from large genomics studies.

Obtaining accurate drug response data in large cohorts of cancer patients is very challenging; thus, most cancer pharmacogenomics discovery is conducted in preclinical studies, typically using cell lines and mouse models. However, these platforms suffer from serious limitations, including small sample sizes. Here, we have developed a novel computational method that allows us to impute drug response in very large clinical cancer genomics data sets, such as The Cancer Genome Atlas (TCGA). The approach works by creating statistical models relating gene expression to drug response in large panels of cancer cell lines and applying these models to tumor gene expression data in the clinical data sets (e.g., TCGA). This yields an imputed drug response for every drug in each patient. These imputed drug response data are then associated with somatic genetic variants measured in the clinical cohort, such as copy number changes or mutations in protein coding genes. These analyses recapitulated drug associations for known clinically actionable somatic genetic alterations and identified new predictive biomarkers for existing drugs.

Integration of Pan-Cancer Cell Line and Single-Cell Transcriptomic Profiles Enables Inference of Therapeutic Vulnerabilities in Heterogeneous Tumors.

Single-cell RNA sequencing (scRNA-seq) greatly advanced the understanding of intratumoral heterogeneity by identifying distinct cancer cell subpopulations. However, translating biological differences into treatment strategies is challenging due to a lack of tools to facilitate efficient drug discovery that tackles heterogeneous tumors. Developing such approaches requires accurate prediction of drug response at the single-cell level to offer therapeutic options to specific cell subpopulations. Here, we developed a transparent computational framework (nicknamed scIDUC) to predict therapeutic efficacies on an individual cell basis by integrating single-cell transcriptomic profiles with large, data-rich pan-cancer cell line screening data sets. This method achieved high accuracy in separating cells into their correct cellular drug response statuses. In three distinct prospective tests covering different diseases (rhabdomyosarcoma, pancreatic ductal adenocarcinoma, and castration-resistant prostate cancer), the predicted results using scIDUC were accurate and mirrored biological expectations. In the first two tests, the framework identified drugs for cell subpopulations that were resistant to standard-of-care (SOC) therapies due to intrinsic resistance or tumor microenvironmental effects, and the results showed high consistency with experimental findings from the original studies. In the third test using newly generated SOC therapy-resistant cell lines, scIDUC identified efficacious drugs for the resistant line, and the predictions were validated with in vitro experiments. Together, this study demonstrates the potential of scIDUC to quickly translate scRNA-seq data into drug responses for individual cells, displaying the potential as a tool to improve the treatment of heterogenous tumors. SIGNIFICANCE: A versatile method that infers cell-level drug response in scRNA-seq data facilitates the development of therapeutic strategies to target heterogeneous subpopulations within a tumor and address issues such as treatment failure and resistance.

Computational Modeling to Identify Drugs Targeting Metastatic Castration-Resistant Prostate Cancer Characterized by Heightened Glycolysis.

Metastatic castration-resistant prostate cancer (mCRPC) remains a deadly disease due to a lack of efficacious treatments. The reprogramming of cancer metabolism toward elevated glycolysis is a hallmark of mCRPC. Our goal is to identify therapeutics specifically associated with high glycolysis. Here, we established a computational framework to identify new pharmacological agents for mCRPC with heightened glycolysis activity under a tumor microenvironment, followed by in vitro validation. First, using our established computational tool, OncoPredict, we imputed the likelihood of drug responses to approximately 1900 agents in each mCRPC tumor from two large clinical patient cohorts. We selected drugs with predicted sensitivity highly correlated with glycolysis scores. In total, 77 drugs predicted to be more sensitive in high glycolysis mCRPC tumors were identified. These drugs represent diverse mechanisms of action. Three of the candidates, ivermectin, CNF2024, and P276-00, were selected for subsequent vitro validation based on the highest measured drug responses associated with glycolysis/OXPHOS in pan-cancer cell lines. By decreasing the input glucose level in culture media to mimic the mCRPC tumor microenvironments, we induced a high-glycolysis condition in PC3 cells and validated the projected higher sensitivity of all three drugs under this condition (p < 0.0001 for all drugs). For biomarker discovery, ivermectin and P276-00 were predicted to be more sensitive to mCRPC tumors with low androgen receptor activities and high glycolysis activities (AR(low)Gly(high)). In addition, we integrated a protein-protein interaction network and topological methods to identify biomarkers for these drug candidates. EEF1B2 and CCNA2 were identified as key biomarkers for ivermectin and CNF2024, respectively, through multiple independent biomarker nomination pipelines. In conclusion, this study offers new efficacious therapeutics beyond traditional androgen-deprivation therapies by precisely targeting mCRPC with high glycolysis.

Integration of Computational Pipeline to Streamline Efficacious Drug Nomination and Biomarker Discovery in Glioblastoma.

Glioblastoma multiforme (GBM) is the deadliest, most heterogeneous, and most common brain cancer in adults. Not only is there an urgent need to identify efficacious therapeutics, but there is also a great need to pair these therapeutics with biomarkers that can help tailor treatment to the right patient populations. We built patient drug response models by integrating patient tumor transcriptome data with high-throughput cell line drug screening data as well as Bayesian networks to infer relationships between patient gene expression and drug response. Through these discovery pipelines, we identified agents of interest for GBM to be effective across five independent patient cohorts and in a mouse avatar model: among them are a number of MEK inhibitors (MEKis). We also predicted phosphoglycerate dehydrogenase enzyme (PHGDH) gene expression levels to be causally associated with MEKi efficacy, where knockdown of this gene increased tumor sensitivity to MEKi and overexpression led to MEKi resistance. Overall, our work demonstrated the power of integrating computational approaches. In doing so, we quickly nominated several drugs with varying known mechanisms of action that can efficaciously target GBM. By simultaneously identifying biomarkers with these drugs, we also provide tools to select the right patient populations for subsequent evaluation.

A Web Application for Predicting Drug Combination Efficacy Using Monotherapy Data and IDACombo.

We recently reported a computational method (IDACombo) designed to predict the efficacy of cancer drug combinations using monotherapy response data and the assumptions of independent drug action. Given the strong agreement between IDACombo predictions and measured drug combination efficacy in vitro and in clinical trials, we believe IDACombo can be of immediate use to researchers who are working to develop novel drug combinations. While we previously released our method as an R package, we have now created an R Shiny application to allow researchers without programming experience to easily utilize this method. The app provides a graphical interface which enables users to easily generate efficacy predictions with IDACombo using provided data from several high-throughput cell line screens or using custom, user-provided data.

Inferring therapeutic vulnerability within tumors through integration of pan-cancer cell line and single-cell transcriptomic profiles.

Single-cell RNA sequencing greatly advanced our understanding of intratumoral heterogeneity through identifying tumor subpopulations with distinct biologies. However, translating biological differences into treatment strategies is challenging, as we still lack tools to facilitate efficient drug discovery that tackles heterogeneous tumors. One key component of such approaches tackles accurate prediction of drug response at the single-cell level to offer therapeutic options to specific cell subpopulations. Here, we present a transparent computational framework (nicknamed scIDUC) to predict therapeutic efficacies on an individual-cell basis by integrating single-cell transcriptomic profiles with large, data-rich pan-cancer cell line screening datasets. Our method achieves high accuracy, with predicted sensitivities easily able to separate cells into their true cellular drug resistance status as measured by effect size (Cohen's d > 1.0). More importantly, we examine our method's utility with three distinct prospective tests covering different diseases (rhabdomyosarcoma, pancreatic ductal adenocarcinoma, and castration-resistant prostate cancer), and in each our predicted results are accurate and mirrored biological expectations. In the first two, we identified drugs for cell subpopulations that are resistant to standard-of-care (SOC) therapies due to intrinsic resistance or effects of tumor microenvironments. Our results showed high consistency with experimental findings from the original studies. In the third test, we generated SOC therapy resistant cell lines, used scIDUC to identify efficacious drugs for the resistant line, and validated the predictions with in-vitro experiments. Together, scIDUC quickly translates scRNA-seq data into drug response for individual cells, displaying the potential as a first-line tool for nuanced and heterogeneity-aware drug discovery.

Simplicity: web-based visualization and analysis of high-throughput cancer cell line screens.

High-throughput drug screens are a powerful tool for cancer drug development. However, the results of such screens are often made available only as raw data, which is intractable for researchers without informatic skills, or as highly processed summary statistics, which can lack essential information for translating screening results into clinically meaningful discoveries. To improve the usability of these datasets, we developed Simplicity, a robust and user-friendly web interface for visualizing, exploring, and summarizing raw and processed data from high-throughput drug screens. Importantly, Simplicity allows for easy recalculation of summary statistics at user-defined drug concentrations. This allows Simplicity's outputs to be used with methods that rely on statistics being calculated at clinically relevant doses. Simplicity can be freely accessed at https://oncotherapyinformatics.org/simplicity/.

Simplicity: Web-Based Visualization and Analysis of High-Throughput Cancer Cell Line Screens.

High-throughput drug screens are a powerful tool for cancer drug development. However, the results of such screens are often made available only as raw data, which is intractable for researchers without informatics skills, or as highly processed summary statistics, which can lack essential information for translating screening results into clinically meaningful discoveries. To improve the usability of these datasets, we developed Simplicity, a robust and user-friendly web interface for visualizing, exploring, and summarizing raw and processed data from high- throughput drug screens. Importantly, Simplicity allows for easy recalculation of summary statistics at user-defined drug concentrations. This allows Simplicity's outputs to be used with methods that rely on statistics being calculated at clinically relevant doses. Simplicity can be freely accessed at https://oncotherapyinformatics.org/simplicity/.

Cooperativity between H3.3K27M and PDGFRA poses multiple therapeutic vulnerabilities in human iPSC-derived diffuse midline glioma avatars.

Diffuse midline glioma (DMG) is a leading cause of brain tumor death in children. In addition to hallmark H3.3K27M mutations, significant subsets also harbor alterations of other genes, such as TP53 and PDGFRA. Despite the prevalence of H3.3K27M, the results of clinical trials in DMG have been mixed, possibly due to the lack of models recapitulating its genetic heterogeneity. To address this gap, we developed human iPSC-derived tumor models harboring TP53R248Q with or without heterozygous H3.3K27M and/or PDGFRAD842V overexpression. The combination of H3.3K27M and PDGFRAD842V resulted in more proliferative tumors when gene-edited neural progenitor (NP) cells were implanted into mouse brains compared to NP with either mutation alone. Transcriptomic comparison of tumors and their NP cells of origin identified conserved JAK/STAT pathway activation across genotypes as characteristic of malignant transformation. Conversely, integrated genome-wide epigenomic and transcriptomic analyses, as well as rational pharmacologic inhibition, revealed targetable vulnerabilities unique to the TP53R248Q; H3.3K27M; PDGFRAD842V tumors and related to their aggressive growth phenotype. These include AREG-mediated cell cycle control, altered metabolism, and vulnerability to combination ONC201/trametinib treatment. Taken together, these data suggest that cooperation between H3.3K27M and PDGFRA influences tumor biology, underscoring the need for better molecular stratification in DMG clinical trials.

Facilitating Drug Discovery in Breast Cancer by Virtually Screening Patients Using In Vitro Drug Response Modeling.

(1) Background: Drug imputation methods often aim to translate in vitro drug response to in vivo drug efficacy predictions. While commonly used in retrospective analyses, our aim is to investigate the use of drug prediction methods for the generation of novel drug discovery hypotheses. Triple-negative breast cancer (TNBC) is a severe clinical challenge in need of new therapies. (2) Methods: We used an established machine learning approach to build models of drug response based on cell line transcriptome data, which we then applied to patient tumor data to obtain predicted sensitivity scores for hundreds of drugs in over 1000 breast cancer patients. We then examined the relationships between predicted drug response and patient clinical features. (3) Results: Our analysis recapitulated several suspected vulnerabilities in TNBC and identified a number of compounds-of-interest. AZD-1775, a Wee1 inhibitor, was predicted to have preferential activity in TNBC (p < 2.2 × 10-16) and its efficacy was highly associated with TP53 mutations (p = 1.2 × 10-46). We validated these findings using independent cell line screening data and pathway analysis. Additionally, co-administration of AZD-1775 with standard-of-care paclitaxel was able to inhibit tumor growth (p < 0.05) and increase survival (p < 0.01) in a xenograft mouse model of TNBC. (4) Conclusions: Overall, this study provides a framework to turn any cancer transcriptomic dataset into a dataset for drug discovery. Using this framework, one can quickly generate meaningful drug discovery hypotheses for a cancer population of interest.

More than fishing for a cure: The promises and pitfalls of high throughput cancer cell line screens.

High-throughput screens in cancer cell lines (CCLs) have been used for decades to help researchers identify compounds with the potential to improve the treatment of cancer and, more recently, to identify genomic susceptibilities in cancer via genome-wide shRNA and CRISPR/Cas9 screens. Additionally, rich genomic and transcriptomic data of these CCLs has allowed researchers to pair this screening data with biological features, enabling efforts to identify biomarkers of treatment response and gene dependencies. In this paper, we review the major CCL screening efforts and the large datasets these screens have made available. We also assess the CCL screens collectively and include a resource with harmonized CCL and compound identifiers to facilitate comparisons across screens. The CCLs in these screens were found to represent a wide range of cancer types, with a strong correlation between the representation of a cancer type and its associated mortality. Patient ages and gender distributions of CCLs were generally as expected, with some notable exceptions of female underrepresentation in certain disease types. Also, ethnicity information, while largely incomplete, suggests that African American and Hispanic patients may be severely underrepresented in these screens. Nearly all genes were targeted in the genetic perturbations screens, but the compounds used for the drug screens target less than half of known cancer drivers, likely reflecting known limitations in our drug design capabilities. Finally, we discuss recent developments in the field and the promise they hold for enabling future screens to overcome previous limitations and lead to new breakthroughs in cancer treatment.