第六章：酶动力学 (Enzyme Kinetics)

6.1 酶动力学：基本酶反应

基本概念

生物化学反应的复杂性使得建立简化模型成为理解现象的必要手段。酶（Enzymes）是活细胞中具有显著催化效率的蛋白质，能选择性作用于称为底物（Substrates）的特定化合物。例如，血红蛋白是一种酶，氧气是其底物。酶在调节生物过程中起着重要作用，可以作为激活剂或抑制剂。

Michaelis-Menten基本酶反应机制

1913年Michaelis和Menten提出的最基本的酶催化反应涉及底物S与酶E反应形成复合物SE，复合物随后转化为产物P和酶。反应机理表示为：

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

其中$k_1$、$k_{-1}$和$k_2$是与反应速率相关的常数参数。

质量作用定律与微分方程组

根据质量作用定律（Law of Mass Action），反应速率与反应物浓度的乘积成正比。用小写字母表示各反应物浓度：

\[ s = [S],\quad e = [E],\quad c = [SE],\quad p = [P] \]

应用质量作用定律得到非线性反应方程组：

\[ \frac{ds}{dt} = -k_1 es + k_{-1}c $$ $$ \frac{de}{dt} = -k_1 es + (k_{-1} + k_2)c $$ $$ \frac{dc}{dt} = k_1 es - (k_{-1} + k_2)c $$ $$ \frac{dp}{dt} = k_2 c \]

初始条件为： $$ s(0) = s_0,\quad e(0) = e_0,\quad c(0) = 0,\quad p(0) = 0 $$

酶守恒定律

酶作为催化剂，其总浓度（游离酶+结合酶）守恒。从方程组可得：

\[ \frac{de}{dt} + \frac{dc}{dt} = 0 \Rightarrow e(t) + c(t) = e_0 \]

代入后得到关于s和c的二维系统：

\[ \frac{ds}{dt} = -k_1 e_0 s + (k_1 s + k_{-1})c $$ $$ \frac{dc}{dt} = k_1 e_0 s - (k_1 s + k_{-1} + k_2)c \]

准稳态近似（Quasi-Steady State Approximation）

假设复合物c的形成初始阶段非常快，之后基本处于平衡态，即$dc/dt \approx 0$，由此得：

\[ c(t) = \frac{e_0 s(t)}{s(t) + K_m},\quad K_m = \frac{k_{-1} + k_2}{k_1} \]

其中$K_m$称为米氏常数（Michaelis Constant）。

代入得底物消耗方程：

\[ \frac{ds}{dt} = -\frac{k_2 e_0 s}{s + K_m} \]

求解得隐式解： $$ s(t) + K_m \ln s(t) = s_0 + K_m \ln s_0 $$

无量纲化分析

引入无量纲变量： $$ \tau = k_1 e_0 t,\quad u(\tau) = \frac{s(t)}{s_0},\quad v(\tau) = \frac{c(t)}{e_0} $$ $$ \lambda = \frac{k_2}{k_1 s_0},\quad K = \frac{k_{-1} + k_2}{k_1 s_0} = \frac{K_m}{s_0},\quad \varepsilon = \frac{e_0}{s_0} $$

得标准准稳态近似的无量纲系统： $$ \frac{du}{d\tau} = -u + (u + K - \lambda)v $$ $$ \varepsilon \frac{dv}{d\tau} = u - (u + K)v $$ 初始条件：$u(0) = 1,\; v(0) = 0$

6.2 瞬态时间估计与无量纲化

两个时间尺度

酶促反应中一个关键特征是酶的催化效能极高，酶浓度远小于底物浓度，即$\varepsilon = e_0/s_0 \ll 1$。

快速瞬态时间尺度$t_c$：由$dc/dt$方程估计，设$s(t) = s_0$：

\[ t_c = \frac{1}{k_1(s_0 + K_m)} \]

慢速时间尺度$t_s$：底物浓度显著变化所需时间：

\[ t_s \approx \frac{s_0}{|ds/dt|_{\max}} \approx \frac{s_0 + K_m}{k_2 s_0} \]

准稳态近似的有效性条件

准稳态近似有效的条件是快速瞬态时间远小于慢速变化时间，即$t_c \ll t_s$，这要求：

\[ \frac{k_2 e_0}{k_1(s_0 + K_m)^2} \ll 1 \]

另一个条件是底物在快速瞬态期间的消耗只是初始量的很小一部分：

\[ \varepsilon = \frac{e_0}{s_0 + K_m} \ll 1 \]

内域与外域的无量纲化

针对快速瞬态（内域）采用时间尺度$\tau = t/t_c = tk_1(s_0 + K_m)$，得到内域方程。

