Edgotype: a fundamental link between genotype and phenotype

Uncovering cell-type-specific immunomodulatory variants and molecular phenotypes in COVID-19 using structurally resolved protein networks

Immunomodulatory variants that lead to the loss or gain of specific protein interactions often manifest only as organismal phenotypes in infectious disease. Here, we propose a network-based approach to integrate genetic variation with a structurally resolved human protein interactome network to prioritize immunomodulatory variants in COVID-19. We find that, in addition to variants that pass genome-wide significance thresholds, variants at the interface of specific protein-protein interactions, even though they do not meet genome-wide thresholds, are equally immunomodulatory. The integration of these variants with single-cell epigenomic and transcriptomic data prioritizes myeloid and T cell subsets as the most affected by these variants across both the peripheral blood and the lung compartments. Of particular interest is a common coding variant that disrupts the OAS1-PRMT6 interaction and affects downstream interferon signaling. Critically, our framework is generalizable across infectious disease contexts and can be used to implicate immunomodulatory variants that do not meet genome-wide significance thresholds.

Pervasive mislocalization of pathogenic coding variants underlying human disorders

Widespread sequencing has yielded thousands of missense variants predicted or confirmed as disease causing. This creates a new bottleneck: determining the functional impact of each variant-typically a painstaking, customized process undertaken one or a few genes and variants at a time. Here, we established a high-throughput imaging platform to assay the impact of coding variation on protein localization, evaluating 3,448 missense variants of over 1,000 genes and phenotypes. We discovered that mislocalization is a common consequence of coding variation, affecting about one-sixth of all pathogenic missense variants, all cellular compartments, and recessive and dominant disorders alike. Mislocalization is primarily driven by effects on protein stability and membrane insertion rather than disruptions of trafficking signals or specific interactions. Furthermore, mislocalization patterns help explain pleiotropy and disease severity and provide insights on variants of uncertain significance. Our publicly available resource extends our understanding of coding variation in human diseases.

Ubiquitin: Not just a one-way ticket to the proteasome, but a therapeutic dial to fine-tune the molecular landscape of disease

Recently, there is a rise in studies that recognize the importance of targeting ubiquitin and related molecular machinery in various therapeutic contexts. Here we briefly discuss the history of ubiquitin, its biological roles in protein degradation and beyond, as well as the current state of ubiquitin-targeting therapeutics across diseases. We conclude that targeting ubiquitin machinery is approaching a renaissance, and tapping its full potential will require embracing a wholistic perspective of ubiquitin's multifaceted roles.

Sliding Window INteraction Grammar (SWING): a generalized interaction language model for peptide and protein interactions

The explosion of sequence data has allowed the rapid growth of protein language models (pLMs). pLMs have now been employed in many frameworks including variant-effect and peptide-specificity prediction. Traditionally, for protein-protein or peptide-protein interactions (PPIs), corresponding sequences are either co-embedded followed by post-hoc integration or the sequences are concatenated prior to embedding. Interestingly, no method utilizes a language representation of the interaction itself. We developed an interaction LM (iLM), which uses a novel language to represent interactions between protein/peptide sequences. Sliding Window Interaction Grammar (SWING) leverages differences in amino acid properties to generate an interaction vocabulary. This vocabulary is the input into a LM followed by a supervised prediction step where the LM's representations are used as features. SWING was first applied to predicting peptide:MHC (pMHC) interactions. SWING was not only successful at generating Class I and Class II models that have comparable prediction to state-of-the-art approaches, but the unique Mixed Class model was also successful at jointly predicting both classes. Further, the SWING model trained only on Class I alleles was predictive for Class II, a complex prediction task not attempted by any existing approach. For de novo data, using only Class I or Class II data, SWING also accurately predicted Class II pMHC interactions in murine models of SLE (MRL/lpr model) and T1D (NOD model), that were validated experimentally. To further evaluate SWING's generalizability, we tested its ability to predict the disruption of specific protein-protein interactions by missense mutations. Although modern methods like AlphaMissense and ESM1b can predict interfaces and variant effects/pathogenicity per mutation, they are unable to predict interaction-specific disruptions. SWING was successful at accurately predicting the impact of both Mendelian mutations and population variants on PPIs. This is the first generalizable approach that can accurately predict interaction-specific disruptions by missense mutations with only sequence information. Overall, SWING is a first-in-class generalizable zero-shot iLM that learns the language of PPIs.

New candidate genes potentially involved in Zika virus teratogenesis

Despite efforts to elucidate Zika virus (ZIKV) teratogenesis, still several issues remain unresolved, particularly on the molecular mechanisms behind the pathogenesis of Congenital Zika Syndrome (CZS). To answer this question, we used bioinformatics tools, animal experiments and human gene expression analysis to investigate genes related to brain development potentially involved in CZS. Searches in databases for genes related to brain development and CZS were performed, and a protein interaction network was created. The expression of these genes was analyzed in a CZS animal model and secondary gene expression analysis (DGE) was performed in human cells exposed to ZIKV. A total of 2610 genes were identified in the databases, of which 1013 were connected. By applying centrality statistics of the global network, 36 candidate genes were identified, which, after selection resulted in nine genes. Gene expression analysis revealed distinctive expression patterns for PRKDC, PCNA, ATM, SMC3 as well as for FGF8 and SHH in the CZS model. Furthermore, DGE analysis altered expression of ATM, PRKDC, PCNA. In conclusion, systems biology are helpful tools to identify candidate genes to be validated in vitro and in vivo. PRKDC, PCNA, ATM, SMC3, FGF8 and SHH have altered expression in ZIKV-induced brain malformations.

Development and validation of a new genotype-phenotype correlation for Niemann-Pick disease type C1

Hundreds of NPC1 variants cause highly heterogeneous phenotypes. This study aims to explore the genotype-phenotype correlation of NPC1, especially for missense variants. In a well-characterized cohort, phenotypes are graded into three clinical forms: mild, intermediate, and severe. Missense residue structural location was stratified into three categories: surface, partially, and fully buried. The association of phenotypes with the topography of the amino acid substitution in the protein structure was investigated in our cohort and validated in two reported cohorts. One hundred six unrelated NPC1 patients were enrolled. A significant correlation of genotype-phenotype was found in 81 classified individuals with two or one (the second was null variant) missense variant (p < 0.001): of 25 patients with at least one missense variant of surface (group A), 19 (76%) mild, six (24%) intermediate, and none severe; of 31 cases with at least one missense variant of partially buried without surface variants (group B), 11 (35%) mild, 16 (52%) intermediate, and four (13%) severe; of the remaining 25 patients with two or one buried missense variants (group C), eight (32%) mild, nine (36%) intermediate, and eight (32%) severe. Additionally, 7-ketocholesterol, the biomarker, was lower in group A than in group B (p = 0.024) and group C (p = 0.029). A model was proposed that accurately predicted phenotypes of 72 of 90 (80%), 73 of85 (86%), and 64 of 69 (93%) patients in our cohort, Italian, and UK cohort, respectively. This study proposed a novel genotype-phenotype correlation in NPC1, linking the underlying molecular pathophysiology with clinical phenotype and aiding genetic counseling and evaluation in clinical practice.

The Extent of Edgetic Perturbations in the Human Interactome Caused by Population-Specific Mutations

