跳转至

第1章：生物化学反应

作者

本章作者为 James P. Keener（犹他大学数学系）与 James Sneyd（奥克兰大学数学系），合著 Mathematical Physiology, Vol. I: Cellular Physiology（Springer, Interdisciplinary Applied Mathematics Vol. 8/I, 第二版, 2009）。本章为全书开篇，其角色是为后续所有章节的生理建模准备"化学-热力学-酶动力学"三件套基础工具：质量作用定律及其与 Gibbs 自由能的联系、酶动力学的两类近似（平衡近似与准稳态近似）、以及协同性与糖酵解振荡这两个非平凡应用示例。1.6 节（数学背景附录）独立给出了全书反复使用的尺度化（nondimensionalization）、相平面与线性稳定性、分岔理论、渐近分析的标准定义与文献指引；该附录事实上是第二卷（系统生理学）乃至本书之外众多生理学建模教科书的共同前置材料。

内容概述

本章围绕"生物化学反应如何用常微分方程定量描述"这一核心问题，按 由热力学约束到动力学方程、再到具体生理示例 的逻辑展开。具体而言：

第一部分（§1.1–§1.3）建立质量作用定律与热力学之间的桥梁：单分子、双分子、三分子反应的速率正比于反应物浓度的乘积（1.2–1.9 式），其与 Gibbs 自由能变化通过 \(K_{eq} = e^{\Delta G^0 / RT}\)（1.18 式）联系起来；详细平衡条件则要求任何闭合反应环中正向速率常数之积等于反向速率常数之积（1.24–1.25 式）。这一节最后以 ATP 水解（\(\Delta G^0 = -31.0\) kJ/mol）作为生理学入口。

第二部分（§1.4）系统讨论酶动力学：从 Michaelis-Menten 1913 模型的二步反应（\(S + E \rightleftharpoons C \to P + E\)）出发，分别用 平衡近似（§1.4.1，设复合物 \(C\) 与 \(S\) 处于瞬时平衡）和 准稳态近似（§1.4.2，Briggs-Haldane 1925，设 \(dc/dt \approx 0\)）推导出两种形式的 Michaelis-Menten 律 \(V = V_{max} s / (K + s)\)，其中 \(K_1 = k_{-1}/k_1\)（平衡近似）vs. \(K_m = (k_{-1} + k_2)/k_1\)（准稳态近似）。作者明确指出二者的近似假设并不相同（前者要求 \(k_{-1} \gg k_2\)，后者要求酶浓度远小于底物浓度），并用无量纲化与相平面语言重新表述了 QSSA，揭示其本质是奇异性摄动问题（\(\epsilon \ll 1\) 时变量 \(x\) 沿"慢流形" \(\sigma - x(\sigma + \kappa) = 0\) 滑动）。

随后 §1.4.3 讨论酶抑制（竞争性 vs. 变构性，分别改变有效 \(K_m\) 或 \(V_{max}\)）、§1.4.4 讨论协同性（Hill 方程（1.70）、MWC 1965 协同模型（1.76）、KNF 1966 序贯模型），§1.4.5 拓展到可逆酶反应（指出可逆情况下两种近似给出分母不同的结果），§1.4.6 介绍Goldbeter-Koshland 开关函数——这是双酶磷酸化-去磷酸化系统实现超灵敏开关响应的核心工具。

第三部分（§1.5）将前述工具应用于糖酵解振荡这一真实生物现象：从 PFK1 酶的变构调控（ATP 抑制、AMP 激活、fructose-2,6-bisphosphate 激活）出发，建立 Sel'kov（1968）二变量模型与 Goldbeter-Lefever（1972）MWC 框架模型；通过准稳态化简得到 (1.107)–(1.108) 形式的二维常微系统，借助线性化与特征方程 (1.115) 证明存在超临界 Hopf 分岔；最后将实验观察到的两个 Hopf 分岔点（20 mM/hr 与 160 mM/hr，周期比 \(T_1/T_2 \approx 2.7\)）与理论预测（\(T_1/T_2 \approx 4.6\)，由 (1.133) 给出）对照。

