Quantitative analysis of protein interaction network dynamics in yeast

Principles of Computation by Competitive Protein Dimerization Networks

Many biological signaling pathways employ proteins that competitively dimerize in diverse combinations. These dimerization networks can perform biochemical computations, in which the concentrations of monomers (inputs) determine the concentrations of dimers (outputs). Despite their prevalence, little is known about the range of input-output computations that dimerization networks can perform (their "expressivity") and how it depends on network size and connectivity. Using a systematic computational approach, we demonstrate that even small dimerization networks (3-6 monomers) are expressive, performing diverse multi-input computations. Further, dimerization networks are versatile, performing different computations when their protein components are expressed at different levels, such as in different cell types. Remarkably, individual networks with random interaction affinities, when large enough (≥8 proteins), can perform nearly all (~90%) potential one-input network computations merely by tuning their monomer expression levels. Thus, even the simple process of competitive dimerization provides a powerful architecture for multi-input, cell-type-specific signal processing.

Dynamic remodeling of Escherichia coli interactome in response to environmental perturbations

Proteins play an essential role in the vital biological processes governing cellular functions. Most proteins function as members of macromolecular machines, with the network of interacting proteins revealing the molecular mechanisms driving the formation of these complexes. Profiling the physiology-driven remodeling of these interactions within different contexts constitutes a crucial component to achieving a comprehensive systems-level understanding of interactome dynamics. Here, we apply co-fractionation mass spectrometry and computational modeling to quantify and profile the interactions of ∼2000 proteins in the bacterium Escherichia coli cultured under 10 distinct culture conditions. The resulting quantitative co-elution patterns revealed large-scale condition-dependent interaction remodeling among protein complexes involved in diverse biochemical pathways in response to the unique environmental challenges. The network-level analysis highlighted interactome-wide biophysical properties and structural patterns governing interaction remodeling. Our results provide evidence of the local and global plasticity of the E. coli interactome along with a rigorous generalizable framework to define protein interaction specificity. We provide an accompanying interactive web application to facilitate the exploration of these rewired networks.

Fit-Seq2.0: An Improved Software for High-Throughput Fitness Measurements Using Pooled Competition Assays

The fitness of a genotype is defined as its lifetime reproductive success, with fitness itself being a composite trait likely dependent on many underlying phenotypes. Measuring fitness is important for understanding how alteration of different cellular components affects a cell's ability to reproduce. Here, we describe an improved approach, implemented in Python, for estimating fitness in high throughput via pooled competition assays.

Network analysis of genome-wide association studies for drug target prioritisation

Over the past decades, genome-wide association studies (GWAS) have led to a dramatic expansion of genetic variants implicated with human traits and diseases. These advances are expected to result in new drug targets but the identification of causal genes and the cell biology underlying human diseases from GWAS remains challenging. Here, we review protein interaction network-based methods to analyse GWAS data. These approaches can rank candidate drug targets at GWAS-associated loci or among interactors of disease genes without direct genetic support. These methods identify the cell biology affected in common across diseases, offering opportunities for drug repurposing, as well as be combined with expression data to identify focal tissues and cell types. Going forward, we expect that these methods will further improve from advances in the characterisation of context specific interaction networks and the joint analysis of rare and common genetic signals.

OUP accepted manuscript

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

An atlas of protein-protein interactions across mouse tissues

Cellular processes arise from the dynamic organization of proteins in networks of physical interactions. Mapping the interactome has therefore been a central objective of high-throughput biology. However, the dynamics of protein interactions across physiological contexts remain poorly understood. Here, we develop a quantitative proteomic approach combining protein correlation profiling with stable isotope labeling of mammals (PCP-SILAM) to map the interactomes of seven mouse tissues. The resulting maps provide a proteome-scale survey of interactome rewiring across mammalian tissues, revealing more than 125,000 unique interactions at a quality comparable to the highest-quality human screens. We identify systematic suppression of cross-talk between the evolutionarily ancient housekeeping interactome and younger, tissue-specific modules. Rewired proteins are tightly regulated by multiple cellular mechanisms and are implicated in disease. Our study opens up new avenues to uncover regulatory mechanisms that shape in vivo interactome responses to physiological and pathophysiological stimuli in mammalian systems.