Until recently, efforts in population genetics have been focused primarily on people of European ancestry. To attenuate this bias, global population studies, such as the 1000 Genomes Project, have revealed differences in genetic variation across ethnic groups. How many of these differences can be attributed to population-specific traits? To answer this question, the mutation data must be linked with functional outcomes. A new "edgotype" concept has been proposed, which emphasizes the interaction-specific, "edgetic", perturbations caused by mutations in the interacting proteins. In this work, we performed systematic in silico edgetic profiling of ~50,000 non-synonymous SNVs (nsSNVs) from the 1000 Genomes Project by leveraging our semi-supervised learning approach SNP-IN tool on a comprehensive set of over 10,000 protein interaction complexes. We interrogated the functional roles of the variants and their impact on the human interactome and compared the results with the pathogenic variants disrupting PPIs in the same interactome. Our results demonstrated that a considerable number of nsSNVs from healthy populations could rewire the interactome. We also showed that the proteins enriched with interaction-disrupting mutations were associated with diverse functions and had implications in a broad spectrum of diseases. Further analysis indicated that distinct gene edgetic profiles among major populations could shed light on the molecular mechanisms behind the population phenotypic variances. Finally, the network analysis revealed that the disease-associated modules surprisingly harbored a higher density of interaction-disrupting mutations from healthy populations. The variation in the cumulative network damage within these modules could potentially account for the observed disparities in disease susceptibility, which are distinctly specific to certain populations. Our work demonstrates the feasibility of a large-scale in silico edgetic study, and reveals insights into the orchestrated play of population-specific mutations in the human interactome.

Neurolysin Knockout Mice in a Diet-Induced Obesity Model

Neurolysin oligopeptidase (E.C.3.4.24.16; Nln), a member of the zinc metallopeptidase M3 family, was first identified in rat brain synaptic membranes hydrolyzing neurotensin at the Pro-Tyr peptide bond. The previous development of C57BL6/N mice with suppression of Nln gene expression (Nln-/-), demonstrated the biological relevance of this oligopeptidase for insulin signaling and glucose uptake. Here, several metabolic parameters were investigated in Nln-/- and wild-type C57BL6/N animals (WT; n = 5-8), male and female, fed either a standard (SD) or a hypercaloric diet (HD), for seven weeks. Higher food intake and body mass gain was observed for Nln-/- animals fed HD, compared to both male and female WT control animals fed HD. Leptin gene expression was higher in Nln-/- male and female animals fed HD, compared to WT controls. Both WT and Nln-/- females fed HD showed similar gene expression increase of dipeptidyl peptidase 4 (DPP4), a peptidase related to glucagon-like peptide-1 (GLP-1) metabolism. The present data suggest that Nln participates in the physiological mechanisms related to diet-induced obesity. Further studies will be necessary to better understand the molecular mechanism responsible for the higher body mass gain observed in Nln-/- animals fed HD.

Latent space search based multimodal optimization with personalized edge-network biomarker for multi-purpose early disease prediction

Considering that cancer is resulting from the comutation of several essential genes of individual patients, researchers have begun to focus on identifying personalized edge-network biomarkers (PEBs) using personalized edge-network analysis for clinical practice. However, most of existing methods ignored the optimization of PEBs when multimodal biomarkers exist in multi-purpose early disease prediction (MPEDP). To solve this problem, this study proposes a novel model (MMPDENB-RBM) that combines personalized dynamic edge-network biomarkers (PDENB) theory, multimodal optimization strategy and latent space search scheme to identify biomarkers with different configurations of PDENB modules (i.e. to effectively identify multimodal PDENBs). The application to the three largest cancer omics datasets from The Cancer Genome Atlas database (i.e. breast invasive carcinoma, lung squamous cell carcinoma and lung adenocarcinoma) showed that the MMPDENB-RBM model could more effectively predict critical cancer state compared with other advanced methods. And, our model had better convergence, diversity and multimodal property as well as effective optimization ability compared with the other state-of-art methods. Particularly, multimodal PDENBs identified were more enriched with different functional biomarkers simultaneously, such as tissue-specific synthetic lethality edge-biomarkers including cancer driver genes and disease marker genes. Importantly, as our aim, these multimodal biomarkers can perform diverse biological and biomedical significances for drug target screen, survival risk assessment and novel biomedical sight as the expected multi-purpose of personalized early disease prediction. In summary, the present study provides multimodal property of PDENBs, especially the therapeutic biomarkers with more biological significances, which can help with MPEDP of individual cancer patients.

Pervasive mislocalization of pathogenic coding variants underlying human disorders

Widespread sequencing has yielded thousands of missense variants predicted or confirmed as disease-causing. This creates a new bottleneck: determining the functional impact of each variant - largely a painstaking, customized process undertaken one or a few genes or variants at a time. Here, we established a high-throughput imaging platform to assay the impact of coding variation on protein localization, evaluating 3,547 missense variants of over 1,000 genes and phenotypes. We discovered that mislocalization is a common consequence of coding variation, affecting about one-sixth of all pathogenic missense variants, all cellular compartments, and recessive and dominant disorders alike. Mislocalization is primarily driven by effects on protein stability and membrane insertion rather than disruptions of trafficking signals or specific interactions. Furthermore, mislocalization patterns help explain pleiotropy and disease severity and provide insights on variants of unknown significance. Our publicly available resource will likely accelerate the understanding of coding variation in human diseases.

Phenotypic heterogeneity in human genetic diseases: ultrasensitivity-mediated threshold effects as a unifying molecular mechanism

Phenotypic heterogeneity is very common in genetic systems and in human diseases and has important consequences for disease diagnosis and treatment. In addition to the many genetic and non-genetic (e.g., epigenetic, environmental) factors reported to account for part of the heterogeneity, we stress the importance of stochastic fluctuation and regulatory network topology in contributing to phenotypic heterogeneity. We argue that a threshold effect is a unifying principle to explain the phenomenon; that ultrasensitivity is the molecular mechanism for this threshold effect; and discuss the three conditions for phenotypic heterogeneity to occur. We suggest that threshold effects occur not only at the cellular level, but also at the organ level. We stress the importance of context-dependence and its relationship to pleiotropy and edgetic mutations. Based on this model, we provide practical strategies to study human genetic diseases. By understanding the network mechanism for ultrasensitivity and identifying the critical factor, we may manipulate the weak spot to gently nudge the system from an ultrasensitive state to a stable non-disease state. Our analysis provides a new insight into the prevention and treatment of genetic diseases.

Homo cerevisiae-Leveraging Yeast for Investigating Protein-Protein Interactions and Their Role in Human Disease

Understanding how genetic variation affects phenotypes represents a major challenge, particularly in the context of human disease. Although numerous disease-associated genes have been identified, the clinical significance of most human variants remains unknown. Despite unparalleled advances in genomics, functional assays often lack sufficient throughput, hindering efficient variant functionalization. There is a critical need for the development of more potent, high-throughput methods for characterizing human genetic variants. Here, we review how yeast helps tackle this challenge, both as a valuable model organism and as an experimental tool for investigating the molecular basis of phenotypic perturbation upon genetic variation. In systems biology, yeast has played a pivotal role as a highly scalable platform which has allowed us to gain extensive genetic and molecular knowledge, including the construction of comprehensive interactome maps at the proteome scale for various organisms. By leveraging interactome networks, one can view biology from a systems perspective, unravel the molecular mechanisms underlying genetic diseases, and identify therapeutic targets. The use of yeast to assess the molecular impacts of genetic variants, including those associated with viral interactions, cancer, and rare and complex diseases, has the potential to bridge the gap between genotype and phenotype, opening the door for precision medicine approaches and therapeutic development.

AD-Syn-Net: systematic identification of Alzheimer's disease-associated mutation and co-mutation vulnerabilities via deep learning

