High-resolution network biology: connecting sequence with function

Broad functional profiling of fission yeast proteins using phenomics and machine learning

Many proteins remain poorly characterized even in well-studied organisms, presenting a bottleneck for research. We applied phenomics and machine-learning approaches with Schizosaccharomyces pombe for broad cues on protein functions. We assayed colony-growth phenotypes to measure the fitness of deletion mutants for 3509 non-essential genes in 131 conditions with different nutrients, drugs, and stresses. These analyses exposed phenotypes for 3492 mutants, including 124 mutants of 'priority unstudied' proteins conserved in humans, providing varied functional clues. For example, over 900 proteins were newly implicated in the resistance to oxidative stress. Phenotype-correlation networks suggested roles for poorly characterized proteins through 'guilt by association' with known proteins. For complementary functional insights, we predicted Gene Ontology (GO) terms using machine learning methods exploiting protein-network and protein-homology data (NET-FF). We obtained 56,594 high-scoring GO predictions, of which 22,060 also featured high information content. Our phenotype-correlation data and NET-FF predictions showed a strong concordance with existing PomBase GO annotations and protein networks, with integrated analyses revealing 1675 novel GO predictions for 783 genes, including 47 predictions for 23 priority unstudied proteins. Experimental validation identified new proteins involved in cellular aging, showing that these predictions and phenomics data provide a rich resource to uncover new protein functions.

The RNA interactome in the Hallmarks of Cancer

Ribonucleic acid (RNA) molecules are indispensable for cellular homeostasis in healthy and malignant cells. However, the functions of RNA extend well beyond that of a protein-coding template. Rather, both coding and non-coding RNA molecules function through critical interactions with a plethora of cellular molecules, including other RNAs, DNA, and proteins. Deconvoluting this RNA interactome, including the interacting partners, the nature of the interaction, and dynamic changes of these interactions in malignancies has yielded fundamental advances in knowledge and are emerging as a novel therapeutic strategy in cancer. Here, we present an RNA-centric review of recent advances in the field of RNA-RNA, RNA-protein, and RNA-DNA interactomic network analysis and their impact across the Hallmarks of Cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

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.

Proteomic profiling reveals the potential mechanisms and regulatory targets of sirtuin 4 in 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine-induced Parkinson's mouse model

Introduction

Parkinson's disease (PD), as a common neurodegenerative disease, currently has no effective therapeutic approaches to delay or stop its progression. There is an urgent need to further define its pathogenesis and develop new therapeutic targets. An increasing number of studies have shown that members of the sirtuin (SIRT) family are differentially involved in neurodegenerative diseases, indicating their potential to serve as targets in therapeutic strategies. Mitochondrial SIRT4 possesses multiple enzymatic activities, such as deacetylase, ADP ribosyltransferase, lipoamidase, and deacylase activities, and exhibits different enzymatic activities and target substrates in different tissues and cells; thus, mitochondrial SIRT4 plays an integral role in regulating metabolism. However, the role and mechanism of SIRT4 in PD are not fully understood. This study aimed to investigate the potential mechanism and possible regulatory targets of SIRT4 in 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced PD mice.

Methods

The expression of the SIRT4 protein in the MPTP-induced PD mouse mice or key familial Parkinson disease protein 7 knockout (DJ-1 KO) rat was compared against the control group by western blot assay. Afterwards, quantitative proteomics and bioinformatics analyses were performed to identify altered proteins in the vitro model and reveal the possible functional role of SIRT4. The most promising molecular target of SIRT4 were screened and validated by viral transfection, western blot assay and reverse transcription quantitative PCR (RT-qPCR) assays.

Results

The expression of the SIRT4 protein was found to be altered both in the MPTP-induced PD mouse mice and DJ-1KO rats. Following the viral transfection of SIRT4, a quantitative proteomics analysis identified 5,094 altered proteins in the vitro model, including 213 significantly upregulated proteins and 222 significantly downregulated proteins. The results from bioinformatics analyses indicated that SIRT4 mainly affected the ribosomal pathway, propionate metabolism pathway, peroxisome proliferator-activated receptor (PPAR) signaling pathway and peroxisome pathway in cells, and we screened 25 potential molecular targets. Finally, only fatty acid binding protein 4 (FABP4) in the PPAR signaling pathway was regulated by SIRT4 among the 25 molecules. Importantly, the alterations in FABP4 and PPARγ were verified in the MPTP-induced PD mouse model.

Discussion

Our results indicated that FABP4 in the PPAR signaling pathway is the most promising molecular target of SIRT4 in an MPTP-induced mouse model and revealed the possible functional role of SIRT4. This study provides a reference for future drug development and mechanism research with SIRT4 as a target or biomarker.

From systems to structure - using genetic data to model protein structures

Understanding the effects of genetic variation is a fundamental problem in biology that requires methods to analyse both physical and functional consequences of sequence changes at systems-wide and mechanistic scales. To achieve a systems view, protein interaction networks map which proteins physically interact, while genetic interaction networks inform on the phenotypic consequences of perturbing these protein interactions. Until recently, understanding the molecular mechanisms that underlie these interactions often required biophysical methods to determine the structures of the proteins involved. The past decade has seen the emergence of new approaches based on coevolution, deep mutational scanning and genome-scale genetic or chemical-genetic interaction mapping that enable modelling of the structures of individual proteins or protein complexes. Here, we review the emerging use of large-scale genetic datasets and deep learning approaches to model protein structures and their interactions, and discuss the integration of structural data from different sources.

Effects of Hydrophobic Residues on the Intracellular Self-Assembly of De Novo Designed Peptide Tags and Their Orthogonality

Protein assemblies forming nano- to micro-sized structures underlie versatile biological events in living systems. For mimicking and engineering these protein assemblies through a bottom-up approach, self-assembling peptides (SAPs) that form nanofibril structures via β-sheets serve as potential practical tags. Nevertheless, the development of SAP tags is still in its infancy, and insight into the relationship between peptide sequences and intracellular self-assembly is limited. In this study, we focused on hydrophobic residues in SAPs and examined the self-assembly of SAP-tagged superfolder GFPs (green fluorescent proteins) in COS-7 cells. Based on XEXK (X; hydrophobic amino acids: F, L, I, V, W, or Y) sequence units, we designed a panel of Xn peptides with different hydrophobic residues (X) and chain lengths (n). We observed that the self-assembly propensity, the size of the assemblies, the influence on protein denaturation, and the subcellular localization differed significantly depending on the hydrophobic amino acid. F9, L9, I7, and V13 peptides formed μm-scaled granules, W13 formed small oligomeric clusters in the cytoplasm, and Y15 formed assemblies in the nucleus. In addition, we investigated the orthogonality of their interaction. Strikingly, W13- and Y15-tagged proteins interacted independently and formed two distinct assemblies in cells. Herein, we have demonstrated the great opportunities for rationalizing artificial protein assemblies and orthogonal structures in an intracellular context using the designed SAPs.

