Phenome-based gene discovery provides information about Parkinson's disease drug targets

Promising biomarkers and therapeutic targets for the management of Parkinson's disease: recent advancements and contemporary research

Parkinson's disease (PD) is one of the progressive neurological diseases which affect around 10 million population worldwide. The clinical manifestation of motor symptoms in PD patients appears later when most dopaminergic neurons have degenerated. Thus, for better management of PD, the development of accurate biomarkers for the early prognosis of PD is imperative. The present work will discuss the potential biomarkers from various attributes covering biochemical, microRNA, and neuroimaging aspects (α-synuclein, DJ-1, UCH-L1, β-glucocerebrosidase, BDNF, etc.) for diagnosis, recent development in PD management, and major limitations with current and conventional anti-Parkinson therapy. This manuscript summarizes potential biomarkers and therapeutic targets, based on available preclinical and clinical evidence, for better management of PD.

Combining phenome-driven drug-target interaction prediction with patients' electronic health records-based clinical corroboration toward drug discovery

Motivation

Predicting drug–target interactions (DTIs) using human phenotypic data have the potential in eliminating the translational gap between animal experiments and clinical outcomes in humans. One challenge in human phenome-driven DTI predictions is integrating and modeling diverse drug and disease phenotypic relationships. Leveraging large amounts of clinical observed phenotypes of drugs and diseases and electronic health records (EHRs) of 72 million patients, we developed a novel integrated computational drug discovery approach by seamlessly combining DTI prediction and clinical corroboration.

Results

We developed a network-based DTI prediction system (TargetPredict) by modeling 855 904 phenotypic and genetic relationships among 1430 drugs, 4251 side effects, 1059 diseases and 17 860 genes. We systematically evaluated TargetPredict in de novo cross-validation and compared it to a state-of-the-art phenome-driven DTI prediction approach. We applied TargetPredict in identifying novel repositioned candidate drugs for Alzheimer’s disease (AD), a disease affecting over 5.8 million people in the United States. We evaluated the clinical efficiency of top repositioned drug candidates using EHRs of over 72 million patients. The area under the receiver operating characteristic (ROC) curve was 0.97 in the de novo cross-validation when evaluated using 910 drugs. TargetPredict outperformed a state-of-the-art phenome-driven DTI prediction system as measured by precision–recall curves [measured by average precision (MAP): 0.28 versus 0.23, P-value < 0.0001]. The EHR-based case–control studies identified that the prescriptions top-ranked repositioned drugs are significantly associated with lower odds of AD diagnosis. For example, we showed that the prescription of liraglutide, a type 2 diabetes drug, is significantly associated with decreased risk of AD diagnosis [adjusted odds ratios (AORs): 0.76; 95% confidence intervals (CI) (0.70, 0.82), P-value < 0.0001]. In summary, our integrated approach that seamlessly combines computational DTI prediction and large-scale patients’ EHRs-based clinical corroboration has high potential in rapidly identifying novel drug targets and drug candidates for complex diseases.

Availability and implementation

nlp.case.edu/public/data/TargetPredict.

A Drug-Side Effect Context-Sensitive Network approach for drug target prediction

SUMMARY

Computational drug target prediction has become an important process in drug discovery. Network-based approaches are commonly used in computational drug-target interaction (DTI) prediction. Existing network-based approaches are limited in capturing the contextual information on how diseases, drugs and genes are connected. Here, we proposed a context-sensitive network (CSN) model for DTI prediction by modeling contextual drug phenotypic relationships. We constructed a Drug-Side Effect Context-Sensitive Network (DSE-CSN) of 139 760 drug-side effect pairs, representing 1480 drugs and 5868 side effects. We also built a protein-protein interaction network (PPIN) of 15 267 gene nodes and 178 972 weighted edges. A heterogeneous network was built by connecting the DSE-CSN and the PPIN through 3684 known DTIs. For each drug on the DSE-CSN, its genetic targets were predicted and prioritized using a network-based ranking algorithm. Our approach was evaluated in both de novo and leave-one-out cross-validation analysis using known DTIs as the gold standard. We compared our DSE-CSN-based model to the traditional similarity-based network (SBN)-based prediction model. The results suggested that the DSE-CSN-based model was able to rank known DTIs highly. In a de novo cross-validation, the area under the receiver operating characteristic (ROC) curve was 0.95. In a leave-one-out cross-validation, the average rank was top 3.2% for known DTIs. When it was compared to the SBN-based model using the Precision-Recall curve, our CSN-based model achieved a higher mean average precision (MAP) (0.23 versus 0.19, P-value<1e-4) in a de novo cross-validation analysis. We further improved the CSN-based DTI prediction by differentially weighting the drug-side effect pairs on the network and showed a significant improvement of the MAP (0.29 versus 0.23, P-value<1e-4). We also showed that the CSN-based model consistently achieved better performances than the traditional SBN-based model across different drug classes. Moreover, we demonstrated that our novel DTI predictions can be supported by published literature. In summary, the CSN-based model, by modeling the context-specific inter-relationships among drugs and side effects, has a high potential in drug target prediction.

