跳转至

第十三章：血液读书笔记

作者

本章作者为 James P. Keener 与 James Sneyd。本章把 Vol I Ch 1（生化反应、酶动力学）和 Ch 7（钙动力学）建立的"分子结合 + 时间延迟"建模工具，扩展到血液系统的系统生理学——血浆化学、红细胞/白细胞/血小板的动力学调控、血红蛋白氧结合曲线、CO2/HCO3⁻/H+ 缓冲系统、白细胞趋化性与炎症反应。本章的数学内容集中在造血周期性疾病（cyclical neutropenia、CML）的延迟微分方程 + Hopf 分岔分析，以及白细胞趋化的 Keller-Segel 方程这两个标志性系统生理建模问题。

内容概述

本章分为 4 大主题：(1) 13.1-13.2 血液基础与造血——血浆化学（理想气体定律、Henry 定律、血红蛋白氧解离曲线）、血细胞分化谱系、造血干细胞分化与红细胞生成素（EPO）反馈、周期性血液疾病（周期中性粒细胞减少症 CN、慢性粒细胞白血病 CML）的延迟微分方程模型；(2) 13.3 红细胞与气体运输——肌红蛋白/血红蛋白的协同结合（MWC、KNF、Adair 中间模型）、Hill 方程、Bohr 效应、Haldane 效应、详细平衡原理；(3) 13.4 白细胞——白细胞趋化性的偏微分方程模型（Keller-Segel）、Zigmond 趋化敏感性、炎症反应；(4) 13.5 凝血（本章中未完整展示）。本章最具数学深度的部分是 13.2.2-13.2.3 的造血延迟微分方程模型，它给出周期为 \(2d + X\) 的自发振荡（与实测 CN 19-21 天周期吻合），并通过外围控制 vs 局部控制的对比分析得出"CN 振荡不是由外围控制系统失稳引起"的关键结论。

核心方程与概念

13.1 血液化学

理想气体定律（13.3）\(P = CRT\)。分压（13.4）\(P_i = x_i P\)。Henry 定律（13.5）\(c_i = \sigma_i P_i\)。Table 13.2 溶解度（Molar/mmHg）：O2 \(1.4 \times 10^{-6}\)、CO2 \(3.3 \times 10^{-5}\)（CO2 溶解度是 O2 的 20 倍）、CO \(1.2 \times 10^{-6}\)。

13.2 造血与周期性疾病

细胞年龄结构模型（13.7-13.9）：细胞年龄密度 \(n(x, t)\) 满足守恒方程 \(\partial n/\partial t + \partial n/\partial x = -\beta n\)，出生边界 \(n(0, t) = F(N(t - d))\)（成熟延迟 \(d\) ≈ 5 天）。稳态解（13.10-13.12）\(F(N_0) = \beta N_0 / (1 - e^{-\beta X})\)，其中 \(X\) 是细胞最大寿命（红细胞 120 天）。

延迟微分方程（13.14，\(\beta = 0\) 极限）\(dN/dt + F(N(t-d-X))e^{-\beta X} - F(N(t-d)) = -\beta N\)。特征方程（13.17）\(F'(N_0)(e^{-\lambda d} - e^{-\lambda(d+X)}) = \lambda\)。Hopf 分岔条件（13.21-13.22）：\(\omega(2d + X) = 2n\pi\)，\(2dF'(N_0) \sin(\omega d) = -\omega d\)。振荡周期（13.24）\(T = 2d + X\)——这是本章最优雅的结论之一：周期 = 2 × 成熟延迟 + 寿命。

关键参数：正常红细胞 \(d = 5\) 天、\(X = 120\) 天，\(X/d = 24\)，远离 Hopf 临界线（图 13.6，\(X/d \approx 14\)）；若 \(X\) 减半（\(X = 60\) 天），\(X/d = 12 < 14\)，触发 Hopf，振荡周期 \(T = 2 \times 5 + 60 = 70\) 天。

13.2.3 周边控制 vs 局部控制

