Next Article in Journal
Dynamic Variations in Multiple Bioactive Constituents under Salt Stress Provide Insight into Quality Formation of Licorice
Previous Article in Journal
Formulation of a Mixture of Plant Extracts for Attenuating Postprandial Glycemia and Diet-Induced Disorders in Rats
Previous Article in Special Issue
Prediction of Disease-related microRNAs through Integrating Attributes of microRNA Nodes and Multiple Kinds of Connecting Edges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Drug Side-Effect Prediction Via Random Walk on the Signed Heterogeneous Drug Network

1
School of Data and Computer Science, Shandong Women’s University, Jinan 250014, China
2
School of Information Science and Engineering, Shandong Normal University, Jinan 250014, China
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(20), 3668; https://doi.org/10.3390/molecules24203668
Submission received: 13 September 2019 / Revised: 8 October 2019 / Accepted: 10 October 2019 / Published: 11 October 2019
(This article belongs to the Special Issue Molecular Computing and Bioinformatics II)

Abstract

:
Drug side-effects have become a major public health concern as they are the underlying cause of over a million serious injuries and deaths each year. Therefore, it is of critical importance to detect side-effects as early as possible. Existing computational methods mainly utilize the drug chemical profile and the drug biological profile to predict the side-effects of a drug. In the utilized drug biological profile information, they only focus on drug–target interactions and neglect the modes of action of drugs on target proteins. In this paper, we develop a new method for predicting potential side-effects of drugs based on more comprehensive drug information in which the modes of action of drugs on target proteins are integrated. Drug information of multiple types is modeled as a signed heterogeneous information network. We propose a signed heterogeneous information network embedding framework for learning drug embeddings and predicting side-effects of drugs. We use two bias random walk procedures to obtain drug sequences and train a Skip-gram model to learn drug embeddings. We experimentally demonstrate the performance of the proposed method by comparison with state-of-the-art methods. Furthermore, the results of a case study support our hypothesis that modes of action of drugs on target proteins are meaningful in side-effect prediction.

1. Introduction

Drug side-effects, or adverse drug reactions (ADRs), can be regarded as undesirable effects that are caused by normal use of drugs. Most are natural pharmacological actions of the drugs and are unavoidable. Side-effects of approved drugs are harmful to patients and can even be fatal. Serious side-effects cause 100,000 deaths per year in the United States and have been one of the leading causes of death [1]. The detection of side-effects is challenging in every stage of drug development.
Traditional side-effect detection methods such as vitro safety profiling and clinical drug safety trials are time-consuming and expensive, and many potential side-effects cannot be detected because there are so many side-effect terms. Many computational methods for analyzing or predicting drug side-effects have been proposed recently [2,3,4,5]. Li et al. [6] combined multiple data sources and proposed an inductive matrix completion method for predicting unknown side-effects. Zhang et al. [7] proposed an integrative label propagation algorithm to predict potential side-effects based on high-order similarity. Zheng et al. [8] built a drug similarity integration framework to measure the similarity between drugs from various perspectives and proposed a highly reliable negative sample selection method to improve the prediction accuracy. Zhang et al. [9] adopted ensemble methods and the feature-selection-based multilabel k-nearest neighbor method (FS-MLKNN) for side-effect prediction. Liu et al. [10] built a series of binary classifiers to determine whether a drug has a specified side-effect. With the advancement of graph embedding technologies [11,12,13], many studies have focused on learning integrated drug embeddings for various prediction tasks. Ma et al. [14] proposed a drug embedding method that is based on multi-view deep autoencoders to predict drug side-effects. Hu et al. [15] proposed a heterogeneous network embedding approach by integrating PPI information into drug embeddings.
The methods that are discussed above mostly utilize chemical profiles (e.g., fingerprints) and biological profiles (e.g., target proteins, pathways, and transporters) of drugs to predict potential side-effects. It has been proved that drug–target interactions play an important role in side-effect prediction [16,17,18]. However, the drug–target interactions that are considered in these works focus only on whether a single drug acts on a target protein and neglect the modes of action of drugs on target proteins. There are various and even opposite action modes when drugs act on target proteins, such as activation, inhibition, agonist, antagonist, potentiator, blocker, inducer, suppressor, and so on. The literature [19] found that modes of action of drugs are strongly related to the therapeutic effects and side-effects of drugs through statistical analysis.
Based on the results of [19], in this paper, we propose utilizing the action mode information of drug–target interactions to learn more effective drug embeddings for drug side-effect prediction. The modes of action of drugs on target proteins can be represented as positive or negative edges in a signed graph (as shown in Table 1). When two drugs act coherently on a common target protein, the two drugs are regarded as similar and are connected with a positive edge. In contrast, when two drugs act incoherently on a common target protein, the two drugs are regarded as dissimilar and are connected with a negative edge. This produces a signed drug network. Signed networks have been frequently used in systems biology [20,21,22,23] and can reveal the deeper complex relations between individuals. The signed drug information and other traditional drug information form a signed heterogeneous information network (signed HIN). A signed HIN is a network with multiple types of nodes and links, and the links can be positive or negative [24]. Multiple types of drug information and every kind of information can be formulated as a signed subnetwork or an unsigned subnetwork. The subnetworks form a signed drug HIN. We projected every drug into a low dimensional vector space and predicted the side-effect of the drug. The low dimension vector representations of drugs synthesized different drug information. The problem of learning drug representations can be formulated as a network embedding task on the signed drug HIN. We propose utilizing random walk [25] on the signed HIN to obtain the drug sequences and utilizing the Skip-gram model [26] to learn drug embeddings. We also experimentally demonstrate that the learned drug embeddings that integrating additional signed drug information can substantially improve the side-effect prediction accuracy.
A flowchart of our side-effect prediction model via random walk on a signed HIN (RW-SHIN) is shown in Figure 1. First, we collected drug profiles from public databases, especially action modes of drug–target interactions, from the DrugBank database [27]. Then, we constructed a signed heterogeneous information network (signed HIN) based on these profiles. The signed HIN contained one signed drug subnetwork and three unsigned drug subnetworks. RW-SHIN performed different random walk procedures on signed subnetworks versus unsigned subnetworks. Here, two biased random walk procedures [28,29], which can effectively incorporate positive and negative relations into a graph, were used to obtain nodes sequences. Then, we utilized the Skip-gram model to learn drug embeddings based on the node sequences. After that, the drug embeddings were input into a fully connected neural network to predict a potential side-effect.