Protein-fragment complementation assays for large-scale analysis of protein-protein interactions

Protein-protein interactions (PPIs) orchestrate nearly all biological processes. They are also considered attractive drug targets for treating many human diseases, including cancers and neurodegenerative disorders. Protein-fragment complementation assays (PCAs) provide a direct and straightforward way to study PPIs in living cells or multicellular organisms. Importantly, PCAs can be used to detect the interaction of proteins expressed at endogenous levels in their native cellular environment. In this review, we present the principle of PCAs and discuss some of their advantages and limitations. We describe their application in large-scale experiments to investigate PPI networks and to screen or profile PPI targeting compounds.

A New Sequential Forward Feature Selection (SFFS) Algorithm for Mining Best Topological and Biological Features to Predict Protein Complexes from Protein-Protein Interaction Networks (PPINs)

Protein-protein interaction plays an important role in the understanding of biological processes in the body. A network of dynamic protein complexes within a cell that regulates most biological processes is known as a protein-protein interaction network (PPIN). Complex prediction from PPINs is a challenging task. Most of the previous computation approaches mine cliques, stars, linear and hybrid structures as complexes from PPINs by considering topological features and fewer of them focus on important biological information contained within protein amino acid sequence. In this study, we have computed a wide variety of topological features and integrate them with biological features computed from protein amino acid sequence such as bag of words, physicochemical and spectral domain features. We propose a new Sequential Forward Feature Selection (SFFS) algorithm, i.e., random forest-based Boruta feature selection for selecting the best features from computed large feature set. Decision tree, linear discriminant analysis and gradient boosting classifiers are used as learners. We have conducted experiments by considering two reference protein complex datasets of yeast, i.e., CYC2008 and MIPS. Human and mouse complex information is taken from CORUM 3.0 dataset. Protein interaction information is extracted from the database of interacting proteins (DIP). Our proposed SFFS, i.e., random forest-based Brouta feature selection in combination with decision trees, linear discriminant analysis and Gradient Boosting Classifiers outperforms other state of art algorithms by achieving precision, recall and F-measure rates, i.e. 94.58%, 94.92% and 94.45% for MIPS, 96.31%, 93.55% and 96.02% for CYC2008, 98.84%, 98.00%, 98.87 % for CORUM humans and 96.60%, 96.70%, 96.32% for CORUM mouse dataset complexes, respectively.

Contextualized Protein-Protein Interactions

Protein-protein interaction (PPI) databases are an important bioinformatics resource, yet existing literature-curated databases usually represent cell-type-agnostic interactions, which is at variance with our understanding that protein dynamics are context specific and highly dependent on their environment. Here, we provide a resource derived through data mining to infer disease- and tissue-relevant interactions by annotating existing PPI databases with cell-contextual information extracted from reporting studies. This resource is applicable to the reconstruction and analysis of disease-centric molecular interaction networks. We have made the data and method publicly available and plan to release scheduled updates in the future. We expect these resources to be of interest to a wide audience of researchers in the life sciences.

•

We present PPI Context: contextualization of existing literature-curated PPIs

•

A resource for filtering PPIs by cell-line information mined from reporting studies

•

A fast and flexible pipeline implementing the presented data mining method

Existing literature-curated protein-protein interaction (PPI) databases usually aggregate cell-type-agnostic interactions, yet PPIs are dependent on environmental conditions. Thus, new methods and resources for inferring the context in which a PPI is reported will extend their application and use in disease-centric modeling. We expect the resource presented in this article to be of high interest to those querying known interactions of proteins of interest, reconstruction and analyses of molecular interaction networks, and multi-omics data integration approaches.

A large accessory protein interactome is rewired across environments

To characterize how protein-protein interaction (PPI) networks change, we quantified the relative PPI abundance of 1.6 million protein pairs in the yeast Saccharomyces cerevisiae across nine growth conditions, with replication, for a total of 44 million measurements. Our multi-condition screen identified 13,764 pairwise PPIs, a threefold increase over PPIs identified in one condition. A few 'immutable' PPIs are present across all conditions, while most 'mutable' PPIs are rarely observed. Immutable PPIs aggregate into highly connected 'core' network modules, with most network remodeling occurring within a loosely connected 'accessory' module. Mutable PPIs are less likely to co-express, co-localize, and be explained by simple mass action kinetics, and more likely to contain proteins with intrinsically disordered regions, implying that environment-dependent association and binding is critical to cellular adaptation. Our results show that protein interactomes are larger than previously thought and contain highly dynamic regions that reorganize to drive or respond to cellular changes.

