Mathematical Modelling for the Role of CD4+T Cells in Tumor-Immune Interactions

Mathematical modelling has been used to study tumor-immune cell interaction. Some models were proposed to examine the effect of circulating lymphocytes, natural killer cells, and CD8+T cells, but they neglected the role of CD4+T cells. Other models were constructed to study the role of CD4+T cells but did not consider the role of other immune cells. In this study, we propose a mathematical model, in the form of a system of nonlinear ordinary differential equations, that predicts the interaction between tumor cells and natural killer cells, CD4+T cells, CD8+T cells, and circulating lymphocytes with or without immunotherapy and/or chemotherapy. This system is stiff, and the Runge–Kutta method failed to solve it. Consequently, the “Adams predictor-corrector” method is used. The results reveal that the patient's immune system can overcome small tumors; however, if the tumor is large, adoptive therapy with CD4+T cells can be an alternative to both CD8+T cell therapy and cytokines in some cases. Moreover, CD4+T cell therapy could replace chemotherapy depending upon tumor size. Even if a combination of chemotherapy and immunotherapy is necessary, using CD4+T cell therapy can better reduce the dose of the associated chemotherapy compared to using combined CD8+T cells and cytokine therapy. Stability analysis is performed for the studied patients. It has been found that all equilibrium points are unstable, and a condition for preventing tumor recurrence after treatment has been deduced. Finally, a bifurcation analysis is performed to study the effect of varying system parameters on the stability, and bifurcation points are specified. New equilibrium points are created or demolished at some bifurcation points, and stability is changed at some others. Hence, for systems turning to be stable, tumors can be eradicated without the possibility of recurrence. The proposed mathematical model provides a valuable tool for designing patients' treatment intervention strategies.

[1]  E. Jaffee,et al.  Role of bone marrow-derived cells in presenting MHC class I-restricted tumor antigens. , 1994, Science.

[2]  D. Jordan,et al.  Nonlinear Ordinary Differential Equations: An Introduction for Scientists and Engineers , 1979 .

[3]  Junde Wu,et al.  Stability analysis of a multi-layer tumor model with free boundary , 2019, Comput. Math. Appl..

[4]  Desmond J. Higham,et al.  Numerical Methods for Ordinary Differential Equations - Initial Value Problems , 2010, Springer undergraduate mathematics series.

[5]  C. N. Brorn,et al.  WHO? , 1896 .

[6]  S. Sastry.  Nonlinear Systems: Analysis, Stability, and Control , 1999 .

[7]  K. Swanson,et al.  A mathematical model for brain tumor response to radiation therapy , 2009, Journal of mathematical biology.

[8]  Jacinto Javier Ibáñez,et al.  Adams-Bashforth and Adams-Moulton methods for solving differential Riccati equations , 2010, Comput. Math. Appl..

[9]  Ricardo Gómez Aiza,et al.  Analysis of a mathematical model of periodically pulsed chemotherapy treatment , 2017 .

[10]  Padmanabhan Seshaiyer,et al.  Mathematical Modeling, Analysis, and Simulation of Tumor Dynamics with Drug Interventions , 2019, Comput. Math. Methods Medicine.

[11]  Antonio Lanzavecchia,et al.  T cell fitness determined by signal strength , 2003, Nature Immunology.

[12]  D. Asemani,et al.  MATHEMATICAL MODELING OF IN-VIVO TUMOR-IMMUNE INTERACTIONS FOR THE CANCER IMMUNOTHERAPY USING MATURED DENDRITIC CELLS , 2018 .

[13]  Junseok Kim,et al.  Mathematical model and its fast numerical method for the tumor growth. , 2015, Mathematical biosciences and engineering : MBE.

[14]  Kazuyuki Aihara,et al.  A Mathematical Model of Prostate Tumor Growth Under Hormone Therapy with Mutation Inhibitor , 2010, J. Nonlinear Sci..

[15]  Yinnian He,et al.  Boundary element method for a free boundary problem modeling three dimensional tumor growth , 2017, Comput. Math. Appl..