2. Materials

2.1. Drug Chemical Profile

We collected the fingerprints of drugs from the PubChem Compound database [30]. Each drug was encoded into an 881-dimensional feature vector, with entries of 0 or 1 representing the presence or absence, respectively, of a chemical substructure term.

2.2. Drug Side-Effect Profile

SIDER [31] is a database that contains the side-effects of drugs on the market that are listed in their package inserts. We downloaded the entire database from http://sideeffects.embl.de/. The information in SIDER is limited because many side-effects may not be detected in clinical trials. OFFSIDES [32] is another side-effect database, which was built by mining the FDA Adverse Event Reporting System (FAERS, http://www.fda.gov/cder/aers/default.htm). We extracted the associations between drugs and side-effects from the two databases. This produced a dataset that consisted of 548 drugs, 1385 side-effect terms, and 41,008 associations between drugs and side-effects, and each drug had 74.8 side-effects on average. An index plot of the number of associated drugs for each side-effect is shown in Figure 2.
Each drug was represented by a 1385-dimensional feature vector, where each element is 1 or 0, which represents the presence or absence, respectively, of a side-effect term for this drug.

2.3. Drug Biological Profile

We crawled the biological profiles of drugs in the DrugBank database [27], which is a widely used public drug information database. We extracted the target proteins of the 548 drugs without considering the action modes of the drugs on the target proteins. Each drug is encoded into a 780-dimensional feature vector, with 0 or 1 representing the presence or absence, respectively, of a target protein for this drug.
Meanwhile, we crawled the action modes of 548 drugs on 780 target proteins and associated every mode with a sign according the procedure that is described in [15] (as presented in Table 1). This produced a signed drug–target association matrix A 1 , 0 , 1 m n , where m is the number of drugs and n is the number of targets.

2.4. Construction of the Signed Drug Heterogeneous Information Network

We use the above three profiles to construct a drug heterogeneous information network (HIN) and learn the drug embeddings based on the constructed drug HIN. A heterogeneous information network is a special type of information network that contains several types of nodes or links. In this work, the drug HIN contains three unsigned subnetworks and one signed drug subnetwork.

2.4.1. Unsigned Drug Subnetworks

We constructed three unsigned drug subnetworks based on drug chemical profiles, drug side-effect profiles and unsigned drug–target associations in the drug biological profiles. The nodes in the three subnetworks correspond to all 548 drugs. The edges are defined by the Jaccard similarities of the feature vectors in every type of profile. The edge weight w i j u between drug pair ( d i , d j ) in an unsigned subnetwork is defined by Equation (1).
w i j u = P r i P r j P r i P r j
where P r i is the feature vector of drug i in a type of drug profile.

2.4.2. Signed Drug Subnetwork

We constructed a signed drug subnetwork G s using the signed drug–target association matrix A 1 , 0 , 1 m n . The node set in G s is the drug set and the edge weight w i j s between drug pair ( d i , d j ) is defined by Equation (2).
w i j s = k = 1 n a i k · a j k k = 1 n a i k · a j k
where a i k is the element of the signed drug–target association matrix A. w i j s ranges from 1 to 1, where w i j s = 1 if drug pair ( d i , d j ) share several targets but have opposite action modes on the shared targets, while w i j s = 1 if drug pair ( d i , d j ) share several targets and have the same action mode on the shared targets and w i j s = 0 if drug pair ( d i , d j ) have no shared targets. The distribution of the edge weights w i j s is plotted in Figure 3.
According to Figure 3, most drug pairs have strong connections for both positive and negative edges. Hence, the drug pairs have only several targets, but the actions on the shared targets are extremely similar or different. The side-effects may differ if the actions on the shared targets differ. The signed drug subnetwork is illustrated in Figure 4.

3. Methods

