Next Article in Journal
Prognostic Factors and Current Treatment Strategies for Renal Cell Carcinoma Metastatic to the Brain: An Overview
Next Article in Special Issue
GATCDA: Predicting circRNA-Disease Associations Based on Graph Attention Network
Previous Article in Journal
LAG-3 Blockade with Relatlimab (BMS-986016) Restores Anti-Leukemic Responses in Chronic Lymphocytic Leukemia
Previous Article in Special Issue
Circular RNA in Chemonaive Lymph Node Negative Colon Cancer Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method to Predict Drug-Target Interactions Based on Large-Scale Graph Representation Learning

1
The Xinjiang Technical Institute of Physics & Chemistry, Chinese Academy of Sciences, Urumqi 830011, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Xinjiang Laboratory of Minority Speech and Language Information Processing, Urumqi 830011, China
4
College of Computer Science and Software Engineering, Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(9), 2111; https://doi.org/10.3390/cancers13092111
Submission received: 27 March 2021 / Revised: 20 April 2021 / Accepted: 22 April 2021 / Published: 27 April 2021
(This article belongs to the Special Issue Circular RNAs: New Insights into the Molecular Biology of Cancer)

Abstract

:

Simple Summary

The traditional process of drug development is lengthy, time-consuming, and costly, whereas very few drugs ever make it to the clinic. The use of computational methods to detect drug side effects greatly reduces the deficiencies in drug clinical trials. Prediction of drug-target interactions is a key step in drug discovery and repositioning. In this article, we proposed a novel method for the prediction of drug-target interactions based on large-scale graph representation learning. This method can be helpful to researchers in clinical trials and drug research and development.

Abstract

Identification of drug-target interactions (DTIs) is a significant step in the drug discovery or repositioning process. Compared with the time-consuming and labor-intensive in vivo experimental methods, the computational models can provide high-quality DTI candidates in an instant. In this study, we propose a novel method called LGDTI to predict DTIs based on large-scale graph representation learning. LGDTI can capture the local and global structural information of the graph. Specifically, the first-order neighbor information of nodes can be aggregated by the graph convolutional network (GCN); on the other hand, the high-order neighbor information of nodes can be learned by the graph embedding method called DeepWalk. Finally, the two kinds of feature are fed into the random forest classifier to train and predict potential DTIs. The results show that our method obtained area under the receiver operating characteristic curve (AUROC) of 0.9455 and area under the precision-recall curve (AUPR) of 0.9491 under 5-fold cross-validation. Moreover, we compare the presented method with some existing state-of-the-art methods. These results imply that LGDTI can efficiently and robustly capture undiscovered DTIs. Moreover, the proposed model is expected to bring new inspiration and provide novel perspectives to relevant researchers.

Graphical Abstract

1. Introduction

Drug repositioning is the process of exploring the new effects of existing drugs except for the original indications for medical treatment. It is a direction with great opportunities and challenges. In addition, it has the advantages of low-cost, short-time and low-risk [1,2]. The drug-target interactions (DTIs) play an important role in drug discovery and drug repositioning. Accurate prediction of DTIs can improve the accuracy of drug clinical trials, thus greatly reducing the risks of experiments. For a long time, the accumulation of a large number of biological experimental data and related literature makes the biological database richer and richer, which provides a favorable condition for the use of computational methods.
Traditional computing methods are mainly divided into two categories: ligand-based methods and structure-based methods. However, structure-based approaches are limited when the 3D structures of the target protein are absent, and ligand-based approaches have low accuracy when there are only a few binding ligands for the target protein [3,4,5,6,7]. In recent years, the widespread recognition of data-driven methods has made machine learning algorithms widely used in biomolecular correlation prediction [8,9,10,11]. There are mainly four related methods of in-silico methods: machine learning-based methods, network-based methods, matrix factor-based methods, and deep learning-based methods [12,13,14]. For example, Ding et al. [15] used substructure fingerprints, physical and chemical properties of organisms, and DTIs as feature extraction methods and input features, and further used SVM for classification. Chen et al. [16] employed gradient boosting decision tree (GBDT) to predict drug-target interactions based on three properties, including IDs of the drug and target, the descriptor of drug and target, DTIs. Luo et al. [17] constructed a heterogeneous network to predict the potential DTIs by integrating the information of multiple drugs. Chen et al. [18] and Ji et al. [19] proposed a multi-molecular network model based on network embedding to predict novel DTIs. Liu et al. [20] proposed a model called NRLMF, which calculates the score of DTIs through logical matrix decomposition, where the properties of the drug and target are expressed in terms of their specificity. Zheng et al. [21] proposed to map the drug and target into a low-rank matrix and to establish the weighted similarity matrix, and solve the problem by using the small square algorithm. Wen et al. [22] used unsupervised learning to extract representations from the original input descriptors to predict DTIs.
Recently, the extensive application of non-Euclidean structured data in graph neural networks has led to various graph-based algorithms [23,24,25,26,27,28,29,30], such as graph convolution networks (GCN), graph attention networks (GAT), graph autoencoders (GAE), graph generative networks, graph spatial-temporal networks, etc. Based on the analysis of biological data, it is found that the biological data network has a good preference for the graph neural network. Gao et al. [31] used long short-term memory (LSTM) and graph convolutional networks (GCN) to represent protein and drug structures, to predict DTIs. Previous work has shown the preferable performance of graph neural network for DTIs [27,32], however, a single understanding of the data relationship between DTIs cannot mine out the hidden information of the graph data well. Therefore, it is necessary to explore the depth information of the drug and target protein through the graph neural network.
In the actual graph, the relationship between two nodes is complex, and the features of each node are usually composed of a variety of attributes. It is necessary to clearly understand the relationship between nodes. Therefore, the extraction of node features should be multi-angle and multi-dimensional. To solve these challenges, we propose a novel method to predict DTIs based on large-scale graph representation learning (LGDTI). Unlike previous graph-based neural network-based approaches, LGDTI aims to gain an in-depth understanding of known drugs and targets association networks through different graph-based representation learning methods. To extract hidden graph features of drugs and targets in a complex biological network, two types of graph representation learning were used to excavate them.