[16]  D. Rowley,et al.  NK cells suppress the generation of Lyt-2+ cytolytic T cells by suppressing or eliminating dendritic cells. , 1986, Journal of immunology.

[17]  Ali Ghaffari,et al.  Suboptimal mixed vaccine and chemotherapy in finite duration cancer treatment: state-dependent Riccati equation control , 2015 .

[18]  A. Radunskaya,et al.  Mathematical Model of Colorectal Cancer with Monoclonal Antibody Treatments , 2013, 1312.3023.

[19]  A. Perelson,et al.  Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis. , 1994, Bulletin of mathematical biology.

[20]  Mustafa Mamat,et al.  Mathematical Model of Cancer Treatments Using Immunotherapy, Chemotherapy and Biochemotherapy , 2013 .

[21]  M. Nazari,et al.  A mixed radiotherapy and chemotherapy model for treatment of cancer with metastasis , 2016 .

[22]  Sara Jabbari,et al.  Mathematical model predicts anti-adhesion–antibiotic–debridement combination therapies can clear an antibiotic resistant infection , 2019, PLoS Comput. Biol..

[23]  H. Frieboes,et al.  Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method. , 2008, Journal of theoretical biology.

[24]  R. Glas,et al.  Recruitment and Activation of Natural Killer (Nk) Cells in Vivo Determined by the Target Cell Phenotype , 2000, The Journal of experimental medicine.

[25]  A. Radunskaya,et al.  Mixed Immunotherapy and Chemotherapy of Tumors: Modeling, Applications and Biological Interpretations , 2022 .

[26]  Damián R. Fernández,et al.  Adjoint method for a tumor growth PDE-constrained optimization problem , 2012, Comput. Math. Appl..

[27]  Kathleen P Wilkie,et al.  A review of mathematical models of cancer-immune interactions in the context of tumor dormancy. , 2013, Advances in experimental medicine and biology.

[28]  Jinchao Xu,et al.  A mathematical model of aortic aneurysm formation , 2017, PloS one.

[29]  O. Olopade,et al.  Treatment of cancer in sub-Saharan Africa. , 2013, The Lancet. Oncology.

[30]  Alain Goriely,et al.  A mathematical model of tumor-immune interactions. , 2012, Journal of theoretical biology.

[31]  Ami Radunskaya,et al.  A delay differential equation model for tumor growth , 2003, Journal of mathematical biology.

[32]  Xiaochuan Hu,et al.  Dynamics of tumor-CD4-cytokine-host cells interactions with treatments , 2018, Appl. Math. Comput..

[33]  Jui-Ling Yu,et al.  Qualitative behavior of systems of tumor–CD4+–cytokine interactions with treatments , 2015 .

[34]  D. Earn,et al.  Interactions Between the Immune System and Cancer: A Brief Review of Non-spatial Mathematical Models , 2011, Bulletin of mathematical biology.

[35]  G. Samanta,et al.  Model Justification and Stratification for Confounding of Chlamydia Trachomatis Disease , 2019, Letters in Biomathematics.

[36]  Lin Shen,et al.  Mathematical Modelling and Analysis of the Tumor Treatment Regimens with Pulsed Immunotherapy and Chemotherapy , 2016, Comput. Math. Methods Medicine.

[37]  D. Kirschner,et al.  Modeling immunotherapy of the tumor – immune interaction , 1998, Journal of mathematical biology.

[38]  Ami Radunskaya,et al.  Modeling Tumor–Immune Dynamics , 2014 .

[39]  M. Caligiuri,et al.  The biology of human natural killer-cell subsets. , 2001, Trends in immunology.

[40]  Rachid Ouifki,et al.  Mathematical model of tumor-immune surveillance. , 2016, Journal of theoretical biology.

[41]  O. Lantz,et al.  CD4 cells can be more efficient at tumor rejection than CD8 cells. , 2005, Blood.

[42]  K. Renee Fister,et al.  Mathematical model creation for cancer chemo-immunotherapy , 2009 .

[43]  O. M. Ogunlaran,et al.  Mathematical Model for an Effective Management of HIV Infection , 2016, BioMed research international.