In our model, side-effects were predicted using the low dimensional vector representations of drugs. The drug representations are learned from the above-mentioned drug profiles. Learning comprehensive drug representations can be formulated as the task of learning node embeddings of the constructed drug HIN. Therefore, we first performed a random walk procedure on every subnetwork in the above-constructed drug HIN to obtain nodes sequences, and then utilized the Skip-gram model to learn drug embeddings based on the obtained node sequences. After that, potential side-effects were predicted using the learned drug embeddings.

3.1. Random Walk on the Signed HIN

The drug HIN that is presented above contains both unsigned subnetworks and signed subnetworks. In the signed subnetworks, the learned node representations must integrate the positive and negative relations. An effective approach is to conduct random walks on the graph to integrate the positive and negative relations into node sequences. In this paper, we use two types of biased random walk procedures [28,29] to obtain many node sequences in each subnetwork. The two biased random walk procedures are conducted on signed graphs versus unsigned graphs.

3.1.1. Biased Random Walk on a Signed Graph

The biased random walk procedure on a signed graph utilizes the direct neighbor relations and common neighbor relations to acquire a node sequence. As illustrated in Figure 5, for a random walker, let t i 1 = s denote the i 1 t h node in the path and let node t i = h be the next-hop node.
If the edge between node s and node h is positive, the walker chooses h as the next hop according to the direct neighbor relation ( s b + ) with transition probability w s h s i n b s ( s ) w s i s .
If the edge between node s and node k is negative, the walker will not choose k as the next hop but may choose k s enemy node j as the next hop, which is guided by the common neighbor relation. The common neighbor relation involves two types of sets: friend sets F ( w ) and enemy sets E ( w ) (as illustrated in Figure 6). If two nodes’ friend sets or enemy sets have a large overlap area, the two nodes share common friends or enemies, respectively, and the relation between them tends to be positive. In contrast, if the enemy set of a node covers a large overlap area of the friend set of another node, the two nodes have very little in common and the relation between them will be weak or even negative.
Here, we define a common neighbor relation similarity l s ( s , h ) , which is used to decide whether the walker jumps to an enemy’s enemy node.
l s ( s , h ) = r F ( s ) F ( h ) w s r s + n E ( s ) F ( h ) w s n s m E ( s ) F ( h ) w s m s p F ( s ) E ( h ) w s p s
If l s ( s , h ) > 0 , h is the enemy’s enemy of node s. The walker will jump to h from s according to the transition probability w s k s · w k h s i n b s ( s ) w s i s · j n e g a t i v e n b s ( k ) w k j s .
The transition probability in the biased random walk procedure is defined in Equation (4).
p s h = p · w s h s i n b s ( s ) w s i s w s h s > 0 q · w s k s · w k h s i n b s ( s ) w s i s · j n e g a t i v e n b s ( k ) w k j s w s k s < 0 a n d w k h s < 0 a n d w s h s = 0 a n d l s ( s , h ) > 0 0 o t h e r w i s e
where p and q are adjustable parameters that guide the walker by selecting one of its friends or an enemy of its enemy as the next hop.
The biased random walk procedure assumes that an enemy’s enemy may be a friend. We investigated the triangle relationships in the signed drug subnetwork to determine whether this assumption holds in our datasets. The results are presented in Table 2.
In the above investigation, we find that an enemy’s enemy is more likely to be a friend in a signed drug subnetwork. Therefore, we can use the biased random walk procedure to combine the positive and negative relations and to obtain node sequences.

3.1.2. Biased Random Walk Procedure on an Unsigned Graph

The biased random walk procedure on an unsigned graph combines two sampling strategies: breadth-first sampling (BFS) and depth-first sampling (DFS). For a random walker, let t i 1 = s denote the ( i 1 )th node and let node t i 2 = t denote the ( i 2 )th node in the path. Node t i = h will be the next-hop node. As illustrated in Figure 7, h 0 or h 1 may be the next-hop node according to BFS and DFS, respectively.
The transition probability in the biased random walk procedure is defined in Equation (5).
p s h = α · w s h u i n b s ( s ) w s i u w t h u 0 β · w s h u i n b s ( s ) w s i u w t h u = 0
where α and β are adjustable parameters that guide the walker selecting the next node following BFS or DFS. t is the last node in the path.

3.2. Learning Drug Embeddings

After obtaining node sequences in the signed HIN, we utilized Skip-gram model [26] to learn node representations. The model of Skip-gram is illustrated in Figure 8 and the optimizer objective function of Skip-gram is presented as Equation (6).
Γ = min f log Pr ( { v i w , , v i 1 , v i + 1 , , v i + w } f ( v i ) )
where f is the node mapping function f : v V R V × d , in which d is the dimension of the embeddings; v i w is the w t h neighbor node of node v i in every node sequence; and w is the window size.

3.3. Prediction Formulation

Once we have obtained the final drug embeddings, we use a fully connected neural network to predict each side-effect, as illustrated the last step in Figure 1. The predictive problem was modeled as a binary classification of each side-effect. For each side-effect, the drugs that are known to cause the side-effect were labeled as positive samples and the remaining drugs were labeled as negative samples. The predictive model utilized the L 2 -norm-regularized logistic regression method [33] to predict each side-effect. The input features of the model were drug embeddings. The loss function is presented as Equation (7).
L ( Z , y , w ) = i C o s t ( y i , f ( z i w + b ) ) + λ w 2 2
where Z is the feature matrix of the drugs and the row vector z i of Z is the embedding of drug i; f is a nonlinear logistic function with L 2 loss; w and b are the weight vector and the bias vector, respectively; and λ is a hyperparameter, which must be learnt from the data.