2. Materials and Methods

2.1. Datasets

In this article, the multi-graph data were collected from DrugBank5.0 [33]. DrugBank5.0 is an open, free, comprehensive database, including drug molecular structures, mechanisms, and drug-target interactions that are constantly being updated. We downloaded 11,396 known DTIs from Drugbank5.0, including 984 drugs and 635 proteins; 11,396 known DTIs are conducted as the benchmark dataset, and in training as the positive sample.

2.2. Drug Attribute Representation

The molecular structure of the drug was extracted from the DrugBank database. The molecular structure is complex and difficult to use directly. To facilitate the calculation of drug molecular structure, it was necessary to vectorize its molecular structure [34]. The molecular fingerprint [35] is an abstract representation of a molecule, which encodes a molecule as a series of bit vectors, in which each bit on the molecular fingerprint corresponds to a molecular fragment, as shown in Figure 1. For the drug data, RDKit [36] was selected to calculate the Morgan fingerprint of the drug molecule.

2.3. Protein Attribute Representation

Protein sequence information was extracted from the STRING database [37]. Proteins are important biological macromolecules. All proteins are polymers formed by the linkage of 20 different amino acids, including (Ala, Val, Leu, Ile, Met, Phe, Trp, Pro), (Gly, Ser, Thr, Cys, Asn, Gln, Tyr), (Arg, Lys, His), and (Asp, Glu). Subsequently, the k-mer method is used [38], and k is set to 3, which translates each protein sequence into a 64-dimensional (4 * 4 * 4) feature vector by calculating the occurrence frequency of each sub-sequence in the entire protein sequence.

2.4. Graph Convolutional Network for Drug-Target Interactions (DTIs)

A graph convolutional network (GCN) [39] is a semi-supervised approach that turns topological associations into topological diagrams. In the algorithm, the input of GCN is the structure of the graph and the characteristics of each node, and the output includes the results at the node level, the results at the graph level, and the pooling information at the node level. Consequently, it is widely used in non-Euclidean spaces.
Let us assume that we have a bipartite graph G = with V = [ v 1 , , v n , , v m + n ] representing n drugs and m proteins, E = [ e i j ] representing the relationship of drug i and protein j . If e i j = 1 , v i and v j has a connection. Furthermore, in the graph the attributes of all nodes X = [ X d X p ] , the attributes of the drug X d = [ x 1 d , , x n d ] T and the attributes of the protein X p = [ x 1 p , , x m p ] T .
In this work, we define the function f ( X , A ) using the spatial method of GCN, where X is the feature set of each node, and A is the adjacency matrix. Therefore, the network communication rules of GCN are as follows:
f ( X , A ) l + 1 = σ ( D ˜ 1 2 A ˜ D ˜ 1 2 X l W l ) ,
in which, A ˜ = A + I n + m is the adjacency matrix added to the self-loop, D ˜ is represented as the degree matrix of A ˜ . W is the weight of the randomly initialized the network. σ represents the activation function of each layer of the neural network, here σ is R e L U ( · ) .
Although GCN has a natural preference for graph data, for DTIs data, we finally determined l = 1 and W is 64 * 64 after analysis and experiment. Then, in the initial training, we found that the algorithm had the problem of over-smoothing. To solve this challenge, we adjusted the defect of the original algorithm for this data. Specifically, after each convolution, we added node features for training, the formula is as follows:
f ( X , A ) l + 1 = Τ ( σ ( D ˜ 1 2 A ˜ D ˜ 1 2 X l W l ) , X ) ,
we adopted this adjusted graph convolution definition in this work.

2.5. Graph Embedding—DeepWalk for DTIs

DeepWalk [40] is a method to learn the potential representation of nodes in a graph and is a widely used algorithm in graph embedding. The main idea of the algorithm is divided into two parts. The first part is to sample the graph based on the random walk and map the node adjacency structure into sequence structure. The second part is to train the Skip-gram model by using the sequences obtained from sampling so that the expression of learning can capture the connectivity between nodes. Let us assume that we have a bipartite graph G = ( V , E ) . V is the set of nodes in the graph, and E is the edge of nodes. Each calculation starts from a given starting point, and then carries out a random walk through the sampled neighbor nodes, repeating the operation until the length of the sampled sequence is equal to the given maximum length, as shown in Algorithm 1.
S i = ( v i | ( v 1 , v 2 , v 3 , , v i 1 ) ) ,
where, S i is the random walk collection sequence, and v i is the random node.
Therefore, in the second part of the algorithm, S is computed by the Skip-gram model. Specifically, a two-layer neural network model is established. The input is the node sequence matrix of S n * m , and the weights in the neural network model are set as W 1 m * h and W 2 h * m respectively. Secondly, through backpropagation, the weight parameters are updated to obtain the representation of the target node, as shown in Algorithm 2.
S = [ S 1 m , S 2 m , , S n m ] T ,
Algorithm 1 DeepWalk ( G , w , d , γ , )
Input: graph G ( V , E )
  windows size w
  representation size d
  epoch γ
  step length t