第四部分（§1.6）是数学背景附录：尺度化（Lin & Segel 1988）、相平面与线性稳定性（Edelstein-Keshet 1988、Boyce & Diprima 1997）、分岔理论（Strogatz 1994、Strogatz 称 Hopf 是全书"出现最频繁的分岔"）、奇异性摄动（含边界层、内部层、过渡层、多尺度方法）以及酶动力学的 QSSA 严格化（Bowen 1963、Heineken 1967、Rubinow 1973、Segel & Slemrod 1989）。本节最后用 QSSA 系统的两时间尺度分析展示了初始层的"快跳跃到慢流形"现象（Fig. 1.14）。

本章的主要数学结果可总结为： 1. 质量作用定律 + Gibbs 自由能 → 平衡常数 \(K_{eq} = e^{\Delta G^0 / RT}\) 2. 准稳态近似 → Michaelis-Menten 律 \(V = V_{max} s / (K_m + s)\)，其本质是奇异性摄动 3. 协同性 → Hill 方程 (1.70)、MWC 饱和函数 (1.76) 4. 双酶循环 → Goldbeter-Koshland 开关函数 (1.92) 5. 糖酵解模型 → Hopf 分岔导致周期振荡

前置知识：常微分方程（特别是相平面分析、线性稳定性）、基本热力学（Gibbs 自由能）、基本酶学（底物-酶-复合物的概念）。1.6 节的奇异性摄动内容对理解 QSSA 的严格性非常关键，但初读时可仅以"小参数 \(\epsilon\) 给出快/慢分离"作为直觉。

核心方程与概念

1. 质量作用定律与热力学

核心方程（1.2）：单分子二聚反应 \(A + B \xrightarrow{k} C\) 的速率

\[ \frac{d[C]}{dt} = k [A] [B] \]

\(k\) 是速率常数，与分子几何、温度有关。注意：质量作用不是不可违反的物理定律，而是有用的模型；在极高/极低浓度下失效（§1.1 注）。

核心方程（1.18）：与热力学的联系

\[ K_{eq} = \frac{k_-}{k_+} = e^{\Delta G^0 / RT} \]

其中 \(\Delta G^0\) 是标准自由能变化（$G^0 = $ 浓度 1 M 时的 Gibbs 自由能，\(P = 1\) atm, \(T = 25^\circ C\)）。重要（§1.2 末段）：热力学只给两个速率常数的比值，不给出各自的绝对大小——所以 \(\Delta G^0\) 告诉我们反应"倾向"哪边进行，但不告诉反应有多快。这是后续动力学分析独立于热力学的根本原因。

生理实例（1.22–1.23）：ATP 水解 \(\Delta G^0 = -31.0\) kJ/mol，故 \(K_{eq} \approx e^{12.5} \approx 2.7 \times 10^5\)。但活细胞中 [ATP]/[ADP][Pi] 比值远高于平衡值——细胞通过不断消耗 ATP 维持远离平衡状态，这是"自由能转导"（free energy transduction）的本质。

2. 详细平衡

核心方程（1.24–1.25）：闭合反应环 \(A \to B \to C \to A\) 中，

\[ k_1 k_2 k_3 = k_{-1} k_{-2} k_{-3} \quad \text{即} \quad K_1 K_2 K_3 = 1 \]

（其中 \(K_i = k_{-i}/k_i\)）。这与浓度无关，是反应网络热力学一致性的必要条件。违反此条件则反应环可作为"自由能引擎"持续产生净流（如 \(\mathrm{F_0 F_1}\)-ATP 合酶正是利用质子梯度驱动反向反应）。

3. 酶动力学：Michaelis-Menten 律的两条推导路线

反应机制（§1.4 开头）：

\[ S + E \underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}} C \xrightarrow{k_2} P + E \]

