第2章：细胞稳态

作者

本章作者与第一章相同——James P. Keener（犹他大学数学系）与 James Sneyd（奥克兰大学数学系）。本章是第 I 卷"细胞生理学"的核心工具章，角色是把第 1 章建立的"化学-热力学-酶动力学"工具链延伸到 跨膜物质转运 这一细胞生理学最普遍的现象。它把全书两个最大的动力学主题——扩散 + 反应 与 膜电位 + 离子转运——首次统一在"守恒律 + 本构关系"框架下。本章 2.6 节建立的 Nernst 方程 + GHK 电流方程 + 膜电容电路模型 三件套是后续第 3 章（膜离子通道）、第 5 章（兴奋性）、第 12 章（心脏电生理）的共同基础。2.9 节的随机过程附录则给第 3 章单通道分析、隐式"状态占据概率"框架做预备。

内容概述

本章围绕"细胞如何维持内环境稳定（细胞稳态）"这一核心问题，按"被动扩散 → 缓冲扩散 → 促进扩散 → 载体介导运输 → 主动运输 → 膜电位 → 渗透与体积控制 → 随机过程"的逻辑展开。

第一部分（§2.1–§2.2）建立扩散的基础。§2.1 用 Table 2.1 给出典型细胞内/外离子浓度（鱿鱼巨大轴突：胞内 [K+]=397 mM、[Na+]=50 mM，胞外 20 mM / 437 mM；红细胞：胞内 [K+]=136 mM，胞外 5 mM；Nernst 电位 VNa=+55–56 mV, VK=-77 至 -101 mV, VCl=-9 至 -86 mV；静息电位 -6 至 -99 mV）。§2.2.1 给出 Fick 定律 $\mathbf J = -D\nabla u$（2.5）和反应-扩散方程 $\partial u/\partial t = D\nabla^2 u + f$（2.7）。§2.2.2 给出 Einstein-Stokes 公式 $D = kT/(6\pi\mu a)$（2.8），并解释为何球形大分子 $DM^{1/3}$ 近似为常数（典型 $D$ 值范围：氢 4.5×10⁻⁵ cm²/s，肌红蛋白 5.1×10⁻⁷ cm²/s，烟草花叶病毒 5.3×10⁻⁸ cm²/s）。§2.2.3 给出膜的"欧姆定律" $J = (D/L)(c_l - c_r)$（2.14）。§2.2.4 给出毛细管扩散的误差函数解 (2.16) 和扩散时间 $t = x^2/D$（Table 2.3：10 nm → 100 ns，1 mm → 16.7 min，1 m → 31.7 年——这是生物学"长距离传输不可能靠纯扩散"的量化论据）。

第二部分（§2.2.5 + §2.3）讨论扩散-反应耦合。§2.2.5 推导出"缓冲扩散"的有效扩散系数 (2.31)：$D_{eff} = (D_h + D_b \cdot w_0/K_{eq}) / (1 + w_0/K_{eq})$。§2.3 用肌红蛋白-氧气的"促进扩散"作为核心应用：通过 $\epsilon_1, \epsilon_2$ 远小于 1 的准稳态化简，导出总流量 $J = (D_s/L)(1 + \mu\rho)(s_0 - s_L)$（2.40）——促进扩散的增强因子 $\mu\rho$ 可以达到 560 倍（用 Wittenberg 1966 的参数 $k_+ = 1.4 \times 10^{10}$ cm³/M/s, $k_- = 11$ s⁻¹, $L = 0.022$ cm, $e_0 = 1.2 \times 10^{-5}$ M/cm³）。§2.3.1 进一步把促进扩散用于肌肉呼吸：稳态方程 (2.43)–(2.44) 积分得 (2.47)，"刚好不出现氧债"的临界外氧浓度 (2.48) 比无促进扩散时显著降低（Fig. 2.4, 2.5）。

第三部分（§2.4）讨论载体介导运输。§2.4.1 葡萄糖转运体：四态 Markov 模型（图 2.6）→ 稳态流量 (2.54)，并演示快平衡近似 $K_d \ll 1$ 的简化（(2.56)–(2.59)）给出与极限相同的答案。§2.4.2 推广到同向（symport）和反向（antiport）转运，给出 (2.62) 和 (2.63)。关键不等式 (2.64)：当 $m$ 个 S 和 $n$ 个 T 反向转运时，$s_e/s_i \cdot m > t_e/t_i \cdot n$ 是 S 净流入的方向。§2.4.3 Na⁺-Ca²⁺ 交换器：3 Na⁺ vs 1 Ca²⁺ 反向转运，每周期转移一个净正电荷（electrogenic）——稳态流 (2.69)，结合电势 (2.75) 得 (2.78) $K_1 K_2 K_3 K_4 = \exp(FV/RT)$，代入 (2.79) 给出电压依赖的通量。重要的数值例子 (2.81)：当 $[Na^+]_e/[Na^+]_i = 8$ 且 $V = -85$ mV 时，$\exp(FV/RT) = 7 \times 10^{-5}$——电生性质让 Na⁺-Ca²⁺ 交换器在静息状态下排 Ca²⁺ 的能力提升 27 倍。

