2025-11-23T07:58:16.677890

Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects

Amster, Rivera, Arredondo
The main purpose of this paper is to analyze the dynamics of the system of time-delay differential equations (DDEs) \begin{equation*} \begin{split} \dot{T}(t)&=T(t) f(t,T(t))-γE(t)T(t),\\ \dot{E}(t)&=σ+ \frac{pE(t)T(t-τ_1)}{g+a T(t-τ_1)}-\frac{mE(t)T(t-τ_2)}{g+a T(t-τ_2)}-ηE(t), \end{split} \end{equation*} where $T=T(t)$ and $E=E(t)$ represent the concentrations of tumor and effector cells at the time $t$. The coefficients $σ$, $μ$, $γ$, and $η$ are all positive, and $f(t, T)$ represents the relative growth rate of tumor cells, corresponding to a generalized logistic growth function that describes periodic time chemotherapeutic effects. The parameter $τ_1 \in \mathbb{R}_{\ge 0}$ is the response time delay of the immune system (mediated by effector cells) to an invasion of tumor cells, while $τ_2 \in \mathbb{R}_{\ge 0}$ represents the time delay of tumor cells in response to the appearance of effector cells.
academic

Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects

基本信息

  • 论文ID: 2510.11135
  • 标题: Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects
  • 作者: Pablo Amster, Andrés Rivera, John A. Arredondo
  • 分类: math.DS (Dynamical Systems)
  • 发表时间: 2025年10月13日
  • 论文链接: https://arxiv.org/abs/2510.11135

摘要

本文主要分析具有时滞的微分方程系统的动力学性质：

