A novel approach to predicting the synergy of anti-cancer drug combinations using document-based feature extraction
BMC Bioinformatics volume 23, Article number: 163 (2022)
To reduce drug side effects and enhance their therapeutic effect compared with single drugs, drug combination research, combining two or more drugs, is highly important. Conducting in-vivo and in-vitro experiments on a vast number of drug combinations incurs astronomical time and cost. To reduce the number of combinations, researchers classify whether drug combinations are synergistic through in-silico methods. Since unstructured data, such as biomedical documents, include experimental types, methods, and results, it can be beneficial extracting features from documents to predict anti-cancer drug combination synergy. However, few studies predict anti-cancer drug combination synergy using document-extracted features.
We present a novel approach for anti-cancer drug combination synergy prediction using document-based feature extraction. Our approach is divided into two steps. First, we extracted documents containing validated anti-cancer drug combinations and cell lines. Drug and cell line synonyms in the extracted documents were converted into representative words, and the documents were preprocessed by tokenization, lemmatization, and stopword removal. Second, the drug and cell line features were extracted from the preprocessed documents, and training data were constructed by feature concatenation. A prediction model based on deep and machine learning was created using the training data. The use of our features yielded higher results compared to the majority of published studies.
Using our prediction model, researchers can save time and cost on new anti-cancer drug combination discoveries. Additionally, since our feature extraction method does not require structuring of unstructured data, new data can be immediately applied without any data scalability issues.
Cancer is a global health threat with a high mortality rate; thus, intense research efforts are ongoing to develop new anti-cancer drugs . The development of new anti-cancer drugs follows two approaches: drug discovery, developing previously undiscovered drugs , and drug repositioning, discovering drugs with anti-cancer effects among drugs used for other diseases . In addition to these approaches, the importance of drug combination studies, combining two or more drugs, is increasing .
Two main aims lie behind drug combinations: minimizing drug side effects by reducing drug dosage and achieving a higher therapeutic effect than that of a single drug . In drug combination research, synergy means that the therapeutic effect is increased by combining two or more drugs . For example, the AZD/crizotinib combination has proven effective in both in-vivo and in-vitro cancer treatment experiments . To discover the optimal drug combination, various conditions must be considered such as the side effects between drugs and the dosage of each drug. The number of drug combinations increases exponentially with regards to the total drug amount and dosages per drug.
Drug combinations are discovered in in-vivo and in-vitro experiments. Recently, high-throughput screening has been applied to drug combination discovery, but discovering and filtering potential cases by this method is time-consuming and expensive. To resolve this problem, researchers used in-silico methods, which analyze large amounts of data to discover hidden and meaningful patterns . In-silico methods discover potential drug combinations by extracting them from biomedical data through machine learning.
Using a chemical and genomic information-based normalization strategy, Preuer et al. identified heterogeneity in the input data and created DeepSynergy, a deep learning based drug synergy model using conical layers . Xia et al. presented a neural network-based computational model to encode molecular feature types (gene expression, microRNA, and proteome) and predicted cell line responses to a subset of drug pairs . Kim et al. presented a drug synergy prediction model based on a multitask deep neural network using molecular and genomic drug features, as well as type and genomic cell line features . Zhang et al presented AuDNNsynergy, a deep learning-based drug combination prediction model with encoded multi-omics data (gene expression, mutation, copy number, etc.) using an autoencoder . Janizek et al. presented TreeCombo, an XGBoost-based drug combination prediction model using chemical and physical drug features and cell line gene expression levels . Celebi et al. proposed a model predicting biologically relevant synergistic drug combinations and cell line features using XGBoost . Jeon et al. presented a personalized drug combination prediction model based on extremely randomized trees(ERT) using genomic information, drug targets, and pharmacological information . Li et al. presented a logistic regression-based drug combination prediction model using cell line features (gene expression and essentiality) and drug features (drug target information) .
In most published studies, drug and cell line features were extracted using structured data. Since structured data is mostly created via manual curation, they cannot be easily applied immediately, making data expansion difficult. Unstructured data, such as biomedical documents, include experimental types, methods, and results . By using this information, the unique characteristics of drugs and cell lines can be expressed, and new relationships between drugs and cell lines can be discovered. Therefore, it can be potentially beneficial to predict anti-cancer drug combination synergy using extracted features from documents. However, few studies predict anti-cancer drug combination synergy by extracting features from documents. To solve this problem, we constructed drug and cell line features using biomedical documents. Since our approach uses unstructured data, few restrictions apply with regards to data scalability.
In this study, we propose a novel approach to predicting the synergy of anti-cancer drug combinations using features extracted from documents. The prediction model was created by applying various deep learning and machine learning methods used in published studies, and the results were compared. The model using our features were achieved higher result than majority of published studies. The results show that our approach has promise compared to other feature extraction methods of published studies. Based on the results, researchers can apply our prediction models to predict whether an anti-cancer drug combination has synergy in a specific cell line before in-vivo or in-vitro experiments.
The process is divided into two steps. First, for document preprocessing, we extracted documents containing keywords of validated anti-cancer drug combinations and target cell lines. Then, the extracted documents were preprocessed. Second, in the prediction model creation process, drug and cell line features were extracted using preprocessed documents, and training data were constructed by concatenating the features. Then, a deep and machine learning-based prediction model was created by performing a 5-fold cross validation of the combinations. Figure 1 shows the overall workflow. The details of each step are as follows.
A large matrix of anti-neoplastic agent combinations(ALMANAC) is the most widely used golden standard data in anti-cancer drug combination research, and it was presented by the National Cancer Institute (NCI) . NCI-ALMANAC has been constructed as a drug dose-response matrix by administering more than 5000 FDA-approved drug combinations to 60 human tumor cell lines. Additionally, ComboScore, the reactivity result of each combination, was calculated based on the bliss  score.
The PubMed database is a repository of biomedical documents. We extracted publications from PubMed containing drug and cell line keywords of the NCI-ALMANAC. As PubMed contains a variety of biomedical publications with no specific format, the same drugs and cell lines are differently marked in the documents. For example, in PubMed id “10026447”, the drug Fluorouracil is written as Fluorouracil, whereas in “10026446” it is written as 5-FU. Fluorouracil and 5-FU represent the same drug, but the labels are different depending on the publication. Therefore, the information on the same drug is dispersed, hindering feature extraction. To resolve this problem, we converted synonyms in documents into representative words. When converting synonyms into representative words, PubChem  and Cellosaurus  were referred to for drug and cell line words, respectively. Next, we separated the words through tokenization in documents, performed lemmatization, and removed stopwords.
Prediction model creation
The Word2vec algorithm is a method for embedding a target word using surrounding words . Therefore, all information in the text can be used for embedding features of drugs and cell line. This is an advantage of using unstructured data such as documents when performing feature extraction. In this step, we extracted document-based features by embedding drugs and cell lines from the preprocessed documents. Table 1 shows the hyperparameter settings used to create the Word2vec model.
Since ComboScore is calculated according to drug concentration in NCI-ALMANAC, we used the average combination ComboScores. Table 2 provides an example of NCI-ALMANAC with the average. We separated the dataset into synergistic and non-synergistic combinations using the average combination ComboScores. We extracted the embedding vectors of drug and cell line terms in NCI-ALMANAC from the Word2vec model, and constructed the training data by concatenating drug and cell line features based on synergistic and non-synergistic combinations. A prediction model based on deep learning and machine learning algorithms that can classify various classes was created using the training data. The deep learning-based model used feed-forward neural network (FFNN)  and autoencoder (AE) , and the machine learning-based model used XGBoost (XGB) , extremely randomized trees (ERT) , and logistic regression (LR) . FFNN is a network model in which several layers of perceptrons are sequentially pasted and connected in the direction of the input layer, the hidden layer, and the output layer. AE is a deep neural network model that compresses and reduces input data, then expands it again to make the resulting data identical to the input data. XGBoost is an ensemble algorithm that combines several decision trees and delivers the best performance in tree boosting. ERT is a modified form of the random forest model, which increases randomness by randomly dividing candidate attributes of each forest tree. LR is a model that classifies data into discrete classes for a given feature. We created a prediction model that can classify the synergy of anticancer drug combinations using each learning model. Kunjie et al. constructed a unified framework that can use computational models or architectures presented in published studies . Generally, direct comparison tests between studies are difficult as the data and algorithms used in each study are different. Using their framework, we directly compared each of the prediction models using the same dataset. In addition, since the framework used the same hyperparameters as the published studies, we created the prediction model using the default settings.
NCI-ALMANAC contains 60 cell lines and 101 drugs. Although synonyms for drugs and cell lines exist in NCI-ALMANAC, in order to extract more abundant documents, we extracted synonyms for drugs and cell lines from PubChem and Cellosaurus, respectively. As a result, the number of words increased about 8 times, and a total of 13,962 words were extracted. We extracted 936,734 documents containing these words from PubMed and converted synonyms into their representative words. We extracted the features of drugs and cell lines based on the Word2vec algorithm using preprocessed documents. A total of 130,180 combinations of 60 cell lines and 68 drugs targeting at least one gene were used in NCI-ALMANAC . Cell line features (expression, mutation, etc.) for the baseline were obtained from CellMinerCDB . In addition, drug target and molecular features were obtained from DrugBank  using the RDKit package in Python. Various deep learning and machine learning algorithms were applied using the constructed drug and cell line features as training data.
Table 3 shows the performance model results according to the algorithm. We evaluated the prediction model performances using the unified framework by Kunjie et al. . Their results were used as baseline, and the performance of the prediction model using our features was compared. The prediction model was evaluated five-times each with 5-fold cross validation. We evaluated the model performance using three indicators: area under the receiver operating characteristic curve (ROC-AUC), area under the precision–recall curve (AUPR) and F1 score. Since we repeated five times to evaluate the prediction model, we used the mean and standard deviation as comparative indicators. The prediction model performance constructed using only our features was higher than the majority of the baseline. Figure 2 shows the ROC curve of the FFNN-based prediction model, which is the best performance in ROC-AUC. Based on these results, we can infer that our feature extraction method was able to extract more drug and cell line information than the baseline method.
We introduced an approach to predicting the synergy of anti-cancer drug combination by using features extracted from documents. Our approach yielded higher prediction performance than the majority of published studies under the same conditions. Our approach makes two contributions to drug combination prediction.
As new drugs are discovered, the number of drug combinations increases dramatically. In practice, testing all drug combination cases proves difficult. Our approach can predict new anti-cancer drug combination synergy, reducing the time and effort invested. In other words, we expect that it would be a useful approach for researchers to predict new anti-drug combination synergy before in-vivo or in-vitro experiments.
Structuring the knowledge from newly published studies is time consuming because it is difficult to conduct without domain experts. Using unstructured data to extract features, our approach could be easily and immediately applied to new research. An additional benefit of our feature extraction method is that the extracted features can be utilized together with features from structure data.
Currently, the present model is exclusive to PubMed abstracts. We used PubMed because it is the largest repository of biomedical publication. In future works, the feature extraction method can be expanded to use more information by including full documents provided by PubMed Central, bioRxiv, and other databases. In the present model, we only extracted the drug and cell line features in NCI-ALMANAC. In future works, we will cover more anti-cancer drug combinations, including various large-scale public anti-cancer drug combination databases such as DrugComb  and DREAM . Finally, our approach could also be used for other purposes. For example, by changing the input data, we can use it to predict the interaction between drugs and genes, and also the effects of drugs and diseases.
Availability of data and materials
The result is available at https://github.com/YongsunShim/anticancer-drug-combination-prediction-using-documents-based-features.git.
The National Cancer Institute
A large matrix of anti-neoplastic agent combinations
Feed-forward neural network
Extremely randomized trees
Area under the receiver operating characteristic curve
Area under the Precision–Recall curve
AstraZeneca-Sanger drug combination prediction DREAM challenge
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34.
Tonkens R. An overview of the drug development process. Phys Exec. 2005;31:48–52.
Li J, Zheng S, Chen B, Butte AJ, Swamidass SJ, Lu Z. A survey of current trends in computational drug repositioning. Brief Bioinform. 2016;17:2–12.
Sun W, Sanderson PE, Zheng W. Drug combination therapy increases successful drug repositioning. Drug Discov Today. 2016;21(7):1189–95.
Cheng F, Kovács IA, Barabási A-L. Network-based prediction of drug combinations. Nat Commun. 2019;10:1–11.
Breitinger H-G. Drug synergy-mechanisms and methods of analysis. In: Acree B, editor. Toxicity and drug testing. InTech; 2012. p. 143–66.
Crystal AS, Shaw AT, Sequist LV, Friboulet L, Niederst MJ, Lockerman EL, Frias RL, Gainor JF, Amzallag A, Greninger P, Dana Lee AK, Gomez-Caraballo M, Elamine L, Howe E, Hur W, Lifshits E, Robinson HE, Katayama R, Faber AC, Awad MM, Ramaswamy S, Mino-Kenudson M, Iafrate AJ, Benes CH, Engelman JA. Patient-derived models of acquired resistance can identify effective drug combinations for cancer. Science. 2014;346:1480–6.
Ekins S, Mestres J, Testa B. In silico pharmacology for drug discovery: methods for virtual ligand screening and profiling. Br J Pharmacol. 2007;152:9–20.
Preuer K, Lewis RPI, Hochreiter S, Bender A, Bulusu KC, Klambauer G. Deepsynergy: predicting anti-cancer drug synergy with deep learning. Bioinformatics. 2018;34:1538–46.
Xia F, Shukla M, Brettin T, Garcia-Cardona C, Cohn J, Allen JE, Maslov S, Holbeck SL, Doroshow JH, Evrard YA, Stahlberg EA, Stevens RL. Predicting tumor cell line response to drug pairs with deep learning. BMC Bioinform. 2018;19:71–9.
Kim Y, Zheng S, Tang J, Zheng WJ, Li Z, Jiang X. Anticancer drug synergy prediction in understudied tissues using transfer learning. J Am Med Inform Assoc. 2021;28:42–51.
Zhang T, Zhang L, Payne PRO, Li F. Synergistic drug combination prediction by integrating multiomics data in deep learning models. In: Markowitz J, editor. Translational bioinformatics for therapeutic development. Springer; 2021. p. 223–38.
Janizek JD, Celik S, Lee S-I. Explainable machine learning prediction of synergistic drug combinations for precision cancer medicine. bioRxiv. 2018.
Celebi R, Bear Don’t Walk O, Movva R, Alpsoy S, Dumontier M. In-silico prediction of synergistic anti-cancer drug combinations using multi-omics data. Sci Rep. 2019;9:1–10.
Jeon M, Kim S, Park S, Lee H, Kang J. In silico drug combination discovery for personalized cancer therapy. BMC Syst Biol. 2018;12:59–67.
Li J, Huo Y, Wu X, Liu E, Zeng Z, Tian Z, Fan K, Stover D, Cheng L, Li L. Essentiality and transcriptome-enriched pathway scores predict drug-combination synergy. Biology. 2020;9:278.
Zhu F, Patumcharoenpol P, Zhang C, Chan YYJ, Meechai A, Vongsangnak W, Shen B. Biomedical text mining and its applications in cancer research. J Biomed Inform. 2013;46:200–11.
Holbeck SL, Camalier R, Crowell JA, Govindharajulu JP, Hollingshead M, Anderson LW, Polley E, Rubinstein L, Srivastava A, Wilsker D, Collins JM, Doroshow JH. The national cancer institute almanac: a comprehensive screening resource for the detection of anticancer drug pairs with enhanced therapeutic activity. Can Res. 2017;77:3564–76.
Bliss CI. The toxicity of poisons applied jointly 1. Ann Appl Biol. 1939;26:585–615.
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, Li Q, Shoemaker BA, Thiessen PA, Yu B, Zaslavsky L, Zhang J, Bolton EE. Pubchem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021;49:1388–95.
Bairoch A. The cellosaurus, a cell-line knowledge resource. J Biomol Tech. 2018;29:25–38.
Mikolov T, Chen K, Corrado G, Dean J. Efficient estimation of word representations in vector space. arXiv preprint arXiv:1301.3781 (2013).
Schmidhuber J. Deep learning in neural networks: an overview. Neural Netw. 2015;61:85–117.
Chen T, Guestrin CE. Xgboost: a scalable tree boosting system. In: Krishnapuram B, editor. Proceedings of the 22nd ACM SIGKDD international conference on knowledge discovery and data mining. 2016. p. 785–94.
Geurts P, Ernst D, Wehenkel L. Extremely randomized trees. In: Fürnkranz J, editor. Machine learning. Springer; 2006. p. 3–42.
Dreiseitl S, Ohno-Machado L. Logistic regression and artificial neural network classification models: a methodology review. J Biomed Inform. 2002;35:352–9.
Fan K, Cheng L, Li L. Artificial intelligence and machine learning methods in predicting anti-cancer drug combination effects. Brief Bioinform. 2021;22:1–12.
Luna A, Elloumi F, Varma S, Wang Y, Rajapakse VN, Aladjem MI, Robert J, Sander C, Pommier Y, Reinhold WC. Cellminer cross-database (CellMinerCDB) version 12: Exploration of patient-derived cancer cell line pharmacogenomics. Nucleic Acids Res. 2021;49:1083–93.
Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, Gautam B, Hassanali M. Drugbank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008;36:901–6.
Zagidullin B, Aldahdooh J, Zheng S, Wang W, Wang Y, Saad J, Malyutina A, Jafari M, Tanoli Z, Pessia A, Tang J. Drugcomb: an integrative cancer drug combination data portal. Nucleic Acids Res. 2019;47:43–51.
Menden MP, Wang D, Mason MJ, Szalai B, Bulusu KC, Guan Y, Yu T, Kang J, Jeon M, Wolfinger R, Nguyen T, Zaslavskiy M, Consortium A-SDCD, Jang IS, Ghazoui Z, Ahsen ME, Vogel R, Neto EC, Norman T, Tang EKY, Garnett MJ, Veroli GYD, Fawell S, Stolovitzky G, Guinney J, Dry JR, Saez-Rodriguez J. Community assessment to advance computational prediction of cancer drug combinations in a pharmacogenomic screen. Nature Communications 2019;10:1–17.
The manuscript was proofread by the Dental Research Institute of Seoul National University. The authors want to thank Hyunwhan Joe for language assistance and help with revising the final manuscript.
This work supported by Korea Institute for Advancement Technology(KIAT) grant funded by the Korea government (Ministry of Trade, Industry and Energy) (No.P0015390, AI-based Drug recommendation system for anti-cancer drug combination discovery) and the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (No.2020R1A4A101942311, 4D Ncleome Dynamics in Skeleto-Muscular Aging).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shim, Y., Lee, M., Kim, PJ. et al. A novel approach to predicting the synergy of anti-cancer drug combinations using document-based feature extraction. BMC Bioinformatics 23, 163 (2022). https://doi.org/10.1186/s12859-022-04698-8
- Anti-cancer drug combination
- Drug synergy
- Document-based feature extraction
- Text mining
- Natural language processing
- Deep learning
- Machine learning