Mathematical modeling of endogenous and exogenously administered T cell recirculation in mouse and its application to pharmacokinetic studies of cell therapies.

Introduction: In vivo T cell migration has been of interest to scientists for the past 60 years. T cell kinetics are important in the understanding of the immune response to infectious agents. More recently, adoptive T cell therapies have proven to be a most promising approach to treating a wide range of diseases, including autoimmune and cancer diseases, whereby the characterization of cellular kinetics represents an important step towards the prediction of therapeutic efficacy. Methods: Here, we developed a physiologically-based pharmacokinetic (PBPK) model that describes endogenous T cell homeostasis and the kinetics of exogenously administered T cells in mouse. Parameter calibration was performed using a nonlinear fixed-effects modeling approach based on published data on T cell kinetics and steady-state levels in different tissues of mice. The Partial Rank Correlation Coefficient (PRCC) method was used to perform a global sensitivity assessment. To estimate the impact of kinetic parameters on exogenously administered T cell dynamics, a local sensitivity analysis was conducted. Results: We simulated the model to analyze cellular kinetics following various T cell doses and frequencies of CCR7+ T cells in the population of infused lymphocytes. The model predicted the effects of T cell numbers and of population composition of infused T cells on the resultant concentration of T cells in various organs. For example, a higher percentage of CCR7+ T cells among exogenously administered T lymphocytes led to an augmented accumulation of T cells in the spleen. The model predicted a linear dependence of T cell dynamics on the dose of adoptively transferred T cells. Discussion: The mathematical model of T cell migration presented here can be integrated into a multi-scale model of the immune system and be used in a preclinical setting for predicting the distribution of genetically modified T lymphocytes in various organs, following adoptive T cell therapies.

Mathematical modeling in autoimmune diseases: from theory to clinical application.

The research & development (R&D) of novel therapeutic agents for the treatment of autoimmune diseases is challenged by highly complex pathogenesis and multiple etiologies of these conditions. The number of targeted therapies available on the market is limited, whereas the prevalence of autoimmune conditions in the global population continues to rise. Mathematical modeling of biological systems is an essential tool which may be applied in support of decision-making across R&D drug programs to improve the probability of success in the development of novel medicines. Over the past decades, multiple models of autoimmune diseases have been developed. Models differ in the spectra of quantitative data used in their development and mathematical methods, as well as in the level of "mechanistic granularity" chosen to describe the underlying biology. Yet, all models strive towards the same goal: to quantitatively describe various aspects of the immune response. The aim of this review was to conduct a systematic review and analysis of mathematical models of autoimmune diseases focused on the mechanistic description of the immune system, to consolidate existing quantitative knowledge on autoimmune processes, and to outline potential directions of interest for future model-based analyses. Following a systematic literature review, 38 models describing the onset, progression, and/or the effect of treatment in 13 systemic and organ-specific autoimmune conditions were identified, most models developed for inflammatory bowel disease, multiple sclerosis, and lupus (5 models each). ≥70% of the models were developed as nonlinear systems of ordinary differential equations, others - as partial differential equations, integro-differential equations, Boolean networks, or probabilistic models. Despite covering a relatively wide range of diseases, most models described the same components of the immune system, such as T-cell response, cytokine influence, or the involvement of macrophages in autoimmune processes. All models were thoroughly analyzed with an emphasis on assumptions, limitations, and their potential applications in the development of novel medicines.

Genotypes of the UCP1 gene polymorphisms and cardiometabolic diseases: A multifactorial study of association with disease probability.

Cardiometabolic diseases (CMDs) are complex disorders with a heterogenous phenotype, which are caused by multiple factors including genetic factors. Single nucleotide polymorphisms (SNPs) rs45539933 (p.Ala64Thr), rs10011540 (c.-112A>C), rs3811791 (c.-1766A>G), and rs1800592 (c.-3826A>G) in the UCP1 gene have been analyzed for association with CMDs in many studies providing controversial results. However, previous studies only considered individual UCP1 SNPs and did not evaluate them in an integrated manner, which is a more powerful approach to uncover genetic component of complex diseases. This study aimed to investigate associations between UCP1 genotype combinations and CMDs or CMD risk factors in the context of non-genetic factors. We performed multiple logistic regression analysis and proposed new methodology of testing different combinations of SNP genotypes. We found that probability of CMDs increased in presence of the three-SNP combination of genotypes with minor alleles of c.-3826A>G and p.Ala64Thr and wild allele of c.-112A>C, with increasing age, body mass index (BMI), body fat percentage (BF%) and may differ between sexes and between countries. The combination of genotypes with c.-3826A>G minor allele and wild homozygotes of c.-112A>C and p.Ala64Thr was associated with increased probability of diabetes. While combination of genotypes with minor alleles of all three SNPs reduced the CMD probability. The present results suggest that age, BMI, sex, and UCP1 three-SNP combinations of genotypes significantly contribute to CMD probability. Varying of c.-112A>C alleles in the genotype combination with minor alleles of c.-3826A>G and p.Ala64Thr markedly changes CMD probability.