跳转至

第9章：神经内分泌细胞

作者

本章作者与 Ch1-Ch8 相同——James P. Keener 与 James Sneyd。本章是 I 卷"细胞生理学"中多振荡机制耦合的典范章——把 Ch5 HH 兴奋性、Ch7 钙振荡、Ch8 突触传递、Ch1 糖酵解振荡整合到神经内分泌细胞的"簇发放（bursting）"。本章角色：（1）胰岛 β 细胞的 Chay-Keizer 1983 模型（§9.1.1）——Rinzel 1985 双稳态 + 相平面分析；（2）Phantom bursting 模型（§9.1.2, Bertram-Sherman 2004）——ER Ca²⁺ 作为慢变量；（3）Compound bursting + 糖酵解（§9.1.3, Wierschem-Bertram 2004）——Ch1.6 糖酵解振荡与 Ch5/7 电生理的耦合；（4）细胞簇的 supercell 模型（§9.1.4, Atwater 1983, Sherman 1988）——Gap junction 介导的同步 + 随机通道共享；（5）Aplysia R-15 的 Plant 1981 抛物型 bursting（§9.1.6）——没有双稳性，两个慢变量；（6）下丘脑-垂体促性腺细胞（§9.2, Li 1997）——空间梯度 + 激素反应。

内容概述

本章围绕"神经内分泌细胞如何产生簇发放"这一核心问题展开。

§9.0 总体框架：神经内分泌细胞的3 大特点——兴奋性、簇发放（不是单动作电位或单周期串）、Ca²⁺ 调节。簇发放的 3 大类（Rinzel 1987, Bertram 1995, Izhikevich 2000）：Type I（双稳性 + 滞后环，β 细胞），Type II（慢振荡驱动，Aplysia R-15），Type III（亚临界 Hopf）。

§9.1.1 Chay-Keizer 1983 β 细胞模型：(9.1)-(9.5) 给出 4 电流 + Ca²⁺ 方程。Rinzel 1985 双稳态分析（§9.1.1 后段）：忽略 m, h 的快动力学 → 简化 3 变量（V, n, c）。核心几何（Fig. 9.3）：快子系统（V, n）的 Z 形稳态曲线。当 c 在中间范围 → 3 稳态（稳定 + 鞍 + 不稳定），其中不稳定被稳定极限环包围。当 c 增大 → 极限环与鞍点碰撞 → 同宿分岔（Fig. 9.3B）→ 没有极限环。分岔图（Fig. 9.4A）：HB (Hopf) + HC (homoclinic) + SN (saddle-node)。

簇发放机制（Fig. 9.4B）：快慢分离——c 慢变 → V 在 Z 形曲线上"切换"两分支。c 上升时：V 在下分支（静息）→ 鞍-结分岔 SN → V 跳到上分支的极限环（active phase）→ 同宿分岔 HC → V 跳回下分支（silent phase）。这是 Type I 簇发放的通用机制。