针对慢速或外域采用时间尺度$T = (1 + \rho)t/t_s$，得到外域方程。

两种无量纲化方式本质相同，只是小参数$\varepsilon$出现在方程的不同位置。

6.3 Michaelis-Menten准稳态分析

奇异摄动问题

当$0 < \varepsilon \ll 1$时，系统(6.13)是一个奇异摄动问题。设$\varepsilon = 0$会降低微分方程组的阶数，因此不能同时满足所有初始条件。

内部（奇异）解

在$\tau = 0$附近，引入伸展变量$\sigma = \tau/\varepsilon$，使$\varepsilon dv/d\tau = dv/d\sigma$。在此尺度下：

\[ \frac{dU_0}{d\sigma} = 0,\quad \frac{dV_0}{d\sigma} = U_0 - (U_0 + K)V_0 $$ $$ U_0(0) = 1,\quad V_0(0) = 0 \]

解为： $$ U_0(\sigma) = 1,\quad V_0(\sigma) = \frac{1}{1+K}(1 - e^{-(1+K)\sigma}) $$

外部（非奇异）解

令$\varepsilon = 0$，从代数方程$v_0 = u_0/(u_0 + K)$得： $$ u_0(\tau) + K \ln u_0(\tau) = 1 - \lambda\tau $$ $$ v_0(\tau) = \frac{u_0(\tau)}{u_0(\tau) + K} $$

匹配原则

当$\sigma \to \infty$（即$\tau \to 0$）时，内部解应趋近于外部解： $$ \lim_{\sigma \to \infty}V_0(\sigma) = \frac{1}{1+K} = \lim_{\tau \to 0}v_0(\tau) $$

一致有效渐近解（O(1)精度）

\[ u(\tau; \varepsilon) = u_0(\tau) + O(\varepsilon) $$ $$ v(\tau; \varepsilon) = \begin{cases} V_0(\sigma) + O(\varepsilon) = \frac{1}{1+K}\left[1 - \exp\left(-(1+K)\frac{\tau}{\varepsilon}\right)\right] + O(\varepsilon), & 0 \leq \tau \ll 1 \\ v_0(\tau) + O(\varepsilon) = \frac{u_0(\tau)}{u_0(\tau) + K}, & \tau = O(1) \end{cases} \]

Michaelis-Menten速率方程

初始反应速率（O(1)精度）： $$ r_0 = \left|\frac{du_0}{d\tau}\right|_{\tau=0} = \frac{\lambda}{1+K} $$

量纲形式： $$ R_0 = \frac{k_2 e_0 s_0}{s_0 + K_m} = \frac{Qs_0}{s_0 + K_m} $$

其中$Q = k_2 e_0$是最大反应速率，$K_m = (k_{-1} + k_2)/k_1$是Michaelis常数。

这就是经典的Michaelis-Menten摄取方程： $$ \frac{ds}{dt} = -\frac{Qs}{K_m + s} $$

6.4 自杀底物动力学（Suicide Substrate Kinetics）

反应机制

自杀底物是一种酶的"自杀性"抑制剂，其机制为：

\[ E + S \underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}} X \xrightarrow{k_2} Y \xrightarrow{\begin{matrix} k_3 \\ k_4 \end{matrix}} \begin{cases} E + P \\ E_i \end{cases} \]

其中$Y$可沿两条途径反应：生成产物P（速率$k_3$）或生成失活酶$E_i$（速率$k_4$）。比率$r = k_3/k_4$称为分配比（Partition Ratio）。

动力学方程

根据质量作用定律： $$ \frac{d[S]}{dt} = -k_1[E][S] + k_{-1}[X] $$ $$ \frac{d[E]}{dt} = -k_1[E][S] + k_{-1}[X] + k_3[Y] $$ $$ \frac{d[X]}{dt} = k_1[E][S] - k_{-1}[X] - k_2[X] $$ $$ \frac{d[Y]}{dt} = k_2[X] - (k_3 + k_4)[Y] $$ $$ \frac{d[E_i]}{dt} = k_4[Y] $$

关键判定参数

Waley条件：关键参数是$r\mu$，其中$\mu = e_0/s_0$
Tatsunami条件：关键参数是$(1+r)\mu$
$(1+r)\mu > 1$：底物耗尽
$(1+r)\mu < 1$：酶全部失活
$(1+r)\mu = 1$：两者同时发生

奇异摄动分析结果

通过Burke等人（1990）的奇异摄动分析，得到两种情况：

Case 1（β < 1）：当$\beta = \psi/(\phi\rho) < 1$时，$e_i \to 1$（酶失活）

Case 2（β > 1）：当$\beta > 1$时，$s \to 0$（底物耗尽）

