A mathematical model is proposed to study the complicated dynamics of tumor-immune interactions with three discrete time delays. The proposed model consists of nine coupled ordinary differential equations (ODEs) whose components are tumor cells, tumor-specific CD8+T cells, dendritic cells, macrophages, regulatory T-cells (Tregs), interleukin-10 (IL-10), IL-12, transforming growth factor-β (TGF-β) and interferon-γ (IFN-γ). By utilizing the quasi-steady-state approximation for cytokine concentrations, we obtain a system of four coupled ODEs. We introduce three discrete time delays into our deterministic system to better understand the tumor-immune dynamics. Basic properties of the system, including existence, uniqueness, positivity, boundedness, and uniform persistence are discussed. Stability analysis of both the delayed and non-delayed system has been performed and the conditions for stability and direction of Hopf bifurcation have been determined. Furthermore, we assess the length of time delay required to maintain the stability of period-1 limit cycle. Parameter estimation techniques are discussed and numerical simulations are provided to support our theoretical analysis. Notably, we observe that time delays do not significantly influence the system behavior for the existing set of parameters value.
Read full abstract