Edgetic perturbation models of human inherited disorders

Ciliopathy-associated missense mutations in IFT140 are tolerated by the inherent resilience of the IFT machinery

Genotype-phenotype correlations of rare diseases are complicated by low patient number, high phenotype variability and compound heterozygosity. Mutations may cause instability of single proteins, and affect protein complex formation or overall robustness of a specific process in a given cell. Ciliopathies offer an interesting case for studying genotype-phenotype correlations as they have a spectrum of severity and include diverse phenotypes depending on different mutations in the same protein. For instance, mutations in the intraflagellar transport protein IFT140 cause a vast spectrum of ciliopathies ranging from isolated retinal dystrophy to severe skeletal abnormalities and multi-organ diseases such as Mainzer-Saldino and Jeune syndrome. Here, the quantitative effects of 23 missense mutations in IFT140, which forms part of the crucial IFT-A complex of the ciliary machinery, were analyzed using affinity purification coupled to mass spectrometry (AP-MS). A subset of 10 mutations led to a significant and domain-specific reduction in IFT140-IFT-A complex interaction indicating complex formation issues and potentially hampering its molecular function. Knockout of IFT140 led to loss of cilia, as shown before. However, phenotypically only mild effects concerning cilia assembly were observed for two out of four tested IFT140 missense mutations. Therefore, our results demonstrate the utility of AP-MS in discerning pathogenic MMs from polymorphisms, and we postulate that reduced function is tolerated by the evolutionarily highly conserved IFT-A system.

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.

The Molecular Landscape of Premature Aging Diseases Defined by Multilayer Network Exploration

Premature Aging (PA) diseases are rare genetic disorders that mimic some aspects of physiological aging at an early age. Various causative genes of PA diseases have been identified in recent years, providing insights into some dysfunctional cellular processes. However, the identification of PA genes also revealed significant genetic heterogeneity and highlighted the gaps in this understanding of PA-associated molecular mechanisms. Furthermore, many patients remain undiagnosed. Overall, the current lack of knowledge about PA diseases hinders the development of effective diagnosis and therapies and poses significant challenges to improving patient care. Here, a network-based approach to systematically unravel the cellular functions disrupted in PA diseases is presented. Leveraging a network community identification algorithm, it is delved into a vast multilayer network of biological interactions to extract the communities of 67 PA diseases from their 132 associated genes. It is found that these communities can be grouped into six distinct clusters, each reflecting specific cellular functions: DNA repair, cell cycle, transcription regulation, inflammation, cell communication, and vesicle-mediated transport. That these clusters collectively represent the landscape of the molecular mechanisms that are perturbed in PA diseases, providing a framework for better understanding their pathogenesis is proposed. Intriguingly, most clusters also exhibited a significant enrichment in genes associated with physiological aging, suggesting a potential overlap between the molecular underpinnings of PA diseases and natural aging.

A perturbation approach for refining Boolean models of cell cycle regulation

Considerable effort is required to build mathematical models of large protein regulatory networks. Utilizing computational algorithms that guide model development can significantly streamline the process and enhance the reliability of the resulting models. In this article, we present a perturbation approach for developing data-centric Boolean models of cell cycle regulation. To evaluate networks, we assign a score based on their steady states and the dynamical trajectories corresponding to the initial conditions. Then, perturbation analysis is used to find new networks with lower scores, in which dynamical trajectories traverse through the correct cell cycle path with high frequency. We apply this method to refine Boolean models of cell cycle regulation in budding yeast and mammalian cells.

Structure-guided Evolutionary Analysis of Interactome Network Rewiring at Single Residue Resolution in Yeasts

Protein-protein interactions (PPIs) are known to rewire extensively during evolution leading to lineage-specific and species-specific changes in molecular processes. However, the detailed molecular evolutionary mechanisms underlying interactome network rewiring are not well-understood. Here, we combine high-confidence PPI data, high-resolution three-dimensional structures of protein complexes, and homology-based structural annotation transfer to construct structurally-resolved interactome networks for the two yeasts S. cerevisiae and S. pombe. We then classify PPIs according to whether they are preserved or different between the two yeast species and compare site-specific evolutionary rates of interfacial versus non-interfacial residues for these different categories of PPIs. We find that residues in PPI interfaces evolve significantly more slowly than non-interfacial residues when using lineage-specific measures of evolutionary rate, but not when using non-lineage-specific measures. Furthermore, both lineage-specific and non-lineage-specific evolutionary rate measures can distinguish interfacial residues from non-interfacial residues for preserved PPIs between the two yeasts, but only the lineage-specific measure is appropriate for rewired PPIs. Finally, both lineage-specific and non-lineage-specific evolutionary rate measures are appropriate for elucidating structural determinants of protein evolution for residues outside of PPI interfaces. Overall, our results demonstrate that unlike tertiary structures of single proteins, PPIs and PPI interfaces can be highly volatile in their evolution, thus requiring the use of lineage-specific measures when studying their evolution. These results yield insight into the evolutionary design principles of PPIs and the mechanisms by which interactions are preserved or rewired between species, improving our understanding of the molecular evolution of PPIs and PPI interfaces at the residue level.

