第三章：膜离子通道读书笔记

书籍信息：Keener & Sneyd, Mathematical Physiology (I: Cellular Physiology), Springer 2009

章节摘要：本章系统介绍了膜离子通道的数学模型与生物物理机制，涵盖电流-电压关系、独立性与饱和性、电扩散模型、壁垒模型、通道门控、单通道分析以及反应速率理论等核心内容。

---

一、电流-电压关系

1.1 基本概念

离子电流的表达式与热力学第一性原理不同，无法从通用理论推导，而必须依赖于描述离子通道的特定模型。本章介绍了两种最常用的离子电流模型：线性模型和Goldman-Hodgkin-Katz（GHK）模型。

线性模型假设离子电流是膜电位的线性函数：

\[I_{Na} = g_{Na}(V - V_{Na}) \qquad(3.1)\]

其中 \(V_{Na} = \frac{RT}{F} \ln \frac{[Na^+]_e}{[Na^+]_i}\) 是钠离子的Nernst势，\(V = V_i - V_e\)。

GHK电流方程（恒定场模型）通过对Nernst-Planck方程在恒定电场假设下积分得到：

\[I_{Na} = P_{Na} \frac{F^2 V}{RT} \frac{[Na^+]_i - [Na^+]_e \exp\left(\frac{-VF}{RT}\right)}{1 - \exp\left(\frac{-VF}{RT}\right)} \qquad(3.2)\]

两种模型的关键区别在于：当存在两种或两种以上离子同时通过膜时，反转电位 \(V_r\) 的形式取决于所采用的模型。

GHK模型的反转电位：

\[V_r = \frac{RT}{F} \ln \frac{P_{Na}[Na^+]_e + P_K[K^+]_e}{P_{Na}[Na^+]_i + P_K[K^+]_i} \qquad(3.3)\]

线性模型的反转电位：

\[V_r = \frac{g_{Na}V_{Na} + g_KV_K}{g_{Na} + g_K} \qquad(3.4)\]

1.2 稳态与瞬时电流-电压曲线

离子通道可以响应膜电位的变化而开启或关闭。测量I-V曲线时需要区分两种不同的I-V曲线：

• 瞬时I-V曲线：电压快速变化后立即测量，反映单个开放通道的特性，假设开放通道的电流变化是瞬时的。
• 稳态I-V曲线：在通道有时间响应电压变化后测量，反映单个通道特性与开放通道比例的乘积。

1.3 离子流动的两类基本模型

描述开放通道中离子流动主要有两类模型：

1. 连续介质模型：将通道描述为连续介质，离子电流由Nernst-Planck电扩散方程决定，通过泊松方程与电场耦合。
2. 壁垒模型（ barrier models）：将通道建模为一系列由壁垒分隔的结合位点，离子通过"跳跃"越过壁垒从一 个结合位点移动到另一个。

---

二、独立性、饱和性与Ussing通量比

2.1 独立性原理

如果离子通过通道的流动与其他离子无关，则称该通道满足独立性原理。设外部浓度为 \(c_e\)，内部浓度为 \(c_i\)，则内向通量满足：

\[J_{in} = k_e c_e \qquad(3.5)\]

外向通量满足：

\[J_{out} = k_i c_i \qquad(3.6)\]

2.2 Ussing通量比

设 \(V_S\) 为离子S的Nernst势，V为跨膜电位差。Ussing通量比定义为内向通量与外向通量的比值：

\[\frac{J_{in}}{J_{out}} = \exp \frac{z(V_S - V)F}{RT} \qquad(3.10)\]

或等价形式：

\[\frac{J_{in}}{J_{out}} = \frac{c_e}{c_i} \exp \frac{-zVF}{RT} \qquad(3.11)\]

重要性质：满足独立性原理的通道必然满足Ussing通量比，但反之不一定成立。许多通道在大离子浓度下表现出对独立性的偏离，但仍满足Ussing通量比。

2.3 单离子孔道的饱和模型

单离子孔道模型假设通道只有一个结合位点，每个结合位点只能结合一个离子。当只有一个结合位点时，稳态通量：

