SIGNOR 3.0, the SIGnaling network open resource 3.0: 2022 update

Connecting the dots: Computational network analysis for disease insight and drug repurposing

Network biology is a powerful framework for studying the structure, function, and dynamics of biological systems, offering insights into the balance between health and disease states. The field is seeing rapid progress in all of its aspects: data availability, network synthesis, network analytics, and impactful applications in medicine and drug development. We review the most recent and significant results in network biomedicine, with a focus on the latest data, analytics, software resources, and applications in medicine. We also discuss what in our view are the likely directions of impactful development over the next few years.

Stage-dependent phosphoproteome remodeling of Parkinson's disease blood cells

The complexity and heterogeneity of PD necessitate advanced diagnostic and prognostic tools to elucidate its molecular mechanisms accurately. In this study, we addressed this challenge by conducting a pilot phospho-proteomic analysis of peripheral blood mononuclear cells (PBMCs) from idiopathic PD patients at varying disease stages to delineate the functional alterations occurring in these cells throughout the disease course and identify key molecules and pathways contributing to PD progression. By integrating clinical data with phospho-proteomic profiles across various PD stages, we identify potential stage-specific molecular signatures indicative of disease progression. This integrative approach allows for the discernment of distinct disease states and enhances our understanding of PD heterogeneity.

From CFTR to a CF signalling network: a systems biology approach to study Cystic Fibrosis

BACKGROUND

Cystic Fibrosis (CF) is a monogenic disease caused by mutations in the gene coding the Cystic Fibrosis Transmembrane Regulator (CFTR) protein, but its overall physio-pathology cannot be solely explained by the loss of the CFTR chloride channel function. Indeed, CFTR belongs to a yet not fully deciphered network of proteins participating in various signalling pathways.

METHODS

We propose a systems biology approach to study how the absence of the CFTR protein at the membrane leads to perturbation of these pathways, resulting in a panel of deleterious CF cellular phenotypes.

RESULTS

Based on publicly available transcriptomic datasets, we built and analyzed a CF network that recapitulates signalling dysregulations. The CF network topology and its resulting phenotypes were found to be consistent with CF pathology.

CONCLUSION

Analysis of the network topology highlighted a few proteins that may initiate the propagation of dysregulations, those that trigger CF cellular phenotypes, and suggested several candidate therapeutic targets. Although our research is focused on CF, the global approach proposed in the present paper could also be followed to study other rare monogenic diseases.

Computational investigation of the functional landscape of the protective role that extra virgin olive oil consumption may have on chronic lymphocytic leukemia

BACKGROUND

The health benefits of the Mediterranean diet are partially attributed to the polyphenols present in extra virgin olive oil (EVOO), which have been shown to have anti-cancer properties. However, the possible effect that EVOO could have on Chronic Lymphocytic Leukemia (CLL) has not been fully explored.

METHODS

 This study investigates the anti-CLL activity of EVOO through a computational multi-level data analysis procedure, focusing on the identification of shared biological functions between them. Specifically, publicly available data from genomics, transcriptomics and proteomics related to EVOO consumption and CLL were collected from several resources and analyzed through a computational pipeline, highlighting common molecular mechanisms and biological processes. Computational verification of a number of the highlighted functional terms associating CLL and EVOO has been performed as well.

RESULTS

Our investigation revealed four molecular pathways and three biological processes that overlap between mechanisms associated with CLL and those impacted by the consumption of EVOO. To further investigate the common biological functions, we focused on AKT1-related terms, aiming to investigate the potential importance of AKT1 in the anti- CLL effects associated with EVOO.

CONCLUSIONS

Overall, the results provide valuable insights into the potential beneficial effect of EVOO in CLL and highlight EVOO's bioactive compounds as promising candidates for future investigations.

A hybrid constrained continuous optimization approach for optimal causal discovery from biological data

MOTIVATION

Understanding causal effects is a fundamental goal of science and underpins our ability to make accurate predictions in unseen settings and conditions. While direct experimentation is the gold standard for measuring and validating causal effects, the field of causal graph theory offers a tantalizing alternative: extracting causal insights from observational data. Theoretical analysis has shown that this is indeed possible, given a large dataset and if certain conditions are met. However, biological datasets, frequently, do not meet such requirements but evaluation of causal discovery algorithms is typically performed on synthetic datasets, which they meet all requirements. Thus, real-life datasets are needed, in which the causal truth is reasonably known. In this work we first construct such a large-scale real-life dataset and then we perform on it a comprehensive benchmarking of various causal discovery methods.

RESULTS

We find that the PC algorithm is particularly accurate at estimating causal structure, including the causal direction which is critical for biological applicability. However, PC does only produces cause-effect directionality, but not estimates of causal effects. We propose PC-NOTEARS (PCnt), a hybrid solution, which includes the PC output as an additional constraint inside the NOTEARS optimization. This approach combines PC algorithm's strengths in graph structure prediction with the NOTEARS continuous optimization to estimate causal effects accurately. PCnt achieved best aggregate performance across all structural and effect size metrics.

AVAILABILITY AND IMPLEMENTATION

https://github.com/zhu-yh1/PC-NOTEARS.

SignalingProfiler 2.0 a network-based approach to bridge multi-omics data to phenotypic hallmarks