Understanding the interplay between dNTP metabolism and genome stability in cancer

The size and composition of the intracellular DNA precursor pool is integral to the maintenance of genome stability, and this relationship is fundamental to our understanding of cancer. Key aspects of carcinogenesis, including elevated mutation rates and induction of certain types of DNA damage in cancer cells, can be linked to disturbances in deoxynucleoside triphosphate (dNTP) pools. Furthermore, our approaches to treat cancer heavily exploit the metabolic interplay between the DNA and the dNTP pool, with a long-standing example being the use of antimetabolite-based cancer therapies, and this strategy continues to show promise with the development of new targeted therapies. In this Review, we compile the current knowledge on both the causes and consequences of dNTP pool perturbations in cancer cells, together with their impact on genome stability. We outline several outstanding questions remaining in the field, such as the role of dNTP catabolism in genome stability and the consequences of dNTP pool expansion. Importantly, we detail how our mechanistic understanding of these processes can be utilised with the aim of providing better informed treatment options to patients with cancer.

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 energetic and allosteric landscape for KRAS inhibition

Thousands of proteins have been validated genetically as therapeutic targets for human diseases1. However, very few have been successfully targeted, and many are considered 'undruggable'. This is particularly true for proteins that function via protein-protein interactions-direct inhibition of binding interfaces is difficult and requires the identification of allosteric sites. However, most proteins have no known allosteric sites, and a comprehensive allosteric map does not exist for any protein. Here we address this shortcoming by charting multiple global atlases of inhibitory allosteric communication in KRAS. We quantified the effects of more than 26,000 mutations on the folding of KRAS and its binding to six interaction partners. Genetic interactions in double mutants enabled us to perform biophysical measurements at scale, inferring more than 22,000 causal free energy changes. These energy landscapes quantify how mutations tune the binding specificity of a signalling protein and map the inhibitory allosteric sites for an important therapeutic target. Allosteric propagation is particularly effective across the central β-sheet of KRAS, and multiple surface pockets are genetically validated as allosterically active, including a distal pocket in the C-terminal lobe of the protein. Allosteric mutations typically inhibit binding to all tested effectors, but they can also change the binding specificity, revealing the regulatory, evolutionary and therapeutic potential to tune pathway activation. Using the approach described here, it should be possible to rapidly and comprehensively identify allosteric target sites in many proteins.

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.

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.

Structure-based prediction of Ras-effector binding affinities and design of "branchegetic" interface mutations

Ras is a central cellular hub protein controlling multiple cell fates. How Ras interacts with a variety of potential effector proteins is relatively unexplored, with only some key effectors characterized in great detail. Here, we have used homology modeling based on X-ray and AlphaFold2 templates to build structural models for 54 Ras-effector complexes. These models were used to estimate binding affinities using a supervised learning regressor. Furthermore, we systematically introduced Ras "branch-pruning" (or branchegetic) mutations to identify 200 interface mutations that affect the binding energy with at least one of the model structures. The impacts of these branchegetic mutants were integrated into a mathematical model to assess the potential for rewiring interactions at the Ras hub on a systems level. These findings have provided a quantitative understanding of Ras-effector interfaces and their impact on systems properties of a key cellular hub.

Toward the analysis of functional proteoforms using mass spectrometry-based stability proteomics

Functional proteomics aims to elucidate biological functions, mechanisms, and pathways of proteins and proteoforms at the molecular level to examine complex cellular systems and disease states. A series of stability proteomics methods have been developed to examine protein functionality by measuring the resistance of a protein to chemical or thermal denaturation or proteolysis. These methods can be applied to measure the thermal stability of thousands of proteins in complex biological samples such as cell lysate, intact cells, tissues, and other biological fluids to measure proteome stability. Stability proteomics methods have been popularly applied to observe stability shifts upon ligand binding for drug target identification. More recently, these methods have been applied to characterize the effect of structural changes in proteins such as those caused by post-translational modifications (PTMs) and mutations, which can affect protein structures or interactions and diversify protein functions. Here, we discussed the current application of a suite of stability proteomics methods, including thermal proteome profiling (TPP), stability of proteomics from rates of oxidation (SPROX), and limited proteolysis (LiP) methods, to observe PTM-induced structural changes on protein stability. We also discuss future perspectives highlighting the integration of top-down mass spectrometry and stability proteomics methods to characterize intact proteoform stability and understand the function of variable protein modifications.

Rewiring Drug Research and Development through Human Data-Driven Discovery (HD3)

