整合CTL免疫应答、潜伏与感染年龄及饱和感染率的HIV年龄结构模型稳定性分析与应用

《Mathematical Biosciences》：Stability analysis of a novel HIV model with CTL immune response, latency age, infection age and a saturated infection rate

本文针对HIV治疗中病毒潜伏库和免疫动力学难题，构建了一个创新性的年龄结构模型，该模型整合了CTL（细胞毒性T淋巴细胞）免疫反应、潜伏年龄、感染年龄和饱和感染率。通过严谨的理论分析，证明了模型解的存在性、有界性，并探究了其渐近平滑性与一致持续性。研究明确了系统存在无病平衡点E0、CTL未激活感染平衡点E1和CTL激活感染平衡点E2三类平衡点，其共存依赖于无免疫激活复制数R1和免疫激活复制数R2两个关键阈值参数。利用特征方程和Lyapunov函数，完成了各平衡点的局部与全局稳定性分析。进一步，模型成功拟合了真实患者病毒载量数据，验证了其生物学合理性，并通过与常微分方程模型对比，凸显了引入年龄结构的重要性。数值模拟抗逆转录病毒治疗，也再次证明了CTL免疫反应在HIV感染过程中的关键作用。

　　

Section snippets

模型推导 (Model derivation)

在t时刻，设C(t)为CTL的数量。设e(t, a)和i(t, b)分别表示潜伏年龄为a的潜伏感染T细胞浓度，以及感染年龄为b的产毒性感染T细胞浓度。未感染T细胞以速率h不断生成，以非线性饱和速率βT(t)V(t)/(1+αV(t))被自由病毒感染，并以自然死亡率μ₁衰减。新感染的T细胞中，一部分比例f进入潜伏感染状态e(t, a)，其余部分则进入产毒性感染状态。

平衡点与复制数 (Equilibria and reproduction numbers)

显然，E₀= (T₀, 0, 0, 0, 0) 是模型的无感染平衡点，其中 T₀= h/μ₁。接下来，我们探讨CTL未激活的感染平衡点E₁= (T₁, e₁(a), i₁(b), V₁, 0) 所满足的方程。从方程组的第一个方程，我们可以推导出 T₁= h(1+αV₁) / [μ₁+ (β+αμ₁)V₁]。从第二个方程，我们得到 e₁(a) 的表达式，它揭示了潜伏细胞池随年龄a的指数衰减规律。

平衡点的局部稳定性 (Local stability of equilibria)

为了探究模型在平衡点附近的动力学行为，我们接下来在平衡点处对系统进行线性化。通过在平衡状态E*附近对原系统进行泰勒展开并忽略高阶无穷小项，我们得到了相应的线性化系统。进一步，通过分析线性化系统的特征方程及其特征根的符号性质，我们可以确定各个平衡点的局部稳定性，从而理解系统在受到微小扰动时的行为是趋向恢复（稳定）还是发散（不稳定）。

数据来源 (Data source)

抗逆转录病毒治疗(ART) 能够迅速抑制病毒复制并显著降低病毒载量，但这同时也可能改变CTL的功能状态，包括其增殖能力、细胞毒性和细胞因子分泌，从而干扰对CTL反应的准确评估。为了确保CTL免疫反应的真实性，并消除ART对病毒动力学分析的混杂影响，我们特意选择了未曾接受过ART治疗的HIV感染者的数据进行分析。

结论 (Conclusion)

在本文中，我们提出了一个新颖的、具有年龄结构的HIV感染动力学模型。这个模型整合了潜伏年龄和感染年龄两个维度，并引入了非线性饱和感染项以及CTL介导的免疫反应机制。该模型旨在捕捉HIV感染过程中的阶段特异性动力学行为：包括潜伏病毒库的形成与激活、被感染细胞的生命周期及其被CTL清除的过程、病毒载量的波动与稳定，以及CTL细胞的增殖与衰减等复杂相互作用。