 Methodology
 Open access
 Published:
Predicting combinations of drugs by exploiting graph embedding of heterogeneous networks
BMC Bioinformatics volumeÂ 23, ArticleÂ number:Â 34 (2022)
Abstract
Background
Drug combination, offering an insight into the increased therapeutic efficacy and reduced toxicity, plays an essential role in the therapy of many complex diseases. Although significant efforts have been devoted to the identification of drugs, the identification of drug combination is still a challenge. The current algorithms assume that the independence of feature selection and drug prediction procedures, which may result in an undesirable performance.
Results
To address this issue, we develop a novel Semisupervised Heterogeneous Network Embedding algorithm (called SeHNE) to predict the combination patterns of drugs by exploiting the graph embedding. Specifically, the ATC similarity of drugs, drugâ€“target, and proteinâ€“protein interaction networks are integrated to construct the heterogeneous networks. Then, SeHNE jointly learns drug features by exploiting the topological structure of heterogeneous networks and predicting drug combination. One distinct advantage of SeHNE is that features of drugs are extracted under the guidance of classification, which improves the quality of features, thereby enhancing the performance of prediction of drugs. Experimental results demonstrate that the proposed algorithm is more accurate than stateoftheart methods on various data, implying that the joint learning is promising for the identification of drug combination.
Conclusions
The proposed model and algorithm provide an effective strategy for the prediction of combinatorial patterns of drugs, implying that the graphbased drug prediction is promising for the discovery of drugs.
Background
Drug discovery is critical for the therapy of complex diseases, particularly for these scare diseases such as cancers [1]. The current algorithms for drug discovery are roughly classified into two classes, i.e., the biochemical and computational strategies. The traditional strategy makes full use of chemical and biological experiments to synthesize the novel drugs and validate the toxicity of drugs. Generally, this strategy is reliable and effective, whereas it is also criticized for the time and finance. Usually, the design of a drug takes more than 10 years with billions of dollars. Therefore, there is a critical need for alternatives for the traditional method. The traditional assumption that drugs can only be applied to specific diseases, which limits the applications of drugs [2]. Actually, drugs for a specific disease may also be a potential alternative for other diseases because some complex diseases with the same or similar underlying mechanisms, requiring the same or similar therapy strategies and drugs [3].
The accumulated identified drugs and corresponding features provide an opportunity for the identification of potential drugs for diseases without drugs for therapy, which significantly alleviates the burden of period and finance. Furthermore, the predicted drug patterns possibly shed light on the products of drugs. On the basis of the assumption that similar drugs have the same or similar performance on therapy, great efforts have been devoted to this issue with an immediate purpose to identify the potential drugs for diseases by exploiting the similarity among drugs [4]. The major difference between these algorithms lies in how to define and infer the similarity by exploiting various features, such as structure, toxicity, and so on. Given the target disease and potential drugs, the vast majority of current algorithms focus on the ranking of potential drugs by measuring the similarity between the known drugs and candidates, ignoring the effects of drugs combination. Actually, the combination of drugs promotes the performance of treatment. For example, the combination of BRAF and V600E significantly improves the therapy of melanoma, and many drug combinations are approved by FDA (Food and Drug Adminstration) [5].
Therefore, it is promising for the identification of combinatorial drugs for diseases. But, it is highly nontrivial to predict the combination of drugs because combinatorial drugs are very likely to produce synergistic, additive, antagonistic, or even suppressive effects [6]. Consequently, the antagonistic or suppressive drugdrug interactions lead to undesirable consequences. For example, drug DDIs, accounting for 3â€“5% of inpatient medication errors, often lead to patient morbidity and mortality [7,8,9]. Even though it is difficult, great efforts have been devoted to the identification of combination drugs largely due to the merits of applications [10,11,12].
Specifically, the most reliable and intuitive strategy for the identification of drug combinations is biochemical experiments, including the typical P450 testing [13] and transporterassociated interactions [14]. Compared to drug discovery, the identification of drug combinations based on biochemical is much more complicated for several reasons. First, the factors involved in drug combinations, such as the positive and negative interactions among drugs, are much more than the traditional single drug discovery, which imposes a great challenge on the design of experiments and protocols since the balance of multiple factors is complicated. Second, the biochemical experimentbased approaches in clinical trials usually are criticized for the expensive cost, intolerable duration, and unpredictable effects of clinical validation. Finally, the selection of candidates of drug combination is usually impractical because the solution space exponentially increases.
Therefore, many computational algorithms have been devoted to the identification of drug combinations to alleviate the burden of biochemical methods by exploiting the similarity of various drugs [15, 16]. On the basis of the computational strategies, the available algorithms are classified into three classes: statistic, feature and networkbased methods. The most intuitive and straightforward statisticbased methods are the Loewe additivity and Bliss independence models [17], which are two commonly used methods for the quantification of synergy between drug combinations. Specifically, the Loewe additivity model assumes that drugs can be combines if the inhibitors have the same or similar mechanism of regulations, whereas Bliss makes use of independence assumption based on probability theory [18]. The advantage of statisticbased methods is simple and easy to implement. However, these approaches are criticized for two limitations. First, the prerequisite of statistic strategy is the large or superlarge scale samples to guarantee the accuracy of prediction, which hampers the application of these algorithms because some drugs cannot be validated on the huge population. Second, the accuracy of prediction is not desirable because the statistic strategy only focuses on the significance of the difference between groups with various responses, neglecting the features of drugs.
To overcome these problems, the featurebased methods predict the combination of drugs by exploiting the machine learning techniques, such as classification, which aim to extract the most discriminative features. The major difference of featurebased algorithms depends on how to select the features of drugs, and what classifiers to choose for prediction. For example, PDCSGB [19] integrates six types of features to predict drug combinations, including the 2dimensional molecular structures, structural similarity, anatomical therapeutic similarity, proteinâ€“protein interaction, chemicalâ€“chemical interaction, and disease pathways, where three classification algorithms to build the drug combination prediction models are proposed. Compared to the statisticbased methods, PDCSGB not only significantly improves the accuracy of prediction of the drug combination, but also ranks the importance of features. Sun et al. [20] predict drug combinations by integrating the gene expression data of multiple drugs, which enhances the performance of algorithms, indicating that gene expression is also a discriminative feature for drug combinations. To validate the role of genomic features, HNAI [21] fuses the drug phenotypic, therapeutic, structural, and genomic similarities to the prediction of drug combinations by using five machine learningbased classifiers.
Even though the featurebased algorithms dramatically outperform the statisticbased methods, the performance is still unsatisfied because the relations among features are ignored, failing to characterize the indirect relations among features. Fortunately, networks (also called graphs) provide an effective and efficient manner to model and characterize complex systems, where vertices denote entities and edges represent interactions among vertices [22, 23]. Thus, many algorithms have been developed by utilizing the networks of features with an immediate purpose to improve the performance of prediction of drug combinations by exploring the indirect relations of features. The key techniques involved in these algorithms concentrate on the network construction and analysis, where network construction determines how to model the features by using similarity of features, and network analysis focus on how to extract the indirect relations from networks to facilitate the prediction of drug combinations. For example, Liu et al. [24] construct a heterogeneous network by exploiting the similarity between drug and protein, in which three types of similarity, i.e., drugâ€“drug, drugâ€“protein, and proteinâ€“protein, are integrated. Then, they perform random walk to extract features of drugs by exploiting the topological structure of the constructed network, which serves as the input of a gradient tree boosting (GTB) classifier to predict drug combinations. To further explore features, NDD [25] employs a nonlinear fusion method for multiple types of similarity to achieve highlevel features, and then predicts drug combinations by using the neural network. Li et al. [26] integrate multiple types of features to improve the accuracy of prediction by using the neighbor recommender strategy in networks. NIMS [27] made use of networks to screen potential drug combinations, where the diseasespecific biological network is treated as a therapeutic target.
To further improve the performance of prediction, integrating features of various entities is also promising. For example, NLLSS [28] predicts the potential synergistic drug combinations by integrating different kinds of information, such as known synergistic drug combinations, drugtarget interactions, and drug chemical structures, which are integrated into a heterogeneous network. Cheng et al. [29] develop a comprehensive drugdrug interaction network incorporating 6946 interactions of 721 approved drugs using data from DrugBank by using the phenotypic similarity, therapeutic similarity, chemical structure similarity, and gene similarity. EPSDC [30] utilizes the ensemble method to predict the drug combinations by integrating multiplesources information, where construct a feature vector for each pair of drugs by exploiting the drug similarity. Then, the rank of drug pairs is performed by analyzing the topological structure of the heterogeneous drugâ€“target network. Finally, EPSDC fulfills the prediction of drug combination by balancing the rank and output of featurebased classifiers.
Even though significant efforts have been devoted to the prediction of drug combinations, many unsolved problems remain. For example, the current networkbased algorithms are timeconsuming, hindering the applications of largescale networks. Furthermore, the accuracy of current methods can be further improved. Finally, most networkbased algorithms construct the heterogeneous network to characterize the interactions among drug and gene/protein. However, these algorithms extract features of drugs by using the topological analysis strategy for homogeneous networks, ignoring heterogeneity of networks. Recently, graph embedding has been applied to heterogeneous networks, aiming to learn features by preserving the topological structure [31,32,33,34,35,36,37,38,39]. Different from homogeneous networks, it is challenging to develop methods for modeling the heterogeneous types of vertices and edges in a unified way. Heterogeneous information network embedding aims to obtain the lowdimensional representation for each vertex by preserving the topological structure of networks. For example, PME [34] utilizes the metric learning to simultaneously preserve the first and secondorder proximity of heterogeneous networks. Dong et al. [35] design the metapathbased random walks to neighborhoods of vertices and then leverages a skipgram model to perform embedding. metapath2vec preserves both the structures and semantics of a given heterogeneous network by simultaneously learning the lowdimensional and latent embedding for vertices. SHINE [38] extracts the latent representations of vertices by preserving the signs of edges by using autoencoder, and there are some deep learning based methods [39]. In this study, we investigate the possibility of predicting drug combinations by fully exploiting the graph embedding in heterogeneous networks, which is one of the major motivations.
To overcome these problems, we develop a semisupervised heterogeneous network embedding algorithm (called SeHNE) to identify drug combination by integrating features of drugs and proteins, which consists of three major components, i.e., network construction, graph embedding for heterogeneous networks, and prediction of drug combinations. To construct the heterogeneous network for drugs, we integrate the drugâ€“drug, proteinâ€“protein interactions, and drugâ€“target associations. The graph embedding for drugs is performed by nonnegative matrix factorization for the drugdrug and drugâ€“target networks, where the basis matrices are fused to generate the heterogeneous features for drug pairs. To incorporate proteinâ€“protein interaction network into feature extraction, we employ the regularization strategy, where the local topological structure of proteins are preserved. Finally, the feature extraction of heterogeneous and prediction of drug combination are jointly learned. In this case, the features are extracted under the guidance of the classifier, thereby improving the discriminative of features. The experimental results demonstrate that the proposed algorithm outperforms stateoftheart methods in terms of various measurements, such as the area under curve (AUC), average precision (AP), and accuracy.
Results
Overview of SeHNE
The overview of SeHNE is depicted in Fig.Â 1, which consists of three major components, i.e., graph embedding in heterogeneous network, and prediction of drug combination. SeHNE jointly learns features of drugs and prediction of drug combination, where graph embedding for heterogeneous networks is selected to obtain features of drugs, and SVM [40] is used for the prediction of drug combination.
To check whether SeHNE is sensitive to the selection of classifier, six typical classifiers, including knearest neighbor (KNN), logistic regression (LR), random forest (RF) [41], gradient boosting tree (GBT) [12], adaboost (ADB) [42] and linear discriminative analysis (LDA) [43], are selected to replace SVM. The performance of SeHNE with various classifiers in terms of different measurements, including accuracy, precision, recall, Fmeasure, MCC, and AUC, is shown in TableÂ 2, where SeHNE obtains a similar performance, implying SeHNE is not sensitive to classifiers.
To check whether SeHNE is sensitive to different similarities, we replace Anatomical Therapeutic Chemical(ATC) with chemical structure similarity (CSS) (Methods Section). The performance of SeHNE is shown in TableÂ 1, where SeHNE is also stable. These results demonstrate that SeHNE is not sensitive to similarity.
Parameter analysis
There are three parameters involved in SeHNE, where parameter \(\lambda _{1}\) and \(\lambda _{2}\) control the importance of drugtarget network and classifier, and k denotes the number of features in graph embedding. We investigate how AUC of SeHNE changes by varying the value of one parameter with values of other parameters fixed. The drug combination data is split into the training and testing data. We use training data for the 10fold crossvalidation to obtain AUC of SeHNE and then utilize training data to construct the model and testing data to measure the accuracy. How AUC of SeHNE changes as parameter k varying from 20 to 180 by fixing \(\lambda _{1}\) and \(\lambda _2\) as 1 is shown in Fig.Â 2A. As k increases from 20 to 140, the accuracy of SeHNE also improves. However, the performance decreases as k keeps increasing. When k is small, the features in embedding are insufficient to characterize drugs. When k is large, the redundancy of features results in undesirable performance. When k equals 140, SeHNE achieves the best performance.
By fixing k=140, and \(\lambda _{2}\)=1, how AUC of SeHNE changes as parameter \(\lambda _{1}\) increases from 0 to 9 with a gap 1 is shown in Fig.Â 2B, where the AUC decreases as \(\lambda _{1}\) increases. The possible reason is that, when \(\lambda _{1}\) is small, the drugâ€“target and drug networks reach a good balance. When \(\lambda _{1}\) is large, the drugâ€“target network dominates the objective function, where features deviate from the drug network. SeHNE obtains the best performance at \(\lambda _{1}\) = 1. FigureÂ 2C shows the AUC of SeHNE by varying \(\lambda _{2}\) from 0 to 9 by setting k = 140, and \(\lambda _{1}\) = 1. As \(\lambda _{2}\) increases from 0 to 1, the performance improves and then keeps stable. Therefore, in the forthcoming experiments, we set \(\lambda _1=\lambda _2=1\), and \(k=140\).
Performance on drug combination prediction
Before presenting the detailed performance of various algorithms, we give an illustrative example as shown in Fig.Â 3, where panel A is the benchmark interactions, and panel B is the predicted DDI of drugs. In Fig.Â 3, colors correspond to the type of interactions among drugs, where the green edges denote antagonistic, and purple edges represent the synergistic relationship. It is obvious the predicted interactions among drugs are highly consistent with the ground truth ones, i.e., the accuracy is 0.7553 on test data. FigureÂ 3 demonstrates that the proposed algorithm is efficient for the prediction of drug combination (TableÂ 2).
In SeHNE, we adopt SVM as the classifier to predict drug combinations. To select the best kernel function for SVM, we compare SeNMF by using various kernel functions, including the linear, logistic regression, and RBF kernels. The AUC and AP scores of SeHNE by using various kernel functions are shown in Fig.Â 4, where panel A is for AUC, and B for average precision score (AP score). These panels show that the polynomial kernel significantly outperforms the others on AUC and AP score. However, the effect of SeHNE with logistic regression is not ideal because the feature space is large. Thus, we select the polynomial kernel in the experiments.
In this study, the graph embedding and prediction are jointly learned. It is natural to ask whether joint learning is promising for the drug combination. We execute SeHNE in two different strategies, i.e., independent learning and joint learning, where independent learning first extracts the graph embedding and then utilizes SVM to predict drub combination. The results are shown in Fig.Â 4A, B, where joint learning is superior to independent learning in terms of AUC and AP scores. These results demonstrate that joint learning is promising for the drug combination.
Finally, we compare SeHNE with stateoftheart to fully validate the performance of various algorithms. Two algorithms are selected for a comparison, including GTB [24] and EPSDC [30]. These algorithms are selected because EPSDC is the newest method that is simultaneously fusing heterogeneous network and ATC similarity. GTB is the typical algorithm drug combination. The 10fold crossvalidation strategy is used to testify the performance of various algorithms. The AUC and AP score are used to quantify the performance of algorithms. The results are shown in Fig.Â 4, where panel C is for AUC and panel D is for AP score.
The result demonstrates that the SeHNE algorithm outperforms the others in terms of AUC score and AP score. The SVM is better than logistic regression for SeHNE. The reason is that SVM is more discriminative since it exploits the critical features for the prediction of the drug combination. There are three reasons why the proposed algorithm outperforms stateoftheart methods: SeHNE extracts features of drugs by exploiting the heterogeneous network, which is more discriminative than current algorithms because the indirect relations are explored. Joint learning improves the quality of features since features are selected under the guidance of classification. Matrix factorization extracts the latent features from heterogeneous networks, which is more comprehensive to depict drug combinations. These results demonstrate that joint learning of heterogeneous networks and classification is promising for the prediction of drug combinations.
Convergence analysis
The proposed algorithm consists of multiple stages, where the objective function is composed of several components, i.e., \({\mathcal {L}}^{[d]}\) loss, \({\mathcal {L}}^{[t]}\) loss, and \({\mathcal {L}}^{[c]}\) loss. To validate the convergence of the proposed algorithm, we check those components changes as the number of iterations changes as shown in Fig.Â 5, where panel A is for \({\mathcal {L}}^{[d]}\) loss, B for \({\mathcal {L}}^{[t]}\) loss, and C for \({\mathcal {L}}^{[c]}\) loss, respectively. From these panels, it is easy to conclude that these subprocedures quickly converge, i.e., they only take 30 iterations to converge.
Finally, we investigate how the objective function of the proposed algorithm changes as the number of iterations, which is shown in Fig.Â 5 D. SeHNE converges within 30 iterations, implying that the proposed algorithm is efficient. There are two reasons why the proposed algorithm quickly converges. First, SeHNE factorizes the drugâ€“target network by regularizing the drugâ€“drug and proteinâ€“protein interaction (PPI) networks, which enhances the efficiency of feature extraction. Second, the heterogeneous features of proteins and drugs serve as prior information, which accelerates the speed of convergence.
Discussion
Drug discovery is critical for the therapy of complex diseases, particularly for these scare cancers. However, the biological experimentbased methods are time and finance consuming, requiring efficient and effective alternatives for this issue. And, the computational approaches provide an alternative for the traditional biochemical strategy by exploiting features of various entities, such as genes, and proteins. Even though great efforts have been devoted to this issue, vast majority of algorithms solely focus on the identification of potential drugs for complex diseases based on the assumption that similar drugs have similar functions. Therefore, current algorithms concentrate on how to define and compute similarity among drugs with various strategies.
Actually, drug combination is also critical needed since therapy of cancers is complicated, where a single drug is insufficient. However, effort for drug combination is really limited largely because the identification of drug combination is much more complicated than the detection of similar drugs. In this study, we present a novel integrative method for drug combination, where drugs, proteins, and interactions are integrated into a heterogeneous network. The proposed algorithm jointly learns the graph embedding and classification. On the one hand, similar to the previous work [24,25,26,27,28,29,30], SeHNE fuses the drugdrug networks, drugprotein networks, and proteinprotein networks into a heterogeneous network and extracts interesting feature for each drug combination from the heterogeneous network. Furthermore, similar to [24, 25], SeHNE takes features of combined drugs as input, and adopts SVM to predict drug combination. On the other hand, different from these works [34, 36, 38], SeHNE joins the procedures of feature extraction and prediction, where matrix factorization is employ to obtain graph embedding as features of drugs.
SeHNE outperforms baselines in terms of accuracy, implying that the joint learning strategy is more accuracy to model and characterize drug combination. There are two reasons explain why the superiority of the proposed algorithm. First, the topological structure of heterogeneous networks provides complemental information for drugs, thereby improving the quality of features of drugs. Second, graph embedding reflects the latent features of drugs by preserving structural information of drugs.
Conclusion
A novel algorithm for the prediction of combination of drugs is proposed, where multiple types entities are integrated to construct heterogeneous networks. Compared with stateoftheart methods, the proposed algorithm fully makes use of the indirect relations among various entities, which provides a better way to characterize the features of drugs. Furthermore, we present joint learning for feature extraction and prediction of drug combinations, where the features of drugs are more discriminative, resulting in an improved performance. The experimental results demonstrate that the proposed methods outperform the current algorithms in terms of accuracy.
Even though the proposed algorithm algorithm is promising for predicting combination of drugs, there are still some unsolved problems for further study, which are listed as

