Next Article in Journal
Mechanisms of Chemopreventive and Therapeutic Proprieties of Ginger Extracts in Cancer
Next Article in Special Issue
Quantum Artificial Neural Network Approach to Derive a Highly Predictive 3D-QSAR Model for Blood–Brain Barrier Passage
Previous Article in Journal
The Intestinal Microbiota: Impacts of Antibiotics Therapy, Colonization Resistance, and Diseases
Previous Article in Special Issue
Ligand-Receptor Interactions and Machine Learning in GCGR and GLP-1R Drug Discovery
Article

Predicting Drug–Target Interactions Based on the Ensemble Models of Multiple Feature Pairs

by 1, 2, 2,* and 1,3,4,*
1
Department of Computer Science & Technology, Tongji University, Shanghai 201804, China
2
Institutes of Physical Science and Information Technology & School of Internet, Anhui University, Hefei 230601, China
3
School of Electrical & Information Engineering, Anhui University of Technology, Ma’anshan 243032, China
4
Key Laboratory of Power Electronics and Motion Control Anhui Education Department, Ma’anshan 243032, China
*
Authors to whom correspondence should be addressed.
Academic Editor: Irina Moreira
Int. J. Mol. Sci. 2021, 22(12), 6598; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22126598
Received: 15 May 2021 / Revised: 9 June 2021 / Accepted: 16 June 2021 / Published: 20 June 2021
(This article belongs to the Special Issue Artificial Intelligence and Machine Learning in Drug Development)

Abstract

Backgroud: The prediction of drug–target interactions (DTIs) is of great significance in drug development. It is time-consuming and expensive in traditional experimental methods. Machine learning can reduce the cost of prediction and is limited by the characteristics of imbalanced datasets and problems of essential feature selection. Methods: The prediction method based on the Ensemble model of Multiple Feature Pairs (Ensemble-MFP) is introduced. Firstly, three negative sets are generated according to the Euclidean distance of three feature pairs. Then, the negative samples of the validation set/test set are randomly selected from the union set of the three negative sets in the validation set/test set. At the same time, the ensemble model with weight is optimized and applied to the test set. Results: The area under the receiver operating characteristic curve (area under ROC, AUC) in three out of four sub-datasets in gold standard datasets was more than 94.0% in the prediction of new drugs. The effectiveness of the proposed method is also shown with the comparison of state-of-the-art methods and demonstration of predicted drug–target pairs. Conclusion: The Ensemble-MFP can weigh the existing feature pairs and has a good prediction effect for general prediction on new drugs.
Keywords: drug–target interactions; ensemble model of Multiple Feature Pairs (Ensemble-MFP); model weight sum; support vector machines drug–target interactions; ensemble model of Multiple Feature Pairs (Ensemble-MFP); model weight sum; support vector machines

1. Introduction

The prediction of drug–target Interaction (DTI) based on machine learning is very important in pharmacology and drug design [1,2,3]. It can also be considered as one direction in chemogenomics, which is a new interdisciplinary subject of biology, chemistry and informatics [4,5]. Traditional DTI methods are time-consuming, costly, and make it difficult to obtain three-dimensional structures of compounds and proteins [6,7,8]. The technology of machine learning accelerates the development of drug–target interactions, especially in reducing the blindness of experiments [9,10,11,12,13].
The characteristics of imbalance datasets in drug–target interaction predictions restrict the development of machine learning [1,14,15,16,17,18]. In the datasets of DTIs, the drug–target pairs with identified interactions which are labeled positive are sparse. At the same time, there are no validated negative samples, that is, non-interaction, in most databases [19]. In other words, the datasets of DTIs cannot provide enough reliable positive and negative samples for machine learning to obtain stable models [20]. To solve this problem, extraction methods of negative samples were studied. The random sampling method is used for negative extraction in various papers, which randomly selects negative samples from unlabeled sets [1,21]. Other negative sampling methods were also discussed. Liu et al. assumed that the negative samples can be extracted by their dissimilar characteristics from positive ones [22]. Hu et al. introduced the method based on Euclidean distance for negative sampling, and obtained better predictions [7]. Moderlet et al. introduced a bootstrap aggregating technique for negative sampling in Positive-Unlabeled (PU) problems [23].
The effective feature pairs selection in DTIs is another problem that restricts machine learning [1,24]. There are many types of features that can describe the characteristics of drugs or target proteins. The feature pair of DTI can be defined as the combination of one or more drug descriptors and one or more target descriptors. The dimensions of drug descriptors and target descriptors can be different or the same according to the feature extraction methods. Researchers have explored many types of feature pairs to predict drug–target interactions. Wei et al. predicted the interactions combined with 881-dimensional drug-descriptors, and target-descriptors of 567-dimensionals and 1449-dimensionals [25]. Bahi et al. combined 193-dimensional drug-descriptors based on the RCPI package and 1290-dimensional target-descriptors from PROFEAT to predict the interactions [11]. Feng et al. proposed the Deep Belief Network (DBN) for DTI based on 6144-dimensional Extended-Connectivity Fingerprints (ECFP) of drugs, and the 8420-dimensional Protein Sequence Composition (PSC) [26].
The prediction based on the Ensemble models of Multiple Feature Pairs (Ensemble-MFP) for new drugs prediction is studied in this work. Negative sampling based on the Euclidean distance, which is used to obtain the most dissimilar samples compared with positive sampling, is highly dependent on the calculated feature pairs, and in particular, the prediction of negative sampling based on different feature pairs is more prone to bias. At the same time, considering that the basic feature pairs of DTI are not clear and it is difficult to discover new feature pairs, an ensemble of the models based on existing feature pairs to have better predictions is necessary. The construction of a validation set and test set is designed to ensure the generalization ability of the algorithm and avoid the problem of overfitting. The final model is the weighted sum of sub-models corresponding to three feature pairs in this work, and the weights are optimized. Finally, the results on the test sets show that the algorithm is effective. Through the prediction of independent datasets by the proposed model, some drug–target pairs with interactions were predicted, which shows that the ensemble model has a good predictive effect on new drugs (see Appendix C). At the same time, we also provide several groups of drug–target pairs that may have interactions for further research in wet-lab. It should be noted that although the research regards the drug–target interaction as a binary classification problem, the actual situation is more complex with the strength of interactions, inhibitor or agonist, and so forth. Therefore, our model has limitations in broader predictions.

