第九章 神经内分泌细胞（Neuroendocrine Cells）
一、概述与生理背景
神经内分泌细胞代表了生命系统中一类独特而重要的细胞类型，它们巧妙地结合了神经元的信息传递特性与内分泌细胞的激素分泌功能。这类细胞能够产生动作电位——这是神经元的标志性特征——同时又能合成和分泌激素，这是内分泌细胞的典型能力。细胞膜上分布着多种离子通道，这些通道的开关状态决定了细胞膜电位的变化。当膜电位发生特定模式的振荡时，钙离子通过电压门控通道进入细胞质，触发激素分泌颗粒的胞吐作用，从而将生物信号传递给下游靶细胞或器官。
在众多神经内分泌细胞中，胰腺β细胞和垂体促性腺激素细胞是研究最为深入的两个范例。胰腺β细胞位于胰岛中，其核心功能是响应血糖浓度变化调节胰岛素分泌。胰岛素是机体维持血糖稳态的关键激素，它促进肝脏、肌肉和脂肪组织对葡萄糖的摄取和利用。当血糖升高时，β细胞增加胰岛素分泌；当血糖下降时，胰岛素分泌停止。这种精密的负反馈调节机制对于维持机体能量代谢平衡至关重要。一旦该调控系统出现紊乱——通常表现为胰岛素分泌不足或胰岛素抵抗——就会导致糖尿病。糖尿病是一种全球范围内广泛流行的慢性代谢性疾病，若控制不当，可引发肾脏衰竭、心血管疾病、视网膜病变等多种严重并发症，甚至危及生命。
垂体促性腺激素细胞则分布在下丘脑和垂体前叶，是下丘脑-垂体-性腺轴的重要组成部分。下丘脑神经元分泌的促性腺激素释放激素（GnRH）通过垂体门脉系统运输至垂体前叶，刺激促性腺激素细胞释放促性腺激素（LH和FSH），后者进一步调控性腺的发育和功能。
本章将详细介绍这两类神经内分泌细胞的数学模型，重点分析其电生理特性和钙离子动力学。贯穿全章的核心概念是快慢时间尺度分离与相互作用——爆发式放电（bursting）这种丰富的动力学行为，本质上源于具有不同时间尺度的变量子系统之间的耦合。
二、胰腺β细胞与Chay-Keizer模型
2.1 爆发式放电的发现与生理意义
胰腺β细胞表现出独特的爆发式放电电活动模式。不同于传统神经元那样产生单个动作电位，β细胞在适当条件下会产生一连串快速动作电位，随后是一段静息期，如此反复，形成"爆发-静息"交替的节律性模式。这种电活动模式与胰岛素分泌紧密相关：在爆发期，细胞产生密集的动作电位，钙离子大量涌入细胞质，高钙浓度触发胰岛素颗粒的胞吐释放；在静息期，钙离子被泵出细胞或储存在内质网中，胰岛素分泌减少或停止。
爆发式放电的周期跨度极大，从数秒到数分钟不等。根据周期长短，通常分为三类：快速爆发式放电（周期约2至5秒）、中等爆发式放电（周期约10至60秒）和慢速爆发式放电（周期约2至4分钟）。
2.2 Chay-Keizer模型的建立
Atwater等人基于大量实验数据识别出了对爆发式放电至关重要的细胞机制。随后，Chay和Keizer（1983）将这些生理知识转化为数学模型，该模型成为后续近二十年研究的基石。Rinzel（1985, 1987）进一步从非线性动力学角度分析了这一模型，揭示了爆发式放电背后的分岔机制。
Chay-Keizer模型包含三种主要离子电流：
（1）钙激活钾通道电流 \(I_{K,Ca}\)：电导受细胞内钙离子浓度调节
（2）电压门控钾通道电流 \(I_K\)：与Hodgkin-Huxley方程形式相同
（3）电压门控钙通道电流 \(I_{Ca}\)：用钠电流形式近似
综合以上电流及漏电流，膜电位方程为：
细胞内游离钙离子浓度方程为：
参数 \(f\) 是缓冲因子（通常约0.007），它使得钙离子方程的演化速度远慢于膜电位方程，从而建立了快慢时间尺度的分离。葡萄糖通过调节 \(k_c\) 来影响细胞行为——葡萄糖浓度升高会加快质膜钙ATPase泵的活性，增加 \(k_c\) 值。
2.3 相平面分析与分岔机制
Rinzel和Lee（1986）通过对Chay-Keizer模型进行简化，建立了爆发式放电的相平面几何框架。忽略m和h门控变量的动力学后，得到快子系统（\(V\)和\(n\)方程）和慢子系统（\(c\)方程）。由于 \(f\) 很小（~0.007），快慢分离自然成立。
快子系统的Z形分岔图：将c视为参数，分析快子系统 \((V, n)\) 的平衡态随c变化的分岔结构，得到特征性的"Z"字形曲线。当c很低时，系统有唯一的稳定平衡点（膜电位较高）；当c很高时，系统也有唯一稳定平衡点（膜电位较低）；当c处于中间范围时，系统存在三个平衡点——下支稳定结点、中支鞍点、上支不稳定结点。在上支不稳定结点周围可能存在稳定极限环。极限环在Hopf分岔点（\(c_{hb}\)）处产生，在同宿分岔点（\(c_{hc}\)）处消失。在 \(c_{hb} < c < c_{hc}\) 范围内，系统呈现双稳态——稳定固定点与稳定极限环共存。
c零cline：由 \(dc/dt = 0\) 得到 \(c = -k_1 I_{Ca}(V)/k_c\)，其位置与 \(k_c\) 成反比。
完整爆发周期：假设初始时 \(c > c_{hc}\)，系统位于Z形曲线下支的稳定平衡点。由于该点位于c零cline下方，c逐渐减小至鞍-结分岔点（SN），系统跳至上支，进入动作电位重复放电的活跃期。膜电位远高于c零cline使c逐渐增加，当c达到 \(c_{hc}\) 时极限环消失，系统返回下支，完成一个爆发周期。
葡萄糖浓度的调控作用：葡萄糖浓度升高使 \(k_c\) 增大，c零cline向左移动。这导致活跃期延长、静息期缩短；\(k_c\) 进一步增大时爆发变为持续放电；极高 \(k_c\) 时则呈完全持续活跃状态。
三、内质网钙离子与幻影爆发式放电模型
3.1 早期模型的困境
Chay-Keizer模型存在两个严重缺陷。第一，模型只能产生快速爆发式放电，无法解释中等和慢速爆发式放电的存在。第二，后来的实验研究表明，细胞内游离钙离子浓度的变化速度实际上相当快，与膜电位的变化几乎同步进行——同时测量膜电位和细胞质钙离子的实验显示，爆发时钙浓度的上升几乎与电压上升一样快。这说明钙离子本身并不是控制爆发周期的那个慢变量。
3.2 内质网的钙库功能与幻影爆发
Chay（1996a,b, 1997）提出，内质网钙离子浓度 \(c_e\) 可能才是真正的慢控制变量。内质网像一个低通滤波器，平滑地反映细胞质钙离子的平均水平，其变化速度远慢于细胞质钙离子。Bertram和Sherman（2004a）进一步发展了这一思想，建立了"幻影爆发式放电"（phantom bursting）模型。
膜电位子模型：
其中 \(I_{Ca} = g_{Ca} m_\infty(V)(V - V_{Ca})\)，\(I_K = g_K n(V - V_K)\)，\(I_{K,Ca} = g_{K,Ca} \omega(c)(V - V_K)\)，\(I_{K,ATP} = g_{K,ATP}(V - V_K)\)。钙激活钾通道的开放比例为：
钙离子子模型：细胞质和内质网钙离子的守恒方程为：
其中 \(J_{serca} = k_{serca} c\)、\(J_{pm} = k_{pm} c\)、\(J_{leak} = k_{leak}(c_e - c)\)。
3.3 快速与中等爆发式放电的机制
快速爆发式放电：当钙激活钾电导较高（\(g_{K,Ca} = 900\) pS）时，\(c_e\) 在整个爆发周期中变化极小，实际上可视为常数。
中等爆发式放电：降低 \(g_{K,Ca}\) 至700 pS后，\(c_e\) 的慢调制作用变得不可或缺。在活跃期，\(c_e\) 逐渐增加，推动c零cline向右移动；当零cline移出Z形曲线的稳定极限环区域时，系统进入静息期。在静息期，\(c_e\) 逐渐减小，零cline向左移动，直到鞍-结分岔再次触发新爆发。活跃期和静息期都被显著延长。
3.4 神经调节剂的作用
乙酰胆碱通过IP3途径增加细胞内钙释放。在模型中加入IP3受体的钙 flux后，当 \([IP3]\) 从0升至0.3 μM时，爆发模式从中等周期转变为快速周期，与实验观察一致。
四、糖酵解与复合爆发式放电
4.1 葡萄糖感知与ATP敏感性钾通道
在低葡萄糖条件下，ATP/ADP比值较低，ATP敏感性钾通道（\(I_{K,ATP}\)）被ADP激活而开放，使细胞超极化；在高葡萄糖条件下，ATP/ADP比值升高，通道被ATP抑制而关闭，细胞去极化，允许爆发式放电发生。
4.2 糖酵解振荡与复合爆发
Wierschem和Bertram（2004）提出，糖酵解途径中代谢物的慢振荡可能是产生慢速和复合爆发式放电的根源。糖酵解振荡子模型为：
其中 \(F([ATP], [ADP]) = [ATP](1 + k_{ADP}[ADP])^2\)。
ATP浓度的振荡会改变Z形分岔曲线的位置。当ATP较低时无爆发；较高时呈持续活跃；中间值时产生间歇性爆发。由于ATP本身在振荡，活跃期的长度呈周期性变化——时而延长、时而缩短，形成"爆发组成的爆发"即复合爆发的特征模式。
五、细胞簇中的同步爆发式放电
5.1 单细胞与胰岛的行为差异
孤立培养的单个β细胞通常不能产生规律、持久的爆发式放电。只有当数千个β细胞聚集形成胰岛并通过缝隙连接电耦合时，才能观察到每个细胞的有规律同步放电。缝隙连接由连接蛋白构成，允许离子直接在相邻细胞间流动。
5.2 通道共享假说
Atwater等人（1983）提出了"通道共享"机制。每个β细胞含有少量高电导的钙激活钾通道，单个通道的随机开放就能对膜电位产生显著扰动。当细胞形成电耦合簇后，每个通道的电流分散到整个网络，单个通道的影响变得微乎其微。每个细胞整合了大量通道的综合效应，电耦合越紧密，行为越规律。
5.3 超细胞模型与有限耦合
Sherman等人（1988）建立了"超细胞"模型，假设所有细胞通过零电阻缝隙连接电耦合。在超细胞模型中，膜电位方程与单细胞形式相同，但 \(p\) 表示簇中开放通道的比例：
当细胞数 \(N_{cell}\) 增加时，\(p\) 的标准差减小，行为越来越规律。
有限耦合模型的单个细胞方程为：
同步爆发式放电的出现需要较强的细胞耦合和较大的簇尺寸。存在一个使爆发周期最大化的最优耦合强度。
六、定性爆发式放电模型与分类
6.1 多项式简化模型
Hindmarsh和Rose（1982, 1984）以及Pernarowski（1994）将爆发式放电模型简化为多项式形式，保留了核心动力学机制。经变量替换后的标准化多项式快子系统为：
该系统存在三个平衡点：左侧稳定结点、中央鞍点（阈值）、右侧被稳定极限环包围的不稳定结点。鞍点的稳定流形起到阈值作用——足够大的扰动使轨迹趋近极限环，较小的扰动衰减回静息态。
6.2 慢变量z的调制
为产生爆发式放电，引入第三个慢变量 \(z\) 调制外加电流：
当 \(r = 0.001\) 和 \(s = 4\) 时，系统表现出与Chay-Keizer模型相同机制的爆发式放电。
6.3 其他细胞类型中的爆发式振荡
爆发式振荡广泛存在于多种细胞中，包括多巴胺能神经元、丘脑皮层中继神经元、下橄榄神经元、Aplysia R15神经元、垂体促性腺激素细胞等。
Aplysia R15神经元的Plant模型与β细胞模型有本质区别。β细胞模型依赖双稳态，爆发终止于同宿分岔，动作电位间期在爆发过程中单调增加；Plant模型无双稳态，爆发由两个慢变量的内在振荡驱动，在爆发开始和结束都穿过同宿分岔，瞬时周期呈抛物线特征。
6.4 爆发式振荡的分类
Rinzel（1987）最早提出分类框架： - I型（滞后型）：由双稳态和滞后环产生，如β细胞模型 - II型（慢振荡型）：由慢变量的内在振荡驱动，如Aplysia R-15神经元 - III型（亚临界Hopf型）：由亚临界Hopf分岔产生
Izhikevich（2000）大大扩展了这一分类，但其中许多类型尚未在真实生物细胞中找到对应实现。
七、下丘脑和垂体促性腺激素细胞
7.1 细胞与激素调控概述
下丘脑-垂体-性腺轴是调控生殖功能的核心内分泌系统。GnRH从下丘脑通过垂体门脉系统运输至垂体前叶，刺激促性腺激素细胞分泌LH和FSH。促性腺激素细胞中，钙离子信号以两种方式产生：自发或诱发的动作电位激活电压门控钙通道导致钙内流；以及GnRH刺激下通过IP3途径从内质网释放钙。
7.2 膜电位与钙离子模型
Li等人（1994, 1995, 1997）建立了促性腺激素细胞的综合模型。膜电位方程为：
钙激活钾通道的电导受膜下钙离子浓度 \(c_R\) 调节：
模型将细胞处理为球对称结构，钙离子浓度随径向距离变化。反应-扩散方程为：
细胞边界处为钙离子跨膜通量。
7.3 模型结果与机制
模型在不同IP3浓度下的响应与实验高度一致。未刺激时呈持续锋电位；添加激动剂后，锋电位频率急剧下降并呈现复杂爆发模式；最高浓度下爆发初期被抑制，随后逐渐恢复。
机制解释如下：IP3诱导的钙释放增加细胞质钙背景水平，激活钙激活钾通道，使动作电位更难产生，锋电位频率下降。持续刺激下质膜钙泵持续工作，内质网钙逐渐耗竭，对钾通道的激活减弱，锋电位频率逐渐恢复。模型展示了ER钙库、细胞质钙和膜电位之间的复杂三方相互作用。
八、公式汇总表与关键参数
8.1 重要公式汇总
| 编号 | 名称/描述 | 公式 |
|---|---|---|
| (9.1) | 钙激活钾通道电导 | \(g_{K,Ca} = \bar{g}_{K,Ca} \frac{c}{K_d + c}\) |
| (9.2) | 电压门控钾通道电导 | \(g_K = \bar{g}_K n^4\) |
| (9.3) | 电压门控钙通道电导 | \(g_{Ca} = \bar{g}_{Ca} m^3 h\) |
| (9.4) | Chay-Keizer膜电位方程 | \(C_m \frac{dV}{dt} = -(g_{K,Ca} + g_K)(V - V_K) - 2g_{Ca}(V - V_{Ca}) - g_L(V - V_L)\) |
| (9.5) | Chay-Keizer钙动力学 | \(\frac{dc}{dt} = f(-k_1 I_{Ca} - k_c c)\) |
| (9.6)-(9.8) | Rinzel-Lee简化快慢系统 | 快：\(V,n\)方程；慢：\(c\)方程 |
| (9.9) | Phantom bursting膜电位 | \(C_m \frac{dV}{dt} = -I_{Ca} - I_K - I_{K,Ca} - I_{K,ATP}\) |
| (9.17) | 钙激活钾通道开放比例 | \(\omega(c) = \frac{c^5}{c^5 + k_D^5}\) |
| (9.21) | 细胞质钙离子方程 | \(\frac{dc}{dt} = f_{cyt}(-\alpha I_{Ca} - J_{pm} + J_{leak} - J_{serca})\) |
| (9.22) | 内质网钙离子方程 | \(\frac{dc_e}{dt} = -\gamma f_{er}(J_{leak} - J_{serca})\) |
| (9.24)-(9.25) | 糖酵解振荡子 | \(\tau_c \frac{d[ATP]}{dt} = v_1 - F\)，\(\tau_c \frac{d[ADP]}{dt} = F - v_2 [ADP]\) |
| (9.26) | 糖酵解速率函数 | \(F([ATP], [ADP]) = [ATP](1 + k_{ADP}[ADP])^2\) |
| (9.29) | 超细胞K通道电导 | \(g_{K,Ca} = \bar{g}_{K,Ca} \frac{c}{K_d + c}\) |
| (9.35) | 超细胞膜电位方程 | \(C_m \frac{dV}{dt} = -I_{Ca} - I_K - \bar{g}_{K,Ca} p(V - V_K)\) |
| (9.44) | 多细胞有限耦合方程 | 含 \(g_c \sum_j d_{ij}(V_i - V_j)\) 耦合项 |
| (9.53)-(9.54) | 多项式快子系统 | \(\frac{dx}{dT} = y - x^3 + 3x^2 + I\)，\(\frac{dy}{dT} = 1 - 5x^2 - y\) |
| (9.55)-(9.57) | 完整多项式爆发模型 | 含慢变量 \(z\) 调制 |
| (9.67) | 促性腺激素细胞膜电位 | \(C \frac{dV}{dt} = -I_{Ca,T} - I_{Ca,L} - I_K - I_{K,Ca} - I_L\) |
| (9.68) | 促性腺激素细胞K通道 | \(I_{K,Ca} = g_{K,Ca} \frac{c_R^4}{c_R^4 + K_c^4} \phi_K(V)\) |
| (9.69)-(9.70) | 空间分布钙模型 | 反应-扩散方程 |
8.2 关键参数表
Chay-Keizer模型参数（表9.1）
| 参数 | 值 | 单位 |
|---|---|---|
| \(C_m\) | 1 | μF/cm² |
| \(\bar{g}_{K,Ca}\) | 0.02 | mS/cm² |
| \(\bar{g}_K\) | 3 | mS/cm² |
| \(\bar{g}_{Ca}\) | 3.2 | mS/cm² |
| \(\bar{g}_L\) | 0.012 | mS/cm² |
| \(V_K\) | -75 | mV |
| \(V_{Ca}\) | 100 | mV |
| \(V_L\) | -40 | mV |
| \(f\) | 0.007 | — |
| \(k_c\) | 0.02 | ms⁻¹ |
Phantom Bursting模型参数（表9.2）
| 参数 | 值 | 单位 |
|---|---|---|
| \(g_{Ca}\) | 1200 | pS |
| \(g_K\) | 3000 | pS |
| \(g_{K,Ca}\) | 300-900 | pS（可变） |
| \(g_{K,ATP}\) | 227.5-350 | pS（可变） |
| \(V_{Ca}\) | 25 | mV |
| \(V_K\) | -75 | mV |
| \(C_m\) | 5300 | fF |
| \(f_{cyt}\) | 0.01 | — |
| \(f_{er}\) | 0.01 | — |
| \(k_{pm}\) | 0.2 | ms⁻¹ |
| \(k_{serca}\) | 0.4 | ms⁻¹ |
| \(\gamma\) | 5 | — |
| \(k_{leak}\) | 0.0005 | ms⁻¹ |
复合爆发式放电模型参数（表9.3）
| 参数 | 值 | 单位 |
|---|---|---|
| \(g_{Ca}\) | 1200 | pS |
| \(g_K\) | 3000 | pS |
| \(g_{K,Ca}\) | 300 | pS |
| \(g_{K,ATP}\) | 350 | pS |
| \(f\) | 0.001 | — |
| \(v_1\) | 10 | mM |
| \(v_2\) | 185 | — |
| \(\tau_c\) | \(1.2 \times 10^{-6}\) | ms |
章节总结
本章系统介绍了神经内分泌细胞的数学建模方法，以胰腺β细胞和促性腺激素细胞为核心实例。
胰腺β细胞部分首先建立了Chay-Keizer模型，通过相平面分析揭示了爆发式放电产生的Z形分岔结构。幻影爆发式放电模型通过引入内质网钙离子作为真正的慢控制变量，成功解释了从快速到中等周期的多种爆发模式。糖酵解振荡与电活动的耦合解释了慢速复合爆发式放电。超细胞模型阐明了单个β细胞与胰岛中同步电活动的联系，说明通道共享和电耦合如何将不规则的噪声转化为规律的节律。
促性腺激素细胞模型展示了膜电位动力学与IP3介导的钙释放之间的复杂相互作用，再现了GnRH刺激下的多相电反应。
贯穿全章的核心概念是快慢时间尺度分离与相互作用。爆发式放电本质上源于具有不同时间尺度的变量子系统之间的耦合。双稳态（平衡态与极限环共存）提供了活跃期存在的底物，而慢变量对分岔参数的调制则驱动系统在双稳态区域内外往复运动。分岔理论（Hopf分岔、同宿分岔、鞍-结分岔）为理解这些复杂动力学提供了精确的几何框架。