In this study, the proposed algorithm only focuses on the combination of drug pairs, rather than the highorder combination, because the space of candidates for combinations of drugs exponentially increases. How to narrow the space of feasible drug combination is the foundation for the exploitation of highorder combination of drugs. The strategy for selecting candidates for highorder combination of drugs is critically needed.

The developed algorithm makes use of the topological structure of heterogeneous networks to extract features of drugs, ignoring the intrinsic features of drugs. How to assign attributes, such as the structure and function of drugs to drugs is also promising for modeling and characterization of drugs.

SeHNE only integrates the information of drugs and proteins without considering the regulation principle of drugs. How to integrate gene expression, drug responses and immune microenvironment are also interesting for the identification of drug combinations.
Methods
In this section, we present the model, procedure, and analysis of the proposed algorithm.
Notations
Before presenting the procedure of SeHNE, we present the notations and formulation of drug combination that are widely used in the forthcoming sections.
Given a group of vertices \(\{v_{1},\ldots ,v_{n}\}\) (n is the number of vertices), a network is denoted by \(G=(V,E)\), where \(E=\{(v_{i},v_{j})\}\) is edge set. The adjacent matrix of G is represented by \(A=(a_{ij})_{n\times n}\) where \(a_{ij}\)=1 if vertex \(v_{i}\) and \(v_{j}\) are connected by an edge, 0 otherwise. A network G is heterogeneous if and only if there are more than one type of vertices in G. For example, there are two types of vertices in the heterogeneous network in Fig.Â 6, where the yellow triangle vertices denote drug, and the blue circle ones are proteins. And, there are three types of interactions, i.e., drugâ€“drug interactions, drugâ€“protein interactions and proteinâ€“protein interactions.
Formally, let \(V^{[d]}=\{v_{1}^{[d]},\ldots ,v_{n_d}^{[d]}\}\), and \(V^{[p]}=\{v_{1}^{[p]},\ldots ,v_{n_{p}}^{[p]}\}\) be the drug, and protein set, respectively. And, \(n_{d}\) and \(n_{p}\) are the number of drugs and proteins, respectively. Among proteins, some proteins are targets of drugs, denoted by \(V^{[t]}=\{v_{1}^{[t]},\ldots ,v_{n_{t}}^{[t]}\}\), where \(V^{[t]}\subset V^{[p]}\). There are three types of interactions, including drugâ€“drug, drugâ€“target (proteins), and proteinâ€“protein interaction, denoted by \(E^{[d]}=\{(v_{i}^{[d]},v_{j}^{[d]})\}\), \(E^{[t]}=\{(v_{i}^{[d]},v_{j}^{[t]})\}\), \(E^{[p]}=\{(v_{i}^{[p]},v_{j}^{[p]})\}\), respectively. For a sake of convince, we use \(G^{[d]}=(V^{[d]},E^{[d]})\), \(G^{[t]}=((V^{[d]},V^{[t]}),E^{[t]})\), and \(G^{[p]}=(V^{[p]},E^{[p]})\), respectively. For the drug network \(G^{[d]}\), the similarity matrix \(S^{[d]}=(s_{ij}^{[d]})\in\)[0,1] is constructed, where \(s_{ij}^{[d]}\) represents the similarity between \(v_{i}^{[d]}\) and \(v_{j}^{[d]}\) in terms of ATC. The incidence matrix for \(G^{[t]}\) is denoted by \(B^{[t]}=(b_{ij}^{[t]})_{n_{d}\times n_{t}}\), where \(b_{ij}^{[t]}\)=1 if protein \(v_{j}^{[t]}\) is the target of drug \(v_{i}^{[d]}\), 0 otherwise. The adjacent matrix for \(G^{[p]}\) is constructed as \(A^{[p]}=(a_{ij}^{[p]})_{n_{p}\times n_{p}}\) with element \(a_{ij}^{[p]}\) as 1 if an interaction between protein \(v_{i}^{[p]}\) and \(v_{j}^{[p]}\) exists, 0 otherwise.
Problem definition
Given the heterogeneous network \(G=(G^{[d]},G^{[t]},G^{[p]})\), \(G^{[d]}\) denotes drugâ€“drug similarity graph, \(G^{[t]}\) denotes drugâ€“target (protein) interaction graph, \(G^{[p]}\) denotes proteinâ€“protein interaction graph. Drug combination aims to construct a prediction function \(\phi\) to predict drugâ€“drug interaction between drug \(v_{i}^{[d]}\) and \(v_{j}^{[d]}\). The prediction function \(\phi\) is defined as
where âˆ’ 1 denotes that the two drugs cannot be combined, and this drug combination may produce an antagonistic or even suppressive effect. \(+\) 1 represents that they can combine. This drug combination may produce a synergistic or additive effect. 0 indicates that they are unrelated. In the study, the observed drug combinations is denoted as C, where \(c_{ij}=0\), if drug combination \((v^{[d]}_i,v_j^{[d]})\) is unobserved or unrelated in medical database, \(\pm 1\) otherwise. \(+\) 1 denotes the synergistic or additive effect, âˆ’ 1 denotes the antagonistic or suppressive effect.
Objective function
As shown in Fig.Â 1, the heterogeneous network is categorized into three classes, i.e., the drugâ€“drug network, drugâ€”protein network, and proteinâ€“protein interaction. Given the similarity matrix of drug network \(S^{[d]}\), nonnegative matrix factorization (NMF) [44] is employed to extract the feature of drugs as
where \(W^{[d]}\) and \(H^{[d]}\) are the basis and feature matrix, respectively. EquationÂ (2) is solved by minimizing the approximation, i.e.,
Here, symmetric NMF (SNMF) [45] is performed to extract the features of drugs by exploiting the topology of \(G^{[d]}\) since the similarity matrix is systematic. In this case, Eq.Â (3) is transformed the optimization problem as
where \(W^{'}\) denotes the transpose of W. By minimizing Eq.Â (4), SeHNE generates a low dimensional feature vector for each drug, i.e., graph embedding.
Graph embedding for drugs in Eq.Â (4) is insufficient to fully characterize the features of drugs since the associations between drug and proteins are neglected. To address this issue, it is wise to extract features from the drugâ€“target network. NMF is employed to extract the features of drugs and proteins, i.e.,
where \(W^{[t]}\) and \(H^{[t]}\) are the graph embedding for drugs and proteins under the drugtarget network, respectively. However, features of proteins \(H^{[t]}\) solely reflects the structure of drugâ€”target networks without exploring the information of proteinâ€“protein interactions. The most intuitive strategy is to obtain features of proteins by factorizing \(S^{[p]}\), i.e.,
Then, we can combine Eqs.Â (5) and (6) to obtain graph embedding for proteins as
However, the size of drugâ€“target network is much less than that of proteinâ€“protein interaction networks, where Eq.Â (7) is dominated by the proteinâ€“protein interaction network. To address the problem, we adopt the regularization strategy to integrate proteinâ€“protein interaction network, where graph embedding for proteins \(H^{[t]}\) must preserve the local topological structure in \(G^{[p]}\). Luckily, the Laplacian regularization meets our expectation [46, 47], which is formulated as
where \(L^{[p]}\) is the Laplacian matrix for \(G^{[p]}\). In this case, the features of proteinâ€“protein interaction network is transformed to graph embedding for drugs.
Notice that there are types of graph embedding for drugs either from drugâ€“drug network or drugâ€“target network. We aggregate them as the embedding of drugs. On the classification of drug combination, the loss function for binary classification is employed by mapping drug pairs with âˆ’ 1 or \(+\) 1, which is formulated as
where \(\ell ^{[hl]}(\phi (E_i,E_j),c_{ij}) = max(0,1c_{ij}\phi (E_i,E_j))\) is the hinge loss, \(\phi (E_i,E_j)=\langle \theta ,K\left( E_i,E_j\right) \rangle\) is the inner product of \(\theta\) and K, \(E_{i}\) is the feature vector of drug \(v_{i}^{[d]}\).
Finally, let \({\mathcal {L}}^{[d]}\) and \({\mathcal {L}}^{[t]}\) denote the loss function of two drug embeddings
By combining Eqs.Â (4), (5), (8) and (9) to construct the joint learning framework for drug combination, we formulate the final objective function of the proposed algorithm as
where \(\lambda _1\), \(\lambda _2\) are parameters.
In the next subsection, we derive the updating rule to minimize the objective function in Eq.Â (11).
Optimization rules
EquationÂ (11) is nonconvex, which can not be directly optimized. The iterationbased strategy is adopted, which updates one variable by fixing the others until the algorithm is convergent.
On the optimization of \(W^{[d]}\), we aim to obtain the optimal matrix \(W^{[d]}\) by fixing \(W^{[t]}\) and \(H^{[t]}\). By removing irrelevant terms to \(W^{[d]}\) and employing the alternating minimization algorithms ANLS [48] to solve symmetric problem, the problem expressed in Eq. (11) is transformed into an optimization problem as
The partial deriative on \(W^{[d]}_i\) is derived as
where \(\nabla _{W^{[d]}_i}\ell (\phi (E_i,E_j),c_{ij})\) is the gradient of hinge loss with respect to \(W^{[d]}_i\).
On the optimization of \(W^{[t]},H^{[t]}\), the problem for \(W^{[t]}\) and \(H^{[t]}\) are deduced as
The partial derivatives for \(W^{[t]}\) and \(H^{[t]}\) are calculated as
and
On the optimization of function \(\phi\), we set it as the same strategy for the softmargin linear SVM. The procedure of SeHNE is illustrated in Algorithm 1.
Algorithm analysis
On the space complexity of the SeHNE algorithm, given a heterogeneous network, the space for ATC similarity network \(O(n_d^2)\), the space for Drugtarget network \(O(n_dn_t)\), and the space for PPI network \(O(n_t^2)\). The space for drug embedding and protein embedding are \(O((n_d+n_t)k)\). Therefore, the overall space complexity is \(O(n_d^2+n_t^2+n_dn_t+(n_d+n_t)k)=O((n_d+n_t)^2+(n_d+n_t)k)=O((n_d+n_t)(n_d+n_t+k))\). In our experiment, \(k<n_d<n_t\). Therefore, \(O((n_t+n_t)(n_t+n_t+n_t))=O(n_t^2)\), demonstrating that the proposed method is efficient in space complexity.
Then, the time complexity is analyzed. For each singletype or bipartite networks, SeHNE consists of three major components: symmetric NMF for similarity network, graph regularized NMF for Drugtarget network and PPIs network, classifier learning for drug embeddings. The time complexity of updating \(W^{[d]}\) is \(O(r(n_d^2k+n^2k))\), where n is the number of known drug combinations, r is the number of iteration. The time complexity of updating \(W^{[t]}\), \(H^{[t]}\) are \(O(r(n_t^2k+n^2k))\). The time complexity of updating SVM classifier is based on the number of support vectors, and much more fast than embedding extraction. So the overall time complexity is \(O(r(n_d^2+n_t^2+n^2)k)\). In our experiment, \(n_d<n_t<n\). Therefore, \(O(r(n_d^2+n_t^2+n^2)k)=O(r(n^2+n^2+n^2)k)=O(rn^2k)\).
Data
The training data for drug combination is downloaded from DrugCombDB [49] with the leukemiarelated Cell lines. The drugtarget interactions are derived from DrugBank, KEGG, and Therapeutic Target Database (TTD) datasets, covering 874 drugs and 1240 targets. Among them, we then removed a set of drugs, which donâ€™t have enough information in the drugâ€“drug network or drugâ€“target network. Last, we obtained a group of 370 drugs and 18,126 drugâ€“drug combinations. In general, there are 18,126 drugdrug combinations, consisting of 5903 synergistic and 12,223 antagonistic drug combinations. We summarize the fundamental properties of the drugâ€“drug combinations network (TableÂ 3). The PPI network of humans covers 15,911 proteins with 217,109 interactions [50].
Criteria
A couple of performance measures were used in our experiment, including Accuracy, Precision, Recall, Fmeasure, Matthews correlation coefficient (MCC), and the area under the receiver operating characteristic curve (AUC). They are formally defined as below:
where P, N, TP, FP, TN, and FN are the numbers of real positives, real negatives, true positives, false positives, true negatives, and false negatives, respectively. The AUC is one of the most popular evaluation metrics [51]. AUC is the area under the receiver operating characteristic (ROC) curve, which plots the true positive rate (TPR) versus the false positive rate (FPR).
Based on the therapeutic organ or system of the drug, drugs are classified in the Anatomical Therapeutic Chemical (ATC) coding system. The ATC similarity between \(v_{i}^{[d]}\) and \(v_{j}^{[d]}\) is defined as:
By observing the similarity matrix obtained from the above formula and the definition of ATC similarity, we can get that the ATC similarity matrix is sparse, which is advantageous for application in largescale networks. On the other hand, the use of ATC as drugsimilarity is limited, as the first level is anatomical so that a drug could be found across different body systems.
We also measured the Chemical Structure Similarity(CSS) of drugs to take different information of drugs into account. We get the smiles of each drug from the Drugbank database and then calculate the Molecular ACCess System (MACCS) fingerprints of drug molecules [52] according to its smiles by RDKit (https://github.com/rdkit/rdkit). MACCS is a binary fingerprint (zeros and ones) that answer 166 fragmentrelated questions. If the explicitly defined fragment exists in the structure, the bit in that position is set to 1, and if not, it is set to 0. Therefore, the drug fingerprint is a binary sequence. Then, the Jaccard similarity method is employed to calculate the chemical structure similarity of drugâ€“drug pairs based on molecular fingerprints. Let A and B represent the counts of bits in the twodrug molecules, respectively, the chemical structure similarity between set (drug) A and B is defined as follows:
To utilize the information of two kinds of similarity matrices, two drugdrug similarities are integrated into a comprehensive similarity measure by the probability disjunction formula as
\(S_{ij}^{ATC}\) is the ATC similarity between \(Drug_{i}\) and \(Drug_{j}\) , \(S_{ij}^{CSS}\) is the chemical structure similarity between \(Drug_{i}\) and \(Drug_{j}\).
Availability of data and materials
The data are publicly available in DrugCombDB (http://drugcombdb.denglab.org/), Drugbank (https://go.drugbank.com/), and TTD (http://db.idrblab.net/ttd/).
Abbreviations
 SVM:

Support vector machine
 TCGA:

The Cancer Genome Atlas
 FPKM:

Fragments Per Kilobase of transcript per Million fragments mapped
 COSMIC:

The Catalogue Of Somatic Mutations In Cancer
 NMF:

Nonnegative matrix factorization
References
Atanasov AG, Supuran CT, Zotchev SB, Dirsch VM. Natural products in drug discovery: advances and opportunities. Nat Rev Drug Discov. 2020;20(3):200â€“16.
Spiro Z, Kovacs IA, Csermely P. Drugtherapy networks and the predictions of novel drug targets. J Biol. 2008;7(6):1â€“5.
Keller Benjamin J, Eichinger Felix, Kretzler Matthias. Formal concept analysis of disease similarity. AMIA Summits Transl Sci Proc. 2012;2012:42â€“51.
Moffat JG, Vincent F, Lee JA, Eder J, Prunotto M. Opportunities and challenges in phenotypic drug discovery: an industry perspective. Nat Rev Drug Discov. 2017;16(8):531â€“43.
Robert C, Karaszewska B, Schachter. Improved overall survival in melanoma with combined dabrafenib and trametinib. N Engl J Med. 2015;372(1):30â€“9.
Jia J, Zhu F, Ma X, Cao Z, Li Y, Chen Y. Mechanisms of drug combinations: interaction and network perspectives. Nat Rev Drug Discov. 2009;8(2):111â€“28.
Leape, Lucian L. Systems analysis of adverse drug events. JAMA. 1995;274(1):35â€“43.
Businaro R. Why we need an efficient and careful pharmacovigilance. J Pharmacovigil. 2013.
Karbownik A, SzaEk E, SobaÅ„ska K, Grabowski T, Wolc A, GrzeKowiak E. Pharmacokinetic drugâ€“drug interaction between erlotinib and paracetamol: a potential risk for clinical practice. Eur J Pharm Sci. 2017;102:55â€“62.
Sun X, Vilar S, Tatonetti NP. Highthroughput methods for combinatorial drug discovery. Sci Transl Med. 2013;5(205):205rv1.
Ali MA, Rizvi S, Syed BA. Trends in the market for antihypertensive drugs. Nat Rev Drug Discov. 2017;16(5):309â€“10.
Friedman J. Greedy function approximation: a gradient boosting machine. Ann Stat. 2001;29(5):1189â€“232.
Veith H, Southall N, Huang R. Comprehensive characterization of cytochrome p450 isozyme selectivity across chemical libraries. Nat Biotechnol. 2009;27(11):1050â€“5.
Huang S, Temple R, Throckmorton DC. Drug interaction studies: study design, data analysis, and implications for dosing and labeling. Clin Pharmacol Ther. 2007;81(2):298â€“304.
Chou TC. Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacol Rev. 2006;58(3):621â€“81.
Li X, Qin G, Yang Q, Chen L, Lu X. Biomolecular networkbased synergistic drug combination discovery. BioMed Res Int. 2016;2016:8518945.
Ryall KA, Tan A. Systems biology approaches for advancing the discovery of effective drug combinations. J Cheminformatics. 2015;7(1):1â€“15.
Fitzgerald J, Schoeberl B, Nielsen U, Sorger P. Systems biology and combination therapy in the quest for clinical efficacy. Nat Chem Biol. 2006;2(9):458â€“66.
Qian X, Xiong Y, Dai H, Kumari KM, Wei DQ. PDCSGB: prediction of effective drug combinations using a stochastic gradient boosting algorithm. J Theor Biol. 2017;417:1â€“7.
Sun Y, Yi X, Xu Q, Wei D. A hadoopbased method to predict potential effective drug combination. Biomed Res Int. 2014;2014:196858.
Cheng F, Zhao Z. Machine learningbased prediction of drugdrug interactions by integrating drug phenotypic, therapeutic, chemical, and genomic properties. J Am Med Inform Assoc. 2014;21(e2):e278â€“86.
Peng J, Guan J, Hui W, Shang X. A novel subnetwork representation learning method for uncovering diseaseâ€“disease relationships. Methods. 2021;192:77â€“84.
Peng J, Xue H, Wei Z, Tuncali I, Shang X. Integrating multinetwork topology for gene function prediction using deep neural networks. Brief Bioinform. 2020;22(2):2096â€“105.
Liu H, Zhang W, Nie L, Ding X, Zou L. Predicting effective drug combinations using gradient tree boosting based on features extracted from drugprotein heterogeneous network. BMC Bioinform. 2019;20(1):645.
Rohani N, Eslahchi C. Drugâ€“drug interaction predicting by neural network using integrated similarity. Sci Rep. 2019;9(1):1â€“11.
Li J, Tong X, Zhu L, Zhang H. A machine learning method for drug combination prediction. Front Genet. 2020;11:1000.
Li S, Zhang B, Zhang N. Network target for screening synergistic drug combinations with application to traditional Chinese medicine. BMC Syst Biol. 2011;5:1â€“13.
Chen X, Ren B, Ming C, Wang Q, Yan G. NLLSS: predicting synergistic drug combinations based on semisupervised learning. PLoS Comput Biol. 2016;12(7):e1004975.
Zhang W, Chen Y, Li D, Yue X. Manifold regularized matrix factorization for drugâ€“drug interaction prediction. J Biomed Inform. 2018;88:90â€“7.
Ding P, Yin R, Luo J, Keong KC. Ensemble prediction of synergistic drug combinations incorporating biological, chemical, pharmacological and network knowledge. IEEE J Biomed Health Inform. 2019;23(3):1336â€“45.
Tang J, Qu M, Mei Q. Pte: predictive text embedding through largescale heterogeneous text networks. In: Proceedings of the 21th ACM SIGKDD international conference on knowledge discovery and data mining. 2015. p. 65â€“1174.
Xu L, Wei X, Cao J, Yu PhilipÂ S. Embedding of embedding (eoe): joint embedding for coupled heterogeneous networks. In: Proceedings of the tenth ACM international conference on web search and data mining. 2017. p. 741â€“749.
Chen T, Sun Y. Taskguided and pathaugmented heterogeneous network embedding for author identification. In: Proceedings of the tenth ACM international conference on web search and data mining. 2017. p. 295â€“304.
Chen H, Yin H, Wang W, Wang H, Nguyen H, Li X. PME: projected metric embedding on heterogeneous networks for link prediction. In: Proceedings of the 24th ACM SIGKDD international conference on knowledge discovery & data mining. 2018. p. 1177â€“1186.
Dong Y, Chawla N, Swami A. Metapath2vec: scalable representation learning for heterogeneous networks. In: Proceedings of the 23rd ACM SIGKDD international conference on knowledge discovery and data mining; 2017. p. 135â€“144.
Yue X, Wang Z, Huang J, Parthasarathy S, Moosavinasab S. Graph embedding on biomedical networks: methods, applications, and evaluations. Bioinformatics. 2020;36(4):1241â€“51.
Wang X, Gong X, Yi J, Zhang W. Predicting genedisease associations from the heterogeneous network using graph embedding. In: 2019 IEEE international conference on bioinformatics and biomedicine (BIBM); 2019. p. 504â€“511.
Wang H, Zhang F, Hou M, Xie X, Guo M, Liu Q. Shine: signed heterogeneous information network embedding for sentiment link prediction. In: Proceedings of the eleventh ACM international conference on web search and data mining; 2018. p. 592â€“600.
Chang Q, Han W, Tang J, Qi G, Aggarwal C, Huang T. Heterogeneous network embedding via deep architectures. In: Proceedings of the 21th ACM SIGKDD international conference on knowledge discovery and data mining; 2015. p. 119â€“128.
Platt JohnÂ C. Sequential minimal optimization: a fast algorithm for training support vector machines. Microsoft Research Technical Report (MSRTR9814); 1998.
Breiman L, Breiman Leo, Cutler RaymondÂ A. Random forests machine learning. J Clin Microbiol. 2001;2:199â€“228.
Freund Y, Schapire R. A desiciontheoretic generalization of online learning and an application to boosting. J Comput Syst Sci. 1997;55(1):119â€“39.
Blei D, Ng A, Jordan M. Latent Dirichlet allocation. J Mach Learn Res. 2003;3:993â€“1022.
Lee DD, Seung HS. Learning the parts of objects by nonnegative matrix factorization. Nature. 1999;401(6755):788â€“91.
Kuang D, Ding C, Park H. Symmetric nonnegative matrix factorization for graph clustering. In: Proceedings of the 2012 SIAM international conference on data mining (SDM); 2012. p. 106â€“117.
Cai D, He X, Wu X, Han J. Nonnegative matrix factorization on manifold. In: 2008 eighth IEEE international conference on data mining; 2008. p. 63â€“72.
Belkin M, Niyogi P. Laplacian Eigenmaps and spectral techniques for embedding and clustering. Adv Neural Inf Process Syst. 2001;14(6):585â€“91.
Zhu Z, Li X, Liu K, Li Q. Dropping symmetry for fast symmetric nonnegative matrix factorization. Adv Neural Inf Process Syst. 2018;31:5154â€“64.
Deng L, Zou B, Zhang W, Liu H. DrugCombDB: a comprehensive database of drug combinations toward network medicine and combination therapy. Nucleic Acids Res. 2020;48(D1):D871â€“81.
Cheng F, KovÃ¡cs I, BarabÃ¡si A. Networkbased prediction of drug combinations. Nat Commun. 2019;10(1):1â€“11.
Fawcett Tom. An introduction to ROC analysis. Pattern Recognit Lett. 2006;27(8):861â€“74.
Vilar S, Uriarte E, Santana L, Lorberbaum T, Hripcsak G, Friedman C, Tatonetti NP. Similaritybased modeling in largescale prediction of drugâ€“drug interactions. Nat Protoc. 2014;9(9):2147â€“63.
Acknowledgements
The authors thank the editor and the referees for their valuable comments and suggestions.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 23 Supplement 1, 2022: Selected articles from the Biological Ontologies and Knowledge bases workshop 2020. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume23supplement1.
Funding
This work was supported by National Natural Science Foundation of China (Nos. 11601431, 61772394), Natural Science Foundation of Qinghai Province (No. 2020ZJ920), and Key Research and Development Program of Gansu (Program No. 21YF5GA063). Publication costs are funded by National Natural Science Foundation of China (Nos. 11601431, 61772394). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
SF, ST and XM design the algorithm, SF, XL and ZD process the data and code the software, XL and XM write the manuscript. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
No ethics approval was required for the study.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Song, F., Tan, S., Dou, Z. et al. Predicting combinations of drugs by exploiting graph embedding of heterogeneous networks. BMC Bioinformatics 23 (Suppl 1), 34 (2022). https://doi.org/10.1186/s12859022045674
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12859022045674