其中$\sigma, \rho, \psi, \phi$是与速率常数相关的无量纲参数。

6.5 合作现象（Cooperative Phenomena）

多结合位点酶

许多酶有多个底物结合位点。例如，血红蛋白有4个氧分子结合位点。当酶的一个位点结合底物后能影响其他位点的活性时，称为别构效应（Allosteric Effect）。若底物增加其他位点的活性则为激活剂，减少则为抑制剂。

双结合位点模型

考虑酶有两个结合位点的情况，反应机理为：

\[ S + E \underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}} C_1 \xrightarrow{k_2} E + P $$ $$ S + C_1 \underset{k_{-3}}{\overset{k_3}{\rightleftharpoons}} C_2 \xrightarrow{k_4} C_1 + P \]

反应速率

对于$e_0/s_0 \ll 1$的情况，Michaelis-Menten反应速率为：

\[ R_0(s_0) = e_0 \cdot \frac{k_2 K_m' + k_4 s_0}{K_m K_m' + K_m' s_0 + s_0^2} \]

其中： $$ K_m = \frac{k_2 + k_{-1}}{k_1},\quad K_m' = \frac{k_4 + k_{-3}}{k_3} $$

Hill方程

合作现象常用Hill方程描述： $$ R_0(S_0) = \frac{Qs_0^n}{K_m + s_0^n} $$

其中$n > 0$为Hill系数： - $n < 1$：负合作性 - $n = 1$：无合作性（Michaelis-Menten） - $n > 1$：正合作性

血红蛋白与氧气的结合呈典型的S形（sigmoidal）曲线，表现为正合作性；而肌红蛋白只有单个结合位点，其摄取曲线为典型的双曲线（Michaelis-Menten形式）。

6.6 自催化、激活与抑制（Autocatalysis, Activation and Inhibition）

自催化反应

自催化是化学物质参与自身生产的过程。简单例子：

\[ A + X \underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}} 2X \]

设A维持恒定浓度$a$，则： $$ \frac{dx}{dt} = k_1 ax - k_{-1}x^2 \Rightarrow x(t) \to x_S = \frac{k_1 a}{k_{-1}} $$

Lotka反应机制

Lotka（1920）提出的假设性振荡器模型： $$ A + X \xrightarrow{k_1} 2X,\quad X + Y \xrightarrow{k_2} 2Y,\quad Y \xrightarrow{k_3} B $$

该系统可转化为Lotka-Volterra方程，其解为时间周期解，但结构不稳定。

Thomas机制

Thomas（1975）机制涉及氧气和尿酸在尿酸酶存在下的反应，其无量纲形式为：

\[ \frac{du}{dt} = a - u - \rho R(u,v) = f(u,v) $$ $$ \frac{dv}{dt} = \alpha(b - v) - \rho R(u,v) = g(u,v) $$ $$ R(u,v) = \frac{uv}{1 + u + Ku^2} \]

该机制表现出底物抑制（Substrate Inhibition）特性。

Gierer-Meinhardt激活剂-抑制剂模型

\[ \frac{du}{dt} = a - bu + \frac{u^2}{v(1 + Ku^2)} = f(u,v) $$ $$ \frac{dv}{dt} = u^2 - v = g(u,v) \]

其中$u$为激活剂（自催化），$v$为抑制剂，且抑制剂抑制激活剂的生产。

一般性质

对于系统$\frac{du}{dt} = f(u,v)$，$\frac{dv}{dt} = g(u,v)$： - 若$\partial g/\partial u > 0$，则$u$激活$v$ - 若$\partial f/\partial v < 0$，则$v$抑制$u$

该系统可出现多稳态、滞后现象和阈值效应。

6.7 多稳态、蘑菇与隔离（Mushrooms and Isolas）

多稳态现象

在某些动力学系统中，参数变化可导致一个、两个或三个正稳态。稳态数目由nullclines的交点数目决定。

滞后回路（Hysteresis Loops）

当存在三个稳态时，其中中间分支不稳定。参数$p$增加和减少时，系统沿不同路径变化，形成滞后回路。

蘑菇结构（Mushrooms）

稳态$u_s$作为参数$p$的函数可呈现"蘑菇"形状，存在两个多稳态区域，每个区域都有滞后现象。

隔离结构（Isolas）

隔离是孤立的闭合解分支。当参数$p$从很小逐渐增加时，$u_s$沿一条路径变化；当$p$从很大逐渐减小时，$u_s$沿另一条不同路径变化。在常规增减过程中不会跨越隔离分支，需要有限扰动才能从一个分支跳到另一个分支。

Gray-Scott模型