4. Results and Discussions

We experimentally compared the performance of the proposed method with several state-of-the-art methods and utilized a case study to support our hypothesis that modes of action of drugs on target proteins are meaningful in side-effect prediction. We also discussed the impact of embedding dimension on the prediction performance to find the best representations of drugs.

4.1. Performance Evaluation Metrics

We used the receiver operating characteristic curve ( R O C curve) to evaluate the performance of each method. The R O C curve is the true-positive rate ( T P R ) as a function of the false-positive rate ( F P R ), which is based on various thresholds, where T P R and F P R are defined in Equations (8) and (9).
T P R = T P T P + F N
F P R = F P F P + T N
where T P , F P , T N , and F N are the numbers of true positives, false positives, true negatives, and false negatives, respectively.
We calculated the area under the R O C curve ( A U R O C ) to evaluate the performance of each method in predicting all side-effects.

4.2. Baselines

We compared our method (RW-SHIN) with several state-of-the-art network embedding algorithms for side-effect prediction. We implemented the following four baselines for comparison:
  • Laplacian eigenmaps [34]: Laplacian eigenmaps is a typical matrix factorization method that has been widely adopted for data analysis of biomedical networks. It aims at factorizing a data matrix of a graph into lower dimensional matrices while preserving the topological properties of the original graph. For the drug HIN, we concatenated the Laplacian eigenmaps of each unsigned subnetwork to construct feature vectors of drugs for side-effect prediction.
  • GCN [35]: GCN is a recently proposed network embedding method that is based on the spectral convolutional operation and realizes state-of-the-art performance on important prediction problems in recommender systems. Here, we linearly integrated the similarity matrices of the unsigned subnetworks in the drug HIN and learned the drug embeddings using GCN.
  • AttSemiGAE [14]: AttSemiGAE utilizes multiview graph autoencoders (GAEs) and adds an attentive mechanism for determining the weights for each view with respect to the corresponding prediction tasks. Here, each unsigned subnetwork in the drug HIN is regarded as a single-view graph in the AttSemiGAE algorithm and the supervised information is the known side-effects of the drugs.
  • RW-HIN: For further validation of the impacts of action modes on the quality of side-effect predictions, we designed a network embedding algorithm that ignores the signed drug information, namely, RW-HIN. The algorithm is based on random walk on an unsigned graph.
These algorithms cannot learn node embeddings on signed graphs and only utilize the following data sources to learn drug embeddings: drug chemical profiles, drug side-effect profiles and unsigned drug biological profiles. We performed the following 5-fold cross-validation procedure for each algorithm: The known drug–side-effect associations are used as a gold-standard set. The drugs in the gold-standard set are split into five equally sized subsets, and each subset is used in turn as the test set, while the remaining four subsets are used as the training set. The final performance is evaluated according to the results on all 5 folds. For accurate comparison, we use the same experimental conditions, namely, the same training drugs and test drugs are used across all methods in each cross-validation fold. The embedding dimension of each algorithm is 32. All parameters in each method (e.g., the regularization parameters, the adjustable parameters, and the embedding dimension) were optimized using grid search with the AUROC score as an objective function.

4.3. Result of Comparison

Figure 9 plots the ROC curves for the five approaches based on the cross-validation experiment. The ROC curve of each method is based on the merged prediction scores of all side-effects. Our method outperforms the other methods, which suggests that the proposed method is more effective than the previous methods. Although RW-SHIN and RW-HIN are based on the same training model, RW-SHIN, which considers the modes of action of drugs on target proteins, outperforms RW-HIN. This suggests that the modes of action of drugs on target proteins are meaningful. In addition, the method with an attentive mechanism outperforms the methods that do not use attentive mechanisms on the same unsigned datasets.
We examined the prediction accuracies for individual side-effects. We calculated the AUROC scores for each side-effect. Figure 10 presents a boxplot that represents the distribution of the resulting AUROC scores for 1385 side-effects for each method. Our signed drug–target information-based method produced the best results. Most high-frequency side-effects (e.g., pseudomembranous colitis, gynecomastia and interstitial nephritis) were predicted with higher accuracy, while most low-frequency side-effects (e.g., diabetic neuropathy, nephrogenic diabetes insipidus and narcolepsy) were predicted with lower accuracy. The ROC-AUC scores are positively correlated with the side-effect frequencies. However, this is not absolute: some high-frequency side-effects (e.g., abdominal pain) were predicted with lower accuracy, while some low-frequency side-effects (e.g., chronic active hepatitis) were predicted with higher accuracy.
We also investigated the side-effect prediction accuracy for every drug by evaluating the T o p K prediction accuracy of each drug. We counted the number of drugs whose T o p 1 , T o p 3 , and T o p 5 prediction results contain at least one known side-effect. If at least one known side-effect is in the T o p K prediction results of a drug, the drug will be counted. The results are presented in Table 3. The proposed signed drug–target association-based method outperforms other methods. The number of drugs whose T o p 1 prediction side-effects are known is 276 (50.3%). The numbers of drugs whose T o p 3 and T o p 5 prediction results contain at least one known side-effect are 354 (64.6%) and 465 (84.9%), respectively.