2. Materials and Methods

2.1. Benchmark Datasets

The benchmark dataset used in this work is the Gold Standard Dataset, which was first introduced by Yamanishi et al. It was collected and constructed in 2008, from KEGG BRITE, DrugBank, BRENDA, and SuperTarget [19,27,28,29,30]. According to the different characteristics of the target protein, it was divided into four sub-datasets: the enzyme, GPCR, ion channel, and nuclear receptor. It is publicly available on http://web.kuicr.kyoto-u.ac.jp/supp/yoshi/drugtarget/, accessed on 1 July 2008. Table 1 shows their statistical information in detail. It can be seen from the Table that the number of positive samples is far less than that of unlabeled samples, that is, the data are seriously imbalanced. It is very important that the prediction research needs reliable negative sample information.

2.2. Evaluation Criteria

The Area Under the Curve for the receiver operating characteristic (Area Under ROC, AUC), is the performance criteria used in this work. The metrics, such as Accuracy, Precision, Recall, and so forth are sufficient in general classification problems, but hold no significance in imbalanced datasets [1]. Some of the parameters used for evaluation are calculated as follows:
A c c u r a c y = T P + T N T P + T N + F P + F N
P r e c i s i o n = T P T P + F P
R e c a l l = T P T P + F N
F 1 _ s c o r e = 2 × P r e c i s i o n × R e c a l l P r e c i s i o n + R e c a l l ,
where TP, TN, FP, and FN are true positive, true negative, false positive and false negative, respectively. In these parameters, “positive” and “negative” represent drug–target pairs labeled as interaction or non-interaction in the benchmark dataset. At the same time, “true” and “false” mean that the prediction of the drug–target pair is right or wrong. ROC curves are drawn according to the True Positive Rate (TPR) and False Positive Rate (FPR) of different thresholds in the classification, and are recommended for comprehensive evaluation, especially in imbalance classification. AUC is the area under the ROC curve, and can be easily compared. It ranges from 0 to 1, and the larger the value, the better the model. Application of AUC can be found in most papers related with classification [1,17,31,32,33].

2.3. Negative Sampling and Data Construction

Negative samples are mainly generated based on the Euclidean distance in this work. Different from the random sampling method, the Euclidean distance-based sampling method holds that the farther the sample is from the positive center, the more reliable the negative sample is [7]. Its formula is as follows:
D i s = ( p o s d , t u n l a b e l d , t ) 2 ,
where p o s d , t denotes the positive samples’ center of the mean calculation. u n l a b e l d , t denotes the unlabeled samples. All unlabeled samples will be sorted according to their distance from the center of positive samples ( D i s ). The larger the D i s is, the more reliable the negative samples are. In order to avoid the negative sample difference caused by different feature units, all features used are firstly normalized. At the same time, Principal Components Analysis (PCA) is performed to avoid the interference of correlation in the calculation of the Euclidean distance. Although the sampling method is effective, it has a high dependence on the selected feature pairs and is difficult to be generalized, especially for the negative samples generated by different feature pairs, or is selected randomly. In order to improve the generalization ability of the model and obtain better prediction results, this work designs an ensemble model method of multiple feature pairs.
Data construction is based on 5-fold cross-validation. In order to make the model reliable for new drugs, the drugs in the dataset are divided with the ratio of 0.6, 0.2, and 0.2, respectively. In other words, DT pairs are divided into the training set, validation set, and test set according to different drugs. This work uses three feature pairs to get three corresponding models, and generates more general negative samples in the validation set and test set. The negative sampling process of the validation set and test set is shown in Figure 1.
Firstly, in the validation set/test set, based on Euclidean distance calculation, three feature pairs are used for negative sampling. Then, the negative samples were combined into U-vali/U-test. Random selection from U-vali/U-test can get more general negative samples for the validation set/test set. For the training set, the three feature pairs are trained respectively by the method of negative sampling based on Euclidean distance, and three models are obtained. According to these three models, the validation set is weighted and optimized to get a better weight vector, which is applied to the test set. The calculation formula of ensemble models is as follows,
d e c = w i × d e c i ,
where w i represents the weight of feature pair i. d e c i represents the decision vectors predicted by m o d e l i . d e c and m o d e l i denote the optimized decision vector in the validation set/test set and the model trained according with feature pair i. The flowchart of the proposed algorithm is shown in Figure 2.

