Mathematical analysis of stochastic models for tumor-immune systems

In this paper we investigate some stochastic models for tumor-immune systems. To describe these models we used a Wiener process, as the noise has a stabilization effect. Their dynamics are studied in terms of stochastic stability in the equilibrium points, by constructing the Lyapunov exponent, depending on the parameters that describe the model. We have studied and analyzed a Kuznetsov-Taylor like stochastic model and a Bell stochastic model for tumor-immune systems. These stochastic models are studied from stability point of view and they were represented using the Euler second order scheme.

💡 Research Summary

The paper addresses the inherent stochasticity of tumor‑immune interactions by formulating two stochastic differential equation (SDE) models that incorporate Wiener‑process noise. Traditional deterministic models fail to capture random fluctuations arising from environmental variability and intracellular signaling noise; the authors therefore augment the classic tumor‑immune dynamics with additive Brownian motion terms. The first model, a Kuznetsov‑Taylor‑type formulation, represents immune activation as a nonlinear function of tumor burden, while the second, a Bell‑type model, includes a probabilistic transition of immune cells upon antigen recognition. Both models consist of coupled equations for tumor cell density (x(t)) and immune cell density (y(t)), with parameters governing growth, death, and interaction rates, plus a noise intensity (\sigma).

Stability analysis proceeds by linearizing the SDEs around equilibrium points ((x^,y^)) and applying Itô calculus to derive a Lyapunov exponent. A negative exponent indicates stochastic asymptotic stability. By sweeping the parameter space, the authors demonstrate that increasing (\sigma) can shift the Lyapunov exponent further into the negative region, revealing a noise‑induced stabilization effect that is absent in the deterministic counterpart. This effect is especially pronounced when the immune response parameter (c) (or the transition rate (k) in the Bell model) exceeds a critical threshold, suggesting that moderate stochastic perturbations may enhance immune surveillance.

Numerical simulations employ a second‑order Euler‑Maruyama scheme, which offers (O(\Delta t^{2})) strong convergence and thus higher fidelity in regimes where noise dominates the dynamics. Simulations across diverse initial conditions illustrate two qualitatively distinct long‑term behaviors: (1) a tumor‑suppressed state where immune cells maintain control and tumor size remains low, and (2) a tumor‑progressive state where immune activity is insufficient and tumor cells proliferate. As the noise intensity grows, trajectories increasingly converge to the suppressed state, confirming the analytical predictions.

The contributions of the work are threefold. First, it provides a rigorous stochastic framework for tumor‑immune dynamics, explicitly modeling environmental and intracellular randomness. Second, it introduces a Lyapunov‑based criterion for stochastic stability, enabling quantitative assessment of how model parameters and noise intensity affect system behavior. Third, it validates the theoretical findings with high‑order numerical experiments, demonstrating that noise can play a constructive role in tumor control. The authors argue that these insights have practical implications for immunotherapy design, where controlled stochastic perturbations (e.g., dosing variability or engineered cytokine fluctuations) might be leveraged to improve treatment outcomes. Future directions include extending the models to incorporate cytokine networks, spatial heterogeneity, and patient‑specific data for personalized predictive simulations.