第四部分（§2.5）讨论主动运输。§2.5.1 简单 ATPase 模型（图 2.10）：六态循环，详细平衡给出 (2.85) $\prod K_i = [L_i]_{eq}/[L_e]_{eq} \cdot [ATP]_{eq}/([ADP]_{eq}[Pi]_{eq}) = \exp(\Delta G^0_{ATP}/RT) = 3.73 \times 10^{-6}$ M⁻¹（由 $\Delta G^0_{ATP} = -31$ kJ/mol 算出），代入 (2.86) 给出 (2.94) 的电压修正版本 (2.97)。§2.5.2 推广到带电离子。§2.5.3 Post-Albers 模型（图 2.11）的 Na⁺-K⁺ ATPase：8 状态循环，包括两个"封闭"（occluded）状态以防止 slippage。细胞能量学（§2.5 开头）：典型动物细胞 1/3 能量用于 Na⁺-K⁺ ATPase，神经细胞可达 2/3。§2.5.4 核质运输：Ran-GTP/GDP 系统，$\Delta G^0_{Ran-GTP binding} = -51$ kJ/mol，水解后下降为 -33 kJ/mol，RanGAP 加速 5×10⁵ 倍。

第五部分（§2.6）建立膜电位理论。§2.6.1 Nernst 方程 (2.104)：$V_S = (RT/zF) \ln([S]_e/[S]_i)$，由化学势 $\Delta G = RT \ln([S]_i/[S]_e) + zFV = 0$ 推出；强调 Nernst 方程只适用于单一离子且独立于离子如何穿越膜。§2.6.2 Gibbs-Donnan 平衡：当膜内有不可透的带电大分子时 ([X] fixed, charge $-z_x$)，(2.109) 给出 $[S]_i$ 的显式解 (2.110)；这是细胞内外电势差的来源之一。§2.6.3 Goldman-Hodgkin-Katz (GHK) 方程：从 Nernst-Planck 方程 (2.114) + 常场近似得单种离子电流 (2.123)；总电流为零得多离子 GHK 电位 (2.126) 或 Na⁺/K⁺/Cl⁻ 三离子特例 (2.127) $V_r = -(RT/F) \ln((P_{Na}[Na^+]_i + P_K[K^+]_i + P_{Cl}[Cl^-]_e)/(P_{Na}[Na^+]_e + P_K[K^+]_e + P_{Cl}[Cl^-]_i))$。重要（§2.6.3 末段）：GHK 方程不是普适的——它依赖常场假设。§2.6.4 电路模型：$C_m dV/dt + I_{ion} = 0$（2.130），$I_S = g_S(V - V_S)$（2.133）。

第六部分（§2.7）渗透与 van't Hoff 定律。从化学势出发推 $\pi_s v = nRT$（2.139）即 van't Hoff 公式（细胞生理浓度下 $\phi$ 修正：葡萄糖/乳糖 $\phi = 1.01$，NaCl $\phi = 0.93$，KCl $\phi = 0.92$）；渗透压可达 19,300 mmHg（37°C 1 osmole/L）——"红细胞在纯水中破裂就是因为这个原因"。

第七部分（§2.8）细胞体积控制。§2.8.1 泵-渗漏模型（pump-leak）：4 维 ODE（Na⁺, K⁺, Cl⁻, 体积）+ 1 个电中性方程 (2.159) → 5 个未知量；引入无量纲变量后稳态化为 (2.165) + (2.166)。关键定理 (2.168)：存在有限细胞体积的充要条件是 $3[Na^+]_e/g_{Na} > 2[K^+]_e/g_K$——即 Na⁺ 通量必须"足够大"，因为 Na⁺-K⁺ ATPase 每周期输出 3 Na⁺ 换 2 K⁺（净输出正电荷）才能维持渗透平衡。真实参数（鱿鱼：$g_K = 0.367$ mS/cm², $g_{Na} = 0.01$ mS/cm²）下此条件成立。§2.8.1 末段讨论两个简化：（1）电导视为常数（实际电压依赖，第 5 章处理）；（2）泵率视为常数。改进模型 (2.169) 让泵率 $P = \rho u^3$ 与内部 Na⁺ 浓度耦合，结果是膜电位对泵率不再敏感（图 2.17 vs 2.16）。§2.8.2 Na⁺ 转运上皮细胞的 Koefoed-Johnsen-Ussing 模型：黏液侧 Na⁺ 通道 + 浆膜侧 Na⁺-K⁺ ATPase；存在 (2.184) 的稳定体积条件 $N_m/N_s < (1+3\rho_n)/(1+2\rho_k)$，超出时细胞破裂。生物学意义：亨利环上皮细胞（Ch17）必须在极高黏液 Na⁺ 下工作，因此采用 §2.8.2 末段所述的"降低 Na⁺ 通道 + 增加 K⁺ 通道"自适应机制。

第八部分（§2.9 附录）随机过程。§2.9.1 Markov 过程：$P(x,t|x_1,t_1,x_2,t_2) = P(x,t|x_2,t_2)$（2.187）。§2.9.2 离散态 Markov 过程：主方程 $d\boldsymbol{\pi}/dt = A^T \boldsymbol{\pi}(t)$（2.199），其中 $A$ 是 (2.200) 形式的生成元矩阵；等待时间 $E(T_i) = 1/K_i$（2.207）。§2.9.3 Gillespie 算法：$T = -(1/K_i) \ln \xi$（2.213）——比朴素离散时间步法效率高得多。§2.9.4 扩散的随机游走解释。§2.9.5 Langevin 方程 $dx = \sqrt{2b(x,t)} dW$（2.225）。§2.9.6 Fokker-Planck 方程 (2.228) + 平均首次穿越时间 (2.236)。§2.9.7 是隐含亮点（Van Milligen 2006 的食物染色+明胶实验）：Fick 定律 $\partial_t c = \partial_x(D \partial_x c)$ 与 Fokker-Planck $\partial_t c = \partial_x^2(D c)$ 在 $D$ 不均匀时给出不同结果；实验支持 FP 形式，因此真实扩散应为 $J = -\nabla(Dc)$（2.243）而非 $J = -D\nabla c$。