Alzheimer's disease (AD) is one of the most challenging neurodegenerative diseases because of its complicated and progressive mechanisms, and multiple risk factors. Increasing research evidence demonstrates that genetics may be a key factor responsible for the occurrence of the disease. Although previous reports identified quite a few AD-associated genes, they were mostly limited owing to patient sample size and selection bias. There is a lack of comprehensive research aimed to identify AD-associated risk mutations systematically. To address this challenge, we hereby construct a large-scale AD mutation and co-mutation framework ('AD-Syn-Net'), and propose deep learning models named Deep-SMCI and Deep-CMCI configured with fully connected layers that are capable of predicting cognitive impairment of subjects effectively based on genetic mutation and co-mutation profiles. Next, we apply the customized frameworks to data sets to evaluate the importance scores of the mutations and identified mutation effectors and co-mutation combination vulnerabilities contributing to cognitive impairment. Furthermore, we evaluate the influence of mutation pairs on the network architecture to dissect the genetic organization of AD and identify novel co-mutations that could be responsible for dementia, laying a solid foundation for proposing future targeted therapy for AD precision medicine. Our deep learning model codes are available open access here: https://github.com/Pan-Bio/AD-mutation-effectors.

Targeting stressor-induced dysfunctions in protein-protein interaction networks via epichaperomes

Diseases are manifestations of complex changes in protein-protein interaction (PPI) networks whereby stressors, genetic, environmental, and combinations thereof, alter molecular interactions and perturb the individual from the level of cells and tissues to the entire organism. Targeting stressor-induced dysfunctions in PPI networks has therefore become a promising but technically challenging frontier in therapeutics discovery. This opinion provides a new framework based upon disrupting epichaperomes - pathological entities that enable dysfunctional rewiring of PPI networks - as a mechanism to revert context-specific PPI network dysfunction to a normative state. We speculate on the implications of recent research in this area for a precision medicine approach to detecting and treating complex diseases, including cancer and neurodegenerative disorders.

Construction and analysis of sample-specific driver modules for breast cancer

BACKGROUND

It is important to understand the functional impact of somatic mutation and methylation aberration at an individual level to implement precision medicine. Recent studies have demonstrated that the perturbation of gene interaction networks can provide a fundamental link between genotype (or epigenotype) and phenotype. However, it is unclear how individual mutations affect the function of biological networks, especially for individual methylation aberration. To solve this, we provided a sample-specific driver module construction method using the 2-order network theory and hub-gene theory to identify individual perturbation networks driven by mutations or methylation aberrations.

RESULTS

Our method integrated multi-omics of breast cancer, including genomics, transcriptomics, epigenomics and interactomics, and provided new insight into the synergistic collaboration between methylation and mutation at an individual level. A common driver pattern of breast cancer was identified from a novel perspective of a driver module, which is correlated to the occurrence and development of breast cancer. The constructed driver module reflects the survival prognosis and degree of malignancy among different subtypes of breast cancer. Additionally, subtype-specific driver modules were identified.

CONCLUSIONS

This study explores the driver module of individual cancer, and contributes to a better understanding of the mechanism of breast cancer driven by the mutations and methylation variations from the point of view of the driver network. This work will help identify new therapeutic combinations of gene mutations and drugs in humans.

Systems Bioinformatics Reveals Possible Relationship between COVID-19 and the Development of Neurological Diseases and Neuropsychiatric Disorders

Coronavirus Disease 2019 (COVID-19) is associated with increased incidence of neurological diseases and neuropsychiatric disorders after infection, but how it contributes to their development remains under investigation. Here, we investigate the possible relationship between COVID-19 and the development of ten neurological disorders and three neuropsychiatric disorders by exploring two pathological mechanisms: (i) dysregulation of host biological processes via virus-host protein-protein interactions (PPIs), and (ii) autoreactivity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epitopes with host "self" proteins via molecular mimicry. We also identify potential genetic risk factors which in combination with SARS-CoV-2 infection might lead to disease development. Our analysis indicated that neurodegenerative diseases (NDs) have a higher number of disease-associated biological processes that can be modulated by SARS-CoV-2 via virus-host PPIs than neuropsychiatric disorders. The sequence similarity analysis indicated the presence of several matching 5-mer and/or 6-mer linear motifs between SARS-CoV-2 epitopes with autoreactive epitopes found in Alzheimer's Disease (AD), Parkinson's Disease (PD), Myasthenia Gravis (MG) and Multiple Sclerosis (MS). The results include autoreactive epitopes that recognize amyloid-beta precursor protein (APP), microtubule-associated protein tau (MAPT), acetylcholine receptors, glial fibrillary acidic protein (GFAP), neurofilament light polypeptide (NfL) and major myelin proteins. Altogether, our results suggest that there might be an increased risk for the development of NDs after COVID-19 both via autoreactivity and virus-host PPIs.

Using transcriptomics to predict and visualize disease status in bighorn sheep (Ovis canadensis)

Increasing risk of pathogen spillover coupled with overall declines in wildlife population abundance in the Anthropocene make infectious disease a relevant concern for species conservation worldwide. While emerging molecular tools could improve our diagnostic capabilities and give insight into mechanisms underlying wildlife disease risk, they have rarely been applied in practice. Here, employing a previously reported gene transcription panel of common immune markers to track physiological changes, we present a detailed analysis over the course of both acute and chronic infection in one wildlife species where disease plays a critical role in conservation, bighorn sheep (Ovis canadensis). Differential gene transcription patterns distinguished between infection statuses over the course of acute infection and differential correlation (DC) analyses identified clear changes in gene co-transcription patterns over the early stages of infection, with transcription of four genes-TGFb, AHR, IL1b and MX1-continuing to increase even as transcription of other immune-associated genes waned. In a separate analysis, we considered the capacity of the same gene transcription panel to aid in differentiating between chronically infected animals and animals in other disease states outside of acute disease events (an immediate priority for wildlife management in this system). We found that this transcription panel was capable of accurately identifying chronically infected animals in the test dataset, though additional data will be required to determine how far this ability extends. Taken together, our results showcase the successful proof of concept and breadth of potential utilities that gene transcription might provide to wildlife disease management, from direct insight into mechanisms associated with differential disease response to improved diagnostic capacity in the field.

Predicting functional consequences of mutations using molecular interaction network features

Variant interpretation remains a central challenge for precision medicine. Missense variants are particularly difficult to understand as they change only a single amino acid in a protein sequence yet can have large and varied effects on protein activity. Numerous tools have been developed to identify missense variants with putative disease consequences from protein sequence and structure. However, biological function arises through higher order interactions among proteins and molecules within cells. We therefore sought to capture information about the potential of missense mutations to perturb protein interaction networks by integrating protein structure and interaction data. We developed 16 network-based annotations for missense mutations that provide orthogonal information to features classically used to prioritize variants. We then evaluated them in the context of a proven machine-learning framework for variant effect prediction across multiple benchmark datasets to demonstrate their potential to improve variant classification. Interestingly, network features resulted in larger performance gains for classifying somatic mutations than for germline variants, possibly due to different constraints on what mutations are tolerated at the cellular versus organismal level. Our results suggest that modeling variant potential to perturb context-specific interactome networks is a fruitful strategy to advance in silico variant effect prediction.

CancerVar: An artificial intelligence-empowered platform for clinical interpretation of somatic mutations in cancer

Several knowledgebases are manually curated to support clinical interpretations of thousands of hotspot somatic mutations in cancer. However, discrepancies or even conflicting interpretations are observed among these databases. Furthermore, many previously undocumented mutations may have clinical or functional impacts on cancer but are not systematically interpreted by existing knowledgebases. To address these challenges, we developed CancerVar to facilitate automated and standardized interpretations for 13 million somatic mutations based on the AMP/ASCO/CAP 2017 guidelines. We further introduced a deep learning framework to predict oncogenicity for these variants using both functional and clinical features. CancerVar achieved satisfactory performance when compared to several independent knowledgebases and, using clinically curated datasets, demonstrated practical utility in classifying somatic variants. In summary, by integrating clinical guidelines with a deep learning framework, CancerVar facilitates clinical interpretation of somatic variants, reduces manual work, improves consistency in variant classification, and promotes implementation of the guidelines.

