跳转至

I 卷总结：细胞生理学（Chapters 1-10）

概述

Keener & Sneyd《Mathematical Physiology》第一卷（Cellular Physiology, Interdisciplinary Applied Mathematics Vol. 8/I, 2nd ed., Springer 2009）系统建立了细胞生理学的数学建模框架。本卷涵盖 10 章内容，从最基础的化学酶动力学出发，层层递进到基因调控与细胞周期，构成一个从分子到细胞的完整数学建模体系。

整体方法论：本卷始终贯穿"理论 + 实验"的双向循环——每个模型都基于实验数据（生化反应、电压钳、Ca²⁺ 测量、基因表达），每个模型都预测新的可测试现象。数学不是目的，而是理解生物系统的工具。

章节摘要

Ch1：生物化学反应（Biochemical Reactions）

  • 质量作用定律 + 热力学详细平衡
  • 酶动力学：Michaelis-Menten、Hill 方程、Goldbeter-Koshland 开关函数
  • 协同性：MWC vs KNF 模型
  • 核心模型：Sel'kov 糖酵解振荡、Goldbeter-Lefever 模型（Hopf 分岔产生振荡）
  • 数学基础：常微分方程（ODE）、准稳态近似（QSSA）、双稳态与 Hopf 分岔

Ch2：细胞稳态（Cellular Homeostasis）

  • 细胞膜结构：脂双层 + 离子通道
  • 扩散方程：Fick 定律、爱因斯坦 \(D = kT/(6\pi\eta r)\)
  • 缓冲扩散：肌红蛋白 O₂ 运输
  • 载体介导运输：葡萄糖（红细胞）、Na/Ca 交换器、ATPase（Post-Albers 模型）
  • 膜电位：Nernst 方程、Goldman-Hodgkin-Katz (GHK) 方程、电缆方程
  • 细胞体积控制：pump-leak 模型、KJU 上皮模型
  • 随机过程附录：Markov、Gillespie、Fokker-Planck

Ch3：膜离子通道（Membrane Ion Channels）

  • 电流-电压关系：线性 vs GHK
  • 独立性原理 + Ussing 通量比
  • 通道模型：2 态 Markov、饱和结合模型、Eyring barrier 模型
  • 多离子通道：Läuger 摩阻模型、long-pore limit
  • 电压门控：HH 风格的门控 + Boltzmann 因子
  • 泵/交换器：Na/Ca 交换器的电压依赖

Ch4：神经元的被动电流传导（Passive Electrical Flow in Neurons）

  • 电缆方程：\(C_m \partial V/\partial t = V/(R_m) + \partial^2 V/(R_c \partial x^2)\)
  • 空间常数 \(\lambda_m\)、时间常数 \(\tau_m\)
  • 4 类边界条件、输入电阻
  • Rall 理论：分支结构 + 等效圆柱 + \(d_0^{3/2} = \sum d_i^{3/2}\)
  • 突触输入的逆问题（逆卷积）
  • 自伴随算子的 Hilbert 空间方法

Ch5：兴奋性（Excitability）

  • Hodgkin-Huxley 1952 模型（核心）：\(C_m dV/dt = -\bar g_K n^4 (V-V_K) - \bar g_{Na} m^3 h (V-V_{Na}) - \bar g_L (V-V_L) + I_{app}\)
  • 4 维 ODE（V, m, h, n）+ 经验拟合函数 \(\alpha_w, \beta_w\)
  • 动作电位的 4 阶段：upstroke → excited → refractory → recovery
  • 相平面分析：快相平面（3 伪稳态）+ 快-慢相平面（cubic 零线）
  • FitzHugh-Nagumo 简化：\(\epsilon dv/dt = f(v) - w + I\), \(dw/dt = v - \gamma w\)
  • Morris-Lecar、Pushchino、van der Pol 变体
  • 振荡 = Hopf 分岔 + 自持

