A Compendium of Co-regulated Protein Complexes in Breast Cancer Reveals Collateral Loss Events

A computational framework for the inference of protein complex remodeling from whole-proteome measurements

Protein complexes are responsible for the enactment of most cellular functions. For the protein complex to form and function, its subunits often need to be present at defined quantitative ratios. Typically, global changes in protein complex composition are assessed with experimental approaches that tend to be time consuming. Here, we have developed a computational algorithm for the detection of altered protein complexes based on the systematic assessment of subunit ratios from quantitative proteomic measurements. We applied it to measurements from breast cancer cell lines and patient biopsies and were able to identify strong remodeling of HDAC2 epigenetic complexes in more aggressive forms of cancer. The presented algorithm is available as an R package and enables the inference of changes in protein complex states by extracting functionally relevant information from bottom-up proteomic datasets.

Proteogenomic data and resources for pan-cancer analysis

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Native top-down proteomics reveals EGFR-ERα signaling crosstalk in breast cancer cells dissociates NUTF2 dimers to modulate ERα signaling and cell growth

Oligomerization of proteins and their modified forms (proteoforms) produces functional protein complexes 1,2. Complexoforms are complexes that consist of the same set of proteins with different proteoforms 3. The ability to characterize these assemblies within cells is critical to understanding the molecular mechanisms involved in disease and to designing effective drugs. An outstanding biological question is how proteoforms drive function and oligomerization of complexoforms. However, tools to define endogenous proteoform-proteoform/ligand interactions are scarce 4. Here, we present a native top-down proteomics (nTDP) strategy that combines size-exclusion chromatography, nano liquid-chromatography in direct infusion mode, field asymmetric ion mobility spectrometry, and multistage mass spectrometry to identify protein assemblies (≤70 kDa) in breast cancer cells and in cells that overexpress EGFR, a resistance model of estrogen receptor-α (ER-α) targeted therapies. By identifying ~104 complexoforms from 17 protein complexes, our nTDP approach revealed several molecular features of the breast cancer proteome, including EGFR-induced dissociation of nuclear transport factor 2 (NUTF2) assemblies that modulate ER activity. Our findings show that the K4 and K55 posttranslational modification sites discovered with nTDP differentially impact the effects of NUTF2 on the inhibition of the ER signaling pathway. By characterizing endogenous proteoform-proteoform/ligand interactions, we reveal the molecular diversity of complexoforms, which allows us to propose a model for ER drug discovery in the context of designing effective inhibitors to selectively bind and disrupt the actions of targeted ER complexoforms.

A NuRD for all seasons

The nucleosome-remodeling and deacetylase (NuRD) complex is an essential transcriptional regulator in all complex animals. All seven core subunits of the complex exist as multiple paralogs, raising the question of whether the complex might utilize paralog switching to achieve cell type-specific functions. We examine the evidence for this idea, making use of published quantitative proteomic data to dissect NuRD composition in 20 different tissues, as well as a large-scale CRISPR knockout screen carried out in >1000 human cancer cell lines. These data, together with recent reports, provide strong support for the idea that distinct permutations of the NuRD complex with tailored functions might regulate tissue-specific gene expression programs.

Protocol to analyze dysregulation of the eIF4F complex in human cancers using R software and large public datasets

Understanding dysregulation of the eukaryotic initiation factor 4F (eIF4F) complex across tumor types is critical to cancer treatment development. We present a protocol and accompanying R package "eIF4F.analysis". We describe analysis of copy number status, gene abundance and stoichiometry, survival probability, expression covariation, correlating genes, mRNA/protein correlation, and protein co-expression. Using publicly available large multi-omics data, eIF4F.analysis permits computationally derived and statistically powerful inferences regarding initiation factor regulation in human cancers and clinical relevance of protein interactions within the eIF4F complex. For complete details on the use and execution of this protocol, please refer to Wu and Wagner (2021).¹.

Experimental reproducibility limits the correlation between mRNA and protein abundances in tumor proteomic profiles

