Dynamic NF-κB and E2F interactions control the priority and timing of inflammatory signalling and cell proliferation
Figures
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig1-v1.tif/full/617,/0/default.jpg)
NF-κB dynamics following TNFα treatment in HeLa and SK-N-AS cells: Mapping the NF-κB response over the cell cycle in synchronized HeLa cells.
(A,B,C and D) The dynamics of RelA-dsRedxp following 10 ng/ml TNFα treatment in transiently transfected SK-N-AS (A), or following 30 pg/ml TNFα treatment in SK-N-AS, and 10 ng/ml TNFα treatment in HeLa cells (C), and at 30 pg/ml for HeLa (D) cells (n=30 cells analysed per condition). (E) The localization of endogenous RelA in different cell cycle phases, observed by immunocytochemistry at 2 hr (G1/S transition), 4 hr (mid S-phase), post-release from double thymidine block and with 15 min TNFα treatment. (F and G) The dynamics of RelA-dsRedxp in transiently transfected HeLa cells synchronized by a double thymidine block, following 10 ng/ml TNFα treatment at G1/S (F), or passing through S-phase (G) (n=20 cells analysed per condition). (H) Western blot of Ser536phopho-RelA (p-RelA), IκBα, and cyclophilin-A (cyclo-A) levels in synchronized HeLa cells harvested at 1 hr time intervals over the G1/S transition following 15 min treatment with TNFα. Also shown are asynchronous, non-stimulated (ASY NST) and asynchronous, stimulated (ASY ST) controls, harvested at t=0.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig2-v1.tif/full/617,/0/default.jpg)
Mapping the NF-κB response over the cell cycle through virtual synchronization.
(A) Selected images from time-lapse imaging of RelA-dsRedxp transiently expressing Hela cells treated with 10 ng/ml TNFα. (B) Virtual synchronization of HeLa cells treated with 10 ng/ml TNFα. Cells were imaged through two successive divisions (M) allowing correlation of cell cycle timing of TNFα treatment (parameter 1) to RelA dynamics (parameters 2, 3 and 4) and cell cycle duration (parameters 1 plus 5). (C) Representative cells of RelA-dsRedxp dynamics following TNFα treatment in asynchronous cells, then virtually synchronized into G1 (n=115), G1/S (n=32), S (n=52) and G2 (n=38) phases.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig2-figsupp1-v1.tif/full/617,/0/default.jpg)
Analysis of cell cycle duration and G1/S timing in HeLa and SK-N-AS cells.
(A) Time series for FUCCI expression in single representative HeLa and SK-N-AS cells. White arrows mark cells before and after the fluorescence levels were detectable. (B) Analysis of cells in (A), showing the G1/S crossing point in fluorescence levels from reporters of SCF (SKP2) (Orange) and APC (Green) E3 ubiquitin ligase activity. (C) Analysis of the G1/S crossing point and cell cycle duration in populations of HeLa and SK-N-AS cells transfected with FUCCI vectors (n ≥ 30 cells for all conditions).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig2-figsupp2-v1.tif/full/617,/0/default.jpg)
Statistical analysis of NF-κB translocation in HeLa cells at inferred cell cycle stages following 10 ng/ml TNFα stimulation.
(A) Analysis of dynamics of initial RelA-dsRedxp translocation with respect to cell cycle phase, using virtual synchronization in HeLa cells. Data were analyzed using nonparametric Anova analysis with Dunn correction for multiple comparisons. Red lines indicate mean normalised amplitude of NF-κB nuclear translocation for different cell cycle phases, and the population average (dotted line). (B) Analysis of nuclear RelA occupancy assessed in non-synchronized cells expressing RelA-dsRedxp following treatment with 10 ng/ml TNFα. Statistical analysis showed significant difference between cell cycle phases with respect to distribution of amplitude of the response (Anova analysis with Dunn correction for multiple comparisons).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig2-figsupp3-v1.tif/full/617,/0/default.jpg)
Statistical analysis of NF-κB translocation in SK-N-AS cells at inferred cell cycle stages following 30 pg/ml TNFα stimulation.
(A) Correlation of estimated cell cycle timing with RelA-dsRedxp N:C peak amplitude following 30 pg/ml TNFα treatment (n=138). (B) Analysis of dynamics of initial RelA-dsRedxp translocation with respect to cell cycle phase. Statistical analysis showed a difference between G1 and S with respect to distribution of amplitude of the response (Anova analysis with Dunn correction for multiple comparisons).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig3-v1.tif/full/617,/0/default.jpg)
Cell cycle length and variability is modified by TNFα addition at G1/S.
Analysis of the timing and variability of mitosis (parameter 1 plus 5 from Figure 2B) following 10 ng/ml TNFα treatment of asynchronous untransfected HeLa cells, compared to subsets of those cells stimulated at late G1- or S-phase. Mean durations were analysed using nonparametric Anova analysis with Dunn correction for multiple comparisons. Variability in the data was analysed using Levene’s test for equality of variance.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig4-v1.tif/full/617,/0/default.jpg)
Physical and functional interaction between NF-κB and E2F-1 systems.
(A) NF-κB-dependent transcription was assessed by luciferase reporter assay (NF-luc), in SK-N-AS cells (n=3, +/- s.d) expressing EGFP-E2F-1, RelA-dsRedxp or both. (B) IκBα and IκBε mRNA levels in SK-N-AS cells (n=3, +/- s.d) following transient expression of EGFP-E2F-1, RelA-DsRedxp or both. (C) E2F-1-dependent transcription as assessed by luciferase reporter assay (CyclinE-luc), in SK-N-AS cells (n=3, +/- s.d) expressing EGFP-E2F-1, RelA-dsRedxp or both. (D) E2F-1 mRNA levels in SK-N-AS cells (n=3, +/- s.d) transiently transfected with RelA-dsRedxp. (E) Representative SK-N-AS cells transiently expressing EGFP-E2F-1 (green), RelA-dsRedxp (red), both fluorescent fusion proteins at different levels, or EGFP-E2F-1, RelA-dsRedxp and IκBα-AmCyan (blue).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig4-figsupp1-v1.tif/full/617,/0/default.jpg)
E2F-1 modulates NF-κB dynamics in the absence of stimulus in SK-N-AS cells.
(A) Time-lapse confocal microscopy of representative SK-N-AS cells transiently transfected with RelA-dsRedxp and EGFP-E2F-1. (B) Trajectories of three representative cells expressing different levels of EGFP-E2F-1. (C) Correlation between RelA-dsRedxp T½ nuclear occupancy (NO) time and EGFP-E2F-1 T½ nuclear degradation time, based on data in (A). (D) Recapitulation of the observed dynamics with an in silico model for physical interaction between RelA (NFkB) and E2F-1 (E2F) (E) Correlation between NF-κB nuclear occupancy time and nuclear E2F-1 degradation time, based on data in (D) (n= 30 cells).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig4-figsupp2-v1.tif/full/617,/0/default.jpg)
E2F-1 modulates NF-κB dynamics in the absence of stimulus in HeLa cells.
(A) Representative HeLa cells transiently transfected with combinations of RelA and E2F-1 fluorescent fusion proteins. (B) Time-lapse confocal microscopy of representative HeLa cells transiently transfected with RelA-dsRedxp and EGFP-E2F-1. (C) Trajectories of three representative cells expressing different levels of EGFP-E2F-1. (D) Correlation between RelA-dsRedxp T½ nuclear occupancy (NO) time and EGFP-E2F-1 T½ nuclear degradation time, based on data in (C) (n=20).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig5-v1.tif/full/617,/0/default.jpg)
Interaction of E2F-1 with RelA.
(A) Representative cell demonstrating co-localisation of E2F1-EGFP and RelA-dsRedxp upon transient transfection. (B) Co-Immunoprecipitation of E2F-1 with RelA pull down in HeLa cells synchronized in late G1 (HeLa cells used for this experiment due to their greater ease of synchronization). (C) FCCS assay between transiently transfected EGFP-E2F-1 and RelA-dsRedxp (red line) or empty-dsRedxp (blue line) fluorescent fusion proteins in single live SK-N-AS cells (+/- s.e.m based on 10 measurements from 10+ cells per condition). (D) Qualitative FRET assay between transiently transfected ECFP-E2F-1 and RelA-EYFP fluorescent fusion proteins in live SK-N-AS cells. First negative control between IkB-ECFP and EYFP-E2F1, and second negative control between free ECFP and EYFP fluorophores expressed in an SK-N-AS cell (shown are average ECFP and EYFP signals (+/- s.e.m based on 20 cells per condition normalised to pre-bleach intensity. p.b. indicates the time point at which photo-bleaching occurred).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig6-v1.tif/full/617,/0/default.jpg)
Mathematical modelling predicts an additional key component for NF-κB - cell cycle interactions: E2F-4 identified as a putative candidate.
(A) Model simulations of RelA-dsRedxp dynamics when co-expressed with EGFP-E2F-1 in cells treated with TNFα. (B) Dynamics analysed in representative SK-N-AS cells treated with 10 ng/ml TNFα expressing RelA-dsRedxp and EGFP-E2F-1 (C) Model simulation of experimental conditions in B, incorporating interactions between NF-κB complexes and a putative E2F-1-induced target protein, subsequently proposed as E2F-4. (D) Analysis of average timing to second peak of NF-κB translocation following TNFα treatment in SK-N-AS cells expressing RelA-dsRedxp alone or with EGFP-E2F-1 (n=20 cells per condition, error bars show s.d.) (E) Assessment of the extent of RelA Ser536 phosphorylation (p-RelA), E2F-4 and IκBα stability by western blot compared to cyclophilin A (cyclo A) amounts in SK-N-AS cells either untreated or treated with 10 ng/ml TNFα and expressing combinations of either untagged or fluorescent RelA-dsRedxp and EGFP-E2F-1. (F) Western blot of E2F-1 and E2F-4 in synchronized HeLa cells, where t=0 is late G1-phase.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig7-v1.tif/full/617,/0/default.jpg)
E2F-4 directly interacts with NF-κB and perturbs RelA dynamics in response to TNFα stimulation.
(A) Single cell trajectories from groups of HeLa cells expressing RelA-dsRedxp and different levels of EGFP-E2F-4 showing the dynamics of RelA-dsRedxp after 10 ng/ml TNFα treatment (n=60 cells). (B) HeLa cells synchronized in S-phase, co-immunoprecipitated with anti-RelA antibody and probed for E2F-4. Also shown are IgG negative controls and whole cell lysate unsynchronized positive control (ctrl). (C) Representative SK-N-AS cells transiently transfected with RelA-dsRedxp and EGFP-E2F-4. (D) FRET assay in live SK-N-AS cells expressing ECFP-E2F-4 and RelA-EYFP fluorescent fusion proteins (shown are average ECFP and EYFP signals (+/- s.e.m) based on 20 cells per condition normalised to pre-bleach intensity. p.b. indicates the point of photo-bleaching). (E) FCCS assay in cells transiently expressing EGFP-E2F-4 and RelA-dsRedxp (red line) or dsRedxp (blue line) fluorescent proteins in single live SK-N-AS cells (+/- s.e.m based on 10 measurements in each of 10+ cells per condition).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig7-figsupp1-v1.tif/full/617,/0/default.jpg)
Analysis of RelA-dsRedxp dynamics in HeLa and SK-N-AS cells co-expressing EGFP-E2F-4 following TNFα stimulation.
The effects of different EGFP-E2F-4 expression levels on the amplitude and timing of the first peak of RelA translocation in HeLa and SK-N-AS cells treated with 10 ng/ml and 30 pg/ml TNFα, respectively. These data indicate how ectopically expressed EGFP-E2F-4 can inhibit the translocation of RelA-dsRedxp in response to TNFα.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig8-v1.tif/full/617,/0/default.jpg)
Effect of cell cycle timing on RelA-dsRedXP translocation in dual BAC HeLa cells (C1-1 line) that co-express E2F-1-Venus fusion protein.
(A) Selected images from time-lapse experiment of dual BAC transfected HeLa stable clone 1-1 showing translocation of RelA-dsRedXP and E2F-1-Venus expression at different cell cycle phases. Cells were treated with 10 ng/ml TNFα. (B) Analysis of the dynamics of initial RelA-dsRedxp translocation in cells ordered at specific cell cycle times with respect to the peak of E2F-1 expression (n = 128). Data were analysed using nonparametric Anova analysis with Dunn correction for multiple comparisons. Red lines indicate mean normalised amplitude of NF-κB nuclear translocation for different cell cycle phases, and the population average (dotted red line). Analysis of nuclear RelA occupancy was assessed in virtually synchronised C 1-1 cells, based on time from cell division and relative to peak E2F-1-Venus expression level. RelA-dsRedxp localization was visualized to allow quantification of translocation, following treatment with 10 ng/ml TNFα. The dotted black line shows the spline fitted level of E2F1 at different times and cell cycle stages (see also Figure 8—figure supplement 1 below). Statistical analysis showed a difference between G1 vs S, and G2 vs S with respect to distribution of amplitude of the RelA translocation response. (C) RelA-dsRedxp dynamics following 10 ng/ml TNFα treatment in asynchronous cells (left panel) and cells virtually synchronised into G1, G1/S, S and G2 phases. The data for each cell was normalised to the amplitude (N:C ratio) at t = 0 min.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig8-figsupp1-v1.tif/full/617,/0/default.jpg)
Virtually synchronized HeLa C 1-1 cells.
Normalised E2F-1-Venus expression at the time of TNFα stimulation of C1-1 cells (data also shown in Figure 8B). E2F-1-Venus expression was normalised to its peak expression. The time axis represents the time of TNFα stimulation relative to the peak of E2F-1-Venus for each cell, where time 0 is the peak of E2F-1-Venus expression. Positive times indicate stimulation after the peak of E2F-1-Venus expression and negative values indicate stimulation events before the peak of E2F-1-Venus expression. The black line shows a spline interpolation of the level of E2F-1-Venus expression. Cell cycle phases were estimated based on measured the E2F1 profile and average cell cycle timing.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig8-figsupp2-v1.tif/full/617,/0/default.jpg)
Physiological and functional expression of E2F-1-Venus in stable BAC-transduced HeLa cells.
(A) HeLa cells stably expressing an E2F-1-Venus fluorescent fusion protein from a 5KB endogenous E2F-1 promoter (Green), transiently transfected with a FUCCI reporter for SCF (SKP-2) activity (Orange). Showing the profile of E2F-1 over two consecutive cell cycles (one parent and two daughter cells), with a peak in late G1. E2F-1 levels dropped during S-phase consistent with rapid rise in SCF (SKP-2) activity and a loss of FUCCI fluorescence. (B) Representative cell from the E2F-1-Venus and RelA-DsRedxp stably transfected population of Hela cells through one full cell cycle.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig8-figsupp3-v1.tif/full/617,/0/default.jpg)
Analysis of the expression of E2F-1-Venus and RelA-DsRedxp translocation in single C1-1 HeLa cells stimulated with 10ng/ml TNFα at different cell cycle phases.
Grey line shows the E2F-1-Venus expression level plotted agains the right y-axis. The red, green blue and orange lines show the timecourse of RelA-dsRedxp localization in exemplar cells in the G1, G1/S, S and G2 phases respectively plotted against the left y-axis. The black vertical line represents the point at which cells were treated with TNFα.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig8-figsupp4-v1.tif/full/617,/0/default.jpg)
Expression and interaction of RelA-dsRedxp and E2F-1-Venus.
(A) Western blot of RelA and α-tubulin levels in dual BAC stable C1-1 and WT HeLa showing exogenous expression of RelA-dsRedxp. (B) Western blot of E2F-1 and cyclophilin A levels in dual BAC stable C1-1 and WT HeLa cells expressing the E2F-1-Venus BAC (C) Total fluorescent molecules per cell for E2F1-Venus at peak expression and RelA-dsRedxp in unstimulated cells (data obtained from Fluorescent Correlation Spectroscopy measurements, and calculated using volume estimates from z-stacked WT HeLa in suspension). (D) FCCS mean correlation curves (+/- s.e.m) between E2F-1-Venus and RelA-dsRedxp (red line, n=46) for TNFα treated BAC stable cells. A comparison to transient empty-dsRedxp is shown (blue line, n=15) (E) Kd determination results using a scatter plot and linear regression (Theil-Sen estimator). The slope of the regression gives the Kd value.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-fig9-v1.tif/full/617,/0/default.jpg)
Schematic representation of NF-κB – E2F interactions.
(A) Predicted mechanisms for NF-κB interaction with E2F proteins over the G1/S transition (B) Model simulations of single cell behaviour.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig1-v1.tif/full/617,/0/default.jpg)
Oscillations in the NF-κB system.
(A) Dynamics of RelA-dsRedxp in transiently transfected SK-N-AS cells following 10 ng/μl TNFα stimulation, plotted over 450 and 150 min respectively. (B) Dynamics of RelA-dsRedxp in transiently transfected HeLa cells following 10 ng/μl TNFα stimulation, plotted over 450 and 150 min respectively.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig2-v1.tif/full/617,/0/default.jpg)
Use of double-Thymidine block to synchronize HeLa cells at G1/S.
Flow cytometric analysis of the distribution of DNA content of non-synchronized HeLa cells and cells harvested at relevant times post-release from Thymidine block.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig3-v1.tif/full/617,/0/default.jpg)
Effect of Double-Thymidine block on tagged and endogenous RelA levels HeLa cells.
Endogenous RelA and tagged RelA-DsRedxp expression levels in unsynchronised and synchronised WT-HeLa and double BAC stable cells. Synchronized fractions at 0, 2 and 4 hr post release of thymidine block. α-Tubulin used as loading control.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig4-v1.tif/full/617,/0/default.jpg)
Negative Co-IP.
(A) Co-Immunoprecipitation of E2F-1 with RelA (pulled down with a RelA antibody) in asynchronous HeLa cells showing no detectable band of E2F-1. (B) Co-Immunoprecipitation of E2F-4 with RelA (pulled down with a RelA antibody) in HeLa cells synchronized in late G1-phase cells showing no detectable band of E2F-4.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig5-v1.tif/full/617,/0/default.jpg)
A typical simulation protocol of the NF-κB:E2F-1 mathematical model.
Simulation protocol of a live cell imaging experiment involving transfection and TNFα stimulation.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig6-v1.tif/full/617,/0/default.jpg)
Simulations from the NF-κB:E2F model.
Blue lines represent E2F1= 0, TR= 1. Black lines represent E2F1= NFkB= 0.1 and TR=0. Red lines represent E2F1= NFkB= 0.1 and TR=1.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig7-v1.tif/full/617,/0/default.jpg)
Cell Cycle length of Clonal HeLa BAC population.
(A) Analysis of cell cycle duration in populations of dual BAC stable cell line (C1-1) with wild type HeLa cells. (B) Analysis of the effects of TNFα treatment in C1-1 and WT HeLa cells.
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig8-v1.tif/full/617,/0/default.jpg)
FCCS control.
FCCS assays between transiently transfected RelA-dsRedxp and IκBα-EGFP (red line) and Empty-DsRedxp and Empty-EGFP (blue line) in single live SK-N-AS cells (+/- s.e.m based on 10 measurements in each of 10+ cells per condition).
![](https://iiif.elifesciences.org/lax/10473%2Felife-10473-app1-fig9-v1.tif/full/617,/0/default.jpg)
FCCS autocorrelation analysis.
Autocorrelation lines for RelA/IkBα, RelA/E2F-1 and RelA/E2F-4 FCCS studies in single live SK-N-AS cells (+/- s.e.m based on 10 measurements in each of 10+ cells per condition).
Tables
A Initial conditions for the NF-κB:E2F mathematical model. Concentrations of NF-κB and E2F-1 added to mimic cell transfection are also included.
Species | Biological name | Initial Conditions Equilibration stage(μM) | Initial Conditions TNFα stimulation (μM) |
---|---|---|---|
NFkB | Cytoplasmic RelA | 0 | 0.004 (+0.1) |
nNFkB | Nuclear RelA | 0 | 0.015 |
E2F1 | Cytoplasmic E2F1 | 0 | 0 (+0.1) |
nE2F1 | Nuclear E2F1 | 0 | 0 |
tIkBa | IκBα mRNA | 0 | 1e-005 |
IkBa | Cytoplasmic IκBα | 0 | 0.017 |
nIkBa | Nuclear IκBα | 0 | 0.004 |
IKKn | Neutral IKK | 0.1 | 0.1 |
IKK | Active IKK | 0 | 0 |
IKKi | Inactive IKKi | 0 | 0 |
tA20 | A20 mRNA | 0 | 1e-005 |
A20 | A20 | 0 | 0.001 |
pIkBa | phospho-IκBα | 0 | 0 |
pIkBaNFkB | phospho-IκBα RelA complex | 0 | 0 |
NFkBE2F1 | cyto. RelA E2F-1 complex | 0 | 0 |
nNFkBE2F1 | nuclear RelA E2F-1 complex | 0 | 0 |
IkBaNFkB | cyto IκBα RelA complex | 0.1 | 0.091 |
nIkBaNFkB | nuclear IκBα RelA complex | 0 | 0.001 |
tE2F4 | E2F-1 target E2F-4 mRNA | 0 | 0 |
E2F4 | E2F-1 target E2F-4 | 0 | 0 |
E2F4NFkB | E2F-4 RelA complex | 0 | 0 |
E2F4IkBaNFkB | E2F-4 IκBα RelA complex | 0 | 0 |
NF-κB:E2F-1 model equations Symbol ‘n’ denotes nuclear variables,‘t’ denotes mRNA transcripts, ‘p’ denotes phosphorylated form of IκBα. Symbols denoting cytoplasmic localisation were omitted.
Model reactions and associated parameters.
Reaction | Symbol | Value | References |
---|---|---|---|
Spatial parameters | |||
Total cell volume | tv | 2700 µm3 | Measured |
C:N ratio | kv | 3.3 | Measured |
Conversion to nuclear volume | nv | ×(kv+1) | - |
Conversion to cytoplasmic volume | cv | ×(1/kv+1) | - |
Initial concentration | |||
Total NF-κB | NF | 0.08 µM | Initialized as cytoplasmic IκBα·NF-κB |
Total IKK | - | 0.08 µM | Initialized as IKKn |
Complex formation & dissociation | |||
IκBα + NF-κB → IκBα·NF-κB nIκBα + nNF-κB → nIκBα·NF-κB | ka1a | 0.5 µM-1s-1 | (Hoffmann et al., 2002) |
IκBα·NF-κB → IκBα + NF-κB nIκBα·nNF-κB → nIκBα + nNF-κB | kd1a | 0.0005 s-1 | (Hoffmann et al., 2002) |
NF-κB + E2F (1 or 4) → NF-κB·E2F nNF-κB + nE2F → nNF-κB·nE2F | ka2e | 0.5 µM-1s-1 | fitted, same as IκBα + NF-κB |
NF-κB·E2F → NF-κB + E2F nNF-κB·nE2F → nNF-κB + nE2F | kd2e | 0.0005 s-1 | fitted, same as IκBα + NF-κB |
NF-κB·E2F1 + IκBα→ IκBα·NF-κB + E2F1 nNF-κB·nE2F1 + nIκBα → nIκBα·NF-κB + nE2F1 | kdis | 0.001 s-1 | fitted |
Transport | |||
NF-κB → nNF-κB | ki1 | 0.0026 s-1 | Measured fitting range: Average 0.0026 ± 0.0018s-1 |
nNF-κB → NF-κB | ke1 | 0.000052 s-1 | ki1/50 (Carlotti et al., 2000) |
E2F1→ nE2F1 | kie | 0.0026 s-1 | fitted, same as NF-κB |
nE2F1 → E2F1 | kee | 0.000052 s-1 | fitted, same as NF-κB |
IκBα → nIκBα | ki3a | 0.00067 s-1 | Measured fitting range: Average 0.00043 ± 0.00024 s-1 |
nIκBα → IκBα | ke3a | 0.000335 s-1 | ki3a/2 (Carlotti et al., 2000) |
nIκBα·nNF-κB → IκBα·NF-κB | ke2a | 0.01 s-1 | Fitted |
NF-κB·E2F1 → nNF-κB·nE2F1 | kine | 0.0026 s-1 | fitted, same as NF-κB |
nNF-κB·nE2F1 → NF-κB·E2F1 | kene | 0.000052 s-1 | fitted, same as NF-κB |
Protein synthesis & degradation | |||
nNF-κB → nNF-κB + tIκBα Order of hill function, h=2 Half-max constant, k=0.065h(fitted) | c1a | 1.4×10-7 µM-1s-1 | Fitted (constrained): 1.07×10-7 – 8.2×10-7µM-1s-1 (Femino et al., 1998); (Cheong et al., 2006) |
tIκBα→ tIκBα + IκBα | c2a | 0.5 s-1 | (Lipniacki et al., 2004) |
NF-κB·IκBα→ NF-κB | c5a | 0.000022 s-1 | (Pando and Verma, 2000; Mathes et al., 2008) |
nNF-κB·nIκBα→ nNF-κB | - | 0 s-1 | Assumed (O'Dea et al., 2007; Mathes et al., 2008) |
nNF-κB → nNF-κB + tA20 Order of hill function, h=2 Half-max constant, k=0.065h | c1 | 1.4×10-7 µM-1s-1 | Assumed to be the same as IκBα |
nE2F1 → nE2F1 + tE2F4 Order of hill function, h=2 Half-max constant, k=0.065h | c1x | 9.8×10-7 µM-1s-1 | Fitted |
tA20→ tA20 + A20 | c2 | 0.5 s-1 | - |
tE2F-4→ tE2F-4 + E2F4 | c2x | 0.5 s-1 | - |
tIκBα→ Sink | c3a | 0.0003 s-1 | Fitted (constrained): 0.00077-0.00029 s-1 (Blattner et al., 2000) |
tA20→ Sink | c3 | 0.00048 s-1 | Fitted, constrained >tIκBα turnover (Ashall et al., 2009) |
tE2F4→ Sink | c3x | 0.00048 s-1 | Fitted |
IκBα→ Sink | c4a | 0.0005 s-1 | Fitted (constrained): 0.000105 – 0.002 s-1 (Pando and Verma, 2000; O'Dea et al., 2007; Mathes et al., 2008) |
A20 → Sink | c4 | 0.0045 s-1 | Fitted |
E2F4 → Sink | c4x | 0.00016 s-1 | Fitted |
E2F1 → Sink | c6e | 0.00016 s-1 | Fitted |
nE2F1 → Sink | c7e | 0.00016 s-1 | Fitted |
NF-κB·E2F1 → Sink | c8ne | 0.00016 s-1 | Fitted |
nNF-κB·nE2F1 → Sink | c9ne | 0.00016 s-1 | Fitted |
TNFα stimulation | |||
TNFα | TR | 1/0 | on/off (Lipniacki et al., 2004) |
IKK parameters | |||
IKKn → IKKa | ka | 0.004 s-1 | Fitted, as above |
IKKa → IKKi | ki | 0.003 s-1 | Fitted, as above |
IKKi → IKKn | kp | 0.0006 s-1 | Fitted |
A20 inhibition rate constant | kbA20 | 0.0018 | Fitted, scales kp dependent on receptor state kbA20×TR |
IKKa + IκBα → pIκBα | kc1a | 0.074 s-1 | Assumed (0.037×2) (Heilker et al., 1999) |
IKKa + IκBα·NF-κB → pIκBα·NF-κB | kc2a | 0.37 s-1 | Assumed (0.037×5×2) (Heilker et al., 1999; Zandi and Karin, 1999) |
pIκBα → Sink | kt1a | 0.1 s-1 | Fitted |
pIκBα·NF-κB → NF-κB | kt2a | 0.1 s-1 | Fitted |
Simulation protocols used throughout the manuscript. TNFα stimulation is invoked via TR=0/1. E2F refers to levels of 'transfection' (in μM). E2F4 off/on refers to whether its transcription is switched on or off.
Figure | Model conditions |
---|---|
3E | TR=0, E2F1= (0.05, 0.1, 0.15), E2F4 off |
4A | TR=1, E2F1 = 0.1, E2F4 off |
4C | TR=1, E2F1 = 0.1, E2F4 on |
4I (G1, G2) | TR=1, E2F1= 0, E2F4 on (but unaffected) |
4I (G1/S) | TR=1, E2F1 = 0.2, E2F4 on |
4I (S) | TR=1, NFkBIkBa = 0.1, E2F1 = 0, E2F4= 0.1 |