 Methodology article
 Open Access
 Published:
The computational prediction of drugdisease interactions using the dualnetwork L_{2,1}CMF method
BMC Bioinformatics volume 20, Article number: 5 (2019)
Abstract
Background
Predicting drugdisease interactions (DDIs) is timeconsuming and expensive. Improving the accuracy of prediction results is necessary, and it is crucial to develop a novel computing technology to predict new DDIs. The existing methods mostly use the construction of heterogeneous networks to predict new DDIs. However, the number of known interacting drugdisease pairs is small, so there will be many errors in this heterogeneous network that will interfere with the final results.
Results
A novel method, known as the dualnetwork L_{2,1}collaborative matrix factorization, is proposed to predict novel DDIs. The Gaussian interaction profile kernels and L_{2,1}norm are introduced in our method to achieve better results than other advanced methods. The network similarities of drugs and diseases with their chemical and semantic similarities are combined in this method.
Conclusions
Cross validation is used to evaluate our method, and simulation experiments are used to predict new interactions using two different datasets. Finally, our prediction accuracy is better than other existing methods. This proves that our method is feasible and effective.
Background
On average, it takes over a dozen years and approximately 1.8 billion dollars to develop a drug [1]. In addition, most drugs have strong side effects or undesirable effects on patients, so these drugs cannot be placed on the market. Therefore, many pharmaceutical companies resort to repositioning of existing drugs on the market [2]. Many known drugs can be found to have new effects for different diseases. In medicine, drug repurposing has two advantages. One advantage is that known drugs have already been approved by the US FDA (Food and Drug Administration) [3]. In other words, these drugs are safe to use. Another advantage is that the side effects of these drugs are known to medical scientists, so these side effects can be better controlled to achieve the desired therapeutic effect. Drug repurposing can help accelerate and facilitate the research and development process in the drug discovery pipeline [4].
The most important factor for drug repositioning is online biological databases. Many public databases, such as KEGG [5], STITCH [6], OMIM [7], DrugBank [8] and ChEMBL [9] store large amounts of information related to drugs and diseases. These databases contain detailed information such as a drug’s chemical structure, side effects, and genomic sequences [10].
In general, the goal of drug repositioning is to discover novel drugdisease interactions (DDIs) using existing drugs. Because a drug is often not specific for one disease, most drugs can treat a variety of diseases. Recently, more methods have been proposed for drug repositioning, such as machine learning [11], text mining [12], network analysis [13] and many other effective methods due to the increasing depth of research [14, 15]. Of course, we can also use the oppositionbased learning particle swarm optimization to predict interactions, such as SNPSNP interactions [16]. For instance, Gottlieb et al. proposed a computational method to discover potential drug indications by constructing drugdrug and diseasedisease similarity classification features [17]. Then, the predicted score of the novel DDIs can be calculated by a logistic regression classifier. Napolitano et al. calculated drug similarities using combined drug datasets [18]. They proposed a multiclass SVM (Support Vector Machine) classifier to predict some novel DDIs. Moreover, some researchers use networkbased models for drug repositioning. The advantage of this network model is that it can fully consider the largescale generation of highthroughput data to build complex biological information interaction networks. Wang et al. proposed a method called TLHGBI to infer novel treatments for diseases [19]. These authors constructed a heterogeneous network and integrated datasets about drugs, diseases and drug targets. Another networkbased prioritization method called DrugNet was proposed by Martinez et al. [20]. This method can predict not only novel drugs but also novel treatments for diseases. Similar to the TLHGBI method, the DrugNet method uses a heterogeneous network to predict novel DDIs using information about drugs, diseases, and targets. Luo et al. developed a computational method to predict novel interactions of known drugs [21]. Furthermore, comprehensive similarity measures and BiRandom Walk (MBiRW) algorithm have been applied to this method. In addition, Luo et al. continued to propose a drug repositioning recommendation system (DRRS) to predict new DDIs by integrating data sources for drugs and diseases [14]. A heterogeneous drugdisease interaction network can be constructed by integrating drugdrug, diseasedisease and drugdisease networks. Moreover, a large drugdisease adjacency matrix can replace the heterogeneous network, including drug pairs, disease pairs, known drugdisease pairs, and unknown drugdisease pairs. A fast and favourable algorithm SVT (Singular Value Thresholding) [22] has been used to complete predicted scores of the drugdisease adjacency matrix for unknown drugdisease pairs. According to previous studies, each method has its own advantages for predicting DDIs. However, after comparing the prediction of these methods, the best method is currently DRRS. The method achieves the highest AUC (area under curve) value and the best prediction [14]. Recently, matrix factorization methods have also been used to identify novel DDIs [23]. The matrix factorization method takes one input matrix and attempts to obtain two other matrices, and then the two matrices are multiplied to approximate the input matrix [23]. Similar to looking for missing interactions in the input matrix, matrix factorization can be used as a good technique to solve the prediction problem. Examples of such matrix factorization methods are the kernel Bayesian matrix factorization method (KBMF2K) [24] and the collaborative matrix factorization method (CMF) [25].
In this work, a simple yet effective matrix factorization model called the DualNetwork L_{2,1}CMF (Dualnetwork L_{2,1}collaborative matrix factorization) is proposed to predict new DDIs based on existing DDIs. However, there are many missing unknown interactions, so a preprocessing step is used to solve this problem. The main purpose of this preprocessing method is to attempt to weight K nearest known neighbours (WKNKN) [26]. Specifically, in the original matrix, WKNKN is used to describe whether there is an interaction between drugdisease pairs, bringing each element closer simply 0 and 1 to a reliable value than. Thus, WKNKN will have a positive impact on the final prediction. Furthermore, unlike the previous matrix factorization methods, L_{2,1}norm [2] and GIP (Gaussian interaction profile) kernels are added to the CMF method. Among them, L_{2,1}norm can avoid overfitting and eliminate some unattached disease pairs [27]. The GIP kernels are used to calculate the drug similarity matrix and the disease similarity matrix [28]. Cross validation is used to evaluate our experimental results. The final experimental results show that after removing some of the interactions, our proposed method is superior to other methods. In addition, a simulation experiment is conducted to predict new interactions.
The results are described in Section 2, including the datasets used in our study and experimental results. The corresponding discussions are presented in Section 3. The conclusion is described in Section 4. Finally, Section 5 describes our proposed method, including specific solution steps and iterative processes.
Results
DDIs datasets
Information about the drugs and diseases was obtained from Gottlieb et al. [17], and the Fdataset comprises multiple data sources. It is the gold standard dataset. This dataset includes 1933 DDIs, 593 drugs and 313 diseases in total. Further information about the drugs and diseases are obtained from Luo et al. [21], and the Cdataset comprises multiple data sources. The Cdataset includes 2353 DDIs, 663 drugs and 409 diseases, including drugs from the DrugBank database and diseases from OMIM (Online Mendelian Inheritance in Man) database [7].
Both datasets contain three matrices: Y ∈ ℝ^{n × m}, S_{D} ∈ ℝ^{n × n} and S_{d} ∈ ℝ^{m × m}. The adjacency matrix Y is proposed to describe the association between drug and disease. In the adjacency matrix, n drugs are represented in rows and m diseases are represented in columns. If drug D(i) is associated with disease d(j), the entity Y(D(i), d(j)) is 1; otherwise it is 0. Sparsity is defined as the ratio of the number of known DDIs to the number of all possible DDIs [14]. Table 1 lists the specific information for these two datasets.
Similarities in the chemical structures of the drugs
The drug similarity matrix is used to predict interactions. The chemical structure information of the drugs constitutes this matrix, S_{D}. The similarity information is derived from the Chemical Development Kit (CDK) [29], and the drugdrug pairs are represented as their 2D chemical fingerprint scores.
Similarities in disease semantics
The disease similarity matrix was used to predict interactions. The matrix S_{d} is represented by the medical descriptions of the diseases. The similarities between diseasedisease pairs were obtained from MimMiner [30]. Therefore, the semantic similarities of the diseases is achieved through text mining. Finally, the meaningful medical information is selected and meaningless data is discarded.
Cross validation experiments
In this study, our experiments are compared to the previous methods (KBMF, HGBI, DrugNet, MBiRW, and DRRS). For each method, 10fold cross validation is repeated ten times. However, before running our method, the preprocessing steps is performed first. The purpose is to solve the problem of missing unknown interactions. This preprocessing step improves the accuracy of the prediction to some extent.
We observe that the interactions between drugs and diseases remain fixed during crossvalidation. In general, the receiver operating characteristic (ROC) curve can be described by changing the true positive rate (TPR, sensitivity) of different levels of the false positive rate (FPR, 1specificity). Moreover, sensitivity and specificity (SPEC) can be written as follows:
where N represents the number of negative samples, TP represents the number of positive samples correctly classified by the classifier and FP represents the number of false positive samples classified by the classifier. Similarly, TN represents the number of negative samples correctly classified by the classifier, and FN represents the number of false negative samples.
A popular evaluation indicator AUC is used to evaluate our approach [31]. AUC is defined as the area under the ROC curve, and it is obvious that the value of this area will not be greater than 1. In general, the value of AUC ranges between 0.5 and 1. The AUC value cannot be less than 0.5. The drugdisease pairs are randomly removed from the interaction matrix Y before running cross validation. This method is called CVp (Cross Validation pairs), and its purpose is to increase the difficulty of the prediction, thereby enabling a more complete assessment of the ability to predict new drugs. In addition, cross validation is performed on the training set to establish the parameters λ_{l}, λ_{d} and λ_{t}. Grid search is used to find the best parameter from the values: λ_{l} ∈ {2^{−2}, 2^{−1}, 2^{0}, 2^{1}}, λ_{d}/λ_{t} ∈ {0, 10^{−4}, 10^{−3}, 10^{−2}, 10^{−1}}.
Prediction of the interaction under CVp
Table 2 lists the experimental results of CVp. The average of the AUC values of the ten cross validation results are taken as the final AUC score. Note that AUC is known to be insensitive to skewed class distributions [32]. The drug disease datasets are highly unbalanced in this study. In other words, there are more negative factors than positive factors. Therefore, the AUC value is a more appropriate measure to evaluate different methods. Table 2 shows the AUC values for different methods, and the best AUC value in each column is shown in bold. Standard deviations are shown in parentheses.
As shown in Table 2, our proposed method, DNL_{2,1}CMF, achieves an AUC of 0.951 on the Cdataset, which is 0.4% higher than DRRS, with an AUC of 0.947. The AUC value of the DrugNet method is the lowest, and our method is 14.7% higher than this value. In addition, our approach also achieves the best results for the Fdataset. Our method achieves an AUC of 0.94, which is 1% higher than DRRS, with an AUC of 0.93. Additionally, the AUC value of the DrugNet method is the lowest, and our method is 16.2% higher than this value. Therefore, our proposed method is better than other existing methods.
In summary, the advantage of our method lies in the introduction of GIP and L_{2,1}norm. GIP can obtain network information on drugs and diseases. L_{2,1}norm can remove undesired drug disease pairs, thus improving prediction accuracy. Some of the previous methods only considered a single drug similarity and a single disease similarity and did not consider their network information. Therefore, our method can achieve better AUC values.
Sensitivity analysis from WKNKN
As mentioned earlier in this paper, because there are some missing unknown interactions in the drug disease interaction matrix Y, a preprocessing method is used to minimize the error. The parameters K and p are fixed. K is the number of nearest known neighbours. p is a decay term where p ≤ 1, and WKNKN is used before running DNL_{2,1}CMF. When K = 5, p = 0.7, the AUC value approaches stability. The sensitivity analysis of these two parameters is shown in Figs. 1 and 2, respectively.
Discussion
Case study
In this subsection, a simulation experiment was conducted. Our method was used to predict the correct drugs in an unknown situation. Therefore, an unknown situation was created by removing some of the DDIs. Y was decomposed into two matrices, A and B, thus the product of these two matrices was used as the final prediction matrix. In this prediction matrix, all elements were no longer 0 and 1. Instead, all elements were close to 0 or 1. Therefore, we compared the elements in Y to determine the final prediction.
On the Cdataset, the seven pairs of interactions related to the drug zoledronic acid (KEGG ID: D01968) were completely removed. The drug was used to prevent skeletal fractures in patients with cancers such as multiple myeloma and prostate cancer. It can also be used to treat the hypercalcemia of malignancy and can be helpful for treating pain from bone metastases. A simulation was conducted to yield the prediction score matrix. Finally, the prediction score matrix counted whether those removed interactions were predicted. At the same time, the new interactions were counted. In other words, the disease most relevant to this drug was found. Among them, all known interactions and three novel interactions were successfully predicted. Table 3 lists the experimental results for the Cdataset. According to the level of relevance, these diseases were sorted from high to low. The known interactions are in bold. It is worth noting that according to our experimental analysis, the eighth disease, osteoporosis, had the strongest interaction with zoledronic acid. More information about the drug is published in DrugBank database.
The complete interactions of the drug hyoscyamine (KEGG ID: D00147) were removed. The drug is mainly used to treat bladder spasm, peptic ulcer disease, diverticulitis, colic, irritable bowel syndrome, cystitis and pancreatitis. This drug is also used to treat certain heart diseases and to control the symptoms of Parkinson’s disease and rhinitis. Fourteen pairs of interactions were removed, and these interactions were still predicted by our method. At the same time, motion sickness was predicted to be related to this drug. More information about the drug is published in https://www.drugbank.ca/drugs/DB00424. Table 4 lists the experimental results.
For the Fdataset, the interactions of the drug cisplatin and the drug dexamethasone were removed, and a simulation experiment was conducted. Table 5 lists the experimental results for cisplatin, and Table 6 lists the experimental results for dexamethasone.
For cisplatin (KEGG ID: D00275), nine interactions were removed. Six known interactions and three novel interactions were successfully predicted. The known interactions are shown in bold. More information about cisplatin is published at https://www.drugbank.ca/drugs/DB00515. For dexamethasone (KEGG ID: D00292), sixteen interactions were removed. Eleven known interactions and four novel interactions were successfully predicted. Moreover, endometriosis can be prevented by dexamethasone. In 2014, the ClinicalTrials.gov database was tested for this disease, and the reliability of this result has been confirmed by clinical trials. Sixtyfour participants were used in the experiment. Detailed experimental results can be found at https://clinicaltrials.gov/ct2/show/study/NCT02056717. Diseases ranked 12, 13, and 14 were not confirmed by ClinicalTrials.gov for treatment with dexamethasone.
According to the above simulation results, our method has good performance for different datasets. According to Table 3 to Table 6, it can be concluded that the advantages of the L_{2,1}norm are increasing the disease matrix sparsity and discarding unwanted disease pairs. This advantage is reflected in the fact that in a drugdisease pair, unwanted noise is removed by the L_{2,1}norm, so the vast majority of known DDIs that have been removed are successfully predicted. Therefore, the addition of GIP kernels and L_{2,1}norm achieved better results than other advanced methods.
Conclusions
In this paper, an effective matrix factorization model is proposed. L_{2,1}norm and GIP kernel are applied in this model. Moreover, the GIP kernel provides more network information for predicting novel DDIs. AUC is used to evaluate the indicators and our method achieves excellent results, so our method is feasible.
It is worth noting that the preprocessing method WKNKN plays an important role in prediction because there are many missing unknown interactions that are addressed by this preprocessing method. This is helpful for the final experimental results. However, the datasets used in this paper still have some limitations. For example, diseasedisease similarity, sequence similarity and GO similarity are not considered. We will collect more similarity information in future work.
In the future, more datasets will be available, and more novel DDIs will be predicted. Of course, we will continue to employ more machine learning methods or deep learning methods to solve drug development problems.
Methods
Problem formalization
Formally, the known interactions Y(D(i), d(j)) of drug D(i) associated with disease d(j) are considered to be a matrix factorization model. The input matrix Y is decomposed into two low rank matrices A and B. These two matrices retain the features of the original matrix. Then, the two matrices are optimized through constraints. Finally, the specific matrices of A and B are obtained. Our mission is to rank all of the drugdisease pairs Y(D(i), d(j)). The most likely interaction pairs have the highest ranking.
Gaussian interaction profile kernel
The method is based on the assumption that diseases that interact with DDIs networks and unrelated drugs in drugdisease networks may show similar interactions with new diseases. D(i) and D(j) represent two drugs, d(i) and d(j) represent two diseases. Their network similarity calculations can be written as:
where γ is a parameter, which is used to adjust the bandwidth of the kernel. In addition, Y(D_{i}) and Y(D_{j}) are the interaction profiles of D_{i} and D_{j}. Similarly, Y(d_{i}) and Y(d_{j}) are the interaction profiles of d_{i} and d_{j}. Then, the two network similarity matrices can be combined with S_{D} and S_{d} to be written as:
where α ∈ [0, 1] is an adjustable parameter. K_{D} is a drug kernel, which represents a linear combination of the drug chemical similarity matrix S_{D} and the drug network similarity matrix GIP_{D}. K_{d} is a disease kernel, which represents a linear combination of the disease semantic similarity matrix S_{d} and the disease network similarity matrix GIP_{d}. Thus, the network information is applied to the prediction of DDIs and performed well in yielding results.
Dualnetwork L_{2,1}collaborative matrix factorization (DNL_{2,1}CMF)
The traditional collaborative matrix factorization (CMF) uses collaborative filtering to predict novel interactions [25]. The objective function of CMF is given as follows:
where ‖⋅‖_{F} is the Frobenius norm and λ_{l}, λ_{d} and λ_{t} are nonnegative parameters.
CMF is an effective method for predicting DDIs. However, this method ignores the network information of drugs and diseases. This problem will reduce the accuracy of the CMF method in predicting novel DDIs.
In this study, an improved collaborative matrix factorization method is used to predict DDIs. The L_{2,1}norm is added to the collaborative matrix factorization method, and drug network information and disease network information are combined with this method. The interaction matrix Y is decomposed into two matrices A and B, where AB^{T} ≈ Y.The dualnetwork L_{2,1}collaborative matrix factorization (DNL_{2,1}CMF) method uses regularization terms to request that the potential feature vectors of similar drugs and similar diseases are similar, and the potential feature vectors of dissimilar drugs and dissimilar diseases are dissimilar [33], where S_{D} ≈ AA^{T} and S_{d} ≈ BB^{T}. Considering that GIP explores kernel network information, the dualnetwork can be interpreted as a drug network and a disease network generated by GIP. Specifically, the interaction profiles can be generated from a drugdisease interaction network. For a classifier, the interaction profiles can be used as feature vectors [34]. Therefore, the kernel method is used, and the kernel can be constructed from the interaction profiles. In summary, because of these advantages, GIP can achieve better results. Therefore, the objective function of DNL_{2,1}CMF method can be written as
where‖⋅‖_{F} is the Frobenius norm and λ_{l}, λ_{d} and λ_{t} are nonnegative parameters. The first term is an approximate model of the matrix Y, whose purpose is to search the latent feature matrices A and B. The Tikhonov regularization is used to minimizes the norms of A, B in the second term, whose purpose is to avoid overfitting. The L_{2,1}norm is applied in B in the third term. The purpose is to increase the sparsity of the disease matrix and discard unwanted disease pairs. For a detailed explanation, please refer to [2]. Based on a previous study [25], the effect of the last two regularization terms is to minimize the squared error between S_{D}(S_{d}) and AA^{T}(BB^{T}).
Initialization of A and B
For the input DDIs matrix Y, the singular value decomposition (SVD) method is used to obtain the initial value of matrix A and matrix B.
where S_{k} is a diagonal matrix and contains the k largest singular values. In addition, the minimization of the objective function is used to predict the outcome of the interactions, but this could lead to unsatisfactory results. Many zeros have not been found, so the WKNKN preprocessing method is used to solve this problem. Figure 3 shows a specific prediction flow chart from the original datasets to the final predicted score matrix.
Optimization algorithm
In this study, the least squares method is used to update A and B. First, L is represented as the objection function of DNL_{2,1}CMF method. Then, ∂L/∂A and ∂L/∂B are set to be 0. According to the alternating least squares method, A and B are updated until convergence. It is worth noting that λ_{l}, λ_{d} and λ_{t} are automatically determined by the cross validation on the training set to the optimal parameter values. Thus, the update rules are as follows:
According to formula (5) and formula (6), K_{D} can be represented by S_{D}, and K_{d} can be represented by S_{d}. These two complete updated rules can be written as:
where D is a diagonal matrix with the ith diagonal element as d_{ii} = 1/2‖(B)^{i}‖_{2}. Therefore, the specific algorithm of DNL_{2,1}CMF is as follows:
Abbreviations
 AUC:

