Tunnelling nanotube formation is driven by Eps8/IRSp53-dependent linear actin polymerization

Tunnelling nanotubes (TNTs) connect distant cells and mediate cargo transfer for intercellular communication in physiological and pathological contexts. How cells generate these actin-mediated protrusions to span lengths beyond those attainable by canonical filopodia remains unknown. Through a combination of micropatterning, microscopy, and optical tweezer-based approaches, we demonstrate that TNTs formed through the outward extension of actin achieve distances greater than the mean length of filopodia and that branched Arp2/3-dependent pathways attenuate the extent to which actin polymerizes in nanotubes, thus limiting their occurrence. Proteomic analysis using epidermal growth factor receptor kinase substrate 8 (Eps8) as a positive effector of TNTs showed that, upon Arp2/3 inhibition, proteins enhancing filament turnover and depolymerization were reduced and Eps8 instead exhibited heightened interactions with the inverted Bin/Amphiphysin/Rvs (I-BAR) domain protein IRSp53 that provides a direct connection with linear actin polymerases. Our data reveals how common protrusion players (Eps8 and IRSp53) form tunnelling nanotubes, and that when competing pathways overutilizing such proteins and monomeric actin in Arp2/3 networks are inhibited, processes promoting linear actin growth dominate to favour tunnelling nanotube formation.

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

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

Analysis of context-specific KRAS-effector (sub)complexes in Caco-2 cells

Ras is a key switch controlling cell behavior. In the GTP-bound form, Ras interacts with numerous effectors in a mutually exclusive manner, where individual Ras-effectors are likely part of larger cellular (sub)complexes. The molecular details of these (sub)complexes and their alteration in specific contexts are not understood. Focusing on KRAS, we performed affinity purification (AP)-mass spectrometry (MS) experiments of exogenously expressed FLAG-KRAS WT and three oncogenic mutants ("genetic contexts") in the human Caco-2 cell line, each exposed to 11 different culture media ("culture contexts") that mimic conditions relevant in the colon and colorectal cancer. We identified four effectors present in complex with KRAS in all genetic and growth contexts ("context-general effectors"). Seven effectors are found in KRAS complexes in only some contexts ("context-specific effectors"). Analyzing all interactors in complex with KRAS per condition, we find that the culture contexts had a larger impact on interaction rewiring than genetic contexts. We investigated how changes in the interactome impact functional outcomes and created a Shiny app for interactive visualization. We validated some of the functional differences in metabolism and proliferation. Finally, we used networks to evaluate how KRAS-effectors are involved in the modulation of functions by random walk analyses of effector-mediated (sub)complexes. Altogether, our work shows the impact of environmental contexts on network rewiring, which provides insights into tissue-specific signaling mechanisms. This may also explain why KRAS oncogenic mutants may be causing cancer only in specific tissues despite KRAS being expressed in most cells and tissues.

Targeted Quantification of Protein Phosphorylation and Its Contributions towards Mathematical Modeling of Signaling Pathways

Post-translational modifications (PTMs) are key regulatory mechanisms that can control protein function. Of these, phosphorylation is the most common and widely studied. Because of its importance in regulating cell signaling, precise and accurate measurements of protein phosphorylation across wide dynamic ranges are crucial to understanding how signaling pathways function. Although immunological assays are commonly used to detect phosphoproteins, their lack of sensitivity, specificity, and selectivity often make them unreliable for quantitative measurements of complex biological samples. Recent advances in Mass Spectrometry (MS)-based targeted proteomics have made it a more useful approach than immunoassays for studying the dynamics of protein phosphorylation. Selected reaction monitoring (SRM)-also known as multiple reaction monitoring (MRM)-and parallel reaction monitoring (PRM) can quantify relative and absolute abundances of protein phosphorylation in multiplexed fashions targeting specific pathways. In addition, the refinement of these tools by enrichment and fractionation strategies has improved measurement of phosphorylation of low-abundance proteins. The quantitative data generated are particularly useful for building and parameterizing mathematical models of complex phospho-signaling pathways. Potentially, these models can provide a framework for linking analytical measurements of clinical samples to better diagnosis and treatment of disease.

Disentangling ERBB Signaling in Breast Cancer Subtypes-A Model-Based Analysis

Simple Summary

Breast cancer subtypes are characterized by the expression and activity of estrogen-, progesterone- and HER2-receptors and differ by the treatment as well as patient prognosis. Tumors of the HER2-subtype overexpress this receptor and are successfully targeted with anti-HER2 therapies. We wanted to know if the HER2-receptor and the downstream signaling network act similarly also in the other subtypes and if this network could potentially be a therapeutic target beyond the HER2-positive subtype. To this end, we quantitatively assessed the wiring of signaling events in the individual subtypes to unravel the characteristics of HER-signaling. Our data along with a model-based analysis suggest that major parts of the intracellular signal transduction network are unchanged between the different breast cancer subtypes and that the clinical differences mostly come from the different levels at which these receptors are present in tumor cells as well as from the particular mutations that are present in individual tumors.

Abstract

