跳转至

第7章：钙动力学

作者

本章作者与 Ch1-Ch6 相同——James P. Keener 与 James Sneyd。本章是 I 卷"细胞生理学"中机制最丰富、应用最广的章节——钙是细胞内最普遍的第二信使，几乎涉及所有生理过程。本章的角色是把 Ch3 的单通道 Markov、Ch4 的扩散方程、Ch5 的兴奋性 + 时间尺度分离、Ch6 的行波传播整合到钙信号这一具体生物问题上。具体内容：（1）钙生理总览（§7.1 开头）——胞外 1 mM、胞内 0.1 μM、4 个数量级梯度；（2）IP3R 8 态模型（DeYoung-Keizer 1992）+ 6 态饱和结合模型（Sneyd-Dufour 2002）（§7.2.5）——8 态 vs 6 态的对比；（3）钙振荡模型——简单 2 变量模型（§7.2.6 末段 Fig. 7.11）+ 牛蛙交感神经元 CICR（§7.2.9 + Fig. 7.15）；（4）钙波 + 螺旋波（§7.3, Fig. 7.3）——Xenopus 卵母细胞的双参数分岔图（Fig. 7.16）；（5）钙缓冲（§7.4）——有效扩散系数 (7.78)、移动 vs 固定缓冲的对比；（6）离散钙源 + fire-diffuse-fire 模型（§7.5）——通道簇 + 跳跃式钙波；（7）随机钙火花（§7.6）——IPR Markov + Gillespie；（8）细胞间钙波（§7.7）——Sneyd 1994 机械刺激 + Höfer 2002 部分再生 + 肝细胞协调振荡；（9）均值场方程（§7.8 附录）——非均匀介质的多尺度均匀化。

内容概述

本章围绕"细胞如何控制 Ca²⁺ 浓度以实现信号转导"这一核心问题展开。

第一部分（§7.1-§7.2）钙生理与单细胞模型。§7.1 给出钙动力学的工具箱：IP3R、RyR、SERCA 泵、PMCA、NCX 交换器、线粒体、缓冲。胞内外梯度 10⁴ 倍（胞外 1 mM vs 胞内 0.1 μM）——这是为什么 Ca²⁺ 是快速信号载体的物理基础。§7.2 给出混合细胞模型（ODE 形式）。

IP3R 8 态模型（§7.2.5, DeYoung-Keizer 1992）：每个亚基有 3 个独立结合位点（IP3 + 激活 Ca²⁺ + 失活 Ca²⁺），3 个亚基都处于 \(S_{110}\) 态时通道开放。关键假设：\(k_5 > k_2, k_4\)（激活快于失活）——这是振荡和波的几何来源。QSSA 化简（7.16）–(7.20)：把 8 态 ODE 化简为 2 态（开放概率 \(\propto pc(1-y)h\) + 失活变量 \(h = 1-y\)）——这是 Ch1.6.3 准稳态在生物物理的经典应用。

IP3R 6 态饱和结合模型（§7.2.5 末段, Sneyd-Dufour 2002）：放弃"8 态 + 3 亚基"的精细结构，改为饱和结合速率（7.21）–(7.25)——这解决了 8 态模型的"动力学范围只跨 1 个数量级"问题（§7.2.5 中段）。

§7.2.6 简单振荡模型：(7.32)-(7.35) 给出 2 变量封闭细胞模型。关键现象（§7.2.6 Fig. 7.11）：\(c\) 零线是 N 形，\(y\) 零线单调 → Hopf 分岔 → 钙振荡。这与 FHN/HH 几何一致（Ch5 末段提到的"通用结构"）。

牛蛙交感神经元 CICR（§7.2.9, Friel 1995）：(7.46)-(7.52) 给出极简模型——Ca²⁺ 激活的释放函数 \(k_3(c)\) 是 Hill 函数。重要发现（§7.2.9 末段）：模型预测的钙振荡时程与实验精确吻合；后续实验验证了模型预测的离子流——这是"理论 + 实验"的双向验证。

§7.2.7 开/闭细胞：(7.43)-(7.45) 把 ODE 重写为 \(c_t\) (总 Ca²⁺) + \(c_e\) (ER Ca²⁺) 形式。重要（§7.2.7 中段）：\(c_t\) 慢变 + \(c, c_e\) 快变 → 两个时间尺度 → 振荡由快动力学调制。

第二部分（§7.3）钙波。Xenopus 螺旋波（Fig. 7.3A, Lechleiter-Clapham 1992）：卵母细胞直径 > 600 μm，足以容纳完整波 + 螺旋。双参数分岔图（Fig. 7.16）：Hopf + 同宿轨道结构（Ch6 通用结构）——\(s, p\) 平面上的 C 形同宿 + U 形 Hopf。行波速度 10-17 μm/s（实验范围）。