Large-scale studies of human proteomes have revealed only a moderate correlation between mRNA and protein abundances. It is unclear to what extent this moderate correlation reflects post-transcriptional regulation and to what extent it reflects measurement error. Here, by analyzing replicate profiles of tumors and cell lines, we show that there is considerable variation in the reproducibility of measurements of transcripts and proteins from individual genes. Proteins with more reproducible measurements tend to have a higher mRNA-protein correlation, suggesting that measurement reproducibility accounts for a substantial fraction of the unexplained variation between mRNA and protein abundances. The reproducibility of individual proteins is somewhat consistent across studies, and we exploit this to develop an aggregate reproducibility score that explains a substantial amount of the variation in mRNA-protein correlations across multiple studies. Finally, we show that pathways previously reported to have a higher-than-average mRNA-protein correlation may simply contain members that can be more reproducibly quantified.

Pan-cancer proteomic map of 949 human cell lines

The proteome provides unique insights into disease biology beyond the genome and transcriptome. A lack of large proteomic datasets has restricted the identification of new cancer biomarkers. Here, proteomes of 949 cancer cell lines across 28 tissue types are analyzed by mass spectrometry. Deploying a workflow to quantify 8,498 proteins, these data capture evidence of cell-type and post-transcriptional modifications. Integrating multi-omics, drug response, and CRISPR-Cas9 gene essentiality screens with a deep learning-based pipeline reveals thousands of protein biomarkers of cancer vulnerabilities that are not significant at the transcript level. The power of the proteome to predict drug response is very similar to that of the transcriptome. Further, random downsampling to only 1,500 proteins has limited impact on predictive power, consistent with protein networks being highly connected and co-regulated. This pan-cancer proteomic map (ProCan-DepMapSanger) is a comprehensive resource available at https://cellmodelpassports.sanger.ac.uk.

Mechanistic insights on cytotoxicity of KOLR, Cinnamomum pauciflorum Nees leaf derived active ingredient, by targeting signaling complexes of phosphodiesterase 3B and rap guanine nucleotide exchange factor 3

Protein signaling complexes play important roles in prevention of several cancer types and can be used for development of targeted therapy. The roles of signaling complexes of phosphodiesterase 3B (PDE3B) and Rap guanine nucleotide exchange factor 3 (RAPGEF3), which are two important enzymes of cyclic adenosine monophosphate (cAMP) metabolism, in cancer have not been fully explored. In the current study, a natural product Kaempferol-3-O-(3'',4''-di-E-p-coumaroyl)-α-L-rhamnopyranoside designated as KOLR was extracted from Cinnamomum pauciflorum Nees leaves. KOLR exhibited higher cytotoxic effects against BxCP-3 pancreatic cancer cell line. In BxPC-3 cells, the KOLR could enhance the formation of RAPGEF 3/ PDE3B protein complex to inhibit the activation of Rap-1 and PI3K-AKT pathway, thereby promoting cell apoptosis and inhibiting cell metastasis. Mutation of RAPGEF3 G557A or low expression of PDE3B inactivated the binding action of KOLR resulting in KOLR resistance. The findings of this study show that PDE3B/RAPGEF3 complex is a potential therapeutic cancer target.

Degron masking outlines degronons, co-degrading functional modules in the proteome

Effective organization of proteins into functional modules (networks, pathways) requires systems-level coordination between transcription, translation and degradation. Whereas the cooperation between transcription and translation was extensively studied, the cooperative degradation regulation of protein complexes and pathways has not been systematically assessed. Here we comprehensively analyzed degron masking, a major mechanism by which cellular systems coordinate degron recognition and protein degradation. For over 200 substrates with characterized degrons (E3 ligase targeting motifs, ubiquitination sites and disordered proteasomal entry sequences), we demonstrate that degrons extensively overlap with protein-protein interaction sites. Analysis of binding site information and protein abundance comparisons show that regulatory partners effectively outcompete E3 ligases, masking degrons from the ubiquitination machinery. Protein abundance variations between normal and cancer cells highlight the dynamics of degron masking components. Finally, integrative analysis of gene co-expression, half-life correlations and functional relationships between interacting proteins point towards higher-order, co-regulated degradation modules (‘degronons’) in the proteome.

Systematic bioinformatics analysis of cooperative degradation of protein complexes indicates that degrons extensively overlap with protein-protein interaction sites, hiding degrons from ubiquitination machinery and suggesting the existence of co-degrading functional modules in the proteome.