Targeted therapies have shown striking success in the treatment of cancer over the last years. However, their specific effects on an individual tumor appear to be varying and difficult to predict. Using an integrative modeling approach that combines mechanistic and regression modeling, we gained insights into the response mechanisms of breast cancer cells due to different ligand–drug combinations. The multi-pathway model, capturing ERBB receptor signaling as well as downstream MAPK and PI3K pathways was calibrated on time-resolved data of the luminal breast cancer cell lines MCF7 and T47D across an array of four ligands and five drugs. The same model was then successfully applied to triple negative and HER2-positive breast cancer cell lines, requiring adjustments mostly for the respective receptor compositions within these cell lines. The additional relevance of cell-line-specific mutations in the MAPK and PI3K pathway components was identified via L₁ regularization, where the impact of these mutations on pathway activation was uncovered. Finally, we predicted and experimentally validated the proliferation response of cells to drug co-treatments. We developed a unified mathematical model that can describe the ERBB receptor and downstream signaling in response to therapeutic drugs targeting this clinically relevant signaling network in cell line that represent three major subtypes of breast cancer. Our data and model suggest that alterations in this network could render anti-HER therapies relevant beyond the HER2-positive subtype.

Integrating protein copy numbers with interaction networks to quantify stoichiometry in clathrin-mediated endocytosis

Proteins that drive processes like clathrin-mediated endocytosis (CME) are expressed at copy numbers within a cell and across cell types varying from hundreds (e.g. auxilin) to millions (e.g. clathrin). These variations contain important information about function, but without integration with the interaction network, they cannot capture how supply and demand for each protein depends on binding to shared and distinct partners. Here we construct the interface-resolved network of 82 proteins involved in CME and establish a metric, a stoichiometric balance ratio (SBR), that quantifies whether each protein in the network has an abundance that is sub- or super-stoichiometric dependent on the global competition for binding. We find that highly abundant proteins (like clathrin) are super-stoichiometric, but that not all super-stoichiometric proteins are highly abundant, across three cell populations (HeLa, fibroblast, and neuronal synaptosomes). Most strikingly, within all cells there is significant competition to bind shared sites on clathrin and the central AP-2 adaptor by other adaptor proteins, resulting in most being in excess supply. Our network and systematic analysis, including response to perturbations of network components, show how competition for shared binding sites results in functionally similar proteins having widely varying stoichiometries, due to variations in both abundance and their unique network of binding partners.

Reconstruction and analysis of a large-scale binary Ras-effector signaling network

BACKGROUND

Ras is a key cellular signaling hub that controls numerous cell fates via multiple downstream effector pathways. While pathways downstream of effectors such as Raf, PI3K and RalGDS are extensively described in the literature, how other effectors signal downstream of Ras is often still enigmatic.

METHODS

A comprehensive and unbiased Ras-effector network was reconstructed downstream of 43 effector proteins (converging onto 12 effector classes) using public pathway and protein-protein interaction (PPI) databases. The output is an oriented graph of pairwise interactions defining a 3-layer signaling network downstream of Ras. The 2290 proteins comprising the network were studied for their implication in signaling crosstalk and feedbacks, their subcellular localizations, and their cellular functions.

RESULTS

The final Ras-effector network consists of 2290 proteins that are connected via 19,080 binary PPIs, increasingly distributed across the downstream layers, with 441 PPIs in layer 1, 1660 in layer 2, and 16,979 in layer 3. We identified a high level of crosstalk among proteins of the 12 effector classes. A class-specific Ras sub-network was generated in CellDesigner (.xml file) and a functional enrichment analysis thereof shows that 58% of the processes have previously been associated to a respective effector pathway, with the remaining providing insights into novel and unexplored functions of specific effector pathways.

CONCLUSIONS

Our large-scale and cell general Ras-effector network is a crucial steppingstone towards defining the network boundaries. It constitutes a 'reference interactome' and can be contextualized for specific conditions, e.g. different cell types or biopsy material obtained from cancer patients. Further, it can serve as a basis for elucidating systems properties, such as input-output relationships, crosstalk, and pathway redundancy. Video Abstract.

Not all RAS mutations are equal: A detailed review of the functional diversity of RAS hot spot mutations

The RAS family of small GTPases are among the most frequently mutated oncogenes in human cancer. Approximately 20% of cancers harbor a RAS mutation, and >150 different missense mutations have been detected. Many of these mutations have mutant-specific biochemical defects that alter nucleotide binding and hydrolysis, effector interactions and cell signaling, prompting renewed efforts in the development of anti-RAS therapies, including the mutation-specific strategies. Previously viewed as undruggable, the recent FDA approval of a KRAS^G12C-selective inhibitor has offered real promise to the development of allele-specific RAS therapies. A broader understanding of the mutational consequences on RAS function must be developed to exploit additional allele-specific vulnerabilities. Approximately 94% of RAS mutations occur at one of three mutational "hot spots" at Gly¹², Gly¹³ and Gln⁶¹. Further, the single-nucleotide substitutions represent >99% of these mutations. Within this scope, we discuss the mutational frequencies of RAS isoforms in cancer, mutant-specific effector interactions and biochemical properties. By limiting our analysis to this mutational subset, we simplify the analysis while only excluding a small percentage of total mutations. Combined, these data suggest that the presence or absence of select RAS mutations in human cancers can be linked to their biochemical properties. Continuing to examine the biochemical differences in each RAS-mutant protein will continue to provide additional breakthroughs in allele-specific therapeutic strategies.