Are transient protein-protein interactions more dispensable?

Protein-protein interactions (PPIs) are key drivers of cell function and evolution. While it is widely assumed that most permanent PPIs are important for cellular function, it remains unclear whether transient PPIs are equally important. Here, we estimate and compare dispensable content among transient PPIs and permanent PPIs in human. Starting with a human reference interactome mapped by experiments, we construct a human structural interactome by building three-dimensional structural models for PPIs, and then distinguish transient PPIs from permanent PPIs using several structural and biophysical properties. We map common mutations from healthy individuals and disease-causing mutations onto the structural interactome, and perform structure-based calculations of the probabilities for common mutations (assumed to be neutral) and disease mutations (assumed to be mildly deleterious) to disrupt transient PPIs and permanent PPIs. Using Bayes' theorem we estimate that a similarly small fraction (<~20%) of both transient and permanent PPIs are completely dispensable, i.e., effectively neutral upon disruption. Hence, transient and permanent interactions are subject to similarly strong selective constraints in the human interactome.

Mapping oncogenic protein interactions for precision medicine

Normal protein‐protein interactions (normPPIs) occur with high fidelity to regulate almost every physiological process. In cancer, this highly organised and precisely regulated network is disrupted, hijacked or reprogrammed resulting in oncogenic protein‐protein interactions (oncoPPIs). OncoPPIs, which can result from genomic alterations, are a hallmark of many types of cancers. Recent technological advances in the field of mass spectrometry (MS)‐based interactomics, structural biology and drug discovery have prompted scientists to identify and characterise oncoPPIs. Disruption of oncoPPI interfaces has become a major focus of drug discovery programs and has resulted in the use of PPI‐specific drugs clinically. However, due to several technical hurdles, studies to build a reference oncoPPI map for various cancer types have not been undertaken. Therefore, there is an urgent need for experimental workflows to overcome the existing challenges in studying oncoPPIs in various cancers and to build comprehensive reference maps. Here, we discuss the important hurdles for characterising oncoPPIs and propose a three‐phase multidisciplinary workflow to identify and characterise oncoPPIs. Systematic identification of cancer‐type‐specific oncogenic interactions will spur new opportunities for PPI‐focused drug discovery projects and precision medicine.

Effect of FKBP12-Derived Intracellular Peptides on Rapamycin-Induced FKBP-FRB Interaction and Autophagy

Intracellular peptides (InPeps) generated by proteasomes were previously suggested as putative natural regulators of protein-protein interactions (PPI). Here, the main aim was to investigate the intracellular effects of intracellular peptide VFDVELL (VFD7) and related peptides on PPI. The internalization of the peptides was achieved using a C-terminus covalently bound cell-penetrating peptide (cpp; YGRKKRRQRRR). The possible inhibition of PPI was investigated using a NanoBiT^® luciferase structural complementation reporter system, with a pair of plasmids vectors each encoding, simultaneously, either FK506-binding protein (FKBP) or FKBP-binding domain (FRB) of mechanistic target of rapamycin complex 1 (mTORC1). The interaction of FKBP-FRB within cells occurs under rapamycin induction. Results shown that rapamycin-induced interaction between FKBP-FRB within human embryonic kidney 293 (HEK293) cells was inhibited by VFD7-cpp (10-500 nM) and FDVELLYGRKKRRQRRR (VFD6-cpp; 1-500 nM); additional VFD7-cpp derivatives were either less or not effective in inhibiting FKBP-FRB interaction induced by rapamycin. Molecular dynamics simulations suggested that selected peptides, such as VFD7-cpp, VFD6-cpp, VFAVELLYGRKKKRRQRRR (VFA7-cpp), and VFEVELLYGRKKKRRQRRR (VFA7-cpp), bind to FKBP and to FRB protein surfaces. However, only VFD7-cpp and VFD6-cpp induced changes on FKBP structure, which could help with understanding their mechanism of PPI inhibition. InPeps extracted from HEK293 cells were found mainly associated with macromolecular components (i.e., proteins and/or nucleic acids), contributing to understanding InPeps' intracellular proteolytic stability and mechanism of action-inhibiting PPI within cells. In a model of cell death induced by hypoxia-reoxygenation, VFD6-cpp (1 µM) increased the viability of mouse embryonic fibroblasts cells (MEF) expressing mTORC1-regulated autophagy-related gene 5 (Atg5), but not in autophagy-deficient MEF cells lacking the expression of Atg5. These data suggest that VFD6-cpp could have therapeutic applications reducing undesired side effects of rapamycin long-term treatments. In summary, the present report provides further evidence that InPeps have biological significance and could be valuable tools for the rational design of therapeutic molecules targeting intracellular PPI.

Deep Mutational Scanning of Protein-Protein Interactions Between Partners Expressed from Their Endogenous Loci In Vivo

Deep mutational scanning (DMS) generates mutants of a protein of interest in a comprehensive manner. CRISPR-Cas9 technology enables large-scale genome editing with high efficiency. Using both DMS and CRISPR-Cas9 therefore allows us to investigate the effects of thousands of mutations inserted directly in the genome. Combined with protein-fragment complementation assay (PCA), which enables the quantitative measurement of protein-protein interactions (PPIs) in vivo, these methods allow for the systematic assessment of the effects of mutations on PPIs in living cells. Here, we describe a method leveraging DMS, CRISPR-Cas9, and PCA to study the effect of point mutations on PPIs mediated by protein domains in yeast.

OUP accepted manuscript

Barcode fusion genetics (BFG) utilizes deep sequencing to improve the throughput of protein–protein interaction (PPI) screening in pools. BFG has been implemented in Yeast two-hybrid (Y2H) screens (BFG-Y2H). While Y2H requires test protein pairs to localize in the nucleus for reporter reconstruction, dihydrofolate reductase protein-fragment complementation assay (DHFR-PCA) allows proteins to localize in broader subcellular contexts and proves to be largely orthogonal to Y2H. Here, we implemented BFG to DHFR-PCA (BFG-PCA). This plasmid-based system can leverage ORF collections across model organisms to perform comparative analysis, unlike the original DHFR-PCA that requires yeast genomic integration. The scalability and quality of BFG-PCA were demonstrated by screening human and yeast interactions for >11 000 bait-prey pairs. BFG-PCA showed high-sensitivity and high-specificity for capturing known interactions for both species. BFG-Y2H and BFG-PCA capture distinct sets of PPIs, which can partially be explained based on the domain orientation of the reporter tags. BFG-PCA is a high-throughput protein interaction technology to interrogate binary PPIs that exploits clone collections from any species of interest, expanding the scope of PPI assays.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.

Mutation Edgotype Drives Fitness Effect in Human

Missense mutations are known to perturb protein-protein interaction networks (known as interactome networks) in different ways. However, it remains unknown how different interactome perturbation patterns (“edgotypes”) impact organismal fitness. Here, we estimate the fitness effect of missense mutations with different interactome perturbation patterns in human, by calculating the fractions of neutral and deleterious mutations that do not disrupt PPIs (“quasi-wild-type”), or disrupt PPIs either by disrupting the binding interface (“edgetic”) or by disrupting overall protein stability (“quasi-null”). We first map pathogenic mutations and common non-pathogenic mutations onto homology-based three-dimensional structural models of proteins and protein-protein interactions in human. Next, we perform structure-based calculations to classify each mutation as either quasi-wild-type, edgetic, or quasi-null. Using our predicted as well as experimentally determined interactome perturbation patterns, we estimate that >∼40% of quasi-wild-type mutations are effectively neutral and the remaining are mostly mildly deleterious, that >∼75% of edgetic mutations are only mildly deleterious, and that up to ∼75% of quasi-null mutations may be strongly detrimental. These estimates are the first such estimates of fitness effect for different network perturbation patterns in any interactome. Our results suggest that while mutations that do not disrupt the interactome tend to be effectively neutral, the majority of human PPIs are under strong purifying selection and the stability of most human proteins is essential to human life.