Pairwise Biological Network Alignment Based on Discrete Bat Algorithm

The development of high-throughput technology has provided a reliable technical guarantee for an increased amount of available data on biological networks. Network alignment is used to analyze these data to identify conserved functional network modules and understand evolutionary relationships across species. Thus, an efficient computational network aligner is needed for network alignment. In this paper, the classic bat algorithm is discretized and applied to the network alignment. The bat algorithm initializes the population randomly and then searches for the optimal solution iteratively. Based on the bat algorithm, the global pairwise alignment algorithm BatAlign is proposed. In BatAlign, the individual velocity and the position are represented by a discrete code. BatAlign uses a search algorithm based on objective function that uses the number of conserved edges as the objective function. The similarity between the networks is used to initialize the population. The experimental results showed that the algorithm was able to match proteins with high functional consistency and reach a relatively high topological quality.

Connecting chemistry and biology through molecular descriptors

Through the representation of small molecule structures as numerical descriptors and the exploitation of the similarity principle, chemoinformatics has made paramount contributions to drug discovery, from unveiling mechanisms of action and repurposing approved drugs to de novo crafting of molecules with desired properties and tailored targets. Yet, the inherent complexity of biological systems has fostered the implementation of large-scale experimental screenings seeking a deeper understanding of the targeted proteins, the disrupted biological processes and the systemic responses of cells to chemical perturbations. After this wealth of data, a new generation of data-driven descriptors has arisen providing a rich portrait of small molecule characteristics that goes beyond chemical properties. Here, we give an overview of biologically relevant descriptors, covering chemical compounds, proteins and other biological entities, such as diseases and cell lines, while aligning them to the major contributions in the field from disciplines, such as natural language processing or computer vision. We now envision a new scenario for chemical and biological entities where they both are translated into a common numerical format. In this computational framework, complex connections between entities can be unveiled by means of simple arithmetic operations, such as distance measures, additions, and subtractions.

The Potential of OMICs Technologies for the Treatment of Immune-Mediated Inflammatory Diseases

Immune-mediated inflammatory diseases (IMIDs), such as inflammatory bowel diseases and inflammatory arthritis (e.g., rheumatoid arthritis, psoriatic arthritis), are marked by increasing worldwide incidence rates. Apart from irreversible damage of the affected tissue, the systemic nature of these diseases heightens the incidence of cardiovascular insults and colitis-associated neoplasia. Only 40–60% of patients respond to currently used standard-of-care immunotherapies. In addition to this limited long-term effectiveness, all current therapies have to be given on a lifelong basis as they are unable to specifically reprogram the inflammatory process and thus achieve a true cure of the disease. On the other hand, the development of various OMICs technologies is considered as “the great hope” for improving the treatment of IMIDs. This review sheds light on the progressive development and the numerous approaches from basic science that gradually lead to the transfer from “bench to bedside” and the implementation into general patient care procedures.

Multi-omics analysis suggests enhanced epileptogenesis in the Cornu Ammonis 3 of the pilocarpine model of mesial temporal lobe epilepsy

Mesial temporal lobe epilepsy (MTLE) is a chronic neurological disorder characterized by the occurrence of seizures, and histopathological abnormalities in the mesial temporal lobe structures, mainly hippocampal sclerosis (HS). We used a multi-omics approach to determine the profile of transcript and protein expression in the dorsal and ventral hippocampal dentate gyrus (DG) and Cornu Ammonis 3 (CA3) in an animal model of MTLE induced by pilocarpine. We performed label-free proteomics and RNAseq from laser-microdissected tissue isolated from pilocarpine-induced Wistar rats. We divided the DG and CA3 into dorsal and ventral areas and analyzed them separately. We performed a data integration analysis and evaluated enriched signaling pathways, as well as the integrated networks generated based on the gene ontology processes. Our results indicate differences in the transcriptomic and proteomic profiles among the DG and the CA3 subfields of the hippocampus. Moreover, our data suggest that epileptogenesis is enhanced in the CA3 region when compared to the DG, with most abnormalities in transcript and protein levels occurring in the CA3. Furthermore, our results show that the epileptogenesis in the pilocarpine model involves predominantly abnormal regulation of excitatory neuronal mechanisms mediated by N-methyl D-aspartate (NMDA) receptors, changes in the serotonin signaling, and neuronal activity controlled by calcium/calmodulin-dependent protein kinase (CaMK) regulation and leucine-rich repeat kinase 2 (LRRK2)/WNT signaling pathways.

mRNAs, proteins and the emerging principles of gene expression control

Gene expression involves transcription, translation and the turnover of mRNAs and proteins. The degree to which protein abundances scale with mRNA levels and the implications in cases where this dependency breaks down remain an intensely debated topic. Here we review recent mRNA-protein correlation studies in the light of the quantitative parameters of the gene expression pathway, contextual confounders and buffering mechanisms. Although protein and mRNA levels typically show reasonable correlation, we describe how transcriptomics and proteomics provide useful non-redundant readouts. Integrating both types of data can reveal exciting biology and is an essential step in refining our understanding of the principles of gene expression control.

Extensive rewiring of the EGFR network in colorectal cancer cells expressing transforming levels of KRASG13D

Protein-protein-interaction networks (PPINs) organize fundamental biological processes, but how oncogenic mutations impact these interactions and their functions at a network-level scale is poorly understood. Here, we analyze how a common oncogenic KRAS mutation (KRASG13D) affects PPIN structure and function of the Epidermal Growth Factor Receptor (EGFR) network in colorectal cancer (CRC) cells. Mapping >6000 PPIs shows that this network is extensively rewired in cells expressing transforming levels of KRASG13D (mtKRAS). The factors driving PPIN rewiring are multifactorial including changes in protein expression and phosphorylation. Mathematical modelling also suggests that the binding dynamics of low and high affinity KRAS interactors contribute to rewiring. PPIN rewiring substantially alters the composition of protein complexes, signal flow, transcriptional regulation, and cellular phenotype. These changes are validated by targeted and global experimental analysis. Importantly, genetic alterations in the most extensively rewired PPIN nodes occur frequently in CRC and are prognostic of poor patient outcomes.

