[论文解读] Analysis of a mathematical model for interactions between T cells and macrophages

Alan D. Rendall|arXiv (Cornell University)|Jun 13, 2010

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 13被引用 41

一句话总结

本文严格分析了由 Lev Bar-Or 提出的基于细胞因子信号传导的 T 细胞与巨噬细胞相互作用的四维常微分方程模型，证明在参数受限条件下系统收敛至平衡态，且巨噬细胞抗原呈递不改变动力学的定性行为。关键发现是，所有观察到的行为——双稳态、收敛性以及 Th1/Th2 反应的主导性——均可通过仅包含 T 细胞的简化两变量模型复现，表明在此框架下巨噬细胞并非核心免疫系统动力学所必需。

ABSTRACT

The aim of this paper is to carry out a mathematical analysis of a system of ordinary differential equations introduced by R. Lev Bar-Or to model the interactions between T cells and macrophages. Under certain restrictions on the parameters of the model, theorems are proved about the number of stationary solutions and their stability. In some cases the existence of periodic solutions or heteroclinic cycles is ruled out. Evidence is presented that the same biological phenomena could be equally well described by a simpler model.

研究动机与目标

• 严格分析 Lev Bar-Or 提出的基于细胞因子的 T 细胞与巨噬细胞相互作用数学模型的长期动力学行为。
• 确定引入巨噬细胞抗原呈递是否会在更简单模型中未出现的情况下，引入定性上新的动力学行为。
• 在何种条件下系统表现出唯一或多个平衡解，并评估其稳定性。
• 研究在完整系统中是否可能出现周期解、异宿环或混沌行为。

提出的方法

• 该模型被表述为一组四阶常微分方程，描述由 T 细胞（C1^T, C2^T）和巨噬细胞（C1^M, C2^M）产生的细胞因子浓度随时间的演化。
• 通过变量替换定义总细胞因子浓度 z1 = C1^T + C1^M 和 z2 = C2^T + C2^M，将系统简化为在不变流形上的二维动力学。
• 使用 Sigmoid 激活函数 g(x) = ½(1 + tanh(x − θ)) 来建模细胞因子对细胞活性的非线性反馈，其中 g’(x) = 2g(x)(1 − g(x))。
• 在参数受限条件下，利用竞争系统理论和雅可比矩阵分析推导出稳定性与收敛性定理。
• 将分析扩展至排除巨噬细胞影响的截断模型，方法是当 i ∈ {1,2} 且 j ∈ {3,4} 时将相互作用系数 a_ij 设为零。
• 将原始模型（Lev Bar-Or, 2006）的数值证据与分析结果进行比较，以评估一致性及模型简化的潜力。

实验结果

研究问题

• RQ1 在何种参数条件下，系统具有唯一的全局吸引平衡解？
• RQ2 引入巨噬细胞抗原呈递是否会导致出现周期解或异宿环等定性上新的动力学行为？
• RQ3 是否必须使用完整的四维模型才能复现原始模型中观察到的 Th1/Th2 双稳态主导现象？还是仅含 T 细胞的简化模型即可复现相同动力学？
• RQ4 是否存在多个稳定平衡解共存的参数区域？其存在性如何依赖于反馈强度与细胞因子阈值？
• RQ5 在一般参数设置下，系统是否表现出复杂行为，如奇异吸引子或同宿轨道？

主要发现

• 当参数满足 C ≤ min{A, B} 时，所有解在 t → ∞ 时收敛至平衡解，且系统为竞争系统，确保无振荡地收敛。
• 当 d1 = d2 = d 且 θ = 0 时，系统退化为标量方程，其动力学可复现所有观察到的行为，包括双稳态。
• 在定理 3 的条件下，平衡解的数量最多为三个，其中最多两个为稳定解，其数量取决于反馈参数的相对强度。
• 尽管原始模型的数值图示暗示可能存在周期解、异宿环或奇异吸引子，但在所分析的参数区域中未发现此类行为的证据。
• 包含巨噬细胞的完整模型的动力学与排除巨噬细胞影响的简化两变量模型在定性上完全一致，表明巨噬细胞抗原呈递并未引入本质性的新动力学。
• 系统表现出双稳态：根据初始条件的不同，解会收敛至 Th1 主导或 Th2 主导的平衡态，与免疫应答极化现象一致。