[44]  Joseph Malinzi,et al.  Mathematical Analysis of a Mathematical Model of Chemovirotherapy: Effect of Drug Infusion Method , 2019, Comput. Math. Methods Medicine.

[45]  D. Kirschner,et al.  A mathematical model of tumor-immune evasion and siRNA treatment , 2003 .

[46]  D. Earn,et al.  Modeling anti-tumor Th1 and Th2 immunity in the rejection of melanoma. , 2010, Journal of theoretical biology.

[47]  L. G. de Pillis,et al.  OPTIMAL CONTROL OF MIXED IMMUNOTHERAPY AND CHEMOTHERAPY OF TUMORS , 2008 .

[48]  G. P. Samanta,et al.  Dynamical Behaviour of a Tumor-Immune System with Chemotherapy and Optimal Control , 2013 .

[49]  M. Kaplan,et al.  CD4+ T‐cell‐mediated anti‐tumor immunity can be uncoupled from autoimmunity via the STAT4/STAT6 signaling axis , 2009, European journal of immunology.

[50]  V. Young,et al.  Adoptive cell therapy with CD4+ T helper 1 cells and CD8+ cytotoxic T cells enhances complete rejection of an established tumour, leading to generation of endogenous memory responses to non-targeted tumour epitopes , 2017, Clinical & translational immunology.

[51]  Kazuyuki Aihara,et al.  A partial differential equation model and its reduction to an ordinary differential equation model for prostate tumor growth under intermittent hormone therapy , 2014, Journal of mathematical biology.

[52]  S. Shabahang,et al.  A Mathematical Model for DC Vaccine Treatment of Type I Diabetes , 2019, Front. Physiol..

[53]  P. Foster,et al.  Immunotherapy of Cytotoxic T Cell–resistant Tumors by T Helper 2 Cells , 2003, The Journal of experimental medicine.

[54]  A. Cesano,et al.  Bringing the Next Generation of Immuno-Oncology Biomarkers to the Clinic , 2018, Biomedicines.

[55]  Lieping Chen,et al.  A Paradigm Shift in Cancer Immunotherapy: From Enhancement to Normalization , 2018, Cell.

[56]  Philipp M. Altrock,et al.  The mathematics of cancer: integrating quantitative models , 2015, Nature Reviews Cancer.

[57]  D. G. Mallet,et al.  Spatial tumor-immune modeling , 2006 .

[58]  I. Melero,et al.  Cytokines in clinical cancer immunotherapy , 2018, British Journal of Cancer.

[59]  G. Samanta,et al.  Analysis of the Dynamics of a Tumor–Immune System with Chemotherapy and Immunotherapy and Quadratic Optimal Control , 2016 .

[60]  Yasmine Abouelseoud,et al.  Innovative Approximate Analytical Solution for Standard Model of Viral Dynamics: Hepatitis C with Direct-Acting Agents as an Implemented Case , 2019 .

[61]  M. Raffeld,et al.  Cancer Regression and Autoimmunity in Patients After Clonal Repopulation with Antitumor Lymphocytes , 2002, Science.

[62]  Mark A. Franklin,et al.  Optimally stable parallel predictors for Adams—Moulton correctors , 1977 .

[64]  Manuela Herman,et al.  An Introduction To Numerical Methods For Differential Equations , 2016 .

[65]  Existence and uniqueness of global solutions for a mathematical model of antiangiogenesis in tumor growth , 2008 .

[66]  Vittorio Cristini,et al.  Three-dimensional multispecies nonlinear tumor growth-II: Tumor invasion and angiogenesis. , 2010, Journal of theoretical biology.

[67]  Raluca Eftimie,et al.  Mathematical Models for Immunology: Current State of the Art and Future Research Directions , 2016, Bulletin of mathematical biology.

[68]  Oliver Lundqvist,et al.  Numerical Methods for Ordinary Differential Equations , 2013, An Introduction to Numerical Methods and Analysis 3e.

[69]  Michelle Ghrist,et al.  Stability ordinates of Adams predictor-corrector methods , 2015 .