\[J = \frac{k_0 k_1 c_0 - k_{-1} k_{-2} c_2}{k_0 c_0 + k_{-2} c_2 + k_{-1} + k_1} \qquad(3.72)\]

该模型满足Ussing通量比但不满足独立性原理，表现为浓度饱和效应。

2.4 多离子孔道与单文件传导

Hodgkin和Keynes（1955）提出的多离子通道机制假设：K⁺离子以单文件形式通过通道，通道内有n个结合位点。对于这种单文件机制，通量比为：

\[\frac{J_{in}}{J_{out}} = \frac{c_e}{c_i} \exp \frac{-VF \cdot n}{RT} \qquad(3.17)\]

实验数据显示K⁺通道的通量比指数约为2.5，表明至少2个K⁺离子同时通过该通道。

---

三、电扩散模型

3.1 Poisson-Nernst-Planck方程

电扩散模型基于Nernst-Planck方程描述离子在浓度梯度和电场作用下的移动。设两种离子S₁和S₂的浓度分别为 \(c_1\) 和 \(c_2\)，电位 \(\phi(x)\) 满足泊松方程：

\[\frac{d^2\phi}{dx^2} = -\frac{zqN_a}{\epsilon}(c_1 - c_2) \qquad(3.19)\]

稳态通量密度满足Nernst-Planck方程：

\[J_1 = -D_1 \frac{dc_1}{dx} + \frac{zF}{RT} c_1 \frac{d\phi}{dx} \qquad(3.20)\]

\[J_2 = -D_2 \frac{dc_2}{dx} - \frac{zF}{RT} c_2 \frac{d\phi}{dx} \qquad(3.21)\]

边界条件：\(c_1(0) = c_i, c_1(L) = c_e, c_2(0) = c_i, c_2(L) = c_e, \phi(0) = V, \phi(L) = 0\)

3.2 短通道或低浓度极限

当通道长度很短或离子浓度很低时（\(\lambda \ll 1\)），令 \(\lambda = 0\) 得到恒定电场假设，电位梯度为常数：

\[\frac{d\psi}{dy} = -v \qquad(3.28)\]

此时得到GHK电流方程：

\[I_1 = \frac{D_1 F^2}{LRT} V \frac{c_i - c_e \exp\left(\frac{-zVF}{RT}\right)}{1 - \exp\left(\frac{-zVF}{RT}\right)} \qquad(3.32)\]

3.3 长通道极限

当通道长度趋于无穷时（\(\eta = 1/\lambda\) 很小），电位满足 \(u_1 = u_2\)，浓度线性分布，电流为线性I-V曲线：

\[J_1 = \frac{D_1 zF}{LRT} \frac{ce - ci}{ln(ce/ci)} (V - \frac{RT}{zF} ln \frac{ce}{ci}) \qquad(3.42)\]

总结：通过取PNP方程的不同极限，可分别得到GHK I-V曲线（短通道/低浓度）或线性I-V曲线（长通道/高浓度）。GHK方程满足独立性和Ussing通量比，而线性I-V曲线两者都不满足。

---

四、壁垒模型

4.1 基本概念

壁垒模型假设离子通过通道的过程是跳过一系列自由能壁垒的过程。势能轮廓图的局部极小值对应结合位点，局部极大值对应阻碍离子流动的壁垒。

速率常数：

\[k_j = \kappa \exp \frac{-G_j}{RT} \qquad(3.43)\]

其中 \(\kappa\) 是具有时间倒数量纲的因子，\(G_j\) 是需要跨越的势垒高度。

4.2 非饱和壁垒模型

最简单的壁垒模型假设每个结合位点的离子浓度可以任意大（通道不饱和）。对于n个壁垒的系统，稳态通量：

\[J = k_0 \lambda \frac{c_0 - c_n e^{-v}}{e^{-n(g+v/n)} - 1} \cdot \frac{e^{-(g+v/n)} - 1}{e^{-(g+v/n)} - 1} \qquad(3.66)\]

均匀膜近似：设 \(g = 0\) 且 \(n \to \infty\)，保持 \(n\lambda = L\) 不变，可得GHK方程：