RoDiCE: robust differential protein co-expression analysis for cancer complexome

MOTIVATION

The full spectrum of abnormalities in cancer-associated protein complexes remains largely unknown. Comparing the co-expression structure of each protein complex between tumor and healthy cells may provide insights regarding cancer-specific protein dysfunction. However, the technical limitations of mass spectrometry-based proteomics, including contamination with biological protein variants, causes noise that leads to non-negligible over- (or under-) estimating co-expression.

RESULTS

We propose a robust algorithm for identifying protein complex aberrations in cancer based on differential protein co-expression testing. Our method based on a copula is sufficient for improving identification accuracy with noisy data compared to conventional linear correlation-based approaches. As an application, we use large-scale proteomic data from renal cancer to show that important protein complexes, regulatory signaling pathways and drug targets can be identified. The proposed approach surpasses traditional linear correlations to provide insights into higher-order differential co-expression structures.

AVAILABILITY AND IMPLEMENTATION

https://github.com/ymatts/RoDiCE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Protein Complex Organization Imposes Constraints on Proteome Dysregulation in Cancer

Protein assembly is a highly dynamic process and proteins can interact in different ways and stoichiometries within a complex. The importance of maintaining protein stoichiometry for complex function and avoiding aggregation of orphan subunits has been demonstrated. However, how exactly the organization of proteins into complexes constrains differential protein abundance in extreme cellular conditions like cancer, where a lot of protein abundance changes occur, has not been systematically investigated. To study this, we collected proteomic data made available by the Clinical Proteomic Tumor Analysis Consortium (CPTAC) to quantify proteomic changes during carcinogenesis and systematically tested five interaction types in complexes to investigate which of these features impact on protein abundance correlation patterns in cancer. We found that higher than expected fraction of protein complex subunits does not show changes in their abundances compared to those in the normal samples. Furthermore, we found that the way proteins interact in complexes indeed constrains their co-abundance patterns. Our results highlight the role of the interactions between the proteins and the need of cancer cells to deal with aberrant changes in protein abundance.

Dual proteome-scale networks reveal cell-specific remodeling of the human interactome

Thousands of interactions assemble proteins into modules that impart spatial and functional organization to the cellular proteome. Through affinity-purification mass spectrometry, we have created two proteome-scale, cell-line-specific interaction networks. The first, BioPlex 3.0, results from affinity purification of 10,128 human proteins-half the proteome-in 293T cells and includes 118,162 interactions among 14,586 proteins. The second results from 5,522 immunoprecipitations in HCT116 cells. These networks model the interactome whose structure encodes protein function, localization, and complex membership. Comparison across cell lines validates thousands of interactions and reveals extensive customization. Whereas shared interactions reside in core complexes and involve essential proteins, cell-specific interactions link these complexes, "rewiring" subnetworks within each cell's interactome. Interactions covary among proteins of shared function as the proteome remodels to produce each cell's phenotype. Viewable interactively online through BioPlexExplorer, these networks define principles of proteome organization and enable unknown protein characterization.

Simplified and Unified Access to Cancer Proteogenomic Data

Comprehensive cancer data sets recently generated by the Clinical Proteomic Tumor Analysis Consortium (CPTAC) offer great potential for advancing our understanding of how to combat cancer. These data sets include DNA, RNA, protein, and clinical characterization for tumor and normal samples from large cohorts of many different cancer types. The raw data are publicly available at various Cancer Research Data Commons. However, widespread reuse of these data sets is also facilitated by easy access to the processed quantitative data tables. We have created a data application programming interface (API) to distribute these processed tables, implemented as a Python package called cptac. We implement it such that users who prefer to work in R can easily use our package for data access and then transfer the data into R for analysis. Our package distributes the finalized processed CPTAC data sets in a consistent, up-to-date format. This consistency makes it easy to integrate the data with common graphing, statistical, and machine-learning packages for advanced analysis. Additionally, consistent formatting across all cancer types promotes the investigation of pan-cancer trends. The data API structure of directly streaming data within a programming environment enhances the reproducibility. Finally, with the accompanying tutorials, this package provides a novel resource for cancer research education. View the software documentation at https://paynelab.github.io/cptac/. View the GitHub repository at https://github.com/PayneLab/cptac.

