Metabolism is fundamental to sustaining human life, with changes in metabolite levels closely related to the occurrence and progression of diseases. The interaction between metabolites and drugs is intricate, encompassing drugs can modulate metabolite concentrations, as well as the metabolites generated through drug metabolism can influence pharmacological toxicity and drug interactions. Currently, a substantial proportion of metabolite-drug associations remains to be fully elucidated, and the data from individual sources are often characterized by incompleteness and noise. Here, we present BioWalk-MDA, a computational framework for large-scale predicting novel interactions between 8354 metabolites and 11 570 drugs. The framework constructs multilayered biomedical knowledge graphs (Multi-BiomedKGs) by integrating biological information across proteins, microbes, and diseases, and incorporated five types of graphs and seven types of associations. It employed random walk and heterogeneous Skip-gram model to extract feature vectors of metabolite-drug pairs and utilized a fully connected neural network (FCNN) to infer novel metabolite-drug associations. The framework demonstrated exceptional performance with an average accuracy of 0.971, an area under the receiver operating characteristic curve (AUROC) value of 0.995, and an area under the precision-recall curve (AUPRC) value of 0.994 in 5-fold cross-validation, surpassing other similar methods. Case studies on three metabolites detectable in blood and three cardiovascular drugs further demonstrated the reliability and efficiency of BioWalk-MDA, and it is anticipated to serve as a valuable tool for exploring metabolite-drug interactions and aiding in drug development and combination strategies.

MOTIVATION: miRNAs play crucial roles in human diseases and newly discovered could be targeted by small molecule (SM) drug compounds. Thus, the identification of small molecule drug compounds (SM) that target dysregulated miRNAs in cancers will provide new insight into cancer biology and accelerate drug discovery for cancer therapy. RESULTS: In this study, we aimed to develop a novel computational method to comprehensively identify associations between SMs and miRNAs. To this end, exploiting multiple molecular interaction databases, we first established an integrated SM-miRNA association network based on 690 561 SM to SM interactions, 291 600 miRNA to miRNA associations, as well as 664 known SM to miRNA targeting pairs. Then, by performing Random Walk with Restart algorithm on the integrated network, we prioritized the miRNAs associated to each of the SMs. By validating our results utilizing an independent dataset we obtained an area under the ROC curve greater than 0.7. Furthermore, comparisons indicated our integrated approach significantly improved the identification performance of those simple modeled methods. This computational framework as well as the prioritized SM-miRNA targeting relationships will promote the further developments of targeted cancer therapies. CONTACT: [email protected], [email protected] or [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The global insight into the relationships between miRNAs and their regulatory influences remains poorly understood. And most of complex diseases may be attributed to certain local areas of pathway (subpathway) instead of the entire pathway. Here, we reviewed the studies on miRNA regulations to pathways and constructed a bipartite miRNAs and subpathways network for systematic analyzing the miRNA regulatory influences to subpathways. We found that a small fraction of miRNAs were global regulators, environmental information processing pathways were preferentially regulated by miRNAs, and miRNAs had synergistic effect on regulating group of subpathways with similar function. Integrating the disease states of miRNAs, we also found that disease miRNAs regulated more subpathways than nondisease miRNAs, and for all miRNAs, the number of regulated subpathways was not in proportion to the number of the related diseases. Therefore, the study not only provided a global view on the relationships among disease, miRNA and subpathway, but also uncovered the function aspects of miRNA regulations and potential pathogenesis of complex diseases. A web server to query, visualize and download for all the data can be freely accessed at http://bioinfo.hrbmu.edu.cn/miR2Subpath.

MOTIVATION: Microarrays datasets frequently contain a large number of missing values (MVs), which need to be estimated and replaced for subsequent data mining. The focus of the paper is to study the effects of different MV treatments for cDNA microarray data on disease classification analysis. RESULTS: By analyzing five datasets, we demonstrate that among three kinds of classifiers evaluated in this study, support vector machine (SVM) classifiers are robust to varied MV imputation methods [e.g. replacing MVs by zero, K nearest-neighbor (KNN) imputation algorithm, local least square imputation and Bayesian principal component analysis], while the classification and regression tree classifiers are sensitive in terms of classification accuracy. The KNNclassifiers built on differentially expressed genes (DEGs) are robust to the varied MV treatments, but the performances of the KNN classifiers based on all measured genes can be significantly deteriorated when imputing MVs for genes with larger missing rate (MR) (e.g. MR > 5%). Generally, while replacing MVs by zero performs relatively poor, the other imputation algorithms have little difference in affecting classification performances of the SVM or KNN classifiers. We further demonstrate the power and feasibility of our recently proposed functional expression profile (FEP) approach as means to handle microarray data with MVs. The FEPs, which are derived from the functional modules that are enriched with sets of DEGs and thus can be consistently identified under varied MV treatments, achieve precise disease classification with better biological interpretation. We conclude that the choice of MV treatments should be determined in context of the later approaches used for disease classification. The suggested exclusion criterion of ignoring the genes with larger MR (e.g. >5%), while justifiable for some classifiers such as KNN classifiers, might not be considered as a general rule for all classifiers.