Ch6：行波在兴奋系统中的传播（Wave Propagation in Excitable Systems）

  • 双稳态方程：\(\partial V/\partial t = \partial^2 V/\partial x^2 + f(V)\)，行波 ODE
  • Shooting 方法证明唯一性
  • 髓鞘化神经：跳跃式传导（saltatory）+ 多发性硬化
  • 离散双稳态方程：耦合阈值 \(D^* = \alpha(1-\alpha)/(2\alpha-1)^2\) → 传播失败
  • 行波脉冲：奇异性摄动分析（外方程 + 内方程）
  • HH 行波：shooting → \(c = 21\) mm/ms（0.5% 误差 vs 实验）
  • 色散曲线 + 运动学理论
  • 高维行波：eikonal-curvature 方程、螺旋波、卷波
  • 应用：心律失常 = 螺旋波、3D 卷波

Ch7：钙动力学（Calcium Dynamics）

  • 钙的"工具箱"：IP3R、RyR、SERCA、PMCA、NCX、线粒体、缓冲
  • 4 个数量级梯度（胞外 1 mM vs 胞内 0.1 μM）
  • IP3R 8 态模型（DeYoung-Keizer 1992）：3 亚基 + 激活/失活
  • IP3R 6 态饱和结合模型（Sneyd-Dufour 2002）
  • 钙振荡模型：2 变量 ODE + Hopf 分岔
  • 牛蛙交感神经元 CICR（Friel 1995）
  • 钙波 + 螺旋波（Xenopus）+ 双参数分岔图
  • 钙缓冲：有效扩散系数、移动 vs 固定缓冲
  • 离散钙源：fire-diffuse-fire 模型
  • 随机钙火花：IPR Markov + Gillespie
  • 细胞间钙波：Sneyd 1994 + Höfer 2002 部分再生 + 肝细胞协调振荡
  • 多尺度均匀化（§7.8 附录）：非均匀介质的有效扩散张量

Ch8：细胞间通讯（Intercellular Communication）

  • 化学突触：量子释放（del Castillo-Katz 1954 + Poisson 极限）
  • Llinás 1976 突触前 5 亚基 Ca²⁺ 通道 + GHK 单通道电流
  • Facilitation 3 个假说：残余自由 Ca²⁺、缓冲饱和、残余结合 Ca²⁺
  • Magleby-Stevens 1972 ACh 动力学 + 两重 QSSA
  • AMPA 6 态 Markov + 简化
  • Gap junction 有效扩散：\(1/D_e = 1/D + 1/(LF)\)（电阻加法）
  • 多尺度均匀化
  • Gap junction 分布的影响：聚集 = 慢转移（"看得见的不等于高通透"）

Ch9：神经内分泌细胞（Neuroendocrine Cells）

  • 胰岛 β 细胞 Chay-Keizer 1983 模型：4 电流 + Ca²⁺ 方程
  • Rinzel 1985 双稳态 + 快慢分离 + Z 形曲线 + 滞后环
  • Phantom bursting（Bertram-Sherman 2004）：ER Ca²⁺ 是真正慢变量
  • Compound bursting + 糖酵解（Wierschem-Bertram 2004）：Ch1 + Ch5 + Ch7 整合
  • Supercell 模型（Sherman 1988）：单细胞随机 + 簇内确定性 + 最优耦合
  • Hindmarsh-Rose 多项式（1984）：通用 bursting 模型
  • Plant 1981 抛物型 bursting（Aplysia R-15）：没有双稳性 + 2 慢变量
  • Li 1997 促性腺细胞球对称模型：Ca²⁺ 空间梯度
  • Bursting 分类：Type I（双稳性）、Type II（两慢变量）、Type III（亚临界 Hopf）

Ch10：细胞功能调控（Regulation of Cell Function）

  • trp 操纵子：负反馈（唯一稳态）
  • lac 操纵子：双稳态（遗传开关，正反馈 + 滞后环）
  • Goodwin 1965 振荡器：级联负反馈，\(n \ge 3, p > 8\) 才振荡
  • Goldbeter 1995 5 变量昼夜节律：Hill 系数 4 + 磷酸化级联
  • Tyson 1999 2 变量相平面：canonical excitable
  • 细胞周期双稳态（Tyson-Novak）：CycB vs Cdh1 互抑 = Start/Finish 不可逆
  • 6 变量细胞周期：S 形曲线 + Hopf 极限环
  • 裂殖酵母模型 + wee1⁻ 突变预测
  • Xenopus 卵母细胞：4 磷酸化状态 + 2 变量极限环