Proteogenomics of Non-smoking Lung Cancer in East Asia Delineates Molecular Signatures of Pathogenesis and Progression

Lung cancer in East Asia is characterized by a high percentage of never-smokers, early onset and predominant EGFR mutations. To illuminate the molecular phenotype of this demographically distinct disease, we performed a deep comprehensive proteogenomic study on a prospectively collected cohort in Taiwan, representing early stage, predominantly female, non-smoking lung adenocarcinoma. Integrated genomic, proteomic, and phosphoproteomic analysis delineated the demographically distinct molecular attributes and hallmarks of tumor progression. Mutational signature analysis revealed age- and gender-related mutagenesis mechanisms, characterized by high prevalence of APOBEC mutational signature in younger females and over-representation of environmental carcinogen-like mutational signatures in older females. A proteomics-informed classification distinguished the clinical characteristics of early stage patients with EGFR mutations. Furthermore, integrated protein network analysis revealed the cellular remodeling underpinning clinical trajectories and nominated candidate biomarkers for patient stratification and therapeutic intervention. This multi-omic molecular architecture may help develop strategies for management of early stage never-smoker lung adenocarcinoma.

Proteome-wide Systems Genetics to Identify Functional Regulators of Complex Traits

Proteomic technologies now enable the rapid quantification of thousands of proteins across genetically diverse samples. Integration of these data with systems-genetics analyses is a powerful approach to identify new regulators of economically important or disease-relevant phenotypes in various populations. In this review, we summarize the latest proteomic technologies and discuss technical challenges for their use in population studies. We demonstrate how the analysis of correlation structure and loci mapping can be used to identify genetic factors regulating functional protein networks and complex traits. Finally, we provide an extensive summary of the use of proteome-wide systems genetics throughout fungi, plant, and animal kingdoms and discuss the power of this approach to identify candidate regulators and drug targets in large human consortium studies.

Divergent organ-specific isogenic metastatic cell lines identified using multi-omics exhibit differential drug sensitivity

Background

Monitoring and treating metastatic progression remains a formidable task due, in part, to an inability to monitor specific differential molecular adaptations that allow the cancer to thrive within different tissue types. Hence, to develop optimal treatment strategies for metastatic disease, an important consideration is the divergence of the metastatic cancer growing in visceral organs from the primary tumor. We had previously reported the establishment of isogenic human metastatic breast cancer cell lines that are representative of the common metastatic sites observed in breast cancer patients.

Methods

Here we have used proteomic, RNAseq, and metabolomic analyses of these isogenic cell lines to systematically identify differences and commonalities in pathway networks and examine the effect on the sensitivity to breast cancer therapeutic agents.

Results

Proteomic analyses indicated that dissemination of cells from the primary tumor sites to visceral organs resulted in cell lines that adapted to growth at each new site by, in part, acquiring protein pathways characteristic of the organ of growth. RNAseq and metabolomics analyses further confirmed the divergences, which resulted in differential efficacies to commonly used FDA approved chemotherapeutic drugs. This model system has provided data that indicates that organ-specific growth of malignant lesions is a selective adaptation and growth process.

Conclusions

The insights provided by these analyses indicate that the rationale of targeted treatment of metastatic disease may benefit from a consideration that the biology of metastases has diverged from the primary tumor biology and using primary tumor traits as the basis for treatment may not be ideal to design treatment strategies.

Development of synthetic lethality in cancer: molecular and cellular classification

Recently, genetically targeted cancer therapies have been a topic of great interest. Synthetic lethality provides a new approach for the treatment of mutated genes that were previously considered unable to be targeted in traditional genotype-targeted treatments. The increasing researches and applications in the clinical setting made synthetic lethality a promising anticancer treatment option. However, the current understandings on different conditions of synthetic lethality have not been systematically assessed and the application of synthetic lethality in clinical practice still faces many challenges. Here, we propose a novel and systematic classification of synthetic lethality divided into gene level, pathway level, organelle level, and conditional synthetic lethality, according to the degree of specificity into its biological mechanism. Multiple preclinical findings of synthetic lethality in recent years will be reviewed and classified under these different categories. Moreover, synthetic lethality targeted drugs in clinical practice will be briefly discussed. Finally, we will explore the essential implications of this classification as well as its prospects in eliminating existing challenges and the future directions of synthetic lethality.