加权成熟时间模型（13.25-13.27）：\(dx/dt = -\alpha x + M(\tilde{x}(t - \tau_m))\)，其中 \(\tilde{x}\) 是 \(x\) 与 gamma 分布 \(g(\tau)\) 的卷积（成熟时间分布）。特征方程（13.31）\(\lambda + \alpha = \mu (a/(\lambda + a))^{k+1} e^{-\lambda \tau_m}\)。稳定性边界（13.38-13.39）\(\alpha(\omega) = -\omega / \tan[\omega \tau_m + (k+1)\tan^{-1}(\omega/a)]\)，\(\mu(\omega) = -\omega R \sin[\omega \tau_m + (k+1)\tan^{-1}(\omega/a)]\)。

关键发现（Hearn et al. 1998）：用 CN 患者参数（\(\tau_m = 1.2\) 天 vs 正常 3.8 天）重算，振荡周期只能在 13.7-14.3 天之间，远低于实测 19-30 天——结论是CN 振荡不能由外围控制系统失稳解释，需要局部控制机制（干细胞动力学）。

13.3 红细胞与气体运输

肌红蛋白结合（13.40-13.41）：\(Y = [\text{O}_2] / (K + [\text{O}_2]) = P_{\text{O}_2} / (K_P + P_{\text{O}_2})\)，\(K_P = 2.6\) mmHg。血红蛋白 4 位点 Adair 模型（13.42-13.47）：\(Y = \sum_{j=0}^{4} j \alpha_j [\text{O}_2]^j / \sum_{j=0}^{4} \alpha_j [\text{O}_2]^j\)，其中 \(\alpha_j = \prod_{i=1}^{j} K_i^{-1}\)，\(K_1 = 45.9\)、\(K_2 = 23.9\)、\(K_3 = 243.1\)、\(K_4 = 1.52\sigma\) mmHg（Roughton et al. 1972）。

关键现象：\(K_4\) 远小于 \(K_1\)-\(K_3\)，血红蛋白倾向于"被填满"——这是正协同性（positive cooperativity）的化学基础。机制（MWC 模型）：血红蛋白存在 R（relaxed，高亲和）和 T（tense，低亲和）两种构象；Perutz 1970 提议 T 态由盐桥稳定，氧结合时盐桥断裂。

Bohr 效应：H+ 是氧结合的变构抑制剂（\(\bar{K}_1 > K_1\)）。详细平衡（13.55）\(K_1/\bar{K}_1 = K_2/\bar{K}_2\)——如果 H+ 抑制 O2 结合，则 O2 也抑制 H+ 结合，这是 Haldane 效应的热力学根源。功能意义：组织 [CO2] 高 → [H+] 升高 → Hb 释放 O2；肺中 [CO2] 低 → [H+] 降低 → Hb 结合 O2。碳酸酐酶催化 CO2 + H2O ⇌ H+ + HCO3⁻，加速 1000 倍（13.54）。

血红蛋白作为 H+ 缓冲（13.59）：[HCO3⁻] = h + n(Z + W) - T_0，每个 Hb 分子可结合 n = 10-20 个 H+，故组织新生的 H+ 几乎全部被 Hb 吸收而不改变血浆 pH（仅从 7.4 降至 7.35）。

13.4 白细胞趋化

双向运动模型（13.65-13.67）：\(J_c = v^+ n^+ - v^- n^-\)，\(\partial n^\pm/\partial t = \mp \partial(vn^\pm)/\partial x + p_\mp n^\mp - p_\pm n^\pm\)。稳态通量（13.69-13.70）\(J_c = -\mu \partial c/\partial x + V_c c - T_p v \partial v/\partial x \cdot c\)，其中随机迁移系数 \(\mu = T_p v^2\)、趋化速度 \(V_c = T_p v(p_- - p_+)\)。