Unraveling how cellular signaling is remodeled upon perturbation is crucial for understanding disease mechanisms and identifying potential drug targets. In this pursuit, computational tools generating mechanistic hypotheses from multi-omics data have invaluable potential. Here, we present a newly implemented version (2.0) of SignalingProfiler, a multi-step pipeline to draw mechanistic hypotheses on the signaling events impacting cellular phenotypes. SignalingProfiler 2.0 derives context-specific signaling networks by integrating proteogenomic data with the prior knowledge-causal network. This is a freely accessible and flexible tool that incorporates statistical, footprint-based, and graph algorithms to accelerate the integration and interpretation of multi-omics data. Through a benchmarking process on three proof-of-concept studies, we demonstrate the tool's ability to generate hierarchical mechanistic networks recapitulating novel and known perturbed signaling and phenotypic outcomes, in both human and mice contexts. In summary, SignalingProfiler 2.0 addresses the emergent need to derive biologically relevant information from complex multi-omics data by extracting interpretable networks.

Unraveling druggable cancer-driving proteins and targeted drugs using artificial intelligence and multi-omics analyses

The druggable proteome refers to proteins that can bind to small molecules with appropriate chemical affinity, inducing a favorable clinical response. Predicting druggable proteins through screening and in silico modeling is imperative for drug design. To contribute to this field, we developed an accurate predictive classifier for druggable cancer-driving proteins using amino acid composition descriptors of protein sequences and 13 machine learning linear and non-linear classifiers. The optimal classifier was achieved with the support vector machine method, utilizing 200 tri-amino acid composition descriptors. The high performance of the model is evident from an area under the receiver operating characteristics (AUROC) of 0.975 ± 0.003 and an accuracy of 0.929 ± 0.006 (threefold cross-validation). The machine learning prediction model was enhanced with multi-omics approaches, including the target-disease evidence score, the shortest pathways to cancer hallmarks, structure-based ligandability assessment, unfavorable prognostic protein analysis, and the oncogenic variome. Additionally, we performed a drug repurposing analysis to identify drugs with the highest affinity capable of targeting the best predicted proteins. As a result, we identified 79 key druggable cancer-driving proteins with the highest ligandability, and 23 of them demonstrated unfavorable prognostic significance across 16 TCGA PanCancer types: CDKN2A, BCL10, ACVR1, CASP8, JAG1, TSC1, NBN, PREX2, PPP2R1A, DNM2, VAV1, ASXL1, TPR, HRAS, BUB1B, ATG7, MARK3, SETD2, CCNE1, MUTYH, CDKN2C, RB1, and SMARCA4. Moreover, we prioritized 11 clinically relevant drugs targeting these proteins. This strategy effectively predicts and prioritizes biomarkers, therapeutic targets, and drugs for in-depth studies in clinical trials. Scripts are available at https://github.com/muntisa/machine-learning-for-druggable-proteins .

Functional implications of glycans and their curation: insights from the workshop held at the 16th Annual International Biocuration Conference in Padua, Italy

Dynamic changes in protein glycosylation impact human health and disease progression. However, current resources that capture disease and phenotype information focus primarily on the macromolecules within the central dogma of molecular biology (DNA, RNA, proteins). To gain a better understanding of organisms, there is a need to capture the functional impact of glycans and glycosylation on biological processes. A workshop titled "Functional impact of glycans and their curation" was held in conjunction with the 16th Annual International Biocuration Conference to discuss ongoing worldwide activities related to glycan function curation. This workshop brought together subject matter experts, tool developers, and biocurators from over 20 projects and bioinformatics resources. Participants discussed four key topics for each of their resources: (i) how they curate glycan function-related data from publications and other sources, (ii) what type of data they would like to acquire, (iii) what data they currently have, and (iv) what standards they use. Their answers contributed input that provided a comprehensive overview of state-of-the-art glycan function curation and annotations. This report summarizes the outcome of discussions, including potential solutions and areas where curators, data wranglers, and text mining experts can collaborate to address current gaps in glycan and glycosylation annotations, leveraging each other's work to improve their respective resources and encourage impactful data sharing among resources. Database URL: https://wiki.glygen.org/Glycan_Function_Workshop_2023.

Overlap in oncogenic and pro-inflammatory pathways associated with areca nut and nicotine exposure

Background

Betel nut/areca nut/Areca catechu is one of the most commonly used psychoactive substance, and is also a major preventable cause of cancer. Unlike other psychoactive substances, such as nicotine, the mechanisms underlying addiction to areca nuts and related oncogenesis remain elusive. Recent reports suggest a possible overlap in the mechanisms of action of nicotine and areca nuts in the human body. Thus, this study aimed to investigate the interactome of human proteins associated with areca nut exposure and the intricate similarities and differences in the effects of the two psychoactive substances on humans.

Methods

A list of proteins associated with areca nut use was obtained from the available literature using terms from Medical Subject Headings (MeSH). Protein-protein interaction (PPI) networks and functional enrichment were analyzed. The results obtained for both psychoactive substances were compared.

Results

Given the limited number of common proteins (36/226, 16%) in the two sets, a substantial overlap (612/1176 nodes, 52%) was observed in the PPI networks, as well as in Gene Ontology. Areca nuts mainly affect signaling pathways through three hub proteins (alpha serine/threonine-protein kinase, tumor protein 53, and interleukin-6), which are common to both psychoactive substances, as well as two unique hub proteins (epidermal growth factor receptor and master regulator of cell cycle entry and proliferative metabolism). Areca nut-related proteins are associated with unique pathways, such as extracellular matrix organization, lipid storage, and metabolism, which are not found in nicotine-associated proteins.

Conclusions

Areca nuts affect regulatory mechanisms, leading to systemic toxicity and oncogenesis. Areca nuts also affect unique pathways that can be studied as potential markers of exposure, as well as targets for anticancer therapeutic agents.

Bioinformatics analysis of the potential biomarkers of Multiple Sclerosis and Guillain-Barré syndrome