Predicted 'wiring landscape' of Ras-effector interactions in 29 human tissues

Ras is a plasma membrane (PM)-associated signaling hub protein that interacts with its partners (effectors) in a mutually exclusive fashion. We have shown earlier that competition for binding and hence the occurrence of specific binding events at a hub protein can modulate the activation of downstream pathways. Here, using a mechanistic modeling approach that incorporates high-quality proteomic data of Ras and 56 effectors in 29 (healthy) human tissues, we quantified the amount of individual Ras-effector complexes, and characterized the (stationary) Ras "wiring landscape" specific to each tissue. We identified nine effectors that are in significant amount in complex with Ras in at least one of the 29 tissues. We simulated both mutant- and stimulus-induced network re-configurations, and assessed their divergence from the reference scenario, specifically discussing a case study for two stimuli in three epithelial tissues. These analyses pointed to 32 effectors that are in significant amount in complex with Ras only if they are additionally recruited to the PM, e.g. via membrane-binding domains or domains binding to activated receptors at the PM. Altogether, our data emphasize the importance of tissue context for binding events at the Ras signaling hub.

The Ins and Outs of RAS Effector Complexes

RAS oncogenes are among the most commonly mutated proteins in human cancers. They regulate a wide range of effector pathways that control cell proliferation, survival, differentiation, migration and metabolic status. Including aberrations in these pathways, RAS-dependent signaling is altered in more than half of human cancers. Targeting mutant RAS proteins and their downstream oncogenic signaling pathways has been elusive. However, recent results comprising detailed molecular studies, large scale omics studies and computational modeling have painted a new and more comprehensive portrait of RAS signaling that helps us to understand the intricacies of RAS, how its physiological and pathophysiological functions are regulated, and how we can target them. Here, we review these efforts particularly trying to relate the detailed mechanistic studies with global functional studies. We highlight the importance of computational modeling and data integration to derive an actionable understanding of RAS signaling that will allow us to design new mechanism-based therapies for RAS mutated cancers.

Systems biology approaches to macromolecules: the role of dynamic protein assemblies in information processing

Macromolecular protein assemblies govern many cellular processes and are disturbed in many diseases including cancer. Often seen as static molecular machines, protein complexes involved in signal transduction networks exhibit intricate dynamics that are critical for their function. Using the RAS-RAF-MEK-ERK pathway as example we discuss recent progress in our understanding of protein complex dynamics achieved through mathematical modelling, computational simulations and structural studies. The emerging picture highlights that both spatial and temporal dynamics cooperate to enable correct signal processing and the fine tuning of timing, duration and strengths of signalling. These dynamic processes are subverted by oncogenic mutations and contribute to tumorigenesis and drug resistance.

Analysis of Ras-effector interaction competition in large intestine and colorectal cancer context

Cancer is the second leading cause of death globally, and colorectal cancer (CRC) is among the five most common cancers. The small GTPase KRAS is an oncogene that is mutated in ~30% of all CRCs. Pharmacological treatments of CRC are currently unsatisfactory, but much hope rests on network-centric approaches to drug development and cancer treatment. These approaches, however, require a better understanding of how networks downstream of Ras oncoproteins are connected in a particular tissue context – here colon and CRC. Previously we have shown that competition for binding to a ‘hub’ protein, such as Ras, can induce a rewiring of signal transduction networks. In this study, we analysed 56 established and predicted effectors that contain a structural domain with the potential ability to bind to Ras oncoproteins and their link to pathways coordinating intestinal homoeostasis and barrier function. Using protein concentrations in colon tissue and Ras-effector binding affinities, a computational network model was generated that predicted how effectors differentially and competitively bind to Ras in colon context. The model also predicted both qualitative and quantitative changes in Ras-effector complex formations with increased levels of active Ras – to simulate its upregulation in cancer – simply as an emergent property of competition for the same binding interface on the surface of Ras. We also considered how the number of Ras-effector complexes at the membrane can be increased by additional domains present in some effectors that are recruited to the membrane in response to specific conditions (inputs/stimuli/growth factors) in colon context and CRC.

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.

The Causes and Consequences of Genetic Interactions (Epistasis)

The same mutation can have different effects in different individuals. One important reason for this is that the outcome of a mutation can depend on the genetic context in which it occurs. This dependency is known as epistasis. In recent years, there has been a concerted effort to quantify the extent of pairwise and higher-order genetic interactions between mutations through deep mutagenesis of proteins and RNAs. This research has revealed two major components of epistasis: nonspecific genetic interactions caused by nonlinearities in genotype-to-phenotype maps, and specific interactions between particular mutations. Here, we provide an overview of our current understanding of the mechanisms causing epistasis at the molecular level, the consequences of genetic interactions for evolution and genetic prediction, and the applications of epistasis for understanding biology and determining macromolecular structures.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Modeling and simulating networks of interdependent protein interactions