Modeling of Multimolecular Complexes

Macromolecular complexes play a key role in cellular function. Predicting the structure and dynamics of these complexes is one of the key challenges in structural biology. Docking applications have traditionally been used to predict pairwise interactions between proteins. However, few methods exist for modeling multi-protein assemblies. Here we present two methods, CombDock and DockStar, that can predict multi-protein assemblies starting from subunit structural models. CombDock can assemble subunits without any assumptions about the pairwise interactions between subunits, while DockStar relies on the interaction graph or, alternatively, a homology model or a cryo-electron microscopy (EM) density map of the entire complex. We demonstrate the two methods using RNA polymerase II with 12 subunits and TRiC/CCT chaperonin with 16 subunits.

Regulostat Inferelator: a novel network biology platform to uncover molecular devices that predetermine cellular response phenotypes

With the emergence of genome editing technologies and synthetic biology, it is now possible to engineer genetic circuits driving a cell's phenotypic response to a stressor. However, capturing a continuous response, rather than simply a binary ‘on’ or ‘off’ response, remains a bioengineering challenge. No tools currently exist to identify gene candidates responsible for predetermining and fine-tuning cell response phenotypes. To address this gap, we devised a novel Regulostat Inferelator (RSI) algorithm to decipher intrinsic molecular devices or networks that predetermine cellular phenotypic responses. The RSI algorithm is designed to extract gene expression patterns from basal transcriptomic data in order to identify ‘regulostat’ constituent gene pairs, which exhibit rheostat-like mode-of-cooperation capable of fine-tuning cellular response. Our proof-of-concept study provides computational evidence for the existence of regulostats and that these networks predetermine cellular response prior to exposure to a stressor or drug. In addition, our work, for the first time, provides evidence of context-specific, drug–regulostat interactions in predetermining drug response phenotypes in cancer cells. Given RSI-inferred regulostat networks offer insights for prioritizing gene candidates capable of rendering a resistant phenotype sensitive to a given drug, we envision that this tool will be of great value in bioengineering and medicine.

Functional geometry of protein interactomes

MOTIVATION

Protein-protein interactions (PPIs) are usually modeled as networks. These networks have extensively been studied using graphlets, small induced subgraphs capturing the local wiring patterns around nodes in networks. They revealed that proteins involved in similar functions tend to be similarly wired. However, such simple models can only represent pairwise relationships and cannot fully capture the higher-order organization of protein interactomes, including protein complexes.

RESULTS

To model the multi-scale organization of these complex biological systems, we utilize simplicial complexes from computational geometry. The question is how to mine these new representations of protein interactomes to reveal additional biological information. To address this, we define simplets, a generalization of graphlets to simplicial complexes. By using simplets, we define a sensitive measure of similarity between simplicial complex representations that allows for clustering them according to their data types better than clustering them by using other state-of-the-art measures, e.g. spectral distance, or facet distribution distance. We model human and baker's yeast protein interactomes as simplicial complexes that capture PPIs and protein complexes as simplices. On these models, we show that our newly introduced simplet-based methods cluster proteins by function better than the clustering methods that use the standard PPI networks, uncovering the new underlying functional organization of the cell. We demonstrate the existence of the functional geometry in the protein interactome data and the superiority of our simplet-based methods to effectively mine for new biological information hidden in the complexity of the higher-order organization of protein interactomes.

AVAILABILITY AND IMPLEMENTATION

Codes and datasets are freely available at http://www0.cs.ucl.ac.uk/staff/natasa/Simplets/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Barcoded DNA nanostructures for the multiplexed profiling of subcellular protein distribution

Massively parallel DNA sequencing is established, yet high-throughput protein profiling remains challenging. Here, we report a barcoding approach that leverages the combinatorial sequence content and the configurational programmability of DNA nanostructures for high-throughput multiplexed profiling of the subcellular expression and distribution of proteins in whole cells. The barcodes are formed by in situ hybridization of tetrahedral DNA nanostructures and short DNA sequences conjugated with protein-targeting antibodies, and by nanostructure-assisted ligation (either enzymatic or chemical) of the nanostructures and exogenous DNA sequences bound to nanoparticles of different sizes (which cause these localization sequences to differentially distribute across subcellular compartments). Compared with linear DNA barcoding, the nanostructured barcodes enhance the signal by more than 100-fold. By implementing the barcoding approach on a microfluidic device for the analysis of rare patient samples, we show that molecular subtypes of breast cancer can be accurately classified and that subcellular spatial markers of disease aggressiveness can be identified.

Unsupervised classification of multi-omics data during cardiac remodeling using deep learning

Integration of multi-omics in cardiovascular diseases (CVDs) presents high potentials for translational discoveries. By analyzing abundance levels of heterogeneous molecules over time, we may uncover biological interactions and networks that were previously unidentifiable. However, to effectively perform integrative analysis of temporal multi-omics, computational methods must account for the heterogeneity and complexity in the data. To this end, we performed unsupervised classification of proteins and metabolites in mice during cardiac remodeling using two innovative deep learning (DL) approaches. First, long short-term memory (LSTM)-based variational autoencoder (LSTM-VAE) was trained on time-series numeric data. The low-dimensional embeddings extracted from LSTM-VAE were then used for clustering. Second, deep convolutional embedded clustering (DCEC) was applied on images of temporal trends. Instead of a two-step procedure, DCEC performes a joint optimization for image reconstruction and cluster assignment. Additionally, we performed K-means clustering, partitioning around medoids (PAM), and hierarchical clustering. Pathway enrichment analysis using the Reactome knowledgebase demonstrated that DL methods yielded higher numbers of significant biological pathways than conventional clustering algorithms. In particular, DCEC resulted in the highest number of enriched pathways, suggesting the strength of its unified framework based on visual similarities. Overall, unsupervised DL is shown to be a promising analytical approach for integrative analysis of temporal multi-omics.

Machine Learning and Integrative Analysis of Biomedical Big Data