Area Under Curve
 CMF:

Collaborative Matrix Factorization
 DDIs:

DrugDisease Interactions
 DNL_{2,1}CMF:

Dualnetwork L_{2,1}Collaborative Matrix Factorization
 DRRS:

Drug Repositioning Recommendation System
 GIP:

Gaussian Interaction Profile
 KBMF2K:

Kernel Bayesian Matrix Factorization
 MBiRW:

Measures and BiRandom Walk
 ROC:

Receive Operating Characteristic
 SVD:

Singular Value Decomposition
 SVT:

Singular Value Thresholding
 TPR:

True Positive Rate; FPR: False Positive Rate
 WKNKN:

Weight K Nearest Known Neighbours
References
Paul SM, Mytelka DS, Dunwiddie CT, Persinger CC, Munos BH, Lindborg SR, Schacht AL. How to improve R&D productivity: the pharmaceutical industry's grand challenge. Nat Rev Drug Discov. 2010;9(3):203–14.
Liu JX, Wang DQ, Zheng CH, Gao YL, Wu SS, Shang JL. Identifying drugpathway association pairs based on L2,1integrative penalized matrix decomposition. BMC Syst Biol. 2017;11(6):119.
Ezzat A, Wu M, Li XL, Kwoh CK. Drugtarget interaction prediction via class imbalanceaware ensemble learning. BMC Bioinformatics. 2016;17(19):509.
Novac N. Challenges and opportunities of drug repositioning. Trends Pharmacol Sci. 2013;34(5):267–72.
Kanehisa M, Goto S, Furumichi M, Mao T, Hirakawa M. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010;38(Database issue):355–60.
Kuhn M, Szklarczyk D, PletscherFrankild S, Blicher TH, Von MC, Jensen LJ, Bork P. STITCH 4: integration of proteinchemical interactions with user data. Nucleic Acids Res. 2014;42(Database issue):401–7.
Amberger J, Bocchini CA, Scott AF, Hamosh A. McKusick’s online Mendelian inheritance in man (OMIM). Nucleic Acids Res. 2009;37(Database issue):793–6.
Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, Pon A, Banco K, Mak C, Neveu V. DrugBank 3.0: a comprehensive resource for ‘omics’ research on drugs. Nucleic Acids Res. 2011;39(Database issue):D1035.
Gaulton A, Bellis LJ, Bento AP, Chambers J, Davies M, Hersey A, Light Y, Mcglinchey S, Michalovich D, Allazikani B. ChEMBL: a largescale bioactivity database for drug discovery. Nucleic Acids Res. 2012;40(Database issue):1100–7.
Banville DL. Mining chemical structural information from the drug literature. Drug Discov Today. 2006;11(1–2):35–42.
Chen X, Yan GY. Semisupervised learning for potential human microRNAdisease associations inference. Sci Rep. 2014;4:5501.
Yang H, Spasic I, Keane JA, Nenadic G. A text mining approach to the prediction of disease status from clinical discharge summaries. J Am Med Inform Assoc. 2009;16(4):596–600.
Oh M, Ahn J, Yoon Y. A networkbased classification model for deriving novel drugdisease associations and assessing their molecular actions. PLoS One. 2014;9(10):e111668.
Luo H, Li M, Wang S, Liu Q, Li Y, Wang J. Computational drug repositioning using lowrank matrix approximation and randomized algorithms. Bioinformatics. 2018;34(11):1904–12.
Zhang L, Xiao M, Zhou J, Yu J. Lineageassociated underrepresented permutations (LAUPs) of mammalian genomic sequences based on a jellyfishbased LAUPs analysis application (JBLA). Bioinformatics. 2018;34(21):3624–30.
Shang J, Sun Y, Li S, Liu JX, Zheng CH, Zhang J. An improved oppositionbased learning particle swarm optimization for the detection of SNPSNP interactions. Biomed Res Int. 2015;2015:524821.
Gottlieb A, Stein GY, Ruppin E, Sharan R. PREDICT: a method for inferring novel drug indications with application to personalized medicine. Mol Syst Biol. 2011;7(1):496.
Napolitano F, Zhao Y, Moreira VM, Tagliaferri R, Kere J, D’Amato M, Greco D. Drug repositioning: a machinelearning approach through data integration. J Cheminform. 2013;5(1):30.
Wang W, Yang S, Zhang X, Li J. Drug repositioning by integrating target information through a heterogeneous network model. Bioinformatics. 2014;30(20):2923–30.
Martínez V, Navarro C, Cano C, Fajardo W, Blanco A. DrugNet: networkbased drugdisease prioritization by integrating heterogeneous data. Artif Intell Med. 2015;63(1):41–9.
Luo H, Wang J, Li M, Luo J, Peng X, Wu FX, Pan Y. Drug repositioning based on comprehensive similarity measures and birandom walk algorithm. Bioinformatics. 2016;32(17):2664.
Cai JF, Cand S, Emmanuel J, Shen Z. A singular value thresholding algorithm for matrix completion. SIAM J Optim. 2008;20(4):1956–82.
Yang J, Li Z, Fan X, Cheng Y. Drug–disease association and drugrepositioning predictions in complex diseases using causal inference–probabilistic matrix factorization. J Chem Inf Model. 2014;54(9):2562–9.
Gönen M. Predicting drug–target interactions from chemical and genomic kernels using Bayesian matrix factorization. Bioinformatics (Oxford, England). 2012;28(18):2304–10.
Shen Z, Zhang YH, Han K, Nandi AK, Honig B, Huang DS. miRNAdisease association prediction with collaborative matrix factorization. Complexity. 2017;2017(9):1–9.
Ezzat A, Zhao P, Wu M, Li XL, Kwoh CK. Drugtarget interaction prediction with graph regularized matrix factorization. IEEE/ACM Trans Comput Biol Bioinformatics. 2017;14(3):646–56.
Liu JX, Wang D, Gao YL, Zheng CH, Shang JL, Liu F, Xu Y. A jointL 2,1 normconstraintbased semisupervised feature extraction for RNASeq data analysis. Neurocomputing. 2017;228(C):263–9.
Song M, Yan Y, Jiang Z. Drugpathway interaction prediction via multiple feature fusion. Mol BioSyst. 2014;10(11):2907–13.
Christoph Steinbeck, †, Yongquan Han, Stefan Kuhn, Oliver Horlacher, Edgar Luttmann A, Willighagen E: The chemistry development kit (CDK): an opensource Java library for chemo and bioinformatics. Cheminform 2003, 34(21):493–500.
Driel MA, Van JB, Gert V, Han G, Brunner LJAM. A textmining analysis of the human phenome. Eur J Hum Genet. 2006;14(5):535–42.
Grau J, Grosse I, Keilwagen J. PRROC: computing and visualizing precisionrecall and receiver operating characteristic curves in R. Bioinformatics. 2015;31(15):2595–7.
Fawcett T. An introduction to ROC analysis. Pattern Recogn Lett. 2006;27(8):861–74.
Ezzat A, Wu M, Li XL, Kwoh CK. Computational prediction of drugtarget interactions using chemogenomic approaches: an empirical survey. Brief Bioinform. 2018;8.
Laarhoven TV, Nabuurs SB, Marchiori E. Gaussian interaction profile kernels for predicting drug–target interaction. Bioinformatics. 2011;27(21):3036–43.
Acknowledgements
Not applicable.
Funding
This work was supported in part by grants from the National Science Foundation of China, Nos. 61872220 and 61572284.
Availability of data and materials
The datasets that support the findings of this study are available in https://github.com/cuizhensdws/drugdiseasedatasets/.
Author information
Authors and Affiliations
Contributions
ZC and JXL jointly contributed to the design of the study. ZC designed and implemented the DNL_{2,1}CMF method, performed the experiments, and drafted the manuscript. JW participated in the design of the study and performed the statistical analysis. JS and LYD contributed to the data analysis. YLG contributed to improving the writing of manuscripts. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Cui, Z., Gao, YL., Liu, JX. et al. The computational prediction of drugdisease interactions using the dualnetwork L_{2,1}CMF method. BMC Bioinformatics 20, 5 (2019). https://doi.org/10.1186/s1285901825756
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1285901825756
Keywords
 Drugdisease interactions
 L_{2,1}norm
 Gaussian interaction profile
 Matrix factorization