Protein interactions are fundamental building blocks of biochemical reaction systems underlying cellular functions. The complexity and functionality of these systems emerge not only from the protein interactions themselves but also from the dependencies between these interactions, as generated by allosteric effects or mutual exclusion due to steric hindrance. Therefore, formal models for integrating and utilizing information about interaction dependencies are of high interest. Here, we describe an approach for endowing protein networks with interaction dependencies using propositional logic, thereby obtaining constrained protein interaction networks ("constrained networks"). The construction of these networks is based on public interaction databases as well as text-mined information about interaction dependencies. We present an efficient data structure and algorithm to simulate protein complex formation in constrained networks. The efficiency of the model allows fast simulation and facilitates the analysis of many proteins in large networks. In addition, this approach enables the simulation of perturbation effects, such as knockout of single or multiple proteins and changes of protein concentrations. We illustrate how our model can be used to analyze a constrained human adhesome protein network, which is responsible for the formation of diverse and dynamic cell-matrix adhesion sites. By comparing protein complex formation under known interaction dependencies versus without dependencies, we investigate how these dependencies shape the resulting repertoire of protein complexes. Furthermore, our model enables investigating how the interplay of network topology with interaction dependencies influences the propagation of perturbation effects across a large biochemical system. Our simulation software CPINSim (for Constrained Protein Interaction Network Simulator) is available under the MIT license at http://github.com/BiancaStoecker/cpinsim and as a Bioconda package (https://bioconda.github.io).

Stoichiometric balance of protein copy numbers is measurable and functionally significant in a protein-protein interaction network for yeast endocytosis

Stoichiometric balance, or dosage balance, implies that proteins that are subunits of obligate complexes (e.g. the ribosome) should have copy numbers expressed to match their stoichiometry in that complex. Establishing balance (or imbalance) is an important tool for inferring subunit function and assembly bottlenecks. We show here that these correlations in protein copy numbers can extend beyond complex subunits to larger protein-protein interactions networks (PPIN) involving a range of reversible binding interactions. We develop a simple method for quantifying balance in any interface-resolved PPINs based on network structure and experimentally observed protein copy numbers. By analyzing such a network for the clathrin-mediated endocytosis (CME) system in yeast, we found that the real protein copy numbers were significantly more balanced in relation to their binding partners compared to randomly sampled sets of yeast copy numbers. The observed balance is not perfect, highlighting both under and overexpressed proteins. We evaluate the potential cost and benefits of imbalance using two criteria. First, a potential cost to imbalance is that ‘leftover’ proteins without remaining functional partners are free to misinteract. We systematically quantify how this misinteraction cost is most dangerous for strong-binding protein interactions and for network topologies observed in biological PPINs. Second, a more direct consequence of imbalance is that the formation of specific functional complexes depends on relative copy numbers. We therefore construct simple kinetic models of two sub-networks in the CME network to assess multi-protein assembly of the ARP2/3 complex and a minimal, nine-protein clathrin-coated vesicle forming module. We find that the observed, imperfectly balanced copy numbers are less effective than balanced copy numbers in producing fast and complete multi-protein assemblies. However, we speculate that strategic imbalance in the vesicle forming module allows cells to tune where endocytosis occurs, providing sensitive control over cargo uptake via clathrin-coated vesicles.

Protein copy numbers are often found to be stoichiometrically balanced for subunits of multi-protein complexes. Imbalance is believed to be deleterious because it lowers complex yield (the dosage balance hypothesis) and increases the risk of misinteractions, but imbalance may also provide unexplored functional benefits. We show here that the benefits of stoichiometric balance can extend to larger networks of interacting proteins. We develop a method to quantify to what degree protein networks are balanced, and apply it to two networks. We find that the clathrin-mediated endocytosis system in yeast is statistically balanced, but not perfectly so, and explore the consequences of imbalance in the form of misinteractions and endocytic function. We also show that biological networks are more robust to misinteractions than random networks when balanced, but are more sensitive to misinteractions under imbalance. This suggests evolutionary pressure for proteins to be balanced and that any conserved imbalance should occur for functional reasons. We explore one such reason in the form of bottlenecking the endocytosis process. Our method can be generalized to other networks and used to identify out-of-balance proteins. Our results provide insight into how network design, expression level regulation, and cell fitness are intertwined.

Interaction Dynamics Determine Signaling and Output Pathway Responses

The understanding of interaction dynamics in signaling pathways can shed light on pathway architecture and provide insights into targets for intervention. Here, we explored the relevance of kinetic rate constants of a key upstream osmosensor in the yeast high-osmolarity glycerol-mitogen-activated protein kinase (HOG-MAPK) pathway to signaling output responses. We created mutant pairs of the Sln1-Ypd1 complex interface that caused major compensating changes in the association (k_on) and dissociation (k_off) rate constants (kinetic perturbations) but only moderate changes in the overall complex affinity (K_d). Yeast cells carrying a Sln1-Ypd1 mutant pair with moderate increases in k_on and k_off displayed a lower threshold of HOG pathway activation than wild-type cells. Mutants with higher k_on and k_off rates gave rise to higher basal signaling and gene expression but impaired osmoadaptation. Thus, the k_on and k_off rates of the components in the Sln1 osmosensor determine proper signaling dynamics and osmoadaptation.

