季节性过敏性鼻炎渐进式常微分方程模型：从核心Th2轴到季节动态与方向性评估的21维多通路整合

《Frontiers in Immunology》：A progressive ODE model of seasonal allergic rhinitis: from core Th2 axis to 21-dimensional multi- pathway integration with seasonal dynamics and directional assessment

　　

季节性过敏性鼻炎（seasonal allergic rhinitis, SAR）影响全球约4亿人口，但目前缺乏一个能够将过敏原暴露至临床症状的完整免疫学级联反应进行定量整合的框架。既往各单独通路（Th2/IgE、TLR4/NF-κB、Th17/Treg）已被分别阐明，但它们的相互作用以及对症状生成的相对贡献仍知之甚少。研究人员构建了一个逐步扩展的常微分方程（ordinary differential equation, ODE）模型，该模型从10维核心Th2/IgE/肥大细胞/嗜酸性粒细胞轴逐步演进为整合TLR4/MyD88/NF-κB先天免疫与Th17/Treg失衡、并包含季节启动与耐受动态的21维系统。模型针对已发表的细胞因子数据进行校准（Leaker 2016, Linton 2025），并针对临床试验药物疗效进行验证（Corren 2021, Bauer 2023）。进行了Sobol全局敏感性分析（5,888个样本，21个参数）。生成了一个包含2,000例患者的虚拟队列，采用对数正态参数分布（变异系数CV = 30%），并通过PCA降维K-means聚类（基于六个输出指标的PC1–PC2）识别出虚拟亚型。通过将基础样本量加倍至N = 512进行收敛性验证：所有S1和ST指数的变化均小于0.013，证实原始N = 256时收敛性充分。核心模型再现了双相症状特征（即刻峰值7.0出现于1.9 h，晚期平台期3.4出现于12–24 h），并揭示了一种机制性二分现象：即刻峰值由肥大细胞-介质轴调控，而总症状负荷则由嗜酸性粒细胞-消退轴控制。19维扩展表明，TLR4与Th17/Treg通路主要调控晚期反应（通过IL-10抑制使AUC降低11.3%）。Sobol分析确定IL-10/Treg轴为总症状负荷的第3–4重要决定因素。季节扩展捕捉到免疫启动现象（IL-5升高46.1%）。三种虚拟亚型映射至已发表的AR（allergic rhinitis，过敏性鼻炎）表型。模型预测Th17/Treg比率与晚期严重程度的定向关联性强于IgE，与已发表临床数据一致。该渐进式ODE模型为理解SAR通路相互作用提供了一个定量框架，并产生了五个可供未来实验验证的可检验假设。虚拟队列中22%的模拟个体表现出严重的IL-17高表型，这一发现支持了Th17/Treg靶向治疗可能使一个生物学上明确界定的SAR患者亚组获益的假设，有望在真实患者队列中得到验证。

季节性过敏性鼻炎（seasonal allergic rhinitis, SAR）是全球性健康问题，影响全球约4亿人，由季节性气传花粉诱导的IgE介导的2型免疫反应所致。其病理生理学涉及过敏原穿透鼻上皮屏障后释放上皮源警报素（TSLP、IL-25、IL-33），激活先天性和适应性2型免疫。经典SAR反应呈双相特征：即刻相（数分钟内）由肥大细胞脱颗粒释放组胺驱动，晚期相（4–24小时）则以嗜酸性粒细胞浸润和持续的Th2细胞因子产生为特征。尽管各单独免疫通路已被广泛表征，但缺乏一个将过敏原暴露到临床症状表达的完整级联进行定量整合的数学框架，以用于理解通路相互作用、识别限速步骤及预测治疗干预结果。常微分方程（ODE）模型在传染病动力学和免疫反应建模中已成功应用，但其在过敏性鼻炎中的构建仍然有限。

研究人员采用渐进式建模策略，首先建立10维核心模型，捕获Th2/IgE/肥大细胞/嗜酸性粒细胞轴，随后扩展至19维以纳入TLR4/NF-κB和Th17/Treg失衡通路，最终扩展至21维并加入季节启动与耐受动态。该方法允许系统评估各通路对症状生成的相对贡献，并确保在每个建模阶段模型复杂性的合理性。模型被明确定位为假设生成工具，从中推导出五个具体的、可通过实验验证的预测，而非声称具备临床预测准确性。