Recent research emphasizes the intricate interplay of genetics and epigenetics in neurological disorders, notably Multiple Sclerosis (MS) and Guillain-Barre Syndrome (GBS), both of which exhibit cardiovascular dysregulation, with GBS often featuring serious bradyarrhythmias requiring prompt recognition and treatment. While cardiovascular autonomic dysfunction in MS is typically less severe, orthostatic intolerance affects around half of MS patients. Their distinction lies in their autoimmune responses, MS is an autoimmune disease affecting the central nervous system, causes demyelination and axon damage, leading to cognitive, ocular, and musculoskeletal dysfunction. In contrast, GBS primarily affects the peripheral nervous system, resulting in paralysis and respiratory complications. Despite their differences, both diseases share environmental risk factors such as viral infections and Vitamin D deficiency. This study aims to explore shared gene expression pathways, functional annotations, and molecular pathways between MS and GBS to enhance diagnostics, pathogenesis understanding, and treatment strategies through molecular analysis techniques. Through the gene expression analysis, five significant genes were found UTS2, TNFSF10, GBP1, VCAN, FOS. Results shows that Common DEGs are linked to apoptosis, bacterial infections, and atherosclerosis. Molecular docking analysis suggests Aflatoxin B1 as a potential therapeutic compound due to its high binding affinity with common differentially expressed proteins.

Comprehensive molecular interaction map of TGFβ induced epithelial to mesenchymal transition in breast cancer

