RIP1-dependent linear and nonlinear recruitments of caspase-8 and RIP3 respectively to necrosome specify distinct cell death outcomes

There remains a significant gap in our quantitative understanding of crosstalk between apoptosis and necroptosis pathways. By employing the SWATH-MS technique, we quantified absolute amounts of up to thousands of proteins in dynamic assembling/de-assembling of TNF signaling complexes. Combining SWATH-MS-based network modeling and experimental validation, we found that when RIP1 level is below ~1000 molecules/cell (mpc), the cell solely undergoes TRADD-dependent apoptosis. When RIP1 is above ~1000 mpc, pro-caspase-8 and RIP3 are recruited to necrosome respectively with linear and nonlinear dependence on RIP1 amount, which well explains the co-occurrence of apoptosis and necroptosis and the paradoxical observations that RIP1 is required for necroptosis but its increase down-regulates necroptosis. Higher amount of RIP1 (>~46,000 mpc) suppresses apoptosis, leading to necroptosis alone. The relation between RIP1 level and occurrence of necroptosis or total cell death is biphasic. Our study provides a resource for encoding the complexity of TNF signaling and a quantitative picture how distinct dynamic interplay among proteins function as basis sets in signaling complexes, enabling RIP1 to play diverse roles in governing cell fate decisions. Electronic supplementary material The online version of this article (10.1007/s13238-020-00810-x) contains supplementary material, which is available to authorized users.


INTRODUCTION
Tumor necrosis factor (TNF) can induce apoptosis or necroptosis depending on cellular contexts, and apoptosis and necroptosis pathways can compete with and convert into each other (Han et al., 2011;Brenner et al., 2015). TNF signaling is initiated by TNF-TNFR1 (TNF receptor 1) ligation that promotes the formation of complex I. Complex I mainly consists of TNFR1-associated death domain protein (TRADD), receptor-interacting protein kinase 1 (RIP1), TNF receptor-associated factor 2 (TRAF2), cellular inhibitor of apoptosis 1/2 (cIAP1/2) and some other proteins. Polyubiquitination of RIP1 occurs in complex I, leading to NF-κB activation. The non-ubiquitinated or de-ubiquitinated RIP1 recruits Fas-associated death domain protein (FADD) and caspase-8 to form complex II, which allows caspase-8 auto-processing, resulting in apoptosis (Brenner et al., 2015). If there is RIP3 expression in cells, RIP3 is recruited to the RIP1-FADD-caspase-8 complex to form necrosome. Autophosphorylation of RIP1 and RIP3 occurs in necrosome and phosphorylated RIP3 recruits mixed lineage kinase domain-like protein (MLKL) (Zhang et al., 2009;Cai et al., 2014), which is then phosphorylated by RIP3 and released from necrosome to execute necroptosis. Caspase-8 exhibits low level protease activity in necrosome which is responsible for the cleavage and inactivation of RIP1 and RIP3 to halt necroptotic process (Tummers et al., 2017;Newton et al., 2019a;Newton et al., 2019b).
Protein interactions are essential in determining signaling consequences. Quantification study of these interactions could clarify how death signals are integrated to direct a specific cell fate. Combination of modeling and experimental approaches has been used to reveal essential features and control circuits of various cellular signaling pathways (Ma et al., 2009;Nakakuki et al., 2010;Yang et al., 2018;Bashor et al., 2019), including NF-κB and apoptosis pathways (Albeck et al., 2008;Shinohara et al., 2014). However, achieving quantitative mechanistic insights is still a major challenge because of limited consistent and comprehensive data in experiments (Kitano, 2005;Cox and Mann, 2011). Recent advances in the quantification capability of proteomics are now able to overcome this limitation (Aebersold and Mann, 2016). To quantitatively analyze the dynamic formation of TNF-induced signaling complexes of cell death, we employed immune-purification coupled with proteomics. SWATH-MS (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry) technique (Ludwig et al., 2018) was used to obtain reproducible and absolute protein amounts across multiple samples, which are currently the most sensitive inputs for mathematical modeling. Using the SWATH-MS-based modeling and experimental analysis, our study provides a quantitative framework for understanding how RIP1 distinctly recruits RIP3 and procaspase-8 into necrosome, thereby determining cell death type through influencing necrosome assembling and deassembling, and also shows a biphasic relationship between RIP1 level and the occurrence of necroptosis or the total cell death.
c Figure 1. Application of IP-SWATH workflow in quantifying TNF signaling complexes. (A) Experimental scheme of IP-SWATH workflow. L929 cells were treated with Flag-TNF or TNF (10 ng/mL) for 5 to 10 timepoints, and TNFR1, RIP1 and RIP3 complexes were immunoprecipitated with anti-Flag agarose beads. Purified heavy isotope amino acid-labeled targeted proteins were spiked into the IP samples. IP samples were digested with trypsin and peptides were analyzed using DDA and SWATH-MS. Group-DIA was employed to analyze SWATH-MS data for generation of pseudo-spectra files. DDA files and pseudospectra files were subjected to database searches, followed by targeted analysis of SWATH-MS. (B) Dynamic range of SWATH-MS. Different amounts of purified recombinant human MLKL protein were spiked into murine MLKL knockout (KO) cell lysates in triplicates, followed by trypsin digestion. The peptide samples were analyzed using SWATH-MS. SWATH-MS data were analyzed using Peakview software, and the product ion peak area of the MLKLspecific peptide LGQLIYEQCEK was directly taken from the Peakview. The sum of top three peak areas of product ions represented the peptide amount. (C) Reproducibility analysis of quantitative protein intensities in TNFR1 dataset. The protein intensities in two samples, each contains three replicates, were compared in TNFR1 dataset, and Pearson correlations were calculated. (D) Heat map of quantitative proteins in TNFR1 dataset. The known TNFR1 complex components were indicated in the right panel. (E) Differential expression analysis of proteins at each timepoint in TNFR1 dataset. Proteins with |Log2(fold change)| > 2 and −Log10(P-value) > 2 were considered significantly changed. Up-regulated proteins were labeled in red, and down-regulated proteins were labeled in blue. (F) XICs of peptides of high-confidence interaction proteins in TNFR1 dataset. Traces in different colors represent different product ions of given peptides. (G) Dynamic range of all identified proteins estimated by molecules per cell (mpc) in TNFR1 dataset. Some representative interaction proteins of TNFR1 were highlighted.