研究人员构建了递进式ODE系统：10维核心模型包含过敏原、上皮警报素、树突状细胞、Th2细胞、Th2细胞因子、IgE、肥大细胞、炎症介质、嗜酸性粒细胞和症状评分共10个状态变量及26个参数；19维扩展新增TLR4/MyD88/NF-κB先天免疫级联和Th17/Treg失衡轴，引入31个新参数；21维季节模型加入启动（P）与耐受（T）两个状态变量以模拟整个花粉季节。参数化采用三阶段策略：基于已发表免疫学文献的参数、基于生物合理性的估计值以及对临床观察症状评分的校准。针对两项已发表的Nasal Allergen Challenge（NAC，鼻过敏原激发试验）研究的细胞因子时序数据进行折叠倍数校准（Leaker 2016与Baumann 2013），采用L-BFGS-B约束优化；针对Linton 2025的绝对浓度数据引入5个缩放因子进行绝对pg/mL校准。所有ODE系统采用LSODA方法求解。执行Sobol全局敏感性分析（Saltelli采样，N = 256，生成5,888个参数组合，涵盖21个关键参数）。构建包含2,000例患者的虚拟队列，56个模型参数和6个初始条件围绕基线值按对数正态分布（变异系数CV = 30%）采样。应用多层感知机（MLP，multi-layer perceptron）代理模型（隐藏层128-64-32，ReLU激活函数，Adam优化器）预测六项输出指标，结合置换重要性分析与主成分分析（PCA，principal component analysis）降维K-means聚类（基于PC1–PC2，K = 3，轮廓系数0.36）识别虚拟亚型。

3.1 核心模型动态（10维）
研究人员通过单次过敏原攻击模拟（A0 = 100 ng/mL），发现核心10维模型再现了SAR的特征性双相反应。在24小时内，过敏原在2小时内清除95%，肥大细胞在0.4小时迅速激活（84%脱颗粒），炎症介质在0.9小时达到峰值（半衰期约10分钟）并快速清除，而嗜酸性粒细胞在14小时达到峰值（174 cells/μL）。这揭示了明确的机制二分现象：肥大细胞-介质轴驱动即刻相反应，而嗜酸性粒细胞轴维持晚期相症状负荷。症状评分表现为即刻峰值（1.9小时为7.0）和随后的晚期平台期（12小时为3.4）。

3.2 扩展模型比较（10维 vs 19维）
在19维模型中，新增通路使总症状负荷降低（AUC 82.7 vs 93.3），晚期相负荷显著降低（52.3 vs 61.9），而即刻相几乎保持不变（9.5 vs 9.6）。消融分析显示，核心+Th17/Treg变体使AUC降低16.2%，而核心+TLR4变体反而使AUC增加3.2%，证实IL-10介导的免疫抑制是晚期症状负荷降低的主要来源，而非TLR4通路本身。

3.3 敏感性分析演变
单因素（OAT）分析显示，即刻峰值由肥大细胞-介质轴参数控制，而总症状负荷由嗜酸性粒细胞-消退轴参数主导。Sobol全局敏感性分析（5,888样本，21参数）对总症状负荷发现k_S2（ST = 0.35）和k_S5（ST = 0.31）为最重要参数，而k_{S_IL10}（ST = 0.13）和k_{IL10_1}（ST = 0.09）分别位列第三和第四，这一贡献是OAT分析完全无法捕捉的。

3.4 药物干预预测
在19维模型中，抗组胺药对即刻峰值的降低最为显著（47.1%），而dupilumab（抗IL-4Rα）提供了最大的总症状负荷降低（AUC降低23.2%）。Treg增强（低剂量IL-2）通过IL-10产生中度但持续的晚期抑制（AUC降低17.7%），抗IL-17提供适度获益（AUC降低3.8%）。参数扫描显示在最强阻断水平（×0.01）下，抗TLR4、抗IL-17及联合用药的AUC最大降幅分别为2.3%、4.3%和5.2%，Bliss独立性检验证实联合效应为次加性而非协同。

3.5 代理建模与虚拟队列探索
MLP代理模型对AUC（R² = 0.79）、Eos_peak（R² = 0.78）、Neu_peak（R² = 0.86）和IL17_peak（R² = 0.88）具有强预测力。置换重要性分析确定k_Med1（0.458）、k_S1（0.390）、k_Med2（0.383）和k_S5（0.267）为前四位最具影响力的特征。PCA降维K-means聚类将2,000例虚拟患者分为三个亚型：亚型1（40%）为高症状Th2-嗜酸性粒细胞反应型；亚型2（38%）为轻度低炎症型；亚型3（22%）为Th17高反应型，其IL-17和中性粒细胞水平显著升高。聚类稳定性通过10个随机种子（ARI = 0.98）和bootstrap重抽样（80%子样本，100次迭代：ARI = 0.94 ± 0.04）验证，K = 3具有最高轮廓系数（0.36）。

3.6 文献校准
折叠倍数校准针对NAC细胞因子数据调整了6个关键参数，使k_PI1降低70%，k_{IL10_1}增加188%，k_C1降低70%。绝对浓度校准针对Linton 2025数据引入5个缩放因子，模型产生的绝对细胞因子浓度落入已发表范围，IL-5和IL-13拟合良好（RMSE分别为4.2 pg/mL和2.7 pg/mL）。