Interactions between RASA2, CADM1, HIF1AN gene polymorphisms and body fatness with breast cancer: a population-based case-control study in China

Genome-wide association studies (GWAS) have indicated that gene polymorphisms in alleles of RAS p21 protein activator 2 (RASA2), cell adhesion molecule 1 (CADM1) and hypoxia inducible factor 1 alpha subunit inhibitor (HIF1AN) are associated with the risk of obesity. In this study, we explored the interactions between candidate SNPs of RASA2 (rs16851483), CADM1 (rs12286929) and HIF1AN (rs17094222) and body fatness for breast cancer risk. Unconditional logistic regression models were applied to measure the associations of related factors with breast cancer by odds ratios (ORs) and 95% confidence intervals (CIs). It was observed that cases had a statistically higher body mass index (BMI ≥ 28 kg/m², OR = 1.77), waist circumference (WC ≥ 90cm, OR = 2.89) and waist-to-hip ratio (WHR ≥ 0.9, OR = 3.41) as compared to controls. Significant differences were also found in the genotype distributions of RASA2 rs16851483 T/T homozygote and CADM1 rs12286929 G/A heterozygote between cases and controls, with an OR of 1.68 (95% CI: 1.10-2.56) and 0.80 (95% CI: 0.64-0.99), respectively. Furthermore, significant interactions were observed between polymorphisms of three genes and body fatness for the risk of breast cancer based on both additive and multiplicative scales. These results of our study suggest that body fatness possibly plays an important role in the development of breast cancer and this risk might be modified by specific genotypes of some potential genes, especially for obese women in China.

Protein-protein binding selectivity and network topology constrain global and local properties of interface binding networks

Protein-protein interactions networks (PPINs) are known to share a highly conserved structure across all organisms. What is poorly understood, however, is the structure of the child interface interaction networks (IINs), which map the binding sites proteins use for each interaction. In this study we analyze four independently constructed IINs from yeast and humans and find a conserved structure of these networks with a unique topology distinct from the parent PPIN. Using an IIN sampling algorithm and a fitness function trained on the manually curated PPINs, we show that IIN topology can be mostly explained as a balance between limits on interface diversity and a need for physico-chemical binding complementarity. This complementarity must be optimized both for functional interactions and against mis-interactions, and this selectivity is encoded in the IIN motifs. To test whether the parent PPIN shapes IINs, we compared optimal IINs in biological PPINs versus random PPINs. We found that the hubs in biological networks allow for selective binding with minimal interfaces, suggesting that binding specificity is an additional pressure for a scale-free-like PPIN. We confirm through phylogenetic analysis that hub interfaces are strongly conserved and rewiring of interactions between proteins involved in endocytosis preserves interface binding selectivity.

Cell signaling as a cognitive process

Cellular identity as defined through morphology and function emerges from intracellular signaling networks that communicate between cells. Based on recursive interactions within and among these intracellular networks, dynamical solutions in terms of biochemical behavior are generated that can differ from those in isolated cells. In this way, cellular heterogeneity in tissues can be established, implying that cell identity is not intrinsically predetermined by the genetic code but is rather dynamically maintained in a cognitive manner. We address how to experimentally measure the flow of information in intracellular biochemical networks and demonstrate that even simple causality motifs can give rise to rich, context-dependent dynamic behavior. The concept how intercellular communication can result in novel dynamical solutions is applied to provide a contextual perspective on cell differentiation and tumorigenesis.

Chromosome-wise Protein Interaction Patterns and Their Impact on Functional Implications of Large-Scale Genomic Aberrations

Gene copy-number changes influence phenotypes through gene-dosage alteration and subsequent changes of protein complex stoichiometry. Human trisomies where gene copy numbers are increased uniformly over entire chromosomes provide generic cases for studying these relationships. In most trisomies, gene and protein level alterations have fatal consequences. We used genome-wide protein-protein interaction data to identify chromosome-specific patterns of protein interactions. We found that some chromosomes encode proteins that interact infrequently with each other, chromosome 21 in particular. We combined the protein interaction data with transcriptome data from human brain tissue to investigate how this pattern of global interactions may affect cellular function. We identified highly connected proteins that also had coordinated gene expression. These proteins were associated with important neurological functions affecting the characteristic phenotypes for Down syndrome and have previously been validated in mouse knockout experiments. Our approach is general and applicable to other gene-dosage changes, such as arm-level amplifications in cancer.

Protein abundance of AKT and ERK pathway components governs cell type-specific regulation of proliferation

