《Journal of Theoretical Biology》：From individual-based models to macroscopic dynamics of antimicrobial resistance
编辑推荐：
本论文引入并探讨了一个动力学框架，用以描述在易感–感染相互作用驱动的微生物感染背景下，被划分为易感、感染、康复和耐药四个室间的群体统计分布的时间演化。主要目标在于量化过度及不当使用抗菌药物所产生的影响——这种滥用在使部分感染个体转入耐药室间的过程中加速了耐药性
本论文引入并探讨了一个动力学框架，用以描述在易感–感染相互作用驱动的微生物感染背景下，被划分为易感、感染、康复和耐药四个室间的群体统计分布的时间演化。主要目标在于量化过度及不当使用抗菌药物所产生的影响——这种滥用在使部分感染个体转入耐药室间的过程中加速了耐药性的传播。该模型由一组玻尔兹曼（Boltzmann）型方程构成，捕捉易感个体与感染个体之间的二元相互作用，并辅以线性重分配算子以表征康复、耐药性形成及再感染过程。本研究表明，在掠碰（grazing-collision）极限下，该玻尔兹曼系统可由一组耦合的福克-普朗克（Fokker–Planck）方程很好地逼近。这一极限描述使得对动力学进行更易处理的数学分析成为可能，其中包括对群体密度长时间行为的刻画。分析结果突显了相互作用项如何驱动系统趋向稳定平衡，并量化了不当抗菌药物使用对耐药个体分布的影响。总体而言，上述结果提供了一个连接动理学理论与经典传染病建模的多尺度视角。
**论文主体解读**
抗菌药物耐药性（antimicrobial resistance, AMR）已成为全球公共卫生领域的持续性重大挑战。相关预测显示，到2050年AMR每年可能导致高达1000万人的死亡。过多且不恰当的抗菌药物使用，特别是抗生素滥用，是加速这一全球性耐药扩散的主要驱动因素。AMR不仅直接威胁人类健康，还严重削弱现代医学的既有成果，如重大手术、癌症化疗和器官移植均因耐药病原体的存在而风险显著增高，同时伴随巨大的社会经济负担。世卫组织（WHO）等多份报告估计，2019年AMR直接导致全球约127万人死亡，并间接关联近495万例死亡，而到2050年全球累计经济成本可能高达100万亿美元。因此，深入理解耐药性产生与传播的机制，并开发有效缓解策略，具有极为紧迫的现实意义。
数学建模在解析AMR的传播机制和评估干预策略方面发挥着关键作用。传统方法主要为基于Kermack–McKendrick经典框架的“仓室模型”，即将群体划分为易感（S）、感染（I）和康复（R）等仓室。然而，这类经典的常微分方程模型通常源自经验或启发式假设，并非从根本上由个体基本相互作用导出，因而往往忽略了个体交互的异质性。为弥补该缺陷，基于动理学理论的建模方法日益受到关注。这类方法通过统计分布描述个体在附加特征变量（如群体密度、接触率或行为特征）上的演化，从而天然地纳入群体异质性。然而，严格地将微观个体轨迹与宏观平均量联系起来的理论结果仍相对稀少。近期，Martalò等研究人员（2026）引入并研究了一种通用仓室流行病模型的均场表述，为理解群体异质性如何影响传染病最终稳态提供了定量框架。
本研究在Martalò等研究人员的框架基础上，结合Bonhoeffer等（1997）、Bergstrom等（2004）及Lehtinen等（2017）对AMR动力学的已有认识，提出了一种新颖的动理学模型。该模型将感染群体划分为两个亚群：感染敏感菌株的个体（I仓室）和感染耐药菌株的个体（J仓室）；同时将接受治疗的群体划分为获得了有效治疗从而康复的个体（T仓室），以及源于耐药感染或敏感感染但接受不当治疗而进入耐药状态的个体（R仓室）。这一结构中，不当治疗导致的由I仓室向R仓室的转移机制构成了本工作的核心特征，为量化不当抗菌药物使用对耐药性出现与传播的影响提供了基础。在该模型框架下，抗菌药物的不当使用对耐药性扩散的影响被证明远比针对性良好的治疗策略的益处更为显著。
为开展这一研究，研究人员建立了由双线性Boltzmann型方程构成的系统，其中包含线性重分配算子与二元相互作用算子。二元相互作用描述了易感与感染个体之间的感染过程，线性算子则刻画了康复、耐药性形成以及再感染等背景驱动的转移过程。研究者进而采用经典的掠碰极限（grazing-collision limit）方法，对模型进行简化，导出了带有时变扩散与漂移项的耦合Fokker–Planck（FP）方程系统。该均场限制性系统允许对群体密度长期行为进行解析研究，包括通过反Gamma型局部平衡态（quasi-equilibrium）的构建及能量距离（Energy distance）估计来刻画各仓室密度向稳定稳态的收敛。此外，通过计算矩系统，研究者获得了各仓室均值和方差的闭合演化方程，其中均值系统非线性的五个常微分方程与经典SIRS型模型结构一致。为分析非线性ODE系统的平衡态稳定性，研究者采用Jacobian矩阵并结合Routh–Hurwitz判据，推导出局部渐近稳定的充分必要条件。在数值验证环节，研究人员结构保持格式（structure preserving scheme）对所得Fokker–Planck方程系统进行了数值近似求积，计算域[0, L]取L=300，划分为N=8001个网格点，并以Δt=Δx/2的时间步长采用半隐式时间推进维持正性保持。模型参数选取满足平衡态存在与稳定的理论条件。
数值实验结果清晰地验证了理论推导。**Test 1. Convergence toward the equilibrium distributions** 中，研究人员数值模拟了系统(33)的演化，在t1=1, t2=10, t3=300三个时刻将各仓室数值解与理论准平衡分布fKq对比，确认长时间后各分布fK收敛至准平衡分布，从而验证了第四部分关于长时间行为的理论结论。**Test 2. Impact of the resistant strain** 中，研究人员系统性变化与不当治疗相关的参数γIR以及因有效治疗向敏感感染恢复/转移的参数γJI（取值范围{0.001, 0.01, 0.10}），考察耐药仓室R的行为变化。结果表明，增大γIR（不当治疗）导致耐药人群均值mR和方差vR均系统性增大，平衡态分布峰值降低、尾部变厚（对数–对数坐标下尤为显著），体现了耐药性扩散增强和群体异质性上升；而增大γJI（有效治疗）则产生互补效应，即均值与方差减小，分布更集中、尾部衰减更快，表明良好治疗能有效控制耐药性。值得强调的是，在相同参数范围内γIR变化对耐药仓室的影响明显强于γJI，这一不对称性提示不当治疗对驱动耐药性的作用强于针对性良好的治疗策略的缓解作用。
综合而言，本研究建立了一个新颖的模型体系，成功搭建了从微观个体相互作用到宏观群体动力学的多尺度桥梁。模型从基本原理出发推导出一个结合敏感与耐药菌株及不同治疗机制的宏观常微分方程模型，同时用动理学框架解析了相应的分布动力学，提供了宏观ODE系统所无法揭示的群体异质性信息。其结论指出，抗生素的滥用不仅直接促进耐药菌群增长，还可能触发涉及原易感个体的新感染途径。这一研究发现对未来系统分析无病平衡态、推导基本再生数及其与治疗参数依赖关系提供了明确的进一步研究方向，并强调了基于真实数据进行参数校准与验证的迫切性。
**研究结论翻译**：本文提出了一种在抗菌药物耐药性存在条件下微生物感染的新模型，旨在刻画感染的长期行为并量化不当抗菌药物使用的影响。研究人员引入了一个由易感–感染相互作用驱动的动力学框架，描述易感、感染、康复及耐药室间群体密度统计分布的时间演化。该模型由描述易感与感染个体二元相互作用的Boltzmann型方程组构成，并辅以涵盖康复与可能再感染动力学的线性重分配算子。模型与经典仓室流行病学模型通过均值演化建立联系，所得均值满足经典SIRS型模型，其相关参数由微观相互作用规则显式导出。采用掠碰极限方法后，Boltzmann系统可被耦合Fokker–Planck方程组良好逼近，从而显著简化了群体密度长时间行为的研究。在此框架内，研究人员证明了解以能量距离收敛于均场系统的唯一稳定平衡态，从而增强了对不当抗菌药物使用如何影响耐药个体分布的量化评估。这些结果为耐药性仓室模型与基于个体建模之间建立了多尺度桥梁，为均场仓室模型的推演以及多智能体系统的建模与控制开辟了新方向。在实际应用层面，动力学描述提供了无法从宏观ODE系统单独恢复的信息——后者仅刻画平均仓室大小，而前者则能表征完整分布及其异质性，包括耐药仓室的平衡分布fR∞。此外，所推导的宏观ODE系统为未来研究提供了若干方向，包括对无病平衡态的系统分析以及基本再生数及其对敏感与耐药感染之间转移率（尤其γIR和γJI）依赖关系的推导。最后，在第五部分第二节所考虑的参数范围内，数值结果表明γIR参数变化对耐药仓室的影响比γJI更显著；两者所量化的治疗相关转移机制之间的明确区分是本文建模的核心特征，上述对比为不当治疗可能比针对性良好的治疗策略对耐药性产生更强影响提供了在所述设置下的定性证据。
