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 recent years, CAR-T cell immunotherapy for hematological malignancies achieved a great response rate in 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 for use as mathematical models. In this work, we develop 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 are considered 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 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 reducing and optimizing the number of in vivo experiments with in silico tests to select specific scenarios that could be tested in experimental research. Such an in silico laboratory is an easy-to-run open-source simulator, built on a Shiny R-based platform called CARTmath. It contains the results of this manuscript as examples and documentation. The developed model together with the CARTmath platform have potential use in assessing different CAR-T cell immunotherapy protocols and its associated efficacy, becoming an accessory for in silico trials.

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.

Mathematical and computational modeling have been increasingly applied in many areas of cancer research, aiming to improve the understanding of tumorigenic mechanisms and to suggest more effective therapy protocols. The mathematical description of the tumor growth dynamics is often made using the exponential, logistic, and Gompertz models. However, recent literature has suggested that the Allee effect may play an important role in the early stages of tumor dynamics, including cancer relapse and metastasis. For a model to provide reliable predictions, it is necessary to have a rigorous evaluation of the uncertainty inherent in the modeling process. In this work, our main objective is to show how a model framework that integrates sensitivity analysis, model calibration, and model selection techniques can improve and systematically characterize model and data uncertainties. We investigate five distinct models with different complexities, which encompass the exponential, logistic, Gompertz, and weak and strong Allee effects dynamics. Using tumor growth data published in the literature, we perform a global sensitivity analysis, apply a Bayesian framework for parameter inference, evaluate the associated sensitivity matrices, and use different information criteria for model selection (First- and Second-Order Akaike Information Criteria and Bayesian Information Criterion). We show that such a wider methodology allows having a more detailed picture of each model assumption and uncertainty, calibration reliability, ultimately improving tumor mathematical description. The used in vivo data suggested the existence of both a competitive effect among tumor cells and a weak Allee effect in the growth dynamics. The proposed model framework highlights the need for more detailed experimental studies on the influence of the Allee effect on the analyzed cancer scenario.

We investigate the feasibility of deriving surrogate ordinary differential equations (ODEs) from individual-based models (IBMs) simulated in the PhysiCell platform. As a preliminary demonstration, we focus on data generated from simulations of several cell cycle models, each with varying complexities and corresponding ODE representations. Our approach begins with an evaluation of the stochastic characteristics, computational requirements, and accuracy of PhysiCell simulations. These simulations heavily depend on parameters such as the time step governing cell phase transitions, initial cell population, and the number of replicates. We then demonstrate the application of a sparse identification machine learning tool, SINDy-SA, to construct surrogate ODE-based models for IBMs. The SINDy-SA framework effectively identified the mathematical structures for all cell cycle models used in our study. Given the ability of surrogate models to reduce computational costs while maintaining accuracy, we believe that the insights from this study will aid researchers in modeling hybrid discrete-continuous multiscale systems using the PhysiCell platform and in developing corresponding surrogate models through the SINDy-SA framework.