Zigmond 趋化敏感性（13.72-13.73）：\(f = (1/2)[1 + \chi_0 \partial N_b/\partial x / (1 + \chi_0 \partial N_b/\partial x)]\)，其中 \(N_b\) 是结合受体数。小梯度极限 \(f \approx 1/2 + (\chi_0/2)(dN_b/da) \partial a/\partial x\)，大梯度极限 \(f \to 1\)。Keller-Segel 极限（13.77）\(J_c = -\mu \partial c/\partial x + \chi c \partial a/\partial x\)，其中 \(\chi = v \chi_0 N_b'(a)\)。

炎症反应三变量模型（13.78-13.79）\(\partial b/\partial t = \mu_b \partial^2 b/\partial y^2 + (k_g - k_d c)b\)（细菌），\(\partial a/\partial t = D \partial^2 a/\partial y^2 + k_p b\)（趋化因子），加上白细胞方程 \(\partial c/\partial t = -\partial J_c/\partial y\)。

关键结论

  1. 红细胞寿命 120 天、每天生成 \(3 \times 10^9\) /kg——骨髓每天替换约 1% 红细胞。
  2. Hopf 分岔周期公式 \(T = 2d + X\)：振荡周期 = 2 × 成熟延迟 + 细胞寿命。这是延迟微分方程理论的"教科书式"结果。
  3. CN 振荡不能用外围控制失稳解释（Hearn 1998）：周边控制模型只能预测 13.7-14.3 天周期，实测 19-30 天。结论：CN 振荡必须由干细胞动力学（局部控制）引起，不是反馈延迟。
  4. Hopf 临界线 \(X/d \approx 14\)：正常红细胞 \(X/d = 24\) 远离不稳定性，但任何缩短寿命的因素（如自身免疫性溶血性贫血、AIHA）可使 \(X/d\) 落入不稳定性区域。
  5. Bohr 效应 ⇔ Haldane 效应（详细平衡）：如果 H+ 抑制 O2 结合，则 O2 必抑制 H+ 结合，二者通过热力学耦合。
  6. 协同性化学基础：Hb 的 \(K_4 \ll K_1, K_2, K_3\)，血红蛋白"倾向被填满"——这是 sigmoidal 氧解离曲线的化学起源。
  7. 肌红蛋白储存、释放给线粒体：在 20 mmHg PO2 下，Hb 35% 饱和但 Mb 90% 饱和——Mb 把 O2 缓冲在肌肉中，与线粒体紧密配合。
  8. CO 毒性（13.3.2）：CO 与 O2 竞争同一结合位点但亲和力高 200 倍——0.5 mmHg CO（无 O2）即可 97% 饱和 Hb，有 O2 时仅 0.7 mmHg CO 即可造成组织缺氧。
  9. 胎儿 Hb 的左移：胎儿 Hb 在低 PO2 下比成人 Hb 多带 30% O2——这是胎儿从母体获取 O2 的关键。
  10. Keller-Segel 趋化方程 \(J_c = -\mu \partial c/\partial x + \chi c \partial a/\partial x\) 是细胞趋化性的标准 PDE 模型，但高维 Keller-Segel 会出现 blow-up（细胞浓度爆聚）——这一现象与肿瘤侵袭、白细胞浸润的病理机制相关。

挑战和开放性问题

  1. CN 振荡的真实机制——周边控制模型预测周期 13.7-14.3 天 vs 实测 19-30 天，4-15 天的差距至今未由延迟微分方程模型弥合。可能涉及：(a) 干细胞周期动力学；(b) 多细胞类型的耦合反馈；(c) 骨髓微环境的物理结构（造血龛 niche）。
  2. 白细胞趋化的分子机制——13.4 节承认"Exactly how this is accomplished, how this should be modeled, and how well the model represents this behavior is the topic of this section"，但 Keller-Segel 是宏观唯象模型，G 蛋白偶联受体下游的信号转导（PI3K、AKT、Rho GTPases）完全缺失。
  3. 炎症反应模型的 blow-up——Keller-Segel 在高维 \(n \geq 2\) 时会出现细胞密度爆聚（finite-time blow-up），对应白细胞过度浸润（如 ARDS、败血症）。本章未讨论爆聚条件、临界 \(\chi/\mu\) 比值的精确阈值。
  4. 造血干细胞分化的随机性——13.2 节明确指出"it has proven difficult to separate the precursors of one cell type from those of another"。单细胞测序（scRNA-seq）2014 年后才给出更清晰的谱系图，但数学模型仍以平均场为主，缺乏随机谱系追踪的 ODE/PDE 公式。
  5. 血凝级联的 ODE 网络——本章未涉及血凝（止血/血栓），但血凝级联是 13.5 节本应包含的 50+ 个 ODE 网络。凝血酶生成的"阈值/开关"行为（少量凝血酶触发大量生成）是生化正反馈的经典例子，需要独立的数学处理。
  6. 红细胞生成的长期反馈——13.2 节说"The control mechanisms that operate when the red blood cell count is too high (a condition called polycythemia or erythrocytosis) are less clear"——多血症（红细胞过多）的反馈机制至今未明，可能涉及低氧诱导因子（HIF）的降解途径。
  7. Bohr 效应的精确模型——(13.55) 是热力学约束，但 \(K_1\)、\(K_2\)、\(\bar{K}_1\)、\(\bar{K}_2\) 的真实值随温度、2,3-BPG 浓度、CO2 浓度的变化是临床血气分析的精确校正依据，但本章只给定了性描述。

个人反思与批判性分析

本章的最大优点在于两个标志性应用：(1) 造血周期 = 成熟延迟 + 寿命——这是延迟微分方程 Hopf 分岔的"教科书级"应用，把临床病理（CN 19-21 天周期）直接与代数公式 \(T = 2d + X\) 对应；(2) Bohr-Haldane 互锁——通过详细平衡原理把两个看似独立的临床现象（Bohr 效应、Haldane 效应）联系起来，是热力学在生物化学中优雅应用的典范。

本章的主要缺陷在于Hearn et al. 1998 的反结论被一笔带过——13.2.3 节末尾明确说 peripheral control 模型不能解释 CN 振荡，但本章没有给出替代的干细胞动力学模型。这是教科书常见的问题：当一个假设被推翻时，缺乏替代理论。这给读者留下"CN 仍未被数学建模"的印象，而事实上有大量后续工作（Colijn & Mackey 2005a,b 等）用耦合干细胞 + 多细胞类型模型重现了 19-21 天周期。

13.3 节的 Bohr 效应 + 详细平衡推导 (13.48-13.55) 是本章最美的一段——它把临床医生习以为常的"CO2 高 → 释放 O2"现象与热力学约束关联起来。这一推导可以被推广到所有变构效应：例如酶的产物抑制、离子通道的电压依赖性失活、细胞信号通路的反馈抑制——核心都是"如果 A 调控 B，则 B 必调控 A"的热力学对称性。

对我自己研究血管生物力学的启发：(1) 细胞年龄结构 PDE 模型 (13.7) 是任何"细胞群体动力学"的标准起点——血管平滑肌细胞的增殖、巨噬细胞在动脉粥样硬化斑块中的聚集、内皮细胞的衰老，都可以用同一框架建模；(2) Bohr 效应的热力学对称性可以推广到血管张力调节——低 O2 引起血管舒张（hypoxic vasodilation），但血管舒张也降低局部 O2 消耗（降低壁剪切应力 → 减少内皮 ATP 释放），二者是否满足某种"详细平衡"原则是一个开放问题；(3) Keller-Segel 爆聚现象在血管生物学中可能与动脉粥样硬化斑块形成相关——单核细胞趋化进入血管壁的密度积累是否会触发斑块不稳定（"爆聚式"破裂）值得用类似数学框架分析。

值得指出，本章没有涵盖：(a) 凝血级联（13.5 节应是独立的一节但未在主流程中）；(b) 血液流变学（剪切率依赖黏度、红细胞聚集）；(c) 免疫学的 T/B 细胞激活网络。这些缺失严重削弱了本章作为"血液系统生理学"的完整性——希望下一版能补充。

重要参考文献

[X1] Mackey, M. C. (1978). A unified hypothesis for the origin of aplastic anemia and periodic haematopoiesis. Blood, 52, 5-15. — 13.2 节造血周期性疾病延迟微分方程模型的原型；周期中性粒细胞减少症（CN）的开创性数学建模。 [X2] Belair, J., Mackey, M. C., & Mahaffy, J. M. (1995). Age-structured and two-delay models for erythropoiesis. Mathematical Biosciences, 128, 317-346. — 13.2.2 节年龄结构 + 延迟模型；\(T = 2d + X\) 周期公式的来源。 [X3] Hearn, T., Haurie, C., Mackey, M. C., et al. (1998). Cyclical neutropenia and the peripheral control of white blood cell production. Journal of Theoretical Biology, 192, 167-181. — 13.2.3 节关键反结论：周边控制模型不能解释 CN 振荡。 [X4] Colijn, C., & Mackey, M. C. (2005a). A mathematical model of hematopoiesis—I. Periodic chronic myelogenous leukemia. Journal of Theoretical Biology, 237, 117-132. — 周期 CML 的耦合干细胞模型。 [X5] Colijn, C., & Mackey, M. C. (2005b). A mathematical model of hematopoiesis—II. Cyclical neutropenia. Journal of Theoretical Biology, 237, 133-146. — CN 的耦合干细胞模型；解释 19-21 天周期。 [X6] Perutz, M. F. (1970). Stereochemistry of cooperative effects in haemoglobin. Nature, 228, 726-739. — 13.3.1 节 MWC 模型的分子机制；T 态盐桥稳定机制。 [X7] Monod, J., Wyman, J., & Changeux, J.-P. (1965). On the nature of allosteric transitions: A plausible model. Journal of Molecular Biology, 12, 88-118. — 13.3.1 节 MWC 协同性模型的原型。 [X8] Adair, G. S. (1925). The hemoglobin system. VI. The oxygen dissociation curve of hemoglobin. Journal of Biological Chemistry, 63, 529-545. — 13.3.1 节血红蛋白 4 结合位点的发现；Adair 中间模型。 [X9] Bohr, C., Hasselbalch, K., & Krogh, A. (1904). Ueber einen in biologischer Beziehung wichtigen Einfluss, den die Kohlensäurespannung des Blutes auf dessen Sauerstoffbindung übt. Skandinavisches Archiv für Physiologie, 16, 402-412. — Bohr 效应的首次报告。 [X10] Roughton, F. J. W., DeLand, E. C., Kernohan, J. C., & Severinghaus, J. W. (1972). Some recent studies of the oxyhemoglobin dissociation curve of human blood under physiological conditions. Respiration Physiology, 16, 378-394. — 13.3.1 节 Hb 4 个 \(K_j\) 值的实验测定。 [X11] Tranquillo, R. T., & Lauffenburger, D. A. (1987). Stochastic model of leukocyte chemosensory movement. Journal of Mathematical Biology, 25, 229-262. — 13.4.1 节白细胞趋化双向运动模型。 [X12] Keller, E. F., & Segel, L. A. (1971). Model for chemotaxis. Journal of Theoretical Biology, 30, 225-234. — 13.4.1 节末尾 Keller-Segel 趋化方程的开创性论文。 [X13] Zigmond, S. H. (1977). Ability of polymorphonuclear leukocytes to orient in gradients of chemotactic factors. Journal of Cell Biology, 75, 606-616. — 13.4.1 节趋化敏感性 \(\chi_0\) 的测定。 [X14] Alt, W., & Lauffenburger, D. A. (1987). Transient behavior of a chemotaxis model. Journal of Mathematical Biology, 24, 691-708. — 13.4.2 节炎症反应三变量模型。 [X15] Haurie, C., Dale, D. C., & Mackey, M. C. (1998). Cyclical neutropenia and other periodic hematological diseases: Evidence for autoinstructive regulation. Blood, 92, 2629-2640. — 周期性血液疾病的自反馈调节假说。 [X16] Alberts, B., Bray, D., Lewis, J., Raff, M., Roberts, K., & Watson, J. D. (1994). Molecular Biology of the Cell (3rd ed.). Garland. — Table 13.1 血细胞类型与浓度数据来源。