核心主题与跨章联系

1. 双稳态（Bistability）作为通用结构

双稳态是本卷最常出现的几何结构： - Ch1：Goldbeter-Koshland 开关函数 + 自催化 - Ch5：FHN 的稳态双稳性 → 兴奋性 - Ch6：双稳态方程的行波存在性 - Ch7：钙振荡的 N 形零线 + 立方 - Ch9：β 细胞簇发放（双稳态 + 慢变量调制） - Ch10：lac 操纵子（遗传开关）、细胞周期（Start/Finish 不可逆）

通用机制：正反馈 + 滞后环 = 双稳态。这是从基因调控到动作电位到簇发放的"通用结构"。

2. 几何与抽象（FHN 范式）

FitzHugh-Nagumo 简化是本卷的"概念框架"： $\(\epsilon \frac{dv}{dt} = f(v) - w + I, \quad \frac{dw}{dt} = v - \gamma w\)$ 其中 \(f(v)\) 是 cubic 零线，\(w\) 是单调恢复零线。 - Ch5：FHN 本身 - Ch6：FHN 行波 - Ch7：钙振荡（N 形 c 零线 + 单调 y 零线 = FHN 几何） - Ch9：β 细胞快子系统（FHN 几何） - Ch10：Tyson 昼夜节律 2 变量（canonical excitable）

这是"几何方法"在生物物理的胜利——所有兴奋系统共享同一几何。

3. 快慢分离（Fast-Slow Separation）

多个时间尺度的分离是本卷的核心方法： - Ch1.6：数学背景（基本技术、渐近分析、准稳态与奇异性摄动） - Ch2：随机 Markov + QSSA - Ch5：HH 4 维 → FHN 2 维 - Ch6：行波 = 外方程（慢） + 内方程（快） - Ch7：IP3R 8 态 → 6 态（QSSA） + 钙振荡（c 慢，V 快） - Ch9：β 细胞（V, n 快，c 慢）+ Aplysia（V 快，x, c 慢） - Ch10：trp / lac（mRNA 快，蛋白质慢）

快慢分离让所有"复杂"现象化简为"低维几何"——这是 I 卷方法论的核心。

4. 随机 vs 确定性

本卷从不回避随机性： - Ch2.9：Markov 过程 + Gillespie 方法 + Fokker-Planck + 首次穿越时间 - Ch3：单通道 Markov 跳变 - Ch7.6：随机钙火花（IPR Markov + Gillespie + Langevin） - Ch9：supercell 模型（随机通道共享 + 簇内平均）

关键发现（Ch7.6.2, Falcke 2003）：在确定性模型非振荡的 [IP3] 区域，随机模型仍有"准周期"波。这意味着"实验观测到振荡 ≠ 确定性模型有极限环"——对系统生物学有深远影响。

5. 多尺度与均匀化

跨尺度的有效参数是本卷处理复杂系统的方法： - Ch7.8：均值场方程（多尺度均匀化），有效扩散张量 - Ch8.2：Gap junction 有效扩散 - Ch9.1.4：细胞簇的 supercell 简化

关键定理（Ch7.8）：\(D_{eff} = (1/V_c) \int (I + \nabla_\xi W) dV\)——有效参数是微观结构的函数。这是 Ch12 心肌组织双域方程的预备。

数学工具清单