设 \(s = [S]\), \(c = [C]\), \(e = [E]\), \(p = [P]\)，守恒 \(e + c = e_0\)（总酶量）。质量作用给出四方程 (1.26)–(1.29)。

路线 A：平衡近似（Michaelis-Menten 1913, §1.4.1）

假设 \(k_1 s e = k_{-1} c\)（复合物与底物瞬时平衡），结合守恒得

\[ c = \frac{e_0 s}{K_1 + s}, \quad V = \frac{dp}{dt} = k_2 c = \frac{V_{max} s}{K_1 + s}, \quad K_1 = \frac{k_{-1}}{k_1}, \ V_{max} = k_2 e_0 \]

路线 B：准稳态近似（Briggs-Haldane 1925, §1.4.2）

引入无量纲变量 \(\sigma = s/s_0\), \(x = c/e_0\), \(\tau = k_1 e_0 t\), \(\kappa = (k_{-1} + k_2)/(k_1 s_0)\), \(\epsilon = e_0/s_0\), \(\alpha = k_{-1}/(k_1 s_0)\)，化简为

\[ \frac{d\sigma}{d\tau} = -\sigma + x(\sigma + \alpha), \quad \epsilon \frac{dx}{d\tau} = \sigma - x(\sigma + \kappa) \]

设 \(\epsilon \ll 1\)（酶浓度远小于底物浓度，典型范围 \(10^{-2}\)–\(10^{-7}\)），令 \(\epsilon \frac{dx}{d\tau} = 0\) 得"慢流形"（slow manifold）\(x = \sigma/(\sigma + \kappa)\)，代入得

\[ \frac{d\sigma}{d\tau} = -q \frac{\sigma}{\sigma + \kappa}, \quad q = \kappa - \alpha = \frac{k_2}{k_1 s_0} \]

转回原变量即

\[ V = \frac{dp}{dt} = \frac{V_{max} s}{s + K_m}, \quad K_m = \frac{k_{-1} + k_2}{k_1} \]

两者关系（§1.4.2 末段）：当 \(k_{-1} \gg k_2\) 时 \(K_m \approx K_1\)，两近似等价。核心区别：QSSA 的有效性来自 \(\epsilon\) 小，而平衡近似的有效性来自 \(k_{-1} \gg k_2\)——这是两种独立但常常同时成立的近似条件。当二者都不成立时（如 \(k_1 e_0\) 不再是 \(k_{-1}, k_2\) 的最大者），需要更细致的奇异性摄动处理（见 Exercise 14 及 Segel 1988）。

实验测定（1.45）：Lineweaver-Burk 图 \(1/V = 1/V_{max} + (K_m/V_{max})(1/s)\) 是经典线性化方法，但实际中常用直接线性图（Eisenthal & Cornish-Bowden 1974）以减小低 \(s\) 处的误差放大。

4. 抑制与协同性

竞争性抑制（§1.4.3 竞争部分）：抑制剂 \(I\) 与底物争夺活性位点

\[ V = \frac{V_{max} s}{s + K_m (1 + i/K_i)}, \quad K_i = \frac{k_{-3}}{k_3} \]

效果：\(K_m\) 增大为 \(K_m(1 + i/K_i)\)，\(V_{max}\) 不变。

变构抑制（§1.4.3 变构部分）：\(I\) 结合于变构位点，4 个状态 \(E, ES, EI, EIS\) 的平衡分析给出

\[ V = \frac{V_{max}}{1 + i/K_i} \cdot \frac{s}{K_s + s} \]

效果：\(V_{max}\) 减小为 \(V_{max}/(1 + i/K_i)\)，\(K_s\) 不变。两种抑制在表观动力学上可以通过 Lineweaver-Burk 图的横/纵截距变化区分。

Hill 方程（1.70，正协同极限）：\(n\) 个独立结合位点 + 极端正协同下

\[ V = \frac{V_{max} s^n}{K_m^n + s^n}, \quad K_m^n = \prod_{i=1}^n K_i \]

取对数得 Hill 系数 \(n\)（通过 \(\ln(V/(V_{max}-V))\) vs. \(\ln s\) 的斜率测定）。注意：\(n\) 通常不是整数，反映多步过程的协同效应而非真正的结合位点数。

MWC 模型（Monod-Wyman-Changeux 1965, 1.72–1.76）：蛋白质有 \(R\)（relaxed）和 \(T\)（tense）两种构象，可有 \(n\) 个结合位点。饱和函数

\[ Y = \frac{s K_1^{-1} (1 + s K_1^{-1})^{n-1} + K_2^{-1} [s K_3^{-1} (1 + s K_3^{-1})^{n-1}]}{(1 + s K_1^{-1})^n + K_2^{-1} (1 + s K_3^{-1})^n} \]

其中 \(K_2 = r_0 / t_0\) 是 \(R \to T\) 构象平衡常数。重要（§1.4.4 MWC 末段）：MWC 模型只能给出正协同，无论 \(K_1 > K_3\) 还是反之。这是因为 \(R\) 构象与 \(T\) 构象之间的协同切换天然是正反馈。

KNF 模型（Koshland-Nemethy-Filmer 1966）：序贯构象变化，可给出负协同。区分 MWC vs. KNF 需要细致实验（如配体平衡时的构象异质性测定）。

5. Goldbeter-Koshland 开关函数

反应（1.83）：\(W \rightleftharpoons W^*\)，\(E_1\) 催化 \(W \to W^*\)，\(E_2\) 催化 \(W^* \to W\)。设两反应均为 Michaelis-Menten 形式

\[ r_1 = \frac{V_1 E_1}{K_1 + W}, \quad r_{-1} = \frac{V_2 E_2}{K_2 + W^*} \]

稳态条件 \(r_1 = r_{-1}\) 给出

\[ \frac{v_1}{v_2} = \frac{(1-y)(\hat K_1 + y)}{y(\hat K_2 + 1 - y)} \]

其中 \(y = W/W_t\), \(\hat K_i = K_i / W_t\), \(v_i = V_i E_i\)。这是 \(y\) 的二次方程，解为 (1.88)–(1.92)：

\[ y = G(v_1, v_2, \hat K_1, \hat K_2) = \beta - \frac{\sqrt{\beta^2 - 4\alpha\gamma}}{2\alpha} \]

\(\alpha, \beta, \gamma\) 见 (1.89)–(1.91)。关键（§1.4.6 末段 + Fig. 1.5）：当 \(\hat K_1, \hat K_2 \ll 1\)（酶反应处于饱和）时，\(y\) 随 \(v_1/v_2\) 呈急剧开关——这是仅靠两个 Michaelis-Menten 反应就能实现"超灵敏"调控的机制，无需大 Hill 系数 \(n\)。该函数在第 10 章（细胞功能调控）的生化网络分析中被广泛使用。

6. 糖酵解振荡的 Sel'kov 与 Goldbeter-Lefever 模型

生理背景（§1.5）：糖酵解三阶段（糖酵解、Krebs 循环、电子传递链），每分子葡萄糖净产 38 ATP。PFK1 酶是关键变构调控点： - ATP 抑制（结合于变构位点） - AMP 激活（与 ATP 竞争变构位点） - fructose-2,6-bisphosphate 激活 - 这构成正反馈（PFK1 活性 → 减少 ATP/AMP 比 → 进一步激活 PFK1）

Sel'kov 模型（1968, 1.97–1.100）：PFK1 (E) 需结合 \(\gamma\) 个 ADP (\(S_2\)) 才激活。无量纲化（1.36）+ 准稳态 (1.105)–(1.106) 后得二维系统

\[ \frac{d\sigma_1}{d\tau} = \nu - f(\sigma_1, \sigma_2), \quad \frac{d\sigma_2}{d\tau} = \alpha f(\sigma_1, \sigma_2) - \eta \sigma_2 \]