Recent developments in high-throughput technologies have accelerated the accumulation of massive amounts of omics data from multiple sources: genome, epigenome, transcriptome, proteome, metabolome, etc. Traditionally, data from each source (e.g., genome) is analyzed in isolation using statistical and machine learning (ML) methods. Integrative analysis of multi-omics and clinical data is key to new biomedical discoveries and advancements in precision medicine. However, data integration poses new computational challenges as well as exacerbates the ones associated with single-omics studies. Specialized computational approaches are required to effectively and efficiently perform integrative analysis of biomedical data acquired from diverse modalities. In this review, we discuss state-of-the-art ML-based approaches for tackling five specific computational challenges associated with integrative analysis: curse of dimensionality, data heterogeneity, missing data, class imbalance and scalability issues.

Use of genetic and chemical synthetic lethality as probes of complexity in bacterial cell systems

Different conditions and genomic contexts are known to have an impact on gene essentiality and interactions. Synthetic lethal interactions occur when a combination of perturbations, either genetic or chemical, result in a more profound fitness defect than expected based on the effect of each perturbation alone. Synthetic lethality in bacterial systems has long been studied; however, during the past decade, the emerging fields of genomics and chemical genomics have led to an increase in the scale and throughput of these studies. Here, we review the concepts of genomics and chemical genomics in the context of synthetic lethality and their revolutionary roles in uncovering novel biology such as the characterization of genes of unknown function and in antibacterial drug discovery. We provide an overview of the methodologies, examples and challenges of both genetic and chemical synthetic lethal screening platforms. Finally, we discuss how to apply genetic and chemical synthetic lethal approaches to rationalize the synergies of drugs, screen for new and improved antibacterial therapies and predict drug mechanism of action.

Mutations in Disordered Regions Can Cause Disease by Creating Dileucine Motifs

Many disease-causing missense mutations affect intrinsically disordered regions (IDRs) of proteins, but the molecular mechanism of their pathogenicity is enigmatic. Here, we employ a peptide-based proteomic screen to investigate the impact of mutations in IDRs on protein-protein interactions. We find that mutations in disordered cytosolic regions of three transmembrane proteins (GLUT1, ITPR1, and CACNA1H) lead to an increased clathrin binding. All three mutations create dileucine motifs known to mediate clathrin-dependent trafficking. Follow-up experiments on GLUT1 (SLC2A1), the glucose transporter causative of GLUT1 deficiency syndrome, revealed that the mutated protein mislocalizes to intracellular compartments. Mutant GLUT1 interacts with adaptor proteins (APs) in vitro, and knocking down AP-2 reverts the cellular mislocalization and restores glucose transport. A systematic analysis of other known disease-causing variants revealed a significant and specific overrepresentation of gained dileucine motifs in structurally disordered cytosolic domains of transmembrane proteins. Thus, several mutations in disordered regions appear to cause "dileucineopathies."

REFINING CELLULAR PATHWAY MODELS USING AN ENSEMBLE OF HETEROGENEOUS DATA SOURCES

Improving current models and hypotheses of cellular pathways is one of the major challenges of systems biology and functional genomics. There is a need for methods to build on established expert knowledge and reconcile it with results of new high-throughput studies. Moreover, the available sources of data are heterogeneous, and the data need to be integrated in different ways depending on which part of the pathway they are most informative for. In this paper, we introduce a compartment specific strategy to integrate edge, node and path data for refining a given network hypothesis. To carry out inference, we use a local-move Gibbs sampler for updating the pathway hypothesis from a compendium of heterogeneous data sources, and a new network regression idea for integrating protein attributes. We demonstrate the utility of this approach in a case study of the pheromone response MAPK pathway in the yeast S. cerevisiae.

Mutations that prevent methylation of cohesin render sensitivity to DNA damage in S. pombe

The canonical role of cohesin is to mediate sister chromatid cohesion. In addition, cohesin plays important roles in processes such as DNA repair and regulation of gene expression. Mounting evidence suggests that various post-translational modifications, including phosphorylation, acetylation and sumoylation regulate cohesin functions. Our mass spectrometry analysis of cohesin purified from Schizosaccharomyces pombe cells revealed that the cohesin subunit Psm1 is methylated on two evolutionarily conserved lysine residues, K536 and K1200. We found that mutations that prevent methylation of Psm1 K536 and K1200 render sensitivity to DNA-damaging agents and show positive genetic interactions with mutations in genes encoding the Mus81-Eme1 endonuclease. Yeast two-hybrid and co-immunoprecipitation assays showed that there were interactions between subunits of the cohesin and Mus81-Eme1 complexes. We conclude that cohesin is methylated and that mutations that prevent methylation of Psm1 K536 and K1200 show synthetic phenotypes with mutants defective in the homologous recombination DNA repair pathway.

Revisiting Antipsychotic Drug Actions Through Gene Networks Associated With Schizophrenia

OBJECTIVE

Antipsychotic drugs were incidentally discovered in the 1950s, but their mechanisms of action are still not understood. Better understanding of schizophrenia pathogenesis could shed light on actions of current drugs and reveal novel "druggable" pathways for unmet therapeutic needs. Recent genome-wide association studies offer unprecedented opportunities to characterize disease gene networks and uncover drug-disease relationships. Polygenic overlap between schizophrenia risk genes and antipsychotic drug targets has been demonstrated, but specific genes and pathways constituting this overlap are undetermined. Risk genes of polygenic disorders do not operate in isolation but in combination with other genes through protein-protein interactions among gene product.

METHOD

The protein interactome was used to map antipsychotic drug targets (N=88) to networks of schizophrenia risk genes (N=328).

RESULTS

Schizophrenia risk genes were significantly localized in the interactome, forming a distinct disease module. Core genes of the module were enriched for genes involved in developmental biology and cognition, which may have a central role in schizophrenia etiology. Antipsychotic drug targets overlapped with the core disease module and comprised multiple pathways beyond dopamine. Some important risk genes like CHRN, PCDH, and HCN families were not connected to existing antipsychotics but may be suitable targets for novel drugs or drug repurposing opportunities to treat other aspects of schizophrenia, such as cognitive or negative symptoms.

CONCLUSIONS

The network medicine approach provides a platform to collate information of disease genetics and drug-gene interactions to shift focus from development of antipsychotics to multitarget antischizophrenia drugs. This approach is transferable to other diseases.

Integrative omics for health and disease