Output: matrix of nodes representation ψ | V | × d
1: Initialization: ψ
2: Build a binary Tree T from V
3: for i = 0 to μ do
4:   V = Shuffle( V )
5: for each v i V do
6:    M v i = R a n d o m W a l k ( G , v i , t )
7:    S k i p G r a m ( ψ , M v i , w )
8: end for
9: end for
Algorithm 2 SkipGram ( ψ , M v i , w )
1: for each v j S v b do
2: for each u k S v i [ j w : j + w ] do
3:    J ( ψ ) = l o g P r ( u k | ψ ( v j ) )
4:    ψ = ψ α × J ψ
5: end for
6: end for

2.6. Construction of the Large-Scale Graph Representation Learning Network

Given a graph G ( V , E ) containing vertices V and edges E , where e i j is regard as a connection of v i and v j . a graph is considered as an adjacency matrix or an incidence matrix [41]. For an adjacency matrix A , A R N × N , is defined as:
A i j = { 1   if   ( v i , v j ) E 0   else   ,
Here, we used an undirected cycled graph, so a i i = 1 . For an incidence matrix B , B R N × M , is defined as:
B i j = { 1   if   v i   and   v j   are   connected 0   else   ,
The function of graph representation learning is to map data from complex graph space to multi-dimensional space. Its form is as follows:
f : V X d ,
where d | V | , V = [ v 1 , v 2 , v 3 , , v n + m ] is the original set of spatial variables and X = [ x 1 , x 2 , x 3 , , x d ] is the projected vector (or the embedded vector) that contains the structural information.
The first-order information is generally used to describe the local similarity between pairs of vertices in a graph [42]. Specifically, if there is an edge between two vertices, the two vertices should be close to each other in the embedded space. If there is no edge connection between two vertices, the first-order proximity between them is 0. Such work usually uses the KL-divergence [43] to calculate the distance by minimizing:
O 1 = ( i , j ) E w i j log p i ( v i , v j ) ,
in which p 1 ( v i , v j ) = 1 1 + exp ( v i T v j ) , v i and v j are the low-dimensional vector representation of the node v i and v j . W i j is the edge weight between node i and j . Although the methods based on the first-order neighbor of nodes are successful in graph embedding, they often fail to combine node substructure and node attributes for optimization. To address this challenge, the advantages of graph convolutional networks in vertex local feature extraction are utilized in Equation (1) to remedy this defect. An example of this algorithm is shown in Figure 2C.
The high-order information is learning the relationship between vertex v i and the other vertices separately [44,45]. Although there is no direct connection between the two vertices in the high-order information, learning that their representation vectors are close means that they should have similar or identical neighbors in the actual relational graph. For example, Figure 2B shows that drug d 1 has a second-order relationship with the target t 2 , drug d 2 and drug d 1 have a shared target t 1 , and target t 3 is a high-order potential candidate for drug d 1 . Then, we abstract high-order information (or global structure information) for each node by the graph embedding method: DeepWalk.
Consequently, we constructed a large-scale graph representation learning network to learn the features of each node, as shown in Figure 2. In which Figure 2A is the drug-target interactions sub-network.

2.7. The Large-Scale Graph Representation Learning DTI (LGDTI) Model Framework

In this study, the proposed LGDTI model contains not only first-order but also high-order graph information. In the first-order graph information, the graph convolutional network is used to capture the first-order neighbor information of the nodes in the graph; in the high-order graph information, the graph embedding algorithm DeepWalk is used to capture the high-order neighbor information of the nodes in the graph. Through these two different methods, the local and global information of each node in the graph is captured by LGDTI. The first-order neighbor information contains the attributes of nodes, which are internal to the node; the high-order neighbor information contains the whole network information of the node, which is called the behavior information. In the end, the two kinds of representation features of nodes obtained from LGDTI are predicted by the random forest classifier. The framework of large-scale graph representation learning as shown in Figure 3. In short, we have three main contributions: (i) we propose to employ specific GCN to learn first-order neighbors’ information (or local structural information) of nodes. (ii) This article proposes to utilize a graph embedding algorithm to learn high-order neighbors’ information (or global structural information) of nodes. (iii) In conclusion, LGDTI can view the DTIs network from multiple perspectives, including three features in the whole feature extraction process: node attributes, node first-order information, and node high-order information.

3. Results and Discussion

3.1. Performance Evaluation of LGDTI Using 5-Fold Cross-Validation

To accurately evaluate the stability and robustness of LGDTI, 5-fold cross-validation was adopted. In detail, the original data set was randomly divided into 5 subsets, among which 4 subsets were selected for each training, and the remaining subsets were used as the test set and repeated 5 times. Additionally, we used five evaluation indicators, including Acc. (Accuracy), MCC. (Matthews’s Correlation Coefficient), Sen. (Sensitivity), Spec. (Specificity), and Perc. (Precision). Moreover, for binary classification, the receiver operating characteristic (ROC) curve can reflect the capability of the model, while the AUC is the area under the ROC curve. The closer the ROC curve is to the upper left corner, the better the performance of the model. Similarly, the value of AUC is also high. The precision-recall (PR) curve contains precision and recall, with recall as the horizontal axis and precision as the vertical axis. On very skewed data sets, the PR curve can give us a comprehensive understanding of the performance of the model. The details of LGDTI under 5-fold cross-validation are shown in Table 1 and Figure 4. The results of each fold AUC, AUPR, and various evaluation criteria show that the proposed method has a better predictive ability. Studying it carefully, the results of each training are close to each other, which shows that the model has preferable stability and robustness.

3.2. Comparison LGDTI with the Different Machine Learning Algorithms

Different machine learning algorithms have different representations of features. By comparing different classification algorithms, including logistic regression (LR), K-nearest neighbor (KNN), gradient boosting decision tree (GBDT), and random forest classifier (RF), we can intuitively see the feature advantages of LGDTI. To make the comparison fairer and more objective, all classification algorithms choose the default parameters. The detailed evaluation results of 5-fold cross-validation are shown in Table 2 and Figure 5.
The results can be explained as follows: (i) for logistic regression, because of the depth and high complexity of input features, it may be difficult to form a linear classification surface, so it is impossible to fit features; (ii) for K-nearest Neighbor, in the characteristics of the sample studied in the early stage, the attributes of the neighboring nodes in the sample have been fused, which makes it impossible to accurately compare K-nearest neighbor; (iii) gradient boosting decision tree and random forest classifier are both ensemble classifiers, which can better solve the shortcomings of a single classifier, especially the random forest classifier, which can achieve preferable results on this dataset.

3.3. Comparison of the Different Feature with Attribute, GF and LGDTI

In summary, LGDTI constructs a graph and combines the first-order and high-order information of the nodes in the graph to denote the characteristics of each node. The first-order graph information aggregates the direct neighbor information of nodes. In graph theory, two nodes have similarities if the structure is similar to the subgraph. The high-order graph information provides a preferable representation of each node’s indirect neighbor information. Therefore, we conducted experiments on the different features of nodes, in which random forest classifier was used, as shown in Table 3 and Figure 6. In Table 3, Attribute has exemplified the feature of drug molecular structure and protein sequence; only first-order graph information is represented as GF; LGDTI includes the first-order and high-order graph information. When only node self-attributes are the worst, while self-attributes of nodes can be enhanced through GCN. Therefore, only the combination of first-order graph information and high-order graph information can better explore the potential features of nodes.

3.4. Compared with Existing State-of-the-Art Prediction Methods

To evaluate the advantage of the proposed method, it is compared with other advanced methods. Although the method proposed by Chen et al. [18] and Ji et al. [19], considers the network information of nodes, it fully expressed the local information of nodes in the network. Then, LGDTI is relatively sufficient for information extraction of nodes, and its high AUROC, AUPR, and ACC are stronger than other methods, as shown in Table 4.
Compared with other methods, node attributes (LGDTI (Only Attribute)), node first-order information (LGDTI (GF)), and the LDGTI model are all better. Among them, in the case of only node attributes, the AUROC, AUPR, and ACC of our model are at least 0.031, 0.0281, and 0.0259 higher respectively. Meanwhile, LGDTI (GF) still has some advantages. Definitively, the AUROC, AUPR, and ACC of the LGDTI model are at least 0.0222, 0.019, and 0.0281 higher than that of Ji et al. methods (Attribute+Behavior), respectively. The first-order neighborhood information aggregation makes node attribute characteristics are enhanced. Furthermore, the integration of first-order information and high-order information of the node will make our method have better prediction ability.

3.5. Case Studies

To test the practical ability of our model, the drugs clozapine and risperidone were exploited to predict potential targets, respectively. Clozapine can be used to treat many types of schizophrenia, and it can directly inhibit the brain stem reticulum up-activation system and has a powerful sedative and hypnotic effect. Risperidone is a psychiatric drug used to treat schizophrenia. In particular, it has an improved effect on the positive and negative symptoms and their accompanying emotional symptoms. It may also reduce the emotional symptoms associated with schizophrenia. In this case study, all known associations in the benchmark dataset were trained by our method, and we sorted the predicted scores of the remaining candidate targets and selected the top 5 targets, as shown in Table 5. The experiment showed that there were 3 targets of the drugs clozapine and risperidone predicted by LGDTI, which could be proved in the SuperTarget database [46]. The remaining unproven targets may be candidates, hopefully, to be explored by medical researchers.

4. Conclusions

Although the accurate and efficient computational model could greatly accelerate the process of identification of DTIs, there is still a huge gap between academia and industry. In this study, we developed a novel method called LGDTI for predicting DTIs. Specifically, the nodes in LGDTI can be represented by 2 kinds of feature including first-order information learned by GCN and high-order information learned by DeepWalk from the graph. in which molecular fingerprint technology was used to extract the attribute of drugs, and the k-mer method was used to extract the attribute of targets. Then, the Random Forest classifier was applied to carry out the relationship prediction task. The presented method obtained the AUC of 0.9455 and the AUPR of 0.9491 under 5-fold cross-validation which is more competitive than several state-of-the-art methods. Moreover, our method can learn three kinds of information about the node, including the node’s attributes, local structure, and global structure. Specifically, LGDTI can integrate attribute information with structural information for learning. The experimental results show that LGDTI has a prominent predictive ability for DTIs. Nevertheless, due to the limitation of the benchmark dataset, the performance of LGDTI cannot be shown collectively in multiple data. Moreover, LGDTI may be greatly improved if two kinds of node information can be better integrated. Consequently, we hope that the proposed model could be utilized to guide drug development and other biological wet experiments.

Author Contributions

B.-W.Z., Z.-H.Y. and L.H. considered the algorithm, arranged the dataset, and performed the analyses. Z.-H.G., L.W. (Lei Wang), Z.-H.C. and L.W. (Leon Wong) wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the grant of the National Key R&D Program of China (2018YFA0902600), and the grants of the National Science Foundation of China, Nos. 61722212, 61861146002 & 61732012.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dickson, M.; Gagnon, J.P. The cost of new drug discovery and development. Discov. Med. 2009, 4, 172–179. [Google Scholar]
  2. DiMasi, J.A.; Hansen, R.W.; Grabowski, H.G. The price of innovation: New estimates of drug development costs. J. Health Econ. 2003, 22, 151–185. [Google Scholar] [CrossRef] [Green Version]
  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. 2016, 17, 2–12. [Google Scholar] [CrossRef] [Green Version]
  4. Napolitano, F.; Zhao, Y.; Moreira, V.M.; Tagliaferri, R.; Kere, J.; D’Amato, M.; Greco, D. Drug repositioning: A machine-learning approach through data integration. J. Cheminform. 2013, 5, 30. [Google Scholar] [CrossRef] [Green Version]
  5. Wu, C.; Gudivada, R.C.; Aronow, B.J.; Jegga, A.G. Computational drug repositioning through heterogeneous network clustering. BMC Syst. Biol. 2013, 7, S6. [Google Scholar] [CrossRef] [Green Version]
  6. Kinnings, S.L.; Liu, N.; Buchmeier, N.; Tonge, P.J.; Xie, L.; Bourne, P.E. Drug discovery using chemical systems biology: Repositioning the safe medicine Comtan to treat multi-drug and extensively drug resistant tuberculosis. PLoS Comput. Biol. 2009, 5, e1000423. [Google Scholar] [CrossRef] [PubMed]
  7. Liu, Z.; Fang, H.; Reagan, K.; Xu, X.; Mendrick, D.L.; Slikker, W., Jr.; Tong, W. In silico drug repositioning–what we need to know. Drug Discov. Today 2013, 18, 110–115. [Google Scholar] [CrossRef]
  8. Bagherian, M.; Sabeti, E.; Wang, K.; Sartor, M.A.; Nikolovska-Coleska, Z.; Najarian, K. Machine learning approaches and databases for prediction of drug-target interaction: A survey paper. Brief. Bioinform. 2021, 22, 247–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Agamah, F.E.; Mazandu, G.K.; Hassan, R.; Bope, C.D.; Thomford, N.E.; Ghansah, A.; Chimusa, E.R. Computational/in silico methods in drug target and lead prediction. Brief. Bioinform. 2020, 21, 1663–1675. [Google Scholar] [CrossRef] [PubMed]
  10. Manoochehri, H.E.; Nourani, M. Drug-target interaction prediction using semi-bipartite graph model and deep learning. BMC Bioinform. 2020, 21, 1–16. [Google Scholar] [CrossRef] [PubMed]
  11. D’Souza, S.; Prema, K.; Balaji, S. Machine learning models for drug-target interactions: Current knowledge and future directions. Drug Discov. Today 2020, 25, 748–756. [Google Scholar] [CrossRef]
  12. Xue, H.; Li, J.; Xie, H.; Wang, Y. Review of drug repositioning approaches and resources. Int. J. Biol. Sci. 2018, 14, 1232. [Google Scholar] [CrossRef] [Green Version]
  13. Luo, H.; Li, M.; Yang, M.; Wu, F.-X.; Li, Y.; Wang, J. Biomedical data and computational models for drug repositioning: A comprehensive review. Brief. Bioinform. 2021, 22, 1604–1619. [Google Scholar] [CrossRef]
  14. Yella, J.K.; Yaddanapudi, S.; Wang, Y.; Jegga, A.G. Changing trends in computational drug repositioning. Pharmaceuticals 2018, 11, 57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ding, Y.; Tang, J.; Guo, F. Identification of drug-target interactions via multiple information integration. Inf. Sci. 2017, 418, 546–560. [Google Scholar] [CrossRef]
  16. Chen, J.; Wang, J.; Wang, X.; Du, Y.; Chang, H. Predicting Drug Target Interactions Based on GBDT. In Proceedings of the International Conference on Machine Learning and Data Mining in Pattern Recognition, New York, NY, USA, 15–19 July 2018; pp. 202–212. [Google Scholar]
  17. Luo, Y.; Zhao, X.; Zhou, J.; Yang, J.; Zhang, Y.; Kuang, W.; Peng, J.; Chen, L.; Zeng, J. A network integration approach for drug-target interaction prediction and computational drug repositioning from heterogeneous information. Nat. Commun. 2017, 8, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Chen, Z.-H.; You, Z.-H.; Guo, Z.-H.; Yi, H.-C.; Luo, G.-X.; Wang, Y.-B. Prediction of Drug-target Interactions from Multi-Molecular Network Based on Deep Walk Embedding Model. Front. Bioeng. Biotechnol. 2020, 8, 338. [Google Scholar] [CrossRef] [PubMed]
  19. Ji, B.-Y.; You, Z.-H.; Jiang, H.-J.; Guo, Z.-H.; Zheng, K. Prediction of drug-target interactions from multi-molecular network based on LINE network representation method. J. Transl. Med. 2020, 18, 1–11. [Google Scholar] [CrossRef]
  20. Liu, Y.; Wu, M.; Miao, C.; Zhao, P.; Li, X.-L. Neighborhood regularized logistic matrix factorization for drug-target interaction prediction. PLoS Comput. Biol. 2016, 12, e1004760. [Google Scholar] [CrossRef] [PubMed]
  21. Zheng, X.; Ding, H.; Mamitsuka, H.; Zhu, S. Collaborative Matrix Factorization with Multiple Similarities for Predicting Drug-Target Interactions. In Proceedings of the 19th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Chicago, IL, USA, 11–14 August 2013; pp. 1025–1033. [Google Scholar]
  22. Wen, M.; Zhang, Z.; Niu, S.; Sha, H.; Yang, R.; Yun, Y.; Lu, H. Deep-learning-based drug-target interaction prediction. J. Proteome Res. 2017, 16, 1401–1409. [Google Scholar] [CrossRef] [PubMed]
  23. Sun, M.; Zhao, S.; Gilvary, C.; Elemento, O.; Zhou, J.; Wang, F. Graph convolutional networks for computational drug development and discovery. Brief. Bioinform. 2020, 21, 919–935. [Google Scholar] [CrossRef]
  24. Zhou, J.; Cui, G.; Zhang, Z.; Yang, C.; Liu, Z.; Wang, L.; Li, C.; Sun, M. Graph neural networks: A review of methods and applications. arXiv 2018, arXiv:1812.08434. [Google Scholar]
  25. Wu, Z.; Pan, S.; Chen, F.; Long, G.; Zhang, C.; Philip, S.Y. A comprehensive survey on graph neural networks. IEEE Trans. Neural Netw. Learn. Syst. 2020, 32, 4–24. [Google Scholar] [CrossRef] [Green Version]
  26. Guo, Z.-H.; You, Z.-H.; Yi, H.-C. Integrative construction and analysis of molecular association network in human cells by fusing node attribute and behavior information. Mol. Ther. Nucleic Acids 2020, 19, 498–506. [Google Scholar] [CrossRef]
  27. Zhao, T.; Hu, Y.; Valsdottir, L.R.; Zang, T.; Peng, J. Identifying drug-target interactions based on graph convolutional network and deep neural network. Brief. Bioinform. 2020, 22, 2141–2150. [Google Scholar] [CrossRef]
  28. Lim, J.; Ryu, S.; Park, K.; Choe, Y.J.; Ham, J.; Kim, W.Y. Predicting drug-target interaction using a novel graph neural network with 3D structure-embedded graph representation. J. Chem. Inf. Modeling 2019, 59, 3981–3988. [Google Scholar] [CrossRef] [PubMed]
  29. Jiang, M.; Li, Z.; Zhang, S.; Wang, S.; Wang, X.; Yuan, Q.; Wei, Z. Drug-target affinity prediction using graph neural network and contact maps. RSC Adv. 2020, 10, 20701–20712. [Google Scholar] [CrossRef]
  30. Yue, X.; Wang, Z.; Huang, J.; Parthasarathy, S.; Moosavinasab, S.; Huang, Y.; Lin, S.M.; Zhang, W.; Zhang, P.; Sun, H. Graph embedding on biomedical networks: Methods, applications and evaluations. Bioinformatics 2020, 36, 1241–1251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Gao, K.Y.; Fokoue, A.; Luo, H.; Iyengar, A.; Dey, S.; Zhang, P. Interpretable Drug Target Prediction Using Deep Neural Representation. In Proceedings of the 27th International Joint Conference on Artificial Intelligence, Stockholm, Sweden, 13–19 July 2018; pp. 3371–3377. [Google Scholar]
  32. Torng, W.; Altman, R.B. Graph convolutional neural networks for predicting drug-target interactions. J. Chem. Inf. Modeling 2019, 59, 4131–4149. [Google Scholar] [CrossRef]
  33. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2017, 46, D1074–D1082. [Google Scholar] [CrossRef]
  34. Li, Y.; Liu, X.z.; You, Z.H.; Li, L.P.; Guo, J.X.; Wang, Z. A computational approach for predicting drug-target interactions from protein sequence and drug substructure fingerprint information. Int. J. Intell. Syst. 2021, 36, 593–609. [Google Scholar] [CrossRef]
  35. Rogers, D.; Hahn, M. Extended-connectivity fingerprints. J. Chem. Inf. Modeling 2010, 50, 742–754. [Google Scholar] [CrossRef] [PubMed]
  36. Landrum, G. Rdkit documentation. Release 2013, 1, 1–79. [Google Scholar]
  37. Szklarczyk, D.; Morris, J.H.; Cook, H.; Kuhn, M.; Wyder, S.; Simonovic, M.; Santos, A.; Doncheva, N.T.; Roth, A.; Bork, P. The STRING database in 2017: Quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2016, 45, D362–D368. [Google Scholar] [CrossRef] [PubMed]
  38. Rizk, G.; Lavenier, D.; Chikhi, R. DSK: K-mer counting with very low memory usage. Bioinformatics 2013, 29, 652–653. [Google Scholar] [CrossRef] [PubMed]
  39. Kipf, T.N.; Welling, M. Semi-supervised classification with graph convolutional networks. arXiv 2016, arXiv:1609.02907. [Google Scholar]
  40. Perozzi, B.; Al-Rfou, R.; Skiena, S. Deepwalk: Online Learning of Social Representations. In Proceedings of the 20th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Association for Computing Machinery, New York, NY, USA, 24–27 August 2014; pp. 701–710. [Google Scholar]
  41. Ding, C.H.; He, X.; Zha, H.; Gu, M.; Simon, H.D. A Min-Max Cut Algorithm for Graph Partitioning and Data Clustering. In Proceedings of the 2001 IEEE International Conference on Data Mining, IEEE, California, CA, USA, 29 November 2001; pp. 107–114. [Google Scholar]
  42. Cavallari, S.; Zheng, V.W.; Cai, H.; Chang, K.C.-C.; Cambria, E. Learning community embedding with community detection and node embedding on graphs. In Proceedings of the 2017 ACM on Conference on Information and Knowledge Management, Association for Computing Machinery, Singapore, 6–10 November 2017; pp. 377–386. [Google Scholar]
  43. Goldberger, J.; Gordon, S.; Greenspan, H. An Efficient Image Similarity Measure Based on Approximations of KL-Divergence between Two Gaussian Mixtures. In Proceedings of the Ninth IEEE International Conference on Computer Vision, Nice, France, 13–16 October 2003; Volume 1, pp. 487–493. [Google Scholar] [CrossRef] [Green Version]
  44. Zhou, C.; Liu, Y.; Liu, X.; Liu, Z.; Gao, J. Scalable graph embedding for asymmetric proximity. Proc. AAAI Conf. Artif. Intell. 2017, 31, 2942–2948. [Google Scholar]
  45. Yang, J.-H.; Chen, C.-M.; Wang, C.-J.; Tsai, M.-F. HOP-rec: High-order proximity for implicit recommendation. In Proceedings of the 12th ACM Conference on Recommender Systems, Vancouver British Columbia, Canada, 2 October 2018; pp. 140–144. [Google Scholar]
  46. Günther, S.; Kuhn, M.; Dunkel, M.; Campillos, M.; Senger, C.; Petsalaki, E.; Ahmed, J.; Urdiales, E.G.; Gewiess, A.; Jensen, L.J. SuperTarget and Matador: resources for exploring drug-target relationships. Nucleic Acids Res. 2007, 36, D919–D922. [Google Scholar] [CrossRef]
Figure 1. A schematic diagram of the drug molecular structure is constructed as bit vectors. A is the structure of a drug molecule, and B, C, and D are all substructures of the drug molecule, corresponding to the converted bit (represented by the small black box), respectively.
Figure 1. A schematic diagram of the drug molecular structure is constructed as bit vectors. A is the structure of a drug molecule, and B, C, and D are all substructures of the drug molecule, corresponding to the converted bit (represented by the small black box), respectively.
Cancers 13 02111 g001
Figure 2. An example of large-scale graph representation learning. (A) The schematic diagram of the relationship between drugs and targets. (B) An example of the graph embedding in drug-target interactions (DTIs). (C) An example of the graph convolutional network.
Figure 2. An example of large-scale graph representation learning. (A) The schematic diagram of the relationship between drugs and targets. (B) An example of the graph embedding in drug-target interactions (DTIs). (C) An example of the graph convolutional network.
Cancers 13 02111 g002
Figure 3. The flowchart of the proposed large-scale graph representation learning DTI (LGDTI). (a) A bipartite graph of DTIs. The solid black line is described as known DTIs, and the dashed red line is described as latent DTIs. (b) Part A constructed an adjacency graph containing a self-loop, in which green nodes are drugs and purple nodes are targets, and the information of first-order neighbors of each node is aggregated through graph convolutional network. Part B represented high-order information of each node in a bipartite graph by DeepWalk. (c) The two kinds of representation features are integrated. (d) Random forest classifier is trained and used for predicting new DTIs.
Figure 3. The flowchart of the proposed large-scale graph representation learning DTI (LGDTI). (a) A bipartite graph of DTIs. The solid black line is described as known DTIs, and the dashed red line is described as latent DTIs. (b) Part A constructed an adjacency graph containing a self-loop, in which green nodes are drugs and purple nodes are targets, and the information of first-order neighbors of each node is aggregated through graph convolutional network. Part B represented high-order information of each node in a bipartite graph by DeepWalk. (c) The two kinds of representation features are integrated. (d) Random forest classifier is trained and used for predicting new DTIs.
Cancers 13 02111 g003
Figure 4. The receiver operating characteristic (ROC) and precision-recall (PR) curves under 5-fold cross-validation.
Figure 4. The receiver operating characteristic (ROC) and precision-recall (PR) curves under 5-fold cross-validation.
Cancers 13 02111 g004
Figure 5. Comparison of the ROC and PR curves performed based on different machine learning classifier.
Figure 5. Comparison of the ROC and PR curves performed based on different machine learning classifier.
Cancers 13 02111 g005
Figure 6. Comparison of the ROC and PR curves performed by random forest classifier based on different features.
Figure 6. Comparison of the ROC and PR curves performed by random forest classifier based on different features.
Cancers 13 02111 g006
Table 1. Five-fold cross-validation results by random forest classifier.
Table 1. Five-fold cross-validation results by random forest classifier.
FoldAcc. (%)MCC (%)Sen. (%)Spec. (%)Prec. (%)AUC (%)
088.3677.1183.2593.4692.7293.93
188.6077.5483.9093.2992.5994.43
288.2276.8982.8593.6092.8394.66
388.4077.2283.1693.6492.9094.51
489.6179.5285.1894.0493.4695.23
Average88.64 ± 0.5677.66 ± 1.0783.67 ± 0.9393.61 ± 0.2892.90 ± 0.3394.55 ± 0.47
Table 2. Comparison of different machine learning classifier.
Table 2. Comparison of different machine learning classifier.
ClassifierAcc. (%)MCC (%)Sen. (%)Spec. (%)Prec. (%)AUC (%)
LR72.54 ± 1.2345.23 ± 2.4776.57 ± 1.2968.51 ± 1.4970.86 ± 1.2178.26 ± 0.78
KNN71.07 ± 1.1546.90 ± 1.8292.99 ± 0.6849.15 ± 2.6964.67 ± 1.0982.63 ± 0.46
GBDT84.98 ± 0.2370.23 ± 0.4180.54 ± 0.6589.41 ± 0.2688.38 ± 0.1991.62 ± 0.38
RF88.64 ± 0.5677.66 ± 1.0783.67 ± 0.9393.61 ± 0.2892.90 ± 0.3394.55 ± 0.47
Table 3. Comparison of different feature using random forest classifier.
Table 3. Comparison of different feature using random forest classifier.
FeatureAcc. (%)MCC (%)Sen. (%)Spec. (%)Prec. (%)AUC (%)
Attribute83.86 ± 0.3267.78 ± 0.6581.62 ± 0.6986.09 ± 0.5685.44 ± 0.4790.89 ± 0.38
GF84.28 ± 0.4668.76 ± 0.9080.67 ± 0.8987.90 ± 0.6286.96 ± 0.5691.41 ± 0.36
LGDTI88.64 ± 0.5677.66 ± 1.0783.67 ± 0.9393.61 ± 0.2892.90 ± 0.3394.55 ± 0.47
Table 4. Compared with existing state-of-the-art prediction methods.
Table 4. Compared with existing state-of-the-art prediction methods.
MethodsDatasetsAUROCAUPRACC
Ji et al. methods (Only Attribute)DrugBank0.87770.88280.8073
Chen et al. methods (Only Attribute)DrugBank0.8779N/A0.8127
LGDTI (Only Attribute)DrugBank0.90890.91090.8386
LGDTI (GF)DrugBank0.91410.91770.8428
Chen et al. methods (Only Behavior)DrugBank0.9206N/A0.8545
Ji et al. methods (Only Behavior)DrugBank0.92180.92860.8575
Ji et al. methods (Attribute+Behavior)DrugBank0.92330.93010.8583
LGDTIDrugBank0.94550.94910.8864
Table 5. New association prediction results for the top 5 targets with clozapine and risperidone.
Table 5. New association prediction results for the top 5 targets with clozapine and risperidone.
Drug NameTarget NameConfirmed
ClozapineAlpha-1D adrenergic receptorSuperTarget
Cytochrome P450 3A5SuperTarget
UDP-glucuronosyltransferase 1A1Unconfirmed
Solute carrier family 22 member 3Unconfirmed
Sodium-dependent serotonin transporterSuperTarget
RisperidoneAlpha-1D adrenergic receptorSuperTarget
Solute carrier family 22 member 8Unconfirmed
Cytochrome P450 2C19Unconfirmed
Sodium-dependent serotonin transporterSuperTarget
Potassium voltage-gated channel subfamily H member 2SuperTarget
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, B.-W.; You, Z.-H.; Hu, L.; Guo, Z.-H.; Wang, L.; Chen, Z.-H.; Wong, L. A Novel Method to Predict Drug-Target Interactions Based on Large-Scale Graph Representation Learning. Cancers 2021, 13, 2111. https://doi.org/10.3390/cancers13092111

AMA Style

Zhao B-W, You Z-H, Hu L, Guo Z-H, Wang L, Chen Z-H, Wong L. A Novel Method to Predict Drug-Target Interactions Based on Large-Scale Graph Representation Learning. Cancers. 2021; 13(9):2111. https://doi.org/10.3390/cancers13092111

Chicago/Turabian Style

Zhao, Bo-Wei, Zhu-Hong You, Lun Hu, Zhen-Hao Guo, Lei Wang, Zhan-Heng Chen, and Leon Wong. 2021. "A Novel Method to Predict Drug-Target Interactions Based on Large-Scale Graph Representation Learning" Cancers 13, no. 9: 2111. https://doi.org/10.3390/cancers13092111

APA Style

Zhao, B. -W., You, Z. -H., Hu, L., Guo, Z. -H., Wang, L., Chen, Z. -H., & Wong, L. (2021). A Novel Method to Predict Drug-Target Interactions Based on Large-Scale Graph Representation Learning. Cancers, 13(9), 2111. https://doi.org/10.3390/cancers13092111

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop