均值回复环境噪声下具有潜伏细胞活化的HIV模型的随机阈值动力学

《Infectious Disease Modelling》：Stochastic threshold dynamics of an HIV model with latent cell activation under mean-reverting environmental noise

　　

受环境变异性对宿主内HIV动力学影响的启发，研究人员开发了一个包含潜伏细胞活化的随机HIV模型。感染率和活化率被建模为对数均值回复Ornstein-Uhlenbeck过程（Ornstein-Uhlenbeck process, OU），确保了参数的正性和生物上合理的波动。研究人员证明了全局正解的存在性和唯一性，并推导了表征持续和灭绝的随机阈值。如果持续阈值超过1，模型存在一个平稳分布（stationary distribution），表明病毒持续存在；相反，即使确定性模型预测持续，足够强的环境噪声也能几乎必然地驱动潜伏细胞、活性感染细胞和游离病毒灭绝。对相关Fokker-Planck方程（Fokker-Planck equation）的分析提供了准平衡附近协方差结构的显式表达式，量化了随机波动。数值模拟支持理论发现并说明了环境噪声如何重塑HIV动力学。研究结果强调了随机变异性在调节病毒持续和潜伏库动力学中的关键作用。

人类免疫缺陷病毒（Human immunodeficiency virus, HIV）是靶向人体免疫系统的逆转录病毒，自发现以来构成重大全球公共卫生挑战。经典数学模型虽能描述病毒复制、免疫反应和治疗效果，但忽略了潜伏感染细胞这一关键生物学特征：部分CD4⁺ T细胞感染后进入长期休眠状态，形成缓慢更新的病毒库，在治疗中断后可迅速释放病毒，且难以被免疫清除和药物根除。现有确定性模型无法解释环境波动（如免疫应答变异性、治疗依从性变化）对病毒动力学的影响，尤其不能预测环境噪声可能导致的病毒灭绝现象。为此，研究人员构建了一个包含潜伏细胞活化的随机宿主内HIV模型，将感染率k(t)和活化率a(t)建模为对数均值回复Ornstein-Uhlenbeck（Ornstein-Uhlenbeck, OU）过程，以同时保证参数正性、均值回复性和时间相关性。该研究旨在分析环境噪声如何影响病毒持续与灭绝，并量化随机波动对潜伏库动态的作用。研究人员通过严格数学分析定义了随机持续阈值R₀^s和随机灭绝阈值R₀^e，证明了当R₀^s>1时模型存在平稳分布（stationary distribution），当R₀^e<1时病毒几乎必然指数灭绝；进一步通过Fokker-Planck方程推导了准平衡点附近的多变量正态分布协方差矩阵。数值模拟验证了理论结果，揭示了噪声强度、回复速率、感染率和活化率对阈值及动力学行为的调节作用。该研究强调了随机变异性在调节病毒持续和潜伏库动力学中的关键作用，为HIV治疗策略提供了理论依据。论文发表在《Infectious Disease Modelling》。

研究人员开展研究用到的主要关键技术方法包括：（1）将感染率k(t)和活化率a(t)建模为对数均值回复Ornstein-Uhlenbeck（OU）过程，引入均值回复率λ_i和噪声强度σ_i，确保参数正性并捕获相关时间结构波动；（2）运用Lyapunov函数和随机分析，通过构造C²函数和紧集，证明平稳分布的存在性；（3）通过分析五维Fokker-Planck方程，利用Hurwitz矩阵和矩阵变换，推导准平衡点附近概率密度函数的显式表达式；（4）采用Milstein高阶数值方法进行离散化模拟，验证理论结果。该研究无样本队列来源，全部基于理论建模与数值模拟。

研究结果按论文主体部分小标题总结如下：

**全局解的存在唯一性（Section 3）**：研究通过构造停时和Lyapunov函数，证明了对任意正初始值，随机模型（1.5）存在唯一的全局正解，且解以概率1保持在正象限内，并进一步确定了不变集Γ（若σ₁²/2λ₁ < θ₁且σ₂²/2λ₂ < θ₂，则解始终在Γ中）。该结果保证了模型的生物学合理性。