Advances in omics technologies - such as genomics, transcriptomics, proteomics and metabolomics - have begun to enable personalized medicine at an extraordinarily detailed molecular level. Individually, these technologies have contributed medical advances that have begun to enter clinical practice. However, each technology individually cannot capture the entire biological complexity of most human diseases. Integration of multiple technologies has emerged as an approach to provide a more comprehensive view of biology and disease. In this Review, we discuss the potential for combining diverse types of data and the utility of this approach in human health and disease. We provide examples of data integration to understand, diagnose and inform treatment of diseases, including rare and common diseases as well as cancer and transplant biology. Finally, we discuss technical and other challenges to clinical implementation of integrative omics.

FunGeneNet: a web tool to estimate enrichment of functional interactions in experimental gene sets

Background

Estimation of functional connectivity in gene sets derived from genome-wide or other biological experiments is one of the essential tasks of bioinformatics. A promising approach for solving this problem is to compare gene networks built using experimental gene sets with random networks. One of the resources that make such an analysis possible is CrossTalkZ, which uses the FunCoup database. However, existing methods, including CrossTalkZ, do not take into account individual types of interactions, such as protein/protein interactions, expression regulation, transport regulation, catalytic reactions, etc., but rather work with generalized types characterizing the existence of any connection between network members.

Results

We developed the online tool FunGeneNet, which utilizes the ANDSystem and STRING to reconstruct gene networks using experimental gene sets and to estimate their difference from random networks. To compare the reconstructed networks with random ones, the node permutation algorithm implemented in CrossTalkZ was taken as a basis. To study the FunGeneNet applicability, the functional connectivity analysis of networks constructed for gene sets involved in the Gene Ontology biological processes was conducted. We showed that the method sensitivity exceeds 0.8 at a specificity of 0.95. We found that the significance level of the difference between gene networks of biological processes and random networks is determined by the type of connections considered between objects. At the same time, the highest reliability is achieved for the generalized form of connections that takes into account all the individual types of connections. By taking examples of the thyroid cancer networks and the apoptosis network, it is demonstrated that key participants in these processes are involved in the interactions of those types by which these networks differ from random ones.

Conclusions

FunGeneNet is a web tool aimed at proving the functionality of networks in a wide range of sizes of experimental gene sets, both for different global networks and for different types of interactions. Using examples of thyroid cancer and apoptosis networks, we have shown that the links over-represented in the analyzed network in comparison with the random ones make possible a biological interpretation of the original gene/protein sets. The FunGeneNet web tool for assessment of the functional enrichment of networks is available at http://www-bionet.sscc.ru/fungenenet/.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4474-7) contains supplementary material, which is available to authorized users.

Dissecting Dynamic and Heterogeneous Proteasome Complexes Using In Vivo Cross-Linking-Assisted Affinity Purification and Mass Spectrometry

Protein-protein interactions are essential for protein complex formation and function. Affinity purification coupled with mass spectrometry (AP-MS) is the method of choice for studying protein-protein interactions at the systems level under different physiological conditions. Although effective in capturing stable protein interactions, transient, weak, and/or dynamic interactors are often lost due to extended procedures during conventional AP-MS experiments. To circumvent this problem, we have recently developed XAP (in vivo cross-linking (X)-assisted affinity purification)-MS strategy to better preserve dynamic protein complexes under native lysis conditions. In addition, we have developed XBAP (in vivo cross-linking (X)-assisted bimolecular tandem affinity purification)-MS method by incorporating XAP with bimolecular affinity purification to define dynamic and heterogeneous protein subcomplexes. Here we describe general experimental protocols of XAP- and XBAP-MS to study dynamic protein complexes and their subcomplexes, respectively. Specifically, we present their applications in capturing and identifying proteasome dynamic interactors and ubiquitin receptor (UbR)-proteasome subcomplexes.

DEEPN as an Approach for Batch Processing of Yeast 2-Hybrid Interactions

We adapted the yeast 2-hybrid assay to simultaneously uncover multiple transient protein interactions within a single screen by using a strategy termed DEEPN (dynamic enrichment for evaluation of protein networks). This approach incorporates high-throughput DNA sequencing and computation to follow competition among a plasmid population encoding interacting partners. To demonstrate the capacity of DEEPN, we identify a wide range of ubiquitin-binding proteins, including interactors that we verify biochemically. To demonstrate the specificity of DEEPN, we show that DEEPN allows simultaneous comparison of candidate interactors across multiple bait proteins, allowing differential interactions to be identified. This feature was used to identify interactors that distinguish between GTP- and GDP-bound conformations of Rab5.

Base-resolution stratification of cancer mutations using functional variomics

A complete understanding of human cancer variants requires new methods to systematically and efficiently assess the functional effects of genomic mutations at a large scale. Here, we describe a set of tools to rapidly clone and stratify thousands of cancer mutations at base resolution. This protocol provides a massively parallel pipeline to achieve high stringency and throughput. The approach includes high-throughput generation of mutant clones by Gateway, confirmation of variant identity by barcoding and next-generation sequencing, and stratification of cancer variants by multiplexed interaction profiling. Compared with alternative site-directed mutagenesis methods, our protocol requires less sequencing effort and enables robust statistical calling of allele-specific effects. To ensure the precision of variant interaction profiling, we further describe two complementary methods-a high-throughput enhanced yeast two-hybrid (HT-eY2H) assay and a mammalian-cell-based Gaussia princeps luciferase protein-fragment complementation assay (GPCA). These independent assays with standard controls validate mutational interaction profiles with high quality. This protocol provides experimentally derived guidelines for classifying candidate cancer alleles emerging from whole-genome or whole-exome sequencing projects as 'drivers' or 'passengers'. For ∼100 genomic mutations, the protocol-including target primer design, variant library construction, and sequence verification-can be completed within as little as 2-3 weeks, and cancer variant stratification can be completed within 2 weeks.

A Systems Survey of Progressive Host-Cell Reorganization during Rotavirus Infection

Pathogen invasion is often accompanied by widespread alterations in cellular physiology, which reflects the hijacking of host factors and processes for pathogen entry and replication. Although genetic perturbation screens have revealed the complexity of host factors involved for numerous pathogens, it has remained challenging to temporally define the progression of events in host cell reorganization during infection. We combine high-confidence genome-scale RNAi screening of host factors required for rotavirus infection in human intestinal cells with an innovative approach to infer the trajectory of virus infection from fixed cell populations. This approach reveals a comprehensive network of host cellular processes involved in rotavirus infection and implicates AMPK in initiating the development of a rotavirus-permissive environment. Our work provides a powerful approach that can be generalized to order complex host cellular requirements along a trajectory of cellular reorganization during pathogen invasion.