Absolute quantification of proteins in TNF signaling complexes by SWATH-MS
We employed MS to systematically analyze the dynamic assembly of three critical complexes-TNFR1, RIP1 and RIP3 complexes of TNF signaling in L929 cell line, which is a well-established cellular model to study necroptosis. The responses of wildtype, Flag-RIP1, and Flag-RIP3 cell lines to TNF stimulation were similar in terms of activation of NF-κB, MAPKs, and MLKL (Fig. S1A). For absolute protein quantification, heavy amino acid-labeled proteins were spiked into purified IP complexes, which were subsequently subjected to DDA and SWATH-MS analyses (Fig. 1A). Targeted extraction of human MLKL peptides (Fig. 1B) and β-galactosidase (Fig. S1B) using Peakview software from SWATH-MS data revealed that excellent quantification accuracy with about 4 orders of dynamic range could be achieved (Table S1). We combined the internal library generated by Group-DIA  and the external library built from DDA for peptide detection ( Fig. S1C and Table S2). Collectively, 1,466, 1,297 and 2,355 proteins were quantified across 18, 30 and 15 IP samples in TNFR1, RIP1 and RIP3 datasets, respectively (Table S2). Protein abundance comparison showed that Pearson correlation coefficients were about 0.89-0.94, 0.81-0.94, and 0.87-0.96 between any two different runs for TNFR1, RIP1 and RIP3 datasets respectively (Figs. 1C and S1D). Coefficients of variation of log2-transformed protein abundance at all time-points were below 10% (Fig. S1E). Thus, the entire experiment is excellent in reproducibility.
Hierarchical clustering of TNFR1, RIP1 and RIP3 interaction proteins indicated a small number of proteins changed in abundance after TNF treatment (Figs. 1D and S1F), suggesting a high level of non-specific background proteins. We therefore distinguished specifically recruited proteins from the background proteins by analyzing the amount change of proteins in complexes at different time periods. Differential expression analysis of each individual protein at different time points was conducted and shown in Figs. 1E and S1G. Manual inspection of all XICs (extracted ion chromatograms) of significantly changed proteins in Figs. 1E and S1G indicated that 11, 9 and 3 proteins were recruited to TNFR1, RIP1 and RIP3 complexes, respectively (Table S2). The representative peptide XICs of RIP1, TRADD and A20 of the 11 proteins recruited to TNFR1 complex were shown in Fig. 1F. These identified proteins are components of the complexes (Brenner et al., 2015), confirming the accuracy of our identification by IP-SWATH.
We quantified absolute amounts of key proteins in TNFR1, RIP1 and RIP3 complexes by using the spiked-in purified heavy amino acids-labeled proteins. The amounts of proteins in the complexes were then estimated by the "TOP3 peptides" approach (Table S3) (Al Shweiki et al., 2017). Absolute amounts of proteins in these three complexes are shown in Figs. 1G and S1H. The significantly changed proteins in TNFR1 complex are highlighted (Fig. 1G). Collectively, our IP-SWATH workflow provides an efficient method to accurately identify and absolutely quantify the dynamically assembled proteins in TNF signaling complexes.

