A full-proteome, interaction-specific characterization of mutational hotspots across human cancers

Rapid accumulation of cancer genomic data has led to the identification of an increasing number of mutational hotspots with uncharacterized significance. Here we present a biologically informed computational framework that characterizes the functional relevance of all 1107 published mutational hotspots identified in approximately 25,000 tumor samples across 41 cancer types in the context of a human 3D interactome network, in which the interface of each interaction is mapped at residue resolution. Hotspots reside in network hub proteins and are enriched on protein interaction interfaces, suggesting that alteration of specific protein–protein interactions is critical for the oncogenicity of many hotspot mutations. Our framework enables, for the first time, systematic identification of specific protein interactions affected by hotspot mutations at the full proteome scale. Furthermore, by constructing a hotspot-affected network that connects all hotspot-affected interactions throughout the whole-human interactome, we uncover genome-wide relationships among hotspots and implicate novel cancer proteins that do not harbor hotspot mutations themselves. Moreover, applying our network-based framework to specific cancer types identifies clinically significant hotspots that can be used for prognosis and therapy targets. Overall, we show that our framework bridges the gap between the statistical significance of mutational hotspots and their biological and clinical significance in human cancers.

Detection of melanogenesis- and anti-apoptosis-associated melanoma factors: Array CGH and PPI mapping integrating study

Background: We investigated melanogenesis- and anti-apoptosis-related melanoma factors in melanoma cells (TXM1, TXM18, A375P, and A375SM). Objective: To find melanoma associated hub factor, high-throughput screening-based techniques integrating with bioinformatics were investigated. Methods: Array CGH analysis was conducted with a commercial system. Total genomic DNAs prepared individually from each cell line with control DNA were properly labeled with Cy3-dCTP and Cy5-dCTP and hybridizations and subsequently performed data treatment by the log2 green (G; test) to red (R; reference) fluorescence ratios (G/R). Gain or loss of copy number was judged by spots with log2-transformed ratios. PPI mapping analysis of detected candidate genes based on the array CGH results was conducted using the human interactome in the STRING database. Energy minimization and a short molecular dynamics (MD) simulation using the implicit solvation model in CHARMM were performed to analyze the interacting residues between YWHAZ and YWHAB. Results: Three genes (BMP-4, BFGF, LEF-1) known to be involved in melanogenesis were found to lose chromosomal copy numbers, and Chr. 6q23.3 was lost in all tested cell lines. Ten hub genes (CTNNB1, PEX13, PEX14, PEX5, IFNG, EXOSC3, EXOSC1, EXOSC8, UBC, and PEX10) were predicted to be functional interaction factors in the network of the 6q23.3 locus. The apoptosis-associated genes E2F1, p50, BCL2L1, and BIRC7 gained, and FGF2 lost chromosomal copy numbers in the tested melanoma cell lines. YWHAB, which gained chromosomal copy numbers, was predicted to be the most important hub protein in melanoma cells. Molecular dynamics simulations for binding YWHAB and YWHAZ were conducted, and the complex was predicted to be energetically and structurally stable through its 3 hydrogen-bond patterns. The number of interacting residues is 27. Conclusion: Our study compares genome-wide screening interactomics predictions for melanoma factors and offers new information for understanding melanogenesis- and anti-apoptosis-associated mechanisms in melanoma. Especially, YWHAB was newly detected as a core factor in melanoma cells.

MOSES: A New Approach to Integrate Interactome Topology and Functional Features for Disease Gene Prediction

Disease gene prediction is to date one of the main computational challenges of precision medicine. It is still uncertain if disease genes have unique functional properties that distinguish them from other non-disease genes or, from a network perspective, if they are located randomly in the interactome or show specific patterns in the network topology. In this study, we propose a new method for disease gene prediction based on the use of biological knowledge-bases (gene-disease associations, genes functional annotations, etc.) and interactome network topology. The proposed algorithm called MOSES is based on the definition of two somewhat opposing sets of genes both disease-specific from different perspectives: warm seeds (i.e., disease genes obtained from databases) and cold seeds (genes far from the disease genes on the interactome and not involved in their biological functions). The application of MOSES to a set of 40 diseases showed that the suggested putative disease genes are significantly enriched in their reference disease. Reassuringly, known and predicted disease genes together, tend to form a connected network module on the human interactome, mitigating the scattered distribution of disease genes which is probably due to both the paucity of disease-gene associations and the incompleteness of the interactome.

Quantitative fragmentomics allow affinity mapping of interactomes

Human protein networks have been widely explored but most binding affinities remain unknown, hindering quantitative interactome-function studies. Yet interactomes rely on minimal interacting fragments displaying quantifiable affinities. Here we measured the affinities of 65,000 interactions involving PDZ domains and their target PDZ-binding motifs (PBM) within a human interactome region particularly relevant for viral infection and cancer. We calculate interactomic distances, identify hot spots for viral interference, generate binding profiles and specificity logos, and explain selected cases by crystallographic studies. Mass spectrometry experiments on cell extracts and literature surveys show that quantitative fragmentomics effectively complement protein interactomics by providing affinities and completeness of coverage, putting a full human interactome affinity survey within realistic reach. Finally, we show that interactome hijacking by the viral PBM of human papillomavirus (HPV) E6 oncoprotein deeply impacts the host cell proteome way beyond immediate E6 binders, illustrating the complex system-wide relationship between interactome and function.

Network neighbors of viral targets and differentially expressed genes in COVID-19 are drug target candidates