Omics Data Complementarity Underlines Functional Cross-Communication in Yeast

Mapping the complete functional layout of a cell and understanding the cross-talk between different processes are fundamental challenges. They elude us because of the incompleteness and noisiness of molecular data and because of the computational intractability of finding the exact answer. We perform a simple integration of three types of baker’s yeast omics data to elucidate the functional organization and lines of cross-functional communication. We examine protein–protein interaction (PPI), co-expression (COEX) and genetic interaction (GI) data, and explore their relationship with the gold standard of functional organization, the Gene Ontology (GO). We utilize a simple framework that identifies functional cross-communication lines in each of the three data types, in GO, and collectively in the integrated model of the three omics data types; we present each of them in our new Functional Organization Map (FOM) model. We compare the FOMs of the three omics datasets with the FOM of GO and find that GI is in best agreement with GO, followed COEX and PPI. We integrate the three FOMs into a unified FOM and find that it is in better agreement with the FOM of GO than those of any omics dataset alone, demonstrating functional complementarity of different omics data.

Architecture of the human interactome defines protein communities and disease networks

The physiology of a cell can be viewed as the product of thousands of proteins acting in concert to shape the cellular response. Coordination is achieved in part through networks of protein-protein interactions that assemble functionally related proteins into complexes, organelles, and signal transduction pathways. Understanding the architecture of the human proteome has the potential to inform cellular, structural, and evolutionary mechanisms and is critical to elucidation of how genome variation contributes to disease^1–3. Here, we present BioPlex 2.0 (Biophysical Interactions of ORFEOME-derived complexes), which employs robust affinity purification-mass spectrometry (AP-MS) methodology⁴ to elucidate protein interaction networks and co-complexes nucleated by more than 25% of protein coding genes from the human genome, and constitutes the largest such network to date. With >56,000 candidate interactions, BioPlex 2.0 contains >29,000 previously unknown co-associations and provides functional insights into hundreds of poorly characterized proteins while enhancing network-based analyses of domain associations, subcellular localization, and co-complex formation. Unsupervised Markov clustering (MCL)⁵ of interacting proteins identified more than 1300 protein communities representing diverse cellular activities. Genes essential for cell fitness^6,7 are enriched within 53 communities representing central cellular functions. Moreover, we identified 442 communities associated with more than 2000 disease annotations, placing numerous candidate disease genes into a cellular framework. BioPlex 2.0 exceeds previous experimentally derived interaction networks in depth and breadth, and will be a valuable resource for exploring the biology of incompletely characterized proteins and for elucidating larger-scale patterns of proteome organization.

The Key Genes of Chronic Pancreatitis which Bridge Chronic Pancreatitis and Pancreatic Cancer Can be Therapeutic Targets

An important question in systems biology is what role the underlying molecular mechanisms play in disease progression. The relationship between chronic pancreatitis and pancreatic cancer needs further exploration in a system view. We constructed the disease network based on gene expression data and protein-protein interaction. We proposed an approach to discover the underlying core network and molecular factors in the progression of pancreatic diseases, which contain stages of chronic pancreatitis and pancreatic cancer. The chronic pancreatitis and pancreatic cancer core network and key factors were revealed and then verified by gene set enrichment analysis of pathways and diseases. The key factors provide the microenvironment for tumor initiation and the change of gene expression level of key factors bridge chronic pancreatitis and pancreatic cancer. Some new candidate genes need further verification by experiments. Transcriptome profiling-based network analysis reveals the importance of chronic pancreatitis genes and pathways in pancreatic cancer development on a system level by computational method and they can be therapeutic targets.

Unified Alignment of Protein-Protein Interaction Networks

Paralleling the increasing availability of protein-protein interaction (PPI) network data, several network alignment methods have been proposed. Network alignments have been used to uncover functionally conserved network parts and to transfer annotations. However, due to the computational intractability of the network alignment problem, aligners are heuristics providing divergent solutions and no consensus exists on a gold standard, or which scoring scheme should be used to evaluate them. We comprehensively evaluate the alignment scoring schemes and global network aligners on large scale PPI data and observe that three methods, HUBALIGN, L-GRAAL and NATALIE, regularly produce the most topologically and biologically coherent alignments. We study the collective behaviour of network aligners and observe that PPI networks are almost entirely aligned with a handful of aligners that we unify into a new tool, Ulign. Ulign enables complete alignment of two networks, which traditional global and local aligners fail to do. Also, multiple mappings of Ulign define biologically relevant soft clusterings of proteins in PPI networks, which may be used for refining the transfer of annotations across networks. Hence, PPI networks are already well investigated by current aligners, so to gain additional biological insights, a paradigm shift is needed. We propose such a shift come from aligning all available data types collectively rather than any particular data type in isolation from others.

Functional variomics and network perturbation: connecting genotype to phenotype in cancer

Proteins interact with other macromolecules in complex cellular networks for signal transduction and biological function. In cancer, genetic aberrations have been traditionally thought to disrupt the entire gene function. It has been increasingly appreciated that each mutation of a gene could have a subtle but unique effect on protein function or network rewiring, contributing to diverse phenotypic consequences across cancer patient populations. In this Review, we discuss the current understanding of cancer genetic variants, including the broad spectrum of mutation classes and the wide range of mechanistic effects on gene function in the context of signalling networks. We highlight recent advances in computational and experimental strategies to study the diverse functional and phenotypic consequences of mutations at the base-pair resolution. Such information is crucial to understanding the complex pleiotropic effect of cancer genes and provides a possible link between genotype and phenotype in cancer.

Genetic interactions and functional analyses of the fission yeast gsk3 and amk2 single and double mutants defective in TORC1-dependent processes

The Target of Rapamycin (TOR) signalling network plays important roles in aging and disease. The AMP-activated protein kinase (AMPK) and the Gsk3 kinase inhibit TOR during stress. We performed genetic interaction screens using synthetic genetic arrays (SGA) with gsk3 and amk2 as query mutants, the latter encoding the regulatory subunit of AMPK. We identified 69 negative and 82 positive common genetic interactors, with functions related to cellular growth and stress. The 120 gsk3-specific negative interactors included genes functioning in translation and ribosomes. The 215 amk2-specific negative interactors included genes functioning in chromatin silencing and DNA damage repair. Both amk2- and gsk3-specific interactors were enriched in phenotype categories related to abnormal cell size and shape. We also performed SGA screen with the amk2 gsk3 double mutant as a query. Mutants sensitive to 5-fluorouracil, an anticancer drug are under-represented within the 305 positive interactors specific for the amk2 gsk3 query. The triple-mutant SGA screen showed higher number of negative interactions than the double mutant SGA screens and uncovered additional genetic network information. These results reveal common and specialized roles of AMPK and Gsk3 in mediating TOR-dependent processes, indicating that AMPK and Gsk3 act in parallel to inhibit TOR function in fission yeast.