主要数学结果： 1. 扩散方程 $\partial_t u = D\nabla^2 u + f$ + Einstein $D = kT/(6\pi\mu a)$（适用于稀薄分子溶液） 2. 促进扩散通量 (2.40)：$J = (D_s/L)(1 + \mu\rho)(s_0 - s_L)$ 3. Na⁺-Ca²⁺ 电生交换器流量 (2.79) 和 (2.81) 的 27 倍增强 4. Nernst 方程 (2.104) + GHK 方程 (2.123)、(2.126) 5. 膜电路方程 (2.130) + (2.133) 6. 细胞体积存在的充要条件 (2.168) 和 (2.184) 7. Gillespie 算法 (2.213) 8. 扩散的 Fokker-Planck 形式 (2.243) 而非 Fick 形式

前置知识：ODE、PDE 基础（特别是扩散方程的分离变量与误差函数解），基本的 Gibbs 自由能化学势（已在 Ch1 介绍），Markov 链基础（若不熟悉可读 Gardiner 2004 入门）。

核心方程与概念

1. 扩散与 Fick 定律

Fick 定律（2.5）：$\mathbf{J} = -D\nabla u$，其中 $D$ 是扩散系数（量纲 cm²/s）。注意：这是本构关系而非物理定律——它在浓度不太高时近似成立（§2.2.1 末段），与 Ohm 定律、Newton 冷却定律同级别。

Einstein 公式（2.8）：$D = kT/(6\pi\mu a)$，$k = R/N_A$。由此球形分子 $D \propto M^{-1/3}$（(2.11) + $\rho$ 近似常数），但小分子例外 $D \propto M^{-1/2}$（Table 2.2 中氢气 4.5×10⁻⁵ 远高于按 $M^{-1/3}$ 预测的 1.5×10⁻⁵）。

膜的欧姆定律（2.14）：$J = (D/L)(c_l - c_r)$，$L/D$ 是有效"电阻"。

扩散时间 $t = x^2/D$（§2.2.4）：Table 2.3 给出关键判断：扩散对 1 mm 量级传输已需要 16.7 min，对 1 m 需要 31.7 年——纯扩散无法解释长距离生理传输，因此循环系统、神经传导都需要主动或对流机制。

2. 缓冲扩散（buffered diffusion）

反应-扩散方程组（2.20）–(2.22)：$H^+ + B \rightleftharpoons HB$ 加上 $H^+$、$B$、$HB$ 的扩散。关键观察 $v + w = w_0$ 不变 → 守恒律 (2.23)。

有效扩散系数（2.31）：当 $u \ll K_{eq}$ 时，

\[ D_{eff} = \frac{D_h + D_b \cdot w_0/K_{eq}}{1 + w_0/K_{eq}} \]

$D_{eff}$ 是 $D_h$ 和 $D_b$ 的凸线性组合。当缓冲剂远多于溶质时（$w_0/K_{eq} \gg 1$）$D_{eff} \to D_b$（缓冲剂限速）；当缓冲剂稀少时 $D_{eff} \to D_h$（自由扩散限速）。

3. 促进扩散（肌红蛋白-O₂）

反应-扩散方程（2.32）–(2.34)：$\partial_t s = D_s s_{xx} - f$, $\partial_t c = D_c c_{xx} + f$（$f$ 是结合/解离速率）。

守恒律（2.36）：$D_s ds/dx + D_c dc/dx = -J$ 给出总通量 $J$ 的积分 (2.37)。

准稳态化简（2.38）+（2.39）：$c = e_0 s/(K + s)$（Michaelis-Menten 形式）代入得 $J = (D_s/L)(1 + \mu\rho)(s_0 - s_L)$（2.40），其中 $\rho = (D_c/D_s)(e_0/K)$, $\mu = K^2/((s_0+K)(s_L+K))$。$\mu\rho$ 是促进扩散的增强因子：低浓度极限 $\mu = 1$, $\rho = 560$（Wittenberg 1966 参数）。

生理意义（§2.3 末段）：自由 O₂ 通量 $J_0 = (D_s/L)(s_0 - s_L)$ 与促进扩散通量之差是"被肌红蛋白存储并缓慢释放"的部分。只有在 $s_L$ 足够低使肌红蛋白释氧时，绑定 O₂ 的通量才有效——这是为什么 $D_c \ll D_s$（560× 慢）但总通量仍可大得多。

§2.3.1 肌肉呼吸：（2.48）给出"刚好不出现氧债"的临界外氧浓度——有肌红蛋白时该临界浓度显著降低（Fig. 2.4, 2.5），即肌红蛋白扩展了肌肉可工作的有氧范围。

4. 载体介导运输（葡萄糖、Na/Ca 交换）

葡萄糖转运体（§2.4.1）：四态 Markov 模型（图 2.6）。稳态流 (2.54)

\[ J = \frac{C_0 (s_e - s_i)}{2K K_d^{-1} [(s_i + K + K_d)(s_e + K + K_d) - K_d^2]} \]

