A Web Application for Predicting Drug Combination Efficacy Using Monotherapy Data and IDACombo
Article Information
Yunong Xia1,†, Alexander L. Ling1,2,†, Weijie Zhang1, Adam Lee1, Mei-Chi Su1, Robert F. Gruener1, Sampreeti Jena1, Yingbo Huang1, Siddhika Pareek1, Yuting Shan1, and R. Stephanie Huang1,*
1Department of Experimental and Clinical Pharmacology, University of Minnesota, Minneapolis, MN 55455, USA
2Harvey Cushing Neuro-oncology Laboratories, Department of Neurosurgery, Hale Building for Transformative Medicine, 4th and 8th floor, Brigham and Women’s Hospital; 60 Fenwood Road, Boston, MA 02116.
*Corresponding Author: R. Stephanie Huang, Department of Experimental and Clinical Pharmacology, University of Minnesota, Minneapolis, MN 55455, USA.
†These authors should be considered as co-first authors.
Received: 15 November 2023; Accepted: 23 November 2023; Published: 00 November 2023.
Citation:
Yunong Xia, Alexander L. Ling, Weijie Zhang, Adam Lee, Mei-Chi Su, Robert F. Gruener, Sampreeti Jena, Yingbo Huang, Siddhika Pareek, Yuting Shan, and R. Stephanie Huang. A Web Application for Predicting Drug Combination Efficacy Using Monotherapy Data and IDACombo. Journal of Cancer Science and Clinical Therapeutics 7 (2023): 253-258.
View / Download Pdf Share at FacebookAbstract
Summary: 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. Availability and Implementation: The R Shiny app itself can be accessed at https://oncotherapyinformatics.org/idacombo/. The source code for the R Shiny app is available on GitHub (https://github.com/yunong-xia/IDACombo-Shiny-App). The R package IDACombo upon which this app is based is also available on GitHub (https://github.com/Alexander-Ling/IDACombo/).
Keywords
IDACombo, Independent drug action, High-throughput drug screens, Drug combinations, Drug repurposing, Cancer, Computational biology.
IDACombo articles IDACombo Research articles IDACombo review articles IDACombo PubMed articles IDACombo PubMed Central articles IDACombo 2023 articles IDACombo 2024 articles IDACombo Scopus articles IDACombo impact factor journals IDACombo Scopus journals IDACombo PubMed journals IDACombo medical journals IDACombo free journals IDACombo best journals IDACombo top journals IDACombo free medical journals IDACombo famous journals IDACombo Google Scholar indexed journals Independent drug action articles Independent drug action Research articles Independent drug action review articles Independent drug action PubMed articles Independent drug action PubMed Central articles Independent drug action 2023 articles Independent drug action 2024 articles Independent drug action Scopus articles Independent drug action impact factor journals Independent drug action Scopus journals Independent drug action PubMed journals Independent drug action medical journals Independent drug action free journals Independent drug action best journals Independent drug action top journals Independent drug action free medical journals Independent drug action famous journals Independent drug action Google Scholar indexed journals High-throughput drug screens articles High-throughput drug screens Research articles High-throughput drug screens review articles High-throughput drug screens PubMed articles High-throughput drug screens PubMed Central articles High-throughput drug screens 2023 articles High-throughput drug screens 2024 articles High-throughput drug screens Scopus articles High-throughput drug screens impact factor journals High-throughput drug screens Scopus journals High-throughput drug screens PubMed journals High-throughput drug screens medical journals High-throughput drug screens free journals High-throughput drug screens best journals High-throughput drug screens top journals High-throughput drug screens free medical journals High-throughput drug screens famous journals High-throughput drug screens Google Scholar indexed journals Drug combinations articles Drug combinations Research articles Drug combinations review articles Drug combinations PubMed articles Drug combinations PubMed Central articles Drug combinations 2023 articles Drug combinations 2024 articles Drug combinations Scopus articles Drug combinations impact factor journals Drug combinations Scopus journals Drug combinations PubMed journals Drug combinations medical journals Drug combinations free journals Drug combinations best journals Drug combinations top journals Drug combinations free medical journals Drug combinations famous journals Drug combinations Google Scholar indexed journals Drug repurposing articles Drug repurposing Research articles Drug repurposing review articles Drug repurposing PubMed articles Drug repurposing PubMed Central articles Drug repurposing 2023 articles Drug repurposing 2024 articles Drug repurposing Scopus articles Drug repurposing impact factor journals Drug repurposing Scopus journals Drug repurposing PubMed journals Drug repurposing medical journals Drug repurposing free journals Drug repurposing best journals Drug repurposing top journals Drug repurposing free medical journals Drug repurposing famous journals Drug repurposing Google Scholar indexed journals Cancer articles Cancer Research articles Cancer review articles Cancer PubMed articles Cancer PubMed Central articles Cancer 2023 articles Cancer 2024 articles Cancer Scopus articles Cancer impact factor journals Cancer Scopus journals Cancer PubMed journals Cancer medical journals Cancer free journals Cancer best journals Cancer top journals Cancer free medical journals Cancer famous journals Cancer Google Scholar indexed journals Computational biology articles Computational biology Research articles Computational biology review articles Computational biology PubMed articles Computational biology PubMed Central articles Computational biology 2023 articles Computational biology 2024 articles Computational biology Scopus articles Computational biology impact factor journals Computational biology Scopus journals Computational biology PubMed journals Computational biology medical journals Computational biology free journals Computational biology best journals Computational biology top journals Computational biology free medical journals Computational biology famous journals Computational biology Google Scholar indexed journals monotherapy articles monotherapy Research articles monotherapy review articles monotherapy PubMed articles monotherapy PubMed Central articles monotherapy 2023 articles monotherapy 2024 articles monotherapy Scopus articles monotherapy impact factor journals monotherapy Scopus journals monotherapy PubMed journals monotherapy medical journals monotherapy free journals monotherapy best journals monotherapy top journals monotherapy free medical journals monotherapy famous journals monotherapy Google Scholar indexed journals statistics articles statistics Research articles statistics review articles statistics PubMed articles statistics PubMed Central articles statistics 2023 articles statistics 2024 articles statistics Scopus articles statistics impact factor journals statistics Scopus journals statistics PubMed journals statistics medical journals statistics free journals statistics best journals statistics top journals statistics free medical journals statistics famous journals statistics Google Scholar indexed journals cancer drug development articles cancer drug development Research articles cancer drug development review articles cancer drug development PubMed articles cancer drug development PubMed Central articles cancer drug development 2023 articles cancer drug development 2024 articles cancer drug development Scopus articles cancer drug development impact factor journals cancer drug development Scopus journals cancer drug development PubMed journals cancer drug development medical journals cancer drug development free journals cancer drug development best journals cancer drug development top journals cancer drug development free medical journals cancer drug development famous journals cancer drug development Google Scholar indexed journals
Article Details
1. Introduction
Despite the vital role of combination drug therapy in cancer treatment, it is impractical to exhaustively screen the huge number of possible drug combinations experimentally. As such, there are concentrated efforts to develop computational algorithms which can accurately predict the efficacy of drug combinations. To this end, we recently developed IDACombo, a computational method which uses pre- clinical measurements of monotherapy cell response data to predict the efficacy of drug combinations under the assumptions of independent drug action (IDA) [1]. IDA hypothesizes that the baseline efficacy of a combination therapy is simply the effect of the single best drug in the combination, and observed clinical efficacies from a large number of trials spanning many different cancer drug combinations were shown to be consistent with IDA rather than drug additivity or synergy [2, 3]. Given the demonstrated clinical relevance of the predictions produced with IDACombo, we previously released the algorithm as an R package. In this work, we created a web-based app for IDACombo which allows researchers with/without computational background to generate drug combination efficacy predictions using a graphical user interface.
2. Underlying software and hosting
The IDACombo app was created using the shiny v1.5.0 package [4] in R v4.0.3 [5] along with the following packages: shinydashboard v0.7.1 [6], shinyhelper v0.3.2 [7], DT v0.16 [8], tidyverse v1.3.0 [9], IDACombo v1.0.2 [1], shinycssloaders v1.0.0 [10, 11], shinyWidgets v0.5.4 [12], rgl v0.100.54 [13], car v3.0-9 [14], data.table v1.13.6 [15], shinyjs v2.0.0 [11], shinybusy v0.2.2 [16], openxlsx v4.2.3 [17], ggplot2 v3.3.3 [9], gridExtra v2.3 [18], promises v1.1.1 [19], future v1.21.0 [20], doFuture v0.12.2 [20], ipc v0.1.3 [21], and memuse v4.1-0 [22, 23]. The app is hosted on virtual machines (VMs) purchased from DigitalOcean (https://www.digitalocean.com/). Each VM is running Ubuntu 20.04.1 LTS (GNU/Linux 5.4.0- 51-generic x86_64) with 8 virtual CPUs, 16 GB RAM, and 100 GB disk space. Traffic is split across VMs using a load balancer. Connections within each VM are handled via Apache server
2.4 and shiny server v1.5.15.953-amd64.
3. Interpreting IDACombo’s Output Metrics
This app provides predicted drug combination efficacies as both summary plots and downloadable tables. The efficacy metrics used in the plots and tabular outputs are briefly described below, with descriptions of the actual plots themselves being provided later in this manuscript alongside descriptions of the analyses to which each plot is relevant. Note that equations and complete descriptions for these metrics can be found in the original IDACombo manuscript [1].
Average Efficacy within a Population: Efficacy metrics for user provided datasets can be any measure of the effect of a drug on the models being tested. For the pre-provided datasets, viability is used as the efficacy metric, with a viability of 0 indicating all cells died when treated with a compound/combination and a viability of 1 indicating all cells remained alive relative to an untreated control. It should be noted that IDACombo generates predictions of efficacy at the population level rather than the individual level, so reported efficacies are averages over all cell lines/models used in a prediction rather than efficacies for individual cell lines/models.
Hazard Ratios (HRs): HRs represent the relative risk of model systems remaining alive following treatment with a drug combination as compared to a control therapy. With provided datasets (or when custom datasets have been uploaded with the “Lower Efficacy Is Better Drug Effect” option selected), HRs are calculated by simply dividing the mean test treatment efficacy across the selected population of cell lines/models by the mean control treatment efficacy. As such, HRs of 1 indicates that the test therapy provides no benefit relative to the control therapy, whereas lower HRs indicate that the test therapy is more effective than the control therapy. When comparing drug combinations which all include a common control therapy (i.e. comparing all of the possible drugs you could combine with 5-fluorouracil to 5-fluorouracil treatment alone), HRs should be the preferred metric of combination efficacy. It is important, however, to look at HRs relative to both the control therapy and the drug(s) being added to that therapy, as it is possible that adding a drug to a therapy significantly improves efficacy relative to the control therapy but is no better than the efficacy achieved using the added drug as a single agent.
Importantly, when datasets have been uploaded without the “Lower Efficacy Is Better Drug Effect” option selected (i.e. when lower efficacy metric values indicate higher model system survival), HRs are calculated by dividing (1 – mean test treatment efficacy) by (1 – mean control treatment efficacy). As such, IDACombo HRs are only well defined when used with efficacy metric that range from 0 to 1, and care should be taken to avoid situations in which mean treatment efficacies become negative.
IDAComboscores: The IDAComboscore provides a metric by which combinations that do not share a common control therapy can be compared. Higher IDAComboscores indicate better combination efficacy, with the aim of maximizing additional cell death caused by a combination relative to its relevant control therapy while also minimizing the HR of a combination relative to its control therapy.
4. User interface:
4.1 Selecting a Dataset to Generate Predictions With:
Users can use the “Dataset Loader” tab to choose whether to generate drug combination predictions using one of the pre-provided drug screening datasets (i.e. CTRPv2, GDSC1, GDSC2, or PRISM-Repurposing, see Table 1) or using their own dataset. Custom datasets should be formatted according to the instructions provided in the app, and a sample file can be downloaded and edited to ensure correct formatting.
Table 1: Pre-provided drug screening datasets datasets available for use with the IDACombo app
Pre-provided datasets were generated using the Simplicity web app (https://oncotherapyinformatics.org/simplicity/) [24] which provides a simple graphical interface with which users can explore and perform calculations with data from high-throughput drug screens in cancer cell lines. While cell line filtering options in the IDACombo app have been kept to a minimum to avoid unnecessary difficulty for users uploading custom datasets to the website, the Simplicity app provides extensive filtering capabilities for drugs and cell lines and can generate datasets formatted for direct use with the IDACombo app.
4.2 Generating Drug Combination Efficacy Predictions
The IDACombo app provides three basic ways to generate drug combination efficacy predictions: (1) 2-Drug predictions, (2) Control Plus One predictions, and (3) Test vs Control predictions. Each of these approaches are contained in a separate tab within the app, with subtabs providing focused or batch functionalities. Brief descriptions for each approach are provided below.
4.2.1. 2-Drug tab: Allows users to generate predictions for 2-drug combinations at a range of concentrations.
- Focused: Generates predictions for the efficacy of a single 2-drug combination across multiple, user selected, concentrations of each drug. Results are provided as an interactive 3D plot of predicted average efficacy across selected cell lines at each drug-concentration combination (Figure 1A) and as a downloadable tabular output.
- Batch Processing: Generates predictions for many 2-drug combinations across all available concentrations for each drug within the selected cell line population. Results are provided as barplots of the 10 combinations predicted to be most efficacious as ranked by IDAComboscore or HR when combined at the maximum concentrations of each drug (Figure 1B). Results are also provided as a downloadable tabular output.
4.2.2. Control Plus One tab: Allows users to generate predictions of the added benefit from adding a single additional drug to a control treatment consisting of one or more drugs combined at single, user specified, and concentrations for each drug.
- Focused: Generates predictions for adding a single specified compound (at multiple, user selected, concentrations) to an existing control treatment. Results are provided as scatterplots showing how combination efficacy (i.e. viability for pre-provided datasets), HRs (vs. control treatment or vs the drug being added), and IDAComboscores are impacted by combining the drug of interest with the control treatment at different concentrations. Results are also provided as a downloadable tabular output.
- Batch Processing: Generates predictions for adding many drugs (one at a time) to an existing control treatment at all available concentration for the drugs being added. Results are provided as barplots of the 10 combinations predicted to be most efficacious as ranked by IDAComboscore or HR when combined at the maximum concentration of each added drug. Results are also provided as a downloadable tabular output.
- Test vs Control tab: Allows user to generate predictions for whether a test therapy consisting of one or more drugs provides more efficacy than a control therapy consisting of one or more drugs. Results are provided as a downloadable tabular output.
All calculation tabs in the app allow users to specify drug concentrations and cell lines to use when generating predictions. Cell lines can be filtered by general cancer type using a drop- down menu, and check-box options can be selected to modify how calculations are performed. Each option has a clickable question mark next to it which provides details for what each option does.
Figure 1: Automatic visualization of IDACombo predictions. A) 3-d plot of measured of measured monotherapy (blue and green points) and predicted combination (purple points) average viability values for the combination of cisplatin + belinostat across a range of concentrations in 800 cancer cell lines. The grey plane represents the best achievable viability by either monotherapy. Plot was produced using the 2-Drug/Focused tab of the IDACombo app.
B) Bar-plots showing top candidates for 2-drug combination with cisplatin. Plots show top IDAComboscores (top panel) and maximum hazard ratios (bottom panel) when predictions were generated using all available GDSC1 cell lines for cisplatin + any all drugs with available Csustained concentrations in the GDSC1 dataset. Predictions were only generated using concentration ranges between 0 and Csustained for each compound.
5. Primary use cases
5.1 Finding the best single drug to combine with a single drug of interest (use case for batch use of 2-Drug function):
A common use case may be when a researcher has a single drug of interest and would like to identify other drugs that could be efficaciously combined with this drug. This can be achieved by using the batch functionality of the “2-Drug” function in our app. After selecting a dataset in the “Dataset Loader” tab, users should navigate to the “Batch Processing” dropdown option under the “2-Drug” tab. In this page, a user can specify a drug of interest and multiple drugs to add. The app will predict efficacies for all 2-drug combinations between the selected drug of interest and drugs to add using all available concentrations for each drug in the loaded dataset. Note that users must also select which cell lines IDACombo should use to generate the predictions. After selecting these options and pressing “RUN”, a table of the prediction results will be generated which can be navigated within the app or downloaded using the “Download DataTable” button. Batch use of the 2-Drug function can be useful not only as a way for finding an effective 2-drug combination therapy which includes a drug of interest, but also as a means of identifying which drug classes (i.e. mechanisms of action) target different cell populations than a user’s drug of interest. This is because, since IDACombo is built upon IDA, the efficacious combinations identified by IDACombo must consist of combinations of drugs which effectively target at least partially non-overlapping cell populations.
5.2 Finding the best drug to add to a standard combination therapy to improve its efficacy (use case for batch processing in Control Plus One function):
Another common use case may be when a researcher would like to identify drugs which can be added to an existing combination therapy to improve its efficacy. This can be achieved using the “Batch Processing” interface under the “Control Plus One” function tab. To perform the analysis, users define the control therapy by selecting the drugs and drug concentrations of the therapy. Users then select which drugs to add to the therapy and which cell lines to use when predicting combination efficacy. Predictions can be generated by pressing “RUN”, and the app will generate efficacy predictions for each combination of the control therapy + additional drug using all available concentrations for each drug that is being added to the control therapy. As with the 2-Drug function, these predictions will be output as a table which can be navigated within the app or downloaded as a tab-delimited text file. Several plots will also be generated to visualize the drugs that are predicted to combine best with the control therapy.
5.3 Estimating HR for control therapy vs test therapy as might be tested in a clinical trial (use case for non-batch use of Test vs Control function):
The last use case we will highlight is when a researcher may wish to compare a test therapy to a control therapy which may or may not share overlapping drugs with the test therapy. In the page of Test Vs Control, a user can specify a test treatment and a control treatment, as well as cell lines which can also be selected based on cell line subgroups as we previously discussed. The HR computed by IDACombo can indicate whether the test therapy is predicted to have improved efficacy compared to the control therapy. Note that this is not yet validated to be useful when the test therapy does not contain all of the drugs in the control therapy (i.e. when comparing two unrelated therapies such that the test therapy is not the control therapy plus one or more additional drugs).
6. Video tutorials
The following video tutorials have been created to facilitate use of the IDACombo app:
- How the IDACombo algorithm works
- How to use this app
- How to generate custom datasets for this app using the Simplicity app
Links to these videos can be found on the introduction page of the IDACombo web app (https://oncotherapyinformatics.org/idacombo/).
Conclusion
The IDACombo web app allows non-programmatic users to predict drug combination efficacy using monotherapy drug-screening data from cancer cell lines. The predictions of this algorithm have been validated against both pre-clinical and clinical studies of drug combination efficacy (Ling and Huang, 2020). The drug screening data necessary to use IDACombo has been built into the app, and users can upload custom datasets to easily convert their own monotherapy drug screening data into clinically meaningful predictions of drug combination efficacy in cancer.
Funding:
This study was supported by NIH/NCI Grants R01CA204856 (R. S. H). R.S.H. also received support from NIH/NCI R01CA229618 and the University of Minnesota (UMN) OACA Faculty Research Development grant. ALL received funding from NIH T32CA079443.
W.Z. received the UMN BICB first year Fellowship, the UMN IDF Fellowship, and the UMN Clinical & Translational Science Institute (CTSI) A-PReP scholarship.
Author Contributions:
Conceptualization and App Development: YX, ALL, RSH
App Beta Testing: YX, ALL, WZ, AL, MCS, RG, SJ, YH, SP, YS
App Maintenance: YX, ALL, WZ, AL, RSH Manuscript Writing: YX, ALL
Manuscript Review and Editing: YX, ALL, WZ, MCS, RSH
References
- Ling A, and Huang RS. Computationally predicting clinical drug combination efficacy with cancer cell line screens and independent drug action. Nat. Commun 11 (2020): 5848.
- Palmer AC and Sorger PK. Combination Cancer Therapy Can Confer Benefit via Patient-to-Patient Variability without Drug Additivity or Synergy. Cell 171 (2017): 1678-1691.
- Palmer AC. et al. Predictable Clinical Benefits without Evidence of Synergy in Trials of Combination Therapies with Immune-Checkpoint Inhibitors. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res 28 (2022): 368-377.
- Winston Chang et al. shiny: Web Application Framework for R (2020).
- R Core Team R: A language and environment for statistical computing (2020).
- Winston Chang and Barbara Borges Ribeiro shinydashboard: Create Dashboards with ‘Shiny’ (2018).
- Chris Mason-Thom shinyhelper: Easily Add Markdown Help Files to ‘shiny’ App Elements (2019).
- Yihui Xie et al. DT: A Wrapper of the JavaScript Library ‘DataTables’ (2020).
- Wickham H. et al. Welcome to the Tidyverse. J. Open Source Softw 4 (2019): 1686.
- Andras Sali and Dean Attali. shinycssloaders: Add Loading Animations to a ‘shiny’ Output While It’s Recalculating (2020).
- Dean Attali shinyjs: Easily Improve the User Experience of Your Shiny Apps in Seconds (2020).
- Victor Perrier et al. shinyWidgets: Custom Inputs Widgets for Shiny. Wickham,H. (2016) ggplot2: Elegant Graphics for Data Analysis Springer (2020).
- Daniel Adler et al. rgl: 3D Visualization Using OpenGL (2020).
- Fox J and Weisberg S. An R Companion to Applied Regression Third. SAGE Publications, Thousand Oaks, CA (2019).
- Matt Dowle and Arun Srinivasan. data.table: Extension of `data.frame` (2020).
- Fanny Meyer and Victor Perrier. shinybusy: Busy Indicator for ‘Shiny’ Applications (2020).
- Philipp Schauberger and Alexander Walker. openxlsx: Read, Write and Edit xlsx Files (2020).
- Baptiste Auguie. gridExtra: Miscellaneous Functions for ‘Grid’ Graphics (2017).
- Joe Cheng promises: Abstractions for Promise-Based Asynchronous Programming (2020).
- Bengtsson H. (2021) A Unifying Framework for Parallel and Distributed Processing in R using Futures. R J 13 (2021): 208.
- Ian E. Fellows ipc: Tools for Message Passing Between Processes (2019).
- Schmidt,D. Guide to the memuse Package (2020a).
- Schmidt,D. Memuse: Memory Estimation Utilities (2020b).
- Alexander L. Ling et al. Simplicity: web-based visualization and analysis of high- throughput cancer cell line screens. Journal of Cancer Science and Clinical Therapeutics 7 (2023), 249-252.
- Basu A. et al. An Interactive Resource to Identify Cancer Genetic and Lineage Dependencies Targeted by Small Molecules. Cell 154 (2013): 1151-1161.
- Seashore-Ludlow B. et al. Harnessing Connectivity in a Large-Scale Small-Molecule Sensitivity Dataset. Cancer Discov 5 (2015): 1210-1223.
- Rees MG. et al. Correlating chemical sensitivity and basal gene expression reveals mechanism of action. Nat. Chem. Biol 12 (2016): 109-116.
- Iorio F. et al. A Landscape of Pharmacogenomic Interactions in Cancer. Cell 166 (2016): 740-754.
- Yang W. et al. Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res 41 (2013): 955-961.
- Garnett MJ. et al. Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature, 483 (2012): 570-575.
- Corsello SM. et al. Discovering the anticancer potential of non-oncology drugs by systematic viability profiling. Nat. Cancer (2020): 1-14.