Cross-Species and Human Inter-Tissue Network Analysis of Genes Implicated in Longevity and Aging Reveal Strong Support for Nutrient Sensing

Intensive research efforts have been undertaken to slow human aging and therefore potentially delay the onset of age-related diseases. These efforts have generated an enormous amount of high-throughput data covering different levels in the physiologic hierarchy, e.g., genetic, epigenetic, transcriptomic, proteomic, and metabolomic, etc. We gathered 15 independent sources of information about genes potentially involved in human longevity and lifespan (N = 5836) and subjected them to various integrated analyses. Many of these genes were initially identified in non-human species, and we investigated their orthologs in three non-human species [i.e., mice (N = 967), fruit fly (N = 449), and worm (N = 411)] for further analysis. We characterized experimentally determined protein-protein interaction networks (PPIN) involving each species' genes from 9 known protein databases and studied the enriched biological pathways among the individually constructed PPINs. We observed three important signaling pathways: FoxO signaling, mTOR signaling, and autophagy to be common and highly enriched in all four species (p-value ≤ 0.001). Our study implies that the interaction of proteins involved in the mechanistic target of rapamycin (mTOR) signaling pathway is somewhat limited to each species or that a "rewiring" of specific networks has taken place over time. To corroborate our findings, we repeated our analysis in 43 different human tissues. We investigated conserved modules in various tissue-specific PPINs of the longevity-associated genes based upon their protein expression. This analysis also revealed mTOR signaling as shared biological processes across four different human tissue-specific PPINs for liver, heart, skeletal muscle, and adipose tissue. Further, we explored our results' translational potential by assessing the protein interactions with all the reported drugs and compounds that have been experimentally verified to promote longevity in the three-comparator species. We observed that the target proteins of the FDA-approved drug rapamycin (a known inhibitor of mTOR) were conserved across all four species. Drugs like melatonin and metformin exhibited shared targets with rapamycin in the human PPIN. The detailed information about the curated gene list, cross-species orthologs, PPIN, and pathways was assembled in an interactive data visualization portal using RStudio's Shiny framework (https://agingnetwork.shinyapps.io/frontiers/).

Rationale and Criteria for a COVID-19 Model Framework

Complex systems are inherently multilevel and multiscale systems. The infectious disease system is considered a complex system resulting from the interaction between three sub-systems (host, pathogen, and environment) organized into a hierarchical structure, ranging from the cellular to the macro-ecosystem level, with multiscales. Therefore, to describe infectious disease phenomena that change through time and space and at different scales, we built a model framework where infectious disease must be considered the set of biological responses of human hosts to pathogens, with biological pathways shared with other pathologies in an ecological interaction context. In this paper, we aimed to design a framework for building a disease model for COVID-19 based on current literature evidence. The model was set up by identifying the molecular pathophysiology related to the COVID-19 phenotypes, collecting the mechanistic knowledge scattered across scientific literature and bioinformatic databases, and integrating it using a logical/conceptual model systems biology. The model framework building process began from the results of a domain-based literature review regarding a multiomics approach to COVID-19. This evidence allowed us to define a framework of COVID-19 conceptual model and to report all concepts in a multilevel and multiscale structure. The same interdisciplinary working groups that carried out the scoping review were involved. The conclusive result is a conceptual method to design multiscale models of infectious diseases. The methodology, applied in this paper, is a set of partially ordered research and development activities that result in a COVID-19 multiscale model.

Integrating 3D structural information into systems biology

Systems biology is a data-heavy field that focuses on systems-wide depictions of biological phenomena necessarily sacrificing a detailed characterization of individual components. As an example, genome-wide protein interaction networks are widely used in systems biology and continuously extended and refined as new sources of evidence become available. Despite the vast amount of information about individual protein structures and protein complexes that has accumulated in the past 50 years in the Protein Data Bank, the data, computational tools, and language of structural biology are not an integral part of systems biology. However, increasing effort has been devoted to this integration, and the related literature is reviewed here. Relationships between proteins that are detected via structural similarity offer a rich source of information not available from sequence similarity, and homology modeling can be used to leverage Protein Data Bank structures to produce 3D models for a significant fraction of many proteomes. A number of structure-informed genomic and cross-species (i.e., virus–host) interactomes will be described, and the unique information they provide will be illustrated with a number of examples. Tissue- and tumor-specific interactomes have also been developed through computational strategies that exploit patient information and through genetic interactions available from increasingly sensitive screens. Strategies to integrate structural information with these alternate data sources will be described. Finally, efforts to link protein structure space with chemical compound space offer novel sources of information in drug design, off-target identification, and the identification of targets for compounds found to be effective in phenotypic screens.

My personal mutanome: a computational genomic medicine platform for searching network perturbing alleles linking genotype to phenotype

Massive genome sequencing data have inspired new challenges in personalized treatments and facilitated oncological drug discovery. We present a comprehensive database, My Personal Mutanome (MPM), for accelerating the development of precision cancer medicine protocols. MPM contains 490,245 mutations from over 10,800 tumor exomes across 33 cancer types in The Cancer Genome Atlas mapped to 94,563 structure-resolved/predicted protein-protein interaction interfaces ("edgetic") and 311,022 functional sites ("nodetic"), including ligand-protein binding sites and 8 types of protein posttranslational modifications. In total, 8884 survival results and 1,271,132 drug responses are obtained for these mapped interactions. MPM is available at https://mutanome.lerner.ccf.org .

Decoding Protein-protein Interactions: An Overview

Drug discovery has focused on the paradigm "one drug, one target" for a long time. However, small molecules can act at multiple macromolecular targets, which serves as the basis for drug repurposing. In an effort to expand the target space, and given advances in X-ray crystallography, protein-protein interactions have become an emerging focus area of drug discovery enterprises. Proteins interact with other biomolecules and it is this intricate network of interactions that determines the behavior of the system and its biological processes. In this review, we briefly discuss networks in disease, followed by computational methods for protein-protein complex prediction. Computational methodologies and techniques employed towards objectives such as protein-protein docking, protein-protein interactions, and interface predictions are described extensively. Docking aims at producing a complex between proteins, while interface predictions identify a subset of residues on one protein that could interact with a partner, and protein-protein interaction sites address whether two proteins interact. In addition, approaches to predict hot spots and binding sites are presented along with a representative example of our internal project on the chemokine CXC receptor 3 B-isoform and predictive modeling with IP10 and PF4.

Dynamic 3D proteomes reveal protein functional alterations at high resolution in situ

Biological processes are regulated by intermolecular interactions and chemical modifications that do not affect protein levels, thus escaping detection in classical proteomic screens. We demonstrate here that a global protein structural readout based on limited proteolysis-mass spectrometry (LiP-MS) detects many such functional alterations, simultaneously and in situ, in bacteria undergoing nutrient adaptation and in yeast responding to acute stress. The structural readout, visualized as structural barcodes, captured enzyme activity changes, phosphorylation, protein aggregation, and complex formation, with the resolution of individual regulated functional sites such as binding and active sites. Comparison with prior knowledge, including other ‘omics data, showed that LiP-MS detects many known functional alterations within well-studied pathways. It suggested distinct metabolite-protein interactions and enabled identification of a fructose-1,6-bisphosphate-based regulatory mechanism of glucose uptake in E. coli. The structural readout dramatically increases classical proteomics coverage, generates mechanistic hypotheses, and paves the way for in situ structural systems biology.

