How to Illuminate the Druggable Genome Using Pharos

Proteome-scale characterisation of motif-based interactome rewiring by disease mutations

Whole genome and exome sequencing are reporting on hundreds of thousands of missense mutations. Taking a pan-disease approach, we explored how mutations in intrinsically disordered regions (IDRs) break or generate protein interactions mediated by short linear motifs. We created a peptide-phage display library tiling ~57,000 peptides from the IDRs of the human proteome overlapping 12,301 single nucleotide variants associated with diverse phenotypes including cancer, metabolic diseases and neurological diseases. By screening 80 human proteins, we identified 366 mutation-modulated interactions, with half of the mutations diminishing binding, and half enhancing binding or creating novel interaction interfaces. The effects of the mutations were confirmed by affinity measurements. In cellular assays, the effects of motif-disruptive mutations were validated, including loss of a nuclear localisation signal in the cell division control protein CDC45 by a mutation associated with Meier-Gorlin syndrome. The study provides insights into how disease-associated mutations may perturb and rewire the motif-based interactome.

Distinguishing heart failure with reduced ejection fraction from heart failure with preserved ejection fraction: A phenomics approach

AIM

Pathophysiological differences between patients with heart failure with preserved (HFpEF) and reduced (HFrEF) ejection fraction (EF) remain unclear. Therefore we used a phenomics approach, integrating selected proteomics data with patient characteristics and cardiac structural and functional parameters, to get insight into differential pathophysiological mechanisms and identify potential treatment targets.

METHODS AND RESULTS

We report data from a representative subcohort of the prospective Singapore Heart Failure Outcomes and Phenotypes (SHOP), including patients with HFrEF (EF <40%, n = 217), HFpEF (EF ≥50%, n = 213), and age- and sex-matched controls without HF (n = 216). We measured 92 biomarkers using a proximity extension assay and assessed cardiac structure and function in all participants using echocardiography. We used multi-block projection to latent structure analysis to integrate clinical, echocardiographic, and biomarker variables. Candidate biomarker targets were cross-referenced with small-molecule and drug databases. The total cohort had a median age of 65 years (interquartile range 60-71), and 50% were women. Protein profiles strongly discriminated patients with HFrEF (area under the curve [AUC] = 0.89) and HFpEF (AUC = 0.94) from controls. Phenomics analyses identified unique druggable inflammatory markers in HFpEF from the tumour necrosis factor receptor superfamily (TNFRSF), which were positively associated with hypertension, diabetes, and increased posterior and relative wall thickness. In HFrEF, interleukin (IL)-8 and IL-6 were possible targets related to lower EF and worsening renal function.

CONCLUSION

We identified pathophysiological mechanisms related to increased cardiac wall thickness parameters and potentially druggable inflammatory markers from the TNFRSF in HFpEF.

Informatic challenges and advances in illuminating the druggable proteome

The understudied members of the druggable proteomes offer promising prospects for drug discovery efforts. While large-scale initiatives have generated valuable functional information on understudied members of the druggable gene families, translating this information into actionable knowledge for drug discovery requires specialized informatics tools and resources. Here, we review the unique informatics challenges and advances in annotating understudied members of the druggable proteome. We demonstrate the application of statistical evolutionary inference tools, knowledge graph mining approaches, and protein language models in illuminating understudied protein kinases, pseudokinases, and ion channels.

FAVA: high-quality functional association networks inferred from scRNA-seq and proteomics data

MOTIVATION

Protein networks are commonly used for understanding how proteins interact. However, they are typically biased by data availability, favoring well-studied proteins with more interactions. To uncover functions of understudied proteins, we must use data that are not affected by this literature bias, such as single-cell RNA-seq and proteomics. Due to data sparseness and redundancy, functional association analysis becomes complex.

RESULTS