The Resistome: A Comprehensive Database of Escherichia coli Resistance Phenotypes

The microbial ability to resist stressful environmental conditions and chemical inhibitors is of great industrial and medical interest. Much of the data related to mutation-based stress resistance, however, is scattered through the academic literature, making it difficult to apply systematic analyses to this wealth of information. To address this issue, we introduce the Resistome database: a literature-curated collection of Escherichia coli genotypes-phenotypes containing over 5,000 mutants that resist hundreds of compounds and environmental conditions. We use the Resistome to understand our current state of knowledge regarding resistance and to detect potential synergy or antagonism between resistance phenotypes. Our data set represents one of the most comprehensive collections of genomic data related to resistance currently available. Future development will focus on the construction of a combined genomic-transcriptomic-proteomic framework for understanding E. coli's resistance biology. The Resistome can be downloaded at https://bitbucket.org/jdwinkler/resistome_release/overview .

Big data mining powers fungal research: recent advances in fission yeast systems biology approaches

Biology research has entered into big data era. Systems biology approaches therefore become the powerful tools to obtain the whole landscape of how cell separate, grow, and resist the stresses. Fission yeast Schizosaccharomyces pombe is wonderful unicellular eukaryote model, especially studying its division and metabolism can facilitate to understanding the molecular mechanism of cancer and discovering anticancer agents. In this perspective, we discuss the recent advanced fission yeast systems biology tools, mainly focus on metabolomics profiling and metabolic modeling, protein-protein interactome and genetic interaction network, DNA sequencing and applications, and high-throughput phenotypic screening. We therefore hope this review can be useful for interested fungal researchers as well as bioformaticians.

Personalized Therapy of Hypertension: the Past and the Future

During the past 20 years, the studies on genetics or pharmacogenomics of primary hypertension provided interesting results supporting the role of genetics, but no actionable finding ready to be translated into personalized medicine. Two types of approaches have been applied: a "hypothesis-driven" approach on the candidate genes, coding for proteins involved in the biochemical machinery underlying the regulation of BP, and an "unbiased hypothesis-free" approach with GWAS, based on the randomness principles of frequentist statistics. During the past 10-15 years, the application of the latter has overtaken the application of the former leading to an enlargement of the number of previously unknown candidate loci or genes but without any actionable result for the therapy of hypertension. In the present review, we summarize the results of our hypothesis-driven approach based on studies carried out in rats with genetic hypertension and in humans with essential hypertension at the pre-hypertensive and early hypertensive stages. These studies led to the identification of mutant adducin and endogenous ouabain as candidate genetic-molecular mechanisms in both species. Rostafuroxin has been developed for its ability to selectively correct Na(+) pump abnormalities sustained by the two abovementioned mechanisms and to selectively reduce BP in rats and in humans carrying the gene variants underlying the mutant adducin and endogenous ouabain (EO) effects. A clinical trial is ongoing to substantiate these findings. Future studies should apply both the candidate gene and GWAS approaches to fully exploit the potential of genetics in optimizing the personalized therapy.

Nonsense-mediated decay regulates key components of homologous recombination

Cells frequently experience DNA damage that requires repair by homologous recombination (HR). Proteins involved in HR are carefully coordinated to ensure proper and efficient repair without interfering with normal cellular processes. In Saccharomyces cerevisiae, Rad55 functions in the early steps of HR and is regulated in response to DNA damage through phosphorylation by the Mec1 and Rad53 kinases of the DNA damage response. To further identify regulatory processes that target HR, we performed a high-throughput genetic interaction screen with RAD55 phosphorylation site mutants. Genes involved in the mRNA quality control process, nonsense-mediated decay (NMD), were found to genetically interact with rad55 phospho-site mutants. Further characterization revealed that RAD55 transcript and protein levels are regulated by NMD. Regulation of HR by NMD extends to multiple targets beyond RAD55, including RAD51, RAD54 and RAD57 Finally, we demonstrate that loss of NMD results in an increase in recombination rates and resistance to the DNA damaging agent methyl methanesulfonate, suggesting this pathway negatively regulates HR under normal growth conditions.

Using biological networks to integrate, visualize and analyze genomics data

Network biology is a rapidly developing area of biomedical research and reflects the current view that complex phenotypes, such as disease susceptibility, are not the result of single gene mutations that act in isolation but are rather due to the perturbation of a gene’s network context. Understanding the topology of these molecular interaction networks and identifying the molecules that play central roles in their structure and regulation is a key to understanding complex systems. The falling cost of next-generation sequencing is now enabling researchers to routinely catalogue the molecular components of these networks at a genome-wide scale and over a large number of different conditions. In this review, we describe how to use publicly available bioinformatics tools to integrate genome-wide ‘omics’ data into a network of experimentally-supported molecular interactions. In addition, we describe how to visualize and analyze these networks to identify topological features of likely functional relevance, including network hubs, bottlenecks and modules. We show that network biology provides a powerful conceptual approach to integrate and find patterns in genome-wide genomic data but we also discuss the limitations and caveats of these methods, of which researchers adopting these methods must remain aware.

Association analysis of the perturbation of interactions in biological pathways and anticancer drug activity

Understanding how different genomic mutational landscapes in patients with cancer lead to different responses to anticancer drugs is an important challenge for realizing precision medicine for cancer. Many studies have analyzed the comprehensive anticancer drug-response profiles and genomic profiles of cancer cell lines to identify the relationship between the anticancer drug response and genomic alternations. However, few studies have focused on interpreting these profiles with a network perspective. In this work, we analyzed genomic alterations in cancer cell lines by considering which interactions in the signaling pathway were perturbed by mutations. With our interaction-centric approach, we identified novel interaction/drug response associations for two drugs (afatinib and ixabepilone) for which no gene-centric association could be found. When we compared the performance of classifiers for predicting the responses to 164 drugs, the classifiers trained with interaction-centric features outperformed the classifiers trained with gene-centric features, despite the smaller number of features (p-value = 2.0 × 10(-3)). By incorporating the interaction information from signaling pathways, we revealed associations between genomic alterations and drug responses that could be missed when using a gene-centric approach.