4.4. Case Studies

We investigated the side-effect prediction results of drugs that are similar in terms of chemical profiles and target proteins but opposite in terms of modes of action on target proteins. We consider Betaxolol (DB00195) and Dobutamine (DB00841) as examples. Betaxolol is a selective beta-1 adrenergic receptor blocker that is used in the treatment of hypertension and glaucoma. Dobutamine is a sympathomimetic drug that is used in the treatment of heart failure and cardiogenic shock. Its primary mechanism is direct stimulation of beta-1 receptors of the sympathetic nervous system.
The chemical similarity of drugs can be calculated using Jaccard similarity, which is shown in Equation (1). The chemical similarity of the two drugs is 0.82. The two drugs have two common target proteins: Beta-1 adrenergic receptor (P08588) and Beta-2 adrenergic receptor (P07550). The unsigned drug biological similarity of the two drugs is 0.67. However, the action modes of the two drugs are opposite. Dobutamine acts as an agonist on the two target proteins, while Betaxolol acts as an antagonist. The top-10 prediction results of the two drugs are presented in Table 4.
The prediction results of the two drugs differed substantially and there was only one common side-effect: nausea. Nausea is a very common side-effect and may be caused by most drugs. This difference is consistent with the known side-effect similarity of the two drugs, which is 0.083. Focusing on the false prediction results of Dobutamine, we found that insomnia is a known side-effect of Pindolol (DB00960), which is an oral beta blocker that is used to treat hypertension. Dobutamine and Pindolol have the same action modes on the common target proteins: the Beta-1 and Beta-2 adrenergic receptors. The chemical similarity of Dobutamine and Pindolol is 0.873. Dobutamine [36] belongs to a class of beta agonists that typically have mild to moderate adverse effects, which include increased heart rate and insomnia.

4.5. Performance Comparison among Embedding Dimensions

To examine the impact of the embedding size on the prediction performance, we compared our two models with various dimensions of drug embeddings in terms of A U R O C . The results are presented in Figure 11. With the increase of the embedding dimension, the side-effect prediction performance initially increased and subsequently decreased. If the embedding dimension is small, the embeddings may lose too much information. If the embedding dimension is very large, the embeddings may incorporate redundant information and tend to overfit. The two methods performed optimally when the embedding dimension was 32. RW-SHIN outperformed RW-HIN at the same embedding dimensions. This supports our assumption that the modes of action of drugs on target proteins can provide useful information for drug embeddings.

5. Conclusions and Future Work

In this work, we proposed an improved drug side-effect prediction method by integrating modes of action of drugs on target proteins into drug information. We proposed a side-effect prediction method that learns more comprehensive drug embeddings based on drug chemical profiles, unsigned drug–target protein associations, signed drug–target protein associations, and drug side-effect profiles. To the best of our knowledge, no previous work considers modes of action of drugs on target proteins in the context of drug side-effect prediction. To formulate this feature of drug information, identical or opposite action modes were defined as positive or negative edges in a signed graph and formed a signed heterogeneous information network with other drug information. We used two bias random walk procedures to obtain drug sequences in the signed HIN and utilized the Skip-gram framework to learn drug embeddings based on the obtained drug sequences. In the cross-validation experiments, all results demonstrate that the prediction performance improved substantially using our method. According to the case study on the drugs Betaxolol and Dobutamine, our method can predict both existing and novel drug–side-effect associations. Moreover, the results of the case study support our hypothesis that action modes are meaningful in side-effect prediction.
In the model, chemical structures, target proteins, action modes, and known side-effects were integrated into a unified framework for learning drug embeddings for side-effect prediction. However, sometimes the modes of action of drugs on target proteins are not always available and complete for all drugs. It has been proven that side-effects of drugs are related to the action modes. Consequently, in our future work, we will utilize known side-effects to predict modes of action of drugs and to further enrich the action mode information. We will also extend the datasets of drugs by mining from electronic medical records [37].

Author Contributions

Conceptualization, B.H.; Data curation, Z.Y.; Investigation, Z.Y.; Methodology, B.H.; Supervision, H.W.

Funding