SWATH-MS-based modeling of TNF-induced cell death signaling network
The amounts of the 11, 9 and 3 proteins recruited to TNFR1, RIP1 and RIP3 complexes respectively ranged from ∼20 to 12,000 mpc as determined by our SWATH-MS (Table S3). The high-abundant proteins are likely to be signaling transducers whereas the others could be catalysts. We therefore simplified illustration of TNF-induced signaling pathway by using these abundant proteins as transducers ( Fig. 2A). Complex I is formed by TRADD and RIP1, competitively binding to TNFR1. Complex II is composed of RIP1, FADD and caspase-8, which has the same components as necrosome except for lacking RIP3. In necrosome, activated caspase-8 cleaves and inactivates RIP1 and RIP3 to inhibit necroptosis. Uncleaved RIP3 phosphorylates and recruits MLKL, eventually leading to necroptosis (Han et al., 2011;Berghe et al., 2014). Experimental data suggest that procaspase-8 can be either activated by TRADD (Mechanism 1) or RIP1 (Mechanism 2) independently or by TRADD and RIP1 together (Mechanism 3) ( Fig The circled letter P denotes phosphorylation state. Simplified schematic diagrams that represent three different regulation mechanisms of pro-caspase-8 (Pro-C8) by TRADD (Mechanism 1), RIP1 (Mechanism 2) or TRADD and RIP1 (Mechanism 3) are shown in the right panel. Mutual inhibition between TRADD and RIP1 is presented to describe the competitive relation between the two proteins binding TNFR1. (B) Kinetic scheme of Models 1, 2 and 3. Lines ended by symbols represent the chemical reactions characterized by Reactions V1-V44 in Table S4. Parameters for the individually numbered reactions are given in Table S6. Ø represents degraded proteins. Model 1 (Mechanism 1) includes the reactions described by black and orange lines, whereas Model 2 (Mechanism 2) involves the reactions described by black and blue lines. Model 3 (Mechanism 3) contains all reactions. (C) Simulation results using Model 1 (lines) and SWATH-MS data (dots) of the time-course responses upon TNF stimulation. Error bars denote standard deviation (SD) for three independent experiments. The unit is molecules per cell (mpc). (D) Deviation (R 2 ) between simulation results of the three models and SWATH-MS data.

RESEARCH ARTICLE
2014; Remijsen et al., 2014). Accordingly, we developed three corresponding mathematical models, Model 1 with Mechanism 1, Model 2 with Mechanism 2 and Model 3 with Mechanism 3, which incorporated protein interactions, phosphorylation, cleavage and enzymatic reactions (Fig. 2B). By comparing the dynamics of these models to experimental data, we can quantitatively explore the emergent properties of TNF signaling complexes. All the biochemical reactions in Fig. 2B (V1-V46) and reaction rates in models are elaborated in Table S4. These models are described by a cast of ordinary differential equations (ODEs), which are represented by compounds and kinetic parameters based on the law of mass action. Complete model descriptions are presented in the Supplementary Materials. To examine the reliability of the models, simulation curves were compared with the SWATH-MS data shown in Figs. 1G and S1H. Comparison results indicated that Models 1, 2 and 3 can well reproduce the amounts of major components in TNFR1, RIP1 and RIP3 complexes obtained by SWATH-MS (Figs. 2C and S2). Squares of the correlation coefficients (R 2 ) between SWATH-MS data and simulation results are all larger than 0.9 for these three models (Fig. 2D). Hence, all of these models exhibit high confidence for further clarifying the transduction mechanism how RIP1 is formed in complexes to exert strategic control, thereby governing cell fate.

