In this paper, we discuss a class of dynamic models of the interaction between tumors and the immune system with antigen delay and Michaelis-Menten type inhibition terms. The Michaelis-Menten type function
and
is used to describe the immune response of effector cells interacting with tumor cells, the linear antigen stimulation term
is used to describe the linear recruitment effect of tumor antigens on effector cells, and the stimulation delay of tumor antigens in the immune system is introduced. The dynamic behavior of the model is studied through qualitative analysis and numerical simulation. Saddle-node bifurcation may occur both in the case with and without time delay. Contrary to the case without time delay, stimulation delay may lead to some complex dynamic behaviors and biological phenomena. In the presence of time delay, the existence condition of Hopf bifurcation at the equilibrium point is obtained. Further discussion shows that the model may exhibit bistability under certain conditions, that is, the growth and development state of the tumor depends on its initial state. Finally, numerical simulation is used to verify the accuracy of the relevant theoretical results, and the corresponding biological significance is briefly discussed.
References
[1]
Wenbo, L. and Wang, J. (2017) Uncovering the Underlying Mechanism of Cancer Tumorigenesis and Development under an Immune Microenvironment from Global Quantification of the Landscape. JournaloftheRoyalSocietyInterface, 14, Article ID: 20170105. https://doi.org/10.1098/rsif.2017.0105
[2]
Schreiber, R.D., Old, L.J. and Smyth, M.J. (2011) Cancer Immunoediting: Integrating Immunity’s Roles in Cancer Suppression and Promotion. Science, 331, 1565-1570. https://doi.org/10.1126/science.1203486
[3]
Sung, W., Hong, T.S., Poznansky, M.C., Paganetti, H. and Grassberger, C. (2022) Mathematical Modeling to Simulate the Effect of Adding Radiation Therapy to Immunotherapy and Application to Hepatocellular Carcinoma. InternationalJournalofRadiationOncology*Biology*Physics, 112, 1055-1062. https://doi.org/10.1016/j.ijrobp.2021.11.008
[4]
Zhang, Z., Li, S., Si, P., Li, X. and He, X. (2022) A Tumor-Immune Model with Mixed Immunotherapy and Chemotherapy: Qualitative Analysis and Optimal Control. JournalofBiologicalSystems, 30, 339-364. https://doi.org/10.1142/s0218339022500127
[5]
Ghaffari Laleh, N., Loeffler, C.M.L., Grajek, J., Staňková, K., Pearson, A.T., Muti, H.S., et al. (2022) Classical Mathematical Models for Prediction of Response to Chemotherapy and Immunotherapy. PLOSComputationalBiology, 18, e1009822. https://doi.org/10.1371/journal.pcbi.1009822
[6]
Pang, L., Shen, L. and Zhao, Z. (2016) Mathematical Modelling and Analysis of the Tumor Treatment Regimens with Pulsed Immunotherapy and Chemotherapy. ComputationalandMathematicalMethodsinMedicine, 2016, Article ID: 6260474. https://doi.org/10.1155/2016/6260474
[7]
Khajanchi, S. and Ghosh, D. (2015) The Combined Effects of Optimal Control in Cancer Remission. AppliedMathematicsandComputation, 271, 375-388. https://doi.org/10.1016/j.amc.2015.09.012
[8]
Ioachim, H.L. (1980) Correlations between Tumor Antigenicity, Malignant Potential, and Local Host Immune Response. In: Witz, I.P. and Hanna, M.G., Eds., InSituExpressionofTumorImmunity, Springer, 213-238. https://doi.org/10.1007/978-1-4684-3677-8_10
[9]
De Boer, R.J., Hogeweg, P., Dullens, H.F., De Weger, R.A. and Den Otter, W. (1985) Macrophage T Lymphocyte Interactions in the Anti-Tumor Immune Response: A Mathematical Model. TheJournalofImmunology, 134, 2748-2758. https://doi.org/10.4049/jimmunol.134.4.2748
[10]
Nishida, N. and Kudo, M. (2016) Immunological Microenvironment of Hepatocellular Carcinoma and Its Clinical Implication. Oncology, 92, 40-49. https://doi.org/10.1159/000451015
[11]
Kuznetsov, V., Makalkin, I., Taylor, M. and Perelson, A. (1994) Nonlinear Dynamics of Immunogenic Tumors: Parameter Estimation and Global Bifurcation Analysis. BulletinofMathematicalBiology, 56, 295-321. https://doi.org/10.1016/s0092-8240(05)80260-5
[12]
Kirschner, D. and Panetta, J.C. (1998) Modeling Immunotherapy of the Tumor—Immune Interaction. JournalofMathematicalBiology, 37, 235-252. https://doi.org/10.1007/s002850050127
[13]
Delisi, C. and Rescigno, A. (1977) Immune Surveillance and Neoplasia—1 A Minimal Mathematical Model. BulletinofMathematicalBiology, 39, 201-221. https://doi.org/10.1016/s0092-8240(77)80008-6
[14]
Adam, J.A. (1996) Effects of Vascularization on Lymphocyte/Tumor Cell Dynamics: Qualitative Features. MathematicalandComputerModelling, 23, 1-10. https://doi.org/10.1016/0895-7177(96)00016-7
[15]
Yang, J., Tang, S. and Cheke, R.A. (2015) Modelling Pulsed Immunotherapy of Tumour–Immune Interaction. MathematicsandComputersinSimulation, 109, 92-112. https://doi.org/10.1016/j.matcom.2014.09.001
[16]
Zhang, G. and Wang, X. (2020) Dynamic Analysis of Tumor-Immune System with Inhibitor Term of Michaelis-Menten Type. Science Technology and Engineering, 20, 7137-7144.
[17]
Gałach, M. (2003) Dynamics of the Tumor-Immune System Competition—The Effect of Time Delay. International Journal of Applied Mathematics Computer Science, 13, 395-406.
[18]
Li, J., Xie, X., Chen, Y. and Zhang, D. (2021) Complex Dynamics of a Tumor-Immune System with Antigenicity. AppliedMathematicsandComputation, 400, Article ID: 126052. https://doi.org/10.1016/j.amc.2021.126052
[19]
Bi, P. and Ruan, S. (2013) Bifurcations in Delay Differential Equations and Applications to Tumor and Immune System Interaction Models. SIAMJournalonAppliedDynamicalSystems, 12, 1847-1888. https://doi.org/10.1137/120887898
[20]
Ruan, S. (2021) Nonlinear Dynamics in Tumor-Immune System Interaction Models with Delays. Discrete&ContinuousDynamicalSystems—B, 26, 541-602. https://doi.org/10.3934/dcdsb.2020282
[21]
Bi, P., Ruan, S. and Zhang, X. (2014) Periodic and Chaotic Oscillations in a Tumor and Immune System Interaction Model with Three Delays. Chaos: AnInterdisciplinaryJournalofNonlinearScience, 24, Article ID: 023101. https://doi.org/10.1063/1.4870363
[22]
Li, J., Xie, X., Zhang, D., Li, J. and Lin, X. (2021) Qualitative Analysis of a Simple Tumor-Immune System with Time Delay of Tumor Action. Discrete&ContinuousDynamicalSystems—B, 26, 5227-5249. https://doi.org/10.3934/dcdsb.2020341
[23]
Li, J., Liu, F., Chen, Y. and Zhang, D. (2022) Dynamic Analysis of a Model on Tumor-Immune System with Regulation of PD-1/PD-L1 and Stimulation Delay of Tumor Antigen. QualitativeTheoryofDynamicalSystems, 21, Article No. 90. https://doi.org/10.1007/s12346-022-00627-5
[24]
Li, J., Ma, X., Chen, Y. and Zhang, D. (2022) Complex Dynamic Behaviors of a Tumor-Immune System with Two Delays in Tumor Actions. DiscreteandContinuousDynamicalSystems—B, 27, 7065-7087. https://doi.org/10.3934/dcdsb.2022033
[25]
Li, J., Chen, Y., Cao, H., Zhang, D. and Zhang, P. (2023) A Simple Model of Tumor-Immune Interaction: The Effect of Antigen Delay. InternationalJournalofBifurcationandChaos, 33, Article ID: 2350129. https://doi.org/10.1142/s0218127423501298
[26]
Niu, B., Guo, Y. and Du, Y. (2018) Hopf Bifurcation Induced by Delay Effect in a Diffusive Tumor-Immune System. InternationalJournalofBifurcationandChaos, 28, Article ID: 1850136. https://doi.org/10.1142/s0218127418501365
[27]
Yu, M., Dong, Y. and Takeuchi, Y. (2016) Dual Role of Delay Effects in a Tumour-Immune System. JournalofBiologicalDynamics, 11, 334-347. https://doi.org/10.1080/17513758.2016.1231347
[28]
Kemwoue, F.F., Deli, V., Edima, H.C., Mendimi, J.M., Gninzanlong, C.L., Dedzo, M.M., et al. (2022) Effects of Delay in a Biological Environment Subject to Tumor Dynamics. Chaos, Solitons&Fractals, 158, Article ID: 112022. https://doi.org/10.1016/j.chaos.2022.112022
[29]
Khajanchi, S. and Banerjee, S. (2014) Stability and Bifurcation Analysis of Delay Induced Tumor Immune Interaction Model. AppliedMathematicsandComputation, 248, 652-671. https://doi.org/10.1016/j.amc.2014.10.009
[30]
Verhulst. F. (1996) Nonlinear Differential Equations and Dynamical Systems. 2th Edition, Springer.