第三部分（§7.4）钙缓冲。线性缓冲（7.81）：\(D_{eff} = D_c K/(K+b_t)\) ——简单比例律。非线性缓冲（7.76）–(7.77)：有效扩散系数 (7.78) \(D_{eff} = (D_c + D_b \theta(c))/(1 + \theta(c))\) ——移动缓冲有"载体"效应（Connelly 1993, Wagner-Keizer 1994）。关键定理（§7.4.2）：固定缓冲对行波存在性无影响（(7.89) 在 \(D_b = 0\) 时退化为无缓冲条件）；移动缓冲可以阻止行波（\(K\) 小时）。

§7.4 末段 BAPTA vs EGTA：BAPTA（快高亲和）vs EGTA（慢高亲和）——两种实验缓冲对钙响应有"截然不同的影响"（Dargan-Parker 2003）。EGTA 让响应"巴尔干化"，BAPTA 让响应"全球化"。这是 §7.4 揭示的"缓冲不只是被动束缚"的反直觉现象。

第四部分（§7.5）离散钙源。§7.5.1 fire-diffuse-fire (FDF) 模型（Pearson-Ponce-Dawson 1998, Keizer 1998）：(7.116)-(7.121) 把通道簇 + 阈值触发 + 瞬时释放 + 扩散结合。关键结果（7.121）：无量纲速度 \(v = L/\tau = D_c / (\eta L)\) 是 \(\eta\) 的函数 → 速度与 \(D_c\) 线性（vs 连续反应-扩散系统的 \(\sqrt{D_c}\)）。为什么是线性：FDF 模型没有 uptake → 钙单调积累。

§7.5 修改 FDF（7.122）：加入线性 uptake \(k_s c\)。关键结果（7.124）：传播失败条件 \(c^* L/\sigma > g_{max}(\beta)\)，\(\beta = \sqrt{k_s L^2/D_c}\)。当 \(D_c\) 足够小 → 传播失败——修复了原 FDF 模型与连续方程的失配。

第五部分（§7.6）随机钙火花。§7.6.1 IPR 随机模型：用 Gillespie 算法（Ch2.9 引入）模拟 IPR Markov 状态图。关键发现（Swillens 1998, 1999）：每个簇约 25 个 IPR，间距 < 12 nm——这是单通道记录推算出来的。Shuai-Jung 2002-2003 用 Langevin 近似代替直接随机模拟——这是 §2.9.5 Langevin 方程在生物物理的应用。

§7.6.2 随机钙波：Falcke 2003a/b 的开创性工作。关键发现：在确定性模型非振荡的 [IP3] 区域，随机模型仍有"准周期"波。这是因为随机触发偶尔激发全簇发放——\(\Delta T_{av}/T_{av}\) 保持常数（Fig. 7.22）。这是"振荡的确定性解释可能错"的重要警示（§7.6.2 末段）。

第六部分（§7.7）细胞间钙波。§7.7.1 机械刺激模型（Sneyd 1994, 1995a, 1995b）：(7.125)-(7.131) 给出 \(p\)（IP3）+ \(c\)（Ca²⁺）+ \(h\)（失活）三变量模型。IP3 通过 gap junction 扩散让胞间波传播。关键预测：\(F\)（IP3 渗透率）需 > 1 μm/s 才能匹配实验。§7.7.2 部分再生（Höfer 2002）：(7.132) 引入 PLCδ 的 Ca²⁺ 正反馈 → 有限范围的波（vs 纯被动扩散的指数衰减 vs 完全再生的无限传播）。§7.7.3 肝细胞协调振荡（Dupont 2000）：相位波——每个细胞独立振荡，IP3 扩散协调相位。

第七部分（§7.8 附录）均值场方程。多尺度均匀化：(7.133)-(7.180) 推导非均匀介质的有效扩散张量。关键结果（7.180）：\(D_{eff} = (1/V_c) \int_{\Omega_c} (I + \nabla_\xi W) dV_\xi\)——有效扩散张量是微观结构的函数。Goel 2006 计算：规则 ER 结构 → 各向同性 \(D_{eff}\)；复杂 ER 结构 → 各向异性（垂直方向扩散 50% 大于其他方向）。这是 Ch12 心肌组织双域方程的预备。

主要数学结果： 1. IP3R 8 态模型 + QSSA 化简（7.16）-（7.20） 2. 钙振荡的 2 变量 ODE + Hopf 分岔（§7.2.6） 3. Friel 1995 CICR 模型（7.46）-（7.52）—— 钙振荡与 [Ca²⁺] 依赖性 4. 双参数分岔图 \(s, p\) 平面（Fig. 7.16）—— C 形同宿 + U 形 Hopf 5. 钙缓冲的有效扩散系数（7.78）—— 移动 vs 固定 6. 行波存在的面积定理（7.89）—— 固定缓冲不改变条件 7. Fire-diffuse-fire 模型（7.121）+ 线性 uptake 修正（7.124） 8. 随机钙火花的 [IP3] 依赖性 + 准周期振荡（§7.6.2） 9. 细胞间波部分再生机制（7.132） 10. 多尺度均匀化的有效扩散张量（7.180）

