MCMC Methods for Parameter Estimation in ODE Systems for CAR-T Cell Cancer Therapy

Chimeric antigen receptor (CAR)-T cell therapy represents a breakthrough in treating resistant hematologic cancers. It is based on genetically modifying T cells transferred from the patient or a donor. Although its implementation has increased over the last few years, CAR-T has many challenges to be addressed, for instance, the associated severe toxicities, such as cytokine release syndrome. To model CAR-T cell dynamics, focusing on their proliferation and cytotoxic activity, we developed a mathematical framework using ordinary differential equations (ODEs) with Bayesian parameter estimation. Bayesian statistics were used to estimate model parameters through Monte Carlo integration, Bayesian inference, and Markov chain Monte Carlo (MCMC) methods. This paper explores MCMC methods, including the Metropolis-Hastings algorithm and DEMetropolis and DEMetropolisZ algorithms, which integrate differential evolution to enhance convergence rates. The theoretical findings and algorithms were validated using Python and Jupyter Notebooks. A real medical dataset of CAR-T cell therapy was analyzed, employing optimization algorithms to fit the mathematical model to the data, with the PyMC library facilitating Bayesian analysis. The results demonstrated that our model accurately captured the key dynamics of CAR-T cell therapy. This conclusion underscores the potential of parameter estimation to improve the understanding and effectiveness of CAR-T cell therapy in clinical settings.