Gray和Scott（1983, 1986）提出的自催化模型可产生蘑菇和隔离现象。该模型描述连续搅拌槽式反应器（CSTR）中的反应：

\[ \frac{dx}{dt} = k_0(x_0 - x) - k_1 xy^2 - k_3 x $$ $$ \frac{dy}{dt} = k_0(y_0 - y) + k_1 xy^2 + k_3 x - k_2 y \]

碘酸-亚砷酸反应

Ganapathisubramanian和Showalter（1984）的实验证明，无机碘酸-亚砷酸反应在适当条件下确实表现出蘑菇和隔离等多稳态行为。实验结果与理论模型预测高度吻合，证明这类复杂动力学行为确实存在于真实化学反应中。

公式汇总表

公式编号	名称	公式
(6.1)	Michaelis-Menten机制	$S + E \rightleftharpoons SE \rightarrow P + E$
(6.9)	Michaelis常数	$K_m = \frac{k_{-1} + k_2}{k_1}$
(6.10)	准稳态底物消耗	$\frac{ds}{dt} = -\frac{k_2 e_0 s}{s + K_m}$
(6.41)	Michaelis-Menten速率	$R_0 = \frac{Qs_0}{s_0 + K_m}$，其中$Q = k_2 e_0$
(6.42)	Michaelis-Menten方程	$\frac{ds}{dt} = -\frac{Qs}{K_m + s}$
(6.15)	快速瞬态时间	$t_c = \frac{1}{k_1(s_0 + K_m)}$
(6.16)	慢速时间尺度	$t_s \approx \frac{s_0 + K_m}{k_2 s_0}$
(6.18)	准稳态有效条件	$\varepsilon = \frac{e_0}{s_0 + K_m} \ll 1$
(6.110)	Hill方程	$R_0(S_0) = \frac{Qs_0^n}{K_m + s_0^n}$
(6.117)	Thomas机制	$\frac{du}{dt} = a - u - \rho\frac{uv}{1+u+Ku^2}$
(6.118)	Gierer-Meinhardt模型	$\frac{du}{dt} = a - bu + \frac{u^2}{v(1+Ku^2)}$
(6.123)	Gray-Scott模型	$\frac{dx}{dt} = k_0(x_0-x) - k_1xy^2 - k_3x$

关键概念回顾

准稳态近似（Quasi-Steady State）：当酶浓度远小于底物浓度时，复合物浓度变化率可近似为零。
奇异摄动分析（Singular Perturbation）：处理含有小参数的微分方程系统，通过内部解与外部解的匹配得到一致有效渐近解。
自催化（Autocatalysis）：产物参与自身的生产过程，是生物振荡器和开关行为的基础。
合作性（Cooperativity）：多结合位点酶显示出的协同效应，导致Hill系数$n \neq 1$。
多稳态与滞后：某些反应系统可具有多个稳定稳态，参数变化导致系统在不同状态间跳跃。
米氏常数$K_m$：衡量酶与底物亲和力的参数，$K_m = (k_{-1} + k_2)/k_1$。

本笔记基于J.D. Murray《Mathematical Biology I》第六章内容整理

公式编号	名称	公式
(6.1)	Michaelis-Menten机制	\(S + E \rightleftharpoons SE \rightarrow P + E\)
(6.9)	Michaelis常数	\(K_m = \frac{k_{-1} + k_2}{k_1}\)
(6.10)	准稳态底物消耗	\(\frac{ds}{dt} = -\frac{k_2 e_0 s}{s + K_m}\)
(6.41)	Michaelis-Menten速率	\(R_0 = \frac{Qs_0}{s_0 + K_m}\)，其中\(Q = k_2 e_0\)
(6.42)	Michaelis-Menten方程	\(\frac{ds}{dt} = -\frac{Qs}{K_m + s}\)
(6.15)	快速瞬态时间	\(t_c = \frac{1}{k_1(s_0 + K_m)}\)
(6.16)	慢速时间尺度	\(t_s \approx \frac{s_0 + K_m}{k_2 s_0}\)
(6.18)	准稳态有效条件	\(\varepsilon = \frac{e_0}{s_0 + K_m} \ll 1\)
(6.110)	Hill方程	\(R_0(S_0) = \frac{Qs_0^n}{K_m + s_0^n}\)
(6.117)	Thomas机制	\(\frac{du}{dt} = a - u - \rho\frac{uv}{1+u+Ku^2}\)
(6.118)	Gierer-Meinhardt模型	\(\frac{du}{dt} = a - bu + \frac{u^2}{v(1+Ku^2)}\)
(6.123)	Gray-Scott模型	\(\frac{dx}{dt} = k_0(x_0-x) - k_1xy^2 - k_3x\)