To address this, we have developed FAVA (Functional Associations using Variational Autoencoders), which compresses high-dimensional data into a low-dimensional space. FAVA infers networks from high-dimensional omics data with much higher accuracy than existing methods, across a diverse collection of real as well as simulated datasets. FAVA can process large datasets with over 0.5 million conditions and has predicted 4210 interactions between 1039 understudied proteins. Our findings showcase FAVA's capability to offer novel perspectives on protein interactions. FAVA functions within the scverse ecosystem, employing AnnData as its input source.

AVAILABILITY AND IMPLEMENTATION

Source code, documentation, and tutorials for FAVA are accessible on GitHub at https://github.com/mikelkou/fava. FAVA can also be installed and used via pip/PyPI as well as via the scverse ecosystem https://github.com/scverse/ecosystem-packages/tree/main/packages/favapy.

Advancing Biomarker Discovery and Therapeutic Targets in Duchenne Muscular Dystrophy: A Comprehensive Review

Mounting evidence underscores the intricate interplay between the immune system and skeletal muscles in Duchenne muscular dystrophy (DMD), as well as during regular muscle regeneration. While immune cell infiltration into skeletal muscles stands out as a prominent feature in the disease pathophysiology, a myriad of secondary defects involving metabolic and inflammatory pathways persist, with the key players yet to be fully elucidated. Steroids, currently the sole effective therapy for delaying onset and symptom control, come with adverse side effects, limiting their widespread use. Preliminary evidence spotlighting the distinctive features of T cell profiling in DMD prompts the immuno-characterization of circulating cells. A molecular analysis of their transcriptome and secretome holds the promise of identifying a subpopulation of cells suitable as disease biomarkers. Furthermore, it provides a gateway to unraveling new pathological pathways and pinpointing potential therapeutic targets. Simultaneously, the last decade has witnessed the emergence of novel approaches. The development and equilibrium of both innate and adaptive immune systems are intricately linked to the gut microbiota. Modulating microbiota-derived metabolites could potentially exacerbate muscle damage through immune system activation. Concurrently, genome sequencing has conferred clinical utility for rare disease diagnosis since innovative methodologies have been deployed to interpret the functional consequences of genomic variations. Despite numerous genes falling short as clinical targets for MD, the exploration of Tdark genes holds promise for unearthing novel and uncharted therapeutic insights. In the quest to expedite the translation of fundamental knowledge into clinical applications, the identification of novel biomarkers and disease targets is paramount. This initiative not only advances our understanding but also paves the way for the design of innovative therapeutic strategies, contributing to enhanced care for individuals grappling with these incapacitating diseases.

Dark kinase annotation, mining, and visualization using the Protein Kinase Ontology

The Protein Kinase Ontology (ProKinO) is an integrated knowledge graph that conceptualizes the complex relationships among protein kinase sequence, structure, function, and disease in a human and machine-readable format. In this study, we have significantly expanded ProKinO by incorporating additional data on expression patterns and drug interactions. Furthermore, we have developed a completely new browser from the ground up to render the knowledge graph visible and interactive on the web. We have enriched ProKinO with new classes and relationships that capture information on kinase ligand binding sites, expression patterns, and functional features. These additions extend ProKinO's capabilities as a discovery tool, enabling it to uncover novel insights about understudied members of the protein kinase family. We next demonstrate the application of ProKinO. Specifically, through graph mining and aggregate SPARQL queries, we identify the p21-activated protein kinase 5 (PAK5) as one of the most frequently mutated dark kinases in human cancers with abnormal expression in multiple cancers, including a previously unappreciated role in acute myeloid leukemia. We have identified recurrent oncogenic mutations in the PAK5 activation loop predicted to alter substrate binding and phosphorylation. Additionally, we have identified common ligand/drug binding residues in PAK family kinases, underscoring ProKinO's potential application in drug discovery. The updated ontology browser and the addition of a web component, ProtVista, which enables interactive mining of kinase sequence annotations in 3D structures and Alphafold models, provide a valuable resource for the signaling community. The updated ProKinO database is accessible at https://prokino.uga.edu.