2.4. Feature Pairs and Algorithm

Feature pairs used in Ensemble-MFP are extracted from the PaDEL-Descriptor and PROFEAT. The PaDEL-Descriptor is a free software for generating drug-descriptors, and is available on https://www.winsite.com/, accessed on 12 October 2010 [34]. PROFEAT is a webservice for calculation protein features, and can be used on http://bidd.group/, accessed on 12 April 2011 [35,36]. Table 2 lists the feature pairs used in the proposed method, which have better predictability in all sub-datasets of the gold standard dataset. In the Table, Estate-FP, MACCS-FP and Sub-FP Count are shorts for Electrotopological State Fingerprints, MACCS Fingerprints, and the Substructure Fingerprints Count, respectively. AAC, APAAC and QSO are short for Amino Acid Composition, Amphiphilic Pseudo-Amino Acid Composition, and Quasi-Sequence-Order descriptors, respectively.
Support vector machines (SVM) and its toolbox Libsvm (version 3.23) are adopted in this work. The Radial Basis Function (RBF) kernel, which can easily process the nonlinear classification problems, is used, and the kernel function only needs to adjust two parameters, c and γ . The parameters are adjusted in the form of an exponent, with the bottom of 2 [37]. Finally, the optimized parameters have good performance in four sub-datasets, that is, c = 2 4 and γ = 2 7 .

3. Results

3.1. Performance on DTIs

The ROC curve is shown in Figure 3, which represents the predictions in the validation set and the test set. All DT-pairs containing these drugs are omitted from the training set. Similarly, predictions on different targets are shown in Figure A1. It is shown that the prediction results of the test sets are very close to those of the validation set, which proves that there are no overfitting problems in this work. More evaluation information about these predictions is listed in Table 3.

3.2. Comparison with State-of-the-Art Methods

Various methods based on the gold standard dataset are compared. Table 4 shows the average results of the proposed method, and compares four feature vector-based algorithms based on the same dataset, such as that by Wang et al., Multi-scale Features Deep Representations (MFDR), Cao et al., and FRnet-DTI [1,9,25,38]. In these methods, the predictions were obtained by 5-fold cross-validation which were the same as our method. Wang et al. used the stacked autoencoder of deep learning based on the drug molecular structure and protein sequence to predict interactions between drugs and targets. Based on the large-scale drug/target features reconstructed by the autoencoder, SVM is used to predict drug–target interactions in the MFDR method. Cao et al. predicted interactions between the drugs and the target proteins according to the MACCS substructure fingerprint of the drug and the amino acid composition, Composition (C), Transformation (T), and Distribution (D) of the target protein. FRnet-DTI is composed of two convolutional neural networks, FRnet-Encode and FRnet-Predict, for feature manipulation and classification. Except for MFDR, other methods only segment the DT-pairs, and do not consider whether there are drugs that have been trained in the test set; we reproduce these models based on the algorithm in their original paper, and test the drug segmentation test set mentioned in this proposed work. In addition, considering that the negative samples of these algorithms are based on random sampling, we also verify the negative samples of the test set (ran-proposed), and the results are shown in the Table 4. The results show that this method has the best prediction effect in GPCR and ion channel. In enzymes, the predicted results of FRnet-DTI were only 1.3% higher than that of the proposed method. Considering the FRnet-DTI algorithm using two convolutional neural networks for feature extraction and prediction, this method is simple to implement and has closed results. For nuclear receptors, the average results are poor with all the compared algorithms, which may be due to the small dataset and lack of enough training information. The lack of information in nuclear receptor also makes the results unstable, as shown in Figure 4. The predictions based on random sampling (ran-proposed) are also comparable with other methods. In Table 4, the best prediction for each sub-dataset is marked as bold.

4. Discussion

4.1. Robustness of Prediction

Robustness of the proposed method is discussed. To show the effectiveness and stability of the proposed algorithm, the experiments were carried out 20 times, and the fluctuations of AUCs are shown in Figure 4. It can be seen that, except for the nuclear receptor, the other three sub-datasets have stable predictions in both the validation set and test set.

4.2. Weight Optimization of Ensemble Models

The weights of different feature pairs are optimized to obtain better predictions. AUC is the evaluation criteria used in the optimization process. In the process of optimization, w 1 ( 0 w 1 1 ) , w 2 ( 0 w 2 ( 1 w 1 ) ) and w 3 ( w 3 = 1 w 1 w 2 ) represent the weights of feature pair 1, feature pair 2, and feature pair 3, respectively. It can be seen from Figure 5 that the prediction results vary with the different weight sequences, which proves the rationality of the Ensemble-MFP algorithm in this work. The maximum predicted results correspond to the optimized weight sequence ( w 1 = 0.1, w 2 = 0.2, w 3 = 0.7).