**平稳分布（Section 4）**：通过构造非负C²函数并利用Young不等式和随机分析，证明了当随机持续阈值R₀^s = k₀a₀αNλ₁λ₂ / [δc(d_L+δ)((λ₁?σ₁²/2)(λ₂?σ₂²/2))] > 1时，模型在不变集Γ上存在至少一个平稳分布g*(·)，表明病毒在分布意义下持续存在。该条件比确定性基本再生数R₀更严格，反映了噪声对持续的抑制作用。

**概率密度函数（Section 5）**：当R₀ > 1时，研究对模型进行线性化变换，求解六维代数方程，获得了准平衡点P*附近状态变量（T, L, T*, V, k, a）的近似多元正态分布密度函数。根据不同的参数条件（如Δ₁、Δ₂、Δ₃、Δ₆是否为零），给出了协方差矩阵Σ的显式表达式（例如Δ₁=0, Δ₂≠0且Δ₆≠0时的形式）。该结果量化了随机波动的幅度和变量间的相关性，允许计算边际概率密度函数。

**灭绝条件（Section 6）**：通过定义R₀^e = (k₀a₀αN) / [δc(d_L+δ)] · exp(σ₁²/λ₁ + σ₂²/λ₂)，证明了当R₀ < 1且R₀^e < 1时，潜伏感染细胞L(t)、活性感染细胞T*(t)和游离病毒V(t)几乎必然指数衰减至零，即疾病灭绝。该条件表明，即使确定性模型预测病毒持续（R₀ > 1），足够强的噪声（使R₀^e < 1）也能驱动病毒灭绝，体现了随机性对动力学的基本改变。

**数值模拟（Section 7）**：采用Milstein高阶数值方法进行离散化模拟，参数值来自表1（如λ=10⁴ ml^?1day^?1, d_T=0.01, δ=1, c=23等）。示例7.1在R₀^s>1条件下展示了各变量的时间序列和频率直方图，验证了平稳分布存在；示例7.2通过4000万迭代点拟合边际密度曲线，确认了概率密度函数表达式与直方图吻合；示例7.3在R₀^e<1条件下展示了L、T*、V迅速衰减至零的灭绝行为；示例7.4分析了噪声强度σ₁和σ₂的影响，发现σ₁（感染率噪声）对各变量的破坏稳定效应大于σ₂（活化率噪声），且回复速率λ_i越大，从初始状态到持久状态的平均首次通过时间越长；示例7.5展示了阈值R₀、R₀^s和R₀^e随感染率k和活化率a的变化趋势，以及三维相图在持续与灭绝条件下的对比。

总结讨论部分：在本研究中，研究人员分析了一个包含潜伏细胞活化的随机宿主内HIV模型，其中感染率和活化率受到由对数均值回复Ornstein-Uhlenbeck过程建模的环境波动的影响。该公式确保了参数的正性，同时捕获了相关的、时间结构化的随机变异性。
- 如果随机持续阈值R₀^s>1，系统在不变集Γ上至少存在一个平稳分布，表明所有状态变量在分布意义下随机持续。
- 如果随机灭绝阈值R₀^e<1，则潜伏感染细胞L(t)、活性感染细胞T*(t)和游离病毒V(t)几乎必然指数衰减，表明即使确定性模型预测持续，足够强的随机扰动也能诱导灭绝。
- 当确定性基本再生数R₀>1时，准平衡点邻域内的解服从近似多元正态分布，其协方差由相关Fokker-Planck方程显式推导，量化了平衡附近的随机波动。

这些结果表明环境变异性从根本上改变了宿主内HIV动力学。作用于关键参数的噪声可以抑制病毒持续或诱导灭绝，强调了在具有潜伏库的模型中纳入随机效应的重要性。在数学上，该研究将确定性阈值理论扩展到具有潜伏感染、乘性均值回复噪声和多尺度动力学的随机框架。不同的随机阈值R₀^s和R₀^e突出了环境噪声如何修改有效繁殖机制。未来工作可以在额外生物学参数中引入随机性，探索非高斯或状态依赖噪声，并建立联系确定性和随机基本再生数的统一框架。