Signaling through the AKT and ERK pathways controls cell proliferation. However, the integrated regulation of this multistep process, involving signal processing, cell growth and cell cycle progression, is poorly understood. Here, we study different hematopoietic cell types, in which AKT and ERK signaling is triggered by erythropoietin (Epo). Although these cell types share the molecular network topology for pro‐proliferative Epo signaling, they exhibit distinct proliferative responses. Iterating quantitative experiments and mathematical modeling, we identify two molecular sources for cell type‐specific proliferation. First, cell type‐specific protein abundance patterns cause differential signal flow along the AKT and ERK pathways. Second, downstream regulators of both pathways have differential effects on proliferation, suggesting that protein synthesis is rate‐limiting for faster cycling cells while slower cell cycles are controlled at the G1‐S progression. The integrated mathematical model of Epo‐driven proliferation explains cell type‐specific effects of targeted AKT and ERK inhibitors and faithfully predicts, based on the protein abundance, anti‐proliferative effects of inhibitors in primary human erythroid progenitor cells. Our findings suggest that the effectiveness of targeted cancer therapy might become predictable from protein abundance.

A Comprehensive View of the β-Arrestinome

G protein-coupled receptors (GPCRs) are membrane receptors critically involved in sensing the environment and orchestrating physiological processes. As such, they transduce extracellular signals such as hormone, neurotransmitters, ions, and light into an integrated cell response. The intracellular trafficking, internalization, and signaling ability of ligand-activated GPCRs are controlled by arrestins, adaptor proteins that they interact with upon ligand binding. β-arrestins 1 and 2 in particular are now considered as hub proteins assembling multiprotein complexes to regulate receptor fate and transduce diversified cell responses. While more than 400 β-arrestin interaction partners have been identified so far, much remains to be learnt on how discrimination between so many binding partners is accomplished. Here, we gathered the interacting partners of β-arrestins through database mining and manual curation of the literature to map the β-arrestin interactome (β-arrestinome). We discussed several parameters that determine compatible (AND) or mutually exclusive (XOR) binding of β-arrestin interactors, such as structural constraints, intracellular abundance, or binding affinity.

HIPPIE v2.0: enhancing meaningfulness and reliability of protein-protein interaction networks