Biphasic roles of RIP1 in RIP3-dependent necroptosis
We next applied these three models to quantitatively dissect the functional roles of RIP1 in necroptosis. RIP3 phosphorylation is used as the marker of necroptosis (Zhang et al., 2009). RIP1 level in wildtype cells is treated as 100% RIP1, which is ∼51,000 ± 3,100 mpc based on three independent SWATH-MS data. By simulation, the reduction of RIP3 phosphorylation is associated with a decrease of RIP1 level in both Model 1 and Model 3 ( To avoid bias, we also used random parameters instead of SWATH-MS data-optimized parameters in these models for analysis. As each of Models 1, 2 and 3 contains more than 50 parameters (Table S6), we reduced parameters to 14, 14 and 16 respectively for simplified Model 1 (SModel 1), Model 2 (SModel 2), and Model 3 (SModel 3) ( Fig. 3B and Supplementary Materials) to make the simulation feasible. Random parameter sets were generated through a Monte Carlo approach. For each SModel, we sampled 10 6 random parameter sets using Latin hypercube sampling method and then for each parameter set determined how RIP1 regulates RIP3 phosphorylation. The positive or negative regulation can be evaluated by the value of Δp-RIP3, which is the level difference of RIP3 phosphorylation between wildtype and the cell containing 50% RIP1 (Fig. 3B). The number of parameter sets that yield negative regulation of RIP3 phosphorylation by RIP1 was calculated for each SModel and the corresponding set number distributions are presented in Fig. 3C. The results indicated that zero parameter set occurs to generate negative regulation with SModel 1 (black line), while certain numbers of parameter sets are obtained with SModels 2 and 3. SModel 2 gives a much larger set number than SModel 3 (Fig. 3C, red and blue lines). Percentage of the summed set numbers in all sampled parameter sets which is referred to as probability, is shown in Fig. 3D. Thus, sampling of random parameters also substantiates that, rather than SModels 1 and 3, SModel 2 is more highly preferred to robustly achieve the negative regulation of RIP3 phosphorylation by RIP1. Gillespie algorithm (Gillespie, 1977) was also performed to evaluate the negative regulation of RIP3 phosphorylation by RIP1 and the result is similar (Supplementary Materials).
To make further comprehensive analysis, RIP3 phosphorylation was simulated under a serial of RIP1 levels (Fig. 3E). As expected, RIP1 down-expression reduces RIP3 phosphorylation in Models 1 and 3, but enhances RIP3 phosphorylation only in Model 2. Model 2 predicted that RIP3 phosphorylation reaches its maximum level as RIP1 decreases to an extremely low level (∼1,000 mpc, e.g., ∼2% of wildtype), and then rapidly declines when RIP1 reduces further, showing a biphasic relationship of RIP3 phosphorylation to RIP1 level.
To test our model prediction, we experimentally knocked down RIP1 to three different expression levels with RIP1specific short hairpin RNA (shRNA), and detected the increase of TNF-induced phosphorylation of RIP3 and MLKL (Fig. 3F). Similar results were obtained when another RIP1 shRNA was used (Fig. S3A). Deletion of RIP1 completely blocked TNF-induced RIP3 phosphorylation (Fig. 3G). The band intensities in Fig. 3F and 3G and a replicate (Fig. S3B) were compared with Model 2 prediction (Fig. 3H), quantitatively confirming the biphasic regulation of RIP3 phosphorylation by RIP1. Therefore, RIP1 functions as an activator in necroptosis within extremely low level range (<∼2% of wildtype) and as an inhibitor at higher level.

Linear and nonlinear recruitments of pro-caspase-8 and RIP3 to necrosome
To dissect the underlying mechanism how RIP1 biphasically regulates RIP3 phosphorylation, we employed Model 2 to examine the dynamic behaviors of RIP1, pro-caspase-8, and Ppm1b , which are known to modulate RIP3 phosphorylation in necrosome (Fig. 4A). When RIP1 level decreased, RIP3 phosphorylation increased, and the amounts of RIP1 and Ppm1b that interacted with RIP3 were barely influenced (Fig. 4B). However, decrease of RIP1 reduced pro-caspase-8 that was associated with RIP3. Since pro-caspase-8 cleaves RIP1/RIP3 in necrosome, we introduced a reduction of pro-caspase-8 level in Model 2 (Figs. 4C and S4). Simulation results indicated that RIP3

RESEARCH ARTICLE
Xiang Li et al. phosphorylation exhibited a gradual increase when procaspase-8 level decreased. However, the trend of RIP1 decrease-associated increase of RIP3 phosphorylation became less pronounced when pro-caspase-8 level decreased. Such increase of RIP3 phosphorylation no longer occurred when pro-caspase-8 expression reduced to zero. Therefore, simulation results suggested that pro-caspase-8 activity is required for the up-regulation of RIP3 phosphorylation by decreasing RIP1 level.
We utilized z-Val-Ala-DL-Asp-fluoromethylketone (zVAD), a pan-caspase inhibitor to experimentally validate the prediction. RIP1 knockdown led to an increase of RIP3 phosphorylation upon TNF treatment (Fig. 4D, top panel). In comparison, RIP1 knockdown barely affected RIP3 phosphorylation when zVAD was included (Fig. 4D, bottom panel). In addition, decrease of RIP1, which enhanced RIP3 phosphorylation in wildtype cells, did not influence RIP3 phosphorylation in caspase-8 KO cells (Fig. 4E). The band intensities of RIP1 level and phosphorylated RIP3 in Western blot were quantified and plotted in Fig. 4F, supporting the idea that pro-caspase-8 recruitment to necrosome is essential for RIP1 negatively regulating RIP3 phosphorylation.
The suppression of RIP3 phosphorylation by pro-caspase-8 requires RIP1 in necrosome (Tummers et al., 2017). To independently address the effects of RIP1 on pro-caspase-8 activity and RIP3 phosphorylation, we made an analysis by removing the inter-inhibition between pro-caspase-8 and RIP3 in Model 2 (Fig. 4G). Simulation with three different RIP1 expression levels suggested that when RIP1 level decreased, the amount of pro-caspase-8 binding to RIP1 was reduced, but RIP3 phosphorylation kept constant (Fig. 4H). Further comprehensive simulations predicted that decrease of RIP1 results in a progressive decrease of procaspase-8 recruitment to necrosome, whereas RIP3 phosphorylation is barely influenced by decreasing RIP1 level from 100% to 2% and then presents a steep reduction (Fig. 4I). Thus, the response of pro-caspase-8 to RIP1 level is mostly linear, whereas RIP3 phosphorylation is not, and presents an ultrasensitive threshold pattern (Fig. 4J).
We next validated the distinct dynamic responses shown in Fig. 4J. Flag-RIP3-expressing RIP3 KO cells pre-treated with RIP1 shRNA were used to examine the behaviors of pro-caspase-8 recruitment and RIP3 phosphorylation in the cells with different RIP1 levels (Fig. 4K). The results suggested that pro-caspase-8 amounts decreased linearly in RIP3 complex when RIP1 decreased while RIP3 phosphorylation almost kept unchanged. The quantified Western blot results shown in Fig. 4L are in good agreement with the predictions in Fig. 4I, confirming the distinct dependence of pro-caspase-8 recruitment to necrosome and RIP3 phosphorylation on RIP1 level.

RIP1-suppressed caspase-8 activation does not cleave RIP1/RIP3
Model 2 well explains RIP1's biphasic roles in necroptosis. However, its manifestation that deletion of RIP1 can completely abolish the involvement of pro-caspase-8 conflicts with the experimental result that TNF induces quick caspase-8 activation and apoptosis in RIP1 KO cells (Fig. 5A) (Vanlangenakker et al., 2011). It is known that in Complex I, TRADD and RIP1 compete for TNFR1 binding and mediate downstream caspase-8 activation (Zheng et al., 2006;Wang et al., 2008). We therefore refined Model 2 by including additional TRADD-dependent caspase-8 activation and named it as Model 4 (Mechanism 4) (Fig. 5B). After the parameter fitting, Model 4 can reproduce not only the MS data with an R 2 value equal to that of Model 2 ( Fig. 5C) but also RIP1's biphasic roles in necroptosis (Fig. 5D, left panel). Importantly, Model 4, but not Model 2, reveals that the activation level of caspase-8 in RIP1 KO cells is higher than that in RIP3 KO cells (right panel in Figs. 5D and S5A), which is consistent with previous studies (Vanlangenakker et al., 2011;Remijsen et al., 2014).
We next sought to dissect the function of RIP1 in apoptosis with the refined model. Simulation demonstrated that the progressive decrease of RIP1 leads to a gradual increase of caspase-8 activation (Fig. 5E). This prediction is supported by our experiment that knocking down of RIP1 promoted caspase-8 activation (Figs. 5F and S3B). Further analysis determined the induced caspase-8 activation was TRADD-dependent. When TRADD expression level was reduced in Model 4, it indeed resulted in progressive decrease of caspase-8 activation at varying RIP1 levels (Figs. 5G and S5B). Consistently, the RIP1 knockdown-induced increase of caspase-8 activation cannot be detected in TRADD KO cells (Fig. 5H) Unexpectedly, topology of Model 4 suggests that unlike RIP1, TRADD-activated caspase-8 neither cleaves RIP1/ RIP3 nor is suppressed by RIP3 (Fig. 5B). Thus, the caspase-8 activation can occur independently of RIP3 phosphorylation in RIP1 knockdown cells (Fig. S3B). RIP1 biphasically regulates necroptosis through the basis set of RIP1, RIP3 and pro-caspase-8, while RIP1 limits apoptosis through the basis set of RIP1, TRADD and pro-caspase-8. Finally, we examined whether Model 4 could precisely quantify RIP1 level-determined total cell death, including apoptosis and necroptosis. We experimentally measured the cell death in RIP1 knockdown and RIP1 KO cells. As shown in Fig. 5I, the predicted death rate curves matched well with the experimental data, exhibiting an n-shaped biphasic relationship between RIP1 expression level and cell death.

Distinct cell death outcomes determined by RIP1 level in TNF signaling complexes
Taken together, we can derive a quantitative picture of TNFinduced complexes formation dynamics, signaling transduction and cell fate decisions under different RIP1 expression levels (Fig. 6A). TNF treatment leads to the recruitment of RIP1 and TRADD to TNFR1 in a competitive manner ( Fig. S6A and S6B) (Zheng et al., 2006). When RIP1 level is extremely low, i.e., < ∼2% of wildtype cells (< ∼1000 mpc), TRADD-mediated caspase-8 activation is predominant, which solely leads to apoptosis (Fig. 6A, right panel). Figure 6B shows the recruitment of pro-caspase-8 to TRADD complex in Flag-TRADD-expressing TRADD KO cell line after RIP1 was knocked down. 48 h DOX treatment significantly reduced overall RIP1 level and RIP1 in TRADD complex, which led to pro-caspase-8 recruitment to TRADD (Fig. 6B). It is known that TRADD and RIP1 might be subsequently released either together or separately from TNFR1 (Fig. 6A, left and middle panel) (Zheng et al., 2006). When RIP1 level is within ∼2%-90% of wildtype (∼1000-46,000 mpc), TRADD-and/or TRADD+-RIP1-dependent apoptosis should occur. Meanwhile, a fraction of RIP1, which is not associated with TRADD, could recruit RIP3 to initiate necroptosis (Fig. 6A, middle panel). This signal branching is supported by the observation that TRADD and RIP3 are located in different complexes, and both complexes contain RIP1 (Fig. 6C). Simultaneous activation of caspase-8 and RIP3 (Fig. S3B) and the suppression of necroptosis by an increase of RIP1 level happen (Fig. 6A, middle panel). Since high RIP1 level eliminates TRADD-induced caspase-8 activation (Fig. 5G), TNF stimulation leads to necroptosis alone when RIP1 level is > ∼90% of wildtype (> ∼46,000 mpc) (Fig. 6A, left panel).
Whether cell undergoes apoptosis or necroptosis depends on which pathway is dominant. As there is no mutual inhibition between the downstream signaling events of apoptosis and necroptosis in RIP1 knockdown cells (Fig. 6A, middle panel), co-occurrence of apoptosis and necroptosis in an individual cell is possible. Immunofluorescence microscopy data supported such an idea (Fig. 6D).
Hence, RIP1 level is a "switch node" in determining the outcomes of cell death. Consistent with the proposed picture shown in Fig. 6A, simulation result indicated that death rate is low at both extremely low and high RIP1 levels, while high death rate is observed within a broad middle range of RIP1 level (Fig. 6E, left panel). We further plotted the contribution of apoptosis and necroptosis to cell death (Fig. 6E, right  panel), showing that the occurrence of apoptosis or necroptosis alone is typical at the range of 0%-2% (0-1000 mpc) or ∼90%-100% (∼46,000-51,000 mpc) RIP1 levels respectively. Simultaneous death outcome occurs at the range of ∼2%-90% RIP1 levels with about equal contribution of apoptosis and necroptosis to total cell death. The equal contribution of apoptosis and necroptosis in TNF-induced cell death is likely to be resulted from the similar increasing tendency of caspase-8 activation (Fig. S4B) and RIP3 phosphorylation (Fig. S5B) to RIP1 level decrease. Counting of apoptotic and necroptotic cells of the experiment shown in Fig. 6D validated the predicted equal contribution in RIP1 knockdown cells (Fig. S6C).

DISCUSSION
Despite the fact that TNF signaling pathway has been extensively studied, quantitative elucidation of the components, their interactions and their distinct dynamic responses in complexes is still lacking. Here we studied TNF-induced cell death from a quantitative aspect. We employed absolute quantification methods to obtain proteomic data of the dynamic assembling and de-assembling of TNF-induced  Xiang Li et al. TNFR1, RIP1, and RIP3 complexes, and then applied these data to mathematical modeling. The SWATH-MS-based model provides an absolutely quantitative picture of TNFinduced cell death signaling, which revealed previously unknown roles of RIP1 level in specifying cell fate decisions.

RESEARCH ARTICLE
As a key upstream regulator in determining cell fate, RIP1 level in cells is modulated by various factors. Chen et al. found that geldanamycin can decrease RIP1 level through decreasing the Hsp90 protein level in primary cortical neurons, which might play a role in the prevention of stroke (Chen et al., 2012). RIP1 expression level is also proved to be elevated by baculovirus transduction in mammalian cells (Wang et al., 2017). Besides, a recent study performed by Someda et al. suggested that retinoic acid enhances the expression of RIP1, which stimulates necroptosis (Someda et al., 2020). Deficiency of selenium, one of the essential trace elements in human and animal body, induces the increase of RIP1 expression and promotes tracheal injury through necroptosis (Wang et al., 2020). Therefore, RIP1 level can be modulated by cellular mechanisms, which could play a role in determining cell viability or death and the mode of cell death.
Sufficient RIP3 expression in a given cell originally had no or low RIP3 expression can switch it death from apoptosis to necroptosis (Zhang et al., 2009). Simulated result suggests that the effect of RIP3 level on death rate and mode of cell death is quite different from that of RIP1 level (Fig. S8A). Necroptosis occurs only when RIP3 level is high, and the death rate is first reduced with the decrease of RIP3 and then increases but the death mode converts to apoptosis (Fig. S8A, lower panel). Further analysis (Fig. S8B) indicates that decrease of RIP3 significantly reduces the contribution of necroptosis to cell death when RIP1 level decreases.
The roles of RIP1 in necroptosis are complicated as RIP1 is not only required for necroptosis but also has an inhibitory effect (Weinlich and Green, 2014). Our study demonstrated that the dual roles are RIP1 level-dependent. RIP1 promotes necroptosis within extremely low level range (< ∼2% of wildtype) and inhibits necroptosis at higher levels, exhibiting a biphasic regulation of necroptosis by RIP1. Further mechanistic insights unraveled that the biphasic relationship is determined by the nonlinear threshold pattern of RIP3 and linear pro-caspase-8 recruitments to necrosome. The identified threshold pattern indicates that full phosphorylation of RIP3 only requires a very low RIP1 level (∼ 2% of wildtype), which is supported by a recent structure study of the RIP1-RIP3 necrosome (Mompean et al., 2018). This pattern is likely because a small amount RIP1-RIP3 hetero-interaction is sufficient to trigger RIP3-RIP3 homo-interaction and the latter can form a higher-order assembly structure of active necrosome (Wu, 2013;Bashor et al., 2019). RIP1 binds to FADD to recruit pro-caspase-8 (Han et al., 2011;Berghe et al., 2014). Here we showed that RIP1 level is linearly associated with the recruitment of pro-caspase-8 to necrosome, which inhibits the signaling toward RIP3 phosphorylation. The result is well in line with previous literature (Oberst et al., 2011;Orozco et al., 2014), suggesting that the proposed mechanism of the triangular relationship among RIP1, RIP3 and caspase-8 in the determination of cell death outcomes might be generally applicable.
The roles of RIP1 in regulating cellular events are known to be regulated by RIP1 protein modifications such as ubiquitylation and phosphorylation (Ofengeim and Yuan, 2013). These modifications affect either RIP1 level or interaction with other proteins or RIP1 kinase activity in different complexes (Dondelinger et al., 2017;Jaco et al., 2017;Xu et al., 2018;Newton et al., 2019a;Newton et al., 2019b). We observed that phosphorylation of RIP1 (Meng et al., 2018) was barely affected when the expression level of RIP1 decreased in cells (Fig. S8C), and the ubiquitination level of RIP1 was significantly reduced with the decrease of RIP1 level in cells (Fig. S6A) but it barely affected NF-κB activation (Fig. S7A). Due to technical limitation, we are unable to dissect the contribution of scaffold and kinase activity of RIP1 in our current study. It is known that the promotion of necroptosis by RIP1 is its kinase activity dependent whereas the suppression is mediated by the scaffold function of RIP1 Kaiser et al., 2014). As both functions of RIP1 are closely associated with its abundance in cells, our analyses presented quantitative relations between necroptosis and the combined effect of RIP1 kinase activity and its scaffold function.
A "mutual inhibition" relationship between apoptosis and necroptosis pathways is generally assumed as caspase-8 and RIP3 can completely block the activation of each other (Han et al., 2011;Berghe et al., 2014). The role of caspase-8 in cleaving RIP1/RIP3 is well known. However, our refined b Figure 6. Distinct cell death outcomes determined by RIP1 level. (A) Proposed scheme of complexes assembly and cell death types that are quantitatively determined by RIP1 level. The thick and thin arrows denote strong and weak signal transduction processes, respectively. (B) Flag-TRADD-expressing TRADD KO cells were treated with TNF. The cell lysates were immunoprecipitated with anti-Flag antibodies. Immunocomplexes and cell lysates were analyzed by Western blot. (C) Flag-TRADD-expressing TRADD KO cells and Flag-RIP3expressing RIP3 KO cells were treated with TNF plus zVAD and then were analyzed by Western blot. (D) RIP1 was knocked down in DOX-inducible MLKL-Flag expressing MLKL KO cells. Cells were induced with DOX and were treated with or without TNF. Cells were then fixed and immunostained for Flag and cleaved-caspase-3, counterstained with hochest. Distributions of MLKL in cells are defined as four types: uniformly diffused with or without cleaved-caspase-3, plasma membrane location with or without cleaved-caspase-3. Each scale bar indicates 10 µm. (E) Two-dimensional diagram simulations of RIP1 decrease induced n-shaped cell death rate upon TNF treatment (left panel). Right panel is the contribution proportions of apoptosis (brown area) and necroptosis (green area) after 6 h TNF treatment. Xiang Li et al. model suggested that TRADD-dependent caspase-8 activation does not cleave RIP1/RIP3. Further experiments ascertained the reason that TRADD and RIP3 are located in two distinct complexes, which provides a regulatory mechanism of the co-occurrence of apoptosis and necroptosis in RIP1 knockdown cells. Simultaneous apoptosis and necroptosis was even observed in individual cells. Thus, the relationship between apoptosis and necroptosis pathways is more like "co-existing", rather than "mutual inhibition". This finding raises a question how cells make fate decisions in "co-existing" state. We propose that there should be a "speed competition" between apoptosis and necroptosis pathways and the cell fate is determined by the pathway which reaches the final destination first.