Genome-Wide Dynamic Evaluation of the UV-Induced DNA Damage Response

Genetic screens in Saccharomyces cerevisiae have allowed for the identification of many genes as sensors or effectors of DNA damage, typically by comparing the fitness of genetic mutants in the presence or absence of DNA-damaging treatments. However, these static screens overlook the dynamic nature of DNA damage response pathways, missing time-dependent or transient effects. Here, we examine gene dependencies in the dynamic response to ultraviolet radiation-induced DNA damage by integrating ultra-high-density arrays of 6144 diploid gene deletion mutants with high-frequency time-lapse imaging. We identify 494 ultraviolet radiation response genes which, in addition to recovering molecular pathways and protein complexes previously annotated to DNA damage repair, include components of the CCR4-NOT complex, tRNA wobble modification, autophagy, and, most unexpectedly, 153 nuclear-encoded mitochondrial genes. Notably, mitochondria-deficient strains present time-dependent insensitivity to ultraviolet radiation, posing impaired mitochondrial function as a protective factor in the ultraviolet radiation response.

Molecular networks in Network Medicine: Development and applications

Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.

Mass-spectrometry-based draft of the Arabidopsis proteome

Plants are essential for life and are extremely diverse organisms with unique molecular capabilities¹. Here we present a quantitative atlas of the transcriptomes, proteomes and phosphoproteomes of 30 tissues of the model plant Arabidopsis thaliana. Our analysis provides initial answers to how many genes exist as proteins (more than 18,000), where they are expressed, in which approximate quantities (a dynamic range of more than six orders of magnitude) and to what extent they are phosphorylated (over 43,000 sites). We present examples of how the data may be used, such as to discover proteins that are translated from short open-reading frames, to uncover sequence motifs that are involved in the regulation of protein production, and to identify tissue-specific protein complexes or phosphorylation-mediated signalling events. Interactive access to this resource for the plant community is provided by the ProteomicsDB and ATHENA databases, which include powerful bioinformatics tools to explore and characterize Arabidopsis proteins, their modifications and interactions.

Next-generation Interactomics: Considerations for the Use of Co-elution to Measure Protein Interaction Networks

Understanding how proteins interact is crucial to understanding cellular processes. Among the available interactome mapping methods, co-elution stands out as a method that is simultaneous in nature and capable of identifying interactions between all the proteins detected in a sample. The general workflow in co-elution methods involves the mild extraction of protein complexes and their separation into several fractions, across which proteins bound together in the same complex will show similar co-elution profiles when analyzed appropriately. In this review we discuss the different separation, quantification and bioinformatic strategies used in co-elution studies, and the important considerations in designing these studies. The benefits of co-elution versus other methods makes it a valuable starting point when asking questions that involve the perturbation of the interactome.

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.

A comprehensive review and evaluation of computational methods for identifying protein complexes from protein–protein interaction networks

Protein complexes are the fundamental units for many cellular processes. Identifying protein complexes accurately is critical for understanding the functions and organizations of cells. With the increment of genome-scale protein–protein interaction (PPI) data for different species, various computational methods focus on identifying protein complexes from PPI networks. In this article, we give a comprehensive and updated review on the state-of-the-art computational methods in the field of protein complex identification, especially focusing on the newly developed approaches. The computational methods are organized into three categories, including cluster-quality-based methods, node-affinity-based methods and ensemble clustering methods. Furthermore, the advantages and disadvantages of different methods are discussed, and then, the performance of 17 state-of-the-art methods is evaluated on two widely used benchmark data sets. Finally, the bottleneck problems and their potential solutions in this important field are discussed.

Neglog: Homology-Based Negative Data Sampling Method for Genome-Scale Reconstruction of Human Protein-Protein Interaction Networks