前置知识：Ch1（酶动力学 + QSSA）、Ch2（守恒 + 扩散 + Markow + Gillespie）、Ch3（单通道 Markov）、Ch5（兴奋性 + Hopf 分岔）、Ch6（行波 + eikonal-curvature）。

核心方程与概念

1. 钙工具箱与 8 态 IP3R 模型

胞内外梯度：\([Ca^{2+}]_{ext} \approx 1\) mM, \([Ca^{2+}]_{cyt} \approx 0.1\) μM——10⁴ 梯度。

IP3R 8 态模型（DeYoung-Keizer 1992, §7.2.5）：每亚基 3 位点（IP3 + 激活 + 失活）→ 3 亚基都开放 → 通道开放。速率常数（Table 7.1）：

\[ \alpha_m = 0.1(25-v)/(e^{(25-v)/10} - 1), \quad \beta_m = 4 e^{-v/18} \]

（这是 Ch5 HH 形式的"复刻"——但 Ch7 用于 IP3R 门控）。

QSSA 化简（7.16）–(7.20)：

\[ x_{110} = \frac{pc}{(p + K_1)(c + K_5)} h, \quad \tau_h \frac{dh}{dt} = h_\infty(c, p) - h \]

\(h\) 失活变量对应 Ch5 HH 的 \(h\)——两个独立的"双时间尺度"结构在细胞激活机制中重复出现。

2. 6 态饱和结合模型

Sneyd-Dufour 2002（§7.2.5 末段）：放弃 8 态结构，引入饱和结合速率 \(\phi(c) = k_2 c / (c + K_1)\)（Michaelis-Menten 形式）。

5 态 ODE（7.26）–(7.30)：\(R + O + A + S + I_1 + I_2 = 1\)——6 态但 5 个独立 ODE（守恒消去一个）。

开放概率（7.31）：\(P_O = (a_1 O + a_2 A)^4\)——4 亚基合作性。

优势（§7.2.5 中段）：动力学范围跨多个数量级（vs 8 态的 1 个数量级）——更符合实验。

3. 钙振荡模型