AVAILABILITY AND IMPLEMENTATION

nlp/case/edu/public/data/DSE/CSN_DTI.

Predicting links between tumor samples and genes using 2-Layered graph based diffusion approach

Background

Determining the association between tumor sample and the gene is demanding because it requires a high cost for conducting genetic experiments. Thus, the discovered association between tumor sample and gene further requires clinical verification and validation. This entire mechanism is time-consuming and expensive. Due to this issue, predicting the association between tumor samples and genes remain a challenge in biomedicine.

Results

Here we present, a computational model based on a heat diffusion algorithm which can predict the association between tumor samples and genes. We proposed a 2-layered graph. In the first layer, we constructed a graph of tumor samples and genes where these two types of nodes are connected by “hasGene” relationship. In the second layer, the gene nodes are connected by “interaction” relationship. We applied the heat diffusion algorithms in nine different variants of genetic interaction networks extracted from STRING and BioGRID database. The heat diffusion algorithm predicted the links between tumor samples and genes with mean AUC-ROC score of 0.84. This score is obtained by using weighted genetic interactions of fusion or co-occurrence channels from the STRING database. For the unweighted genetic interaction from the BioGRID database, the algorithms predict the links with an AUC-ROC score of 0.74.

Conclusions

We demonstrate that the gene-gene interaction scores could improve the predictive power of the heat diffusion model to predict the links between tumor samples and genes. We showed the efficient runtime of the heat diffusion algorithm in various genetic interaction network. We statistically validated our prediction quality of the links between tumor samples and genes.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3056-2) contains supplementary material, which is available to authorized users.

The Alzheimer's comorbidity phenome: mining from a large patient database and phenome-driven genetics prediction

Objective

Alzheimer’s disease (AD) is a severe neurodegenerative disorder and has become a global public health problem. Intensive research has been conducted for AD. But the pathophysiology of AD is still not elucidated. Disease comorbidity often associates diseases with overlapping patterns of genetic markers. This may inform a common etiology and suggest essential protein targets. US Food and Drug Administration (FDA) Adverse Event Reporting System (FAERS) collects large-scale postmarketing surveillance data that provide a unique opportunity to investigate disease co-occurrence pattern. We aim to construct a heterogeneous network that integrates disease comorbidity network (DCN) from FAERS with protein–protein interaction (PPI) to prioritize the AD risk genes using network-based ranking algorithm.

Materials and Methods

We built a DCN based on indication data from FAERS using association rule mining. DCN was further integrated with PPI network. We used random walk with restart ranking algorithm to prioritize AD risk genes.

Results

We evaluated the performance of our approach using AD risk genes curated from genetic association studies. Our approach achieved an area under a receiver operating characteristic curve of 0.770. Top 500 ranked genes achieved 5.53-fold enrichment for known AD risk genes as compared to random expectation. Pathway enrichment analysis using top-ranked genes revealed that two novel pathways, ERBB and coagulation pathways, might be involved in AD pathogenesis.

Conclusion

We innovatively leveraged FAERS, a comprehensive data resource for FDA postmarket drug safety surveillance, for large-scale AD comorbidity mining. This exploratory study demonstrated the potential of disease-comorbidities mining from FAERS in AD genetics discovery.

Combining mechanism-based prediction with patient-based profiling for psoriasis metabolomics biomarker discovery