其中 $K = k_-/k_+$, $K_d = k/k_+$。$K_d$ 反映构象变化速率（类似扩散的"温度"）。注意（§2.4.1 中段）：有两种解释——(a) 群体模型（$p_i + p_e + c_i + c_e = C_0$，每个是浓度）；(b) 单分子模型（守恒 = 1，每个是概率）——两个给出相同的 ODE 但物理意义不同。Markov 解释（§2.4.1 末段注）：第 (b) 种解释下，每个变量是单分子处于该状态的概率，速率常数 $k_{ij}$ 是分子在状态 $i$ 的等待时间倒数。

快平衡简化（2.56）–(2.59)：$K_d \ll 1$ 时构象变化是限速步骤；可消去 $c_i, c_e, p_i, p_e$ 中的两个，得 $J = (kK/2)(s_i - s_e)/[(s_i + K)(s_e + K)]$——这正是 (2.54) 在 $K_d \to 0$ 的极限。

Symport/Antiport（§2.4.2）：$m$ 个 S 和 $n$ 个 T 同/反向转运。(2.62) symport 流量；(2.63) antiport 流量。重要不等式 (2.64)：antiport S 净流入需要 $(s_e/s_i)^m > (t_e/t_i)^n$。

Na⁺-Ca²⁺ 电生交换器（§2.4.3）：3 Na⁺ vs 1 Ca²⁺，每周期净转移 1 个正电荷出细胞。稳态流 (2.69) 和含电压版本 (2.79)

\[ J = \frac{k_1 k_2 k_3 k_4 (c_i n_e^3 - e^{FV/RT} c_e n_i^3)}{16 \text{ 个正项}} \]

关键 (2.81)：当 $n_e/n_i = 8$, $V = -85$ mV 时 $\exp(FV/RT) = 7 \times 10^{-5}$，因此电生性质让交换器在 $c_i/c_e > 7 \times 10^{-5}$ 时仍能净排 Ca²⁺（实际细胞内/外 Ca²⁺ 比为 10⁴ 量级）——即电生性质让交换器排 Ca²⁺ 能力提升 27 倍。

5. 主动运输（ATPase）

简单 ATPase（§2.5.1, 图 2.10）：六态循环。详细平衡约束 (2.85)

\[ \prod_{i=1}^{6} K_i = \frac{[L_i]_{eq}}{[L_e]_{eq}} \cdot \frac{[ATP]_{eq}}{[ADP]_{eq}[Pi]_{eq}} = e^{\Delta G^0_{ATP}/RT} = 3.73 \times 10^{-6} \text{ M}^{-1} \]

稳态流 (2.86) 表明通量正比于"远离平衡的距离"——细胞通过维持 [ATP] 远高于平衡值实现主动运输。重要物理结论（§2.5.1 中段）："如果 $[L_e]$ 足够高，可迫使泵反转，让 L 从外向内流同时合成 ATP"——这是 ATP 合酶（Ch7）的原理。

Na⁺-K⁺ ATPase（§2.5.3）：Post-Albers 模型 8 状态（图 2.11），含两个"封闭"（occluded）状态以防止 slippage（§2.5.3 中段）。稳态流 (2.100) 形式与 (2.86) 类似。能量学（§2.5 开头）：典型动物细胞 ~1/3 能量用于此泵，神经细胞 ~2/3。

Na⁺-K⁺ ATPase 与强心苷（§2.8.1 末段注）：digitalis（digoxin）抑制 30-40% Na⁺-K⁺ ATPase → 胞内 Na⁺ 升高 → Na⁺-Ca²⁺ 交换器受抑 → 胞内 Ca²⁺ 升高 → 心肌收缩力增强（治疗心衰）——这是 1980s 美国第四大处方药的药理机制。

6. 膜电位理论

Nernst 方程（2.104）

\[ V_S = \frac{RT}{zF} \ln \frac{[S]_e}{[S]_i} = \frac{kT}{zq} \ln \frac{[S]_e}{[S]_i} \]

关键注（§2.6.1 末段）：Nernst 方程只对单一离子成立；多离子情况下没有"普适"的反转电位，必须明确电流-电压模型。27°C 下 $RT/F = 25.8$ mV。

Nernst-Planck 方程（2.114）

\[ \mathbf{J} = -D \left( \nabla c + \frac{zF}{RT} c \nabla \phi \right) \]

GHK 单种离子电流（2.123）

\[ I_S = P_S \frac{z^2 F^2}{RT} V \frac{c_i - c_e \exp(-zFV/RT)}{1 - \exp(-zFV/RT)} \]

GHK 多离子零电流电位（2.126）

\[ V = -\frac{RT}{F} \ln \frac{\sum_{z=-1} P_j c_{j,e} + \sum_{z=1} P_j c_{j,i}}{\sum_{z=-1} P_j c_{j,i} + \sum_{z=1} P_j c_{j,e}} \]

Na⁺/K⁺/Cl⁻ 三离子特例（2.127）：$V_r = -(RT/F) \ln \frac{P_{Na}[Na^+]_i + P_K[K^+]_i + P_{Cl}[Cl^-]_e}{P_{Na}[Na^+]_e + P_K[K^+]_e + P_{Cl}[Cl^-]_i}$。重要（§2.6.3 末段）：GHK 方程依赖常场假设——其他模型给出不同的 $I_{ion}(V, c)$ 关系（Ch3 讨论）。