•

Dynamic structural proteomic screens detect functional changes at high resolution

•

Detect enzyme activity, phosphorylation, and molecular interactions in situ

•

Generate new molecular hypotheses and increase functional proteomics coverage

•

Enabled discovery of a regulatory mechanism of glucose uptake in E. coli

Single-cell network biology for resolving cellular heterogeneity in human diseases

Understanding cellular heterogeneity is the holy grail of biology and medicine. Cells harboring identical genomes show a wide variety of behaviors in multicellular organisms. Genetic circuits underlying cell-type identities will facilitate the understanding of the regulatory programs for differentiation and maintenance of distinct cellular states. Such a cell-type-specific gene network can be inferred from coregulatory patterns across individual cells. Conventional methods of transcriptome profiling using tissue samples provide only average signals of diverse cell types. Therefore, reconstructing gene regulatory networks for a particular cell type is not feasible with tissue-based transcriptome data. Recently, single-cell omics technology has emerged and enabled the capture of the transcriptomic landscape of every individual cell. Although single-cell gene expression studies have already opened up new avenues, network biology using single-cell transcriptome data will further accelerate our understanding of cellular heterogeneity. In this review, we provide an overview of single-cell network biology and summarize recent progress in method development for network inference from single-cell RNA sequencing (scRNA-seq) data. Then, we describe how cell-type-specific gene networks can be utilized to study regulatory programs specific to disease-associated cell types and cellular states. Moreover, with scRNA data, modeling personal or patient-specific gene networks is feasible. Therefore, we also introduce potential applications of single-cell network biology for precision medicine. We envision a rapid paradigm shift toward single-cell network analysis for systems biology in the near future.

Gene regulatory networks reconstructed from single-cell RNA sequencing datasets are allowing researchers to better understand the molecular circuits and cell states that contribute to complex human disease. Junha Cha and Insuk Lee from Yonsei University in Seoul, South Korea, review the concept of ‘single-cell network biology’, which involves using computational algorithms on genetic expression data from thousands of cells to infer functional interactions in various biological contexts. This systems biology approach to analyzing the profiles of messenger RNA in single cells is helping researchers discover new signaling pathways that could serve as disease biomarkers or therapeutic targets. In the future, patient-specific models of personal gene networks could explain why certain genetic variants affect disease risk. This research could also eventually lead to new types of individualized medical treatments.

Inferring the molecular and phenotypic impact of amino acid variants with MutPred2

Identifying pathogenic variants and underlying functional alterations is challenging. To this end, we introduce MutPred2, a tool that improves the prioritization of pathogenic amino acid substitutions over existing methods, generates molecular mechanisms potentially causative of disease, and returns interpretable pathogenicity score distributions on individual genomes. Whilst its prioritization performance is state-of-the-art, a distinguishing feature of MutPred2 is the probabilistic modeling of variant impact on specific aspects of protein structure and function that can serve to guide experimental studies of phenotype-altering variants. We demonstrate the utility of MutPred2 in the identification of the structural and functional mutational signatures relevant to Mendelian disorders and the prioritization of de novo mutations associated with complex neurodevelopmental disorders. We then experimentally validate the functional impact of several variants identified in patients with such disorders. We argue that mechanism-driven studies of human inherited disease have the potential to significantly accelerate the discovery of clinically actionable variants.

Amino acid homorepeats in proteins

Amino acid homorepeats, or homorepeats, are polypeptide segments found in proteins that contain stretches of identical amino acid residues. Although abnormal homorepeat expansions are linked to pathologies such as neurodegenerative diseases, homorepeats are prevalent in eukaryotic proteomes, suggesting that they are important for normal physiology. In this Review, we discuss recent advances in our understanding of the biological functions of homorepeats, which range from facilitating subcellular protein localization to mediating interactions between proteins across diverse cellular pathways. We explore how the functional diversity of homorepeat-containing proteins could be linked to the ability of homorepeats to adopt different structural conformations, an ability influenced by repeat composition, repeat length and the nature of flanking sequences. We conclude by highlighting how an understanding of homorepeats will help us better characterize and develop therapeutics against the human diseases to which they contribute.

Guilt-by-Association - Functional Insights Gained From Studying the LRRK2 Interactome

The Parkinson's disease-associated Leucine-rich repeat kinase 2 (LRRK2) is a complex multi-domain protein belonging to the Roco protein family, a unique group of G-proteins. Variants of this gene are associated with an increased risk of Parkinson's disease. Besides its well-characterized enzymatic activities, conferred by its GTPase and kinase domains, and a central dimerization domain, it contains four predicted repeat domains, which are, based on their structure, commonly involved in protein-protein interactions (PPIs). In the past decades, tremendous progress has been made in determining comprehensive interactome maps for the human proteome. Knowledge of PPIs has been instrumental in assigning functions to proteins involved in human disease and helped to understand the connectivity between different disease pathways and also significantly contributed to the functional understanding of LRRK2. In addition to an increased kinase activity observed for proteins containing PD-associated variants, various studies helped to establish LRRK2 as a large scaffold protein in the interface between cytoskeletal dynamics and the vesicular transport. This review first discusses a number of specific LRRK2-associated PPIs for which a functional consequence can at least be speculated upon, and then considers the representation of LRRK2 protein interactions in public repositories, providing an outlook on open research questions and challenges in this field.

Computer-aided biomarker discovery for precision medicine: data resources, models and applications

Biomarkers are a class of measurable and evaluable indicators with the potential to predict disease initiation and progression. In contrast to disease-associated factors, biomarkers hold the promise to capture the changeable signatures of biological states. With methodological advances, computer-aided biomarker discovery has now become a burgeoning paradigm in the field of biomedical science. In recent years, the 'big data' term has accumulated for the systematical investigation of complex biological phenomena and promoted the flourishing of computational methods for systems-level biomarker screening. Compared with routine wet-lab experiments, bioinformatics approaches are more efficient to decode disease pathogenesis under a holistic framework, which is propitious to identify biomarkers ranging from single molecules to molecular networks for disease diagnosis, prognosis and therapy. In this review, the concept and characteristics of typical biomarker types, e.g. single molecular biomarkers, module/network biomarkers, cross-level biomarkers, etc., are explicated on the guidance of systems biology. Then, publicly available data resources together with some well-constructed biomarker databases and knowledge bases are introduced. Biomarker identification models using mathematical, network and machine learning theories are sequentially discussed. Based on network substructural and functional evidences, a novel bioinformatics model is particularly highlighted for microRNA biomarker discovery. This article aims to give deep insights into the advantages and challenges of current computational approaches for biomarker detection, and to light up the future wisdom toward precision medicine and nation-wide healthcare.

RNA Polymerase II CTD phosphatase Rtr1 fine-tunes transcription termination

RNA Polymerase II (RNAPII) transcription termination is regulated by the phosphorylation status of the C-terminal domain (CTD). The phosphatase Rtr1 has been shown to regulate serine 5 phosphorylation on the CTD; however, its role in the regulation of RNAPII termination has not been explored. As a consequence of RTR1 deletion, interactions within the termination machinery and between the termination machinery and RNAPII were altered as quantified by Disruption-Compensation (DisCo) network analysis. Of note, interactions between RNAPII and the cleavage factor IA (CF1A) subunit Pcf11 were reduced in rtr1Δ, whereas interactions with the CTD and RNA-binding termination factor Nrd1 were increased. Globally, rtr1Δ leads to decreases in numerous noncoding RNAs that are linked to the Nrd1, Nab3 and Sen1 (NNS) -dependent RNAPII termination pathway. Genome-wide analysis of RNAPII and Nrd1 occupancy suggests that loss of RTR1 leads to increased termination at noncoding genes. Additionally, premature RNAPII termination increases globally at protein-coding genes with a decrease in RNAPII occupancy occurring just after the peak of Nrd1 recruitment during early elongation. The effects of rtr1Δ on RNA expression levels were lost following deletion of the exosome subunit Rrp6, which works with the NNS complex to rapidly degrade a number of noncoding RNAs following termination. Overall, these data suggest that Rtr1 restricts the NNS-dependent termination pathway in WT cells to prevent premature termination of mRNAs and ncRNAs. Rtr1 facilitates low-level elongation of noncoding transcripts that impact RNAPII interference thereby shaping the transcriptome.