3.7 临床试验比较：校准验证与独立验证
针对dupilumab的独立验证（Corren 2021），模型预测AUC降低27.8% vs 临床试验组间差异33.9%，峰值症状降低仅0.3% vs 临床24.4%。针对鼻内糖皮质激素（INCS，intranasal corticosteroids）的校准验证（Bauer 2023），模型预测7小时IL-13降低93.0%，与临床目标完全一致，TNSS AUC_0–30min降低27.2% vs 临床32.0%。INCS症状降幅被低估，因为糖皮质激素不直接抑制肥大细胞脱颗粒，而临床观察到的症状降低可能包含当前模型未捕捉的机制，如鼻血管通透性降低和黏膜水肿。

3.8 慢性季节暴露（21维模型）
56天季节模拟产生的症状轨迹遵循花粉包络线，症状峰值达7.58（第11天）。花粉-症状相关性为r = 0.67（Spearman），处于de Weger报道的范围（r = 0.51–0.65）内。在匹配花粉水平下，季节早期与晚期的症状差异为?0.37（模型）vs ?0.30至?0.41（de Weger 2011）。重复NAC启动模拟（5次挑战）显示IL-5从第1次到第5次挑战增加46.1%（189.9 → 277.4 pg/mL），与Orban 2020一致。

3.9 H4方向一致性评估：Th17/Treg作为晚期预测因子
三个虚拟亚型映射至已发表的AR表型。在2,000例虚拟队列中，基线Th17/Treg比率与晚期AUC相关（ρ = 0.079, p = 4×10^?4），而IgE无显著相关（ρ = 0.025, p = 0.26）。这种方向性排序（Th17/Treg > IgE）与Cai 2024的间接证据方向一致，但绝对效应量较弱（ROC AUC = 0.530 vs 0.519），表明H4得到部分支持：方向正确但效应极小（p值由大样本量N = 2,000驱动），需要在真实患者队列中进一步验证。

讨论部分指出，渐进式建模策略在不同复杂性阶段揭示了关键发现：10维核心模型确立了SAR的根本机制二分法，解释了为何H1抗组胺药能有效缓解急性症状，而靶向2型细胞因子的生物制剂对慢性症状提供持续获益。19维扩展表明TLR4与Th17/Treg通路主要调控晚期反应，而Th17/Treg比率在24小时内从0.20升至1.62（稳定于1.07），由IL-6依赖性Th17分化和相互Treg抑制驱动。21维季节扩展捕捉到免疫启动（IL-5 +46.1%）和早晚期症状差异（方向正确但幅度被低估）。Sobol全局敏感性分析揭示IL-10/Treg轴参数对总症状负荷的显著贡献是单因素扰动分析无法发现的，凸显了全局方法在多通路模型中的重要性。虚拟队列方法遵循定量系统药理学（QSP，quantitative systems pharmacology）范式，代表了过敏性鼻炎领域首次此类虚拟队列分析。22%的IL-17高表型亚群（亚型3）支持了抗IL-17治疗可能使标准Th2聚焦治疗应答不足的生物学定义亚组获益的假设。局限性包括：超过半数的模型参数（51%）仍为估计值；药物机制被简化；虚拟队列完全基于合成数据（对数正态分布，CV = 30%）；Th17/Treg相互拮抗以线性项建模而非RORγt/Foxp3转录因子动力学；未进行正式的参数可识别性分析；以及在极端IL-10抑制下症状评分可降至?0.48（影响<3%队列，已在后处理中clamp至(0, 10)区间）。

本文提出了一个逐步扩展的SAR ODE模型，该模型从10维核心Th2/IgE/肥大细胞/嗜酸性粒细胞轴演变为整合核心Th2轴、TLR4/NF-κB先天免疫和Th17/Treg失衡三条免疫学通路，以及启动与耐受两个动态过程的21维系统，涵盖急性24小时NAC与慢性8周季节两个时间尺度。模型经绝对细胞因子浓度校准和临床试验药物疗效验证，其关于Th17/Treg失衡作为晚期严重程度标志的预测与已发表临床数据方向一致。对2,000例患者虚拟队列的人工智能分析识别出三个映射至已发表AR表型的虚拟亚型。模型生成了五个具有具体验证方法的假设，定位为假设生成工具而非预测性临床模型。虚拟队列中22%模拟个体表现为严重IL-17高表型的发现，支持Th17/Treg靶向治疗可能使标准Th2聚焦治疗应答不足的SAR生物学定义亚组获益的假设，有望在真实患者队列中得到验证。