\[J = P_S \frac{c_0 - c_n e^{-v}}{1 - e^{-v}} \qquad(3.68)\]

4.3 饱和壁垒模型：单离子孔道

最简饱和模型：具有单个结合位点的单离子孔道，通过动力学机制：

\[X + S_i \underset{k_{-1}}{\overset{k_0}{\rightleftharpoons}} XS \underset{k_{-2}}{\overset{k_1}{\rightleftharpoons}} X + S_e \qquad(3.69)\]

稳态通量由公式（3.72）给出，表现出饱和特性：当 \(c_0\) 较小时，通量近似线性；当 \(c_0\) 增大时，通量趋向最大值 \(k_1\)。

Ussing通量比验证：单离子饱和模型满足Ussing通量比，但不满足独立性原理。

多结合位点：当通道有n个壁垒（n-1个结合位点）时，稳态通量的一般形式为：

\[J = \frac{k_0 c_0 - k_{-n} c_n \pi_{n-1}}{\phi_{n-1} + \beta k_0 c_0 + k_{-n} c_n(\alpha \phi_{n-1} - \beta \phi_{n-2})} \qquad(3.86)\]

4.4 多离子孔道

多离子孔道模型允许通道内同时存在多个离子。对于两个结合位点的最简多离子模型，有4个可能的通道状态。其分析更为复杂，通常需要数值方法求解。

4.5 电生泵与交换体

通道的电压依赖性可能来源于转运蛋白构象变化时电荷跨膜移动所涉及的能垒。模型将构象变化速率常数设为：

\[k_2 = \bar{k}_2 \exp \frac{FV}{RT}, \quad k_{-2} = \bar{k}_{-2} \exp \frac{-FV}{RT} \qquad(3.97,3.99)\]

---

五、通道门控

5.1 两种状态的K⁺通道模型

最简单的K⁺通道门控模型假设通道可以处于闭合（C）或开放（O）状态：

\[C \underset{\beta(V)}{\overset{\alpha(V)}{\rightleftharpoons}} O \qquad(3.102)\]

设 \(g\) 为开放通道的比例，则：

\[\frac{dg}{dt} = \alpha(V)(1-g) - \beta(V)g \qquad(3.103)\]

或写成标准形式：

\[\tau_g(V) \frac{dg}{dt} = g_\infty(V) - g \qquad(3.104)\]

其中 \(g_\infty(V) = \alpha/(\alpha + \beta)\)，\(\tau_g(V) = 1/(\alpha + \beta)\)。

由自由能论证，平衡常数为：

\[\frac{\beta}{\alpha} = k_0 \exp \frac{aFV}{RT} \qquad(3.106)\]

\[g_\infty(V) = \frac{1}{1 + k_0 \exp(\frac{aFV}{RT})} = \frac{1}{2} + \frac{1}{2} \tanh(b(V - V_0)) \qquad(3.107,3.108)\]

5.2 多亚基模型

当通道由k个相同的独立亚基组成时，每个亚基可独立开启或关闭。通道电导与 \(n^k\) 成正比，其中n满足：

\[\frac{dn}{dt} = \alpha(1-n) - \beta n \qquad(3.112)\]

这种多亚基门控模型是Hodgkin-Huxley兴奋性模型的基础。

5.3 钠通道

Na⁺通道既激活又失活，需要更复杂的模型。假设通道有一个h亚基和两个m亚基，每 个亚基可开或闭，通道仅当所有亚基都开放时才导电。

状态 \(S_{ij}\)（i个开放的m亚基，j个开放的h亚基）的动力学由六个微分方程描述，但系统具有稳定的流形不变集：

\[x_{ij} = \text{二项分布} \cdot \text{因子} \]

其中m和h分别满足：

\[\frac{dm}{dt} = \alpha(1-m) - \beta m \qquad(3.117)\]

\[\frac{dh}{dt} = \gamma(1-h) - \delta h \qquad(3.118)\]

5.4 激动剂控制的离子通道

del Castillo和Katz（1957）提出的模型将激动剂结合步骤与门控步骤明确分开：