Psoriasis is a chronic, debilitating skin condition that affects approximately 125 million individuals worldwide. The cause of psoriasis appears multifactorial, and no unified mitigating signal or single antigenic target has been identified to date. Metabolomic studies hold great potential for explaining disease mechanism, facilitating early diagnosis, and identifying potential therapeutic areas. Here, we present an integrated disease metabolomic biomarker discovery strategy that combines mechanism-based biomarker discovery with clinical sample-based metabolomic profiling. We applied this strategy in identifying and understanding metabolite biomarkers for psoriasis. The key innovation of our strategy is a novel mechanism-based metabolite prediction system, mmPredict, which assimilates vast amounts of existing knowledge of diseases and metabolites. mmPredict first constructed a psoriasis-specific mouse mutational phenotype profile. It then constructed phenotype profiles for a total of 259,170 chemicals/metabolites using known chemical genetics and human metabolomic data. Metabolites were then prioritized based on the phenotypic similarities between disease- and metabolites. We evaluated mmPredict using 150 metabolites identified using our in-house metabolome profiling study of psoriasis patient samples. mmPredict found 96 of the 150 metabolites and ranked them highly (recall: 0.64, mean ranking: 8.73%, median ranking: 2.33%, p-value: 4.75E-44). These results show that mmPredict is consistent with, as well as a complement to, traditional human metabolomic profiling studies. We then developed a strategy to combine outputs from both systems and found that the oxidative product of linoleic acid, 13(S)-hydroxy-9Z,11E-octadecadienoic acid (13- HODE), ranked highly by both mmPredict and our in-house experiments. Our integrated analysis indicates that 13- HODE may be a mechanistic link between psoriasis and cardiovascular comorbidities associated with psoriasis. In summary, we developed an integrated metabolomic prediction system that combines both human metabolomic studies and mechanism-based prediction and demonstrated its application in the skin disease psoriasis. Our system is highly general and can be applied to other diseases when patient-based metabolomic profiling data becomes more increasingly available. Data is publicly available at: http://nlp.

CASE

edu/public/data/mmPredict_PSO.

Using a novel computational drug-repositioning approach (DrugPredict) to rapidly identify potent drug candidates for cancer treatment

Computation-based drug-repurposing/repositioning approaches can greatly speed up the traditional drug discovery process. To date, systematic and comprehensive computation-based approaches to identify and validate drug-repositioning candidates for epithelial ovarian cancer (EOC) have not been undertaken. Here, we present a novel drug discovery strategy that combines a computational drug-repositioning system (DrugPredict) with biological testing in cell lines in order to rapidly identify novel drug candidates for EOC. DrugPredict exploited unique repositioning opportunities rendered by a vast amount of disease genomics, phenomics, drug treatment, and genetic pathway and uniquely revealed that non-steroidal anti-inflammatories (NSAIDs) rank just as high as currently used ovarian cancer drugs. As epidemiological studies have reported decreased incidence of ovarian cancer associated with regular intake of NSAIDs, we assessed whether NSAIDs could have chemoadjuvant applications in EOC and found that (i) NSAID Indomethacin induces robust cell death in primary patient-derived platinum-sensitive and platinum- resistant ovarian cancer cells and ovarian cancer stem cells and (ii) downregulation of β-catenin is partially driving effects of Indomethacin in cisplatin-resistant cells. In summary, we demonstrate that DrugPredict represents an innovative computational drug- discovery strategy to uncover drugs that are routinely used for other indications that could be effective in treating various cancers, thus introducing a potentially rapid and cost-effective translational opportunity. As NSAIDs are already in routine use in gynecological treatment regimens and have acceptable safety profile, our results will provide with a rationale for testing NSAIDs as potential chemoadjuvants in EOC patient trials.

Predict drug permeability to blood-brain-barrier from clinical phenotypes: drug side effects and drug indications

Motivation

Blood-Brain-Barrier (BBB) is a rigorous permeability barrier for maintaining homeostasis of Central Nervous System (CNS). Determination of compound's permeability to BBB is prerequisite in CNS drug discovery. Existing computational methods usually predict drug BBB permeability from chemical structure and they generally apply to small compounds passing BBB through passive diffusion. As abundant information on drug side effects and indications has been recorded over time through extensive clinical usage, we aim to explore BBB permeability prediction from a new angle and introduce a novel approach to predict BBB permeability from drug clinical phenotypes (drug side effects and drug indications). This method can apply to both small compounds and macro-molecules penetrating BBB through various mechanisms besides passive diffusion.

Results

We composed a training dataset of 213 drugs with known brain and blood steady-state concentrations ratio and extracted their side effects and indications as features. Next, we trained SVM models with polynomial kernel and obtained accuracy of 76.0%, AUC 0.739, and F 1 score (macro weighted) 0.760 with Monte Carlo cross validation. The independent test accuracy was 68.3%, AUC 0.692, F 1 score 0.676. When both chemical features and clinical phenotypes were available, combining the two types of features achieved significantly better performance than chemical feature based approach (accuracy 85.5% versus 72.9%, AUC 0.854 versus 0.733, F 1 score 0.854 versus 0.725; P  < e -90 ). We also conducted de novo prediction and identified 110 drugs in SIDER database having the potential to penetrate BBB, which could serve as start point for CNS drug repositioning research.

