Signaling networks are critical for virtually all cell functions. Our current knowledge of cell signaling has been summarized in signaling pathway databases, which, while useful, are highly biased toward well-studied processes, and do not capture context specific network wiring or pathway cross-talk. Mass spectrometry-based phosphoproteomics data can provide a more unbiased view of active cell signaling processes in a given context, however, it suffers from low signal-to-noise ratio and poor reproducibility across experiments. While progress in methods to extract active signaling signatures from such data has been made, there are still limitations with respect to balancing bias and interpretability. Here we present phuEGO, which combines up-to-three-layer network propagation with ego network decomposition to provide small networks comprising active functional signaling modules. PhuEGO boosts the signal-to-noise ratio from global phosphoproteomics datasets, enriches the resulting networks for functional phosphosites and allows the improved comparison and integration across datasets. We applied phuEGO to five phosphoproteomics data sets from cell lines collected upon infection with SARS CoV2. PhuEGO was better able to identify common active functions across datasets and to point to a subnetwork enriched for known COVID-19 targets. Overall, phuEGO provides a flexible tool to the community for the improved functional interpretation of global phosphoproteomics datasets.

Autism spectrum disorder (ASD) is a highly heritable complex disease that affects 1% of the population, yet its underlying molecular mechanisms are largely unknown. Here we study the problem of predicting causal genes for ASD by combining genome-scale data with a network propagation approach. We construct a predictor that integrates multiple omic data sets that assess genomic, transcriptomic, proteomic, and phosphoproteomic associations with ASD. In cross validation our predictor yields mean area under the ROC curve of 0.87 and area under the precision-recall curve of 0.89. We further show that it outperforms previous gene-level predictors of autism association. Finally, we show that we can use the model to predict genes associated with Schizophrenia which is known to share genetic components with ASD.

Kinase fusion genes are the most active fusion gene group in human cancer fusion genes. To help choose the clinically significant kinase so that the cancer patients that have fusion genes can be better diagnosed, we need a metric to infer the assessment of kinases in pan-cancer fusion genes rather than relying on the sample frequency expressed fusion genes. Most of all, multiple studies assessed human kinases as the drug targets using multiple types of genomic and clinical information, but none used the kinase fusion genes in their study. The assessment studies of kinase without kinase fusion gene events can miss the effect of one of the mechanisms that enhance the kinase function in cancer. To fill this gap, in this study, we suggest a novel way of assessing genes using a network propagation approach to infer how likely individual kinases influence the kinase fusion gene network composed of ~5K kinase fusion gene pairs. To select a better seed of propagation, we chose the top genes via dimensionality reduction like a principal component or latent layer information of six features of individual genes in pan-cancer fusion genes. Our approach may provide a novel way to assess of human kinases in cancer.

BACKGROUND: Genome-wide association studies (GWAS) have enabled large-scale analysis of the role of genetic variants in human disease. Despite impressive methodological advances, subsequent clinical interpretation and application remains challenging when GWAS suffer from a lack of statistical power. In recent years, however, the use of information diffusion algorithms with molecular networks has led to fruitful insights on disease genes.
RESULTS: We present an overview of the design choices and pitfalls that prove crucial in the application of network propagation methods to GWAS summary statistics. We highlight general trends from the literature, and present benchmark experiments to expand on these insights selecting as case study three diseases and five molecular networks. We verify that the use of gene-level scores based on GWAS P-values offers advantages over the selection of a set of \'seed\' disease genes not weighted by the associated P-values if the GWAS summary statistics are of sufficient quality. Beyond that, the size and the density of the networks prove to be important factors for consideration. Finally, we explore several ensemble methods and show that combining multiple networks may improve the network propagation approach.

Combination therapy is a promising strategy for cancers, increasing therapeutic options and reducing drug resistance. Yet, systematic identification of efficacious drug combinations is limited by the combinatorial explosion caused by a large number of possible drug pairs and diseases. At present, machine learning techniques have been widely applied to predict drug combinations, but most studies rely on the response of drug combinations to specific cell lines and are not entirely satisfactory in terms of mechanism interpretability and model scalability. Here, we proposed a novel network propagation-based machine learning framework to predict synergistic drug combinations. Based on the topological information of a comprehensive drug-drug association network, we innovatively introduced an affinity score between drug pairs as one of the features to train machine learning models. We applied network-based strategy to evaluate their therapeutic potential to different cancer types. Finally, we identified 17 specific-, 21 general- and 40 broad-spectrum antitumor drug combinations, in which 69% drug combinations were validated by vitro cellular experiments, 83% drug combinations were validated by literature reports and 100% drug combinations were validated by biological function analyses. By quantifying the network relationships between drug targets and cancer-related driver genes in the human protein-protein interactome, we show the existence of four distinct patterns of drug-drug-disease relationships. We also revealed that 32 biological pathways were correlated with the synergistic mechanism of broad-spectrum antitumor drug combinations. Overall, our model offers a powerful scalable screening framework for cancer treatments.

Inflammatory cytokines released by synovium after trauma disturb the gene regulatory network and have been implicated in the pathophysiology of osteoarthritis. A mechanistic understanding of how aging perturbs this process can help identify novel interventions. Here, we introduced network paradigms to simulate cytokine-mediated pathological communication between the synovium and cartilage. Cartilage-specific network analysis of injured young and aged murine knees revealed aberrant matrix remodeling as a transcriptomic response unique to aged knees displaying accelerated cartilage degradation. Next, network-based cytokine inference with pharmacological manipulation uncovered IL6 family member, Oncostatin M (OSM), as a driver of the aberrant matrix remodeling. By implementing a phenotypic drug discovery approach, we identified that the activation of OSM recapitulated an \"inflammatory\" phenotype of knee osteoarthritis and highlighted high-value targets for drug development and repurposing. These findings offer translational opportunities targeting the inflammation-driven osteoarthritis phenotype.