In an era of unparalleled technical advancement, the pharmaceutical industry is struggling to transform data into increased research and development efficiency, and, as a corollary, new drugs for patients. Here, we briefly review some of the commonly discussed issues around this counterintuitive innovation crisis. Looking at both industry- and science-related factors, we posit that traditional preclinical research is front-loading the development pipeline with data and drug candidates that are unlikely to succeed in patients. Applying a first principles analysis, we highlight the critical culprits and provide suggestions as to how these issues can be rectified through the pursuit of a Human Data-driven Discovery (HD3) paradigm. Consistent with other examples of disruptive innovation, we propose that new levels of success are not dependent on new inventions, but rather on the strategic integration of existing data and technology assets. In support of these suggestions, we highlight the power of HD3, through recently published proof-of-concept applications in the areas of drug safety analysis and prediction, drug repositioning, the rational design of combination therapies and the global response to the COVID-19 pandemic. We conclude that innovators must play a key role in expediting the path to a largely human-focused, systems-based approach to drug discovery and research.

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.

A systematic review of computational approaches to understand cancer biology for informed drug repurposing

Cancer is the second leading cause of death globally, trailing only heart disease. In the United States alone, 1.9 million new cancer cases and 609,360 deaths were recorded for 2022. Unfortunately, the success rate for new cancer drug development remains less than 10%, making the disease particularly challenging. This low success rate is largely attributed to the complex and poorly understood nature of cancer etiology. Therefore, it is critical to find alternative approaches to understanding cancer biology and developing effective treatments. One such approach is drug repurposing, which offers a shorter drug development timeline and lower costs while increasing the likelihood of success. In this review, we provide a comprehensive analysis of computational approaches for understanding cancer biology, including systems biology, multi-omics, and pathway analysis. Additionally, we examine the use of these methods for drug repurposing in cancer, including the databases and tools that are used for cancer research. Finally, we present case studies of drug repurposing, discussing their limitations and offering recommendations for future research in this area.

Formate overflow drives toxic folate trapping in MTHFD1 inhibited cancer cells

Cancer cells fuel their increased need for nucleotide supply by upregulating one-carbon (1C) metabolism, including the enzymes methylenetetrahydrofolate dehydrogenase-cyclohydrolase 1 and 2 (MTHFD1 and MTHFD2). TH9619 is a potent inhibitor of dehydrogenase and cyclohydrolase activities in both MTHFD1 and MTHFD2, and selectively kills cancer cells. Here, we reveal that, in cells, TH9619 targets nuclear MTHFD2 but does not inhibit mitochondrial MTHFD2. Hence, overflow of formate from mitochondria continues in the presence of TH9619. TH9619 inhibits the activity of MTHFD1 occurring downstream of mitochondrial formate release, leading to the accumulation of 10-formyl-tetrahydrofolate, which we term a 'folate trap'. This results in thymidylate depletion and death of MTHFD2-expressing cancer cells. This previously uncharacterized folate trapping mechanism is exacerbated by physiological hypoxanthine levels that block the de novo purine synthesis pathway, and additionally prevent 10-formyl-tetrahydrofolate consumption for purine synthesis. The folate trapping mechanism described here for TH9619 differs from other MTHFD1/2 inhibitors and antifolates. Thus, our findings uncover an approach to attack cancer and reveal a regulatory mechanism in 1C metabolism.

Genetics of mitochondrial diseases: Current approaches for the molecular diagnosis

Mitochondrial diseases are a genetically and phenotypically variable set of monogenic disorders. The main characteristic of mitochondrial diseases is a defective oxidative phosphorylation. Both nuclear and mitochondrial DNA encode the approximately 1500 mitochondrial proteins. Since identification of the first mitochondrial disease gene in 1988 a total of 425 genes have been associated with mitochondrial diseases. Mitochondrial dysfunctions can be caused both by pathogenic variants in the mitochondrial DNA or the nuclear DNA. Hence, besides maternal inheritance, mitochondrial diseases can follow all modes of Mendelian inheritance. The maternal inheritance and tissue specificity distinguish molecular diagnostics of mitochondrial disorders from other rare disorders. With the advances made in the next-generation sequencing technology, whole exome sequencing and even whole-genome sequencing are now the established methods of choice for molecular diagnostics of mitochondrial diseases. They reach a diagnostic rate of more than 50% in clinically suspected mitochondrial disease patients. Moreover, next-generation sequencing is delivering a constantly growing number of novel mitochondrial disease genes. This chapter reviews mitochondrial and nuclear causes of mitochondrial diseases, molecular diagnostic methodologies, and their current challenges and perspectives.

Mitotic MTH1 Inhibitors in Treatment of Cancer

The DNA damage response (DDR) protein MTH1 is sanitising the oxidized dNTP pool and preventing incorporation of oxidative damage into DNA and has an emerging role in mitosis. It is a stress-induced protein and often found to be overexpressed in cancer. Mitotic MTH1 inhibitors arrest cells in mitosis and result in incorporation of oxidative damage into DNA and selective killing of cancer cells. Here, I discuss the leading mitotic MTH1 inhibitor TH1579 (OXC-101, karonudib), now being evaluated in clinical trials, and describe its dual effect on mitosis and incorporation of oxidative DNA damage in cancer cells. I describe why MTH1 inhibitors that solely inhibits the enzyme activity fail to kill cancer cells and discuss if MTH1 is a valid target for cancer treatment. I discuss emerging roles of MTH1 in regulating tubulin polymerisation and mitosis and the necessity of developing the basic science insights along with translational efforts. I also give a perspective on how edgetic perturbation is making target validation difficult in the DDR field.