工具 章节 应用
质量作用 + 详细平衡 Ch1, Ch2, Ch7, Ch8 化学反应
Michaelis-Menten / Hill 方程 Ch1, Ch3, Ch7, Ch9, Ch10 酶动力学
Goldbeter-Koshland 开关 Ch1, Ch9, Ch10 双向调控
QSSA + 奇异性摄动 Ch1, Ch2, Ch5, Ch6, Ch7, Ch9, Ch10 多时间尺度
准稳态 + Markov Ch2, Ch3, Ch7, Ch8 单通道
Gillespie / Langevin Ch2, Ch7 随机模拟
Nernst / GHK 方程 Ch2, Ch3, Ch8 膜电位
电缆方程 Ch2, Ch4, Ch5, Ch6 神经元
FHN 几何 Ch5, Ch6, Ch7, Ch9, Ch10 兴奋性
Hopf 分岔 Ch5, Ch7, Ch9, Ch10 振荡
同宿分岔 Ch5, Ch6, Ch9 簇发放
行波 ODE + shooting Ch6 传播
奇异性摄动（外/内方程） Ch6 行波分析
eikonal-curvature 方程 Ch6 高维行波
相平面分析 Ch4, Ch5, Ch7, Ch9, Ch10 稳定性
离散 → 连续均匀化 Ch7.8, Ch8.2 多尺度
自伴随算子 Ch4 Hilbert 空间方法
泊松分布 Ch8 量子释放
双稳态 + 滞后环 Ch7, Ch9, Ch10 不可逆转换
守恒律 + 物质平衡 Ch2, Ch3, Ch7 钙动力学
配体-门控（Ch3.5.4） Ch1, Ch3, Ch7, Ch8 通道动力学

关键洞见（10 条）

  1. Hodgkin-Huxley 1952 是生理学最成功的生物物理模型（Ch5）——4 维 ODE 精确预测动作电位的形状、传导速度、阈值、不应期——0.5% 误差 vs 实验。70 年持续验证。

  2. FitzHugh-Nagumo 是兴奋性的"通用模型"（Ch5）——2 维 ODE 提取了兴奋性的核心几何（cubic 零线 + 时间尺度分离）。所有兴奋系统（HH, FHN, 钙振荡, 昼夜节律, 簇发放）共享此结构。

  3. 钙振荡 = 兴奋系统（Ch7）——N 形 c 零线 + 单调 y 零线 = FHN 几何。双参数分岔图（s, p 平面）显示 C 形同宿 + U 形 Hopf——与 FHN 完全相同。

  4. 随机可以"创造"振荡（Ch7.6.2, Falcke 2003）——确定性模型非振荡区域，随机模型仍有准周期。这意味着"实验观测到振荡 ≠ 确定性模型有极限环"。

  5. 双稳态 = 不可逆转换（Ch9, Ch10）——正反馈 + 滞后环 = 双稳态。从 lac 操纵子到细胞周期，这是"通用开关结构"。

  6. 快慢分离 = 高维到低维（Ch5-Ch10）——HH 4 维 → FHN 2 维、细胞周期 6 维 → 2 维慢变量、Xenopus 9 维 → 2 维极限环。这是 I 卷方法论的核心。

  7. Gap junction 分布"反直觉"（Ch8.2.4）——聚集的 plaque 反而降低通透——"看得见的不等于高通透"。这是"集体行为不能从个体推断"的范例。

  8. 细胞簇 = 统计平均 = 确定性（Ch9.1.4）——N_cell = 1 → 随机，N_cell → ∞ → 确定性。这是大数定律在细胞生物学的应用。

  9. 多尺度均匀化 = 非均匀介质的有效参数（Ch7.8）——\(D_{eff}\) 是微观结构的函数。这是 Ch12 心肌组织双域方程的预备。

  10. "理论 + 实验"双向循环贯穿全卷——每个模型都基于实验数据，每个模型都预测新的可测试现象。这是生物数学建模的核心。

与 II 卷的衔接

I 卷建立的所有工具，在 II 卷（系统生理学）中有直接应用：

  • Ch11（循环系统）：心输出量调控涉及基因表达（Ch10）+ Ca²⁺ 信号（Ch7）+ 膜电位（Ch5）。
  • Ch12（心脏电生理）：心肌组织是Ch5 HH + Ch6 行波 + Ch7 钙火花 + Ch8 gap junction + Ch9 bursting + Ch10 基因的综合。Ch7.8 多尺度均匀化是 Ch12 双域方程的预备。
  • Ch13-15（代谢、肾、肌肉）：代谢酶的基因调控是Ch10 的直接应用。
  • Ch16+（内分泌、生殖）：激素信号通过Ch10 基因表达 + Ch7 Ca²⁺ 实现长期效应。
  • Ch17-19（神经生理）：神经可塑性、LTP/LTD 涉及Ch10 基因表达 + Ch8 突触 + Ch5 动作电位。
  • Ch20（视觉系统）：光感受器适应涉及Ch5 HH 风格 + Ch7 Ca²⁺ + Ch10 基因。

