A multiscale model of the CD8 T cell immune response structured by intracellular content

During the primary CD8 T cell immune response, CD8 T cells undergo proliferation and continuous differentiation, acquiring cytotoxic abilities to address the infection and generate an immune memory. At the end of the response, the remaining CD8 T cells are antigen-specific memory cells that will respond stronger and faster in case they are presented this very same antigen again. We propose a nonlinear multiscale mathematical model of the CD8 T cell immune response describing dynamics of two inter-connected physical scales. At the intracellular scale, the level of expression of key proteins involved in proliferation, death, and differentiation of CD8 T cells is modeled by a delay differential system whose dynamics define maturation velocities of CD8 T cells. At the population scale, the amount of CD8 T cells is represented by a discrete density and cell fate depends on their intracellular content. We introduce the model, then show essential mathematical properties (existence, uniqueness, positivity) of solutions and analyse their asymptotic behavior based on the behavior of the intracellular regulatory network. We numerically illustrate the model's ability to qualitatively reproduce both primary and secondary responses, providing a preliminary tool for investigating the generation of long-lived CD8 memory T cells and vaccine design.