4.3. Comparison between Ensemble Models and Individual Model

It is shown that the prediction results based on the Ensemble models of multiple feature pairs are better than the individual feature pair model in the test set. For each sub-dataset in Figure 6, Ensemble represents the predictions based on the Ensemble-MFP method, and Fea-1, Fea-2, and Fea-3 represent the results based on only feature pair 1, feature pair 2, and feature pair 3, respectively. In order to make the comparison reliable, all the positive and negative samples used in the training set, validation set, and test set in this part are the same. It is shown that the result of ensemble models is better than that of the individual model with each feature pair, which proves the superiority of the ensemble design. In addition, even if the multiple feature pairs used in the ensemble model are connected to form longer features with weights, better prediction results cannot be obtained, because the ensemble model can simulate more general negative samples (Figure A2).

4.4. External Validation

External validation dataset is used to demonstrate the effectiveness of the proposed method. The datasets used in DeepDTI [21], which was extracted based on DrugBank, is used for external validation. At the same time, the independent dataset extracted from the Drug Mechanism of ChEMBL, retaining the inhibitors and Homo sapiens, is tested [39]. After removing the same drugs of the gold standard dataset in the model training, two external datasets were tested with random negative samples. In Table 5, “DeepDTI” denotes the results in their original paper, “proposed-DeepDTI” and “proposed-ChEMBL”, which represent the results on the two external validation datasets based on the proposed method. The results in Table 5 shows the effectiveness of our proposed method, and TPR and TNR represent the True Positive Rate (TP/(TP + FN)) and True Negative Rate (TN/(TN + FP)), respectively. In addition, two predicted drug–target pairs were demonstrated as interactions (Lysine (DB00194) interacts with SLC7A4 (O43246) [40,41], and Micafungin (DB01141) interacts with FKSA (A2QLK4)) [42,43].

5. Conclusions

In this work, an algorithm based on the Ensemble models of Multiple Feature Pairs (Ensemble-MFP) is proposed for drug–target interaction predictions. Three models are obtained through three feature pairs, and the weights of the models are optimized on the validation set and applied on the test set. In order to make the model more general, the negative samples in the validation set/test set are collected randomly from three negative sets, which are extracted based on the Euclidean distance of three feature pairs. It is shown that, compared with the individual model of the single feature pair on the test set in the algorithm, the prediction effects of the Ensemble-MFP are better, which proves the effectiveness of the method. In addition, according to the external validation and demonstration results of the predicted DT pairs, the proposed method has a significance contribution on the drug design.
The algorithm can be further extended based on the details of more feature pairs. For the sake of simplicity, only three feature pairs are studies in this work. In addition, more feature pairs can be added to the algorithm. At the same time, according to the drug–target pairs predicted, we believe that our algorithm will supply more potential DT-pairs for wet-lab people, and motivate more researchers to study DTI in depth. Finally, the binary classification method restricts the further development of DTI to a certain extent, which will be the development direction in the future.

Author Contributions

Conceptualization, methodology and software, C.W.; writing—original draft preparation, C.W.; writing—review and editing, J.Z.; supervision, P.C. and B.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Nos. 61472282, and 61672035), Educational Commission of Anhui Province (No. KJ2019ZD05), Anhui Province Funds for Excellent Youth Scholars in Colleges (gxyqZD2016068), and Anhui Scientific Research Foundation for Returned Scholars.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The benchmark dataset used in the method was first introduced by Yamanishi et al., and was collected based on KEGG BRITE, DrugBank, BRENDA and SuperTarget. It can be downloaded freely at http://web.kuicr.kyoto-u.ac.jp/supp/yoshi/drugtarget/, accessed on 1 July 2008, and “Adjacency matrix of the gold standard drug–target interaction data” are used in the proposed method. The gold standard data set had been used in many studies related to drug–target interactions, and it is easy to compare the predicted results. PaDEL-Descriptor is a free software to extract drug descriptors, which is produced by the Pharmaceutical Data Exploration Laboratory and can be downloaded from https://www.winsite.com/, accessed on 12 October 2010. The method adopts version 2.21. PROFEAT is a webservice for computing the characteristics of related proteins on http://bidd.group/, accessed on 12 April 2011. Libsvm version 3.23 is used in this method. It was produced by Lin et al and can be downloaded in https://www.csie.ntu.edu.tw/~cjlin/libsvm/, accessed on 5 May 2003. The index parameter adjustment method with 2 at the bottom is based on the 2016 libsvm practical guide. The algorithm is implemented on MacOS High Sierra (10.13.6) platform by MATLAB R2018b. And the code is available on GitHub (https://github.com/Wangcheng0788/Ensemble-MFP, accessed on 16 June 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DTIdrug–target Interaction
Ensemble-MFPFusion of Multiple Feature Pairs
AUCArea Under the Curve for ROC
ROCReceiver Operating Characteristics
PUPositive-Unlabeled problems
GPCRG Protein-Coupled Receptors
SVMSupport Vector Machines
RBFRadial Basis Function

Appendix A. Predictions for New Targets