特征方程 (1.115) 的根实部符号由

\[ H = \alpha f_2 - \eta - f_1 \]

决定（\(f_i = \partial f/\partial \sigma_i\)）。\(H(\nu) = (1-\nu)[(1+y)(\eta\gamma + (\nu-1)y)] - \eta\) 在 \(\gamma > 1\) 时必有 Hopf 分岔点。定量预测（Fig. 1.9）：\(\nu\) 超过阈值后出现稳定极限环，\(\nu\) 减小时振幅迅速扩张直到无界——与实验"高、低流速下稳定"不完全吻合。

Goldbeter-Lefever 模型（1972, 1.118–1.134）：用 MWC 框架建模 PFK1 的二聚体。14 种化学态，12 个内部复合物在准稳态下化简为 (1.107)–(1.108) 形式但 \(f(\sigma_1, \sigma_2)\) 不同；在额外假设 \(k_3 = 0\), \(k_1 \gg k_{-1}\), \(k \gg k_{-2}\) 下进一步化简为 \(f(\sigma_1, \sigma_2) = \sigma_1(1+\sigma_2)^2\)。稳定性条件 (1.130) 化为三次多项式 (1.131)

\[ \eta y^3 - y + 2 = 0, \quad y = 1 + \nu/\eta \]

由此可同时解出两个 Hopf 分岔点。与实验对比（1.132–1.134）：设 \(\nu \propto\) 葡萄糖供给，要求 \(\nu_2/\nu_1 = 160/20 = 8\)，解得 \(\eta = 116.7, y_1 = 2.08, y_2 = 9.61\)，对应周期比 \(T_1/T_2 = \sqrt{y_2/y_1} = 2.15\)（用 (1.134) 的简化形式）——实验值为 2.7，理论预测 2.7（书中给出 4.6，但用 (1.134) 精确计算应得 \(\sqrt{9.61/2.08} \approx 2.15\)，与 (1.134) 表达式一致——见 §1.5 末段与 Fig. 1.11）。

7. 奇异性摄动与初始层