The work is partially supported by the National Natural Science Foundation of China ( N o s . 61672329 , 61373149 , 61472233 , 61572300 , 81273704 ), Project of Shandong Province Higher Educational Science and Technology Program( N o . J 18 K A 370 ), Discipline Talent Team Cultivation Program of Shandong Women’s University ( N o . 209201904 ) and High-level Scientific Research Project Cultivation Fund of Shandong Women’s University ( N o . 2018 G S P S J 09 ). We also gratefully acknowledge the support of NVIDIA Corporation with the donation of the TITAN X GPU used for this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Giacomini, K.M.; Krauss, R.M.; Dan, M.R.; Eichelbaum, M.; Hayden, M.R.; Nakamura, Y. When good drugs go bad. Nature 2007, 446, 975–977. [Google Scholar] [CrossRef] [PubMed]
  2. Yamanishi, Y.; Pauwels, E.; Kotera, M. Drug side-effect prediction based on the integration of chemical and biological spaces. J. Chem. Inf. Model. 2012, 52, 3284–3292. [Google Scholar] [CrossRef]
  3. Li, J.; Zheng, S.; Chen, B.; Butte, A.J.; Swamidass, S.J.; Lu, Z. A survey of current trends in computational drug repositioning. Brief. Bioinform. 2015, 17, 2–12. [Google Scholar] [CrossRef] [Green Version]
  4. Xu, B.; Shi, X.; Zhao, Z.; Zheng, W. Leveraging biomedical resources in bi-lstm for drug–drug interaction extraction. IEEE Access 2018, 6, 33432–33439. [Google Scholar] [CrossRef]
  5. Vilar, S.; Tatonetti, N.P.; Hripcsak, G. 3D pharmacophoric similarity improves multi adverse drug event identification in pharmacovigilance. Sci. Rep. 2015, 5, 8809. [Google Scholar] [CrossRef] [PubMed]
  6. Rong, L.; Dong, Y.; Kuang, Q.; Wu, Y.; Li, Y.; Min, Z.; Li, M. Inductive matrix completion for predicting adverse drug reactions (adrs) integrating drug–target interactions. Chemom. Intell. Lab. Syst. 2015, 144, 71–79. [Google Scholar]
  7. Zhang, P.; Wang, F.; Hu, J.; Sorrentino, R. Label propagation prediction of drug–drug interactions based on clinical side effects. Sci. Rep. 2015, 5, 12339. [Google Scholar] [CrossRef] [PubMed]
  8. Zheng, Y.; Peng, H.; Ghosh, S.; Lan, C.; Li, J. Inverse similarity and reliable negative samples for drug side-effect prediction. BMC Bioinform. 2019, 19, 554. [Google Scholar] [CrossRef]
  9. Zhang, W.; Liu, F.; Luo, L.; Zhang, J. Predicting drug side effects by multi-label learning and ensemble learning. BMC Bioinform. 2015, 16, 365. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, M.; Wu, Y.; Chen, Y.; Sun, J.; Zhao, Z.; Chen, X.W.; Matheny, M.E.; Xu, H. Large-scale prediction of adverse drug reactions using chemical, biological, and phenotypic properties of drugs. J. Am. Med. Inform. Assoc. 2012, 19, e28–e35. [Google Scholar] [CrossRef] [Green Version]
  11. Yan, S.; Xu, D.; Zhang, B.; Zhang, H.J.; Yang, Q.; Lin, S. Graph embedding and extensions: A general framework for dimensionality reduction. IEEE Trans. Pattern Anal. Mach. Intell. 2006, 29, 40–51. [Google Scholar] [CrossRef]
  12. Cao, S. Deep Neural Network foR Learning Graph Representations. In Proceedings of the Thirtieth Aaai Conference on Artificial Intelligence, Phoenix, AZ, USA, 12–17 February 2016. [Google Scholar]
  13. Huang, Z.; Mamoulis, N. Heterogeneous information network embedding for meta path based proximity. arXiv 2017, arXiv:1701.05291. [Google Scholar]
  14. Ma, T.; Xiao, C.; Zhou, J.; Wang, F. Drug similarity integration through attentive multi-view graph auto-encoders. arXiv 2018, arXiv:1804.10850. [Google Scholar]
  15. Hu, B.; Wang, H.; Wang, L.; Yuan, W. Adverse Drug Reaction Predictions Using Stacking Deep Heterogeneous Information Network Embedding Approach. Molecules 2018, 23, 3193. [Google Scholar] [CrossRef]
  16. Campillos, M.; Kuhn, M.; Gavin, A.C.; Jensen, L.J.; Bork, P. Drug target identification using side-effect similarity. Science 2008, 321, 263–266. [Google Scholar] [CrossRef] [PubMed]
  17. Mizutani, S.; Pauwels, E.; Stoven, V.; Goto, S.; Yamanishi, Y. Relating drug–protein interaction network with drug side effects. Bioinformatics 2012, 28, i522–i528. [Google Scholar] [CrossRef]
  18. Yamanishi, Y.; Kotera, M.; Moriya, Y.; Sawada, R.; Goto, S. Dinies: Drug–target interaction network inference engine based on supervised analysis. Nucleic Acids Res. 2014, 42, W39–W45. [Google Scholar] [CrossRef]
  19. Torres, N.B.; Altafini, C. Drug combinatorics and side effect estimation on the signed human drug–target network. BMC Syst. Biol. 2015, 10, 74. [Google Scholar] [CrossRef]
  20. Iacono, G.; Altafini, C. Monotonicity, frustration, and ordered response: An analysis of the energy landscape of perturbed large-scale biological networks. BMC Syst. Biol. 2010, 4, 83. [Google Scholar] [CrossRef]
  21. Iacono, G.; Ramezani, F.; Soranzo, N.; Altafini, C. Determining the distance to monotonicity of a biological network: A graph-theoretical approach. IET Syst. Biol. 2010, 4, 223. [Google Scholar] [CrossRef]
  22. Milo, R.; Shen-Orr, S.; Itzkovitz, S.; Kashtan, N.; Chklovskii, D.; Alon, U. Network motifs: Simple building blocks of complex networks. Science 2002, 298, 824–827. [Google Scholar] [CrossRef] [PubMed]
  23. Sontag, E.D. Monotone and near-monotone biochemical networks. Syst. Synth. Biol. 2007, 1, 59–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wang, H.; Zhang, F.; Min, H.; Xing, X.; Guo, M.; Qi, L. Shine: Signed heterogeneous information network embedding for sentiment link prediction. arXiv 2017, arXiv:1712.00732. [Google Scholar]
  25. Lovász, L. Random walks on graphs: A survey. In Combinatorics, Paul Erdos Is Eighty; János Bolyai Mathematical Society: Hungary, Budapest, 1993; pp. 353–397. [Google Scholar]
  26. Mikolov, T.; Sutskever, I.; Chen, K.; Corrado, G.; Dean, J. Distributed representations of words and phrases and their compositionality. Adv. Neural Inf. Process. Syst. 2013, 26, 3111–3119. [Google Scholar]
  27. Knox, C.; Law, V.; Jewison, T.; Liu, P.; Wishart, D.S. Drugbank 3.0: A comprehensive resource for ‘omics’ research on drugs. Nucleic Acids Res. 2010, 39, D1035–D1041. [Google Scholar] [CrossRef] [PubMed]
  28. Hu, B.; Wang, H.; Yu, X.; Yuan, W.; He, T. Sparse network embedding for community detection and sign prediction in signed social networks. J. Ambient Intell. Hum. Comput. 2019, 10, 175–186. [Google Scholar] [CrossRef]
  29. Grover, A.; Leskovec, J. Node2vec: Scalable feature learning for networks. arXiv 2016, arXiv:1607.00653. [Google Scholar]
  30. Wang, Y.; Xiao, J.; Suzek, T.O.; Zhang, J.; Wang, J.; Bryant, S.H. Pubchem: A public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009, 37, W623–W633. [Google Scholar] [CrossRef]
  31. Kuhn, M.; Campillos, M.; Letunic, I.; Jensen, L.J.; Bork, P. A side effect resource to capture phenotypic effects of drugs. Mol. Syst. Biol. 2010, 6, 343. [Google Scholar] [CrossRef]
  32. Tatonetti, N.P.; Ye, P.P.; Daneshjou, R.; Altman, R.B. Data-driven prediction of drug effects and interactions. Sci. Transl. Med. 2012, 4, 125ra31. [Google Scholar] [CrossRef]
  33. Hastie, T.; Tibshirani, R.; Friedman, J.; Franklin, J. The elements of statistical learning: Data mining, inference and prediction. Math. Intell. 2005, 27, 83–85. [Google Scholar]
  34. Belkin, M. Laplacian eigenmaps and spactral techniques for embedding and clustering. Adv. Neural Inf. Process. Syst. 2002, 14, 585–591. [Google Scholar]
  35. Kipf, T.N.; Welling, M. Semi-supervised classification with graph convolutional networks. arXiv 2016, arXiv:1609.02907. [Google Scholar]
  36. Ishide, T. Denopamine, a selective beta1-receptor agonist and a new coronary vasodilator. Curr. Med. Res. Opin. 2002, 18, 407–413. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, L.; Wang, H.; Song, Y.; Wang, Q. MCPL-Based FT-LSTM: Medical Representation Learning-Based Clinical Prediction Model for Time Series Events. IEEE Accesss 2019, 7, 70253–70264. [Google Scholar] [CrossRef]