CORUM: the comprehensive resource of mammalian protein complexes-2019

CORUM is a database that provides a manually curated repository of experimentally characterized protein complexes from mammalian organisms, mainly human (67%), mouse (15%) and rat (10%). Given the vital functions of these macromolecular machines, their identification and functional characterization is foundational to our understanding of normal and disease biology. The new CORUM 3.0 release encompasses 4274 protein complexes offering the largest and most comprehensive publicly available dataset of mammalian protein complexes. The CORUM dataset is built from 4473 different genes, representing 22% of the protein coding genes in humans. Protein complexes are described by a protein complex name, subunit composition, cellular functions as well as the literature references. Information about stoichiometry of subunits depends on availability of experimental data. Recent developments include a graphical tool displaying known interactions between subunits. This allows the prediction of structural interconnections within protein complexes of unknown structure. In addition, we present a set of 58 protein complexes with alternatively spliced subunits. Those were found to affect cellular functions such as regulation of apoptotic activity, protein complex assembly or define cellular localization. CORUM is freely accessible at http://mips.helmholtz-muenchen.de/corum/.

Biomolecular Databases and Subnetwork Identification Approaches of Interest to Big Data Community: An Expert Review

Next-generation sequencing approaches and genome-wide studies have become essential for characterizing the mechanisms of human diseases. Consequently, many researchers have applied these approaches to discover the genetic/genomic causes of common complex and rare human diseases, generating multiomics big data that span the continuum of genomics, proteomics, metabolomics, and many other system science fields. Therefore, there is a significant and unmet need for biological databases and tools that enable and empower the researchers to analyze, integrate, and make sense of big data. There are currently large number of databases that offer different types of biological information. In particular, the integration of gene expression profiles and protein-protein interaction networks provides a deeper understanding of the complex multilayered molecular architecture of human diseases. Therefore, there has been a growing interest in developing methodologies that integrate and contextualize big data from molecular interaction networks to identify biomarkers of human diseases at a subnetwork resolution as well. In this expert review, we provide a comprehensive summary of most popular biomolecular databases for molecular interactions (e.g., Biological General Repository for Interaction Datasets, Kyoto Encyclopedia of Genes and Genomes and Search Tool for The Retrieval of Interacting Genes/Proteins), gene-disease associations (e.g., Online Mendelian Inheritance in Man, Disease-Gene Network, MalaCards), and population-specific databases (e.g., Human Genetic Variation Database), and describe some examples of their usage and potential applications. We also present the most recent subnetwork identification approaches and discuss their main advantages and limitations. As the field of data science continues to emerge, the present analysis offers a deeper and contextualized understanding of the available databases in molecular biomedicine.

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

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

Promoter Activity Buffering Reduces the Fitness Cost of Misregulation

Organisms regulate gene expression through changes in the activity of transcription factors (TFs). In yeast, the response of genes to changes in TF activity is generally assumed to be encoded in the promoter. To directly test this assumption, we chose 42 genes and, for each, replaced the promoter with a synthetic inducible promoter and measured how protein expression changes as a function of TF activity. Most genes exhibited gene-specific TF dose-response curves not due to differences in mRNA stability, translation, or protein stability. Instead, most genes have an intrinsic ability to buffer the effects of promoter activity. This can be encoded in the open reading frame and the 3' end of genes and can be implemented by both autoregulatory feedback and by titration of limiting trans regulators. We show experimentally and computationally that, when misexpression of a gene is deleterious, this buffering insulates cells from fitness defects due to misregulation.

Multi-omics Characterization of Interaction-mediated Control of Human Protein Abundance levels