膜电路模型（2.130）+$I_S = g_S(V - V_S)$（2.133）

\[ C_m \frac{dV}{dt} + \sum_S I_S = 0 \quad \text{且} \quad I_S = g_S(V - V_S) \]

电容（2.129）：$C_m = k\epsilon_0/d$。典型细胞膜 $C_m = 1.0$ μF/cm²，由此反推介电常数 $k \approx 8.5$（介于真空 1 和油 3 之间）。

7. 渗透与 van't Hoff 定律

核心公式（2.139）$\pi_s v = nRT$（理想溶液）；修正 (2.141) $\pi_s v = \phi nRT$。37°C 1 osmole/L 对应 19,300 mmHg ≈ 250 m 水柱——"红细胞在纯水中破裂的原因"。

8. 细胞体积控制（pump-leak 模型）

4 维 ODE + 1 个代数约束（2.150）–(2.159）：变量 $[Na^+]_i, [K^+]_i, [Cl^-]_i, w, V$，约束是胞内电中性 (2.159)。

关键定理 (2.168)：存在有限稳态体积的充要条件是

\[ 3 \frac{[Na^+]_e}{g_{Na}} > 2 \frac{[K^+]_e}{g_K} \]

物理意义：Na⁺-K⁺ ATPase 每周期净输出 1 个正电荷，只有当 $g_{Na}$ 相对 $g_K$ 足够大时，漏入的 Na⁺ 才能被泵出，从而避免 Na⁺ 持续积累导致渗透压升高。

鳞虾轴突数值例子（§2.8.1 末段）：$g_K = 0.367$ mS/cm², $g_{Na} = 0.01$ mS/cm², $[Na^+]_e = 437$ mM, $[K^+]_e = 20$ mM → $3 \times 437/0.01 = 1.31 \times 10^5 > 2 \times 20/0.367 = 109$ ✓。

改进模型 (2.169)：$P = \rho u^3$（$u = [Na^+]_i/[Na^+]_e$）——泵率随内部 Na⁺ 浓度的 3 次方变化。结果（图 2.17 vs 2.16）：膜电位对泵率不敏感（自调节），细胞体积也更稳定。

KJU 上皮细胞模型（§2.8.2）：黏液侧 Na⁺ 通道 + 浆膜侧 Na⁺-K⁺ ATPase。稳定体积的条件 (2.184) $N_m/N_s < (1+3\rho_n)/(1+2\rho_k)$。细胞实际策略：降低黏液侧 Na⁺ 电导 + 增加体积激活的 K⁺ 电导（(2.180) 之后段落）——这让细胞能在远高于 (2.184) 阈值的黏液 Na⁺ 下工作。

9. 随机过程附录

Markov 主方程（2.199）$d\boldsymbol{\pi}/dt = A^T \boldsymbol{\pi}(t)$，$A$ 是 (2.200) 形式的生成元。注意：这是与 ODE 形式相同的方程——确定性 ODE 的"浓度"和随机 ODE 的"占据概率"是同一组方程。

等待时间（2.207）$E(T_i) = 1/K_i$。Gillespie 算法（2.213）$T = -(1/K_i) \ln \xi$。

Fokker-Planck 方程（2.228）$\partial_t p = -\partial_x(a p) + \partial_x^2(b p)$，对应 Langevin (2.225)。

平均首次穿越时间方程（2.236）$a(y) T_y + b(y) T_{yy} = -1$，边界条件：反射壁 $T_y = 0$，吸收壁 $T = 0$。例子（2.237）–(2.238)：$T(x) = (L^2 - x^2)/(2D)$。

Fick vs Fokker-Planck（§2.9.7, Van Milligen 2006 食物染色实验）：当 $D$ 不均匀时

Fick 形式：$\partial_t c = \partial_x(D c_x)$ → 稳态 $D c_x = 0$
Fokker-Planck 形式：$\partial_t c = \partial_x^2(D c)$ → 稳态 $D c = $ 常数

实验支持 Fokker-Planck：真实扩散应为 $J = -\nabla(D c)$ 而非 $J = -D\nabla c$——这是一个常被忽略的细节，对细胞内异质介质（如细胞骨架、膜蛋白拥挤区）有重要影响。

关键结论

扩散的生物学边界（§2.2.4 + Table 2.3）：$t = x^2/D$，1 mm 需要 16.7 min，1 m 需要 31.7 年。纯扩散无法解释长距离生理传输——这是循环系统（Ch11）、神经传导（Ch5）、轴突运输（Ch8）必须依赖主动机制的物理原因。
促进扩散的"反直觉"机制（§2.3）：肌红蛋白（$D_c \ll D_s$）反而能增强 O₂ 通量 560 倍，原因是总通量 $J$ 包含自由 O₂ 通量 + 绑定 O₂ 通量——绑定部分虽慢但量大，且只有当 $s_L$ 足够低使 Mb 释氧时生效。
电生性质让 Na⁺-Ca²⁺ 交换器提升 27 倍（§2.4.3 + (2.81)）：膜电位 $V = -85$ mV 下 $\exp(FV/RT) = 7 \times 10^{-5}$，这让交换器在 $c_i/c_e > 7 \times 10^{-5}$ 仍能净排 Ca²⁺——实际细胞 $c_i/c_e \sim 10^{-4}$ 是临界值。电生性是 Na⁺-Ca²⁺ 交换器在静息状态下能工作的根本原因。
细胞体积存在的充要条件 (2.168)：$3[Na^+]_e/g_{Na} > 2[K^+]_e/g_K$。这一条件不是显然成立的——$g_{Na}$ 远小于 $g_K$ 是细胞相对普遍的状态，关键是 $[Na^+]_e$ 远高于 $[K^+]_e$（10-20 倍）来补偿。当泵被 ouabain 抑制时 Na⁺-K⁺ ATPase 失效，$g_{Na}$ 渗漏持续，胞内 Na⁺ 持续上升 → 渗透压升高 → 细胞肿胀破裂（Exercise 25）。
GHK 方程不是普适的（§2.6.3 末段）：Nernst 方程来自热力学约束，是普适的；但 GHK 电流方程依赖常场假设，多离子反转电位依赖渗透性参数。对真实离子通道（特别是高电场、复杂几何）需要其他模型（Ch3）。
膜电容与典型参数（§2.6.4）：$C_m = 1$ μF/cm²，对应 $\epsilon \approx 8.5$，膜厚 ~7.5 nm。一个 $500$ μm 半径的鱿鱼轴突，$V = 100$ mV 时总电荷 $Q = \pi \times 10^{-8}$ C/cm，相比胞内 K⁺ 总电荷（$0.1\pi$ C/cm）只有 $10^{-7}$——细胞内电中性假设是极好的近似。
Post-Albers 模型的"封闭态"防 slippage（§2.5.3 中段）：Na⁺ 释放到胞外、K⁺ 释放到胞内时都经过封闭态（结合但不可接触任一侧），这保证每个循环只向前走、不后退。这是真实 ATPase 高效率的物理基础。
Fick 定律在非均匀 $D$ 下失效（§2.9.7 + Van Millegen 2006）：实验证明真实扩散是 $J = -\nabla(Dc)$（Fokker-Planck 形式）而非 $J = -D\nabla c$（Fick 形式）。这一区别对细胞内拥挤介质（细胞骨架、膜筏）的建模有重要影响——大多数教材仍用 Fick 形式，但 Keener-Sneyd 在此明确指出。
Gillespie 算法是离散态 Markov 模拟的最优方法（§2.9.3 + (2.213)）：$T = -\ln\xi/K_i$ 直接采样等待时间，比离散时间步法快得多、且无 $dt \to 0$ 的限制。对单离子通道记录的分析（Ch3）和生化网络模拟有核心重要性。

挑战和开放性问题

促进扩散的"低氧增强"是反直觉的（§2.3 末段）：肌红蛋白 $D_c \ll D_s$，单看 $D$ 应减慢通量；通量增强来自绑定 O₂ 的"大容量存储"。但什么条件下促进扩散会失效？当 $s_L$ 过高使 Mb 不释氧时，增强因子会消失——但 §2.3 末段没有给出 $s_L$ 临界值的显式条件。
Na⁺-Ca²⁺ 交换器模型的"细节耦合"（§2.4.3 末段）：原书假设 3 Na⁺ 和 1 Ca²⁺ 同时结合释放。真实（Hilgemann 2004）是逐步结合；Kang-Hilgemann 2004 已发展详细动力学模型，但本书模型与此的偏差未量化。
Post-Albers 模型的"电压依赖分布"（§2.5.3 末段）：Smith-Crampin 2004 把电压依赖放在 $K_2, K_5$（Na⁺ 结合/解离），但 §2.5.2 末段承认"模型可以有不同的 $\phi(V)$，只要 $\prod K_i$ 满足 (2.96)"——如何从实验上区分这些等价模型？原书没有给出。
GHK 方程的边界（§2.6.3 末段）：常场假设在 $V$ 大或 $D$ 不均匀时失效。真实离子通道（Ch3）有复杂几何（selectivity filter、gate）和高电场（~10⁷ V/m），需要 Poisson-Nernst-Planck 方程或其他更细模型。GHK 仍是好的零级近似，但定量预测差。
细胞体积控制模型的简化（§2.8.1 末段注）：(a) $g_{Na}, g_K$ 当作常数（实际电压依赖，Ch5）；(b) 泵率当常数（实际 ATP 依赖）。虽然 §2.8.1 末段给出了耦合 ATP 浓度的改进 (2.169)，但没有给出能自洽计算 ATP 浓度的闭合模型。Exercise 26 把电势依赖加入泵率，但分析也只到定性层面。
KJU 上皮模型的有限性（§2.8.2 末段）：基本模型在 (2.184) 阈值外细胞破裂，但亨利环（Ch17）必须在极高黏液 Na⁺ 下工作。真实细胞采用什么机制扩展工作范围？原书只说"降低 Na⁺ 电导 + 增加 K⁺ 电导"，没有给出分子机制（是哪些通道、什么反馈）。
§2.9.7 Fick vs Fokker-Planck 的实际后果（隐含问题）：作者指出了 Van Milligen 2006 的实验证据，但对细胞生理学的实际后果没有展开。细胞内（细胞骨架、膜蛋白、内质网等）都是强非均匀介质；用 Fick 而非 FP 形式会引入多大误差？这是值得深入研究的问题。
核质运输的数学模型（§2.5.4 末段）：原书只给出了定性的生物物理描述（Ran-GTP 亲和力差 10⁴ 倍，RanGAP 加速 5×10⁵ 倍），没有像 Na⁺-Ca²⁺ 那样给出完整的 Markov 模型和稳态流量。Exercise 14 留给读者。
随机 vs 确定性的边界（§2.9 开头 + 整个附录隐含）：离子通道、基因表达、钙火花等过程的"随机性"何时不可忽略？原书未给出统一的判定标准。这是当前系统生物学的重要开放问题。

个人反思与批判性分析

建模哲学层面：

Keener-Sneyd 在 Ch2 的"由守恒律出发 → 选择本构关系 → 求解具体问题"方法论是生物物理建模的标准范式。具体而言： 1. §2.2 给出"反应-扩散"框架 (2.6)，第 1 章给出"反应"框架，两者结合就是 Ch6 行波、Ch7 钙波的数学基础。 2. §2.4 的"载体介导运输"和 §2.5 的"主动运输"使用相同的 Markov 状态图范式——这是第 3 章单离子通道 Markov 模型、Ch5 HH 门控、Ch7 IP3R、Ch8 突触囊泡释放、Ch9 簇发放、Ch10 转录调控的共同技术。 3. §2.6 的"膜电路模型"是 Ch5 HH、Ch12 心脏动作电位、Ch20 耳蜗放大器等所有兴奋性问题的核心。

值得专门提出的是：§2.4.1 末段"两种解释（群体 vs 单分子）"是 Ch3 单通道分析的伏笔——单通道记录（patch clamp）直接观察单分子状态转换，必须用概率解释而非浓度解释。Ch2 → Ch3 的过渡是把 ODE 框架拓展到随机 ODE 的关键。

与第 1 章的衔接：

§2.7 渗透压公式 $\pi_s v = nRT$ 在 Ch1 的化学势语言下完全自然（(2.135)–(2.138)），但 Ch2 把它用到细胞生理学（红细胞破裂、肾近端小管重吸收 Ch17）。Ch1 末段 QSSA + Ch2 末段 Markov 等待时间分布在 Ch3 合并为"准稳态化简的单通道动力学"——这是学习曲线中最陡的一段。

与后续章节的衔接：

Ch3（膜离子通道）：GHK 电流方程 (2.123) → Poisson-Nernst-Planck；载体介导 (2.54) → 单通道 Markov；Na⁺-Ca²⁺ 交换器 (2.79) → 心脏电生理 (Ch12)。
Ch5（兴奋性）：膜电路 (2.130) → HH 方程；体积控制 (2.168) → 阈值与发放。
Ch7（钙动力学）：促进扩散 (2.40) → 钙缓冲；Na⁺-Ca²⁺ 交换器 → 钙外排；Ca²⁺ ATPase → SERCA 泵。
Ch11（循环）：Nernst-Planck (2.114) → 动脉血流的氧扩散。
Ch12（心脏）：Na⁺-K⁺ ATPase (2.100) → 心脏代谢需求。
Ch17（肾）：渗透 (2.139) → 逆流倍增机制。

值得重新推导的方程：

(2.40) 促进扩散增强因子 $\mu\rho$：从 (2.32)–(2.34) 的守恒律 (2.36) 出发，代入 QSSA (2.39) 后重新推导；体会"两个守恒律 + QSSA"如何给出 $1 + \mu\rho$ 的结构。
(2.79) Na⁺-Ca²⁺ 电生交换器流量：从 4 态 Markov 出发，结合详细平衡 (2.78)，重新推导带电压的稳态流。注意 (2.78) 来自 §2.4.3 末段的详细平衡——这是热力学一致性的具体体现。
(2.97) 带电压的 ATPase 流量：从 §2.5.1 推导 + 自由能分解 (2.95) 出发。注意 (2.96) 是"通用"的——任何 ATPase 都满足，关键是电压依赖可以分配到不同步骤（§2.5.2 末段）。
(2.165)–(2.166) pump-leak 模型的稳态解：从 4 维 ODE + 电中性 (2.159) 出发，引入无量纲后化简为 2 个代数方程。这是用无量纲化把高维 ODE 化简为代数方程的范例。
(2.213) Gillespie 算法：从等待时间分布 (2.203) 出发做反变换积分。这是 5 行推导但能解决"如何高效模拟 Markov 过程"——值得手算一次。
(2.238) 平均首次穿越时间：从 (2.237) 边界值问题出发，积分两次即可。对"分子到达靶点需要多久"这类问题有普适价值。

对作者建模选择的评价：

优点：（1）把"细胞稳态"分解为"扩散 + 转运 + 渗透"三件事，每件用守恒律 + 本构关系处理，逻辑清晰；（2）§2.4–2.5 的 Markov 模型展示了"用统一框架处理被动/主动/载体介导"；（3）§2.6 的"GHK 方程 + 膜电路"是电生理学的标准入门，简洁完整；（4）§2.8 的 pump-leak 模型优雅地展示了"为什么 Na⁺/K⁺ 比例不同时体积可以稳定"。
不足：（1）§2.7 渗透的"van't Hoff 公式仅是理想近似"在 §2.7 中段承认，但之后直接用理想公式——没有给出对真实细胞修正的完整讨论；（2）§2.8.1 的电中性假设是基于"Debye 层电荷远小于离子总电荷"（$10^{-7}$ 量级），但没有讨论在强外加电场下这个假设会失效；（3）§2.9 附录虽然提到 Fokker-Planck 但主体仍是确定性——随机建模只在 (2.218)–(2.226) 出现，与主体内容衔接不够；（4）没有给出 Ch3 单通道 patch clamp 数据的随机分析，这是 Ch2 → Ch3 的最大缺口。

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

促进扩散的"低氧增强"机制：肌红蛋白的存在让 O₂ 通量增强 560 倍，但 $D_c \ll D_s$——是否有"最大增强因子"的理论上限？给定 $D_s, D_c, e_0, K, s_0, s_L$ 的范围？
Na⁺-Ca²⁺ 交换器"3:1 化学计量"的生理约束：为什么是 3 Na⁺ 换 1 Ca²⁺ 而不是 2:1 或 4:1？是否与电生性质的最优性有关？Läuger 1991 等已讨论过这个。
GHK vs Eyring（单通道 barrier model）的边界：给定 $V$ 和离子浓度，何时 GHK 足够、何时必须用 PNP 或 Eyring？这是一个没有简单答案的工程判断。
细胞体积控制模型的"压力驱动"缺失：原书假设"机械压力差为零"且"弹性回复力可忽略"，但真实细胞有细胞骨架的弹性恢复（特别是红细胞）。当这两个假设失效时，(2.153) 应该改成什么形式？
§2.9.7 Fick vs Fokker-Planck 的实际后果：细胞内 $D$ 不均匀是普遍现象。用 Fick 而非 FP 形式会引入多大误差？是否有可观测的实验后果？

重要参考文献

[X1] Fick A. On liquid diffusion. Philos Mag 1855;10:30–39. （Fick 定律的原始推导。）

[X2] Einstein A. On the theory of the Brownian motion. Ann Phys 1906;19:371–381. （Einstein 公式 $D = kT/f$ 的推导。）

[X3] Nernst W. Die elektromotorische Wirksamkeit der Ionen. Z Phys Chem 1889;4:129–181. （Nernst 方程。）

[X4] Goldman DE. Potential, impedance, and rectification in membranes. J Gen Physiol 1943;27:37–60. （GHK 电流方程。）

[X5] Hodgkin AL, Katz B. The effect of sodium ions on the electrical activity of the giant axon of the squid. J Physiol 1949;108:37–77. （Hodgkin-Katz 方程——GHK 方程的拓展。）

[X6] Planck M. Ueber die Erregung von Electricität und Wärme in Electrolyten. Ann Phys Chem 1890;39:161–186. （Nernst-Planck 方程的早期推导。）

[X7] Wyman J. Facilitated diffusion and the possible role of myoglobin as a transport mechanism. J Biol Chem 1966;241:115–121. （肌红蛋白促进扩散的实验证据。）

[X8] Wittenberg BA, Wittenberg JB, Caldwell PRB. Role of myoglobin in the oxygen supply to red skeletal muscle. J Biol Chem 1975;250:9038–9043. （促进扩散参数 $k_+, k_-$ 的实验测定——§2.3 的数值来源。）

[X9] Koefoed-Johnsen V, Ussing HH. The nature of the frog skin potential. Acta Physiol Scand 1958;42:298–308. （上皮细胞 Na⁺ 转运与体积控制的经典模型——KJU 模型。）

[X10] Tosteson DC, Hoffman JF. Regulation of cell volume by active cation transport in high and low potassium sheep red cells. J Gen Physiol 1960;44:169–194. （红细胞体积控制——pump-leak 模型的早期工作。）

[X11] Hilgemann DW. New insights into the molecular mechanism of the electrogenic Na/Ca exchange. Ann N Y Acad Sci 2004;976:127–138. （Na⁺-Ca²⁺ 交换器详细动力学。）

[X12] Smith NP, Crampin EJ. Development of models of active ion transport for whole-cell modelling: cardiac sodium-potassium pump as a case study. Prog Biophys Mol Biol 2004;85:387–405. （Post-Albers Na⁺-K⁺ ATPase 的现代数学实现。）

[X13] Apell HJ. Structure-function relationship in P-type ATPases—a biophysical approach. Annu Rev Biophys Biophys Chem 1989;18:221–246.

[X14] Post RL, Hegyvary C, Kume S. Activation by adenosine triphosphate in the phosphorylation kinetics of sodium and potassium ion transport adenosine triphosphatase. J Biol Chem 1972;247:6530–6540. （Post-Albers 模型的早期形式。）

[X15] Gillespie DT. Exact stochastic simulation of coupled chemical reactions. J Phys Chem 1977;81:2340–2361. （Gillespie 算法原始论文。）

[X16] Van Milligen BP, Bons PD, Carreras BA, Sánchez R. On the applicability of Fick's law to diffusion in inhomogeneous systems. Eur J Phys 2006;26:913–925. （§2.9.7 的实验基础——证明 Fick 在 $D$ 不均匀时失效。）

[X17] Segel IH. Biochemical Calculations (2nd ed.). Wiley, 1976. （van't Hoff 公式和缓冲扩散的标准参考。）

[X18] Hille B. Ion Channels of Excitable Membranes (3rd ed.). Sinauer, 2001. （Ch3 单通道分析的经典教科书——§2.4.1 Markov 模型部分的延伸。）

[X19] Murray JD. Mathematical Biology (3rd ed., 2 vols). Springer, 2002/2003. （Ch1 已引；Ch2 的反应-扩散和上皮细胞部分与之有重叠。）

[X20] Gardiner CW. Handbook of Stochastic Methods (3rd ed.). Springer, 2004. （§2.9 附录的完整参考。）