§9.1.2 Phantom Bursting（Bertram-Sherman 2004）：实验观察（Fig. 9.5, Zhang 2003）：Ca²⁺ 上升几乎与 V 一样快——Chay-Keizer 的"c 是慢变量"假设失效。新假设（Chay 1996）：ER Ca²⁺ \(c_e\) 是慢变量——\(c_e\) 在 active phase 上升，在 silent phase 下降——低通滤波 cytoplasmic Ca²⁺。模型（9.9)-(9.22）：4 电流（Ca²⁺ + K⁺ + K(Ca) + K(ATP)）+ 2 Ca²⁺ 变量（c + \(c_e\)）+ 2 门控（n + Ca²⁺-敏感 K 通道激活 \(\omega\)）。ATP 敏感 K⁺ 通道（9.14）是葡萄糖调节器——ATP/ADP 比低 → K(ATP) 开放 → 超极化 → 无 burst；高 → 关闭 → burst。

Fast Bursting（Fig. 9.6）：\(c_e\) 变化小，c 是有效慢变量——等效 Chay-Keizer。Medium Bursting（Fig. 9.9）：\(c_e\) 变化大，\(c_e\) 是慢变量——c 零线随 \(c_e\) 移动 → burst period 拉长。

§9.1.3 Compound Bursting + 糖酵解（Wierschem-Bertram 2004）：关键观察（§9.1.3 开头）：簇发放周期可2-4 分钟（比 medium 10-60 s 还长）——需要"更慢"的振荡源。假说：糖酵解振荡（Ch1.6 提到的 Goldbeter-Lefever 模型）→ ATP/ADP 慢振荡 → \(g_{K,ATP}\) 调制 → burst period 调制。模型（9.24)-(9.28）：糖酵解 2 变量（[ATP], [ADP]）+ 简化 Goldbeter-Lefever + Chay-Keizer 简化。Compound bursting（Fig. 9.10）：burst of bursts——每个 active phase 本身由多个子 burst 组成。机制：[ATP] 在 [ATP] 低时 Z 形零线在 lower branch（无 burst）→ [ATP] 中时 Z 形零线在中分支（burst）→ [ATP] 高时 Z 形零线在 upper branch（持续 active）。[ATP] 振荡 → burst 周期调制。

§9.1.4 Cell Cluster 的"Supercell"模型（Atwater 1983, Sherman 1988, Sherman-Rinzel 1991）：实验事实（Fig. 9.11）：孤立 β 细胞不规则 burst（Fig. 9.11C），成簇 β 细胞规则 burst（Fig. 9.11A, B）。机制：每个 K(Ca) 通道的"大电导"（比单通道 ~10×）+ 小通道数（每个细胞 ~100-1000）→ 单一通道的开放对 V 影响显著 → 随机通道开口 → V 随机扰动。簇内耦合：N_cell 个细胞 + gap junction → V 同步 + 通道电流被 N_cell 平均化 → 每通道的相对影响 1/N_cell → 宏观确定性。

Supercell 公式（9.42)-(9.43）：\(N_{cell} C_m dV/dt = -N_{cell}(I_{Ca} + I_K) - \hat g \sum N_i^o (V - V_K)\) → 化简为 1 大细胞，有效 K(Ca) 电导降低 \(1/N_{cell}\)。多细胞模型（9.44）：有限 gap junction 电导 \(g_c\) → (9.46) 显示 supercell 极限 \(g_c \to \infty\)。

关键发现（§9.1.4 末段 Fig. 9.14）：存在一个最优耦合强度让 burst period 最大化（Sherman 1994 详细分析）——强耦合和弱耦合都不如中间最优。

§9.1.5 Hindmarsh-Rose 多项式模型（9.47)-(9.57）：简化 FHN + 第三慢变量 z：\(dx/dT = y - x^3 + 3x^2 + I - z\), \(dy/dT = 1 - 5x^2 - y\), \(dz/dT = r[s(x-x_1) - z]\)。关键参数 r = 0.001, s = 4 → bursting。Liénard 形式（Exercise 8）：\(d^2x/dT^2 + F(x) dx/dT + G(x, z) = -\epsilon H(x, z)\)——van der Pol 振子 + 慢变量。

§9.1.6 Aplysia R-15 抛物型 bursting（Plant 1981）：关键差异 vs Chay-Keizer：没有双稳性，两个慢变量（x 和 c）。抛物型：burst period 在 burst 中间最短（远离同宿分岔），burst 头尾最长（接近同宿分岔）——period 是 t 的抛物线。分类（§9.1.6 末段）：Rinzel 1987 把 bursting 分为 3 类，Izhikevich 2000 扩展到大量可能机制。

§9.2 下丘脑-垂体促性腺细胞（Li 1997）：关键创新（§9.2 末段）：球对称空间梯度——Ca²⁺ 在膜附近高，胞内低。膜 + IP3R + 缓冲的完整模型（9.67)-(9.72）。关键预测（Fig. 9.22）：GnRH 刺激 → IP3 增加 → 大 Ca²⁺ 振荡 → \(I_{K,Ca}\) 强烈激活 → V spike 频率下降（phase T-1）→ Ca²⁺ 整体下降（pmp 漏）→ \(c_e\) 下降 → Ca²⁺ 振荡变小变快（phase SS-2）→ V spike 频率回升 → Ca²⁺ 内流 → \(c_e\) 恢复 → 基础状态。

主要数学结果： 1. Chay-Keizer 1983 模型 + Rinzel 1985 快慢分离（9.6)-(9.8） 2. Z 形稳态 + HB + HC + SN 三分岔（图 9.4） 3. Phantom bursting（Bertram-Sherman 2004）——ER Ca²⁺ 慢变量 4. Compound bursting（Wierschem-Bertram 2004）——糖酵解振荡 + 简化 FHN 5. Supercell 公式（9.42)-(9.43）—— 通道共享 6. Sherman 1994 最优耦合强度 7. Hindmarsh-Rose 多项式 + Liénard 形式 8. Plant 1981 抛物型 bursting（没有双稳性） 9. Li 1997 球对称 Ca²⁺ 梯度 + IP3R 整合 10. bursting 分类（Rinzel 1987, Izhikevich 2000）

前置知识：Ch1.6（Goldbeter-Lefever 糖酵解）、Ch5（HH, FHN）、Ch7（IP3R, ER Ca²⁺, SERCA）、Ch8（gap junction）。

核心方程与概念

1. Chay-Keizer 1983 β 细胞模型

4 电流 + Ca²⁺ 方程（9.4）–(9.5）：

\[ C_m \frac{dV}{dt} = -(g_{K,Ca} + g_K)(V - V_K) - 2g_{Ca}(V - V_{Ca}) - g_L(V - V_L) \]
\[ \frac{dc}{dt} = f(-k_1 I_{Ca} - k_c c) \]

离子电导（9.1）–(9.3）： - \(g_{K,Ca} = \bar g_{K,Ca} c/(K_d + c)\)——Ca²⁺ 激活 K⁺ - \(g_K = \bar g_K n^4\)——电压门控 K⁺（Ch5 HH \(n\) 门控） - \(g_{Ca} = \bar g_{Ca} m^3 h\)——电压门控 Ca²⁺（用 HH Na⁺ 形式）

关键参数（Table 9.1）：\(\bar g_{K,Ca} = 0.02\), \(\bar g_K = 3\), \(\bar g_{Ca} = 3.2\), \(\bar g_L = 0.012\) mS/cm², \(K_d = 1\) μM, \(f = 0.007\)（f 很小让 Ca²⁺ 慢）, \(k_c = 0.02\) ms⁻¹。

2. Rinzel 1985 双稳态分析

简化模型（9.6）–(9.8）：忽略 m, h 动力学 → \(I_{Ca} = \bar g_{Ca} m_\infty^3 h_\infty (V - V_{Ca})\)（代数而非 ODE）。

Z 形稳态曲线（Fig. 9.4A）： - 下分支：稳定（静息，\(V\) 低） - 中分支：不稳定（鞍点） - 上分支：不稳定节点（周围有稳定极限环）

分岔： - Hopf (HB)：\(c_{HB}\) 处上分支失稳 → 稳定极限环出现 - Homoclinic (HC)：\(c_{HC}\) 处极限环与鞍点碰撞 → 极限环消失 - Saddle-Node (SN)：Z 形曲线在 fold 处出现/消失稳态

Bursting 机制（Fig. 9.4B）：c 慢变 → V 在 Z 形曲线上切换。c 减小时：下分支 → SN → 上分支（极限环）→ active；c 增大时：上分支（极限环）→ HC → 下分支 → silent。

3. Phantom Bursting（Bertram-Sherman 2004）

4 电流 + 2 Ca²⁺ 变量（9.9）–(9.22）：

\[ I_{Ca} = g_{Ca} m_\infty(V)(V - V_{Ca}), \quad I_{K,Ca} = g_{K,Ca} \omega(c)(V - V_K) \]
\[ I_{K,ATP} = g_{K,ATP}(V - V_K), \quad \omega(c) = \frac{c^5}{c^5 + k_D^5} \]
\[ \frac{dc}{dt} = f_{cyt}(-\alpha I_{Ca} - J_{pm} + J_{leak} - J_{serca}) \]
\[ \frac{dc_e}{dt} = -\gamma f_{er}(J_{leak} - J_{serca}) \]

K(Ca) 通道的 Hill 形式（9.17）\(\omega(c) = c^5/(c^5 + k_D^5)\)——5 个 Ca²⁺ 合作结合。

K(ATP) 通道（9.14）：ATP 调节——ATP 高 → K(ATP) 关闭 → 去极化 → burst；ATP 低 → K(ATP) 开放 → 超极化 → 无 burst。这是"葡萄糖传感器"的机制。

关键（§9.1.2 Fig. 9.7）：fast bursting：\(c_e\) 变化小，c 是有效慢变量。Medium bursting（Fig. 9.8）：\(c_e\) 变化大，\(c_e\) 是慢变量，c 零线随 \(c_e\) 移动 → burst period 拉长。

4. Compound Bursting + 糖酵解（Wierschem-Bertram 2004）

完整模型（9.24）–(9.28）：

糖酵解 2 变量（9.24）–(9.25）：

\[ \tau_c \frac{d[ATP]}{dt} = v_1 - F([ATP], [ADP]), \quad \tau_c \frac{d[ADP]}{dt} = F([ATP], [ADP]) - v_2 [ADP] \]

其中 \(F([ATP], [ADP]) = [ATP](1 + k_{ADP} [ADP])^2\)（Ch1.6 简化 Goldbeter-Lefever）。

Compound bursting（Fig. 9.10）：[ATP] 振荡 → Z 形稳态曲线的水平移动（§9.1.3 Fig. 9.10a, b, c）。当 [ATP] 低：Z 形曲线 → c 零线在下分支 → 无 burst（a）。当 [ATP] 中：Z 形曲线在中分支 → c 零线穿越稳定极限环 → burst（b）。当 [ATP] 高：Z 形曲线在上分支 → 持续 active（c）。Compound bursting = [ATP] 在三个区域的循环。

5. Supercell 通道共享（Atwater 1983, Sherman 1988）

单细胞 vs 簇：(9.30) 通道 \(C \xrightleftharpoons[k_-]{k_+}\) \(O\)（Hill 形式 (9.29)）。关键参数（Table 9.4）：单 K(Ca) 通道电导 \(\hat g = 30\) nS，比 \(g_K = 2.5\), \(g_{Ca} = 1.4\) nS 大 10 倍。

Supercell 公式（9.42）–(9.43）：

\[ C_m \frac{dV}{dt} = -(I_{Ca} + I_K) - \bar g_{K,Ca} p (V - V_K), \quad p = \frac{1}{N_{cell} \bar N} \sum_{i=1}^{N_{cell}} N_i^o \]

——p 是全簇的开放比例。方差：\(Var(p) = Var(N^o)/(N_{cell} \bar N)\)——随 \(N_{cell}\) 减小。

最优耦合强度（§9.1.4 末段, Sherman 1994）：存在一个最优 \(g_c\) 让 burst period 最大化——强耦合和弱耦合都不如。

6. Hindmarsh-Rose 多项式

(9.55)–(9.57）：

\[ \frac{dx}{dT} = y - x^3 + 3x^2 + I - z, \quad \frac{dy}{dT} = 1 - 5x^2 - y, \quad \frac{dz}{dT} = r[s(x - x_1) - z] \]

关键参数：\(r = 0.001\)（慢变量 z 的时间尺度）, \(s = 4\)。

Liénard 形式（Exercise 8）：

\[ \frac{d^2x}{dT^2} + F(x) \frac{dx}{dT} + G(x, z) = -\epsilon H(x, z), \quad \frac{dz}{dT} = \epsilon H(x, z) \]

——van der Pol 振子 + 慢调制。

7. Plant 1981 抛物型 bursting

3 慢变量（9.58）–(9.62）：

快子系统（V, h, n）：标准 HH Na⁺ + K⁺ + K(Ca) + 漏。

慢子系统（x, c）：\(\tau_x dx/dt = x_\infty(V) - x\), \(dc/dt = f(k_1 x (V_{Ca} - V) - c)\)。

关键差异（Fig. 9.19）：没有双稳性——同宿分岔 HC 与鞍-结分岔 SN1 重合。机制：两个慢变量（x, c）的振荡 → 慢轨道穿越 HC 线 → burst。

Period 抛物线（§9.1.6 Fig. 9.20B）：burst 头尾 period 长（接近 HC），burst 中间 period 短（远离 HC）。

8. Li 1997 促性腺细胞球对称模型

完整方程（9.67）–(9.72）：

\[ \frac{\partial c}{\partial t} = D \nabla^2 c + J_{IPR} - J_{serca}, \quad \frac{\partial c_e}{\partial t} = D_e \nabla^2 c_e - \gamma(J_{IPR} - J_{serca}) \]

边界条件（9.71）：

\[ D \frac{\partial c}{\partial r}\bigg|_{r=R} = -\alpha(I_{Ca,T} + I_{Ca,L}) - J_{pm} \]

关键（Fig. 9.22）：5 阶段（SS-1, T-1, T-2, SS-2, basal）——ATP/ADP + IP3 + Ca²⁺ + V 的多变量耦合。

关键结论

  1. 簇发放 = 快慢分离 + 双稳态（§9.1.1, Rinzel 1985）：2 个快变量（V, n）+ 1 个慢变量（c）。c 慢变 → V 在 Z 形曲线上切换——Type I 簇发放的通用机制。

  2. Chay-Keizer 1983 是 β 细胞簇发放的"标准模型"（§9.1.1）：4 电流 + Ca²⁺ 方程——Rinzel 1985 的相平面分析揭示几何机制。这是"快慢分离"在电生理的经典胜利。

  3. Phantom Bursting 揭示 ER Ca²⁺ 是真正慢变量（§9.1.2, Bertram-Sherman 2004）：实验发现 cytoplasmic Ca²⁺ 几乎与 V 一样快（Zhang 2003）→ Chay-Keizer 的"c 慢"假设失效 → ER Ca²⁺ 才是慢变量。低通滤波 cytoplasmic Ca²⁺是机制。

  4. Compound Bursting = 糖酵解振荡 + 简化 FHN（§9.1.3, Wierschem-Bertram 2004）：最慢的 burst 周期对应糖酵解振荡（Ch1.6 Goldbeter-Lefever）。Ch1 + Ch5 + Ch7 的"3 章合一"。

  5. K(ATP) 通道是"葡萄糖传感器"（§9.1.2 末段）：ATP 高 → K(ATP) 关闭 → burst；ATP 低 → K(ATP) 开放 → 无 burst。这是 β 细胞响应血糖的核心机制——临床降糖药（磺酰脲类）作用于 K(ATP)。

  6. Supercell 通道共享解释"单细胞不规则 + 簇内规则"（§9.1.4, Atwater 1983, Sherman 1988）：单 K(Ca) 通道电导大 10× + 通道数小→ 单细胞随机扰动大 → 不规则 burst。簇内 N_cell 平均化→ 规则 burst。这是"统计平均 = 确定性"在生物系统的范例。

  7. 存在最优耦合强度（§9.1.4 末段, Sherman 1994）：强耦合和弱耦合都不如中间最优——这是"非单调耦合"现象。真实胰岛的耦合强度可能就在这个最优点附近。

  8. Plant 1981 抛物型 bursting 没有双稳性（§9.1.6）：2 个慢变量（x, c）→ 慢振荡穿越 HC 线 → burst。Period 是 t 的抛物线——"parabolic bursting"名称的来源。这是 Type II 簇发放的经典模型。

  9. Hindmarsh-Rose 多项式是"通用 bursting 模型"（§9.1.5）：3 变量 + Liénard 形式。通过参数 r, s, I 的调整可以产生多种 bursting 模式（square-wave, tapered, 等等）——Izhikevich 2000 给出系统分类。

  10. Bursting 分类（§9.1.6 末段, Rinzel 1987, Izhikevich 2000）：Type I（双稳性 + 滞后）、Type II（两个慢变量）、Type III（亚临界 Hopf）——3 大类，但具体生物实现可能用更复杂的机制。"通用模型 vs 现象模型"的张力。

挑战和开放性问题

  1. ATP 敏感 K⁺ 通道的精确模型（§9.1.2 末段）：Ch9 把 K(ATP) 简化为常数电导——真实 K(ATP) 有复杂的 ATP 依赖动力学。这能否解释 Type III 簇发放？

  2. Compound bursting 的代谢耦合（§9.1.3 末段, Bertram 2007）：Smolen 模型 + Magnus-Keizer 线粒体代谢——比 Wierschem-Bertram 复杂得多。糖酵解 vs 线粒体哪个对 β 细胞簇发放更重要？

  3. 最优耦合强度的机制（§9.1.4 末段, Sherman 1994）：为何 burst period 与 \(g_c\) 关系是非单调的？这是"同步 vs 异步"的权衡吗？

  4. Aplysia R-15 抛物型 bursting 在其他细胞中的普适性（§9.1.6 末段）：只有 Plant 1981 一个详细模型——是否需要更多实现？

  5. Bursting 频率编码（§9.1 开头）：burst 频率携带信息——不同频率触发不同基因表达、不同分泌速率。这与 spike 频率编码的关系是什么？

  6. 多振荡机制的精确耦合（§9.1.3 + §9.2）：糖酵解 + 线粒体 + ER Ca²⁺ + cytoplasmic Ca²⁺ + V——5 个时间尺度。何时哪个主导？这是"多时间尺度"建模的核心挑战。

  7. 下丘脑-垂体其他细胞类型（§9.2 开头）：促肾上腺皮质细胞（LeBeau 1997, 1998）、加压素神经元（Roper 2003, 2004）——这些模型的精确比较？

  8. Spatial 梯度的作用（§9.2.1 末段, Li 1997）：c 膜附近高，胞内低——这让 \(I_{K,Ca}\) 用 \(c_R\)（膜附近）而非平均 c。这个"空间局部化"在简单模型（β 细胞用均匀 c）中被忽略——这是否是 Aplysia R-15 vs β 细胞的本质差异？

  9. 随机簇发放（§9.1.4 + §9.1.6）：随机通道 + 簇内耦合让"单细胞 → 簇"是确定性增强。但单细胞本身的随机性如何？这能否解释 β 细胞的"不规则 burst"？

  10. Clinical 转化（§9.1.2 末段）：糖尿病的根因之一是β 细胞功能失调——簇发放异常。磺酰脲类作用于 K(ATP) 通道。更精细的"代谢-电生理"模型能否指导个体化治疗？

个人反思与批判性分析

建模哲学层面：

Ch9 是 I 卷"建模应用的高峰"——把 Ch1（糖酵解）+ Ch5（HH/FHN）+ Ch7（钙振荡）+ Ch8（gap junction）整合到一个具体的生物系统（神经内分泌细胞）。几个重要特点：

  1. Rinzel 1985 的快慢分离是 Ch5 思想的直接应用（§9.1.1）：HH 4 维 ODE → 简化为 3 维（V, n, c）→ 进一步快慢分离到 (V, n) + c。这是"Ch5 FHN 简化"思想在生物物理的具体胜利。

  2. Bertram-Sherman 2004 phantom bursting 是"实验驱动的建模"（§9.1.2）：Zhang 2003 的新实验数据让 Chay-Keizer 模型失败 → 重新假设 ER Ca²⁺ 慢变量 → 新模型。这是"模型必须跟随实验"的范例。

  3. Wierschem-Bertram 2004 compound bursting 是"跨章整合"（§9.1.3）：Ch1 糖酵解 + Ch5/7 簇发放 = 3 章合一。这是 I 卷"工具组合"的典范。

  4. Supercell 模型是"统计平均 = 确定性"在生物系统的范例（§9.1.4）：N_cell = 1 → 随机，N_cell → ∞ → 确定性。这是大数定律在细胞生物学的应用。

  5. Plant 1981 抛物型 bursting 揭示"没有双稳性也能 burst"（§9.1.6）：Chay-Keizer 不是唯一的簇发放机制——两个慢变量可以替代双稳性。这是"机制多元性"在神经内分泌细胞的体现。

  6. Hindmarsh-Rose 模型是"通用簇发放模型"（§9.1.5）：3 变量 + Liénard 形式 → 通过参数调整可产生多种 bursting。这是"模型简化的极限"。

与后续章节的衔接：

  • Ch10（细胞功能调控）：基因调控网络的双稳态 = 慢变量的"立方"几何（与 FHN 同源）。Ch9 的"双稳性 bursting"和 Ch10 的"双稳性基因表达"是同一数学结构在不同生物尺度的应用。
  • Ch12（心脏电生理）：心肌细胞的早期后去极化 (EAD)和延迟后去极化 (DAD)是另一种 bursting。Ch9 的 bursting 分析工具直接应用到 Ch12。
  • Ch16+（神经系统）：神经元的 bursting 普遍存在（Fig. 9.18）——丘脑-皮层、嗅觉、下丘脑等。Ch9 的模型是"通用神经簇发放"的基础。
  • Ch18（细胞间耦合）：胰腺 β 细胞的"簇内同步"是细胞间耦合的范例——Ch9 + Ch8 整合。

值得重新推导的方程：

  • (9.6)–(9.8) Rinzel 简化：从 Chay-Keizer (9.1)-(9.5) 出发，假设 m, h 在瞬时平衡 → 3 变量。关键是 m_∞, h_∞——Ch5.1.2 的拟合函数。
  • (9.4) Chay-Keizer 膜方程：4 电流求和 → 1 阶 ODE。关键是电流符号（Ca²⁺ 是内向电流 → \(V_{Ca} > 0\)，K⁺ 是外向 → \(V_K < 0\)）。
  • (9.42)–(9.43) Supercell 公式：从 (9.30) 单通道 Markov 出发，对 N_cell 个细胞求和 → 1 维宏观变量。关键是 \(p = \sum N_i^o / (N_{cell} \bar N)\)。
  • (9.17) K(Ca) Hill 形式：5 个 Ca²⁺ 合作结合 → \(\omega(c) = c^5/(c^5 + k_D^5)\)。Ch1.4.4 Hill 方程的简化。
  • (9.55)–(9.57) Hindmarsh-Rose：从 FHN (Ch5.2) 出发 + 慢变量 z 调制 I → 3 变量。关键是 I 替换为 I - z。
  • (9.58)–(9.62) Plant 模型：从标准 HH + K(Ca) + Ca²⁺ 慢激活 x → 抛物型 bursting。关键是 x 的慢动力学。

对作者建模选择的评价：

  • 优点：（1）Rinzel 1985 的快慢分离 + 相平面分析是经典；（2）Bertram-Sherman 2004 phantom bursting 完美展示"实验驱动建模"；（3）Wierschem-Bertram 2004 compound bursting 跨章整合清晰；（4）Sherman 1988 supercell 模型揭示"统计平均 = 确定性"；（5）Hindmarsh-Rose 多项式简洁优雅；（6）Plant 1981 抛物型 bursting 揭示"机制多元性"；（7）Li 1997 球对称 Ca²⁺ 梯度是空间建模的范例；（8）bursting 分类（Izhikevich 2000）是数学严密 + 物理直观的结合。
  • 不足：（1）Chay-Keizer 1983 在 2009 年已经"过时"——更现代的模型包括 Smolen + Bertram 2004 等；（2）Phantom bursting 的"实验驱动"虽然精彩，但ER Ca²⁺ 慢变量的精确机制仍不清楚；（3）Compound bursting 的代谢耦合是 §9.1.3 末段提到但没详细展开；（4）最优耦合强度的机制（Sherman 1994）§9.1.4 末段只一句话提及，没有详细分析；（5）Aplysia R-15 vs β 细胞的对比应该给出更系统的"机制分类"；（6）下丘脑-垂体促性腺细胞（§9.2）的 5 阶段（Fig. 9.22）虽然精彩，但没有详细方程推导。

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

  1. ER Ca²⁺ 慢变量的精确机制（§9.1.2）：phantom bursting 假设"c 是快速 + c_e 是慢"——这是否在所有 β 细胞中都成立？不同物种（小鼠 vs 人）有差异吗？
  2. Compound bursting 的代谢耦合（§9.1.3）：糖酵解 vs 线粒体哪个对簇发放更重要？Smolen 模型 + Magnus-Keizer 模型能给出代谢-电生理的完整耦合吗？
  3. 最优耦合强度（§9.1.4 末段, Sherman 1994）：为何 burst period 与 g_c 是非单调的？这是"同步增强但 burst 缩短"的权衡吗？真实胰岛的 g_c 处于最优点吗？
  4. Aplysia R-15 vs β 细胞的本质差异（§9.1.6）：球对称 Ca²⁺ 梯度（§9.2）是否让"双稳性 + 单慢变量"能产生更丰富的 bursting？这是"空间 = 额外维度"的核心论点吗？
  5. Bursting 频率编码（§9.1 开头）：burst 频率 vs spike 频率——哪个携带更多信息？真实细胞（如 GnRH 神经元）如何区分？

重要参考文献

[X1] Chay TR, Keizer J. Minimal model for membrane oscillations in the pancreatic β-cell. Biophys J 1983;42:181–190. （Chay-Keizer β 细胞模型——§9.1.1 原始论文。）

[X2] Rinzel J. Bursting oscillations in an excitable membrane model. In: Ordinary and Partial Differential Equations, Sleeman BD, Jarvis RJ (eds.), Springer Lecture Notes in Mathematics 1151, 1985, pp. 304–316. （Rinzel 1985 快慢分离——§9.1.1 的相平面分析。）

[X3] Rinzel J, Lee YS. On different mechanisms for membrane potential bursting. In: Bifurcation Analysis and Nonlinear Eigenvalue Problems, Gurel K, Rössler OE (eds.), Springer, 1986.

[X4] Rinzel J, Lee YS. Dissection of a model for neuronal parabolic bursting. J Math Biol 1987;25:653–675. （Aplysia 抛物型 bursting 的数学分析。）

[X5] Chay TR. Electrical bursting and luminal Ca²⁺ oscillation in excitable cell models. Biol Cybern 1996;75:419–431. （ER Ca²⁺ 作为慢变量——phantom bursting 假说的早期工作。）

[X6] Bertram R, Sherman A. Negative calcium feedback: the road to phantom bursting. J Biol Phys 2004a;30:133–152. （Phantom bursting 的完整论文——§9.1.2。）

[X7] Bertram R, Sherman A. A calcium-based phantom bursting model for pancreatic islets. Bull Math Biol 2004b;66:1313–1344.

[X8] Zhang M, Goforth P, Bertram R, Sherman A, Satin L. The Ca²⁺ dynamics of isolated mouse β-cells and islets: implications for mathematical models. Biophys J 2003;84:2852–2870. （c 不是慢变量的关键实验——§9.1.2 Fig. 9.5。）

[X9] Wierschem A, Bertram R. Complex bursting in pancreatic islets: a role for glycolytic oscillations. J Theor Biol 2004;228:513–521. （Compound bursting + 糖酵解——§9.1.3。）

[X10] Bertram R, Satin L, Pedersen MG, Luciani DS, Sherman A. Interaction of glycolysis and mitochondrial respiration in metabolic oscillations of pancreatic islets. Biophys J 2007a;92:1544–1555.

[X11] Atwater I, Dawson CM, Scott A, Eddlestone G, Rojas E. The nature of the oscillatory behaviour in electrical activity for pancreatic β-cell. J Horm Metab Res 1980;suppl. 10:100–107. （最早的 β 细胞 bursting 模型。）

[X12] Atwater I, Rinzel J. The β-cell electrical bursting: underlying mechanism and gap-junction coupling. J Mol Endocrinol 1986;1:59–64. （Cell cluster 的早期工作。）

[X13] Sherman A, Rinzel J, Keizer J. Emergence of organized bursting in clusters of pancreatic β-cells by channel sharing. Biophys J 1988;54:411–425. （Supercell 模型——§9.1.4。）

[X14] Sherman A, Rinzel J. Model for synchronization of pancreatic β-cells by gap junction coupling. Biophys J 1991;59:547–559. （Multicell 模型——有限 g_c。）

[X15] Sherman A. Anti-phase, asymmetric and aperiodic oscillations in excitable cells—I. Coupled bursters. Bull Math Biol 1994;56:801–835. （最优耦合强度——§9.1.4 末段。）

[X16] Tsaneva-Atanasova K, Sherman A, Li Y, Roper P, Bertram R. Mechanism of synchronization in electrically coupled pancreatic β-cells. J Theor Biol 2006;241:173–184.

[X17] Plant RE. Bifurcation and resonance in a model for bursting nerve cells. J Math Biol 1981;11:15–32. （Aplysia R-15 抛物型 bursting——§9.1.6。）

[X18] Hindmarsh JL, Rose RM. A model of neuronal bursting using three coupled first order differential equations. Proc R Soc Lond B 1984;221:87–102. （Hindmarsh-Rose 多项式——§9.1.5。）

[X19] Pernarowski M. Fast subsystem bifurcation analysis of a model of bursting nerve cells. J Math Biol 1994;32:473–480. （Hindmarsh-Rose 分岔分析。）

[X20] de Vries G. Multiple bifurcations and stability in a neural network model. PhD thesis, Utrecht University, 1995. （bursting 综述。）

[X21] Izhikevich EM. Neural excitability, spiking and bursting. Int J Bifurc Chaos 2000;10:1171–1266. （bursting 系统分类——§9.1.6 末段。）

[X22] Wang X-J, Rinzel J. Oscillatory and bursting properties of neurons. In: The Handbook of Brain Theory and Neural Networks, Arbib MA (ed.), MIT Press, 1995, pp. 686–691. （Fig. 9.18 bursting 细胞类型汇总。）

[X23] Li Y-X, Stojilković SS, Keizer J, Rinzel J. Excitation of pituitary gonadotrophs by hypothalamic releasing hormones. Cell Calcium 1994;15:39–50. （促性腺细胞模型——§9.2.1。）

[X24] Li Y-X, Rinzel J, Keizer J, Stojilković SS. Calcium oscillations in pituitary gonadotrophs: comparison of experiment and theory. Proc Natl Acad Sci USA 1995;92:11261–11265.

[X25] Li Y-X, Rinzel J, Vergara L, Stojilković SS. Spontaneous electrical and calcium oscillations in a pituitary cell line. J Biol Chem 1997;272:24310–24316. （球对称 Ca²⁺ 梯度模型——§9.2.1。）

[X26] LeBeau AP, Van Goor F, Stojilković SS, Sherman A. Modeling of membrane excitability in gonadotropin-releasing hormone-secreting hypothalamic neurons regulated by Ca²⁺-dependent and Ca²⁺-independent Ca²⁺-sensitive K⁺ channels. Biol Cybern 1997;77:167–180. （GnRH 神经元模型。）

[X27] Roper P, Callaway JC, Armstrong WE. Burst initiation and termination in phasic vasopressin cells of the rat supraoptic nucleus: a combined mathematical, electrical, and Ca²⁺ fluorescence study. J Neurosci 2004;24:4818–4831. （加压素神经元。）

[X28] Stojilković SS, Catt KJ. Calcium oscillations in anterior pituitary cells. Endocr Rev 1992;13:256–280. （垂体 Ca²⁺ 振荡综述。）

[X29] Sherman A. Theoretical aspects of synchronized bursting in β-cells. In: Pacemaker Activity and Intercellular Communication, Huizinga JD (ed.), CRC Press, 1995.

[X30] Keizer J, Smolen P. Bursting electrical activity in pancreatic β cells caused by Ca²⁺- and voltage-inactivated Ca²⁺ channels. Proc Natl Acad Sci USA 1991;88:3897–3901. （β 细胞 bursting 的非 Ca²⁺ 慢变量假说。）

[X31] Smolen P, Keizer J. Slow voltage inactivation of Ca²⁺ currents and bursting mechanisms for the pancreatic β-cell. J Membr Biol 1992;127:9–19.

[X32] Magnus G, Keizer J. Minimal model of β-cell mitochondrial Ca²⁺ handling. Am J Physiol 1997;273:C717–C733. （线粒体代谢模型。）

[X33] Bertram R, Pedersen MG, Luciani DS, Sherman A. A simplified model for mitochondrial ATP production. J Math Biol 2006;53:325–338.

[X34] Atwater I, Rosario L, Rojas E. Properties of the Ca-activated K⁺ channel in pancreatic β-cells. Cell Calcium 1983;4:561–570. （K(Ca) 通道——§9.1.4 supercell 假设的基础。）

[X35] Sneyd J, et al. A model of the calcium oscillations in exocrine pancreas. Cell Calcium 1991;12:151–158.

[X36] Kopell N, Ermentrout GB. Symmetry and phase-locking in chains of weakly coupled oscillators. Commun Pure Appl Math 1986;39:623–660. （弱耦合振子同步——§9.1.6 Plant 模型引用。）

[X37] Minorsky N. Nonlinear Oscillations. Van Nostrand, 1962. （Liénard 形式——Exercise 8 引用。）

[X38] Stoker JJ. Nonlinear Vibrations in Mechanical and Electrical Systems. Wiley, 1950. （van der Pol 振子——Exercise 8 引用。）

[X39] Tsaneva-Atanasova K, et al. Quantification of autocrine/paracrine effects in pituitary gonadotrophs using a compartmentalized mathematical model. J Cell Sci 2005;118:4485–4500.