I 卷的 10 章是II 卷所有应用章的细胞基础。

建议阅读顺序

对于想快速入门的读者： 1. Ch1.1-1.5：酶动力学 + Hill 方程 2. Ch2.1-2.5：扩散 + 膜电位 3. Ch5：HH 模型（最经典） 4. Ch6.1-6.2：行波基础 5. Ch7.1-7.3：钙动力学 6. Ch10.3：细胞周期（双稳态的精彩应用）

对于想深入建模的读者： 1. Ch1.6：奇异性摄动 + Goldbeter-Lefever 2. Ch2.6-2.9：载体介导运输 + 随机 Markov 3. Ch3 全部：通道 Markov 的完整工具 4. Ch5 + Ch6：HH + 行波（电生理核心） 5. Ch7 全部：钙信号 + 缓冲 + 多尺度 6. Ch9 全部：神经内分泌（多振荡机制耦合） 7. Ch10 全部：基因调控（双稳态 + 滞后环）

笔记与原书关系

本笔记严格遵循 Keener-Sneyd 2009 第二版的章节结构，每个章节的笔记都对应原书的一章。笔记的 7 节结构严格遵循 SOP §2.2：

  1. 作者（Author）
  2. 内容概述（Content Overview）
  3. 核心方程与概念（Core Equations and Concepts）
  4. 关键结论（Key Conclusions）
  5. 挑战和开放性问题（Challenges and Open Questions）
  6. 个人反思与批判性分析（Personal Reflections and Critical Analysis）
  7. 重要参考文献（Key References）

笔记的 7 节 + 概述补充 5 节 + 内容概述补充 5 节结构让每章笔记既有 SOP 合规性（7 节标准），又有足够的"个人化"分析（反思、批判、对后续章节的展望）。

总结

Keener-Sneyd 第一卷是细胞生理学数学建模的"百科全书"。从 Ch1 的化学酶动力学到 Ch10 的细胞周期，本卷建立了从分子到细胞的完整数学工具链。核心洞见是：

  • 几何方法（FHN 范式、cubic 零线、N 形零线）比具体参数更重要
  • 快慢分离让高维系统化简为低维几何
  • 双稳态是"不可逆转换"和"开关"的核心
  • 随机性与确定性互补
  • 多尺度均匀化是处理非均匀介质的统一方法
  • "理论 + 实验"的双向循环贯穿全卷

I 卷的 10 章为 II 卷的所有应用章奠定了数学基础。这是数学建模在生物物理学的经典教科书。

重要参考文献

[X1] Keener JP, Sneyd J. Mathematical Physiology (2nd ed., Vol. I: Cellular Physiology). Springer Interdisciplinary Applied Mathematics Vol. 8/I, 2009. （I 卷全书。）

[X2] Hodgkin AL, Huxley AF. A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol 1952;117:500–544. （HH 模型原始论文——Ch5 核心。）

[X3] FitzHugh R. Impulses and physiological states in theoretical models of nerve membrane. Biophys J 1961;1:445–466. （FHN 简化——Ch5 简化基础。）

[X4] Rinzel J. Bursting oscillations in an excitable membrane model. In: Ordinary and Partial Differential Equations, Springer Lecture Notes in Mathematics 1151, 1985, pp. 304–316. （β 细胞快慢分离——Ch9 核心。）

[X5] Chay TR, Keizer J. Minimal model for membrane oscillations in the pancreatic β-cell. Biophys J 1983;42:181–190. （Chay-Keizer β 细胞模型——Ch9 起点。）

[X6] De Young GW, Keizer J. Single inositol 1,4,5-trisphosphate receptor channel properties. Proc Natl Acad Sci USA 1992;89:100–115. （IP3R 8 态模型——Ch7 标准。）

[X7] Sneyd J, Dufour J-F. A dynamic model of the type-2 inositol trisphosphate receptor. Proc Natl Acad Sci USA 2002;99:2398–2403. （IP3R 6 态饱和结合模型——Ch7 改进。）

