基于病毒载量与ISG表达结构化的新型病毒感染模型揭示干扰素信号调控感染进程与严重性的机制

本研究开发了一种整合细胞内病毒-干扰素刺激基因(ISG)动态的病毒感染数学模型，通过将感染细胞按病毒载量(v)和ISG表达(i)进行结构化，实现了从分子互作到组织水平的跨尺度模拟。研究发现增强的ISG诱导会通过抑制I型干扰素(IFN-I)产生和延缓健康细胞耗竭而延长感染进程，同时证实ISG介导的病毒抑制效应是控制感染严重程度的关键。该模型为解析先天免疫中病毒-宿主相互作用提供了计算高效的灵活框架，可用于评估靶向ISG通路的治疗策略。

病毒感染一直是全球公共卫生面临的重大挑战，从快速传播的流感病毒和SARS-CoV-2到持久潜伏的HIV和丙型肝炎病毒，不同类型的病毒给人类健康带来多样化威胁。病毒感染究竟会走向急性自限还是慢性持续，很大程度上取决于病毒与宿主免疫系统的复杂博弈。在先天免疫应答中，I型干扰素（Type I interferons, IFN-I）扮演着核心角色，它能够激活数百个干扰素刺激基因（Interferon-stimulated genes, ISGs），从而在细胞内建立抗病毒状态。然而，ISG与病毒在细胞内的分子互作如何影响感染的整体进程和严重程度，至今仍有许多未解之谜。

传统基于常微分方程（ODEs）的病毒动力学模型虽能模拟组织水平的感染动态，却忽略了细胞内部的异质性。而多尺度模型和基于个体的模型虽然能刻画细胞分子过程，但计算成本高昂且难以进行数学分析。为此，研究人员开发了一个新型数学模型，首次将感染细胞按照细胞内病毒RNA浓度和ISG表达水平进行结构化，从而巧妙地将细胞内病毒-ISG动态整合到组织水平的病毒感染框架中。

该研究基于先前SARS-CoV-2先天免疫模型构建，保留了病毒动力学核心参数，但将细胞类型简化为健康细胞，使模型更适用于体外细胞培养系统并易于推广到其他疾病。模型通过偏微分方程（PDE）描述感染细胞在病毒-ISG参数空间中的分布变化，跟踪四个关键变量：健康细胞浓度C(t)、结构化感染细胞浓度C_n(v,i,t)、细胞外自由病毒浓度V(t)和细胞外IFN-I浓度I(t)。

在技术方法层面，研究者采用有限差分法数值求解PDE方程，时间步长dt=0.001天，空间步长dx=0.05，使用迎风格式处理输运项，梯形法近似积分。ODE部分采用显式欧拉法求解。参数估值主要来自文献报道的病毒复制、降解和分泌动力学数据，包括病毒复制率α=2.77 day^-1、ISG由IFN-I的激活率β=5×10^-7 pg/(mL·day)、病毒分泌率ρ=0.53 day^-1等关键参数。模型使用Python实现，计算30天感染过程仅需52秒，保证了高效的系统性参数分析。

3.1. Validation using COVID-19 and HIV patient data

研究人员使用pre-alphaCOVID-19毒株和HIV患者的病毒载量数据对模型进行了验证。通过将初始病毒接种量调整为V₀=1.2 copies/mL并提高健康细胞浓度至1,530，模型准确再现了pre-alpha患者病毒载量的时间动态和峰值特征。对于HIV感染，通过调整ISG激活率β=1.4×10^-7 pg/(mL·day)、病毒产生率k₃=730 copies/(mL·day)等参数，模型成功拟合了慢性感染患者的病毒持久化模式，证明了模型在急性和慢性感染中的适用性。

3.2. Elevated ISG induction by type I IFN prolongs infection by suppressing type I IFN production in infected cells

研究发现增强的ISG激活虽然降低了病毒载量峰值，却意外地延长了感染持续时间。当ISG激活率较高时，感染30天后仍可检测到显著病毒载量，而低ISG激活则导致感染在20天内清除。机制分析表明，高ISG表达抑制了感染细胞的IFN-I产生，同时通过限制病毒分泌减缓了健康细胞的耗竭，从而为病毒持久存在创造了条件。感染细胞表型分析揭示，高ISG激活情况下，高病毒/高ISG表型细胞比例达到43%，这些细胞病毒分泌能力低但可持续存在。

3.3. Effective suppression of viral replication by ISGs is essential for controlling infection severity

通过系统调整ISG抑制病毒复制的效能（efficacy，eff），研究发现当eff<67%时，感染面积 under the curve (AUC)与ISG效能呈线性负相关；而当eff>67%时，这种关系转为指数型，小幅提升效能即可大幅降低病毒负担。峰值病毒载量与ISG效能的关系类似：eff<80%时呈线性负相关，超过该阈值后转为指数关系。这些结果表明，高效能的ISG介导的病毒复制抑制是控制感染严重程度的关键，且存在明显的效能阈值效应。

3.4. Moderate viral secretion by infected cells maximizes viral load production

病毒分泌率ρ对感染进程呈现非单调影响。较低分泌率（ρ=0.2 day^-1)导致病毒在细胞内滞留，促进细胞凋亡，从而限制感染传播；中等分泌率（ρ=0.8 day^-1)使病毒产量最大化；而过高分泌率（ρ=2.0 day^-1)反而减少细胞内病毒复制，降低总体病毒产量。这种钟形曲线关系表明病毒进化可能面临分泌策略的权衡，适度分泌最有利于病毒传播。

研究讨论部分指出，该模型通过将感染细胞按病毒载量和ISG表达结构化，首次在计算高效的框架下实现了细胞内病毒-ISG互作与组织水平感染动态的整合。模型揭示了IFN-I信号的双重角色：适度激活可控制感染严重程度，而过强激活反而通过调节病毒分泌和防止健康细胞耗竭促进病毒持久性。这为理解慢性病毒感染中持续的ISG表达提供了理论解释。

研究发现ISG抑制病毒复制的效能存在临界阈值（67-80%），只有超过该阈值才能显著降低病毒负担。这一发现对针对ISG通路的治疗策略具有重要启示：可能需要实现高度有效的ISG激活才能获得临床收益。此外，病毒分泌率的非单调效应提示，调控病毒释放可能成为新的治疗思路，这与当前主要针对病毒复制的策略形成互补。

该研究的优势在于模型兼具计算效率与生物学真实性，参数大多来自实验数据，且适用于急性和慢性感染建模。局限性包括未考虑适应性免疫应答、病毒株进化、空间效应以及IFN-I对健康细胞的保护作用。未来工作将扩展模型纳入T细胞、B细胞应答，以及抗病毒治疗与外源IFN-I给药的组合效应，为针对ISG-病毒互作的精准治疗策略提供理论指导。

论文发表在《Journal of Thermal Biology》，通过创新的数学建模方法揭示了细胞内病毒-ISG动态平衡如何决定感染结局，为理解病毒感染免疫调控提供了新视角，也为开发靶向先天免疫通路的抗病毒策略奠定了理论基础。