Disease-associated H58Y mutation affects the nuclear dynamics of human DNA topoisomerase IIβ

DNA topoisomerase II (TOP2) is an enzyme that resolves DNA topological problems and plays critical roles in various nuclear processes. Recently, a heterozygous H58Y substitution in the ATPase domain of human TOP2B was identified from patients with autism spectrum disorder, but its biological significance remains unclear. In this study, we analyzed the nuclear dynamics of TOP2B with H58Y (TOP2B H58Y). Although wild-type TOP2B was highly mobile in the nucleus of a living cell, the nuclear mobility of TOP2B H58Y was markedly reduced, suggesting that the impact of H58Y manifests as low protein mobility. We found that TOP2B H58Y is insensitive to ICRF-187, a TOP2 inhibitor that halts TOP2 as a closed clamp on DNA. When the ATPase activity of TOP2B was compromised, the nuclear mobility of TOP2B H58Y was restored to wild-type levels, indicating the contribution of the ATPase activity to the low nuclear mobility. Analysis of genome-edited cells harboring TOP2B H58Y showed that TOP2B H58Y retains sensitivity to the TOP2 poison etoposide, implying that TOP2B H58Y can undergo at least a part of its catalytic reactions. Collectively, TOP2 H58Y represents a unique example of the relationship between a disease-associated mutation and perturbed protein dynamics.

Targeting the DNA damage response and repair in cancer through nucleotide metabolism

The exploitation of the DNA damage response and DNA repair proficiency of cancer cells is an important anticancer strategy. The replication and repair of DNA are dependent upon the supply of deoxynucleoside triphosphate (dNTP) building blocks, which are produced and maintained by nucleotide metabolic pathways. Enzymes within these pathways can be promising targets to selectively induce toxic DNA lesions in cancer cells. These same pathways also activate antimetabolites, an important group of chemotherapies that disrupt both nucleotide and DNA metabolism to induce DNA damage in cancer cells. Thus, dNTP metabolic enzymes can also be targeted to refine the use of these chemotherapeutics, many of which remain standard of care in common cancers. In this review article, we will discuss both these approaches exemplified by the enzymes MTH1, MTHFD2 and SAMHD1. © 2022 The Authors. Molecular Oncology published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.

In silico analyses of Wnt1 nsSNPs reveal structurally destabilizing variants, altered interactions with Frizzled receptors and its deregulation in tumorigenesis

Wnt1 is the first mammalian Wnt gene, which is discovered as proto-oncogene and in human the gene is located on the chromosome 12q13. Mutations in Wnt1 are reported to be associated with various cancers and other human diseases. The structural and functional consequences of most of the non-synonymous SNPs (nsSNPs), present in the human Wnt1 gene, are not known. In the present work, extensive bioinformatics analyses are used to screen 292 nsSNPs of Wnt1 for predicting pathogenic and harmless polymorphisms. We have identified 10 highly deleterious nsSNPs among which 7 are located within the highly conserved areas. These 10 nsSNPs are also predicted to affect the post-translational modifications of Wnt1. Further, structure based stability analyses of these 10 highly deleterious nsSNPs revealed 8 variants as highly destabilizing. These 8 highly destabilizing variants were shown to have high BC score and high RMSIP score from normal mode analyses. Based on the deformation energies, obtained from the normal mode analyses, variants like G169A, G169S, G331R and G331S were found to be unstable. Molecular Dynamics (MD) simulations revealed structural stability and fluctuation of WT Wnt1 and its prioritized variants. RMSD remained fluctuating mostly between 4 and 5 Å and occasionally between 3.5 and 5.5 Å ranges. RMSF in the CTD region (residues 330-360) of the binding pocket were lower compared to that of WT. Studying the impacts of nsSNPs on the binding interface of Wnt1 and seven Frizzled receptors have predicted substitutions which can stabilize or destabilize the binding interface. We have found that Wnt1 and FZD8-CRD is the best docked complex in our study. MD simulation based analyses of wild type Wnt1-FZD8-CRD complex and the 8 prioritized variants revealed that RMSF was higher in the unstructured regions and RMSD remained fluctuating in the region of 5 Å ± 1 Å. We have also observed differential Wnt1 gene expression pattern in normal, tumor and metastatic conditions across different tissues. Wnt1 gene expression was significantly higher in metastatic tissues of lungs, colon and skin; and was significantly lower in metastatic tissues of breast, esophagus and kidney. We have also found that Wnt1 deregulation is associated with survival outcome in patients with gastric and breast cancer. Furthermore, these computationally screened highly deleterious nsSNPs of Wnt1 can be analyzed in population based genetic studies and may help understand the Wnt1 associated diseases.