Exploring DrugCentral: from molecular structures to clinical effects

DrugCentral, accessible at https://drugcentral.org , is an open-access online drug information repository. It covers over 4950 drugs, incorporating structural, physicochemical, and pharmacological details to support drug discovery, development, and repositioning. With around 20,000 bioactivity data points, manual curation enhances information from several major digital sources. Approximately 724 mechanism-of-action (MoA) targets offer updated drug target insights. The platform captures clinical data: over 14,300 on- and off-label uses, 27,000 contraindications, and around 340,000 adverse drug events from pharmacovigilance reports. DrugCentral encompasses information from molecular structures to marketed formulations, providing a comprehensive pharmaceutical reference. Users can easily navigate basic drug information and key features, making DrugCentral a versatile, unique resource. Furthermore, we present a use-case example where we utilize experimentally determined data from DrugCentral to support drug repurposing. A minimum activity threshold t should be considered against novel targets to repurpose a drug. Analyzing 1156 bioactivities for human MoA targets suggests a general threshold of 1 µM: t = 6 when expressed as - log[Activity(M)]). This applies to 87% of the drugs. Moreover, t can be refined empirically based on water solubility (S): t = 3 - logS, for logS < - 3. Alongside the drug repurposing classification scheme, which considers intellectual property rights, market exclusivity protections, and market accessibility, DrugCentral provides valuable data to prioritize candidates for drug repurposing programs efficiently.

Identification of drug repurposing candidates for the treatment of anxiety: A genetic approach

Anxiety disorders are a group of prevalent and heritable neuropsychiatric diseases. We previously conducted a genome-wide association study (GWAS) which identified genomic loci associated with anxiety; however, the biological consequences underlying the genetic associations are largely unknown. Integrating GWAS and functional genomic data may improve our understanding of the genetic effects on intermediate molecular phenotypes such as gene expression. This can provide an opportunity for the discovery of drug targets for anxiety via drug repurposing. We used the GWAS summary statistics to determine putative causal genes for anxiety using MAGMA and colocalization analyses. A transcriptome-wide association study was conducted to identify genes with differential genetically regulated levels of gene expression in human brain tissue. The genes were integrated with a large drug-gene expression database (Connectivity Map), discovering compounds that are predicted to "normalise" anxiety-associated expression changes. The study identified 64 putative causal genes associated with anxiety (35 genes upregulated; 29 genes downregulated). Drug mechanisms adrenergic receptor agonists, sigma receptor agonists, and glutamate receptor agonists gene targets were enriched in anxiety-associated genetic signal and exhibited an opposing effect on the anxiety-associated gene expression signature. The significance of the project demonstrated genetic links for novel drug candidates to potentially advance anxiety therapeutics.

Illuminate the Functions of Dark Proteins Using the Reactome-IDG Web Portal

Understudied or dark proteins have the potential to shed light on as-yet undiscovered molecular mechanisms that underlie phenotypes and suggest innovative therapeutic approaches for many diseases. The Reactome-IDG (Illuminating the Druggable Genome) project aims to place dark proteins in the context of manually curated, highly reliable pathways in Reactome, the most comprehensive, open-source biological pathway knowledgebase, facilitating the understanding functions and predicting therapeutic potentials of dark proteins. The Reactome-IDG web portal, deployed at https://idg.reactome.org, provides a simple, interactive web page for users to search pathways that may functionally interact with dark proteins, enabling the prediction of functions of dark proteins in the context of Reactome pathways. Enhanced visualization features implemented at the portal allow users to investigate the functional contexts for dark proteins based on tissue-specific gene or protein expression, drug-target interactions, or protein or gene pairwise relationships in the original Reactome's systems biology graph notation (SBGN) diagrams or the new simplified functional interaction (FI) network view of pathways. The protocols in this chapter describe step-by-step procedures to use the web portal to learn biological functions of dark proteins in the context of Reactome pathways. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Search for interacting pathways of a protein Support Protocol: Interacting pathway results for an annotated protein Alternate Protocol: Use individual pairwise relationships to predict interacting pathways of a protein Basic Protocol 2: Using the IDG pathway browser to study interacting pathways Basic Protocol 3: Overlaying tissue-specific expression data Basic Protocol 4: Overlaying protein/gene pairwise relationships in the pathway context Basic Protocol 5: Visualizing drug/target interactions.