Proteogenomic studies of cancer samples have shown that copy-number variation can be attenuated at the protein level for a large fraction of the proteome, likely due to the degradation of unassembled protein complex subunits. Such interaction-mediated control of protein abundance remains poorly characterized. To study this, we compiled genomic, (phospho)proteomic and structural data for hundreds of cancer samples and find that up to 42% of 8,124 analyzed proteins show signs of post-transcriptional control. We find evidence of interaction-dependent control of protein abundance, correlated with interface size, for 516 protein pairs, with some interactions further controlled by phosphorylation. Finally, these findings in cancer were reflected in variation in protein levels in normal tissues. Importantly, expression differences due to natural genetic variation were increasingly buffered from phenotype differences for highly attenuated proteins. Altogether, this study further highlights the importance of posttranscriptional control of protein abundance in cancer and healthy cells.

Overdosage of Balanced Protein Complexes Reduces Proliferation Rate in Aneuploid Cells

Cells with complex aneuploidies display a wide range of phenotypic abnormalities. However, the molecular basis for this has been mainly studied in trisomic (2n + 1) and disomic (n + 1) cells. To determine how karyotype affects proliferation in cells with complex aneuploidies, we generated 92 2n + x yeast strains in which each diploid cell has between 3 and 12 extra chromosomes. Genome-wide and, for individual protein complexes, proliferation defects are caused by the presence of protein complexes in which all subunits are balanced at the 3-copy level. Proteomics revealed that over 50% of 3-copy members of imbalanced complexes were expressed at only 2n protein levels, whereas members of complexes in which all subunits are stoichiometrically balanced at 3 copies per cell had 3n protein levels. We validated this finding using orthogonal datasets from yeast and from human cancers. Taken together, our study provides an explanation of how aneuploidy affects phenotype.

MOLI: multi-omics late integration with deep neural networks for drug response prediction

MOTIVATION

Historically, gene expression has been shown to be the most informative data for drug response prediction. Recent evidence suggests that integrating additional omics can improve the prediction accuracy which raises the question of how to integrate the additional omics. Regardless of the integration strategy, clinical utility and translatability are crucial. Thus, we reasoned a multi-omics approach combined with clinical datasets would improve drug response prediction and clinical relevance.

RESULTS

We propose MOLI, a multi-omics late integration method based on deep neural networks. MOLI takes somatic mutation, copy number aberration and gene expression data as input, and integrates them for drug response prediction. MOLI uses type-specific encoding sub-networks to learn features for each omics type, concatenates them into one representation and optimizes this representation via a combined cost function consisting of a triplet loss and a binary cross-entropy loss. The former makes the representations of responder samples more similar to each other and different from the non-responders, and the latter makes this representation predictive of the response values. We validate MOLI on in vitro and in vivo datasets for five chemotherapy agents and two targeted therapeutics. Compared to state-of-the-art single-omics and early integration multi-omics methods, MOLI achieves higher prediction accuracy in external validations. Moreover, a significant improvement in MOLI's performance is observed for targeted drugs when training on a pan-drug input, i.e. using all the drugs with the same target compared to training only on drug-specific inputs. MOLI's high predictive power suggests it may have utility in precision oncology.

AVAILABILITY AND IMPLEMENTATION

https://github.com/hosseinshn/MOLI.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Understanding Disease Variants through the Lens of Protein Interactions

High-density interaction mapping of mitochondrial proteins provides clues to molecular mechanisms implicated in the progression of neurological disorders.

Disentangling Genetic and Environmental Effects on the Proteotypes of Individuals

Proteotypes, like genotypes, have been found to vary between individuals in several studies, but consistent molecular functional traits across studies remain to be quantified. In a meta-analysis of 11 proteomics datasets from humans and mice, we use co-variation of proteins in known functional modules across datasets and individuals to obtain a consensus landscape of proteotype variation. We find that individuals differ considerably in both protein complex abundances and stoichiometry. We disentangle genetic and environmental factors impacting these metrics, with genetic sex and specific diets together explaining 13.5% and 11.6% of the observed variation of complex abundance and stoichiometry, respectively. Sex-specific differences, for example, include various proteins and complexes, where the respective genes are not located on sex-specific chromosomes. Diet-specific differences, added to the individual genetic backgrounds, might become a starting point for personalized proteotype modulation toward desired features.

•

Benchmarking of datasets on human and mouse proteotypes

•

Consistent co-variation landscape of functional modules across individuals

•