Biological interacting units identified in human protein networks reveal tissue-functional diversification and its impact on disease

•

Biological processes are exerted by groups of physically interacting proteins.

•

Proteins display variable biological roles depending on tissue-interactomic context.

•

Tissue-specific protein-protein interaction networks reveal functional diversification.

•

Most disease associated genes/proteins display tissue-specific phenotypes.

•

Protein interaction network analysis is a valuable resource to identify disease genes.

Protein-protein interactions (PPI) play an essential role in the biological processes that occur in the cell. Therefore, the dissection of PPI networks becomes decisive to model functional coordination and predict pathological de-regulation. Cellular networks are dynamic and proteins display varying roles depending on the tissue-interactomic context. Thus, the use of centrality measures in individual proteins fall short to dissect the functional properties of the cell. For this reason, there is a need for more comprehensive, relational, and context-specific ways to analyze the multiple actions of proteins in different cells and identify specific functional assemblies within global biomolecular networks. Under this framework, we define Biological Interacting units (BioInt-U) as groups of proteins that interact physically and are enriched in a common Gene Ontology. A search strategy was applied on 33 tissue-specific (TS) PPI networks to generate BioInt libraries associated with each particular human tissue. The cross-tissue comparison showed that housekeeping assemblies incorporate different proteins and exhibit distinct network properties depending on the tissue. Furthermore, disease genes (DGs) of tissue-associated pathologies preferentially accumulate in units in the expected tissues, which in turn were more central in the TS networks. Overall, the study reveals a tissue-specific functional diversification based on the identification of specific protein units and suggests vulnerabilities specific of each tissue network, which can be applied to refine protein-disease association methods.

Structural Determinants of Yeast Protein-Protein Interaction Interface Evolution at the Residue Level

Interfaces of contact between proteins play important roles in determining the proper structure and function of protein-protein interactions (PPIs). Therefore, to fully understand PPIs, we need to better understand the evolutionary design principles of PPI interfaces. Previous studies have uncovered that interfacial sites are more evolutionarily conserved than other surface protein sites. Yet, little is known about the nature and relative importance of evolutionary constraints in PPI interfaces. Here, we explore constraints imposed by the structure of the microenvironment surrounding interfacial residues on residue evolutionary rate using a large dataset of over 700 structural models of baker's yeast PPIs. We find that interfacial residues are, on average, systematically more conserved than all other residues with a similar degree of total burial as measured by relative solvent accessibility (RSA). Besides, we find that RSA of the residue when the PPI is formed is a better predictor of interfacial residue evolutionary rate than RSA in the monomer state. Furthermore, we investigate four structure-based measures of residue interfacial involvement, including change in RSA upon binding (ΔRSA), number of residue-residue contacts across the interface, and distance from the center or the periphery of the interface. Integrated modeling for evolutionary rate prediction in interfaces shows that ΔRSA plays a dominant role among the four measures of interfacial involvement, with minor, but independent contributions from other measures. These results yield insight into the evolutionary design of interfaces, improving our understanding of the role that structure plays in the molecular evolution of PPIs at the residue level.

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.

Mapping the energetic and allosteric landscapes of protein binding domains

Allosteric communication between distant sites in proteins is central to biological regulation but still poorly characterized, limiting understanding, engineering and drug development^1-6. An important reason for this is the lack of methods to comprehensively quantify allostery in diverse proteins. Here we address this shortcoming and present a method that uses deep mutational scanning to globally map allostery. The approach uses an efficient experimental design to infer en masse the causal biophysical effects of mutations by quantifying multiple molecular phenotypes-here we examine binding and protein abundance-in multiple genetic backgrounds and fitting thermodynamic models using neural networks. We apply the approach to two of the most common protein interaction domains found in humans, an SH3 domain and a PDZ domain, to produce comprehensive atlases of allosteric communication. Allosteric mutations are abundant, with a large mutational target space of network-altering 'edgetic' variants. Mutations are more likely to be allosteric closer to binding interfaces, at glycine residues and at specific residues connecting to an opposite surface within the PDZ domain. This general approach of quantifying mutational effects for multiple molecular phenotypes and in multiple genetic backgrounds should enable the energetic and allosteric landscapes of many proteins to be rapidly and comprehensively mapped.

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.

The properties of human disease mutations at protein interfaces

The assembly of proteins into complexes and their interactions with other biomolecules are often vital for their biological function. While it is known that mutations at protein interfaces have a high potential to be damaging and cause human genetic disease, there has been relatively little consideration for how this varies between different types of interfaces. Here we investigate the properties of human pathogenic and putatively benign missense variants at homomeric (isologous and heterologous), heteromeric, DNA, RNA and other ligand interfaces, and at different regions in proteins with respect to those interfaces. We find that different types of interfaces vary greatly in their propensity to be associated with pathogenic mutations, with homomeric heterologous and DNA interfaces being particularly enriched in disease. We also find that residues that do not directly participate in an interface, but are close in three-dimensional space, show a significant disease enrichment. Finally, we observe that mutations at different types of interfaces tend to have distinct property changes when undergoing amino acid substitutions associated with disease, and that this is linked to substantial variability in their identification by computational variant effect predictors.

