第七章：钙动力学读书笔记
书籍信息
书名：Mathematical Physiology (I: Cellular Physiology) 作者：James Keener & John Sneyd 出版社：Springer, 2009（第二版） 章节：第七章——Calcium Dynamics（钙动力学）
一、概述与生理背景
钙离子（Ca²⁺）是细胞内最重要的第二信使之一，其浓度变化展现出极其复杂的时空行为模式。许多类型的细胞在受到激动剂刺激后会表现出Ca²⁺振荡现象。根据振荡产生机制的不同，可将Ca²⁺振荡分为两大类：第一类依赖于细胞膜电位的周期性波动以及经由电压门控Ca²⁺通道进入细胞的Ca²⁺（例如心脏细胞）；第二类发生在电压钳制条件下。本章主要关注第二类振荡，并进一步根据Ca²⁺通量是通过雷诺定受体（ryanodine receptors, RyR）还是IP₃受体（IP₃ receptors, IPR）进行分类讨论。
Ca²⁺振荡在众多细胞控制过程中发挥着重要作用。Berridge等人（2003）的综述指出，Ca²⁺振荡涉及卵子受精时的激活、轴突生长、细胞迁移、基因表达、植物根毛节结的形成、肌肉发育以及上皮细胞释放细胞因子等多种生理过程。在许多情况下，振荡作为一种频率编码信号，使细胞能够利用Ca²⁺作为第二信使，同时避免长时间高Ca²⁺浓度带来的毒性效应。例如，已知促性腺激素细胞中的胞吐作用依赖于Ca²⁺振荡的频率，基因表达也可以被Ca²⁺峰频率所调控。然而，仍有许多例子中Ca²⁺振荡所携带的信号尚未被明确解码。
IP₃依赖性Ca²⁺振荡的周期从几秒到几分钟不等。在某些细胞类型中，这些振荡发生在IP₃浓度恒定的情况下；而在另一些细胞类型中，它们似乎由[IP₃]的振荡所驱动。其他因素（如磷酸酶和激酶）对IP₃受体的调节也在设定振荡周期中起重要作用，而细胞外Ca²⁺的内流是另一个重要的调节机制。因此，控制Ca²⁺振荡的机制多种多样，期望一个单一模型能够捕捉所有细胞类型的重要行为是不现实的。然而，大多数模型有很多共同之处，通过研究少量模型可以学到很多关于整体方法的知识。Berridge等人（2003）提出了"细胞Ca²⁺信号工具包"的概念，该工具包包含大量组件（受体、G蛋白、通道、缓冲剂、泵、交换器等）；通过表达所需的组件，每个细胞可以微调其细胞内Ca²⁺的时空特性。Ca²⁺动力学模型遵循类似的模块化方法，每个细胞工具包成员都有相应的建模工具包模块。
图7.2展示了来自多种细胞类型的典型Ca²⁺振荡，包括：肝细胞受血管加压素刺激后的振荡、大鼠腮腺受卡巴胆碱刺激后的振荡、促性腺激素细胞受促性腺激素释放激素刺激后的振荡、金色仓鼠卵受精后的振荡，以及两种不同浓度卡巴胆碱刺激胰岛素瘤细胞后的振荡。
二、均匀混合细胞模型：钙振荡
2.1 基本模型框架
假设细胞是均匀混合的，则每种物质的浓度在整个细胞中均匀分布。用\(c\)表示细胞质中游离Ca²⁺离子的浓度，\(c_e\)表示内质网（ER）中Ca²⁺的均匀浓度。Ca²⁺的守恒方程为：
其中\(v\)是区室的体积，\(\tilde{J}_{net}\)是进入该区室的净Ca²⁺通量（单位：摩尔每秒）。通常假设细胞及其内部区室的体积是恒定的，因此守恒方程写为：
其中\(J_{net}\)是进入细胞质的净通量（每细胞质体积），单位是浓度每秒。
具体而言，根据图7.1中的通量，有：
对于ER中Ca²⁺的方程，必须考虑细胞质体积与ER体积的不同：
其中\(\gamma = v_{cyt}/v_{ER}\)是细胞质体积与ER体积的比率。
2.2 钙内流
Ca²⁺从细胞外进入细胞的内流通常依赖于电压。然而，当Ca²⁺振荡发生在恒定电压下时（非兴奋性细胞的典型情况），电压依赖性不重要。内流还依赖于许多其他因素，包括Ca²⁺、IP₃、花生四烯酸和ER中的Ca²⁺浓度。有证据表明，在某些细胞类型中，ER的耗竭会导致经由储存操纵通道（store-operated channels, SOCs）的Ca²⁺内流增加。一个常见的简化假设是假设\(J_{leak}\)是IP₃浓度\(p\)的线性递增函数：
2.3 线粒体
线粒体在Ca²⁺处理方面是一个高度复杂的过程。线粒体的功能之一是大量摄取和释放Ca²⁺，但速度相对较慢。因此，线粒体倾向于调节波的后缘、降低波幅，并改变长期的振荡行为。
2.4 钙泵和交换器
钙ATP酶：早期Ca²⁺ATPase泵的模型使用Hill方程类型。例如，来自Lytton等人（1992）的数据显示，通过ATPase的通量大约是\(c\)的S形函数，Hill系数约为2。因此，一个常见的模型是：
MacLennan等人（1997）构建了一个更详细的模型，该泵可以处于两种基本构象之一：E₁和E₂。在E₁构象中，泵从细胞质中结合两个Ca²⁺离子，然后暴露一个磷酸化位点。一旦磷酸化，泵切换到E₂构象，此时Ca²⁺结合位点暴露于ER管腔，且对Ca²⁺的亲和力低得多。由此产生的稳态通量为：
钙交换器：另一种重要的Ca²⁺从细胞质中移除的方式是通过Na⁺-Ca²⁺交换器，它们以牺牲三个Na⁺进入为代价从细胞质中移除一个Ca²⁺。
三、IP₃受体模型
3.1 受体的基本特性
IP₃受体的基本特性是它们对Ca²⁺或IP₃的阶跃变化表现出时间依赖性响应。因此，对IP₃或Ca²⁺的阶跃增加，受体开放概率首先上升到峰值，然后下降到较低的平台。这种下降称为受体的适应，因为它是针对维持的Ca²⁺或IP₃浓度的适应。
适应现在被认为部分是由于Ca²⁺不仅刺激自身释放，而且在较慢的时间尺度上抑制自身释放这一事实。假设Ca²⁺对IP₃受体的连续激活和失活是IP₃依赖性Ca²⁺振荡和波的基础。
3.2 八态IP₃受体模型
De Young和Keizer（1992）的模型假设IP₃受体由三个等价且独立的亚基组成，所有亚基必须处于导电状态才能有Ca²⁺流动。每个亚基有一个IP₃结合位点、一个激活Ca²⁺结合位点和一个失活Ca²⁺结合位点，每个位点可以处于占据或未占据状态，因此每个亚基可以处于八种状态之一。微分方程为：
其中：
IP₃受体的稳态开放概率作为Ca²⁺的函数呈现钟形——在低Ca²⁺时，增加Ca²⁺会增加受体的开放概率；而在高Ca²⁺时，增加Ca²⁺会降低开放概率。
3.3 八态模型的简化
八态受体模型的复杂性为寻找保留其基本特性的更简单模型提供了充分动机。由于IP₃快速结合到其结合位点，Ca²⁺快速结合到激活位点，我们可以放弃这些结合过程的瞬态细节，假设受体相对于IP₃结合和Ca²⁺激活处于准稳态。简化后的模型可以写成：
开放概率为：
3.4 具有饱和结合速率的模型
Sneyd和Dufour（2002）构建了一个基于Taylor及其同事定性模型的IP₃受体模型，该模型具有饱和结合速率。主要特点包括：
- IP₃受体可以在没有Ca²⁺的情况下被IP₃打开，但电导较低
- IP₃受体可以在没有IP₃的情况下被Ca²⁺失活
- 一旦IP₃结合，IP₃受体可以自发失活（到关闭状态S）
- 一旦IP₃结合，IP₃受体也可以结合Ca²⁺来激活受体
- 一旦IP₃受体被Ca²⁺结合激活，可以通过结合额外的Ca²⁺被失活
- IP₃和Ca²⁺的结合是顺序的
四、简单钙动力学模型
4.1 封闭细胞模型
假设只有两个通量，IPR和SERCA通量：
由于唯一的通量在ER和细胞质之间，有：
因此\(c + \gamma c_e = c_t\)不变。这种模型称为封闭细胞模型。
使用Hill函数模拟SERCA泵：
4.2 更复杂的模型
如果忽略RyR和线粒体通量，使用IP₃受体的饱和结合模型，使用SERCA泵的四态Markov模型，使用质膜泵的Hill函数，并假设Ca²⁺缓冲是快速且线性的，则得到：
4.3 开放和封闭细胞模型
引入新变量：
这个新变量\(c_t\)是细胞中Ca²⁺的摩尔总数（细胞质和ER）除以细胞质体积，是细胞总Ca²⁺含量的度量。重新排列方程得到：
当\(\delta = 0\)时，Ca²⁺既不进入也不离开细胞，\(c_t\)保持恒定，这称为封闭细胞模型。当\(\delta \neq 0\)时，我们有开放细胞模型，其中细胞中Ca²⁺的总量可以随时间变化。
五、IP₃动力学
在上述模型中，Ca²⁺振荡发生在恒定IP₃浓度下。IP₃的作用是激活IPR；一旦受体被激活，Ca²⁺反馈接管，振荡周期由Ca²⁺对IPR反馈的动力学以及膜和ER通量的相互作用控制。
然而，这是一个过度简化。已知IP₃的产生速率依赖于Ca²⁺，在某些细胞类型中，Ca²⁺振荡伴随着[IP₃]的振荡。目前尚不清楚[IP₃]振荡是否是Ca²⁺振荡所必需的，因为前者可能只是被动的跟随后者。
Ca²⁺影响IP₃浓度的两种主要方式：首先，Ca²⁺激活PLC已知在许多细胞类型中发生；其次，IP₃以两种主要方式降解，每种都受到反馈调节——它被IP₃ 5-磷酸酶去磷酸化为肌醇1,4-二磷酸（IP₂），或被IP₃ 3-激酶磷酸化为肌醇1,3,4,5-四磷酸（IP₄）。
六、雷诺定受体
RyR是Ca²⁺从细胞内储存释放的第二种主要方式。RyR共享许多与IP₃受体相似的结构和功能特性。Ca²⁺可以通过Ca²⁺诱导的Ca²⁺释放（CICR）从肌浆网或内质网通过雷诺定受体释放。
牛蛙交感神经元中的Ca²⁺振荡：交感神经元对咖啡因或轻度去极化反应产生稳健且可重复的Ca²⁺振荡。尽管这些振荡依赖于外部Ca²⁺，但它们发生在固定的膜电位下，并涉及经由RyR从ER释放Ca²⁺。CICR的简单模型（Friel, 1995）为牛蛙交感神经元中这些振荡的行为提供了出色的定量描述：
七、钙波
7.1 空间分布模型
在某些细胞类型中，Ca²⁺振荡几乎均匀地发生在整个细胞中。在这种情况下，测量细胞中任何一点的Ca²⁺浓度都会得到相同的时间进程，均匀混合模型是合适的。然而，更常见的是，每次振荡表现为在细胞中移动的波——这些细胞内"振荡"实际上是周期性细胞内波。要模拟和理解这种空间分布行为，需要包含Ca²⁺扩散。
在细胞质和ER中，Ca²⁺假设通过正常扩散并与缓冲剂反应。因此：
均质化技术用于推导有效扩散系数的平均场方程。
7.2 螺旋波模拟
Xenopus卵母细胞中波传播的常见实验程序是在细胞内光释放一剂IP₃并观察随后的Ca²⁺活动。当在域中间释放一剂IP₃S₃时，它在释放部位引起大量Ca²⁺释放。然后IP₃S₃扩散穿过细胞，过程中释放Ca²⁺。释放的Ca²⁺对IP₃受体的激活可以导致储存的周期性Ca²⁺释放，Ca²⁺在IP₃受体之间的扩散有助于稳定波，产生规则的周期性行波。
7.3 行波方程和分岔分析
使用行波变量\(\xi = x + st\)，其中\(s\)是波速，可以将反应扩散方程转化为两个常微分方程：
行波、 行波前缘和周期波分别对应于行波方程的同宿轨道、异宿轨道和极限周期。
八、钙缓冲
8.1 基本缓冲模型
Ca²⁺在所有细胞中都被大量缓冲，至少99%的可用Ca²⁺与大的Ca²⁺结合蛋白结合。缓冲的基本化学反应可以表示为：
简单的Ca²⁺缓冲模型为：
8.2 快速缓冲或过量缓冲
如果缓冲剂具有快速动力学，其对细胞内Ca²⁺动力学的影响可以简单地分析。如果\(k^-\)和\(k^+\)相对于Ca²⁺反应的时间常数较大，我们取准稳态：
因此：
其中\(K = k^-/k^+\)。
有效扩散系数为：
8.3 缓冲波的存在性
缓冲波的存在条件为：
immobile缓冲剂（即\(Db = 0\)）对双稳态方程中行波的存在没有影响，而mobile缓冲剂在足够量存在时可以消除行波。
九、离散钙源
9.1 离散释放模型
在所有先前给出的模型中，Ca²⁺从ER的释放被假设为空间均匀的。实际上，情况并非如此。在Xenopus卵母细胞中，IPR排列成簇，密度约为每30μm²一个簇，每个簇包含约25个IPR。
假设Ca²⁺从离散释放位点释放但被连续移除：
9.2 火-弥散-火模型
在火-弥散-火（fire-diffuse-fire）模型中，一旦Ca²⁺达到释放位点的阈值\(c^*\)，该位点立即释放固定量\(\sigma\)的Ca²⁺。因此，Ca²⁺波通过释放位点的顺序 firing 来传播：
十、钙Puff与随机建模
10.1 随机IPR模型
在所有上述模型中，Ca²⁺释放被建模为确定性的。然而，现在已知这并不总是合适的。实际上，每个Ca²⁺振荡或波都是由称为puff的随机基本释放事件构建的，每个puff对应于单个或少量IPR的Ca²⁺释放。
在低IP₃浓度下，单个簇的点状释放发生；在高[IP₃]下，这些局部释放事件协调成全局细胞内波。
随机IPR模型的基本假设是每个IPR释放事件可以建模为随机Markov过程，而Ca²⁺的扩散和其他泵和Ca²⁺通量的作用可以确定性地建模。
10.2 随机波模型
Falcke（2003a,b）研究了随着[IP₃]增加Xenopus卵母细胞中从puff到波的转变。在低[IP₃]下，只能观察到puff；每个簇释放的Ca²⁺不足以刺激相邻簇的Ca²⁺释放。然而，随着[IP₃]增加，每个IPR的敏感性和释放的Ca²⁺量都增加。这允许从成核位点发展出全局波。
特别有趣的是，在模拟中，周期性全局波发生在IP₃浓度下，对于这些浓度，模型的确定性版本是非振荡的。换句话说，振荡波不一定是振荡动力学的结果。它们可以由随机过程产生。
十一、细胞间钙波
11.1 机械刺激的细胞间Ca²⁺波
Sanderson及其同事提出了一种上皮细胞中机械刺激细胞间Ca²⁺波的模型。他们提出机械刺激在受刺激细胞中产生大量IP₃，这种IP₃通过被动扩散移动，通过间隙连接从一个细胞到另一个细胞。由于IP₃从ER释放Ca²⁺，IP₃从细胞到细胞的扩散导致相应的细胞间Ca²⁺波。
模型中，上皮细胞培养被建模为方形细胞网格。假设IP₃通过被动扩散移动并以饱和动力学降解。因此，如果\(p\)表示[IP₃]：
11.2 部分再生
关于先前模型的主要问题之一是一个IP₃分子是否能够扩散通过多个细胞而不被降解。一个IP₃分子不太可能存活足够长的时间来引起传播距离达数百微米的细胞间波。然而，如果IP₃的产生是再生的，波可以潜在地无限传播。
Höfer等人（2002）表明，部分再生机制可以传播比简单IP₃扩散远得多的细胞间Ca²⁺波，但并非无限；波最终仍然终止。
十二、公式汇总表
| 公式编号 | 公式内容 | 描述 |
|---|---|---|
| (7.1) | \(\frac{d(vc)}{dt} = \tilde{J}_{net}\) | Ca²⁺守恒基本方程 |
| (7.2) | \(\frac{dc}{dt} = J_{net} = \tilde{J}_{net}/v\) | 浓度形式的守恒方程 |
| (7.3) | \(\frac{dc}{dt} = J_{IPR} + J_{RyR} + J_{in} - J_{pm} - J_{serca} - J_{on} + J_{off} + J_{uni} - J_{mito}\) | 细胞质Ca²⁺变化方程 |
| (7.5) | \(J_{in} = \alpha_1 + \alpha_2 p\) | 内流模型（IP₃线性函数） |
| (7.6) | \(J_{serca} = \frac{V_p c^2}{K_p^2 + c^2}\) | SERCA泵Hill函数模型 |
| (7.7) | \(J_{serca} = \frac{c^2 - K_1 K_2 K_3 K_4 K_5 K_6 c_e^2}{\alpha_1 c^2 + \alpha_2 c_e^2 + \alpha_3 c^2 c_e^2 + \alpha_4}\) | MacLennan SERCA模型稳态通量 |
| (7.16) | \(\frac{dy}{dt} = \frac{(k_{-4} K_1 K_2 + k_{-2} pK_4)c}{K_4 K_2(p + K_1)}(1 - y) - \frac{k_{-2} p + k_{-4} K_3}{p + K_3}y\) | IP₃受体简化模型方程 |
| (7.17) | \(\tau_y(c,p) \frac{dy}{dt} = y_\infty(c,p) - y\) | IP₃受体准稳态形式 |
| (7.19) | \(x_{110} = \frac{pc}{(p+K_1)(c+K_5)} h\) | 简化IP₃受体开放概率 |
| (7.20) | \(\tau_h(c,p) \frac{dh}{dt} = h_\infty(c,p) - h\) | IP₃受体失活变量方程 |
| (7.32) | \(\frac{dc}{dt} = (k_f P_O + J_{er})(c_e - c) - J_{serca}\) | 封闭细胞Ca²⁺模型 |
| (7.34) | \(J_{serca} = \frac{V_p c^2}{K_p^2 + c^2}\) | Hill函数SERCA模型 |
| (7.36) | \(\frac{dc}{dt} = k_f P_O + J_{er}(c_e - c) - J_{serca} + J_{in} - J_{pm}\) | 开放细胞Ca²⁺模型 |
| (7.43) | \(c_t = c + \frac{c_e}{\gamma}\) | 总Ca²⁺含量定义 |
| (7.44) | \(\frac{dc_t}{dt} = \delta(J_{in} - J_{pm})\) | 总Ca²⁺变化方程 |
| (7.46) | \(\frac{dc}{dt} = J_{L1} - J_{P1} + J_{L2} - J_{P2}\) | CICR模型细胞质方程 |
| (7.47) | \(\frac{dc_e}{dt} = \gamma(-J_{L2} + J_{P2})\) | CICR模型ER方程 |
| (7.53) | \(\frac{\partial c}{\partial t} = \nabla \cdot (D_c \nabla c) - J_{on} + J_{off}\) | 空间分布细胞质方程 |
| (7.70) | \(\frac{\partial c}{\partial t} = D_c \nabla^2 c + f(c) + k^- b - k^+ c(b_t - b)\) | 缓冲Ca²⁺动力学 |
| (7.73) | \(b = \frac{b_t c}{K + c}\) | 快速缓冲准稳态 |
| (7.78) | \(D_{eff} = \frac{D_c + D_b \theta(c)}{1 + \theta(c)}\) | 有效扩散系数 |
| (7.93) | \(\frac{\partial c}{\partial t} = D_c \frac{\partial^2 c}{\partial x^2} - k_s c + \sum_n \delta(x - nL) f(c)\) | 离散释放模型 |
| (7.116) | \(\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.125) | \(\frac{\partial p}{\partial t} = D_p \nabla^2 p - \frac{V_p p}{k_p + p}\) | IP₃扩散降解模型 |
十三、关键概念总结
-
Ca²⁺信号工具包：Berridge等人提出的概念，包含多种Ca²⁺信号组分，每个细胞选择表达特定组件以微调其Ca²⁺时空特性。
-
适应现象：IP₃受体对持续刺激的反应先上升后下降的特性，由Ca²⁺的激活和失活时间尺度差异造成。
-
CICR机制：Ca²⁺诱导的Ca²⁺释放，是RyR和IPR工作的基本机制。
-
均匀混合 vs 空间分布模型：当Ca²⁺波存在时需要空间分布模型和扩散项。
-
火-弥散-火模型：描述离散释放位点通过Ca²⁺扩散通讯的简化模型。
-
随机性：Ca²⁺动力学本质上具有随机性，puff是构建全局事件的基本单元。
-
细胞间波传播：通过间隙连接扩散IP₃或Ca²⁺，或通过细胞外 messenger（如ATP）。
十四、参考文献
- Berridge MJ, Galione A (1988) Cytosolic calcium oscillators. FASEB J 2:3074-3082
- Berridge MJ, Bootman MD, Roderick HL (2003) Calcium signalling: dynamics, homeostasis and remodelling. Nat Rev Mol Cell Biol 4:517-529
- De Young GW, Keizer J (1992) A single-pool inositol 1,4,5-trisphosphate-receptor-based model for agonist-stimulated oscillations in Ca²⁺ concentration. Proc Natl Acad Sci USA 89:9895-9899
- Falcke M (2004) Reading the patterns in living cells – the physics of Ca²⁺ signaling. Adv Phys 53:255-440
- Sneyd J, Falcke M (2005) Models of the inositol trisphosphate receptor. Prog Biophys Mol Biol 89:173-191
- Sneyd J, Dufour JF (2002) A dynamic model of the type 2 inositol trisphosphate receptor. Proc Natl Acad Sci USA 99:2398-2403
- Friel DD (1995) [Ca²⁺]i oscillations in sympathetic neurons: an experimental test of a theoretical model. Biophys J 68:1752-1766
读书笔记完成日期：2026年5月11日
备注：本笔记涵盖了Keener & Sneyd《数学生理学》第七章的全部主要内容，包括Ca²⁺振荡和波的生理背景、均匀混合细胞模型、IP₃受体模型、Ca²⁺波的空间分布建模、Ca²⁺缓冲、离散钙源、火-弥散-火模型、随机建模以及细胞间Ca²⁺波等内容。所有LaTeX公式均以\(...\)格式呈现。