Rapid reconstruction of genome-scale protein-protein interaction (PPI) networks is instrumental in understanding the cellular processes and disease pathogenesis and drug reactions. However, lack of experimentally verified negative data (i.e., pairs of proteins that do not interact) is still a major issue that needs to be properly addressed in computational modeling. In this study, we take advantage of the very limited experimentally verified negative data from Negatome to infer more negative data for computational modeling. We assume that the paralogs or orthologs of two non-interacting proteins also do not interact with high probability. We coin an assumption as "Neglog" this assumption is to some extent supported by paralogous/orthologous structure conservation. To reduce the risk of bias toward the negative data from Negatome, we combine Neglog with less biased random sampling according to a certain ratio to construct training data. L₂-regularized logistic regression is used as the base classifier to counteract noise and train on a large dataset. Computational results show that the proposed Neglog method outperforms pure random sampling method with sound biological interpretability. In addition, we find that independent test on negative data is indispensable for bias control, which is usually neglected by existing studies. Lastly, we use the Neglog method to validate the PPIs in STRING, which are supported by gene ontology (GO) enrichment analyses.

Transcription Is Just the Beginning of Gene Expression Regulation: The Functional Significance of RNA-Binding Proteins to Post-transcriptional Processes in Plants

Plants have developed sophisticated mechanisms to compensate and respond to ever-changing environmental conditions. Research focus in this area has recently shifted towards understanding the post-transcriptional mechanisms that contribute to RNA transcript maturation, abundance and function as key regulatory steps in allowing plants to properly react and adapt to these never-ending shifts in their environments. At the center of these regulatory mechanisms are RNA-binding proteins (RBPs), the functional mediators of all post-transcriptional processes. In plants, RBPs are becoming increasingly appreciated as the critical modulators of core cellular processes during development and in response to environmental stimuli. With the majority of research on RBPs and their functions historically in prokaryotic and mammalian systems, it has more recently been unveiled that plants have expanded families of conserved and novel RBPs compared with their eukaryotic counterparts. To better understand the scope of RBPs in plants, we present past and current literature detailing specific roles of RBPs during stress response, development and other fundamental transition periods. In this review, we highlight examples of complex regulation coordinated by RBPs with a focus on the diverse mechanisms of plant RBPs and the unique processes they regulate. Additionally, we discuss the importance for additional research into understanding global interactions of RBPs on a systems and network-scale, with genome mining and annotation providing valuable insight for potential uses in improving crop plants in order to maintain high-level production in this era of global climate change.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

Analysis of networks of host proteins in the early time points following HIV transduction

Background

Utilization of quantitative proteomics data on the network level is still a challenge in proteomics data analysis. Currently existing models use sophisticated, sometimes hard to implement analysis techniques. Our aim was to generate a relatively simple strategy for quantitative proteomics data analysis in order to utilize as much of the data generated in a proteomics experiment as possible.

Results

In this study, we applied label-free proteomics, and generated a network model utilizing both qualitative, and quantitative data, in order to examine the early host response to Human Immunodeficiency Virus type 1 (HIV-1). A weighted network model was generated based on the amount of proteins measured by mass spectrometry, and analysis of weighted networks and functional sub-networks revealed upregulation of proteins involved in translation, transcription, and DNA condensation in the early phase of the viral life-cycle.

Conclusion

A relatively simple strategy for network analysis was created and applied to examine the effect of HIV-1 on host cellular proteome. We believe that our model may prove beneficial in creating algorithms, allowing for both quantitative and qualitative studies of proteome change in various biological and pathological processes by quantitative mass spectrometry.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2990-3) contains supplementary material, which is available to authorized users.

Co-complex protein membership evaluation using Maximum Entropy on GO ontology and InterPro annotation

Motivation

Protein–protein interactions (PPI) play a crucial role in our understanding of protein function and biological processes. The standardization and recording of experimental findings is increasingly stored in ontologies, with the Gene Ontology (GO) being one of the most successful projects. Several PPI evaluation algorithms have been based on the application of probabilistic frameworks or machine learning algorithms to GO properties. Here, we introduce a new training set design and machine learning based approach that combines dependent heterogeneous protein annotations from the entire ontology to evaluate putative co-complex protein interactions determined by empirical studies.