Many cellular RNAs including those that encode for proteins are produced by the enzyme RNA Polymerase II. In this work, we have defined a new role for the phosphatase Rtr1 in the regulation of RNA Polymerase II progression from the start of transcription to the 3’ end of the gene where the nascent RNA from protein-coding genes is typically cleaved and polyadenylated. Deletion of the gene that encodes RTR1 leads to changes in the interactions between RNA polymerase II and the termination machinery. Rtr1 loss also causes early termination of RNA Polymerase II at many of its target gene types, including protein coding genes and noncoding RNAs. Evidence suggests that the premature termination observed in RTR1 knockout cells occurs through the termination factor and RNA binding protein Nrd1 and its binding partner Nab3. Deletion of RRP6, a known component of the Nrd1-Nab3 termination coupled RNA degradation pathway, is epistatic to RTR1 suggesting that Rrp6 is required to terminate and/or degrade many of the noncoding RNAs that have increased turnover in RTR1 deletion cells. These findings suggest that Rtr1 normally promotes elongation of RNA Polymerase II transcripts through prevention of Nrd1-directed termination.

Domain-mediated interactions for protein subfamily identification

Within a protein family, proteins with the same domain often exhibit different cellular functions, despite the shared evolutionary history and molecular function of the domain. We hypothesized that domain-mediated interactions (DMIs) may categorize a protein family into subfamilies because the diversified functions of a single domain often depend on interacting partners of domains. Here we systematically identified DMI subfamilies, in which proteins share domains with DMI partners, as well as with various functional and physical interaction networks in individual species. In humans, DMI subfamily members are associated with similar diseases, including cancers, and are frequently co-associated with the same diseases. DMI information relates to the functional and evolutionary subdivisions of human kinases. In yeast, DMI subfamilies contain proteins with similar phenotypic outcomes from specific chemical treatments. Therefore, the systematic investigation here provides insights into the diverse functions of subfamilies derived from a protein family with a link-centric approach and suggests a useful resource for annotating the functions and phenotypic outcomes of proteins.

Genetics' Piece of the PI: Inferring the Origin of Complex Traits and Diseases from Proteome-Wide Protein-Protein Interaction Dynamics

How do common and rare genetic polymorphisms contribute to quantitative traits or disease risk and progression? Multiple human traits have been extensively characterized at the genomic level, revealing their complex genetic architecture. However, it is difficult to resolve the mechanisms by which specific variants contribute to a phenotype. Recently, analyses of variant effects on molecular traits have uncovered intermediate mechanisms that link sequence variation to phenotypic changes. Yet, these methods only capture a fraction of genetic contributions to phenotype. Here, in reviewing the field, it is proposed that complex traits can be understood by characterizing the dynamics of biochemical networks within living cells, and that the effects of genetic variation can be captured on these networks by using protein-protein interaction (PPI) methodologies. This synergy between PPI methodologies and the genetics of complex traits opens new avenues to investigate the molecular etiology of human diseases and to facilitate their prevention or treatment.

Enriching Human Interactome with Functional Mutations to Detect High-Impact Network Modules Underlying Complex Diseases

Rapid progress in high-throughput -omics technologies moves us one step closer to the datacalypse in life sciences. In spite of the already generated volumes of data, our knowledge of the molecular mechanisms underlying complex genetic diseases remains limited. Increasing evidence shows that biological networks are essential, albeit not sufficient, for the better understanding of these mechanisms. The identification of disease-specific functional modules in the human interactome can provide a more focused insight into the mechanistic nature of the disease. However, carving a disease network module from the whole interactome is a difficult task. In this paper, we propose a computational framework, Discovering most IMpacted SUbnetworks in interactoMe (DIMSUM), which enables the integration of genome-wide association studies (GWAS) and functional effects of mutations into the protein-protein interaction (PPI) network to improve disease module detection. Specifically, our approach incorporates and propagates the functional impact of non-synonymous single nucleotide polymorphisms (nsSNPs) on PPIs to implicate the genes that are most likely influenced by the disruptive mutations, and to identify the module with the greatest functional impact. Comparison against state-of-the-art seed-based module detection methods shows that our approach could yield modules that are biologically more relevant and have stronger association with the studied disease. We expect for our method to become a part of the common toolbox for the disease module analysis, facilitating the discovery of new disease markers.

bZIP edgetic mutations: at the frontier of plant metabolism, development and stress trade-off

This article comments on: Garg A, Kirchler T, Fillinger S, Wanke F, Stadelhofer B, Stahl M, Chaban C. 2019. Targeted manipulation of bZIP53 DNA-binding properties influences Arabidopsis metabolism and growth. Journal of Experimental Botany 70, 5659–5671.

The pathobiology of perturbed mutant huntingtin protein-protein interactions in Huntington's disease

Mutations are at the root of many human diseases. Still, we largely do not exactly understand how they trigger pathogenesis. One, more recent, hypothesis has been that they comprehensively perturb protein-protein interaction (PPI) networks and significantly alter key biological processes. Under this premise, many rare genetic disorders with Mendelian inheritance, like Huntington's disease and several spinocerebellar ataxias, are likely to be caused by complex genotype-phenotype relationships involving abnormal PPIs. These altered PPI networks and their effects on cellular pathways are poorly understood at the molecular level. In this review, we focus on PPIs that are perturbed by the expanded pathogenic polyglutamine tract in huntingtin (HTT), the protein which, in its mutated form, leads to the autosomal dominant, neurodegenerative Huntington's disease. One aspect of perturbed mutant HTT interactions is the formation of abnormal protein species such as fibrils or large neuronal inclusions as a result of homotypic and heterotypic aberrant molecular interactions. This review focuses on abnormal PPIs that are associated with the assembly of mutant HTT aggregates in cells and their potential relevance in disease. Furthermore, the mechanisms and pathobiological processes that may contribute to phenotype development, neuronal dysfunction and toxicity in Huntington's disease brains are also discussed. This article is part of the Special Issue "Proteomics".

Systematic domain-based aggregation of protein structures highlights DNA-, RNA- and other ligand-binding positions

Domains are fundamental subunits of proteins, and while they play major roles in facilitating protein–DNA, protein–RNA and other protein–ligand interactions, a systematic assessment of their various interaction modes is still lacking. A comprehensive resource identifying positions within domains that tend to interact with nucleic acids, small molecules and other ligands would expand our knowledge of domain functionality as well as aid in detecting ligand-binding sites within structurally uncharacterized proteins. Here, we introduce an approach to identify per-domain-position interaction ‘frequencies’ by aggregating protein co-complex structures by domain and ascertaining how often residues mapping to each domain position interact with ligands. We perform this domain-based analysis on ∼91000 co-complex structures, and infer positions involved in binding DNA, RNA, peptides, ions or small molecules across 4128 domains, which we refer to collectively as the InteracDome. Cross-validation testing reveals that ligand-binding positions for 2152 domains are highly consistent and can be used to identify residues facilitating interactions in ∼63–69% of human genes. Our resource of domain-inferred ligand-binding sites should be a great aid in understanding disease etiology: whereas these sites are enriched in Mendelian-associated and cancer somatic mutations, they are depleted in polymorphisms observed across healthy populations. The InteracDome is available at http://interacdome.princeton.edu.