RESEARCH ARTICLE
Although the biphasic relationship between RIP1 level and cell death is demonstrated in vitro but its relevance to RIP1's function in vivo is unclear. Nonetheless, a few in vivo evidences support this finding. The effect of RIP1 knockdown in mouse was evaluated by Suda et al. Decrease of RIP1 markedly exacerbates liver injury through massive apoptosis but not necroptosis (Suda et al., 2016), which could be resulted from TRADD dependent apoptosis presented in our Model 4 (Fig. 5G). Kaiser et al. showed that RIP1-deficient animals are more sensitive to necroptosis and apoptosis in response to diverse stimuli (Kaiser et al., 2014), which could resemble the situation that decrease of RIP1 level increased caspase-8 activation and RIP3 phosphorylation in our Model 4 (Fig. 6E). We believe that due to cell type differences and tissue organ specificities the in vivo role of RIP1 should be very complicated, and the downregulation of RIP1 level added one more layer to the complicity. In conclusion, our study presented comprehensive multi-approach analyses for understanding the roles of RIP1 on cell death determination. We suggest that precise quantitative description of signaling transduction with SWATH-MS-based modeling can provide a deeper insight into signaling dynamics that encode biological decision making.
To generate RIP1 N-terminal Flag knock-in L929 cells, homologous recombination strategy was employed. The rAAV (Adeno-associated virus) targeting vector was constructed by insertion of left homologous arm (about 1kb genomic DNA sequences upstream of the start codon of ripk1) and right homologous arm (about 1 kb genomic DNA sequences downstream of the star codon of ripk1) into the rAAV-Neo-Lox P-flag KI vector. Targeting rAAV virues were packaged in 293T cells. L929 cells were infected with the targeting rAAV virus and then selected for neomycin-resistant clones 24 h post infection. Those clones were then screened for homologous recombination by genomic PCR and the positive clones were infected with adenovirus expressing Cre-recombinase to excise the neomycin gene cassette. The final successful flag knock-in clones were confirmed by genomic PCR.
Purification of heavy amino acid labeled proteins 293T cells were cultured in amino acid deficient DMEM (Thermo) supplemented with "heavy" 13C6, 15N2 L-Lysine and 13C6, 15N4 L-Arginine (Sigma), 10% dialyzed FBS (Thermo), 0.1 mg/mL streptomycin, and 0.2 U/mL penicillin. Cells were grown for six doubling times prior to use in order for the cells to be fully incorporated with the labeled amino acids. Since overexpression of RIP3 in 293T causes extensive cell death, we synthesize the DNA sequence encoding the peptides of RIP3 that are frequently detected by mass spectrometry. RIP3 fragment DNA and the full-length TRADD DNA sequences were cloned into overexpression vectors that harbor 6× His tag. Overexpression plasmids were transfected into 293T. After 48 h, the cells were lysed with 8 mol/L urea and the target proteins were purified using Ni-NTA column and eluted with 250 mmol/L imidazole. The purified proteins were separated on SDS-PAGE for purity verification (Fig. S1I). The absolute amount of purified proteins was determined by AAA-MS (Table S3).
Ten 15-cm dish cells were collected for each time point experiment, and three biological replicates were carried out. After 3× Flag-TNF or TNF treatment, cells were immediately washed twice with ice-cold PBS and harvested by scraping. The harvested cells were washed with PBS and lysed for 30 min on ice in HBS lysis buffer (12.5 mmol/L HEPES, 150 mmol/L NaCl, 1% Nonidet P-40, pH 7.5) with the protease inhibitor cocktail. Cell lysates were then spun down at 20,000 ×g for 30 min. The soluble fraction was collected, and immunoprecipitated overnight with anti-Flag M2 antibody-conjugated agarose at 4°C. Resins containing protein complexes were washed three times with HBS lysis buffer. Proteins were subsequently eluted twice with 0.2 mg/mL of 3× flag peptides in HBS lysis buffer for 30 min each time, and elution was pooled for a final volume of 300 μL. Proteins in the elution were precipitated with 20% RIP1-dependent linear and nonlinear recruitments RESEARCH ARTICLE trichloroacetic acid (TCA) and the pellet was washed two times with 1-mL cold acetone and dried in speedvac.
TCA-precipitated proteins were re-suspended in 50 μL 1% SDC (sodium deoxycholate) in 10 mmol/L TCEP (Tris (2-carboxyethyl) phosphine hydrochloride), 40 mmol/L CAA (chloroacetamide), 100 mmol/L Tris-HCl pH 8.5. The appropriate amount of purified heavy amino acid labeled TRADD was spiked into the TNFR1 and RIP1 IP samples, and purified heavy amino acid labeled RIP3 was spiked into RIP3 IP samples. IP samples were subsequently shaken at 37°C for 30 min. Protein concentration was measured with 660 nm Protein Assay Reagent (Piece, Thermo). 1% SDC was then diluted to 0.5% SDC with water, and trypsin (Sigma) was added at the ratio of 1:50 (trypsin: protein). The digestion was performed at 37°C overnight. 1% TFA was added to stop the reactions followed by centrifugation at 20,000 ×g for 10 min. The supernatants were subsequently transferred to C18 STAGEtips for desalting. Tryptic peptides were eluted with 70% ACN/1% FA and dried in speedvac.

