Chimeric Antigen Receptor (CAR)-T cell immunotherapy revolutionized cancer treatment and consists of the genetic modification of T lymphocytes with a CAR gene, aiming to increase their ability to recognize and kill antigen-specific tumor cells. The dynamics of CAR-T cell responses in patients presents a multiphasic kinetics with distribution, expansion, contraction, and persistence phases. The characteristics and duration of each phase depend on the tumor type, the infused product, and on patient-specific characteristics. We present a mathematical model which describes the multiphasic CAR-T cell dynamics resulting from the interplay between CAR-T and tumor cells, considering patient and product heterogeneities. The CAR-T cell population is divided into functional (distributed and effector), memory, and exhausted CAR-T cell phenotypes. The model is able to describe the diversity of CAR-T cell dynamic behaviors in different patients and hematological cancers as well as their therapy outcomes. Our results indicate that the joint assessment of the area under the concentration-time curve in the first 28 days and the corresponding fraction of non-exhausted CAR-T cells may be considered as potential markers to classify therapy responses. Overall, the analysis of different CAR-T cell phenotypes can be a key aspect for a better understanding of the whole CAR-T cell dynamics.

In this manuscript we use an exactly solvable stochastic binary model for regulation of gene expression to analyse the dynamics of response to a treatment aiming to modulate the number of transcripts of RKIP gene. We demonstrate the usefulness of our method simulating three treatment scenarios aiming to reestablish RKIP gene expression dynamics towards pre-cancerous state: i. to increase the promoter’s ON state duration; ii. to increase the mRNAs’ synthesis rate; iii. to increase both rates. We show that the pre-treatment kinetic rates of ON and OFF promoter switching speeds and mRNA synthesis and degradation will affect the heterogeneity and time for treatment response. Hence, we present a strategy for reducing drug dosage by simultaneously targeting multiple kinetic rates. That enables a reduction of treatment response time and heterogeneity which in principle diminishes the chances of emergence of resistance to treatment. This approach may be useful for inferring kinetic constants related to expression of antimetastatic genes or oncogenes and on the design of multi-drug therapeutic strategies targeting master regulatory genes.

Immunotherapy has gained great momentum with chimeric antigen receptor T cell (CAR-T) therapy, in which patient’s T lymphocytes are genetically manipulated to recognize tumor-specific antigens increasing tumor elimination efficiency. In the last years, CAR-T cell immunotherapy for hematological malignancies achieved a great response rate on patients and is a very promising therapy for several other malignancies. Each new CAR design requires a preclinical proof-of-concept experiment using immunodeficient mouse models. The absence of a functional immune system in these mice makes them simple and suitable to be mathematically modeled. In this work, we developed a three population mathematical model to describe tumor response to CAR-T cell immunotherapy in immunodeficient mouse models, encompassing interactions between a non-solid tumor and CAR-T cells (effector and long-term memory). We account for several phenomena, such as tumor-induced immunosuppression, memory pool formation, and conversion of memory into effector CAR-T cells in the presence of new tumor cells. Individual donor and tumor specificities were considered as uncertainties in the model parameters. Our model is able to reproduce several CAR-T cell immunotherapy scenarios, with different CAR receptors and tumor targets reported in the literature. We found that therapy effectiveness mostly depends on some specific parameters such as the differentiation of effector to memory CAR-T cells, CAR-T cytotoxic capacity, tumor growth rate, and tumor-induced immunosuppression. In summary, our model can contribute to reduce and optimize the number of in vivo experiments with in silico tests to select specific scenarios that could be tested in experimental research. Such in silico laboratory was made available in a Shiny R-based platform called CARTmath. It is an open-source, easy to run simulator, available at github.com/tmglncc/CARTmath or directly on the webpage cartmath.lncc.br, containing this manuscript results as examples and documentation. The developed model, together with the CARTmath platform, provides potential use for assessing different CAR-T cell immunotherapy protocols and associated efficacy, becoming an accessory towards in silico trials.