Consensus molecular environment of schizophrenia risk genes in coexpression networks shifting across age and brain regions

Schizophrenia is a neurodevelopmental brain disorder whose genetic risk is associated with shifting clinical phenomena across the life span. We investigated the convergence of putative schizophrenia risk genes in brain coexpression networks in postmortem human prefrontal cortex (DLPFC), hippocampus, caudate nucleus, and dentate gyrus granule cells, parsed by specific age periods (total N = 833). The results support an early prefrontal involvement in the biology underlying schizophrenia and reveal a dynamic interplay of regions in which age parsing explains more variance in schizophrenia risk compared to lumping all age periods together. Across multiple data sources and publications, we identify 28 genes that are the most consistently found partners in modules enriched for schizophrenia risk genes in DLPFC; twenty-three are previously unidentified associations with schizophrenia. In iPSC-derived neurons, the relationship of these genes with schizophrenia risk genes is maintained. The genetic architecture of schizophrenia is embedded in shifting coexpression patterns across brain regions and time, potentially underwriting its shifting clinical presentation.

Artificial Intelligence and Machine Learning for Lead-to-Candidate Decision-Making and Beyond

The use of artificial intelligence (AI) and machine learning (ML) in pharmaceutical research and development has to date focused on research: target identification; docking-, fragment-, and motif-based generation of compound libraries; modeling of synthesis feasibility; rank-ordering likely hits according to structural and chemometric similarity to compounds having known activity and affinity to the target(s); optimizing a smaller library for synthesis and high-throughput screening; and combining evidence from screening to support hit-to-lead decisions. Applying AI/ML methods to lead optimization and lead-to-candidate (L2C) decision-making has shown slower progress, especially regarding predicting absorption, distribution, metabolism, excretion, and toxicology properties. The present review surveys reasons why this is so, reports progress that has occurred in recent years, and summarizes some of the issues that remain. Effective AI/ML tools to derisk L2C and later phases of development are important to accelerate the pharmaceutical development process, ameliorate escalating development costs, and achieve greater success rates.

The International Mouse Phenotyping Consortium: comprehensive knockout phenotyping underpinning the study of human disease