[X8] del Castillo J, Katz B. Quantal components of the end-plate potential. J Physiol 1954;124:560–573. （量子释放——Ch8 基础。）

[X9] Magleby KL, Stevens CF. A quantitative description of end-plate currents. J Physiol 1972;223:173–197. （ACh 动力学——Ch8 标准。）

[X10] Bertram R, Sherman A. Negative calcium feedback: the road to phantom bursting. J Biol Phys 2004a;30:133–152. （Phantom bursting——Ch9 重要更新。）

[X11] 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 模型——Ch9 关键。）

[X12] Plant RE. Bifurcation and resonance in a model for bursting nerve cells. J Math Biol 1981;11:15–32. （Aplysia 抛物型 bursting——Ch9 重要变体。）

[X13] Tyson JJ, Novak B. Regulation of the eukaryotic cell cycle: molecular antagonism, hysteresis, and irreversible transitions. J Theor Biol 2001;210:249–263. （细胞周期双稳态——Ch10 核心。）

[X14] Tyson JJ, et al. A detailed mathematical model for the mammalian circadian clock. Biophys J 1999;77:2411–2417. （昼夜节律 2 变量——Ch10 简化典范。）

[X15] Jacob F, Monod J. Genetic regulatory mechanisms in the synthesis of proteins. J Mol Biol 1961;3:318–356. （lac / trp 操纵子——Ch10 基础。）

[X16] Goodwin BC. Oscillatory behavior in enzymatic control processes. In: Advances in Enzyme Regulation, Weber G (ed.), Pergamon, 1965, pp. 425–438. （Goodwin 振荡器——Ch10 起点。）

[X17] Novak B, Tyson JJ. Numerical analysis of a comprehensive model of M-phase control in Xenopus oocyte extracts and intact embryos. J Cell Sci 1993;106:1153–1168. （Xenopus 卵母细胞 9 变量——Ch10 详细模型。）

[X18] Watson JD, Crick FHC. A structure for deoxyribose nucleic acid. Nature 1953;171:737–738. （DNA 双螺旋——Ch10 分子生物学基础。）

[X19] Alberts B, et al. Molecular Biology of the Cell, 3rd ed. Garland, 1994. （分子生物学标准教材——Ch10 推荐。）

[X20] Murray JD. Mathematical Biology (3rd ed., 2 vols). Springer, 2002/2003. （行波传播的标准教材——Ch6 推荐。）

[X21] Strogatz SH. Nonlinear Dynamics and Chaos. Addison-Wesley, 1994. （非线性动力学入门——Ch5 Ch7 Ch9 推荐。）

[X22] Fall CP, Marland ES, Wagner JM, Tyson JJ (eds.). Cellular Biophysics. MIT Press, 2002. （Ch10 节选的细胞生物物理教材。）

[X23] Hille B. Ion Channels of Excitable Membranes, 3rd ed. Sinauer, 2001. （通道物理学的标准教材——Ch3 Ch5 推荐。）

[X24] Ermentrout GB, Kopell N. Introduction to dynamical systems and bifurcations. In: Methods in Neuronal Modeling, 2nd ed., Koch C, Segev I (eds.), MIT Press, 1998. （分岔理论在神经科学——Ch5 Ch6 推荐。）

[X25] Sneyd J, Falcke M, et al. A comparative study of mathematical models of the inositol trisphosphate receptor. Prog Biophys Mol Biol 2004;89:136–143. （IP3R 模型对比——Ch7 §7.2.5 引用。）

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

[X27] 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. （Fig. 9.18 bursting 细胞类型汇总——Ch9 §9.1.6 引用。）

[X28] Harris-Warrick RM, et al. (eds.). Dynamic Biological Networks: The Stomatogastric Nervous System. MIT Press, 1992. （多振荡机制耦合——Ch9 隐含引用。）

[X29] Keener JP, Sneyd J. Mathematical Physiology (2nd ed., Vol. II: Systems Physiology). Springer, 2009. （II 卷——I 卷的应用延伸。）

[X30] Tyson JJ, Novak B. A model for restriction point control of the mammalian cell cycle. J Theor Biol 2004;230:563–579. （哺乳动物细胞周期——Ch10 §10.3.4 引用。）