Fig. 1

Reproduction of intracellular and cell dynamics data characterizing an early CD8+ T cell immune response. a Kinetics of IL-2, IL-2R, IL-2•IL-2R complex, T-bet, Fas* and cleaved Caspase from 72 h pi to 120 h pi. Molecular concentration is represented by squares (in vivo mRNA data, corresponding to the right-hand y-axis; n = 3; Caspase data are Caspase-3 data from ImmGen) or solid lines (simulation results, corresponding to the left-hand y-axis; mean value is indicated by a solid line, max and min values are upper and lower boundaries of the shaded areas respectively; n = 10 simulations). b CD8⁺ T cell dynamics, from 72 h pi to 120 h pi. Fold change of CD8⁺ T cell number is represented by cross marks for in vivo data (mean ± SD; cells were collected from 4 mice at each time point) and by a solid line for simulation results (mean value is indicated by a solid line, max and min values are upper and lower boundaries of the shaded area respectively; n = 10 simulations). c Proportion of the phenotypes within a simulated CD8⁺ T cell population (mean; n = 10 simulations). d Representative in silico simulation snapshots of 72 h, 96 h and 120 h pi. Color-coded are APCs (green), naïve cells (light blue), pre-activated cells (dark blue), activated cells (orange), effector cells (red) and apoptotic cells (black)