Protein complexes vary in their stoichiometry across individuals

•

Quantifying effects of genetic sex and specific diets on complexes

Quantitative proteomic analysis of prostate tissue specimens identifies deregulated protein complexes in primary prostate cancer

Background

Prostate cancer (PCa) is the most frequently diagnosed non-skin cancer and a leading cause of mortality among males in developed countries. However, our understanding of the global changes of protein complexes within PCa tissue specimens remains very limited, although it has been well recognized that protein complexes carry out essentially all major processes in living organisms and that their deregulation drives the pathogenesis and progression of various diseases.

Methods

By coupling tandem mass tagging-synchronous precursor selection-mass spectrometry/mass spectrometry/mass spectrometry with differential expression and co-regulation analyses, the present study compared the differences between protein complexes in normal prostate, low-grade PCa, and high-grade PCa tissue specimens.

Results

Globally, a large downregulated putative protein–protein interaction (PPI) network was detected in both low-grade and high-grade PCa, yet a large upregulated putative PPI network was only detected in high-grade but not low-grade PCa, compared with normal controls. To identify specific protein complexes that are deregulated in PCa, quantified proteins were mapped to protein complexes in CORUM (v3.0), a high-quality collection of 4274 experimentally verified mammalian protein complexes. Differential expression and gene ontology (GO) enrichment analyses suggested that 13 integrin complexes involved in cell adhesion were significantly downregulated in both low- and high-grade PCa compared with normal prostate, and that four Prothymosin alpha (ProTα) complexes were significantly upregulated in high-grade PCa compared with normal prostate. Moreover, differential co-regulation and GO enrichment analyses indicated that the assembly levels of six protein complexes involved in RNA splicing were significantly increased in low-grade PCa, and those of four subcomplexes of mitochondrial complex I were significantly increased in high-grade PCa, compared with normal prostate.

Conclusions

In summary, to the best of our knowledge, the study represents the first large-scale and quantitative, albeit indirect, comparison of individual protein complexes in human PCa tissue specimens. It may serve as a useful resource for better understanding the deregulation of protein complexes in primary PCa.

Electronic supplementary material

The online version of this article (10.1186/s12014-019-9236-2) contains supplementary material, which is available to authorized users.

Network Medicine in the Age of Biomedical Big Data

Network medicine is an emerging area of research dealing with molecular and genetic interactions, network biomarkers of disease, and therapeutic target discovery. Large-scale biomedical data generation offers a unique opportunity to assess the effect and impact of cellular heterogeneity and environmental perturbations on the observed phenotype. Marrying the two, network medicine with biomedical data provides a framework to build meaningful models and extract impactful results at a network level. In this review, we survey existing network types and biomedical data sources. More importantly, we delve into ways in which the network medicine approach, aided by phenotype-specific biomedical data, can be gainfully applied. We provide three paradigms, mainly dealing with three major biological network archetypes: protein-protein interaction, expression-based, and gene regulatory networks. For each of these paradigms, we discuss a broad overview of philosophies under which various network methods work. We also provide a few examples in each paradigm as a test case of its successful application. Finally, we delineate several opportunities and challenges in the field of network medicine. We hope this review provides a lexicon for researchers from biological sciences and network theory to come on the same page to work on research areas that require interdisciplinary expertise. Taken together, the understanding gained from combining biomedical data with networks can be useful for characterizing disease etiologies and identifying therapeutic targets, which, in turn, will lead to better preventive medicine with translational impact on personalized healthcare.

An integrative systems genetic analysis of mammalian lipid metabolism

Dysregulation of lipid homeostasis is a precipitating event in the pathogenesis and progression of hepatosteatosis and metabolic syndrome. These conditions are highly prevalent in developed societies and currently have limited options for diagnostic and therapeutic intervention. Here, using a proteomic and lipidomic-wide systems genetic approach, we interrogated lipid regulatory networks in 107 genetically distinct mouse strains to reveal key insights into the control and network structure of mammalian lipid metabolism. These include the identification of plasma lipid signatures that predict pathological lipid abundance in the liver of mice and humans, defining subcellular localization and functionality of lipid-related proteins, and revealing functional protein and genetic variants that are predicted to modulate lipid abundance. Trans-omic analyses using these datasets facilitated the identification and validation of PSMD9 as a previously unknown lipid regulatory protein. Collectively, our study serves as a rich resource for probing mammalian lipid metabolism and provides opportunities for the discovery of therapeutic agents and biomarkers in the setting of hepatic lipotoxicity.