Absolute quantification of RIP1 protein amount in individual cells
To obtain the absolute amount of RIP1 in individual L929 cell, the RIP1 peptide (HQAIFDNTTSLTDEHLNPIR) containing "heavy" 13C6, 15N2 L-Lysine and 13C6, 15N4 L-Arginine was purchased from the Go Top Peptide Biotech Company (Hangzhou, China). The number of L929 cells was counted and collected followed by threetime PBS wash. Three biological experiments were performed. Cells were lysed with 1% SDC in 10 mmol/L TCEP (Tris(2-carboxyethyl) phosphine hydrochloride), 40 mmol/L CAA (chloroacetamide), 100 mmol/L Tris-HCl pH 8.5, followed by 37°C for 30 min. The lysate was cleared by centrifugation. Protein concentration was assayed with 660 nm protein Assay Reagent (Piece, Thermo). 100 μg proteins were subjected to trypsin digestion. The digestion was performed at 37°C overnight. Heavy amino acid labeled peptide was added into the reactions. 1% TFA was added to stop the reactions followed by centrifugation at 20,000 ×g for 10 min. The supernatants were subsequently transferred to C18 STAGEtips for desalting. Tryptic peptides were eluted with 70% ACN/1% FA and dried in SpeedyVac. The L929 cell lysate peptides were dissolved in 0.1% FA and about 5 μg peptides (we assume digestion efficiency is 100%) were measured on mass spectrometry in SWATH-MS mode. The LC gradient time is 180 min, and SWATH-MS setting is 100 VW (See below for detail).
For DDA, 250-ms MS1 scan was performed in the range of 350-1,250 m/z, and up to 20 most intense precursors with charge state Group-DIA analysis SWATH wiff files were converted to profile mzXML files with MSconvert (Proteowizard, Version 3.0.4472). The mzXML files were split into 1 MS1 mzML and 32 or 100 MS2 mzML files using in-house scripts. These mzML files were input into Group-DIA analysis. 18 runs in TNFR1 dataset, 30 runs in RIP1 dataset and 17 runs in RIP3 dataset were separately analyzed by Group-DIA. Group-DIA composed of four modules, "alignment", "analysis", "identification" and "validation". For internal spectral library building, only "alignment" and "analysis" modules were performed. MS1 profiles in multiple runs were first aligned, and XICs of precursors and product ions were concatenated for multiple runs. Similarity comparison was performed between the concatenated XICs of precursors and product ions. The precursors and product ions pairs were extracted and stored in mgf and mzML formats.

