Abstract:Mathematical modeling provides a rigorous way to quantify immunological processes and discriminate between alternative mechanisms driving specific biological phenomena. It is typical that mathematical models of immunological phenomena are developed by modelers to explain specific sets of experimental data after the data have been collected by experimental collaborators. Whether the available data are sufficient to accurately estimate model parameters or to discriminate between alternative models is not typically investigated. While previously collected data may be sufficient to guide development of alternative models and help estimating model parameters, such data often do not allow to discriminate between alternative models. As a case study, we develop a series of power analyses to determine optimal sample sizes that allow for accurate estimation of model parameters and for discrimination between alternative models describing clustering of CD8 T cells around Plasmodium liver stages. In our typical experiments, mice are infected intravenously with Plasmodium sporozoites that invade hepatocytes (liver cells), and then activated CD8 T cells are transferred into the infected mice. The number of T cells found in the vicinity of individual infected hepatocytes at different times after T cell transfer is counted using intravital microscopy. We previously developed a series of mathematical models aimed to explain highly variable number of T cells per parasite; one of such models, the density-dependent recruitment (DDR) model, fitted the data from preliminary experiments better than the alternative models, such as the density-independent exit (DIE) model. Here, we show that the ability to discriminate between these alternative models depends on the number of parasites imaged in the analysis; analysis of about [Formula: see text] parasites at 2, 4, and 8 h after T cell transfer will allow for over 95% probability to select the correct model. The type of data collected also has an impact; following T cell clustering around individual parasites over time (called as longitudinal (LT) data) allows for a more precise and less biased estimates of the parameters of the DDR model than that generated from a more traditional way of imaging individual parasites in different liver areas/mice (cross-sectional (CS) data). However, LT imaging comes at a cost of a need to keep the mice alive under the microscope for hours which may be ethically unacceptable. We finally show that the number of time points at which the measurements are taken also impacts the precision of estimation of DDR model parameters; in particular, measuring T cell clustering at one time point does not allow accurately estimating all parameters of the DDR model. Using our case study, we propose a general framework on how mathematical modeling can be used to guide experimental designs and power analyses of complex biological processes.

A Nonconventional Analysis of CD$4^{+}$ and CD$8^{+}$ T Cell Responses During and After Acute Lymphocytic Choriomeningitis Virus Infection

A mathematical study of CD8+ T cell responses calibrated with human data

Dynamical Modeling and Qualitative Analysis of a Delayed Model for CD8 T Cells in Response to Viral Antigens

Mathematical Modeling to Guide Experimental Design: T Cell Clustering as a Case Study

IL-2 Sensitivity and Exogenous IL-2 Concentration Gradient Tune the Productive Contact Duration of CD8 + T Cell-Apc: a Multiscale Modeling Study

Dynamics of Activation and Regulation of the Immune Response to Attack by Viral Pathogens Using Mathematical Modeling

Functional Characterization of Ly49<SUP>+</SUP>CD8 T-Cells in Both Normal Condition and During Anti-Viral Response

Dynamical Analysis of Antigen-Driven T-cell Infection Model with Multiple Delays

Measles Infection Dose Responses: Insights from Mathematical Modeling

Analyzing Mouse B Cell Responses Specific to LCMV Infection.

Novel Lymphocytic Choriomeningitis Virus Strain Sustains Abundant Exhausted Progenitor CD8 T Cells without Systemic Viremia

Modelling cross-reactivity and memory in the cellular adaptive immune response to influenza infection in the host

Mathematical modeling of antigenicity for HIV dynamics

Modeling interleukin-2-based immunotherapy in AIDS pathogenesis

A statistical analysis of memory CD8 T cell differentiation: An application of a hierarchical state space model to a short time course microarray experiment

Accessing Early Differentiation of Virus-Specific Follicular Helper CD4+ T Cell in Acute LCMV-Infected Mice

The Impact of the Cd8+Cell Non-Cytotoxic Antiviral Response (cnar) and Cytotoxic T Lymphocyte (ctl) Activity in A Cell-To-Cell Spread Model for Hiv-1 with A Time Delay

Immunodominance Analysis Through Interactions of CD8+ T Cells and DCs in Lymph Nodes

Dynamic Modeling of Antibody Repertoire Reshaping in Response to Viral Infections

Viral dynamics model with CTL immune response incorporating antiretroviral therapy

In silico model of infection of a CD4(+) T-cell by a human immunodeficiency type 1 virus, and a mini-review on its molecular pathophysiology