UNASSIGNED: Lead identification is a fundamental step to prioritize candidate compounds for downstream drug discovery process. Machine learning (ML) and deep learning (DL) approaches are widely used to identify lead compounds using both chemical property and experimental information. However, ML or DL methods rarely consider compound similarity information directly since ML and DL models use abstract representation of molecules for model construction. Alternatively, data mining approaches are also used to explore chemical space with drug candidates by screening undesirable compounds. A major challenge for data mining approaches is to develop efficient data mining methods that search large chemical space for desirable lead compounds with low false positive rate.
UNASSIGNED: In this work, we developed a network propagation (NP) based data mining method for lead identification that performs search on an ensemble of chemical similarity networks. We compiled 14 fingerprint-based similarity networks. Given a target protein of interest, we use a deep learning-based drug target interaction model to narrow down compound candidates and then we use network propagation to prioritize drug candidates that are highly correlated with drug activity score such as IC50. In an extensive experiment with BindingDB, we showed that our approach successfully discovered intentionally unlabeled compounds for given targets. To further demonstrate the prediction power of our approach, we identified 24 candidate leads for CLK1. Two out of five synthesizable candidates were experimentally validated in binding assays. In conclusion, our framework can be very useful for lead identification from very large compound databases such as ZINC.

Introduction: Network-based methods are promising approaches in systems toxicology because they can be used to predict the effects of drugs and chemicals on health, to elucidate the mode of action of compounds, and to identify biomarkers of toxicity. Over the years, the network biology community has developed a wide range of methods, and users are faced with the task of choosing the most appropriate method for their own application. Furthermore, the advantages and limitations of each method are difficult to determine without a proper standard and comparative evaluation of their performance. This study aims to evaluate different network-based methods that can be used to gain biological insight into the mechanisms of drug toxicity, using valproic acid (VPA)-induced liver steatosis as a benchmark. Methods: We provide a comprehensive analysis of the results produced by each method and highlight the fact that the experimental design (how the method is applied) is relevant in addition to the method specifications. We also contribute with a systematic methodology to analyse the results of the methods individually and in a comparative manner. Results: Our results show that the evaluated tools differ in their performance against the benchmark and in their ability to provide novel insights into the mechanism of adverse effects of the drug. We also suggest that aggregation of the results provided by different methods provides a more confident set of candidate genes and processes to further the knowledge of the drug\'s mechanism of action. Discussion: By providing a detailed and systematic analysis of the results of different network-based tools, we aim to assist users in making informed decisions about the most appropriate method for systems toxicology applications.

Alzheimer\'s disease (AD) poses a profound human, social, and economic burden. Previous studies suggest that extra virgin olive oil (EVOO) may be helpful in preventing cognitive decline. Here, we present a network machine learning method for identifying bioactive phytochemicals in EVOO with the highest potential to impact the protein network linked to the development and progression of the AD. A balanced classification accuracy of 70.3 ± 2.6% was achieved in fivefold cross-validation settings for predicting late-stage experimental drugs targeting AD from other clinically approved drugs. The calibrated machine learning algorithm was then used to predict the likelihood of existing drugs and known EVOO phytochemicals to be similar in action to the drugs impacting AD protein networks. These analyses identified the following ten EVOO phytochemicals with the highest likelihood of being active against AD: quercetin, genistein, luteolin, palmitoleate, stearic acid, apigenin, epicatechin, kaempferol, squalene, and daidzein (in the order from the highest to the lowest likelihood). This in silico study presents a framework that brings together artificial intelligence, analytical chemistry, and omics studies to identify unique therapeutic agents. It provides new insights into how EVOO constituents may help treat or prevent AD and potentially provide a basis for consideration in future clinical studies.

OBJECTIVE: If lesions in sporadic amyotrophic lateral sclerosis (ALS) originate from a single focal onset site and spread contiguously by prion-like cell-to-cell propagation at a constant speed, the lesion spread time should be proportional to the anatomical distance. We verify this model in the patients.
METHODS: In 29 sporadic ALS patients with hand onset followed by spread to shoulder and leg, we retrospectively evaluated the inter/intra-regional spread time ratio: time interval of symptoms from hand-to-leg divided by that from hand-to-shoulder. We also obtained the corresponding inter-/intra-regional distance ratios of spinal cord from magnetic resonance imaging of 12 patients, and those of primary motor cortex from coordinates using neuroimaging software.
RESULTS: Inter-/intra-regional spread time ratios ranged from 0.29 to 6.00 (median 1.20). Distance ratios ranged from 1.85 to 2.86 in primary motor cortex and from 5.79 to 8.67 in spinal cord. Taken together with clinical manifestations, of 27 patients with the requisite information available, lesion spreading was consistent with the model in primary motor cortex in 4 (14.8%) patients, and in spinal cord in only 1 (3.7%) patient. However, in more patients (12 of 29 patients: 41.4%), the inter-regional spread times in a long anatomical distance of hand-to-leg were shorter than or equal to the intra-regional spread times in a short anatomical distance of hand-to-shoulder.
CONCLUSIONS: Contiguous cell-to-cell propagation at a constant speed might not play a major role at least in distant lesion spreading of ALS. Several mechanisms can be responsible for progression in ALS.