Similar prediction results are obtained when the target proteins do not appear repeatedly in training set, validation set and test set, which proves the effectiveness on new target protein prediction of the Ensemble-MFP. The prediction results are shown in Figure A1, and Ensemble represents the results of proposed method, Fea-1, Fea-2, Fea-3 denote the predictions based on the model of feature pair 1, feature pair 2 and feature pair 3 respectively.
Figure A1. The predictions for new targets. The experiments are based on different targets of training set, validation set and test set.
Figure A1. The predictions for new targets. The experiments are based on different targets of training set, validation set and test set.
Ijms 22 06598 g0a1

Appendix B. Comparison between Simple Connection of Feature Pairs and Ensemble-MFP

The prediction effect of three feature pairs on simple combination is worse than that of different models trained by three features. The difference between the two methods lies in the construction of negative samples. The simple connection of multiple feature pairs is similar to the result of a single feature pair, and cannot predict more general negative samples. In contrast, the method proposed in this work, on the one hand, simulates more general negative samples to a certain extent by combining and randomly selecting three groups of negative samples; on the other hand, by optimizing the weights, the model can get better prediction in general samples. We test and demonstrate several different experimental situations, including: (1) the case of simple feature connection, we test the negative samples based on Euclidean distance screening (Sim); (2) the case of simple feature connection, the more general negative samples designed in this work are tested (Gener-Sim); (3) the proposed results (Proposed). For the sake of fairness, in case (2), we also optimize the weights of each simply connected features. It can be seen from the Figure A2 that although the result of negative samples based on single feature extraction is very good (Sim), it is difficult to achieve good prediction for more general negative samples (Gener-Sim) with longer feature forms. In contrast, the weighted method mentioned in this work is better.
Figure A2. Comparison between simple connection of feature pairs and Ensemble-MFP.
Figure A2. Comparison between simple connection of feature pairs and Ensemble-MFP.
Ijms 22 06598 g0a2

Appendix C. Predicted Drug–Target Pairs

The top ten predicted drug–target pairs in GPCR are listed in Table A1. We use the proposed algorithm to predict the drug–target pairs in the test set and rank them according to the decision values. After several times of algorithm prediction, 10 groups of drug–target pairs with high decision scores were selected, and the information was input into the drug database (DrugBank) for query. Through the query, we found that the two groups of predicted drug–target pairs recorded as unlabeled drug–target pairs in 2008 had interaction records in the database. In addition, Trimipramine (D00394) interacts with HTR1A and HTR1B, so there may be interaction between D00394 and HTR1F (hsa3355) [44,45,46,47]. In the Table A1, 8 out of 10 pairs of predicted drug targets can be further demonstrated by wet-lab people. The results show that this method can effectively predict new drugs and is of great significance for drug development.
Table A1. Top ten drug–target pairs predicted in GPCR.
Table A1. Top ten drug–target pairs predicted in GPCR.
GPCRDrugDrug NameTargetGene NameRecord of DatabaseScore
1D00394Trimipraminehsa3355HTR1F-1.0973
2D00563Mirtazapinehsa3355HTR1F-1.0623
3D00394Trimipraminehsa154ADRB2-1.0365
4D02566Maprotilinehsa3355HTR1F-1.0350
5D00563Mirtazapinehsa154ADRB2-0.9933
6D00483Propranololhsa3355HTR1F-0.9854
7D00394Trimipraminehsa147ADRA1BDrugBank0.9768
8D00394Trimipraminehsa1128CHRM1-0.9765
9D00563Mirtazapinehsa147ADRA1B-0.9696
10D00394Trimipraminehsa3350HTR1ADrugBank0.9666