Results

PPI annotations are built combinatorically using corresponding GO terms and InterPro annotation. We use a S.cerevisiae high-confidence complex dataset as a positive training set. A series of classifiers based on Maximum Entropy and support vector machines (SVMs), each with a composite counterpart algorithm, are trained on a series of training sets. These achieve a high performance area under the ROC curve of ≤0.97, outperforming go2ppi—a previously established prediction tool for protein-protein interactions (PPI) based on Gene Ontology (GO) annotations.

Availability and implementation

https://github.com/ima23/maxent-ppi

Supplementary information

Supplementary data are available at Bioinformatics online.

Differential analysis of combinatorial protein complexes with CompleXChange

BACKGROUND

Although a considerable number of proteins operate as multiprotein complexes and not on their own, organism-wide studies so far are only able to quantify individual proteins or protein-coding genes in a condition-specific manner for a sizeable number of samples, but not their assemblies. Consequently, there exist large amounts of transcriptomic data and an increasing amount of data on proteome abundance, but quantitative knowledge on complexomes is missing. This deficiency impedes the applicability of the powerful tool of differential analysis in the realm of macromolecular complexes. Here, we present a pipeline for differential analysis of protein complexes based on predicted or manually assigned complexes and inferred complex abundances, which can be easily applied on a whole-genome scale.

RESULTS

We observed for simulated data that results obtained by our complex abundance estimation algorithm were in better agreement with the ground truth and physicochemically more reasonable compared to previous efforts that used linear programming while running in a fraction of the time. The practical usability of the method was assessed in the context of transcription factor complexes in human monocyte and lymphoblastoid samples. We demonstrated that our new method is robust against false-positive detection and reports deregulated complexomes that can only be partially explained by differential analysis of individual protein-coding genes. Furthermore we showed that deregulated complexes identified by the tool potentially harbor significant yet unused information content.

CONCLUSIONS

CompleXChange allows to analyze deregulation of the protein complexome on a whole-genome scale by integrating a plethora of input data that is already available. A platform-independent Java binary, a user guide with example data and the source code are freely available at https://sourceforge.net/projects/complexchange/ .

iSeq 2.0: A Modular and Interchangeable Toolkit for Interaction Screening in Yeast

We developed a flexible toolkit for combinatorial screening in Saccharomyces cerevisiae, which generates large libraries of cells, each uniquely barcoded to mark a combination of DNA elements. This interaction sequencing platform (iSeq 2.0) includes genomic landing pads that assemble combinations through sequential integration of plasmids or yeast mating, 15 barcoded plasmid libraries containing split selectable markers (URA3_AI, KanMX_AI, HphMX_AI, and NatMX_AI), and an array of ∼24,000 "double-barcoder" strains that can make existing yeast libraries iSeq compatible. Various DNA elements are compatible with iSeq: DNA introduced on integrating plasmids, engineered genomic modifications, or entire genetic backgrounds. DNA element libraries are modular and interchangeable, and any two libraries can be combined, making iSeq capable of performing many new combinatorial screens by short-read sequencing.

Genome-wide analysis of yeast expression data based on a priori generated co-regulation cliques

DNA microarrays are highly sensitive tools to evaluate the gene expression status of organismic samples and standardized array formats exist for many different sample types. Differential expression studies usually utilize the strongest upor downregulated genes to generate networks visualizing the relationships among these genes. To include all yeast genes in one analysis and to get broader information on all cellular responses, we test a priori input of predefined genome-wide expression cliques and subsequent statistical analysis of the expression data. To this end, we generate a set of 72 co-regulation cliques using the information from 3196 microarray experiments. The obtained cliques performed highly significant in gene ontology and transcription factor enrichment analyses. We then tested the clique set on individual microarray experiments reporting on responses to pheromone, glycerol versus glucose based growth and the cellular response to heat. In all cases a highly significant determination of affected expression cliques was possible based on their average expression differences, the positions of their genes within hit rankings (UpRegScore) or the enrichment of the Top200 hits in certain cliques. The 72 cliques were finally used to compare experiments, which reported on the transcriptional response to polyglutamine proteins of different lengths. Using the predefined clique set it is possible to identify with high sensitivity and good significance sample and condition specific changes to gene expression. We thus conclude that an analysis, starting with these 72 preformed expression cliques, can complement traditional microarray analyses by visualizing the entire response on a static genome-wide gene set.