Breast cancer is one of the prevailing cancers globally, with a high mortality rate. Metastatic breast cancer (MBC) is an advanced stage of cancer, characterised by a highly nonlinear, heterogeneous process involving numerous singling pathways and regulatory interactions. Epithelial-mesenchymal transition (EMT) emerges as a key mechanism exploited by cancer cells. Transforming Growth Factor-β (TGFβ)-dependent signalling is attributed to promote EMT in advanced stages of breast cancer. A comprehensive regulatory map of TGFβ induced EMT was developed through an extensive literature survey. The network assembled comprises of 312 distinct species (proteins, genes, RNAs, complexes), and 426 reactions (state transitions, nuclear translocations, complex associations, and dissociations). The map was developed by following Systems Biology Graphical Notation (SBGN) using Cell Designer and made publicly available using MINERVA ( http://35.174.227.105:8080/minerva/?id=Metastatic_Breast_Cancer_1 ). While the complete molecular mechanism of MBC is still not known, the map captures the elaborate signalling interplay of TGFβ induced EMT-promoting MBC. Subsequently, the disease map assembled was translated into a Boolean model utilising CaSQ and analysed using Cell Collective. Simulations of these have captured the known experimental outcomes of TGFβ induced EMT in MBC. Hub regulators of the assembled map were identified, and their transcriptome-based analysis confirmed their role in cancer metastasis. Elaborate analysis of this map may help in gaining additional insights into the development and progression of metastatic breast cancer.

The landscape of cancer-rewired GPCR signaling axes

We explored the dysregulation of G-protein-coupled receptor (GPCR) ligand systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes. Multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes and are associated to specific transcriptional programs and to patient survival patterns. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including, e.g., muscarinic, adenosine, 5-hydroxytryptamine, and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens, which we further validated. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Inference of drug off-target effects on cellular signaling using interactome-based deep learning

Many diseases emerge from dysregulated cellular signaling, and drugs are often designed to target specific signaling proteins. Off-target effects are, however, common and may ultimately result in failed clinical trials. Here we develop a computer model of the cell's transcriptional response to drugs for improved understanding of their mechanisms of action. The model is based on ensembles of artificial neural networks and simultaneously infers drug-target interactions and their downstream effects on intracellular signaling. With this, it predicts transcription factors' activities, while recovering known drug-target interactions and inferring many new ones, which we validate with an independent dataset. As a case study, we analyze the effects of the drug Lestaurtinib on downstream signaling. Alongside its intended target, FLT3, the model predicts an inhibition of CDK2 that enhances the downregulation of the cell cycle-critical transcription factor FOXM1. Our approach can therefore enhance our understanding of drug signaling for therapeutic design.

Unveiling the signaling network of FLT3-ITD AML improves drug sensitivity prediction

Currently, the identification of patient-specific therapies in cancer is mainly informed by personalized genomic analysis. In the setting of acute myeloid leukemia (AML), patient-drug treatment matching fails in a subset of patients harboring atypical internal tandem duplications (ITDs) in the tyrosine kinase domain of the FLT3 gene. To address this unmet medical need, here we develop a systems-based strategy that integrates multiparametric analysis of crucial signaling pathways, and patient-specific genomic and transcriptomic data with a prior knowledge signaling network using a Boolean-based formalism. By this approach, we derive personalized predictive models describing the signaling landscape of AML FLT3-ITD positive cell lines and patients. These models enable us to derive mechanistic insight into drug resistance mechanisms and suggest novel opportunities for combinatorial treatments. Interestingly, our analysis reveals that the JNK kinase pathway plays a crucial role in the tyrosine kinase inhibitor response of FLT3-ITD cells through cell cycle regulation. Finally, our work shows that patient-specific logic models have the potential to inform precision medicine approaches.

Exploration of functional relations among differentially co-expressed genes identifies regulators in glioblastoma

The conventional computational approaches to investigating a disease confront inherent constraints as they often need to improve in delving beyond protein functional associations and grasping their deeper contextual significance within the disease framework. Such context-specificity can be explored using clinical data by evaluating the change in interaction between the biological entities in different conditions by investigating the differential co-expression relationships. We believe that the integration and analysis of differential co-expression and the functional relationships, primarily focusing on the source nodes, will open novel insights about disease progression as the source proteins could trigger signaling cascades, mostly because they are transcription factors, cell surface receptors, or enzymes that respond instantly to a particular stimulus. A thorough contextual investigation of these nodes could lead to a helpful beginning point for identifying potential causal linkages and guiding subsequent scientific investigations to uncover mechanisms underlying observed associations. Our methodology includes functional protein-protein Interaction (PPI) data and co-expression information and filters functional linkages through a series of critical steps, culminating in the identification of a robust set of regulators. Our analysis identified eleven key regulators-AKT1, BRCA1, CAMK2G, CUL1, FGFR3, KIF3A, NUP210, PRKACB, RAB8A, RPS6KA2 and TGFB3-in glioblastoma. These regulators play a pivotal role in disease classification, cell growth control, and patient survivability and exhibit associations with immune infiltrations and disease hallmarks. This underscores the importance of assessing correlation towards causality in unraveling complex biological insights.

Transomics2cytoscape: an automated software for interpretable 2.5-dimensional visualization of trans-omic networks

Biochemical network visualization is one of the essential technologies for mechanistic interpretation of omics data. In particular, recent advances in multi-omics measurement and analysis require the development of visualization methods that encompass multiple omics data. Visualization in 2.5 dimension (2.5D visualization), which is an isometric view of stacked X-Y planes, is a convenient way to interpret multi-omics/trans-omics data in the context of the conventional layouts of biochemical networks drawn on each of the stacked omics layers. However, 2.5D visualization of trans-omics networks is a state-of-the-art method that primarily relies on time-consuming human efforts involving manual drawing. Here, we present an R Bioconductor package 'transomics2cytoscape' for automated visualization of 2.5D trans-omics networks. We confirmed that transomics2cytoscape could be used for rapid visualization of trans-omics networks presented in published papers within a few minutes. Transomics2cytoscape allows for frequent update/redrawing of trans-omics networks in line with the progress in multi-omics/trans-omics data analysis, thereby enabling network-based interpretation of multi-omics data at each research step. The transomics2cytoscape source code is available at https://github.com/ecell/transomics2cytoscape .

Camptothecin structure simplification elaborated new imidazo[2,1-b]quinazoline derivative as a human topoisomerase I inhibitor with efficacy against bone cancer cells and colon adenocarcinoma

Camptothecin is a pentacyclic natural alkaloid that inhibits the hTop1 enzyme involved in DNA transcription and cancer cell growth. Camptothecin structure pitfalls prompted us to design new congeners using a structure simplification strategy to reduce the ring extension number from pentacyclic to tetracyclic while maintaining potential stacking of the new compounds with the DNA base pairs at the Top1-mediated cleavage complex and aqueous solubility, as well as minimizing compound-liver toxicity. The principal axis of this study was the verification of hTop1 inhibiting activity as a possible mechanism of action and the elaboration of new simplified inhibitors with improved pharmacodynamic and pharmacokinetic profiling using three structure panels (A-C) of (isoquinolinoimidazoquinazoline), (imidazoquinazoline), and (imidazoisoquinoline), respectively. DNA relaxation assay identified five compounds as hTop1 inhibitors belonging to the imidazoisoquinolines 3a,b, the imidazoquinazolines 12, and the isoquinolinoimidazoquinazolines 7a,b. In an MTT cytotoxicity assay against different cancer cell lines, compound 12 was the most potent against HOS bone cancer cells (IC50 = 1.47 μM). At the same time, the other inhibitors had no detectable activity against any cancer cell type. Compound (12) demonstrated great penetrating power in the HOS cancer cells' 3D-multicellular tumor spheroid model. Bioinformatics research of the hTop1 gene revealed that the TP53 cell proliferative gene is in the network of hTop1. The finding is confirmed empirically using the gene expression assay that proved the increase in p53 expression. The impact of structure simplification on compound 12 profile, characterized by the absence of acute oral liver toxicity when compared to Doxorubicin as a standard inhibitor, the lethal dose measured on Swiss Albino female mice and reported at LD50 = 250 mg/kg, and therapeutic significance in reducing colon adenocarcinoma tumor volume by 75.36 % after five weeks of treatment with compound 12. The molecular docking solutions of the active CPT-based derivative 12 and the inactive congener 14 into the active site of hTop1 and the activity cliffing of such MMP directed us to recommend the addition of HBD and HBA variables to compound 12 imidazoquinazoline core scaffold to enhance the potency via hydrogen bond formation with the major groove amino acids (Asp533, Lys532) as well as maintaining the hydrogen bond with the minor groove amino acid Arg364.

CardiOmics signatures reveal therapeutically actionable targets and drugs for cardiovascular diseases

Cardiovascular diseases are the leading cause of death worldwide, with heart failure being a complex condition that affects millions of individuals. Single-nucleus RNA sequencing has recently emerged as a powerful tool for unraveling the molecular mechanisms behind cardiovascular diseases. This cutting-edge technology enables the identification of molecular signatures, intracellular networks, and spatial relationships among cardiac cells, including cardiomyocytes, mast cells, lymphocytes, macrophages, lymphatic endothelial cells, endocardial cells, endothelial cells, epicardial cells, adipocytes, fibroblasts, neuronal cells, pericytes, and vascular smooth muscle cells. Despite these advancements, the discovery of essential therapeutic targets and drugs for precision cardiology remains a challenge. To bridge this gap, we conducted comprehensive in silico analyses of single-nucleus RNA sequencing data, functional enrichment, protein interactome network, and identification of the shortest pathways to physiological phenotypes. This integrated multi-omics analysis generated CardiOmics signatures, which allowed us to pinpoint three therapeutically actionable targets (ADRA1A1, PPARG, and ROCK2) and 15 effective drugs, including adrenergic receptor agonists, adrenergic receptor antagonists, norepinephrine precursors, PPAR receptor agonists, and Rho-associated kinase inhibitors, involved in late-stage cardiovascular disease clinical trials.

PathDIP 5: improving coverage and making enrichment analysis more biologically meaningful

Pathway Data Integration Portal (PathDIP) is an integrated pathway database that was developed to increase functional gene annotation coverage and reduce bias in pathway enrichment analysis. PathDIP 5 provides multiple improvements to enable more interpretable analysis: users can perform enrichment analysis using all sources, separate sources or by combining specific pathway subsets; they can select the types of sources to use or the types of pathways for the analysis, reducing the number of resulting generic pathways or pathways not related to users' research question; users can use API. All pathways have been mapped to seven representative types. The results of pathway enrichment can be summarized through knowledge-based pathway consolidation. All curated pathways were mapped to 53 pathway ontology-based categories. In addition to genes, pathDIP 5 now includes metabolites. We updated existing databases, included two new sources, PathBank and MetabolicAtlas, and removed outdated databases. We enable users to analyse their results using Drugst.One, where a drug-gene network is created using only the user's genes in a specific pathway. Interpreting the results of any analysis is now improved by multiple charts on all the results pages. PathDIP 5 is freely available at https://ophid.utoronto.ca/pathDIP.

omicSynth: An open multi-omic community resource for identifying druggable targets across neurodegenerative diseases

Treatments for neurodegenerative disorders remain rare, but recent FDA approvals, such as lecanemab and aducanumab for Alzheimer disease (MIM: 607822), highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use summary-data-based Mendelian randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer disease, 3 amyotrophic lateral sclerosis (MIM: 105400), 5 Lewy body dementia (MIM: 127750), 46 Parkinson disease (MIM: 605909), and 9 progressive supranuclear palsy (MIM: 601104) target genes passing multiple test corrections (pSMR_multi < 2.95 × 10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics, classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these, 69.8% are expressed in the disease-relevant cell type from single-nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development, and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as riluzole in Alzheimer disease. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community.

The landscape of SETBP1 gene expression and transcription factor activity across human tissues

The SET binding protein 1 (SETBP1) gene encodes a transcription factor (TF) involved in various cellular processes. Variants in SETBP1 can result in three different diseases determined by the introduction (germline vs. somatic) and location of the variant. Germline variants cause the ultra-rare pediatric Schinzel Giedion Syndrome (SGS) and SETBP1 haploinsufficiency disorder (SETBP1-HD), characterized by severe multisystemic abnormalities with neurodegeneration or a less severe brain phenotype accompanied by hypotonia and strabismus, respectively. Somatic variants in SETBP1 are associated with hematological malignancies and cancer development in other tissues in adults. To better understand the tissue-specific mechanisms involving SETBP1, we analyzed publicly available RNA-sequencing (RNA-seq) data from the Genotype-Tissue Expression (GTEx) project. We found SETBP1 and its known target genes were widely expressed across 31 adult human tissues. K-means clustering identified three distinct expression patterns of SETBP1 targets across tissues. Functional enrichment analysis (FEA) of each cluster revealed gene sets related to transcriptional regulation, DNA binding, and mitochondrial function. TF activity analysis of SETBP1 and its target TFs revealed tissue-specific TF activity, underscoring the role of tissue context-driven regulation and suggesting its impact in SETBP1-associated disease. In addition to uncovering tissue-specific molecular signatures of SETBP1 expression and TF activity, we provide a Shiny web application to facilitate exploring TF activity across human tissues for 758 TFs. This study provides insight into the landscape of SETBP1 expression and TF activity across 31 non-diseased human tissues and reveals tissue-specific expression and activity of SETBP1 and its targets. In conjunction with the web application we constructed, our framework enables researchers to generate hypotheses related to the role tissue backgrounds play with respect to gene expression and TF activity in different disease contexts.

Curation of causal interactions mediated by genes associated with autism accelerates the understanding of gene-phenotype relationships underlying neurodevelopmental disorders

Autism spectrum disorder (ASD) comprises a large group of neurodevelopmental conditions featuring, over a wide range of severity and combinations, a core set of manifestations (restricted sociality, stereotyped behavior and language impairment) alongside various comorbidities. Common and rare variants in several hundreds of genes and regulatory regions have been implicated in the molecular pathogenesis of ASD along a range of causation evidence strength. Despite significant progress in elucidating the impact of few paradigmatic individual loci, such sheer complexity in the genetic architecture underlying ASD as a whole has hampered the identification of convergent actionable hubs hypothesized to relay between the vastness of risk alleles and the core phenotypes. In turn this has limited the development of strategies that can revert or ameliorate this condition, calling for a systems-level approach to probe the cross-talk of cooperating genes in terms of causal interaction networks in order to make convergences experimentally tractable and reveal their clinical actionability. As a first step in this direction, we have captured from the scientific literature information on the causal links between the genes whose variants have been associated with ASD and the whole human proteome. This information has been annotated in a computer readable format in the SIGNOR database and is made freely available in the resource website. To link this information to cell functions and phenotypes, we have developed graph algorithms that estimate the functional distance of any protein in the SIGNOR causal interactome to phenotypes and pathways. The main novelty of our approach resides in the possibility to explore the mechanistic links connecting the suggested gene-phenotype relations.

Knowledge Graphs and Their Applications in Drug Discovery

Knowledge graphs represent information in the form of entities and relationships between those entities. Such a representation has multiple potential applications in drug discovery, including democratizing access to biomedical data, contextualizing or visualizing that data, and generating novel insights through the application of machine learning approaches. Knowledge graphs put data into context and therefore offer the opportunity to generate explainable predictions, which is a key topic in contemporary artificial intelligence. In this chapter, we outline some of the factors that need to be considered when constructing biomedical knowledge graphs, examine recent advances in mining such systems to gain insights for drug discovery, and identify potential future areas for further development.

Specifying cellular context of transcription factor regulons for exploring context-specific gene regulation programs

Understanding the role of transcription and transcription factors in cellular identity and disease, such as cancer and autoimmunity, is essential. However, comprehensive data resources for cell line-specific transcription factor-to-target gene annotations are currently limited. To address this, we developed a straightforward method to define regulons that capture the cell-specific aspects of TF binding and transcript expression levels. By integrating cellular transcriptome and transcription factor binding data, we generated regulons for four common cell lines comprising both proximal and distal cell line-specific regulatory events. Through systematic benchmarking involving transcription factor knockout experiments, we demonstrated performance on par with state-of-the-art methods, with our method being easily applicable to other cell types of interest. We present case studies using three cancer single-cell datasets to showcase the utility of these cell-type-specific regulons in exploring transcriptional dysregulation. In summary, this study provides a valuable tool and a resource for systematically exploring cell line-specific transcriptional regulations, emphasizing the utility of network analysis in deciphering disease mechanisms.

Precise Therapeutic Targeting of Distinct NRXN1+/- Mutations

As genetic studies continue to identify risk loci that are significantly associated with risk for neuropsychiatric disease, a critical unanswered question is the extent to which diverse mutations--sometimes impacting the same gene-- will require tailored therapeutic strategies. Here we consider this in the context of rare neuropsychiatric disorder-associated copy number variants (2p16.3) resulting in heterozygous deletions in NRXN1, a pre-synaptic cell adhesion protein that serves as a critical synaptic organizer in the brain. Complex patterns of NRXN1 alternative splicing are fundamental to establishing diverse neurocircuitry, vary between the cell types of the brain, and are differentially impacted by unique (non-recurrent) deletions. We contrast the cell-type-specific impact of patient-specific mutations in NRXN1 using human induced pluripotent stem cells, finding that perturbations in NRXN1 splicing result in divergent cell-type-specific synaptic outcomes. Via distinct loss-of-function (LOF) and gain-of-function (GOF) mechanisms, NRXN1+/- deletions cause decreased synaptic activity in glutamatergic neurons, yet increased synaptic activity in GABAergic neurons. Stratification of patients by LOF and GOF mechanisms will facilitate individualized restoration of NRXN1 isoform repertoires; towards this, antisense oligonucleotides knockdown mutant isoform expression and alters synaptic transcriptional signatures, while treatment with β-estradiol rescues synaptic function in glutamatergic neurons. Given the increasing number of mutations predicted to engender both LOF and GOF mechanisms in brain disease, our findings add nuance to future considerations of precision medicine.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

Prechordal structures act cooperatively in early trabeculae development of gnathostome skull

The vertebrate skull is formed by mesoderm and neural crest (NC) cells. The mesoderm contributes to the skull chordal domain, with the notochord playing an essential role in this process. The NC contributes to the skull prechordal domain, prompting investigation into the embryonic structures involved in prechordal neurocranium cartilage formation. The trabeculae cartilage, a structure of the prechordal neurocranium, arises at the convergence of prechordal plate (PCP), ventral midline (VM) cells of the diencephalon, and dorsal oral ectoderm. This study examines the molecular participation of these embryonic structures in gnathostome trabeculae development. PCP-secreted SHH induces its expression in VM cells of the diencephalon, initiating a positive feedback loop involving SIX3 and GLI1. SHH secreted by the VM cells of the diencephalon acts on the dorsal oral ectoderm, stimulating condensation of NC cells to form trabeculae. SHH from the prechordal region affects the expression of SOX9 in NC cells. BMP7 and SHH secreted by PCP induce NKX2.1 expression in VM cells of the diencephalon, but this does not impact trabeculae formation. Molecular cooperation between PCP, VM cells of the diencephalon, and dorsal oral ectoderm is crucial for craniofacial development by NC cells in the prechordal domain.

Monocytes subsets altered distribution and dysregulated plasma hsa-miR-21-5p and hsa-miR-155-5p in HCV-linked liver cirrhosis progression to hepatocellular carcinoma

PURPOSE

The authors aim to investigate the altered monocytes subsets distribution in liver cirrhosis (LC) and subsequent hepatocellular carcinoma (HCC) in association with the expression level of plasma Homo sapiens (has)-miR-21-5p and hsa-miR-155-5p. A step toward non-protein coding (nc) RNA precision medicine based on the immune perturbation manifested as altered monocytes distribution, on top of LC and HCC.

METHODS

Seventy-nine patients diagnosed with chronic hepatitis C virus (CHCV) infection with LC were enrolled in the current study. Patients were sub-classified into LC group without HCC (n = 40), LC with HCC (n = 39), and 15 apparently healthy controls. Monocyte subsets frequencies were assessed by flow cytometry. Real-time quantitative PCR was used to measure plasma hsa-miR-21-5p and hsa-miR-155-5p expression.

RESULTS

Hsa-miR-21-5p correlated with intermediate monocytes (r = 0.30, p = 0.007), while hsa-miR-155-5p negatively correlated with non-classical monocytes (r = - 0.316, p = 0.005). ROC curve analysis revealed that combining intermediate monocytes frequency and hsa-miR-21 yielded sensitivity = 79.5%, specificity = 75%, and AUC = 0.84. In comparison, AFP yielded a lower sensitivity = 69% and 100% specificity with AUC = 0.85. Logistic regression analysis proved that up-regulation of intermediate monocytes frequency and hsa-miR-21-5p were independent risk factors for LC progression to HCC, after adjustment for co-founders.

CONCLUSION

Monocyte subsets differentiation in HCC was linked to hsa-miR-21-5p and hsa-miR-155-5p. Combined up-regulation of intermediate monocytes frequency and hsa-miR-21-5p expression could be considered a sensitive indicator of LC progression to HCC. Circulating intermediate monocytes and hsa-miR-21-5p were independent risk factors for HCC evolution, clinically and in silico proved.

A framework for considering prior information in network-based approaches to omics data analysis

For decades, molecular biologists have been uncovering the mechanics of biological systems. Efforts to bring their findings together have led to the development of multiple databases and information systems that capture and present pathway information in a computable network format. Concurrently, the advent of modern omics technologies has empowered researchers to systematically profile cellular processes across different modalities. Numerous algorithms, methodologies, and tools have been developed to use prior knowledge networks (PKNs) in the analysis of omics datasets. Interestingly, it has been repeatedly demonstrated that the source of prior knowledge can greatly impact the results of a given analysis. For these methods to be successful it is paramount that their selection of PKNs is amenable to the data type and the computational task they aim to accomplish. Here we present a five-level framework that broadly describes network models in terms of their scope, level of detail, and ability to inform causal predictions. To contextualize this framework, we review a handful of network-based omics analysis methods at each level, while also describing the computational tasks they aim to accomplish.

The Cyclin-dependent kinase 1: more than a cell cycle regulator

The Cyclin-dependent kinase 1, as a serine/threonine protein kinase, is more than a cell cycle regulator as it was originally identified. During the last decade, it has been shown to carry out versatile functions during the last decade. From cell cycle control to gene expression regulation and apoptosis, CDK1 is intimately involved in many cellular events that are vital for cell survival. Here, we provide a comprehensive catalogue of the CDK1 upstream regulators and substrates, describing how this kinase is implicated in the control of key 'cell cycle-unrelated' biological processes. Finally, we describe how deregulation of CDK1 expression and activation has been closely associated with cancer progression and drug resistance.

The landscape of SETBP1 gene expression and transcription factor activity across human tissues

Background

The SET binding protein 1 (SETBP1) gene encodes a transcription factor (TF) involved in various cellular processes. Distinct SETBP1 variants have been linked to three different diseases. Germline variants cause the ultra-rare pediatric Schinzel Giedion Syndrome (SGS) and SETBP1 haploinsufficiency disorder (SETBP1-HD), characterized by severe multisystemic abnormalities with neurodegeneration or a less severe brain phenotype accompanied by hypotonia and strabismus, respectively. Somatic variants in SETBP1 are associated with hematological malignancies and cancer development in other tissues in adults.

Results

To better understand the tissue-specific mechanisms involving SETBP1, we analyzed publicly available RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project. We found SETBP1, and its known target genes were widely expressed across 31 adult human tissues. K-means clustering identified three distinct expression patterns of SETBP1 targets across tissues. Functional enrichment analysis (FEA) of each cluster revealed gene sets related to transcription regulation, DNA binding, and mitochondrial function. TF activity analysis of SETBP1 and its target TFs revealed tissue-specific TF activity, underscoring the role of tissue context-driven regulation and suggesting its impact in SETBP1-associated disease. In addition to uncovering tissue-specific molecular signatures of SETBP1 expression and TF activity, we provide a Shiny web application to facilitate exploring TF activity across human tissues for 758 TFs.

Conclusions

This study provides insight into the landscape of SETBP1 expression and TF activity across 31 non-diseased human tissues and reveals tissue-specific expression and activity of SETBP1 and its targets. In conjunction with the web application we constructed, our framework enables researchers to generate hypotheses related to the role tissue backgrounds play with respect to gene expression and TF activity in different disease contexts.

The landscape of cancer rewired GPCR signaling axes

We explored the dysregulation of GPCR ligand signaling systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes, which revealed that multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes. We showed that biosynthetic pathway enrichment from enzyme expression recapitulated pathway activity signatures from metabolomics datasets, providing valuable surrogate information for GPCRs responding to organic ligands. We found that several GPCRs signaling components were significantly associated with patient survival in a cancer type-specific fashion. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many pairs involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including e.g., muscarinic, adenosine, 5-hydroxytryptamine and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

A Network of 17 Microtubule-Related Genes Highlights Functional Deregulations in Breast Cancer

A wide panel of microtubule-associated proteins and kinases is involved in coordinated regulation of the microtubule cytoskeleton and may thus represent valuable molecular markers contributing to major cellular pathways deregulated in cancer. We previously identified a panel of 17 microtubule-related (MT-Rel) genes that are differentially expressed in breast tumors showing resistance to taxane-based chemotherapy. In the present study, we evaluated the expression, prognostic value and functional impact of these genes in breast cancer. We show that 14 MT-Rel genes (KIF4A, ASPM, KIF20A, KIF14, TPX2, KIF18B, KIFC1, AURKB, KIF2C, GTSE1, KIF15, KIF11, RACGAP1, STMN1) are up-regulated in breast tumors compared with adjacent normal tissue. Six of them (KIF4A, ASPM, KIF20A, KIF14, TPX2, KIF18B) are overexpressed by more than 10-fold in tumor samples and four of them (KIF11, AURKB, TPX2 and KIFC1) are essential for cell survival. Overexpression of all 14 genes, and underexpression of 3 other MT-Rel genes (MAST4, MAPT and MTUS1) are associated with poor breast cancer patient survival. A Systems Biology approach highlighted three major functional networks connecting the 17 MT-Rel genes and their partners, which are centered on spindle assembly, chromosome segregation and cytokinesis. Our studies identified mitotic Aurora kinases and their substrates as major targets for therapeutic approaches against breast cancer.

Inferring single-cell gene regulatory network by non-redundant mutual information

Gene regulatory network plays a crucial role in controlling the biological processes of living creatures. Deciphering the complex gene regulatory networks from experimental data remains a major challenge in system biology. Recent advances in single-cell RNA sequencing technology bring massive high-resolution data, enabling computational inference of cell-specific gene regulatory networks (GRNs). Many relevant algorithms have been developed to achieve this goal in the past years. However, GRN inference is still less ideal due to the extra noises involved in pseudo-time information and large amounts of dropouts in datasets. Here, we present a novel GRN inference method named Normi, which is based on non-redundant mutual information. Normi manipulates these problems by employing a sliding size-fixed window approach on the entire trajectory and conducts average smoothing strategy on the gene expression of the cells in each window to obtain representative cells. To further alleviate the impact of dropouts, we utilize the mixed KSG estimator to quantify the high-order time-delayed mutual information among genes, then filter out the redundant edges by adopting Max-Relevance and Min Redundancy algorithm. Moreover, we determined the optimal time delay for each gene pair by distance correlation. Normi outperforms other state-of-the-art GRN inference methods on both simulated data and single-cell RNA sequencing (scRNA-seq) datasets, demonstrating its superiority in robustness. The performance of Normi in real scRNA-seq data further reveals its ability to identify the key regulators and crucial biological processes.

omicSynth: an Open Multi-omic Community Resource for Identifying Druggable Targets across Neurodegenerative Diseases

Treatments for neurodegenerative disorders remain rare, although recent FDA approvals, such as Lecanemab and Aducanumab for Alzheimer's Disease, highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use Summary-data-based Mendelian Randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer's disease, 3 amyotrophic lateral sclerosis, 5 Lewy body dementia, 46 Parkinson's disease, and 9 Progressive supranuclear palsy target genes passing multiple test corrections (pSMR_multi < 2.95×10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics - classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these 69.8% are expressed in the disease relevant cell type from single nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as Riluzole in AD. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community [https://nih-card-ndd-smr-home-syboky.streamlit.app/].

Identifying miRNA biomarkers for breast cancer and ovarian cancer: a text mining perspective

BACKGROUND

microRNA (miRNAs) are small, non-coding RNAs that mediate post-transcriptional gene silencing. Numerous studies have demonstrated the critical role of miRNAs in the development of breast cancer and ovarian cancer. To reduce potential bias from individual studies, a more comprehensive approach of exploring miRNAs in cancer research is essential. This study aims to explore the role of miRNAs in the development of breast cancer and ovarian cancer.

METHODS

Abstracts of the publications were tokenized and the biomedical terms (miRNA, gene, disease, species) were identified and extracted for vectorization. Predictive analyses were conducted with four machine learning models: K-Nearest Neighbors (KNN), Support Vector Machines (SVM), Random Forest (RF), and Naïve Bayes. Both holdout validation and cross-validation were utilized. Feature importance will be identified for miRNA-cancer networks construction.

RESULTS

We found that miR-182 is highly specific to female cancers. miR-182 targets different genes in regulating breast cancer and ovarian cancer. Naïve Bayes provided a promising prediction model for breast cancer and ovarian cancer with miRNAs and genes combination, with an accuracy score greater than 60%. Feature importance identified miR-155 and miR-199 are critical for breast cancer and ovarian cancer prediction, with miR-155 being highly related to breast cancer, whereas miR-199 being more associated with ovarian cancer.

CONCLUSION

Our approach effectively identified potential miRNA biomarkers associated with breast cancer and ovarian cancer, providing a solid foundation for generating novel research hypotheses and guiding future experimental studies.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

New trends in synthetic drugs and natural products targeting 20S proteasomes in cancers

Cancer is a global health challenge that remains to be a field of extensive research aiming to find new anticancer therapeutics. The 20S proteasome complex is one of the targets of anticancerdrugs, as it is correlated with several cancer types. Herein, we aim to discuss the 20S proteasome subunits and investigatethe currently studied proteasome inhibitors targeting the catalytically active proteasome subunits. In this review, we summarize the proteindegradation mechanism of the 20S proteasome complex and compareit with the 26S proteasome complex. Afterwards, the localization of the 20S proteasome is summarized as well as its use as a diagnosticandprognostic marker. The FDA-approved proteasome inhibitors (PIs) under clinical trials are summarized and their current limited use in solid tumors is also reviewed in addition to the expression of theβ5 subunit in differentcell lines. The review discusses in-silico analysis of the active subunit of the 20S proteasome complex. For development of new proteasome inhibitor drugs, the natural products inhibiting the 20S proteasome are summarized, as well as novel methodologies and challenges for the natural product discovery and current information about the biosynthetic gene clusters encoding them. We herein briefly summarize some resistancemechanismsto the proteasomeinhibitors. Additionally, we focus on the three main classes of proteasome inhibitors: 1] boronic acid, 2] beta-lactone and 3] epoxide inhibitor classes, as well as other PI classes, and their IC₅₀ values and their structure-activity relationship (SAR). Lastly,we summarize several future prospects of developing new proteasome inhibitors towards the treatment of tumors, especially solid tumors.

Recent advances in kinase signaling network profiling by mass spectrometry

Mass spectrometry-based phosphoproteomics is currently the leading methodology for the study of global kinase signaling. The scientific community is continuously releasing technological improvements for sensitive and fast identification of phosphopeptides, and their accurate quantification. To interpret large-scale phosphoproteomics data, numerous bioinformatic resources are available that help understanding kinase network functional role in biological systems upon perturbation. Some of these resources are databases of phosphorylation sites, protein kinases and phosphatases; others are bioinformatic algorithms to infer kinase activity, predict phosphosite functional relevance and visualize kinase signaling networks. In this review, we present the latest experimental and bioinformatic tools to profile protein kinase signaling networks and provide examples of their application in biomedicine.