The International Mouse Phenotyping Consortium (IMPC; https://www.mousephenotype.org/) web portal makes available curated, integrated and analysed knockout mouse phenotyping data generated by the IMPC project consisting of 85M data points and over 95,000 statistically significant phenotype hits mapped to human diseases. The IMPC portal delivers a substantial reference dataset that supports the enrichment of various domain-specific projects and databases, as well as the wider research and clinical community, where the IMPC genotype-phenotype knowledge contributes to the molecular diagnosis of patients affected by rare disorders. Data from 9,000 mouse lines and 750 000 images provides vital resources enabling the interpretation of the ignorome, and advancing our knowledge on mammalian gene function and the mechanisms underlying phenotypes associated with human diseases. The resource is widely integrated and the lines have been used in over 4,600 publications indicating the value of the data and the materials.

Pharos 2023: an integrated resource for the understudied human proteome

The Illuminating the Druggable Genome (IDG) project aims to improve our understanding of understudied proteins and our ability to study them in the context of disease biology by perturbing them with small molecules, biologics, or other therapeutic modalities. Two main products from the IDG effort are the Target Central Resource Database (TCRD) (http://juniper.health.unm.edu/tcrd/), which curates and aggregates information, and Pharos (https://pharos.nih.gov/), a web interface for fusers to extract and visualize data from TCRD. Since the 2021 release, TCRD/Pharos has focused on developing visualization and analysis tools that help reveal higher-level patterns in the underlying data. The current iterations of TCRD and Pharos enable users to perform enrichment calculations based on subsets of targets, diseases, or ligands and to create interactive heat maps and UpSet charts of many types of annotations. Using several examples, we show how to address disease biology and drug discovery questions through enrichment calculations and UpSet charts.

Human oncogenic viruses: an overview of protein biomarkers in viral cancers and their potential use in clinics

INTRODUCTION

Although the idea that carcinogenesis might be caused by viruses was first voiced about 100 years ago, today's data disappointingly show that we have not made much progress in preventing and/or treating viral cancers in a century. According to recent studies, infections are responsible for approximately 13% of cancer development in the world. Today, it is accepted and proven by many authorities that Epstein-Barr virus (EBV), Hepatitis B virus (HBV), Hepatitis C virus (HCV), Human Herpesvirus 8 (HHV8), Human T-cell Lymphotropic virus 1 (HTLV1) and highly oncogenic Human Papillomaviruses (HPVs) cause or/and contribute to cancer development in humans.

AREAS COVERED

Considering the insufficient prevention and/or treatment strategies for viral cancers, in this review we present the current knowledge on protein biomarkers of oncogenic viruses. In addition, we aimed to decipher their potential for clinical use by evaluating whether the proposed biomarkers are expressed in body fluids, are druggable, and act as tumor suppressors or oncoproteins.

EXPERT OPINION

Consequently, we believe that this review will shed light on researchers and provide a guide to find remarkable solutions for the prevention and/or treatment of viral cancers.

Inhibition of CDK8/19 Mediator kinase potentiates HER2-targeting drugs and bypasses resistance to these agents in vitro and in vivo

Breast cancers (BrCas) that overexpress oncogenic tyrosine kinase receptor HER2 are treated with HER2-targeting antibodies (such as trastuzumab) or small-molecule kinase inhibitors (such as lapatinib). However, most patients with metastatic HER2+ BrCa have intrinsic resistance and nearly all eventually become resistant to HER2-targeting therapy. Resistance to HER2-targeting drugs frequently involves transcriptional reprogramming associated with constitutive activation of different signaling pathways. We have investigated the role of CDK8/19 Mediator kinase, a regulator of transcriptional reprogramming, in the response of HER2+ BrCa to HER2-targeting drugs. CDK8 was in the top 1% of all genes ranked by correlation with shorter relapse-free survival among treated HER2+ BrCa patients. Selective CDK8/19 inhibitors (senexin B and SNX631) showed synergistic interactions with lapatinib and trastuzumab in a panel of HER2+ BrCa cell lines, overcoming and preventing resistance to HER2-targeting drugs. The synergistic effects were mediated in part through the PI3K/AKT/mTOR pathway and reduced by PI3K inhibition. Combination of HER2- and CDK8/19-targeting agents inhibited STAT1 and STAT3 phosphorylation at S727 and up-regulated tumor suppressor BTG2. The growth of xenograft tumors formed by lapatinib-sensitive or -resistant HER2+ breast cancer cells was partially inhibited by SNX631 alone and strongly suppressed by the combination of SNX631 and lapatinib, overcoming lapatinib resistance. These effects were associated with decreased tumor cell proliferation and altered recruitment of stromal components to the xenograft tumors. These results suggest potential clinical benefit of combining HER2- and CDK8/19-targeting drugs in the treatment of metastatic HER2+ BrCa.

Identification of novel γ-globin inducers among all potential erythroid druggable targets

Human γ-globin is predominantly expressed in fetal liver erythroid cells during gestation from 2 nearly identical genes, HBG1 and HBG2, that are both perinatally silenced. Reactivation of these fetal genes in adult red blood cells can ameliorate many symptoms associated with the inherited β-globinopathies, sickle cell disease, and Cooley anemia. Although promising genetic strategies to reactivate the γ-globin genes to treat these diseases have been explored, there are significant barriers to their effective implementation worldwide; alternatively, pharmacological induction of γ-globin synthesis could readily reach the majority of affected individuals. In this study, we generated a CRISPR knockout library that targeted all erythroid genes for which prospective or actual therapeutic compounds already exist. By probing this library for genes that repress fetal hemoglobin (HbF), we identified several novel, potentially druggable, γ-globin repressors, including VHL and PTEN. We demonstrate that deletion of VHL induces HbF through activation of the HIF1α pathway and that deletion of PTEN induces HbF through AKT pathway stimulation. Finally, we show that small-molecule inhibitors of PTEN and EZH induce HbF in both healthy and β-thalassemic human primary erythroid cells.

Considerations and challenges for sex-aware drug repurposing

Sex differences are essential factors in disease etiology and manifestation in many diseases such as cardiovascular disease, cancer, and neurodegeneration [33]. The biological influence of sex differences (including genomic, epigenetic, hormonal, immunological, and metabolic differences between males and females) and the lack of biomedical studies considering sex differences in their study design has led to several policies. For example, the National Institute of Health's (NIH) sex as a biological variable (SABV) and Sex and Gender Equity in Research (SAGER) policies to motivate researchers to consider sex differences [204]. However, drug repurposing, a promising alternative to traditional drug discovery by identifying novel uses for FDA-approved drugs, lacks sex-aware methods that can improve the identification of drugs that have sex-specific responses [7, 11, 14, 33]. Sex-aware drug repurposing methods either select drug candidates that are more efficacious in one sex or deprioritize drug candidates based on if they are predicted to cause a sex-bias adverse event (SBAE), unintended therapeutic effects that are more likely to occur in one sex. Computational drug repurposing methods are encouraging approaches to develop for sex-aware drug repurposing because they can prioritize sex-specific drug candidates or SBAEs at lower cost and time than traditional drug discovery. Sex-aware methods currently exist for clinical, genomic, and transcriptomic information [1, 7, 155]. They have not expanded to other data types, such as DNA variation, which has been beneficial in other drug repurposing methods that do not consider sex [114]. Additionally, some sex-aware methods suffer from poorer performance because a disproportionate number of male and female samples are available to train computational methods [7]. However, there is development potential for several different categories (i.e., data mining, ligand binding predictions, molecular associations, and networks). Low-dimensional representations of molecular association and network approaches are also especially promising candidates for future sex-aware drug repurposing methodologies because they reduce the multiple hypothesis testing burden and capture sex-specific variation better than the other methods [151, 159]. Here we review how sex influences drug response, the current state of drug repurposing including with respect to sex-bias drug response, and how model organism study design choices influence drug repurposing validation.

Path to Clonal Theranostics in Luminal Breast Cancers

Integrating tumor heterogeneity in the drug discovery process is a key challenge to tackle breast cancer resistance. Identifying protein targets for functionally distinct tumor clones is particularly important to tailor therapy to the heterogeneous tumor subpopulations and achieve clonal theranostics. For this purpose, we performed an unsupervised, label-free, spatially resolved shotgun proteomics guided by MALDI mass spectrometry imaging (MSI) on 124 selected tumor clonal areas from early luminal breast cancers, tumor stroma, and breast cancer metastases. 2868 proteins were identified. The main protein classes found in the clonal proteome dataset were enzymes, cytoskeletal proteins, membrane-traffic, translational or scaffold proteins, or transporters. As a comparison, gene-specific transcriptional regulators, chromatin related proteins or transmembrane signal receptor were more abundant in the TCGA dataset. Moreover, 26 mutated proteins have been identified. Similarly, expanding the search to alternative proteins databases retrieved 126 alternative proteins in the clonal proteome dataset. Most of these alternative proteins were coded mainly from non-coding RNA. To fully understand the molecular information brought by our approach and its relevance to drug target discovery, the clonal proteomic dataset was further compared to the TCGA breast cancer database and two transcriptomic panels, BC360 (nanoString^®) and CDx (Foundation One^®). We retrieved 139 pathways in the clonal proteome dataset. Only 55% of these pathways were also present in the TCGA dataset, 68% in BC360 and 50% in CDx. Seven of these pathways have been suggested as candidate for drug targeting, 22 have been associated with breast cancer in experimental or clinical reports, the remaining 19 pathways have been understudied in breast cancer. Among the anticancer drugs, 35 drugs matched uniquely with the clonal proteome dataset, with only 7 of them already approved in breast cancer. The number of target and drug interactions with non-anticancer drugs (such as agents targeting the cardiovascular system, metabolism, the musculoskeletal or the nervous systems) was higher in the clonal proteome dataset (540 interactions) compared to TCGA (83 interactions), BC360 (419 interactions), or CDx (172 interactions). Many of the protein targets identified and drugs screened were clinically relevant to breast cancer and are in clinical trials. Thus, we described the non-redundant knowledge brought by this clone-tailored approach compared to TCGA or transcriptomic panels, the targetable proteins identified in the clonal proteome dataset, and the potential of this approach for drug discovery and repurposing through drug interactions with antineoplastic agents and non-anticancer drugs.

Target 2035 - update on the quest for a probe for every protein

Twenty years after the publication of the first draft of the human genome, our knowledge of the human proteome is still fragmented. The challenge of translating the wealth of new knowledge from genomics into new medicines is that proteins, and not genes, are the primary executers of biological function. Therefore, much of how biology works in health and disease must be understood through the lens of protein function. Accordingly, a subset of human proteins has been at the heart of research interests of scientists over the centuries, and we have accumulated varying degrees of knowledge about approximately 65% of the human proteome. Nevertheless, a large proportion of proteins in the human proteome (∼35%) remains uncharacterized, and less than 5% of the human proteome has been successfully targeted for drug discovery. This highlights the profound disconnect between our abilities to obtain genetic information and subsequent development of effective medicines. Target 2035 is an international federation of biomedical scientists from the public and private sectors, which aims to address this gap by developing and applying new technologies to create by year 2035 chemogenomic libraries, chemical probes, and/or biological probes for the entire human proteome.

Knowledge graph analytics platform with LINCS and IDG for Parkinson's disease target illumination

Background

LINCS, "Library of Integrated Network-based Cellular Signatures", and IDG, "Illuminating the Druggable Genome", are both NIH projects and consortia that have generated rich datasets for the study of the molecular basis of human health and disease. LINCS L1000 expression signatures provide unbiased systems/omics experimental evidence. IDG provides compiled and curated knowledge for illumination and prioritization of novel drug target hypotheses. Together, these resources can support a powerful new approach to identifying novel drug targets for complex diseases, such as Parkinson's disease (PD), which continues to inflict severe harm on human health, and resist traditional research approaches.

Results

Integrating LINCS and IDG, we built the Knowledge Graph Analytics Platform (KGAP) to support an important use case: identification and prioritization of drug target hypotheses for associated diseases. The KGAP approach includes strong semantics interpretable by domain scientists and a robust, high performance implementation of a graph database and related analytical methods. Illustrating the value of our approach, we investigated results from queries relevant to PD. Approved PD drug indications from IDG’s resource DrugCentral were used as starting points for evidence paths exploring chemogenomic space via LINCS expression signatures for associated genes, evaluated as target hypotheses by integration with IDG. The KG-analytic scoring function was validated against a gold standard dataset of genes associated with PD as elucidated, published mechanism-of-action drug targets, also from DrugCentral. IDG's resource TIN-X was used to rank and filter KGAP results for novel PD targets, and one, SYNGR3 (Synaptogyrin-3), was manually investigated further as a case study and plausible new drug target for PD.

Conclusions

The synergy of LINCS and IDG, via KG methods, empowers graph analytics methods for the investigation of the molecular basis of complex diseases, and specifically for identification and prioritization of novel drug targets. The KGAP approach enables downstream applications via integration with resources similarly aligned with modern KG methodology. The generality of the approach indicates that KGAP is applicable to many disease areas, in addition to PD, the focus of this paper.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04530-9.

Getting Started with the IDG KMC Datasets and Tools

The Illuminating the Druggable Genome (IDG) consortium is a National Institutes of Health (NIH) Common Fund program designed to enhance our knowledge of under-studied proteins, more specifically, proteins unannotated within the three most commonly drug-targeted protein families: G-protein coupled receptors, ion channels, and protein kinases. Since 2014, the IDG Knowledge Management Center (IDG-KMC) has generated several open-access datasets and resources that jointly serve as a highly translational machine-learning-ready knowledgebase focused on human protein-coding genes and their products. The goal of the IDG-KMC is to develop comprehensive integrated knowledge for the druggable genome to illuminate the uncharacterized or poorly annotated portion of the druggable genome. The tools derived from the IDG-KMC provide either user-friendly visualizations or ways to impute the knowledge about potential targets using machine learning strategies. In the following protocols, we describe how to use each web-based tool to accelerate illumination in under-studied proteins. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Interacting with the Pharos user interface Basic Protocol 2: Accessing the data in Harmonizome Basic Protocol 3: The ARCHS4 resource Basic Protocol 4: Making predictions about gene function with PrismExp Basic Protocol 5: Using Geneshot to illuminate knowledge about under-studied targets Basic Protocol 6: Exploring under-studied targets with TIN-X Basic Protocol 7: Interacting with the DrugCentral user interface Basic Protocol 8: Estimating Anti-SARS-CoV-2 activities with DrugCentral REDIAL-2020 Basic Protocol 9: Drug Set Enrichment Analysis using Drugmonizome Basic Protocol 10: The Drugmonizome-ML Appyter Basic Protocol 11: The Harmonizome-ML Appyter Basic Protocol 12: GWAS target illumination with TIGA Basic Protocol 13: Prioritizing kinases for lists of proteins and phosphoproteins with KEA3 Basic Protocol 14: Converting PubMed searches to drug sets with the DrugShot Appyter.

The Playbooks of Medicinal Chemistry Design Moves

Large databases of biologically relevant molecules, such as ChEMBL, SureChEMBL, or compound collections of pharmaceutical or agrochemical companies, are invaluable sources of medicinal chemistry information, albeit implicit. We developed a modified matched molecular pair approach to systematically and exhaustively extract the transformations in these databases and distill them into snippets of explicit design knowledge that are easily interpretable and directly applicable. The resulting "playbooks of medicinal chemistry design moves" capture the collective pharmaceutical and agrochemical research expertise across multiple chemists, companies, targets, and projects. They can be queried in an automated fashion for systematic prospective design and compound generation. The ChEMBL playbook and an application to exploit it are available at https://github.com/mahendra-awale/medchem_moves.

TCRD and Pharos 2021: mining the human proteome for disease biology

In 2014, the National Institutes of Health (NIH) initiated the Illuminating the Druggable Genome (IDG) program to identify and improve our understanding of poorly characterized proteins that can potentially be modulated using small molecules or biologics. Two resources produced from these efforts are: The Target Central Resource Database (TCRD) (http://juniper.health.unm.edu/tcrd/) and Pharos (https://pharos.nih.gov/), a web interface to browse the TCRD. The ultimate goal of these resources is to highlight and facilitate research into currently understudied proteins, by aggregating a multitude of data sources, and ranking targets based on the amount of data available, and presenting data in machine learning ready format. Since the 2017 release, both TCRD and Pharos have produced two major releases, which have incorporated or expanded an additional 25 data sources. Recently incorporated data types include human and viral-human protein-protein interactions, protein-disease and protein-phenotype associations, and drug-induced gene signatures, among others. These aggregated data have enabled us to generate new visualizations and content sections in Pharos, in order to empower users to find new areas of study in the druggable genome.