Figure 1. Flowchart of our side-effects prediction model (random walk on a signed heterogeneous information network (RW-SHIN)).
Figure 1. Flowchart of our side-effects prediction model (random walk on a signed heterogeneous information network (RW-SHIN)).
Molecules 24 03668 g001
Figure 2. Index plot of the number of associated drugs for each side-effect.
Figure 2. Index plot of the number of associated drugs for each side-effect.
Molecules 24 03668 g002
Figure 3. Distribution of the edge weights w i j s . (The number of edges is 4544. The number of negative edges with weights w i j s < 0.5 is 1074, while the number of positive edges with weights w i j s > 0.5 is 3369.)
Figure 3. Distribution of the edge weights w i j s . (The number of edges is 4544. The number of negative edges with weights w i j s < 0.5 is 1074, while the number of positive edges with weights w i j s > 0.5 is 3369.)
Molecules 24 03668 g003
Figure 4. Signed drug sub-network. (Nodes are drugs and the number of edges is 4544. There are 3416 positive edges and 1128 negative edges. Lines in red color represent positive links while lines in green color represent negative links.)
Figure 4. Signed drug sub-network. (Nodes are drugs and the number of edges is 4544. There are 3416 positive edges and 1128 negative edges. Lines in red color represent positive links while lines in green color represent negative links.)
Molecules 24 03668 g004
Figure 5. Illustration of the biased random walk procedure in a signed graph. (The walker is now resting on node s and evaluating its next node h. Green dotted lines and red solid lines represent negative edges and positive edges, respectively. Blue lines with arrows indicate the next-hop node h.)
Figure 5. Illustration of the biased random walk procedure in a signed graph. (The walker is now resting on node s and evaluating its next node h. Green dotted lines and red solid lines represent negative edges and positive edges, respectively. Blue lines with arrows indicate the next-hop node h.)
Molecules 24 03668 g005
Figure 6. Illustration of the enemy sets and friend sets. (The overlap area is the intersection of sets.)
Figure 6. Illustration of the enemy sets and friend sets. (The overlap area is the intersection of sets.)
Molecules 24 03668 g006
Figure 7. Illustration of the biased random walk procedure in an unsigned graph. (The walker is now resting on node s and evaluating its next node h. t is the last step in the path. Red solid lines represent edges in the graph. Blue lines with arrows represent the next possible steps.)
Figure 7. Illustration of the biased random walk procedure in an unsigned graph. (The walker is now resting on node s and evaluating its next node h. t is the last step in the path. Red solid lines represent edges in the graph. Blue lines with arrows represent the next possible steps.)
Molecules 24 03668 g007
Figure 8. Skip-gram model.
Figure 8. Skip-gram model.
Molecules 24 03668 g008
Figure 9. Receiver operator characteristic (ROC) curves based on the 5-fold cross-validation experiment.
Figure 9. Receiver operator characteristic (ROC) curves based on the 5-fold cross-validation experiment.
Molecules 24 03668 g009
Figure 10. Boxplots of the area under the ROC curve (AUROC) scores for every side-effect.
Figure 10. Boxplots of the area under the ROC curve (AUROC) scores for every side-effect.
Molecules 24 03668 g010
Figure 11. Impact of the embedding dimension on the prediction performance.
Figure 11. Impact of the embedding dimension on the prediction performance.
Molecules 24 03668 g011
Table 1. Modes of action of drugs and the corresponding edge signs [15].
Table 1. Modes of action of drugs and the corresponding edge signs [15].
Edge SignAction Modes in DrugBank
Positive (+)agonist; partial agonist; activator; stimulator; inducer; positive allosteric modulator; potentiator; positive modulator
Negative (−)Inhibitor; inhibitory allosteric modulator; inhibitor competitive; antagonist; partial antagonist; negative modulator; inverse agonist; blocker; suppressor; desensitize the target; neutralizer; reducer
Not classifiable (0)antibody; cofactor; modulator; binder; chaperone; cleavage; metabolizer; ligand; product of; component of; chelator; cross-linking/alkylation; intercalation; adduct; acetylation; allosteric modulator
Table 2. Number of triads in the signed drug subnetwork.
Table 2. Number of triads in the signed drug subnetwork.
Number of TriadsNumber of T1 (+ + +)Number of T2 (− − +)Number of T3 (+ + −)Number of T4 (− − −)
33,47019,700 (58.86%)12,266 (36.65%)1216 (3.63%)288 (0.86%)
Table 3. Numbers of drugs whose t o p K prediction results contain one known side-effect. (Bold values are the best prediction results.)
Table 3. Numbers of drugs whose t o p K prediction results contain one known side-effect. (Bold values are the best prediction results.)
Method T o p 1 T o p 3 T o p 5
random20 ± 2.4852 ± 3.4697 ± 3.44
Laplacian eigenmaps228 ± 8.08316 ± 7.23429 ± 6.45
GCN249 ± 7.28343 ±6.39441 ± 5.27
AttSemiGAE264 ± 4.79365 ± 5.81460± 6.42
RW-HIN258 ±8.12328 ± 7.47454 ± 7.18
RW-SHIN276 ±5.48354 ± 6.32465 ± 8.93
Table 4. Prediction of the o p 10 side-effects of Betaxolol and Dobutamine based on RW-SHIN.
Table 4. Prediction of the o p 10 side-effects of Betaxolol and Dobutamine based on RW-SHIN.
Top K BetaxololConfirmedDobutamineConfirmed
K = 1 dysphasiayeshypokalemiayes
K = 2 perforated gastric ulcernosneezingno
K = 3 nauseayespruritusyes
K = 4 headacheyesnauseayes
K = 5 mood disordersyesinsomniano
K = 6 shoulder painnostreptococcal pharyngitisno
K = 7 hypercholesterolemiayesheartburnyes
K = 8 otitis externayesskin necrosisyes
K = 9 tumornoadrenal diseaseno
K = 10 ear painyesdyskinesiano

Share and Cite

MDPI and ACS Style

Hu, B.; Wang, H.; Yu, Z. Drug Side-Effect Prediction Via Random Walk on the Signed Heterogeneous Drug Network. Molecules 2019, 24, 3668. https://doi.org/10.3390/molecules24203668

AMA Style

Hu B, Wang H, Yu Z. Drug Side-Effect Prediction Via Random Walk on the Signed Heterogeneous Drug Network. Molecules. 2019; 24(20):3668. https://doi.org/10.3390/molecules24203668

Chicago/Turabian Style

Hu, Baofang, Hong Wang, and Zhenmei Yu. 2019. "Drug Side-Effect Prediction Via Random Walk on the Signed Heterogeneous Drug Network" Molecules 24, no. 20: 3668. https://doi.org/10.3390/molecules24203668

APA Style

Hu, B., Wang, H., & Yu, Z. (2019). Drug Side-Effect Prediction Via Random Walk on the Signed Heterogeneous Drug Network. Molecules, 24(20), 3668. https://doi.org/10.3390/molecules24203668

Article Metrics

Back to TopTop