Availability and Implementation

https://github.com/bioinformatics-gao/CASE-BBB-prediction-Data.

Contact

rxx@case.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Parkinson's Disease: From Pathogenesis to Pharmacogenomics

Parkinson's disease (PD) is the second most important age-related neurodegenerative disorder in developed societies, after Alzheimer's disease, with a prevalence ranging from 41 per 100,000 in the fourth decade of life to over 1900 per 100,000 in people over 80 years of age. As a movement disorder, the PD phenotype is characterized by rigidity, resting tremor, and bradykinesia. Parkinson's disease -related neurodegeneration is likely to occur several decades before the onset of the motor symptoms. Potential risk factors include environmental toxins, drugs, pesticides, brain microtrauma, focal cerebrovascular damage, and genomic defects. Parkinson's disease neuropathology is characterized by a selective loss of dopaminergic neurons in the substantia nigra pars compacta, with widespread involvement of other central nervous system (CNS) structures and peripheral tissues. Pathogenic mechanisms associated with genomic, epigenetic and environmental factors lead to conformational changes and deposits of key proteins due to abnormalities in the ubiquitin-proteasome system together with dysregulation of mitochondrial function and oxidative stress. Conventional pharmacological treatments for PD are dopamine precursors (levodopa, l-DOPA, l-3,4 dihidroxifenilalanina), and other symptomatic treatments including dopamine agonists (amantadine, apomorphine, bromocriptine, cabergoline, lisuride, pergolide, pramipexole, ropinirole, rotigotine), monoamine oxidase (MAO) inhibitors (selegiline, rasagiline), and catechol-O-methyltransferase (COMT) inhibitors (entacapone, tolcapone). The chronic administration of antiparkinsonian drugs currently induces the "wearing-off phenomenon", with additional psychomotor and autonomic complications. In order to minimize these clinical complications, novel compounds have been developed. Novel drugs and bioproducts for the treatment of PD should address dopaminergic neuroprotection to reduce premature neurodegeneration in addition to enhancing dopaminergic neurotransmission. Since biochemical changes and therapeutic outcomes are highly dependent upon the genomic profiles of PD patients, personalized treatments should rely on pharmacogenetic procedures to optimize therapeutics.

Disease gene prioritization by integrating tissue-specific molecular networks using a robust multi-network model

BACKGROUND

Accurately prioritizing candidate disease genes is an important and challenging problem. Various network-based methods have been developed to predict potential disease genes by utilizing the disease similarity network and molecular networks such as protein interaction or gene co-expression networks. Although successful, a common limitation of the existing methods is that they assume all diseases share the same molecular network and a single generic molecular network is used to predict candidate genes for all diseases. However, different diseases tend to manifest in different tissues, and the molecular networks in different tissues are usually different. An ideal method should be able to incorporate tissue-specific molecular networks for different diseases.

RESULTS

In this paper, we develop a robust and flexible method to integrate tissue-specific molecular networks for disease gene prioritization. Our method allows each disease to have its own tissue-specific network(s). We formulate the problem of candidate gene prioritization as an optimization problem based on network propagation. When there are multiple tissue-specific networks available for a disease, our method can automatically infer the relative importance of each tissue-specific network. Thus it is robust to the noisy and incomplete network data. To solve the optimization problem, we develop fast algorithms which have linear time complexities in the number of nodes in the molecular networks. We also provide rigorous theoretical foundations for our algorithms in terms of their optimality and convergence properties. Extensive experimental results show that our method can significantly improve the accuracy of candidate gene prioritization compared with the state-of-the-art methods.

CONCLUSIONS

In our experiments, we compare our methods with 7 popular network-based disease gene prioritization algorithms on diseases from Online Mendelian Inheritance in Man (OMIM) database. The experimental results demonstrate that our methods recover true associations more accurately than other methods in terms of AUC values, and the performance differences are significant (with paired t-test p-values less than 0.05). This validates the importance to integrate tissue-specific molecular networks for studying disease gene prioritization and show the superiority of our network models and ranking algorithms toward this purpose. The source code and datasets are available at http://nijingchao.github.io/CRstar/ .