简单 2 变量封闭细胞（7.32）–(7.35）：

\[ \frac{dc}{dt} = (k_f P_O + J_{er})(c_e - c) - J_{serca}, \quad J_{serca} = \frac{V_p c^2}{K_p^2 + c^2} \]

\(c\) 零线：N 形（来自 \(P_O\) 包含 \(\min\) 函数）。\(y\) 零线：单调。单交点 + Hopf 分岔 → 振荡。几何与 FHN 完全一致（Ch5）——钙振荡 = 兴奋系统。

Friel 1995 牛蛙 CICR（7.46）–(7.52）：

\[ k_3 = \kappa_1 + \frac{\kappa_2 c^n}{K_d^n + c^n} \]

—— Hill 形式的 Ca²⁺ 激活释放。关键（§7.2.9 末段）：仅 Ca²⁺ 激活（无失活）即可产生生理学振荡——IPR 的失活是"可选项"。

4. 钙波 + 行波

Xenopus 螺旋波（Fig. 7.3A）：卵母细胞 > 600 μm → 完整波 + 螺旋。波长 150 μm，周期 8 s。

双参数分岔图（Fig. 7.16）：\(s\) (速度), \(p\) (IP3) 平面。C 形同宿 + U 形 Hopf——与 Ch6 的 FHN 行波结构完全相同。这是钙波 = 行波的进一步证据。

关键预测（§7.3.2 末段）：波速 10-17 μm/s（实验范围），随 \(p\) 增加而增加。

5. 钙缓冲

线性缓冲 + 准稳态（7.72）–(7.81）：

\[ b = \frac{b_t c}{K + c}, \quad \frac{\partial c}{\partial t} = \frac{K}{K + b_t} (D_c \nabla^2 c + f(c)) \]

——有效扩散系数 \(D_{eff} = D_c K/(K+b_t)\)——简单的比例缩放。

非线性缓冲（7.76）–(7.78)：\(D_{eff} = (D_c + D_b \theta(c))/(1 + \theta(c))\) ——移动缓冲的"载体"效应（Wagner-Keizer 1994）。

BAPTA vs EGTA 关键实验（§7.4 末段, Dargan-Parker 2003）：快高亲和（BAPTA）→ 全球化；慢高亲和（EGTA）→ 巴尔干化。这是 §7.4 揭示的"缓冲不只是被动束缚"的反直觉现象。

6. 离散钙源 + Fire-Diffuse-Fire

§7.5 FDF 模型（7.116）–(7.121）：

\[ \frac{\partial c}{\partial t} = D_c \frac{\partial^2 c}{\partial x^2} + \sigma \sum_n \delta(x - nL) \delta(t - t_n) \]

速度条件（7.121）：

\[ \frac{c^* L}{\sigma} = g(\eta), \quad \eta = \frac{D_c \tau}{L^2} \]

关键（§7.5.1 末段）：速度与 \(D_c\) 线性（vs 连续反应-扩散系统的 \(\sqrt{D_c}\)）——因为 FDF 没有 uptake。

§7.5 修改 FDF（7.124）：\(g_\beta(\eta)\) 加 \(\beta = \sqrt{k_s L^2/D_c}\) 项 → 当 \(D_c\) 足够小 → 传播失败——修复与连续方程的失配。

7. 随机钙火花

IPR 随机模型（§7.6.1）：用 Gillespie 算法（Ch2.9）模拟 IPR Markov。关键发现（Swillens 1998, 1999）：每个簇约 25 个 IPR，间距 < 12 nm。

Langevin 近似（Shuai-Jung 2002-2003）：用随机 ODE 代替直接模拟——Ch2.9.5 Langevin 方程在 IPR 系统的应用。

关键发现（§7.6.2, Falcke 2003a/b）：在确定性模型非振荡的 [IP3] 区域，随机模型仍有"准周期"波。\(\Delta T_{av}/T_{av}\) 保持常数（Fig. 7.22）。这是"振荡的确定性解释可能错"的重要警示。

8. 细胞间钙波

机械刺激（§7.7.1, Sneyd 1994）：\(p\) 通过 gap junction 扩散 → 触发邻居细胞 Ca²⁺ 释放 → 胞间波。关键预测：\(F\)（IP3 渗透率）需 > 1 μm/s。

部分再生（§7.7.2, Höfer 2002）：(7.132) 引入 PLCδ 的 Ca²⁺ 正反馈 → 有限范围的波（vs 完全再生的无限传播）。

肝细胞相位波（§7.7.3, Dupont 2000）：每个细胞独立振荡，IP3 扩散协调相位——相位波 vs 行波的对比。

9. 均值场方程（§7.8）

多尺度均匀化（7.162）–(7.180）：

\[ \frac{\partial u_0}{\partial t} = D_{eff} \frac{\partial^2 u_0}{\partial z^2} + f(u_0), \quad D_{eff} = \frac{1}{\bar R} \]

有效扩散张量（7.180）：

\[ D_{eff} = \frac{1}{V_c} \int_{\Omega_c} (I + \nabla_\xi W) dV_\xi \]

Goel 2006 计算：规则 ER → 各向同性；复杂 ER → 各向异性（垂直方向 50% 大）。

这是 Ch12 心肌组织双域方程的预备——非均匀介质的扩散必须用有效参数。

关键结论

  1. IP3R 8 态模型是钙信号传导的标准模型（§7.2.5）：Ca²⁺ 激活 + Ca²⁺ 失活 = 双时间尺度 = 振荡的通用结构（与 HH 的 \(m, h\) 同源）。QSSA 化简到 2 变量让所有后续分析可行。

  2. 钙振荡 = 兴奋系统（§7.2.6 + Ch5）：\(c\) 零线 N 形 + 单调恢复零线 + Hopf 分岔——与 HH / FHN 的几何完全相同。这是"通用兴奋性结构"在钙系统的具体实现。

  3. Friel 1995 模型展示"理论预测 → 实验验证"（§7.2.9 末段）：模型预测振荡时程 → 实验吻合；模型预测离子流 → 实验验证。这是 Ch1-7 一直强调的"建模 + 实验"双向循环的典范。

  4. 双参数分岔图是钙波的"指纹"（Fig. 7.16, §7.3.2）：C 形同宿 + U 形 Hopf——与 FHN 行波（Ch6）、HH 振荡（Ch5）完全相同的结构。这是"跨模型普遍性"的强证据。

  5. 缓冲不只是被动束缚（§7.4）：BAPTA（快高亲和）让响应"全球化"，EGTA（慢高亲和）让响应"巴尔干化"。缓冲的"动力学"决定它的功能——这是 §7.4 的核心洞见。

  6. Fire-diffuse-fire 模型揭示离散钙源的独特性（§7.5）：速度与 \(D_c\) 线性（vs 连续方程 \(\sqrt{D_c}\)）——因为 FDF 没有 uptake。当加入 uptake 后（(7.124)），\(D_c\) 小 → 传播失败——修复与连续方程的失配。

  7. 随机模型可以"创造"振荡（§7.6.2, Falcke 2003）：在确定性模型非振荡的 [IP3] 区域，随机模型仍有"准周期"波。这意味着"实验观测到振荡 ≠ 确定性模型有极限环"。对系统生物学有深远影响。

  8. 细胞间钙波有三种机制（§7.7）：被动扩散（Sneyd 1994，\(F > 1\) μm/s 必要条件）、部分再生（Höfer 2002，PLCδ 反馈）、相位波（Dupont 2000，IP3 协调）。这三种机制可能同时存在（§7.7.3 末段）。

  9. 多尺度均匀化是非均匀介质的关键工具（§7.8）：\(D_{eff} = (1/V_c) \int (I + \nabla_\xi W) dV\)——有效参数是微观结构的函数。这是 Ch12 心肌组织双域方程的预备。

  10. 钙的"工具箱"非常丰富（§7.1 末段, Berridge 2003）：每个细胞选择"工具箱"中的特定组件来精细调控时空动力学。建模也是"工具箱"式——每个模型组件对应一个生物组件。

挑战和开放性问题

  1. IP3R 8 态 vs 6 态的边界（§7.2.5）：两个模型在数学上不严格可比——8 态的"快激活 + 慢失活"假设 vs 6 态的"饱和结合"假设。何时用哪个？ 取决于实验数据的时间分辨率与动力学范围。

  2. 钙火花 vs 钙波的定量关系（§7.6.2 末段）：Falcke 2003 的随机模型虽然重现了 [IP3] 依赖性，但\(\Delta T_{av}/T_{av}\) 保持常数的机制——是 [IP3] 改变的"触发概率"还是"重置时间"？没有明确解析。

  3. 细胞间钙波的传播距离（§7.7.2 末段）：Höfer 2002 的部分再生机制让波传播更远，但不能无限传播。什么决定了"部分再生"的上限？这与组织异质性的关系是什么？

  4. BAPTA vs EGTA 的"巴尔干化 vs 全球化"机制（§7.4 末段）：Dargan-Parker 2003 提供了实验现象，但精确的模型解释还没有。这需要 EPR 局部 Ca²⁺ 微域的精细建模（§7.6.1 末段提到的"局部高浓度"）。

  5. 随机 vs 确定性的边界（§7.6.2 末段）：Falcke 2003 的"随机振荡"是否在某些细胞类型中实际上比确定性极限环更准确？真实细胞的 IPR 数（25/簇）是"小到需要随机 + 大到需要 Markov"的中间地带。

  6. 多尺度均匀化的有效性（§7.8）：\(D_{eff}\) 的显式计算（Goel 2006）依赖周期性微观结构。真实细胞的 ER 是非周期、非各向同性的——是否需要重新定义 \(D_{eff}\)？

  7. 细胞类型差异（§7.2.6 末段）：不同细胞类型（肝细胞、胰腺腺泡、气道平滑肌）的钙振荡模型形式相同但参数不同。是否存在通用的"参数标度律"？这能让一个模型的预测推广到其他细胞。

  8. 缓冲的多缓冲动力学（§7.4 末段 + Exercise 12）：真实细胞有 200+ Ca²⁺ 缓冲蛋白（Carafoli 2001）。多缓冲的复杂动力学如何影响振荡和波？是否需要 200 个 ODE？

  9. IP3 动力学 vs Ca²⁺ 动力学的相对重要性（§7.2.8 末段, Sneyd 2006）：不同细胞类型的实验给出不同结论。胰腺腺泡 = IP3 振荡必要，气道平滑肌 = 恒定 IP3 即可。这些差异的机制是什么？

  10. 2D vs 3D 钙波（§7.3.1 末段）：Xenopus 卵母细胞是 2D 螺旋波的良好平台。真实细胞（如肝细胞、胰腺腺泡）是 3D——3D 螺旋波 + 卷波（Ch6 的工具）如何应用于钙信号？这尚未充分研究。

个人反思与批判性分析

建模哲学层面：

Ch7 是 I 卷中应用最广的章节——钙信号传导涉及几乎所有生理过程。建模哲学上有几个特点：

  1. "工具箱"式建模（§7.1 末段, Berridge 2003）：每个细胞选择"工具箱"中的特定组件，建模也是组合式——这是模块化建模的早期范例。

  2. 多时间尺度 + QSSA（§7.2.5 末段 + §7.2.7）：从 8 态 ODE 到 2 变量 ODE，从 ODE 到 \(c_t\) + \(c\) 双时间尺度——反复使用 Ch1.6.3 + Ch2 的工具。

  3. 随机 + 确定性双轨（§7.6）：Falcke 2003 揭示确定性极限环 ≠ 唯一机制——随机触发可以"创造"准周期振荡。这是对"系统生物学 = ODE"的根本挑战。

  4. 多尺度均匀化（§7.8）：从微观结构（μm）到宏观现象（mm）需要有效参数。这是 Ch12 心肌组织双域方程的数学预备。

与后续章节的衔接：

  • Ch8（细胞间通讯）：细胞间钙波 + gap junction 介导的扩散——是 §7.7 的延伸。
  • Ch9（神经内分泌）：胰岛 β 细胞的簇发放（bursting）= Ch5 兴奋性 + Ch7 慢 Ca²⁺ 反馈。\(c_t\) 慢 + \(c\) 快的双时间尺度（§7.2.7）在这里是关键。
  • Ch10（细胞功能调控）：基因调控网络的双稳态 = 慢变量立方几何——与 §7.2 钙振荡的"通用结构"同源。
  • Ch11（循环系统）：血管内皮的钙信号传导——§7 的所有工具。
  • Ch12（心脏电生理）：心肌细胞 L 型 Ca²⁺ 通道 + RyR + SR 释放 = §7.1 的工具箱 + §7.5 离散钙源 + §7.8 多尺度均匀化。Ch12 是 Ch7 的直接应用。
  • Ch15（肌肉）：肌细胞的兴奋-收缩耦合 = HH 动作电位 → Ca²⁺ 释放 → 横桥循环。
  • Ch17（肾）：肾近端小管的钙信号。
  • Ch19（视网膜）：光感受器适应 = Ca²⁺ 反馈。

值得重新推导的方程：

  • (7.16) IP3R QSSA 化简：从 8 态 ODE (7.11) 出发，对激活和 IP3 结合做 QSSA，得到 2 变量系统。这是 Ch1.6.3 QSSA 在生物物理的经典应用。
  • (7.32) 简单钙振荡模型：从 (7.3) 出发选 IPR + SERCA 两个通量 + Hill 形式 (7.34) + 8 态 QSSA (7.35)。2 变量系统有 N 形零线 + Hopf。
  • (7.78) 有效扩散系数：从 (7.70) 出发对快速缓冲做 QSSA → 非线性反应-扩散方程 → 有效扩散 (7.78)。关键是移动缓冲的"载体"效应。
  • (7.89) 行波存在条件：从 (7.87) 出发用 Ch6 的双稳态方程理论 → 固定缓冲（\(D_b = 0\)）不影响行波存在性。
  • (7.121) FDF 速度公式：从 (7.116) 出发用 Green 函数 (7.118) → 阈值条件 (7.119) → 周期性稳态 → 速度条件 (7.121)。
  • (7.180) 有效扩散张量：从 (7.166) 出发做两尺度渐近展开 → 一致性条件 → 有效张量。关键是周期性微观结构。

对作者建模选择的评价：

  • 优点：（1）IP3R 8 态 + 6 态的对比清晰——展示了建模的演进（8 态太精细 / 6 态更适合工程）；（2）钙振荡的相平面分析与 Ch5 HH 完美呼应——"通用结构"的强证据；（3）Friel 1995 模型的"理论 + 实验"双向循环是 Ch1-7 一直强调的典范；（4）钙缓冲（§7.4）的"动力学"洞见是 §7.4 的核心；（5）FDF 模型（§7.5）展示了离散钙源的独特性；（6）随机钙火花（§7.6）提出了"振荡的确定性解释可能错"的重要警示；（7）细胞间波（§7.7）的三种机制（被动/部分再生/相位波）清晰分类；（8）多尺度均匀化（§7.8）是非均匀介质的关键工具。
  • 不足：（1）§7.2.5 末段对"8 态 vs 6 态"的边界讨论不充分——何时用哪个没有清晰判据；（2）§7.6.2 末段对"随机 vs 确定性"的边界承认"implications have yet to be fully digested"——这是当前最深刻的开放问题；（3）§7.4 末段 BAPTA vs EGTA 的实验现象没有详细模型解释；（4）§7.5 FDF 模型与连续方程的失配在原版明显（速度 ∝ \(D_c\) 而非 \(\sqrt{D_c}\)），加入 uptake 才修复——这暗示 FDF 模型的原始形式有根本缺陷；（5）§7.7 末段"机制争议"承认"much remains to be discovered"——细胞间波机制仍未完全理解。

若能与作者对话，我会问：

  1. 8 态 vs 6 态的边界：在什么实验数据下，应该选 8 态（精细门控）还是 6 态（饱和结合）？是否有定量的判据？
  2. 随机 vs 确定性的边界：Falcke 2003 的"随机振荡"是否在所有细胞类型中都重要？还是只在"小簇"细胞（如早期胚胎细胞）中？
  3. BAPTA vs EGTA 的"巴尔干化 vs 全球化"：Dargan-Parker 2003 的实验现象是否已被完整建模？如果有，关键参数是什么？
  4. 细胞间钙波与神经系统信息处理：Ch7 §7.7 末段提到"细胞间 Ca²⁺ 波在大脑信息处理中起重要作用"。这个论断的实验证据有多强？建模是否能验证？
  5. 多尺度均匀化的现代发展：§7.8 附录给出了经典的多尺度均匀化。现代发展（如随机均匀化、非局部均匀化、机器学习辅助的微观结构重构）是否能给出更准确的 \(D_{eff}\)？

重要参考文献

[X1] Berridge MJ, Bootman MD, Roderick HL. Calcium signalling: dynamics, homeostasis and remodelling. Nat Rev Mol Cell Biol 2003;4:517–529. （§7.1 钙工具箱的综述。）

[X2] Berridge MJ. The AM and FM of calcium signalling. Nature 1997;386:759–760.

[X3] Carafoli E. Calcium signaling: A historical account. Biol Res 2002;35:107–111.

[X4] Hofer AM, Brown EM. Extracellular calcium sensing and signalling. Nat Rev Mol Cell Biol 2003;4:530–538.

[X5] Webb SE, Miller AL. Calcium signalling during embryonic development. Nat Rev Mol Cell Biol 2003;4:539–551.

[X6] Orrenius S, Zhivotovsky B, Nicotera P. Regulation of cell death: the calcium-apoptosis link. Nat Rev Mol Cell Biol 2003;4:552–565.

[X7] MacLennan DH, Kranias EG. Calcium: a life and death signal. Nature 2003;408:641–645.

[X8] De Young GW, Keizer J. Single inositol 1,4,5-trisphosphate receptor channel properties. Proc Natl Acad Sci USA 1992;89:100–115. （8 态 IP3R 模型原始论文——§7.2.5。）

[X9] Sneyd J, Dufour J-F. A dynamic model of the type-2 inositol trisphosphate receptor. Proc Natl Acad Sci USA 2002;99:2398–2403. （6 态饱和结合模型。）

[X10] Li Y-X, Rinzel J. Equations for InsP3 receptor-mediated [Ca²⁺]i oscillations derived from a detailed kinetic model: a Hodgkin-Huxley like formalism. J Theor Biol 1994;166:461–473. （IP3R QSSA 化简的另一种形式。）

[X11] Keizer J, De Young GW. Simplification of a realistic model of IP3-induced Ca²⁺ oscillations. J Theor Biol 1994;166:431–442.

[X12] Tang Y, Stephenson JL, Othmer HG. Simplification and analysis of models of calcium oscillations based on IP3-sensitive calcium channel kinetics. Biophys J 1996;70:246–263.

[X13] Sneyd J, Falcke M. Insights into the regulation of calcium oscillations by the cADPR pathway. J Theor Biol 2005;233:51–61.

[X14] Sneyd J, et al. A comparison of three models of the inositol trisphosphate receptor. Prog Biophys Mol Biol 2004;89:136–143. （8 态 vs 6 态 vs 简化模型的对比。）

[X15] MacLennan DH, et al. A model of the SERCA pump. J Biol Chem 1997;272:28815–28818. （SERCA 6 态泵模型——§7.2.4。）

[X16] Lytton J, et al. Functional comparisons between isoforms of the sarcoplasmic or endoplasmic reticulum family of calcium pumps. J Biol Chem 1992;267:14483–14489.

[X17] Friel DD. [Ca²⁺]i oscillations in sympathetic neurons: an experimental test of a theoretical model. Biophys J 1995;68:1752–1766. （牛蛙交感神经元 CICR 模型——§7.2.9。）

[X18] Bezprozvanny I, Watras J, Ehrlich BE. Bell-shaped calcium-response curves of Ins(1,4,5)P3-gated and calcium-gated channels from endoplasmic reticulum of cerebellum. Nature 1991;351:751–754. （8 态模型的实验数据来源。）

[X19] Atri A, Amundson J, Clapham D, Sneyd J. A single-pool model for intracellular calcium oscillations and waves in Xenopus laevis oocytes. Biophys J 1993;65:1727–1739. （Xenopus 螺旋波模型——Fig. 7.3B。）

[X20] Lechleiter JD, Clapham DE. Molecular mechanisms of intracellular calcium excitability in Xenopus oocytes. Cell 1992;69:283–294. （Xenopus 钙波的实验发现——Fig. 7.3A。）

[X21] Marchant J, Parker I. Role of elementary Ca²⁺ puffs in generating repetitive Ca²⁺ transients. EMBO J 2001;20:65–77. （钙火花的实验基础。）

[X22] Cheng H, Lederer WJ, Cannell MB. Calcium sparks: elementary events underlying excitation-contraction coupling in heart muscle. Science 1993;262:740–744. （钙火花的发现——心肌 RyR。）

[X23] Smith GD, et al. Asymptotic analysis of buffered calcium diffusion near a point source. SIAM J Appl Math 2001;61:1816–1838. （钙缓冲的渐近分析。）

[X24] Wagner J, Keizer J. Effects of rapid buffers on Ca²⁺ diffusion and Ca²⁺ oscillations. Biophys J 1994;67:447–456. （移动缓冲的"载体"效应。）

[X25] Irving M, et al. A new method of measuring oxygen diffusion and consumption in skeletal muscle. News Physiol Sci 1990;5:269–272.

[X26] Falcke M. Reading the patterns in living cells — the physics of Ca²⁺ signaling. Adv Phys 2004;53:255–440. （钙动力学的最完整综述（200 页）——Falcke 本人评论。）

[X27] Pearson JE, Ponce-Dawson S. Simple model of a fire-diffuse-fire model. SIAM J Appl Math 1998;58:1667–1681. （fire-diffuse-fire 模型原始论文。）

[X28] Keizer J, Smith GD, Ponce-Dawson S, Pearson JE. Saltatory propagation of Ca²⁺ waves by Ca²⁺ sparks. Biophys J 1998;75:595–600.

[X29] Coombes S. The effect of noise on phase waves. SIAM J Appl Dyn Syst 2001;1:417–438.

[X30] Falcke M, Hudson JL, Camacho P, Lechleiter JD. Impact of mitochondrial Ca²⁺ cycling on pattern formation and stability. Biophys J 1999;77:37–44.

[X31] Sneyd J, et al. A model of calcium waves in pancreatic acinar cells. Biophys J 1995b;69:1718–1728. （胰腺腺泡钙波。）

[X32] Sneyd J, et al. A model of the mechanism of intercellular calcium wave propagation in epithelial cells. J Theor Biol 1994;170:131–146. （§7.7.1 机械刺激模型。）

[X33] Sneyd J, et al. Intercellular calcium waves mediated by diffusion of inositol (1,4,5)-trisphosphate: a two-dimensional model. Am J Physiol 1995a;268:C1537–C1545.

[X34] Sanderson MJ, Charles AC, Boitano S, Dirksen ER. Mechanical stimulation and intercellular communication increases intracellular Ca²⁺ in epithelial cells. Cell Regul 1990;1:585–596. （§7.7 胞间波的早期工作。）

[X35] Höfer T, Politi A, Heinrich R. Intercellular Ca²⁺ wave propagation mediated by gap junctions. BioSystems 2002;67:67–74. （§7.7.2 部分再生。）

[X36] Dupont G, et al. Calcium dynamics in hepatocytes: experiments and models. J Biol Phys 2000;26:317–326. （§7.7.3 肝细胞相位波。）

[X37] Höfer T. Model of intercellular calcium signalling in epithelial and nonexcitable cells. BioSystems 1999;50:73–82.

[X38] Dupont G, Swillens S, Tordjmann T, Combettes L. Hierarchy of mechanisms involved in generating Ca²⁺ oscillations. Biochim Biophys Acta 2007;1773:1311–1322.

[X39] Falcke M. Deterministic and stochastic models of intracellular Ca²⁺ waves. N J Phys 2003a;5:96.1–96.28. （§7.6.2 随机钙波。）

[X40] Shuai JW, Jung P. Optimal intracellular calcium signaling. Phys Rev Lett 2002a;88:068102.

[X41] Fox RF, Lu YN. Emergent collective behavior in large numbers of globally coupled independently stochastic ion channels. Phys Rev E 1994;49:3421–3431. （Hodgkin-Huxley 的 Langevin 扩展——Ch7 §7.6.1 的方法学参考。）

[X42] Keener J. Stochastic and deterministic methods in mathematical biology. Stoch Dyn 2006;6:465–484. （§7.6.2 末段的引用。）

[X43] Goldbeter A, Dupont G, Berridge MJ. Minimal model for signal-induced Ca²⁺ oscillations and for their frequency encoding through protein phosphorylation. Proc Natl Acad Sci USA 1990;87:1461–1465. （Exercise 3 的两池模型。）

[X44] Tsien RW, Tsien RY. Calcium channels, stores, and oscillations. Annu Rev Cell Biol 1990;6:715–760. （早期综述。）

[X45] Dupont G, Erneux C. Simulations of the effects of inositol 1,4,5-trisphosphate 3-kinase and 5-phosphatase activities on Ca²⁺ oscillations. Cell Calcium 1997;22:321–331.

[X46] Politi A, et al. Models of IP3 and Ca²⁺ oscillations: frequency encoding and specificity. Cell Calcium 2006;40:157–165.

[X47] Bennett MVL, et al. Gap junctions: new tools, new answers, new questions. Neuron 1991;6:305–320.

[X48] Goel P, Sneyd J, Nirody MK. A homogenization approach to effective diffusion in heterogeneous excitable media. Phys Rev E 2006;74:021911. （§7.8 有效扩散张量的具体计算。）

[X49] Bensoussan A, Lions J-L, Papanicolaou G. Asymptotic Analysis for Periodic Structures. North-Holland, 1978. （多尺度均匀化的数学基础——§7.8 附录。）

[X50] Wagner J, Keizer J, Bhatt S. A theoretical analysis of calcium waves in dictyostelium discoideum. J Cell Biol 1990;111:2715–2727.

[X51] MacLennan DH, Rice WJ, Green NM. The mechanism of Ca²⁺ transport by sarco(endo)plasmic reticulum Ca²⁺-ATPases. J Biol Chem 1997;272:28815–28818.