§1.6.2–1.6.3 是全章的方法论核心。三类奇异性摄动问题： - 边界层（boundary layer）：快变仅在边界附近，如 \(\epsilon u'' + u' + u = 0\) 的解 \(u(x;\epsilon) = (1-e)e^{-x/\epsilon} + e^{1-x} + O(\epsilon)\) - 内部层/过渡层（interior/transition layer）：形如 \(\tan((x-x_0)/\epsilon)\)，对兴奋性系统（第 5 章）至关重要 - 多尺度（multiple scales）：如心肌动作电位（细胞结构远小于波前）、耳蜗膜（波长远短于膜长度变化）

QSSA 的严格化（§1.6.3）：用两时间尺度展开 \(\sigma = \sigma_0 + \epsilon \sigma_1 + \cdots\), \(x = x_0 + \epsilon x_1 + \cdots\)。在慢时间 \(\tau\) 上，\(x_0 = \sigma_0/(\sigma_0 + \kappa)\)（即 (1.149)），但不满足初值 \(x(0) = 0\)。切换到快时间 \(\eta = \tau/\epsilon\) 后 \(\tilde x_0 = (1-e^{-(1+\kappa)\eta})/(1+\kappa)\)（即 (1.155)），匹配得全时间有效解 \(x(\tau) = (\sigma/(\sigma + \kappa))(1 - e^{-(1+\kappa)\tau/\epsilon})\)（即 (1.158)）。这是匹配渐近展开（matched asymptotic expansion）的标准范例。

类比（§1.6.3 末段）：平衡近似对应的慢变量不是 \(\sigma\) 而是 \(\sigma + \alpha x\)（\(\alpha = e_0/s_0\)），即"快平衡"约束了不同的线性组合。

关键结论

  1. 质量作用定律 + Gibbs 自由能 = 平衡约束（1.18）：仅给出速率常数之比，不给出绝对速率。这是热力学与动力学的根本分工。

  2. Michaelis-Menten 律的两种近似本质不同（§1.4.1–1.4.2）：平衡近似假设 \(k_{-1} \gg k_2\)，QSSA 假设 \(\epsilon = e_0/s_0 \ll 1\)。当 \(k_1 e_0\) 不再小于 \(k_{-1}, k_2\) 时（Exercise 14），需要重新构造无量纲方案。

  3. QSSA 是奇异性摄动（§1.6.3 + Fig. 1.14）：解先在快时间尺度"跳到"慢流形 \(\sigma - x(\sigma + \kappa) = 0\)，再沿慢流形滑动到稳态。初始层时间常数约为 \(\epsilon/(1+\kappa)\)，典型 \(10^{-7}\) 量级（\(e_0/s_0\) 在 \(10^{-2}\)–\(10^{-7}\) 之间）——这就是为什么大多数生化反应中初始层"生理上不重要"。

  4. MWC vs. KNF 协同模型的根本差异（§1.4.4）：MWC 是"协同"（concerted）模型——所有亚基同步构象变化——只能给出正协同；KNF 是"序贯"（sequential）模型——单亚基独立变化——可给出正/负两种协同。对氧合血红蛋白的实验证据支持二者兼具。

  5. Goldbeter-Koshland 函数实现"零协同超灵敏"（§1.4.6 + (1.92)）：仅需两个 Michaelis-Menten 酶反应（不需要大 \(n\) 的 Hill 函数），在 \(\hat K_1, \hat K_2 \ll 1\) 时即可获得接近 Heaviside 的开关响应。这是生化网络（特别第 10 章）构建"分子开关"的基础工具。

  6. 糖酵解振荡 = Hopf 分岔（§1.5）：Sel'kov 模型在 \(\gamma > 1\) 时有超临界 Hopf 分岔；Goldbeter-Lefever 模型在 \(\eta\) 足够大时存在两个 Hopf 分岔点，区间内稳定周期解——更接近 Hess-Boiteux 1973 的双分岔实验。理论与实验周期比的偏差（理论 2.7，实验 4.6）暴露了简化的局限。

  7. 分岔理论工具是全书的"通用语"（§1.6.1 末段 + Fig. 1.9, 1.11）：本书所有分岔图均用 XPPAUT (Ermentrout 2002) 绘制，它是 AUTO 的便捷前端。读者应熟悉稳态分岔、Hopf 分岔、鞍结分岔、同宿分岔四种基本类型；其中 Hopf 是本书出现最频繁的。

挑战和开放性问题

  1. "广泛适用性"的隐含局限（§1.1 末段注）：质量作用在极高/极低浓度下失效——但细胞内的局部浓度常处于这两种极端（mM 级别局部富集、信号转导中的 nM 级稀溶液）。原书未给出统一的高浓度/低浓度修正框架。

  2. 细胞内非理想性（§1.2 中段）：非理想溶液应使用化学活度 \(a\) 替代浓度 \(c\)；但活度系数的精确计算在细胞内（高浓度大分子拥挤环境）仍是开放问题。这是连接本章与第 2 章（细胞稳态）"膜内大分子"问题的桥梁，但本章未深入。

  3. QSSA 有效性的边界（Exercise 14 + Segel 1988）：当 \(e_0/s_0\) 不再是远小于 1（如视网膜视锥细胞中 cAMP 水解），原方案失效。已有替代推导，但缺乏统一判据。这是 Segel-Slemrod 1989 后续工作的起点，本书未充分覆盖。

  4. MWC vs. KNF 的实验区分（§1.4.4 末段）：对氧合血红蛋白"两种机制都参与"，但对绝大多数协同酶无法干净地区分。本书未给出"在哪个实验观察上必定有差异"的判据。

  5. Goldbeter-Koshland 的"超灵敏"是否真有优势（§1.4.6 隐含）：函数 (1.92) 在 \(\hat K_1, \hat K_2 \ll 1\) 时是"开关"，但这要求酶永远处于饱和状态——这与"酶工作在线性区间"的传统工程直觉冲突。细胞如何平衡调控灵敏度与代谢通量利用率？本章未讨论。

  6. 糖酵解振荡的起源（§1.5 中段脚注）：PFK 反馈是经典解释，但 hexose 转运动力学（Bier 1996）、ATP 自催化（Reijenga 2002）、整体网络效应都曾被提出作为替代机制。实验上完整细胞 vs. 提取物的振荡行为差异（前者正弦、后者弛豫）也未被单一模型解释。

  7. 初始层生理重要性的判断准则（§1.6.3 末段隐含）："大多数生化反应初始层太快所以不重要"是经验法则。但对于信号转导（如钙火花、IP3 释放）的快速局部动力学，初始层恰恰是核心。本书将这一问题推迟到第 5、7 章，但未给出一个统一的时间尺度判定方法。

个人反思与批判性分析

建模哲学层面：

Keener-Sneyd 的"以奇异性摄动为中心"视角与 Murray 等教科书有显著差异。Murray 的 Mathematical Biology（2002）倾向于"先写出方程再用标准 ODE 理论分析"；Keener-Sneyd 始终强调无量纲化 + QSSA + 多尺度匹配三件套。这在 Ch1 末尾（§1.6.3）的 QSSA 严格化中达到极致——把一个看似平凡的"酶动力学"问题转化为边界层 + 慢流形的完整摄动问题。这种风格贯穿全卷（第 5、6、9、10 章的振荡与脉冲问题均如此处理）。

值得专门提出的是：QSSA 严格化在 Ch1 §1.6.3 的"双时间尺度匹配"（(1.149) + (1.155) → (1.158)）事实上是第 6 章（行波传播）和第 5 章（兴奋性）"快速-慢系统相平面"分析的预演。读者如果在 Ch1 把这个匹配过程彻底理解，Ch5/6 的几何直觉会非常顺。

与第 2 章的衔接：

Ch1 的"ATP 水解 + PFK 反馈"在 Ch2（细胞稳态）会被进一步用于 离子主动转运（Na/K-ATPase）和 膜电位（Nernst 方程）。两个 QSSA 工具（MWC 模型、Goldbeter-Koshland 函数）会在 Ch7（钙动力学）被直接用于 IP3 受体的多状态门控，在 Ch10（细胞功能调控）被用于基因调控网络的双稳态。

值得重新推导的方程：

  • (1.32) vs. (1.43)：QSSA 推导中的无量纲化与"快变量近似"步骤（(1.36)–(1.44)），是奇异性摄动应用的最简范例。
  • (1.76) MWC 饱和函数：从图 1.4 的状态图出发，自己重新推一遍 6 态 / 一般 \(n\) 态，会发现"二项式系数从哪里来"。
  • (1.92) Goldbeter-Koshland 函数：从稳态方程 (1.86) 出发解二次方程，注意根号前的负号是物理选择（\(y\) 应在 \([0,1]\) 范围内），见 Exercise 12 的"Why negative square root"问题。
  • (1.124) Goldbeter-Lefever 化简：从 14 态到 (1.107)–(1.108) 的步骤在 Exercise 18 中让读者重做；这是 AUTO/XPPPAUT 配合手算的典型工作量。

对作者建模选择的评价：

  • 优点：把"准稳态 + 多尺度"作为统一语言贯穿始终，使得跨章节的方法论一致；MWC 和 Goldbeter-Koshland 介绍早且有具体生理应用。
  • 不足：（1）没有足够强调 QSSA 的失败模式——只在 Exercise 14 隐含；（2）糖酵解模型对实验的双 Hopf 周期比预测仍有 ~50% 误差（2.7 vs. 4.6），作者未讨论这一偏差的来源（可能是 \(\nu \propto\) 葡萄糖供给这个假设过强）；（3）§1.6 数学背景过于浓缩——奇异性摄动只用了一段话，但后续章节反复依赖它；建议非 ODE 专业背景的读者预备 Lin-Segel 1988 或 Murray 2002 作为外部参考。

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

  1. §1.4.5 末段说"可逆情况下平衡近似和 QSSA 给出分母不同的结果（(1.80) vs. (1.82)）"——这两种结果在什么实验条件下能被区分？\(\Delta G\) 测量的精度够吗？
  2. MWC 模型的"5 个独立参数"（\(K_1, K_2, K_3, n\)）中，哪些能从单分子 FRET 实验单独标定？目前似乎只能从饱和曲线整体拟合。
  3. Hopf 分岔点 \(\nu_1, \nu_2\) 的精确预测是金标准；Goldbeter-Lefever 的 \(T_1/T_2 = 4.6\) 误差 70%（vs. 实验 2.7）——这是模型结构问题（缺少 fructose-2,6-bisphosphate 调控）还是参数标定问题？

重要参考文献

[X1] Michaelis L, Menten ML. Die Kinetik der Invertinwirkung. Biochem Z 1913;49:333–369. （Michaelis-Menten 二步反应机制与平衡近似。）

[X2] Briggs GE, Haldane JBS. A note on the kinetics of enzyme action. Biochem J 1925;19:338–339. （准稳态近似。）

[X3] Eisenthal R, Cornish-Bowden A. The direct linear plot: a new graphical procedure for estimating enzyme kinetic parameters. Biochem J 1974;139:715–720. （Lineweaver-Burk 的稳健替代。）

[X4] Cornish-Bowden A, Eisenthal R. Statistical considerations in the estimation of enzyme kinetic parameters by the direct linear plot and other methods. Biochem J 1974;139:721–730.

[X5] Monod J, Wyman J, Changeux JP. On the nature of allosteric transitions: a plausible model. J Mol Biol 1965;12:88–118. （MWC 协同模型。）

[X6] Koshland DE, Nemethy G, Filmer D. Comparison of experimental binding data and theoretical models in proteins containing subunits. Biochemistry 1966;5:365–385. （KNF 序贯模型。）

[X7] Goldbeter A, Koshland DE. An amplified sensitivity arising from covalent modification in biological systems. Proc Natl Acad Sci USA 1981;78:6840–6844. （Goldbeter-Koshland 开关函数。）

[X8] Sel'kov EE. Self-oscillations in glycolysis. Eur J Biochem 1968;4:79–86. （糖酵解二变量模型。）

[X9] Goldbeter A, Lefever R. Dissipative structures for an allosteric model. Biophys J 1972;12:1302–1315. （MWC 框架下的糖酵解振荡。）

[X10] Hess B, Boiteux A. Heterogeneous, metabolic oscillations in a cell-free extract of S. carlsbergensis. Hoppe-Seyler's Z Physiol Chem 1973;354:17–20. （糖酵解振荡的原始实验观察。）

[X11] Strogatz SH. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering. Addison-Wesley, 1994. （分岔理论的入门标准教材。）

[X12] Lin CC, Segel LA. Mathematics Applied to Deterministic Problems in the Natural Sciences. SIAM, 1988. （尺度化与奇异性摄动的标准参考。）

[X13] Segel LA, Slemrod M. The quasi-steady-state assumption: a case study in perturbation. SIAM Rev 1989;31:446–477. （QSSA 严格化的经典综述。）

[X14] Ermentrout B. Simulating, Analyzing, and Animating Dynamical Systems: A Guide to XPPAUT for Researchers and Students. SIAM, 2002. （XPPAUT 软件手册——本书所有分岔图绘制工具。）

[X15] Stryer L. Biochemistry (3rd ed.). WH Freeman, 1988. （糖酵解途径的标准生物化学参考——本节部分参数来源。）