The COVID-19 pandemic is raging. It revealed the importance of rapid scientific advancement towards understanding and treating new diseases. To address this challenge, we adapt an explainable artificial intelligence algorithm for data fusion and utilize it on new omics data on viral–host interactions, human protein interactions, and drugs to better understand SARS-CoV-2 infection mechanisms and predict new drug–target interactions for COVID-19. We discover that in the human interactome, the human proteins targeted by SARS-CoV-2 proteins and the genes that are differentially expressed after the infection have common neighbors central in the interactome that may be key to the disease mechanisms. We uncover 185 new drug–target interactions targeting 49 of these key genes and suggest re-purposing of 149 FDA-approved drugs, including drugs targeting VEGF and nitric oxide signaling, whose pathways coincide with the observed COVID-19 symptoms. Our integrative methodology is universal and can enable insight into this and other serious diseases.

Assessment of community efforts to advance computational prediction of protein-protein interactions

Comprehensive insights from the human protein-protein interaction (PPI) network, known as the human interactome, can provide important insights into the molecular mechanisms of complex biological processes and diseases. Despite the remarkable experimental efforts undertaken to date to determine the structure of the human interactome, many PPIs remain unmapped. Computational approaches, especially network-based methods, can facilitate the identification of new PPIs. Many such approaches have been proposed. However, a systematic evaluation of existing network-based methods in predicting PPIs is still lacking. Here, we report community efforts initiated by the International Network Medicine Consortium to benchmark the ability of 24 representative network-based methods to predict PPIs across five different interactomes, including a synthetic interactome generated by the duplication-mutation-complementation model, and the interactomes of four different organisms: A. thaliana, C. elegans, S. cerevisiae, and H. sapiens. We selected the top-seven methods through a computational validation on the human interactome. We next experimentally validated their top-500 predicted PPIs (in total 3,276 predicted PPIs) using the yeast two-hybrid assay, finding 1,177 new human PPIs (involving 633 proteins). Our results indicate that task-tailored similarity-based methods, which leverage the underlying network characteristics of PPIs, show superior performance over other general link prediction methods. Through experimental validation, we confirmed that the top-ranking methods show promising performance externally. For example, from the top 500 PPIs predicted by an advanced similarity-base method [MPS(B&T)], 430 were successfully tested by Y2H with 376 testing positive, yielding a precision of 87.4%. These results establish advanced similarity-based methods as powerful tools for the prediction of human PPIs.

Algorithmically Deduced FREM2 Molecular Pathway Is a Potent Grade and Survival Biomarker of Human Gliomas

Gliomas are the most common malignant brain tumors with high mortality rates. Recently we showed that the FREM2 gene has a role in glioblastoma progression. Here we reconstructed the FREM2 molecular pathway using the human interactome model. We assessed the biomarker capacity of FREM2 expression and its pathway as the overall survival (OS) and progression-free survival (PFS) biomarkers. To this end, we used three literature and one experimental RNA sequencing datasets collectively covering 566 glioblastomas (GBM) and 1097 low-grade gliomas (LGG). The activation level of deduced FREM2 pathway showed strong biomarker characteristics and significantly outperformed the FREM2 expression level itself. For all relevant datasets, it could robustly discriminate GBM and LGG (p < 1.63 × 10−13, AUC > 0.74). High FREM2 pathway activation level was associated with poor OS in LGG (p < 0.001), and low PFS in LGG (p < 0.001) and GBM (p < 0.05). FREM2 pathway activation level was poor prognosis biomarker for OS (p < 0.05) and PFS (p < 0.05) in LGG with IDH mutation, for PFS in LGG with wild type IDH (p < 0.001) and mutant IDH with 1p/19q codeletion(p < 0.05), in GBM with unmethylated MGMT (p < 0.05), and in GBM with wild type IDH (p < 0.05). Thus, we conclude that the activation level of the FREM2 pathway is a potent new-generation diagnostic and prognostic biomarker for multiple molecular subtypes of GBM and LGG.

Predicting anticancer hyperfoods with graph convolutional networks

Background Recent efforts in the field of nutritional science have allowed the discovery of disease-beating molecules within foods based on the commonality of bioactive food molecules to FDA-approved drugs. The pioneering work in this field used an unsupervised network propagation algorithm to learn the systemic-wide effect on the human interactome of 1962 FDA-approved drugs and a supervised algorithm to predict anticancer therapeutics using the learned representations. Then, a set of bioactive molecules within foods was fed into the model, which predicted molecules with cancer-beating potential.The employed methodology consisted of disjoint unsupervised feature generation and classification tasks, which can result in sub-optimal learned drug representations with respect to the classification task. Additionally, due to the disjoint nature of the tasks, the employed approach proved cumbersome to optimize, requiring testing of thousands of hyperparameter combinations and significant computational resources.To overcome the technical limitations highlighted above, we represent each drug as a graph (human interactome) with its targets as binary node features on the graph and formulate the problem as a graph classification task. To solve this task, inspired by the success of graph neural networks in graph classification problems, we use an end-to-end graph neural network model operating directly on the graphs, which learns drug representations to optimize model performance in the prediction of anticancer therapeutics. Results The proposed model outperforms the baseline approach in the anticancer therapeutic prediction task, achieving an F1 score of 67.99%±2.52% and an AUPR of 73.91%±3.49%. It is also shown that the model is able to capture knowledge of biological pathways to predict anticancer molecules based on the molecules’ effects on cancer-related pathways. Conclusions We introduce an end-to-end graph convolutional model to predict cancer-beating molecules within food. The introduced model outperforms the existing baseline approach, and shows interpretability, paving the way to the future of a personalized nutritional science approach allowing the development of nutrition strategies for cancer prevention and/or therapeutics.