References

  1. Rayhan, F.; Ahmed, S.; Mousavian, Z.; Farid, D.M.; Shatabda, S. FRnet-DTI: Deep convolutional neural network for drug–target interaction prediction. Heliyon 2020, 6, e03444. [Google Scholar] [CrossRef]
  2. Parsons, A.B.; Brost, R.L.; Ding, H.; Li, Z.; Zhang, C.; Sheikh, B.; Brown, G.W.; Kane, P.M.; Hughes, T.R.; Boone, C. Integration of chemical-genetic and genetic interaction data links bioactive compounds to cellular target pathways. Nat. Biotechnol. 2004, 22, 62. [Google Scholar] [CrossRef] [PubMed]
  3. He, Z.; Zhang, J.; Shi, X.H.; Hu, L.L.; Kong, X.; Cai, Y.D.; Chou, K.C. Predicting drug–target interaction networks based on functional groups and biological features. PLoS ONE 2010, 5, e9603. [Google Scholar] [CrossRef] [PubMed]
  4. Claes, R.A.; Mats, G.G.; Helena, S. Quantitative Chemogenomics: Machine-Learning Models of Protein-Ligand Interaction. Curr. Top. Med. Chem. 2011, 11, 1978–1993. [Google Scholar]
  5. Bredel, M.; Jacoby, E. Chemogenomics: An emerging strategy for rapid target and drug discovery. Nat. Rev. Genet. 2004, 5, 262–275. [Google Scholar] [CrossRef]
  6. Alaimo, S.; Pulvirenti, A.; Giugno, R.; Ferro, A. drug–target interaction prediction through domain-tuned network-based inference. Bioinformatics 2013, 29, 2004–2008. [Google Scholar] [CrossRef]
  7. Hu, S.; Xia, D.N.; Su, B.; Chen, P.; Li, J. A Convolutional Neural Network System to Discriminate drug–target Interactions. IEEE/ACM Trans. Comput. Biol. Bioinform. 2019. [Google Scholar] [CrossRef] [PubMed]
  8. Jacob, L.; Hoffmann, B.; Stoven, V.; Vert, J.P. Virtual screening of GPCRs: An in silico chemogenomics approach. BMC Bioinform. 2008, 9, 363. [Google Scholar] [CrossRef]
  9. Wang, L.; You, Z.H.; Chen, X.; Xia, S.X.; Liu, F.; Yan, X.; Zhou, Y.; Song, K.J. A Computational-Based Method for Predicting drug–target interactions by Using Stacked Autoencoder Deep Neural Network. J. Comput. Biol. 2018, 25, 361–373. [Google Scholar] [CrossRef]
  10. Cheng, F.; Liu, C.; Jiang, J.; Lu, W.; Li, W.; Liu, G.; Zhou, W.; Huang, J.; Tang, Y. Prediction of drug–target interactions and drug repositioning via network-based inference. PLoS Comput. Biol. 2012, 8, e1002503. [Google Scholar] [CrossRef]
  11. Bahi, M.; Batouche, M. Drug–target Interaction Prediction in Drug Repositioning Based on Deep Semi-Supervised Learning. In IFIP Advances in Information and Communication Technology; Springer: London, UK, 2018; Volume 522, pp. 302–313. [Google Scholar] [CrossRef]
  12. Gove, R.; Faytong, J. Machine Learning and Event-Based Software Testing: Classifiers for Identifying Infeasible GUI Event Sequences. Adv. Comput. 2012, 86, 109–135. [Google Scholar]
  13. Bleakley, K.; Yamanishi, Y. Supervised prediction of drug–target interactions using bipartite local models. Bioinformatics 2009, 25, 2397–2403. [Google Scholar] [CrossRef]
  14. Bing, W.; Fang, A.; Xue, S.; Kim, S.; Xiang, Z. DISCO2: A Comprehensive Peak Alignment Algorithm for Two-Dimensional Gas Chromatography Time-of-Flight Mass Spectrometry. In Lecture Notes in Computer Science, Proceedings of the Bio-Inspired Computing and Applications—7th International Conference on Intelligent Computing, ICIC 2011, Zhengzhou, China, 11–14 August 2011; Revised Selected Papers; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  15. Chen, H.; Zhang, Z. A semi-supervised method for drug–target interaction prediction with consistency in networks. PLoS ONE 2013, 8, e62975. [Google Scholar] [CrossRef]
  16. Mousavian, Z.; Khakabimamaghani, S.; Kavousi, K.; Masoudi-Nejad, A. drug–target interaction prediction from PSSM based evolutionary information. J. Pharmacol. Toxicol. Methods 2016, 78, 42–51. [Google Scholar] [CrossRef] [PubMed]
  17. Rayhan, F.; Ahmed, S.; Shatabda, S.; Farid, D.M.; Mousavian, Z.; Dehzangi, A.; Rahman, M.S. iDTI-ESBoost: Identification of Drug Target Interaction Using Evolutionary and Structural Features with Boosting. Sci. Rep. 2017, 7, 17731. [Google Scholar] [CrossRef] [PubMed]
  18. Ezzat, A.; Zhao, P.; Wu, M.; Li, X.L.; Kwoh, C.K. drug–target Interaction Prediction with Graph Regularized Matrix Factorization. IEEE/ACM Trans. Comput. Biol. Bioinform. 2017, 14, 646–656. [Google Scholar] [CrossRef]
  19. Yamanishi, Y.; Araki, M.; Gutteridge, A.; Honda, W.; Kanehisa, M. Prediction of drug–target interaction networks from the integration of chemical and genomic spaces. Bioinformatics 2008, 24, i232–i240. [Google Scholar] [CrossRef]
  20. Bleakley, K.; Biau, G.; Vert, J.P. Supervised reconstruction of biological networks with local models. Bioinformatics 2007, 23, i57–i65. [Google Scholar] [CrossRef]
  21. 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]
  22. Liu, H.; Sun, J.; Guan, J.; Zheng, J.; Zhou, S. Improving compound-protein interaction prediction by building up highly credible negative samples. Bioinformatics 2015, 31, i221–i229. [Google Scholar] [CrossRef] [PubMed]
  23. Mordelet, F.; Vert, J.P. A bagging SVM to learn from positive and unlabeled examples. Pattern Recognit. Lett. 2013, 37, 201–209. [Google Scholar] [CrossRef]
  24. Wang, C.; Wang, W.; Lu, K.; Zhang, J.; Wang, B. Predicting drug–target interactions with Electrotopological State Fingerprints and Amphiphilic Pseudo Amino Acid Composition. Int. J. Mol. Sci. 2020, 21, 5694. [Google Scholar] [CrossRef]
  25. Hu, P.W.; Chan, K.C.C.; You, Z.H. Large-scale prediction of drug–target interactions from deep representations. In Proceedings of the 2016 International Joint Conference on Neural Networks (IJCNN), Vancouver, BC, Canada, 24–29 July 2016. [Google Scholar]
  26. Feng, Q.; Dueva, E.; Cherkasov, A.; Ester, M. PADME: A Deep Learning-based Framework for drug–target Interaction Prediction. arXiv 2018, arXiv:1807.09741. [Google Scholar]
  27. Kanehisa, M. From genomics to chemical genomics: New developments in KEGG. Nucleic Acids Res. 2006, 34, D354–D357. [Google Scholar] [CrossRef] [PubMed]
  28. Wishart, D.S.; Knox, C.; Guo, A.C.; Cheng, D.; Shrivastava, S.; Tzur, D.; Gautam, B.; Hassanali, M. DrugBank: A knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008, 36, D901–D906. [Google Scholar] [CrossRef] [PubMed]
  29. Stefan, G.; Michael, K.; Mathias, D.; Monica, C.; Christian, S.; Evangelia, P.; Jessica, A.; Garcia, U.E.; Andreas, G.; Juhl, J.L. SuperTarget and Matador: Resources for exploring drug–target relationships. Nuclc Acids Res. 2008, 36, D919–D922. [Google Scholar]
  30. Ida, S.; Chang, A.; Christian, E.; Marion, G.; Christian, H.; Gregor, H.; Dietmar, S. BRENDA, the enzyme database: Updates and major new developments. Nucleic Acids Res. 2004, 32, D431–D433. [Google Scholar]
  31. Ezzat, A.; Wu, M.; Li, X.L.; Kwoh, C.K. drug–target interaction prediction using ensemble learning and dimensionality reduction. Methods 2017, 129, 81–88. [Google Scholar] [CrossRef]
  32. Lee, I.; Nam, H. Identification of drug–target interaction by a random walk with restart method on an interactome network. BMC Bioinform. 2018, 19, 208. [Google Scholar] [CrossRef]
  33. Ozturk, H.; Ozkirimli, E.; Ozgur, A. A comparative study of SMILES-based compound similarity functions for drug–target interaction prediction. BMC Bioinform. 2016, 17, 128. [Google Scholar] [CrossRef]
  34. Yap, C.W. PaDEL-descriptor: An open source software to calculate molecular descriptors and fingerprints. J. Comput. Chem. 2011, 32, 1466–1474. [Google Scholar] [CrossRef] [PubMed]
  35. Li, Z.R.; Lin, H.H.; Han, L.Y.; Jiang, L.; Chen, X.; Chen, Y.Z. PROFEAT: A web server for computing structural and physicochemical features of proteins and peptides from amino acid sequence. Nucleic Acids Res. 2006, 34, W32–W37. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, P.; Tao, L.; Zeng, X.; Qin, C.; Chen, S.Y.; Zhu, F.; Yang, S.Y.; Li, Z.R.; Chen, W.P.; Chen, Y.Z. PROFEAT Update: A Protein Features Web Server with Added Facility to Compute Network Descriptors for Studying Omics-Derived Networks. J. Mol. Biol. 2017, 429, 416–425. [Google Scholar] [CrossRef]
  37. Hsu, C.; Chang, C.; Lin, C. A practical guide to support vector classification. Bju Int. 2008, 101, 1396–1400. [Google Scholar]
  38. Cao, D.S.; Liu, S.; Xu, Q.S.; Lu, H.M.; Huang, J.H.; Hu, Q.N.; Liang, Y.Z. Large-scale prediction of drug–target interactions using protein sequences and drug topological structures. Anal. Chim. Acta 2012, 752, 1–10. [Google Scholar] [CrossRef]
  39. Anna, G.; Anne, H.; Michał, N.; Patrícia, B.A.; Jon, C.; David, M.; Prudence, M.; Francis, A.; Bellis, L.J.; Elena, C.U. The ChEMBL database in 2017. Nucleic Acids Res. 2016, 45, D945–D954. [Google Scholar]
  40. Collins, J.E.; Wright, C.L.; Edwards, C.A.; Davis, M. A genome annotation-driven approach to cloning the human ORFeome. Genome Biol. 2004, 5, 1–11. [Google Scholar] [CrossRef] [PubMed]
  41. Gerhard, D.S.; Wagner, L.; Feingold, E.A.; Shenmen, C.M.; Grouse, L.H.; Schuler, G.; Klein, S.L.; Old, S.; Rasooly, R.; Good, P.; et al. The status, quality, and expansion of the NIH full-length cDNA project: The Mammalian Gene Collection (MGC). Genome Res. 2004, 14, 2121–2127. [Google Scholar]
  42. Pel, H.J.; Winde, J.D.; Archer, D.B.; Dyer, P.S.; Hofmann, G.; Schaap, P.J.; Turner, G.; Vries, R.D.; Al Ba Ng, R.; Albermann, K. Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88. Nat. Biotechnol. 2007, 25, 221–231. [Google Scholar] [CrossRef]
  43. Damveld, R.A.; van Kuyk, P.A.; Arentshorst, M.; Klis, F.M.; van den Hondel, C.A.M.J.J.; Ram, A.F.J. Expression of agsA, one of five 1,3-α-d-glucan synthase-encoding genes in Aspergillus niger, is induced in response to cell wall stress. Fungal Genet. Biol. 2005, 42, 165–177. [Google Scholar] [CrossRef]
  44. Kawanishi, Y.; Harada, S.; Tachikawa, H.; Okubo, T.; Shiraishi, H. Novel mutations in the promoter and coding region of the human 5-HT1A receptor gene and association analysis in schizophrenia. Am. J. Med. Genet. 2010, 81, 434–439. [Google Scholar] [CrossRef]
  45. Nakhai, B.; Nielsen, D.; Linnoila, M.; Goldman, D. 2 Naturally Occurring Amino Acid Substitutions in the Human 5-HT1A Receptor: Glycine 22 to Serine 22 and Isoleucine 28 to Valine 28. Biochem. Biophys. Res. Commun. 1995, 210, 530–536. [Google Scholar] [CrossRef]
  46. Wright, C.D.; Chen, Q.; Baye, N.L.; Huang, Y.; Healy, C.L.; Kasinathan, S.; O’Connell, T.D. Nuclear alpha1-adrenergic receptors signal activated ERK localization to caveolae in adult cardiac myocytes. Circ. Res. 2008, 103, 992–1000. [Google Scholar] [CrossRef] [PubMed]
  47. Wright, C.D.; Wu, S.C.; Dahl, E.F.; Sazama, A.J.; O’Connell, T.D. Nuclear Localization Drives α1-Adrenergic Receptor Oligomerization and Signaling in Cardiac Myocytes. Cell. Signal. 2012, 24, 794–802. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The negative sampling process of the validation set and test set.