The increasing number of experimentally detected interactions between proteins makes it difficult for researchers to extract the interactions relevant for specific biological processes or diseases. This makes it necessary to accompany the large-scale detection of protein–protein interactions (PPIs) with strategies and tools to generate meaningful PPI subnetworks. To this end, we generated the Human Integrated Protein–Protein Interaction rEference or HIPPIE (http://cbdm.uni-mainz.de/hippie/). HIPPIE is a one-stop resource for the generation and interpretation of PPI networks relevant to a specific research question. We provide means to generate highly reliable, context-specific PPI networks and to make sense out of them. We just released the second major update of HIPPIE, implementing various new features. HIPPIE grew substantially over the last years and now contains more than 270 000 confidence scored and annotated PPIs. We integrated different types of experimental information for the confidence scoring and the construction of context-specific networks. We implemented basic graph algorithms that highlight important proteins and interactions. HIPPIE's graphical interface implements several ways for wet lab and computational scientists alike to access the PPI data.

Network biology concepts in complex disease comorbidities

The co-occurrence of diseases can inform the underlying network biology of shared and multifunctional genes and pathways. In addition, comorbidities help to elucidate the effects of external exposures, such as diet, lifestyle and patient care. With worldwide health transaction data now often being collected electronically, disease co-occurrences are starting to be quantitatively characterized. Linking network dynamics to the real-life, non-ideal patient in whom diseases co-occur and interact provides a valuable basis for generating hypotheses on molecular disease mechanisms, and provides knowledge that can facilitate drug repurposing and the development of targeted therapeutic strategies.

Protein-protein interactions generate hidden feedback and feed-forward loops to trigger bistable switches, oscillations and biphasic dose-responses

Protein-protein interactions (PPIs) defined as reversible association of two proteins to form a complex, are undoubtedly among the most common interaction motifs featured in cells. Recent large-scale proteomic studies have revealed an enormously complex interactome of the cell, consisting of tens of thousands of PPIs with numerous signalling hubs. PPIs have functional roles in regulating a wide range of cellular processes including signal transduction and post-translational modifications, and de-regulation of PPIs is implicated in many diseases including cancers and neuro-degenerative disorders. Despite the ubiquitous appearance and physiological significance of PPIs, our understanding of the dynamic and functional consequences of these simple motifs remains incomplete, particularly when PPIs occur within large biochemical networks. We employ quantitative, dynamic modelling to computationally analyse salient dynamic features of the PPI motifs and PPI-containing signalling networks varying in topological architecture. Our analyses surprisingly reveal that simple reversible PPI motifs, when being embedded into signalling cascades, could give rise to extremely rich and complex regulatory dynamics in the absence of explicit positive and negative feedback loops. Our work represents a systematic investigation of the dynamic properties of PPIs in signalling networks, and the results shed light on how this simple event may potentiate diverse and intricate behaviours in vivo.

Tuneable endogenous mammalian target complementation via multiplexed plasmid-based recombineering

Understanding the quantitative functional consequences of human disease mutations requires silencing of endogenous genes and expression of mutants at close to physiological levels. Changing protein levels above or below these levels is also important for system perturbation and modelling. Fast design optimization demands flexible interchangeable cassettes for endogenous gene silencing and tuneable expression. Here, we introduce ‘TEMTAC’, a multigene recombineering and delivery system for simultaneous siRNA-based knockdown and regulated mutant (or other variant) expression with different dynamic ranges. We show its applicability by confirming known phenotypic effects for selected mutations for BRAF, HRAS, and SHP2.

Advances in template-based protein docking by utilizing interfaces towards completing structural interactome

The increase in the number of structurally determined protein complexes strengthens template-based docking (TBD) methods for modelling protein-protein interactions (PPIs). These methods utilize the known structures of protein complexes as templates to predict the quaternary structure of the target proteins. The templates may be partial or complete structures. Interface based (partial) methods have recently gained interest due in part to the observation that the interface regions are reusable. We describe how available template interfaces can be used to obtain the structural models of protein interactions. Despite the agreement that a majority of the protein complexes can be modelled using the available Protein Data Bank (PDB) structures, a handful of studies argue that we need more template proteins to increase the structural coverage of PPIs. We also discuss the performance of the interface TBD methods at large scale, and the significance of capturing multiple conformations for improving accuracy.

Recurrent inactivating RASA2 mutations in melanoma

Analysis of 501 melanoma exomes identified RASA2, encoding a RasGAP, as a tumor-suppressor gene mutated in 5% of melanomas. Recurrent loss-of-function mutations in RASA2 were found to increase RAS activation, melanoma cell growth and migration. RASA2 expression was lost in ≥30% of human melanomas and was associated with reduced patient survival. These findings identify RASA2 inactivation as a melanoma driver and highlight the importance of RasGAPs in cancer.

The Architecture of the TIR Domain Signalosome in the Toll-like Receptor-4 Signaling Pathway

Activated Toll-like receptors (TLRs) cluster in lipid rafts and induce pro- and anti-tumor responses. The organization of the assembly is critical to the understanding of how these key receptors control major signaling pathways in the cell. Although several models for individual interactions were proposed, the entire TIR-domain signalosome architecture has not been worked out, possibly due to its complexity. We employ a powerful algorithm, crystal structures and experimental data to model the TLR4 and its cluster. The architecture that we obtain with 8 MyD88 molecules provides the structural basis for the MyD88-templated myddosome helical assembly and receptor clustering; it also provides clues to pro- and anti-inflammatory signaling pathways branching at the signalosome level to Mal/MyD88 and TRAM/TRIF pro- and anti-inflammatory pathways. The assembly of MyD88 death domain (DD) with TRAF3 (anti-viral/anti-inflammatory) and TRAF6 (pro-inflammatory) suggest that TRAF3/TRAF6 binding sites on MyD88 DD partially overlap, as do IRAK4 and FADD. Significantly, the organization illuminates mechanisms of oncogenic mutations, demonstrates that almost all TLR4 parallel pathways are competitive and clarifies decisions at pathway branching points. The architectures are compatible with the currently-available experimental data and provide compelling insights into signaling in cancer and inflammation pathways.

Unraveling the molecular architecture of a G protein-coupled receptor/β-arrestin/Erk module complex

β-arrestins serve as signaling scaffolds downstream of G protein-coupled receptors, and thus play a crucial role in a plethora of cellular processes. Although it is largely accepted that the ability of β-arrestins to interact simultaneously with many protein partners is key in G protein-independent signaling of GPCRs, only the precise knowledge of these multimeric arrangements will allow a full understanding of the dynamics of these interactions and their functional consequences. However, current experimental procedures for the determination of the three-dimensional structures of protein-protein complexes are not well adapted to analyze these short-lived, multi-component assemblies. We propose a model of the receptor/β-arrestin/Erk1 signaling module, which is consistent with most of the available experimental data. Moreover, for the β-arrestin/Raf1 and the β-arrestin/ERK interactions, we have used the model to design interfering peptides and shown that they compete with both partners, hereby demonstrating the validity of the predicted interaction regions.

Dissecting the calcium-induced differentiation of human primary keratinocytes stem cells by integrative and structural network analyses

The molecular details underlying the time-dependent assembly of protein complexes in cellular networks, such as those that occur during differentiation, are largely unexplored. Focusing on the calcium-induced differentiation of primary human keratinocytes as a model system for a major cellular reorganization process, we look at the expression of genes whose products are involved in manually-annotated protein complexes. Clustering analyses revealed only moderate co-expression of functionally related proteins during differentiation. However, when we looked at protein complexes, we found that the majority (55%) are composed of non-dynamic and dynamic gene products ('di-chromatic'), 19% are non-dynamic, and 26% only dynamic. Considering three-dimensional protein structures to predict steric interactions, we found that proteins encoded by dynamic genes frequently interact with a common non-dynamic protein in a mutually exclusive fashion. This suggests that during differentiation, complex assemblies may also change through variation in the abundance of proteins that compete for binding to common proteins as found in some cases for paralogous proteins. Considering the example of the TNF-α/NFκB signaling complex, we suggest that the same core complex can guide signals into diverse context-specific outputs by addition of time specific expressed subunits, while keeping other cellular functions constant. Thus, our analysis provides evidence that complex assembly with stable core components and competition could contribute to cell differentiation.

An analysis of critical factors for quantitative immunoblotting

Immunoblotting (also known as Western blotting) combined with digital image analysis can be a reliable method for analyzing the abundance of proteins and protein modifications, but not every immunoblot-analysis combination produces an accurate result. I illustrate how sample preparation, protocol implementation, detection scheme, and normalization approach profoundly affect the quantitative performance of immunoblotting. This study implemented diagnostic experiments that assess an immunoblot-analysis workflow for accuracy and precision. The results showed that ignoring such diagnostics can lead to pseudoquantitative immunoblot data that markedly overestimate or underestimate true differences in protein abundance.

Versatile communication strategies among tandem WW domain repeats

Interactions mediated by short linear motifs in proteins play major roles in regulation of cellular homeostasis since their transient nature allows for easy modulation. We are still far from a full understanding and appreciation of the complex regulation patterns that can be, and are, achieved by this type of interaction. The fact that many linear-motif-binding domains occur in tandem repeats in proteins indicates that their mutual communication is used extensively to obtain complex integration of information toward regulatory decisions. This review is an attempt to overview, and classify, different ways by which two and more tandem repeats cooperate in binding to their targets, in the well-characterized family of WW domains and their corresponding polyproline ligands.

Attenuation of the suppressive activity of cellular splicing factor SRSF3 by Kaposi sarcoma-associated herpesvirus ORF57 protein is required for RNA splicing

Kaposi sarcoma-associated herpesvirus (KSHV) ORF57 is a multifunctional post-transcriptional regulator essential for viral gene expression during KSHV lytic infection. ORF57 requires interactions with various cellular proteins for its function. Here, we identified serine/arginine-rich splicing factor 3 (SRSF3, formerly known as SRp20) as a cellular cofactor involved in ORF57-mediated splicing of KSHV K8β RNA. In the absence of ORF57, SRSF3 binds to a suboptimal K8β intron and inhibits K8β splicing. Knockdown of SRSF3 promotes K8β splicing, mimicking the effect of ORF57. The N-terminal half of ORF57 binds to the RNA recognition motif of SRSF3, which prevents SRSF3 from associating with the K8β intron RNA and therefore attenuates the suppressive effect of SRSF3 on K8β splicing. ORF57 also promotes splicing of heterologous non-KSHV transcripts that are negatively regulated by SRSF3, indicating that the effect of ORF57 on SRSF3 activity is independent of RNA target. SPEN proteins, previously identified as ORF57-interacting partners, suppress ORF57 splicing activity by displacing ORF57 from SRSF3-RNA complexes. In summary, we have identified modulation of SRSF3 activity as the molecular mechanism by which ORF57 promotes RNA splicing.

Visualizing and manipulating temporal signaling dynamics with fluorescence-based tools

The use of genome-wide proteomic and RNA interference approaches has moved our understanding of signal transduction from linear pathways to highly integrated networks centered on core nodes. However, probing the dynamics of flow of information through such networks remains technically challenging. In particular, how the temporal dynamics of an individual pathway can elicit distinct outcomes in a single cell type and how multiple pathways may interact sequentially or synchronously to influence cell fate remain open questions in many contexts. The development of fluorescence-based reporters and optogenetic regulators of pathway activity enables the analysis of signaling in living cells and organisms with unprecedented spatiotemporal resolution and holds the promise of addressing these key questions. We present a brief overview of the evidence for the importance of temporal dynamics in cellular regulation, introduce these fluorescence-based tools, and highlight specific studies that leveraged these tools to probe the dynamics of information flow through signaling networks. In particular, we highlight two studies in Caenorhabditis elegans sensory neurons and cultured mammalian cells that demonstrate the importance of signal dynamics in determining cellular responses.

Analyzing protein-protein interactions in the post-interactomic era. Are we ready for the endgame?

Mapping protein-protein interactions in genome-wide scales revealed thousands of novel binding partners in each of the explored model organisms. Organizing these hits in comprehensive ways is becoming increasingly important for systems biology approaches to understand complex cellular processes and diseases. However, proteome wide interaction techniques and their resulting global networks are not revealing the topologies of networks that are truly operating in the cell. In this short review I will discuss which prerequisites have to be fulfilled and which experimental methods might be practicable to translate primary protein interaction data into network presentations that help in understanding cellular processes.

2013: Signaling breakthroughs of the year

The editorial staff and distinguished scientists in the field of cell signaling nominated diverse research as advances for 2013. Breakthroughs in understanding how spatial and temporal signals control cellular behavior ranged from filtering high-frequency stimuli to interpreting circadian inputs. This year's nominations also highlight the importance of understanding cell signaling in the context of physiology and disease, such as links between the nervous system and cancer. Furthermore, the application of new techniques to study cell signaling--such as optogenetics, DNA editing with CRISPR-Cas9, and sequencing untranslated regions of transcripts--continues to expand the realm and impact of signaling research.