Nearly all proteins interact with other molecules as part of their biological function. For example, proteins can interact with other copies of the same type of protein, with different proteins, with DNA, or with small ligand molecules. Many mutations at protein interfaces, the regions of proteins that interact with other molecules, are known to cause human genetic disease. In this study, we first investigate how different types of protein interfaces have different tendencies to be associated with disease. We also show that the closer a mutation is to an interface, the more likely it is to cause disease. Finally, we study how mutations at different types of interfaces tend to be associated with different changes in amino acid properties, which appears to influence our ability to computationally predict the effects of mutations. Ultimately, we hope that consideration of protein interface properties will eventually improve our ability to identify new disease-causing mutations.

Decoding the cellular effects of genetic variation through interaction proteomics

It is often unclear how genetic variation translates into cellular phenotypes, including how much of the coding variation can be recovered in the proteome. Proteogenomic analyses of heterogenous cell lines revealed that the genetic differences impact mostly the abundance and stoichiometry of protein complexes, with the effects propagating post-transcriptionally via protein interactions onto other subunits. Conversely, large scale binary interaction analyses of missense variants revealed that loss of interaction is widespread and caused by about 50% disease-associated mutations, while deep scanning mutagenesis of binary interactions identified thousands of interaction-deficient variants per interaction. The idea that phenotypes arise from genetic variation through protein-protein interaction is therefore substantiated by both forward and reverse interaction proteomics. With improved methodologies, these two approaches combined can close the knowledge gap between nucleotide sequence variation and its functional consequences on the cellular proteome.

Extracting Significant Comorbid Diseases from MeSH Index of PubMed

Text mining is an important research area to be explored in terms of understanding disease associations and have an insight in disease comorbidities. The reason for comorbid occurrence in any patient may be genetic or molecular interference from any other processes. Comorbidity and multimorbidity may be technically different, yet still are inseparable in studies. They have overlapping nature of associations and hence can be integrated for a more rational approach. The association rule generally used to determine comorbidity may also be helpful in novel knowledge prediction or may even serve as an important tool of assessment in surgical cases. Another approach of interest may be to utilize biological vocabulary resources like UMLS/MeSH across a patient health information and analyze the interrelationship between different health conditions. The protocol presented here can be utilized for understanding the disease associations and analyze at an extensive level.

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.

Annotating Cancer-Related Variants at Protein-Protein Interface with Structure-PPi

A comprehensive analysis of germline and somatic variants requires complex computational approaches that combine next-generation sequencing (NGS)-based omics data with curated annotations from public repositories. Here, we describe Structure-PPi, which facilitates the analysis of cancer-related variants onto protein 3D structures, interaction interfaces, and other important functional sites (i.e., catalytic, ligand-binding, posttranslational modification). Our approach relies on features extracted from Interactome3D, UniProtKB, InterPro, APPRIS, dbNSFP, and COSMIC databases and provides complementary information to pathogenicity prediction methods. Thus, Structure-PPi helps in the discrimination of false-positive predictions and adds both mechanistic and biological insights into the role of variants in a given cancer. An online version of the tools is available at https://rbbt.bsc.es/StructurePPI/ .

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.

GPCRs Are Optimal Regulators of Complex Biological Systems and Orchestrate the Interface between Health and Disease

GPCRs arguably represent the most effective current therapeutic targets for a plethora of diseases. GPCRs also possess a pivotal role in the regulation of the physiological balance between healthy and pathological conditions; thus, their importance in systems biology cannot be underestimated. The molecular diversity of GPCR signaling systems is likely to be closely associated with disease-associated changes in organismal tissue complexity and compartmentalization, thus enabling a nuanced GPCR-based capacity to interdict multiple disease pathomechanisms at a systemic level. GPCRs have been long considered as controllers of communication between tissues and cells. This communication involves the ligand-mediated control of cell surface receptors that then direct their stimuli to impact cell physiology. Given the tremendous success of GPCRs as therapeutic targets, considerable focus has been placed on the ability of these therapeutics to modulate diseases by acting at cell surface receptors. In the past decade, however, attention has focused upon how stable multiprotein GPCR superstructures, termed receptorsomes, both at the cell surface membrane and in the intracellular domain dictate and condition long-term GPCR activities associated with the regulation of protein expression patterns, cellular stress responses and DNA integrity management. The ability of these receptorsomes (often in the absence of typical cell surface ligands) to control complex cellular activities implicates them as key controllers of the functional balance between health and disease. A greater understanding of this function of GPCRs is likely to significantly augment our ability to further employ these proteins in a multitude of diseases.

De novo individualized disease modules reveal the synthetic penetrance of genes and inform personalized treatment regimens