\dot{T}(t) &= T(t)f(t,T(t)) - \gamma E(t)T(t) \\ \dot{E}(t) &= \sigma + \frac{pE(t)T(t-\tau_1)}{g+aT(t-\tau_1)} - \frac{mE(t)T(t-\tau_2)}{g+aT(t-\tau_2)} - \eta E(t) \end{align}$$ 其中$T(t)$和$E(t)$分别表示时刻$t$肿瘤细胞和效应细胞的浓度；$f(t,T)$表示肿瘤细胞的相对增长率，对应描述周期性化疗效应的广义logistic增长函数；$\tau_1$是免疫系统对肿瘤细胞入侵的响应时滞，$\tau_2$是肿瘤细胞对效应细胞出现的响应时滞。 ## 研究背景与动机 ### 问题背景 1. **肿瘤-免疫系统相互作用**: 肿瘤是机体组织的异常增长，正常情况下免疫系统可以控制，但肿瘤的形成往往与免疫系统缺陷有关 2. **时滞效应的重要性**: 生物系统中普遍存在时滞现象，免疫系统对肿瘤的识别和响应需要时间，肿瘤对免疫攻击的反应也存在延迟 3. **化疗的周期性效应**: 化疗通常采用周期性给药方案，需要考虑其对肿瘤增长的时变影响 ### 研究动机 1. **现有模型的局限性**: 经典的肿瘤-免疫模型多采用logistic增长，但实验数据显示Gompertz模型拟合效果更好 2. **时滞参数的影响**: 需要深入分析不同时滞参数对系统动力学的影响 3. **周期解的存在性**: 探索在化疗周期性作用下系统周期解的存在条件 ## 核心贡献 1. **建立了广义Richards型增长的肿瘤-免疫时滞模型**，包含两个不同的时滞参数$\tau_1$和$\tau_2$ 2. **完整分析了平衡点的存在性和稳定性**，给出了不同参数条件下平衡点的分类 3. **确定了稳定性切换的临界时滞值**，提供了显式的稳定区间 4. **证明了周期解的存在性**，通过延拓方法建立了从自治系统到非自治系统的周期解分支 5. **给出了Hopf分岔的条件**，分析了系统从稳定平衡点向周期解的转换 ## 方法详解 ### 任务定义 研究具有双时滞的肿瘤-免疫竞争系统的动力学性质，包括： - 平衡点的存在性和稳定性分析 - 时滞参数对系统稳定性的影响 - 周期解的存在性证明 - Hopf分岔分析 ### 模型架构 #### 基本模型框架 系统(2)描述了肿瘤细胞T(t)和效应细胞E(t)的相互作用： **肿瘤细胞方程**： $$\dot{T}(t) = T(t)f(t,T(t)) - \gamma E(t)T(t)$$ - 第一项：肿瘤细胞的自然增长 - 第二项：效应细胞对肿瘤细胞的杀伤 **效应细胞方程**： $$\dot{E}(t) = \sigma + \frac{pE(t)T(t-\tau_1)}{g+aT(t-\tau_1)} - \frac{mE(t)T(t-\tau_2)}{g+aT(t-\tau_2)} - \eta E(t)$$ - $\sigma$：效应细胞的恒定输入 - 第二项：肿瘤细胞对效应细胞的激活（Holling II型功能反应） - 第三项：效应细胞的耗竭 - 第四项：效应细胞的自然死亡 #### 增长函数设计 采用Richards型增长函数： $$f(t,T) = r\left[\frac{(1-T^\beta)}{\beta} - b(t)\right]$$ 其中： - $\beta \in (0,1]$：形状参数 - $b(t)$：周期性化疗函数 - 当$\beta=1$时退化为logistic增长 - 当$\beta \to 0$时趋近于Gompertz增长 ### 技术创新点 1. **双时滞机制**：区分免疫激活时滞$\tau_1$和免疫耗竭时滞$\tau_2$，更准确反映生物过程 2. **广义增长模型**：通过参数$\beta$统一logistic和Gompertz增长模式 3. **参数缩放技术**：通过变换(14)简化参数分析 4. **稳定性切换曲线**：给出了$(\tau_1,\tau_2)$平面上的稳定性边界 ## 理论分析 ### 平衡点分析 #### 无瘤平衡点 对于所有参数值，系统总存在无瘤平衡点$(0,\sigma/\eta)$。 **稳定性条件**（定理10）： 定义$\Delta := \gamma\sigma - rb\eta$，则： - 若$\Delta < 0$：无瘤平衡点不稳定 - 若$\Delta = 0$：无瘤平衡点稳定 - 若$\Delta > 0$：无瘤平衡点渐近稳定 #### 共存平衡点 当$a=0$时，共存平衡点$(T^*,E^*)$满足方程(17)： $$\mu(T^{*\beta} - b)T^* + T^{*\beta} = b - \sigma$$ **存在性条件**（命题8）： - 条件†：$\sigma < b$且$\mu \leq 0$ - 条件‡：$\sigma \leq b$且$0 < \mu \leq \mu_c$ 其中$\mu_c := \left(\frac{1}{b}\left(\frac{1-\beta}{1+\beta}\right)^{\beta-1}\right)^{1/\beta}$ ### 稳定性分析 #### 小时滞情况 **定理12**：对于充分小的$a,\tau_1,\tau_2 \geq 0$： 1. 在条件†下，$(T^*,E^*)$局部渐近稳定 2. 在条件‡下，若满足额外条件则$(T^*,E^*)$稳定 #### 一般时滞情况 **定理14**：当$\tau_1 = \tau_2 = \tau$时，存在临界值$\tau_c$使得： - $\tau < \tau_c$：平衡点稳定 - $\tau > \tau_c$：平衡点失稳并发生Hopf分岔 临界值由以下方程确定： $$G(y) = y^4 + \left[(\lambda_1^*)^2 + (\lambda_2^*)^2\right]y^2 + (\lambda_1^*\lambda_2^*)^2 - N^2 = 0$$ ### 周期解存在性 **定理16**（延拓定理）：设$(T^*,E^*)$是$b \equiv b_0$时的可行平衡点，则存在常数$\phi^*,\tau_j^* > 0$和$r^* > 0$，使得当$\|\phi\|_\infty < \phi^*$，$\tau_1 < \tau_1^*$，$\tau_2 < \tau_2^*$时，系统(2)至少有一个$\omega$-周期解$(T,E)$满足$\|(T-T^*,E-E^*)\|_\infty < r^*$。 **证明思路**： 1. 验证线性化系统无非平凡周期解（引理4） 2. 应用隐函数定理和Fredholm择一定理 3. 建立从自治无时滞系统到非自治时滞系统的光滑分支 ## 数值验证 ### 参数设置 论文采用文献中的典型参数值： - $g \approx 10^6$（相对较大） - 缩放后$a \to a/g = O(10^{-6})$（相对较小） ### 稳定性切换曲线 **定理15**给出了$(\tau_1,\tau_2)$平面上的稳定性切换曲线： $$C_{s,k} = \{(\tau_1^{\pm}(y), \tau_2^{\pm}(y)) \in \mathbb{R}_+^2 : y \in I, s,k \in \mathbb{Z}\}$$ 其中$\tau_1^{\pm}$和$\tau_2^{\pm}$由方程(35)和(36)给出。 ## 相关工作 ### 肿瘤增长模型 1. **经典模型**：指数增长、logistic增长、Gompertz增长 2. **实验验证**：多项研究表明Gompertz模型拟合实验数据效果最佳 3. **参数相关性**：Gompertz模型的两个参数间存在显著相关性 ### 肿瘤-免疫相互作用模型 1. **Kuznetsov等(1994)**：建立了基础的肿瘤-免疫模型框架 2. **Li等(2021,2022)**：研究了单时滞和双时滞的影响 3. **本文贡献**：扩展到广义Richards增长和周期性化疗效应 ### 时滞微分方程理论 1. **稳定性分析**：特征方程方法、Lyapunov函数方法 2. **分岔理论**：Hopf分岔、稳定性切换 3. **周期解**：延拓方法、不动点理论 ## 结论与讨论 ### 主要结论 1. **平衡点分类完整**：给出了所有可能平衡点的存在条件和分类 2. **稳定性条件明确**：提供了显式的稳定性判据和临界时滞值 3. **周期解存在性**：证明了在化疗周期作用下周期解的存在 4. **分岔现象清晰**：揭示了系统从稳定到振荡的转换机制 ### 局限性 1. **参数简化**：$a=0$的假设虽有生物学依据但限制了一般性 2. **数值验证不足**：理论结果缺乏充分的数值模拟验证 3. **生物学验证**：模型预测需要更多实验数据支持 ### 未来方向 1. **三维扩展**：考虑更多细胞类型（如调节性T细胞） 2. **空间异质性**：引入扩散项研究空间模式 3. **随机效应**：考虑噪声对系统动力学的影响 4. **优化控制**：基于模型设计最优化疗方案 ## 深度评价 ### 优点 1. **理论严谨**：数学分析完整，证明过程清晰 2. **模型合理**：生物学背景明确，参数含义清楚 3. **方法创新**：双时滞机制和广义增长函数的结合 4. **结果实用**：为化疗方案设计提供理论指导 ### 不足 1. **数值部分薄弱**：第6节"数值验证"标注为"构建中" 2. **参数估计缺失**：未提供实际参数估计方法 3. **模型验证不足**：缺乏与临床数据的对比 ### 影响力 1. **理论贡献**：丰富了时滞微分方程在生物数学中的应用 2. **应用价值**：为肿瘤治疗的数学建模提供新思路 3. **方法论意义**：延拓方法在周期解分析中的应用具有一般性 ### 适用场景 1. **理论研究**：时滞微分方程、动力系统理论 2. **生物数学**：肿瘤增长模型、免疫系统建模 3. **医学应用**：化疗方案优化、治疗效果预测 ## 参考文献 论文引用了27篇相关文献，主要包括： - 肿瘤增长模型的经典文献 - 时滞微分方程稳定性理论 - 肿瘤-免疫相互作用的数学模型 - 分岔理论和周期解的存在性定理 --- **总体评价**：这是一篇理论性很强的数学生物学论文，在时滞微分方程和肿瘤-免疫系统建模方面做出了重要贡献。理论分析详尽，方法创新，但数值验证部分有待完善。对理解肿瘤-免疫系统的复杂动力学具有重要意义。