Unbiased Fitness Estimation of Pooled Barcode or Amplicon Sequencing Studies

Standard practice for phenotyping complex cell pools is to measure the fold enrichment of genotype-specific amplicons after a period of competitive growth. Here, we show that fold-enrichment measures cannot be compared across genotype pools with different fitness distributions. We develop a method to calculate an unbiased estimate of relative fitness by tracking abundances over several time points and show how to optimize experimental protocols to minimize fitness measurement error.

Context-specific interactions in literature-curated protein interaction databases

BACKGROUND

Databases of literature-curated protein-protein interactions (PPIs) are often used to interpret high-throughput interactome mapping studies and estimate error rates. These databases combine interactions across thousands of published studies and experimental techniques. Because the tendency for two proteins to interact depends on the local conditions, this heterogeneity of conditions means that only a subset of database PPIs are interacting during any given experiment. A typical use of these databases as gold standards in interactome mapping projects, however, assumes that PPIs included in the database are indeed interacting under the experimental conditions of the study.

RESULTS

Using raw data from 20 co-fractionation experiments and six published interactomes, we demonstrate that this assumption is often false, with up to 55% of purported gold standard interactions showing no evidence of interaction, on average. We identify a subset of CORUM database complexes that do show consistent evidence of interaction in co-fractionation studies, and we use this subset as gold standards to dramatically improve interactome mapping as judged by the number of predicted interactions at a given error rate.

CONCLUSIONS

We recommend using this CORUM subset as the gold standard set in future co-fractionation studies. More generally, we recommend using the subset of literature-curated PPIs that are specific to the experimental context whenever possible.

Doxorubicin induces an extensive transcriptional and metabolic rewiring in yeast cells

Doxorubicin is one of the most effective chemotherapy drugs used against solid tumors in the treatment of several cancer types. Two different mechanisms, (i) intercalation of doxorubicin into DNA and inhibition of topoisomerase II leading to changes in chromatin structure, (ii) generation of free radicals and oxidative damage to biomolecules, have been proposed to explain the mode of action of this drug in cancer cells. A genome-wide integrative systems biology approach used in the present study to investigate the long-term effect of doxorubicin in Saccharomyces cerevisiae cells indicated the up-regulation of genes involved in response to oxidative stress as well as in Rad53 checkpoint sensing and signaling pathway. Modular analysis of the active sub-network has also revealed the induction of the genes significantly associated with nucleosome assembly/disassembly and DNA repair in response to doxorubicin. Furthermore, an extensive re-wiring of the metabolism was observed. In addition to glycolysis, and sulfate assimilation, several pathways related to ribosome biogenesis/translation, amino acid biosynthesis, nucleotide biosynthesis, de novo IMP biosynthesis and one-carbon metabolism were significantly repressed. Pentose phosphate pathway, MAPK signaling pathway biological processes associated with meiosis and sporulation were found to be induced in response to long-term exposure to doxorubicin in yeast cells.

Quantitative analysis of protein interaction network dynamics in yeast

Many cellular functions are mediated by protein–protein interaction networks, which are environment dependent. However, systematic measurement of interactions in diverse environments is required to better understand the relative importance of different mechanisms underlying network dynamics. To investigate environment‐dependent protein complex dynamics, we used a DNA‐barcode‐based multiplexed protein interaction assay in Saccharomyces cerevisiae to measure in vivo abundance of 1,379 binary protein complexes under 14 environments. Many binary complexes (55%) were environment dependent, especially those involving transmembrane transporters. We observed many concerted changes around highly connected proteins, and overall network dynamics suggested that “concerted” protein‐centered changes are prevalent. Under a diauxic shift in carbon source from glucose to ethanol, a mass‐action‐based model using relative mRNA levels explained an estimated 47% of the observed variance in binary complex abundance and predicted the direction of concerted binary complex changes with 88% accuracy. Thus, we provide a resource of yeast protein interaction measurements across diverse environments and illustrate the value of this resource in revealing mechanisms of network dynamics.