\[A + T \underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}} AT \underset{\alpha}{\overset{\beta}{\rightleftharpoons}} AR \qquad(3.127)\]

该模型将亲和力（激动剂与通道的结合）和效能（结合后通道开放的倾向）分开，是理解激动剂控制通道的核心框架。

5.5 药物与毒素

许多药物通过阻断特定离子通道发挥作用。通道阻滞剂模型假设未结合通道P和被阻滞通道B之间存在平衡：

\[P + D \underset{k_-}{\overset{k_+}{\rightleftharpoons}} B \qquad(3.131)\]

稳态下：

\[\frac{P}{P_0} = \frac{K_d}{K_d + D} \qquad(3.133)\]

其中 \(K_d = k_-/k_+\) 是解离平衡常数。

---

六、单通道分析

6.1 基本概念

自1970年代末以来，膜片钳记录技术的发展使测量含单个或少量离子通道的膜片上的离子电流成为可能。

单通道记录的核心数学理论基于离散空间连续时间马尔可夫过程。

单通道数据可提供宏观测量无法获得的信息：开放时间分布、闭合时间分布、通道在失活前的开关次数等。

6.2 钠通道的单通道分析

对于三态Na⁺通道模型（C→O→I），闭合时间分布为 \(\alpha e^{-\alpha t}\)，开放时间分布为 \((\beta + \gamma) e^{-(\beta + \gamma)t}\)。

首次转移到开放状态的概率为 \(A = \alpha/(\alpha + \delta)\)。

通道在最终失活前开放k次的概率：

\[P[N=k] = (AB)^k \frac{1-AB}{1-B} \cdot \frac{1}{B} \qquad(3.136)\]

6.3 激动剂控制通道的单通道分析

对于激动剂控制的三态模型（3.127），由于结合和解离转换不可直接观测，形成"隐藏马尔可夫过程"。

开时间分布为 \(\alpha e^{-\alpha t}\)（均值 \(1/\alpha\)）。

闭时间分布是两种指数分布的和，由公式（3.140）给出。

6.4 与实验数据的比较

实验数据通常以开放和闭合时间的列表形式给出，显示在直方图中。

使用对数时间变换可以更清楚地区分多个指数成分：对于概率密度函数 \(\phi(t)\)，绘制 \(e^x \phi(e^x)\) 关于 \(x = \ln(t)\) 的图。

---

七、附录：反应速率

7.1 Boltzmann分布

Arrhenius方程指出反应速率 \(k\) 与 \(\exp(-B/T)\) 成正比，其中 \(B = G/R\) 或 \(G/k\)。

Boltzmann分布的推导基于概率论：对于n个粒子分布在k个能量状态 \(\{U_i\}\) 的系统，最可能的分布为：

\[n_i = \alpha e^{-\beta U_i} \qquad(3.157)\]

当应用于反应速率时，处于活化能垒Ur的粒子比例为 \(e^{-\beta\Delta U}\)，因此：

\[k \propto e^{-\beta \Delta U} \qquad(3.159)\]

7.2 Fokker-Planck方程方法

对于单个分子行为，使用 Langevin 方程描述：