Current understandings of individual disease etiology and therapeutics are limited despite great need. To fill the gap, we propose a novel computational pipeline that collects potent disease gene cooperative pathways to envision individualized disease etiology and therapies. Our algorithm constructs individualized disease modules de novo, which enables us to elucidate the importance of mutated genes in specific patients and to understand the synthetic penetrance of these genes across patients. We reveal that importance of the notorious cancer drivers TP53 and PIK3CA fluctuate widely across breast cancers and peak in tumors with distinct numbers of mutations and that rarely mutated genes such as XPO1 and PLEKHA1 have high disease module importance in specific individuals. Furthermore, individualized module disruption enables us to devise customized singular and combinatorial target therapies that were highly varied across patients, showing the need for precision therapeutics pipelines. As the first analysis of de novo individualized disease modules, we illustrate the power of individualized disease modules for precision medicine by providing deep novel insights on the activity of diseased genes in individuals.

Identification of discriminative gene-level and protein-level features associated with pathogenic gain-of-function and loss-of-function variants

Identifying whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity is an important step toward understanding disease mechanisms because they may result in markedly different clinical phenotypes. Here, we generated an extensive database of documented germline GOF and LOF pathogenic variants by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. We then investigated various gene- and protein-level features of GOF and LOF variants and applied machine learning and statistical analyses to identify discriminative features. We found that GOF variants were enriched in essential genes, for autosomal-dominant inheritance, and in protein binding and interaction domains, whereas LOF variants were enriched in singleton genes, for protein-truncating variants, and in protein core regions. We developed a user-friendly web-based interface that enables the extraction of selected subsets from the GOF/LOF database by a broad set of annotated features and downloading of up-to-date versions. These results improve our understanding of how variants affect gene/protein function and may ultimately guide future treatment options.

Network-based protein-protein interaction prediction method maps perturbations of cancer interactome

The perturbations of protein-protein interactions (PPIs) were found to be the main cause of cancer. Previous PPI prediction methods which were trained with non-disease general PPI data were not compatible to map the PPI network in cancer. Therefore, we established a novel cancer specific PPI prediction method dubbed NECARE, which was based on relational graph convolutional network (R-GCN) with knowledge-based features. It achieved the best performance with a Matthews correlation coefficient (MCC) = 0.84±0.03 and an F1 = 91±2% compared with other methods. With NECARE, we mapped the cancer interactome atlas and revealed that the perturbations of PPIs were enriched on 1362 genes, which were named cancer hub genes. Those genes were found to over-represent with mutations occurring at protein-macromolecules binding interfaces. Furthermore, over 56% of cancer treatment-related genes belonged to hub genes and they were significantly related to the prognosis of 32 types of cancers. Finally, by coimmunoprecipitation, we confirmed that the NECARE prediction method was highly reliable with a 90% accuracy. Overall, we provided the novel network-based cancer protein-protein interaction prediction method and mapped the perturbation of cancer interactome. NECARE is available at: https://github.com/JiajunQiu/NECARE.

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.

Protein domain-based approaches for the identification and prioritization of therapeutically actionable cancer variants

The tremendous number of cancer variants that can be detected by NGS analyses has required the development of computational approaches to prioritize mutations on the basis of their biological and clinical significance. Standard strategies take a gene-centric approach to the problem, allowing exclusively the identification of highly frequent variants. On the contrary, protein domain (PD)-based approaches allow to identify functionally relevant low frequency variants by searching for mutations that recur on analogous residues across homologous proteins (i.e. containing the same PD). Such approaches enable to transfer information about the effects and druggability from one known mutation to unknown ones. Here we describe how PD-based strategies work, and discuss how they could be exploited for mutation prioritization. The principle that mutations clustered on specific residues of PDs have the same functional consequences and are therapeutically actionable in a similar manner could help the choice of patient-specific targeted drugs, eventually improving the management of cancer patients.

Exploring absent protein function in yeast: assaying post translational modification and human genetic variation

Yeast is a valuable eukaryotic model organism that has evolved many processes conserved up to humans, yet many protein functions, including certain DNA and protein modifications, are absent. It is this absence of protein function that is fundamental to approaches using yeast as an in vivo test system to investigate human proteins. Functionality of the heterologous expressed proteins is connected to a quantitative, selectable phenotype, enabling the systematic analyses of mechanisms and specificity of DNA modification, post-translational protein modifications as well as the impact of annotated cancer mutations and coding variation on protein activity and interaction. Through continuous improvements of yeast screening systems, this is increasingly carried out on a global scale using deep mutational scanning approaches. Here we discuss the applicability of yeast systems to investigate absent human protein function with a specific focus on the impact of protein variation on protein-protein interaction modulation.

Selection of functional EPHB2 genotypes from ENU mutated grass carp treated with GCRV

Background

N-ethyl-N-nitrosourea (ENU) mutagenesis is a useful method for the genetic engineering of plants, and the production of functional mutants in animal models including mice and zebrafish. Grass carp reovirus (GCRV) is a haemorrhagic disease of grass carp which has caused noteworthy losses in fingerlings over the last few years. To overcome this problem, we used ENU mutant grass carp in an attempt to identify functional resistance genes for future hereditary rearing projects in grass carp.

