Modeling the Nonmonotonic Immune Response in a Tumor–Immune System Interaction

Abstract

Tumor–immune system interactions are very complicated, being highly nonlinear and not well understood. A large number of tumors can potentially weaken the immune system through various mechanisms such as secreting cytokines that suppress the immune response. In this paper, we propose a tumor–immune system interaction model with a nonmonotonic immune response function and adoptive cellular immunotherapy (ACI). The model has a tumor-free equilibrium and at most three tumor-presence equilibria (low, moderate and high ones). The stability of all equilibria is studied by analyzing their characteristic equations. The consideration of nonmonotonic immune response results in a series of bifurcations such as the saddle-node bifurcation, transcritical bifurcation, Hopf bifurcation and Bogdanov–Takens bifurcation. In addition, numerical simulation results show the coexistence of periodic orbits and homoclinic orbits. Interestingly, along with various bifurcations, we also found two bistable scenarios: the coexistence of a stable tumor-free as well as a high-tumor-presence equilibrium and the coexistence of a stable-low as well as a high-tumor-presence equilibrium, which can show symmetric and antisymmetric properties in a range of model parameters and initial cell concentrations. The new findings indicate that under ACI, patients can possibly reach either a stable tumor-free state or a low-tumor-presence state in the presence of nonmonotonic immune response once the immune system is activated.