Estimating dispensable content in the human interactome

Protein-protein interaction (PPI) networks (interactome networks) have successfully advanced our knowledge of molecular function, disease and evolution. While much progress has been made in quantifying errors and biases in experimental PPI datasets, it remains unknown what fraction of the error-free PPIs in the cell are completely dispensable, i.e., effectively neutral upon disruption. Here, we estimate dispensable content in the human interactome by calculating the fractions of PPIs disrupted by neutral and non-neutral mutations. Starting with the human reference interactome determined by experiments, we construct a human structural interactome by building homology-based three-dimensional structural models for PPIs. Next, we map common mutations from healthy individuals as well as Mendelian disease-causing mutations onto the human structural interactome, and perform structure-based calculations of how these mutations perturb the interactome. Using our predicted as well as experimentally-determined interactome perturbation patterns by common and disease mutations, we estimate that <~20% of the human interactome is completely dispensable.

The fraction of protein-protein interactions (PPIs) that can be disrupted without fitness effect is unknown. Here, the authors model how disease-causing mutations and common mutations carried by healthy people perturb the interactome, and estimate that <20% of human PPIs are completely dispensable.

Applications of in Silico Methods for Design and Development of Drugs Targeting Protein-Protein Interactions

BACKGROUND

Identification of Protein-Protein Interactions (PPIs) is a major challenge in modern molecular biology and biochemistry research, due to the unquestionable role of proteins in cells, biological process and pathological states. Over the past decade, the PPIs have evolved from being considered a highly challenging field of research to being investigated and examined as targets for pharmacological intervention.

OBJECTIVE

Comprehension of protein interactions is crucial to known how proteins come together to build signalling pathways, to carry out their functions, or to cause diseases, when deregulated. Multiplicity and great amount of PPIs structures offer a huge number of new and potential targets for the treatment of different diseases.

METHODS

Computational techniques are becoming predominant in PPIs studies for their effectiveness, flexibility, accuracy and cost. As a matter of fact, there are effective in silico approaches which are able to identify PPIs and PPI site. Such methods for computational target prediction have been developed through molecular descriptors and data-mining procedures.

RESULTS

In this review, we present different types of interactions between protein-protein and the application of in silico methods for design and development of drugs targeting PPIs. We described computational approaches for the identification of possible targets on protein surface and to detect of stimulator/ inhibitor molecules.

CONCLUSION

A deeper study of the most recent bioinformatics methodologies for PPIs studies is vital for a better understanding of protein complexes and for discover new potential PPI modulators in therapeutic intervention.

Protein ensembles link genotype to phenotype

Classically, phenotype is what is observed, and genotype is the genetic makeup. Statistical studies aim to project phenotypic likelihoods of genotypic patterns. The traditional genotype-to-phenotype theory embraces the view that the encoded protein shape together with gene expression level largely determines the resulting phenotypic trait. Here, we point out that the molecular biology revolution at the turn of the century explained that the gene encodes not one but ensembles of conformations, which in turn spell all possible gene-associated phenotypes. The significance of a dynamic ensemble view is in understanding the linkage between genetic change and the gained observable physical or biochemical characteristics. Thus, despite the transformative shift in our understanding of the basis of protein structure and function, the literature still commonly relates to the classical genotype-phenotype paradigm. This is important because an ensemble view clarifies how even seemingly small genetic alterations can lead to pleiotropic traits in adaptive evolution and in disease, why cellular pathways can be modified in monogenic and polygenic traits, and how the environment may tweak protein function.

Pathway-specific protein domains are predictive for human diseases

Protein domains are basic functional units of proteins. Many protein domains are pervasive among diverse biological processes, yet some are associated with specific pathways. Human complex diseases are generally viewed as pathway-level disorders. Therefore, we hypothesized that pathway-specific domains could be highly informative for human diseases. To test the hypothesis, we developed a network-based scoring scheme to quantify specificity of domain-pathway associations. We first generated domain profiles for human proteins, then constructed a co-pathway protein network based on the associations between domain profiles. Based on the score, we classified human protein domains into pathway-specific domains (PSDs) and non-specific domains (NSDs). We found that PSDs contained more pathogenic variants than NSDs. PSDs were also enriched for disease-associated mutations that disrupt protein-protein interactions (PPIs) and tend to have a moderate number of domain interactions. These results suggest that mutations in PSDs are likely to disrupt within-pathway PPIs, resulting in functional failure of pathways. Finally, we demonstrated the prediction capacity of PSDs for disease-associated genes with experimental validations in zebrafish. Taken together, the network-based quantitative method of modeling domain-pathway associations presented herein suggested underlying mechanisms of how protein domains associated with specific pathways influence mutational impacts on diseases via perturbations in within-pathway PPIs, and provided a novel genomic feature for interpreting genetic variants to facilitate the discovery of human disease genes.

Protein domains are basic functional units of proteins, yet domain-based pathway annotations for proteins are challenging tasks because many domains are pervasive among diverse pathways. Therefore, we developed a network-based scoring scheme to measure pathway specificity of domains, and then used it to identify pathway-specific domains. Surprisingly, we observed substantially more disease mutations in pathway-specific domains than non-specific domains. We found evidences that mutations of pathway-specific domains tend to perturb pathway integrity via disrupting within-pathway protein-protein interactions. We also demonstrated prediction capacity of pathway-specific domains for complex diseases with experimental validations. Our study demonstrated the usefulness of pathway information for protein domains in interpreting non-random distribution of disease mutations among domains and identification of disease genes and variants.

RMut: R package for a Boolean sensitivity analysis against various types of mutations

There have been many in silico studies based on a Boolean network model to investigate network sensitivity against gene or interaction mutations. However, there are no proper tools to examine the network sensitivity against many different types of mutations, including user-defined ones. To address this issue, we developed RMut, which is an R package to analyze the Boolean network-based sensitivity by efficiently employing not only many well-known node-based and edgetic mutations but also novel user-defined mutations. In addition, RMut can specify the mutation area and the duration time for more precise analysis. RMut can be used to analyze large-scale networks because it is implemented in a parallel algorithm using the OpenCL library. In the first case study, we observed that the real biological networks were most sensitive to overexpression/state-flip and edge-addition/-reverse mutations among node-based and edgetic mutations, respectively. In the second case study, we showed that edgetic mutations can predict drug-targets better than node-based mutations. Finally, we examined the network sensitivity to double edge-removal mutations and found an interesting synergistic effect. Taken together, these findings indicate that RMut is a flexible R package to efficiently analyze network sensitivity to various types of mutations. RMut is available at https://github.com/csclab/RMut.

A network-centric approach to drugging TNF-induced NF-κB signaling

Target-centric drug development strategies prioritize single-target potency in vitro and do not account for connectivity and multi-target effects within a signal transduction network. Here, we present a systems biology approach that combines transcriptomic and structural analyses with live-cell imaging to predict small molecule inhibitors of TNF-induced NF-κB signaling and elucidate the network response. We identify two first-in-class small molecules that inhibit the NF-κB signaling pathway by preventing the maturation of a rate-limiting multiprotein complex necessary for IKK activation. Our findings suggest that a network-centric drug discovery approach is a promising strategy to evaluate the impact of pharmacologic intervention in signaling.

Chemical perturbation of specific protein–protein interactions is notoriously difficult, yet necessary when complete inhibition of a signalling pathway is detrimental to the cell. Here, the authors use a systems approach and identify two first-in-class small molecules that specifically inhibit TNF-induced NF-κB activation.