Improving the Accuracy of Fitted Atomic Models in Cryo-EM Density Maps of Protein Assemblies Using Evolutionary Information from Aligned Homologous Proteins

Cryo-Electron Microscopy (cryo-EM) has become an important technique to obtain structural insights into large macromolecular assemblies. However the resolution of the density maps do not allow for its interpretation at atomic level. Hence they are combined with high resolution structures along with information from other experimental or bioinformatics techniques to obtain pseudo-atomic models. Here, we describe the use of evolutionary conservation of residues as obtained from protein structures and alignments of homologous proteins to detect errors in the fitting of atomic structures as well as improve accuracy of the protein-protein interfacial regions in the cryo-EM density maps.

The post-genomic era of biological network alignment

Biological network alignment aims to find regions of topological and functional (dis)similarities between molecular networks of different species. Then, network alignment can guide the transfer of biological knowledge from well-studied model species to less well-studied species between conserved (aligned) network regions, thus complementing valuable insights that have already been provided by genomic sequence alignment. Here, we review computational challenges behind the network alignment problem, existing approaches for solving the problem, ways of evaluating their alignment quality, and the approaches' biomedical applications. We discuss recent innovative efforts of improving the existing view of network alignment. We conclude with open research questions in comparative biological network research that could further our understanding of principles of life, evolution, disease, and therapeutics.

Cell-based screens and phenomics with fission yeast

Next-generation sequencing approaches have considerably advanced our understanding of genome function and regulation. However, the knowledge of gene function and complex cellular processes remains a challenge and bottleneck in biological research. Phenomics is a rapidly emerging area, which seeks to rigorously characterize all phenotypes associated with genes or gene variants. Such high-throughput phenotyping under different conditions can be a potent approach toward gene function. The fission yeast Schizosaccharomyces pombe (S. pombe) is a proven eukaryotic model organism that is increasingly used for genomewide screens and phenomic assays. In this review, we highlight current large-scale, cell-based approaches used with S. pombe, including computational colony-growth measurements, genetic interaction screens, parallel profiling using barcodes, microscopy-based cell profiling, metabolomic methods and transposon mutagenesis. These diverse methods are starting to offer rich insights into the relationship between genotypes and phenotypes.

A Pan-Cancer Catalogue of Cancer Driver Protein Interaction Interfaces

Despite their importance in maintaining the integrity of all cellular pathways, the role of mutations on protein-protein interaction (PPI) interfaces as cancer drivers has not been systematically studied. Here we analyzed the mutation patterns of the PPI interfaces from 10,028 proteins in a pan-cancer cohort of 5,989 tumors from 23 projects of The Cancer Genome Atlas (TCGA) to find interfaces enriched in somatic missense mutations. To that end we use e-Driver, an algorithm to analyze the mutation distribution of specific protein functional regions. We identified 103 PPI interfaces enriched in somatic cancer mutations. 32 of these interfaces are found in proteins coded by known cancer driver genes. The remaining 71 interfaces are found in proteins that have not been previously identified as cancer drivers even that, in most cases, there is an extensive literature suggesting they play an important role in cancer. Finally, we integrate these findings with clinical information to show how tumors apparently driven by the same gene have different behaviors, including patient outcomes, depending on which specific interfaces are mutated.

Until now, most efforts in cancer genomics have focused on identifying genes and pathways driving tumor development. Although this has been unquestionably a success, as evidenced by the fact that we now have an extensive catalogue of cancer driver genes and pathways, there is still a poor understanding of why patients with the same affected driver genes may have different disease outcomes or drug responses. This is precisely the aim of this work-to show how by considering proteins as multifunctional factories instead of monolithic black boxes, it is possible to identify novel cancer driver genes and propose molecular hypotheses to explain such heterogeneity. To that end we have mapped the mutation profiles of 5,989 cancer patients from TCGA to more than 10,000 protein structures, leading us to identify 103 protein interaction interfaces enriched in somatic mutations. Finally, we have integrated clinical annotations as well as proteomics data to show how tumors apparently driven by the same gene can display different behaviors, including patient outcomes, depending on which specific interfaces are mutated.

Comprehensive assessment of cancer missense mutation clustering in protein structures

Large-scale tumor sequencing projects enabled the identification of many new cancer gene candidates through computational approaches. Here, we describe a general method to detect cancer genes based on significant 3D clustering of mutations relative to the structure of the encoded protein products. The approach can also be used to search for proteins with an enrichment of mutations at binding interfaces with a protein, nucleic acid, or small molecule partner. We applied this approach to systematically analyze the PanCancer compendium of somatic mutations from 4,742 tumors relative to all known 3D structures of human proteins in the Protein Data Bank. We detected significant 3D clustering of missense mutations in several previously known oncoproteins including HRAS, EGFR, and PIK3CA. Although clustering of missense mutations is often regarded as a hallmark of oncoproteins, we observed that a number of tumor suppressors, including FBXW7, VHL, and STK11, also showed such clustering. Beside these known cases, we also identified significant 3D clustering of missense mutations in NUF2, which encodes a component of the kinetochore, that could affect chromosome segregation and lead to aneuploidy. Analysis of interaction interfaces revealed enrichment of mutations in the interfaces between FBXW7-CCNE1, HRAS-RASA1, CUL4B-CAND1, OGT-HCFC1, PPP2R1A-PPP2R5C/PPP2R2A, DICER1-Mg2+, MAX-DNA, SRSF2-RNA, and others. Together, our results indicate that systematic consideration of 3D structure can assist in the identification of cancer genes and in the understanding of the functional role of their mutations.

Characteristics Analyses and Comparisons of the Protein Structure Networks Constructed by Different Methods

Protein structure networks (PSNs) were widely used in analyses of protein structure and function. In this work, we analyzed and compared the characters of PSNs by different methods. The degrees of the different types of the nodes were found to be associated with the amino acid characters, including SAS, secondary structure, hydropathy and the volume of amino acids. It showed that PSNs by the methods of CA10, SC10 and AT5 inherited more amino acid characters and had higher correlations with the original protein structures. And PSNs by these three methods would be powerful tools in understanding the characters of protein structures.