Figure 1. The negative sampling process of the validation set and test set.
Ijms 22 06598 g001
Figure 2. The flowchart of the Ensemble-MFP algorithm.
Figure 2. The flowchart of the Ensemble-MFP algorithm.
Ijms 22 06598 g002
Figure 3. ROC-curves predicted on the validation set and the test set of each sub-dataset. The experiments are based on different drugs of the training set, validation set, and test set.
Figure 3. ROC-curves predicted on the validation set and the test set of each sub-dataset. The experiments are based on different drugs of the training set, validation set, and test set.
Ijms 22 06598 g003
Figure 4. Fluctuations of AUCs in four sub-datasets.
Figure 4. Fluctuations of AUCs in four sub-datasets.
Ijms 22 06598 g004
Figure 5. Heatmap of various weights.
Figure 5. Heatmap of various weights.
Ijms 22 06598 g005
Figure 6. Comparison between Ensemble Models and Individual Model.
Figure 6. Comparison between Ensemble Models and Individual Model.
Ijms 22 06598 g006
Table 1. Statistics of gold-standard datasets.
Table 1. Statistics of gold-standard datasets.
EnzymeGPCRIon ChannelNuclear Receptor
Drugs44522321054
Targets6649520426
DTIs2926635147690
unlabeled DT-pairs292,55420,55041,3641314
Table 2. Feature Pairs used in Ensemble-MFP.
Table 2. Feature Pairs used in Ensemble-MFP.
Drug DescriptorDimensionTarget DescriptorDimension
Feature Pair 1Estate-FP79AAC20
Feature Pair 2MACCS-FP166APAAC80
Feature Pair 3Sub-FP Count307QSO160
Table 3. Prediction results of the proposed method.
Table 3. Prediction results of the proposed method.
EnzymeGPCRIon ChannelNuclear Receptor
Accuracy (%)89.92 ± 0.93 #96.50 ± 0.7085.01 ± 1.6884.32 ± 12.44
Precision (%)90.37 ± 0.9398.89 ± 0.1684.90 ± 1.6891.29 ± 13.70
Recall (%)100 ± 0.0097.14 ± 0.96100.00 ± 0.0089.68 ± 17.01
F1-scores (%)94.94 ± 0.7298.01 ± 0.2391.83 ± 0.9890.48 ± 8.72
AUC (%)95.92 ± 0.3994.32 ± 0.5795.97 ± 0.2683.87 ± 7.38
# The value in the Table means the average ± standard deviation.
Table 4. Comparison with state-of-the-art methods on the gold standard dataset. Proposed and ran-proposed represent the predictions with proposed work and random sampling, respectively.
Table 4. Comparison with state-of-the-art methods on the gold standard dataset. Proposed and ran-proposed represent the predictions with proposed work and random sampling, respectively.
AUCEnzymeGPCRIon ChannelNuclear Receptor
Wang et al.0.9160.8970.9070.775
MFDR0.9690.9040.9330.886
Cao et al.0.9380.8390.8750.809
FRnet-DTI0.9720.9120.9430.872
Proposed0.9590.9430.9600.839
ran-proposed0.9330.9080.9250.821
Table 5. Predictions on external validation datasets.
Table 5. Predictions on external validation datasets.
TPR (%)TNR (%)Accuracy (%)AUC (%)
DeepDTI82.2789.5385.8891.58
proposed-DeepDTI86.2388.6591.6993.01
proposed-ChEMBL90.0993.1890.5792.78
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Back to TopTop