CoExpresso: assess the quantitative behavior of protein complexes in human cells

Background

Translational and post-translational control mechanisms in the cell result in widely observable differences between measured gene transcription and protein abundances. Herein, protein complexes are among the most tightly controlled entities by selective degradation of their individual proteins. They furthermore act as control hubs that regulate highly important processes in the cell and exhibit a high functional diversity due to their ability to change their composition and their structure. Better understanding and prediction of these functional states demands methods for the characterization of complex composition, behavior, and abundance across multiple cell states. Mass spectrometry provides an unbiased approach to directly determine protein abundances across different cell populations and thus to profile a comprehensive abundance map of proteins.

Results

We provide a tool to investigate the behavior of protein subunits in known complexes by comparing their abundance profiles across up to 140 cell types available in ProteomicsDB. Thorough assessment of different randomization methods and statistical scoring algorithms allows determining the significance of concurrent profiles within a complex, therefore providing insights into the conservation of their composition across human cell types as well as the identification of intrinsic structures in complex behavior to determine which proteins orchestrate complex function. This analysis can be extended to investigate common profiles within arbitrary protein groups. CoExpresso can be accessed through http://computproteomics.bmb.sdu.dk/Apps/CoExpresso.

Conclusions

With the CoExpresso web service, we offer a potent scoring scheme to assess proteins for their co-regulation and thereby offer insight into their potential for forming functional groups like protein complexes.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2573-8) contains supplementary material, which is available to authorized users.

Synthetic lethal therapies for cancer: what's next after PARP inhibitors?

The genetic concept of synthetic lethality has now been validated clinically through the demonstrated efficacy of poly(ADP-ribose) polymerase (PARP) inhibitors for the treatment of cancers in individuals with germline loss-of-function mutations in either BRCA1 or BRCA2. Three different PARP inhibitors have now been approved for the treatment of patients with BRCA-mutant ovarian cancer and one for those with BRCA-mutant breast cancer; these agents have also shown promising results in patients with BRCA-mutant prostate cancer. Here, we describe a number of other synthetic lethal interactions that have been discovered in cancer. We discuss some of the underlying principles that might increase the likelihood of clinical efficacy and how new computational and experimental approaches are now facilitating the discovery and validation of synthetic lethal interactions. Finally, we make suggestions on possible future directions and challenges facing researchers in this field.

E-Cadherin/ROS1 Inhibitor Synthetic Lethality in Breast Cancer

The cell adhesion glycoprotein E-cadherin (CDH1) is commonly inactivated in breast tumors. Precision medicine approaches that exploit this characteristic are not available. Using perturbation screens in breast tumor cells with CRISPR/Cas9-engineered CDH1 mutations, we identified synthetic lethality between E-cadherin deficiency and inhibition of the tyrosine kinase ROS1. Data from large-scale genetic screens in molecularly diverse breast tumor cell lines established that the E-cadherin/ROS1 synthetic lethality was not only robust in the face of considerable molecular heterogeneity but was also elicited with clinical ROS1 inhibitors, including foretinib and crizotinib. ROS1 inhibitors induced mitotic abnormalities and multinucleation in E-cadherin-defective cells, phenotypes associated with a defect in cytokinesis and aberrant p120 catenin phosphorylation and localization. In vivo, ROS1 inhibitors produced profound antitumor effects in multiple models of E-cadherin-defective breast cancer. These data therefore provide the preclinical rationale for assessing ROS1 inhibitors, such as the licensed drug crizotinib, in appropriately stratified patients.Significance: E-cadherin defects are common in breast cancer but are currently not targeted with a precision medicine approach. Our preclinical data indicate that licensed ROS1 inhibitors, including crizotinib, should be repurposed to target E-cadherin-defective breast cancers, thus providing the rationale for the assessment of these agents in molecularly stratified phase II clinical trials. Cancer Discov; 8(4); 498-515. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 371.