Spectral library building
To enable in-depth exploration of SWATH data, we combined internal spectral library and DDA spectral library. Group-DIA-generated mgf files were converted to mzXML files using TPP msconvert (Trans Proteomics Pipeline, Version 4.8). DDA wiff files were converted to centroided mzXML files with qtofpeakpicker (Trans Proteomics Pipeline, Version 4.8). These mzXML files were searched with Comet (Version 2017.01) and X!tandem (Version 2013.06.15.1,

RESEARCH ARTICLE
Xiang Li et al. native and k-score) against the full non-redundant, canonical mouse genome as annotated by UniprotKB/Swiss-Prot (downloaded in September, 2014) appendant with common contaminants and reversed sequence decoys (33,864 sequences includes decoys). The search parameters were set as followed, parent monoisotopic tolerance 50 ppm, modification 57.021464@C, potential modification 15.994915@M and maximum missed cleavage sites 2. The pep.xml files were validated with PeptideProphet and combined with iProphet. The Mayu (version 1.07) was utilized for FDR estimation. The peptide ions filtered at 1% protein FDR were imported into spectraST for spectral library building. The retention time of peptides in sptxt file was replaced with iRT time using spectrast2spectrast_irt. py script, where ciRT peptides were used for retention time normalization. Subsequently, the sptxt files were made consensus nonabundant sptxt files with spectraST. Sptxt files were converted to Peakview software compatible txt files using spectrast2tsv.py script.

Targeted analysis of SWATH-MS data using Peakview software
The spectral libraries were converted to Peakview-compatible file. All b-ions from heavy amino acid labeled peptides were filtered out from the library. The parameters of Peakveiw software (Version 2.2) were set as followed. The peptides from the bait proteins and iRT peptides were used for retention time alignment between runs. "Number of transitions per peptide" was "6", "False Discovery Rate" was "1%", "XIC Extraction window" was "10 min", and "XIC width (Da)" was "0.05". cultures were examined, and similar results were obtained in at least three independent experiments.