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 statistical analysis of memory CD8 T cell differentiation: An application of a hierarchical state space model to a short time course microarray experiment

Model-based assessment of the Role of Uneven Partitioning of Molecular Content on Heterogeneity and Regulation of Differentiation in CD8 T-cell Immune Responses

Inference of Spatiotemporal Effects on Cellular State Transitions from Time-Lapse Microscopy

A hidden spatial-temporal Markov random field model for network-based analysis of time course gene expression data

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

A variational deep-learning approach to modeling memory T cell dynamics

State-space modeling of dynamic genetic networks

A data-driven Boolean model explains memory subsets and evolution in CD8+ T cell exhaustion

Functional Diversity of Memory CD8 T Cells is Spatiotemporally Imprinted

modelling of CD8 T cell immune response links genetic regulation to population dynamics

Genetic analysis of differentiation of T-helper lymphocytes

Computational investigations into the orgins of 'short term' biochemical memory in T cell activation

Time-variant Clustering Model for Understanding Cell Fate Decisions.

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

Hierarchical Control of State Transitions in Dense Associative Memories

Bayesian analysis of immune response dynamics with sparse time series data

Methods for Time Series Analysis of RNA-seq Data with Application to Human Th17 Cell Differentiation.

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

Jointly Modeling Longitudinal and Event Time Data With Application to Acquired Immunodeficiency Syndrome

Network estimation in State Space Model with L1-regularization constraint

The architectural design of CD8+ T cell responses in acute and chronic infection: Parallel structures with divergent fates