\[m\frac{d^2x}{dt^2} + \nu \frac{dx}{dt} + U'(x) = \sqrt{2\nu kT} W(t) \qquad(3.161)\]

忽略惯性项，得到Fokker-Planck方程：

\[\nu \frac{\partial p}{\partial t} = \frac{\partial}{\partial x}(U'(x)p) + kT \frac{\partial^2 p}{\partial x^2} \qquad(3.163)\]

稳态解为连续形式的Boltzmann分布：

\[p(x) = \frac{1}{A} \exp \frac{-U(x)}{kT} \qquad(3.164)\]

7.3 Kramers结果

反应速率可由Fokker-Planck方程的平均首次退出时间模型得出。对于势能为 \(U(x) = Gu(Lx)\) 的系统，平均首次通过时间为：

\[\tau(0) \approx \frac{\pi \nu L^2}{\sqrt{G u''(0)|u''(1)|}} e^{G/kT} \qquad(3.179)\]

由此得到Arrhenius速率表达式，其中：

\[\kappa = \frac{\sqrt{G u''(0)|u''(1)|}}{\pi \nu L^2} \qquad(3.180)\]

这就是著名的Kramers公式，表明前因子 \(\kappa\) 与温度无关（当势垒较高时），但与势垒高度G有关。

---

公式汇总表

编号	公式名称	公式形式	适用场景
(3.1)	线性I-V关系	\(I_{Na} = g_{Na}(V - V_{Na})\)	枪乌贼轴突Na⁺通道
(3.2)	GHK电流方程	\(I_{Na} = P_{Na} \frac{F^2 V}{RT} \frac{[Na^+]_i - [Na^+]_e \exp(-VF/RT)}{1 - \exp(-VF/RT)}\)	脊椎动物轴突
(3.10)	Ussing通量比	\(\frac{J_{in}}{J_{out}} = \exp\frac{z(V_S - V)F}{RT}\)	独立性原理验证
(3.17)	单文件通量比	\(\frac{J_{in}}{J_{out}} = \frac{c_e}{c_i} \exp\frac{-VF \cdot n}{RT}\)	多离子通道
(3.66)	非饱和壁垒模型	\(J = k_0\lambda \frac{c_0 - c_n e^{-v}}{e^{-n(g+v/n)} - 1}\)	无饱和离子通道
(3.72)	单离子饱和模型	\(J = \frac{k_0 k_1 c_0 - k_{-1} k_{-2} c_2}{k_0 c_0 + k_{-2} c_2 + k_{-1} + k_1}\)	饱和离子通道
(3.86)	多结合位点通量	\(J = \frac{k_0 c_0 - k_{-n} c_n \pi_{n-1}}{\phi_{n-1} + \beta k_0 c_0 + k_{-n} c_n(\alpha\phi_{n-1} - \beta\phi_{n-2})}\)	多结合位点通道
(3.103)	两态门控动力学	\(\frac{dg}{dt} = \alpha(V)(1-g) - \beta(V)g\)	电压门控通道
(3.107)	稳态开放概率	\(g_\infty(V) = \frac{1}{1 + k_0 \exp(aFV/RT)}\)	门控电荷移动
(3.112)	多亚基开放动力学	\(\frac{dn}{dt} = \alpha(1-n) - \beta n\)	多亚基通道
(3.130)	激动剂稳态开放概率	\(\bar{\gamma}_1 = \frac{1}{1 + \frac{\beta}{\alpha} + \frac{\alpha}{k_{-1}}a}\)	激动剂控制通道
(3.148)	Arrhenius方程	\(k_i = \kappa \exp\frac{-G_i}{RT}\)	反应速率理论
(3.179)	Kramers公式	$\tau(0) \approx \frac{\pi \nu L^2}{\sqrt{G u''(0)	u''(1)

---

章节要点总结

1. 电流-电压关系的模型依赖性：Nernst方程是通用的，但离子电流表达式不是通用的。线性模型和GHK模型是最常用的两种近似。

2. 独立性原理与Ussing通量比：独立性原理意味着通量与浓度成正比；Ussing通量比是更基本的约束，可用于验证通道机制。

3. PNP方程的两种极限：短通道/低浓度→GHK方程；长通道/高浓度→线性I-V曲线。

4. 壁垒模型的核心思想：离子通过跳跃越过势能壁垒，速率常数由Arrhenius方程给出。

5. 单离子vs多离子通道：单离子通道满足Ussing通量比但不满足独立性；多离子通道可能违反Ussing通量比。

6. 通道门控的马尔可夫性质：通道在离散状态之间随机转换，可用马尔可夫过程描述。

7. 单通道数据的丰富信息：单通道记录包含宏观测量无法获得的动力学信息。

8. Kramers反应速率理论：基于Fokker-Planck方程的分子穿越势垒时间模型，是Arrhenius方程的统计力学基础。

---

参考文献：本章主要参考文献包括Hodgkin & Huxley (1952a), Hodgkin & Keynes (1955), Ussing (1949), Läuger (1973), Hille (2001), Colquhoun & Hawkes (1977, 1981, 1982), Kramers (1940)等经典著作。