Results

This study used ENU-mutated grass carp to identify genetic markers associated with resistance to the haemorrhagic disease caused by GCRV. Bulked segregant analysis (BSA) was performed on two homozygous gynogenetic ENU grass carp groups who were susceptible or resistant to GCRV. This analysis identified 466,162 SNPs and 197,644 InDels within the genomes of these mixed pools with a total of 170 genes annotated in the associated region, including 49 genes with non-synonymous mutations at SNP sites and 25 genes with frame shift mutations at InDel sites. Of these 170 mutated genes, 5 randomly selected immune-related genes were shown to be more strongly expressed in the resistant group as compared to the susceptible animals. In addition, we found that one immune-related gene, EPHB2, presented with two heterozygous SNP mutations which altered the animal’s responded to GCRV disease. These SNPs were found in the intron region of EPHB2 at positions 5859 (5859^G > A) and 5968 (5968^G > A) and were significantly (p = 0.002, 0.003) associated with resistance to GCRV. These SNP sites were also shown to correlate with the GCRV-resistant phenotype in these ENU grass carp. We also evaluated the mortality of the different ENU fish genotypes in response to GCRV and the SNPs in EPHB2. The outcomes of these evaluations will be useful in future selections of GCRV-resistant genes for genetic breeding in grass carp.

Conclusion

Our results provide a proof of concept for the application of BSA-sequence analysis in detecting genes responsible for specific functional genotypes and may help to develop better methods for marker-assisted selection, especially for disease resistance in response to GCRV.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07858-x.

Yeast proteins do not practice social distancing as species hybridize

With the current COVID-19 pandemic, we all realized how important interactions are. Interactions are everywhere. At the cellular level, protein interactions play a key role and their ensemble, also called interactome, is often referred as the basic building blocks of life. Given its importance, the maintenance of the integrity of the interactome is a real challenge in the cell. Many events during evolution can disrupt interactomes and potentially result in different characteristics for the organisms. However, the molecular underpinnings of changes in interactions at the cellular level are still largely unexplored. Among the perturbations, hybridization puts in contact two different interactomes, which may lead to many changes in the protein interaction network of the hybrid, including gains and losses of interactions. We recently investigated the fate of the interactomes after hybridization between yeast species using a comparative proteomics approach. A large-scale conservation of the interactions was observed in hybrids, but we also noticed the presence of proteostasis-related changes. This suggests that, despite a general robustness, small differences may accumulate in hybrids and perturb their protein physiology. Here, we summarize our work with a broader perspective on the importance of interactions.

Differential Network Analysis: A Statistical Perspective

Networks effectively capture interactions among components of complex systems, and have thus become a mainstay in many scientific disciplines. Growing evidence, especially from biology, suggest that networks undergo changes over time, and in response to external stimuli. In biology and medicine, these changes have been found to be predictive of complex diseases. They have also been used to gain insight into mechanisms of disease initiation and progression. Primarily motivated by biological applications, this article provides a review of recent statistical machine learning methods for inferring networks and identifying changes in their structures.

e-MutPath: computational modeling reveals the functional landscape of genetic mutations rewiring interactome networks

Understanding the functional impact of cancer somatic mutations represents a critical knowledge gap for implementing precision oncology. It has been increasingly appreciated that the interaction profile mediated by a genomic mutation provides a fundamental link between genotype and phenotype. However, specific effects on biological signaling networks for the majority of mutations are largely unknown by experimental approaches. To resolve this challenge, we developed e-MutPath (edgetic Mutation-mediated Pathway perturbations), a network-based computational method to identify candidate ‘edgetic’ mutations that perturb functional pathways. e-MutPath identifies informative paths that could be used to distinguish disease risk factors from neutral elements and to stratify disease subtypes with clinical relevance. The predicted targets are enriched in cancer vulnerability genes, known drug targets but depleted for proteins associated with side effects, demonstrating the power of network-based strategies to investigate the functional impact and perturbation profiles of genomic mutations. Together, e-MutPath represents a robust computational tool to systematically assign functions to genetic mutations, especially in the context of their specific pathway perturbation effect.

Insulin's Discovery: New Insights on Its Hundredth Birthday: From Insulin Action and Clearance to Sweet Networks

In 2021, the 100th anniversary of the isolation of insulin and the rescue of a child with type 1 diabetes from death will be marked. In this review, we highlight advances since the ingenious work of the four discoverers, Frederick Grant Banting, John James Rickard Macleod, James Bertram Collip and Charles Herbert Best. Macleoad closed his Nobel Lecture speech by raising the question of the mechanism of insulin action in the body. This challenge attracted many investigators, and the question remained unanswered until the third part of the 20th century. We summarize what has been learned, from the discovery of cell surface receptors, insulin action, and clearance, to network and precision medicine.

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.

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.