个人学习笔记

  • 首页
  • NONMEM
  • PBPK
  • R语言
  • 其他
  • 常用链接
  • 博客简介
  • 隐私政策
学海无涯勤可渡,书山万仞志能攀
  1. 首页
  2. 酶动力学
  3. 正文

Michaelis–Menten 方程怎么推导和应用:从八步公式到抑制剂与药物发现

2026年9月25日 14点热度 0人点赞 0条评论

摘要

本笔记根据 Easton Henderson 于 2026 年 1 月 12 日发表的英文综述《Mastering the Michaelis-Menten Equation》整理。文章从 E + S ⇌ ES → E + P 出发,串起公式推导、实验测定、非线性拟合、抑制剂分析和药物发现应用。下面按“理论 → 实验 → 拟合 → 机制解释 → 排错”的顺序改写,便于直接作为实验清单使用。

原文:Mastering the Michaelis-Menten Equation: A Comprehensive Guide to Derivation, Assumptions, and Modern Applications in Drug Discovery。

Michaelis-Menten 酶动力学示意图
图 1 原文配图,用于定位 E、S、ES 和 P 的关系。

一、模型从哪里开始

经典模型处理单底物、单产物反应:E + S ⇌ ES → E + P。k1 表示结合,k−1 表示解离,k2(即 kcat)表示催化。实验测量的是反应刚开始时的初始速率 v0,而不是整条进程曲线的平均速率。

二、八步推导 Michaelis–Menten 方程

  1. 定义产物速率: v = d[P]/dt = k2[ES]。
  2. 写出 ES 微分方程: d[ES]/dt = k1[E][S] − k−1[ES] − k2[ES]。
  3. 应用稳态近似:过渡期结束后令 d[ES]/dt ≈ 0,得到 k1[E][S] = (k−1+k2)[ES]。
  4. 应用酶守恒: [E]T = [E]+[ES],所以 [E]=[E]T−[ES]。
  5. 代入并解出 [ES]: [ES] = k1[E]T[S] / ((k−1+k2)+k1[S])。
  6. 定义 Km: Km=(k−1+k2)/k1,于是 [ES]=[E]T[S]/(Km+[S])。
  7. 定义 Vmax:酶饱和时 [ES]≈[E]T,所以 Vmax=k2[E]T。
  8. 得到最终式: v0=Vmax[S]/(Km+[S])。当 [S]=Km 时,速率等于 Vmax/2。

这条式子给出矩形双曲线:低底物区近似一级动力学,高底物区逐渐饱和。稳态推导不要求快速平衡;只有当 k2 远小于 k−1 时,Km 才可近似为结合解离常数 Kd。

三、参数怎么读

参数 关系式 解释 研发用途
Km (k−1+kcat)/k1 达到半最大速率的底物浓度,不总等于亲和力 设计底物梯度、解释竞争和 IC50
Vmax kcat[E]T 给定酶量下的最大速率 比较催化能力、检查酶活性
kcat k2 单个活性位点每秒的周转数 比较酶效率
kcat/Km 低底物区表观二级速率常数 综合表示捕获并转化底物的能力 比较底物选择性

四、实验前的假设检查

  • 稳态:只在线性初始区取速率;若出现快速 burst,应改用 pre-steady-state 分析。
  • 酶守恒:反应期间酶不明显失活、聚集或吸附。
  • 底物不耗竭:让 [S]0 远大于 [E]T,通常把转化率控制在 5% 以下;当 [E]T/[S]0<0.01 时耗竭通常小于 1%,0.01–0.05 时应考虑积分方程,超过 0.05 时应拟合完整时间过程。
  • 产物可忽略:初始区内 [P] 近似为零,逆反应和产物抑制尚未明显发生。
  • 单底物、无协同性:双底物、多个结合位点或变构调节要改用相应机制模型。

五、从配液到初始速率的操作步骤

  1. 准备能维持 pH、离子强度和酶稳定性的缓冲液及辅因子(如 Mg2+、ATP、NADH/NADPH)。
  2. 记录酶浓度、活性批次、底物纯度和稳定性。
  3. 以预计 Km 为中心设计约 0.2×–5× Km 的 6–8 个底物浓度,设置重复。
  4. 加入无酶/热灭活酶、无底物、阳性酶和参考抑制剂对照。
  5. 用多道移液器或自动化方式同步启动反应。
  6. 用吸光、荧光、放射性或 LC–MS 追踪产物,先确认检测信号在线性范围。
  7. 对每条产物-时间曲线只在线性早期做回归,斜率就是 v0。
  8. 检查不同酶量下速率是否成比例、曲线是否线性;若不满足,先排查失活、耗竭、混合延迟和背景。

六、Km 和 Vmax 的拟合

推荐直接对原始 v0–[S] 数据做非线性回归:v0=Vmax[S]/(Km+[S]),报告估计值、95% 置信区间、残差和重复实验一致性。GraphPad Prism、SigmaPlot、R、Python 或专用酶动力学软件都可以完成。不要把 Lineweaver–Burk 作为主要估计方法,因为倒数会放大低底物区误差;它更适合作为诊断图。

诊断式为 1/v=(Km/Vmax)(1/[S])+1/Vmax:y 截距是 1/Vmax,斜率是 Km/Vmax,x 截距是 −1/Km。

七、抑制剂实验如何从 IC50 走到 Ki

先做无抑制剂底物梯度,得到基线 Km 和 Vmax;再构建“多档底物 × 多档抑制剂”的二维矩阵,每个组合做重复,从每条曲线提取 v0,最后对全部数据做全局拟合。

类型 表观 Vmax 表观 Km 机制
竞争性 不变 升高 抑制剂与底物竞争游离酶
非竞争性 降低 理想情况下不变 可结合游离酶和 ES
反竞争性 降低 降低 优先结合 ES
混合型 降低 可升可降 对游离酶和 ES 的亲和力不同

IC50 是特定底物浓度下的表观指标,不能直接当作 Ki;只有在机制、底物浓度和模型都明确后,才可使用 Cheng–Prusoff 关系等校正。

八、出现偏离时怎样排错

  • 低底物区残差大:检查检测限、背景扣除和移液误差,必要时使用加权回归。
  • 高底物区速率下降:考虑底物抑制、溶解度、黏度或变构效应。
  • 曲线先快后慢:检查酶失活、产物抑制、底物耗竭和反应可逆性。
  • 不同酶量下速率不成比例:检查吸附、聚集、活性位点计数和底物耗竭。
  • 跨实验日 Km 漂移:核对温度、pH、辅因子新鲜度、底物批次和仪器线性范围。
  • 双底物体系:改用有序、随机或 ping-pong 机制模型,不要把表观 Km 当成单底物真值。

九、与药物发现的连接

经过机制验证的 Km、Vmax、kcat 和 Ki 可用于比较候选物的抑制强度、区分结合效应与催化效应,并进入 IVIVE、药物相互作用和系统生物学模型。解释这些参数时必须同时记录底物浓度、酶浓度、温度、pH、辅因子和检测方式,因为它们决定结果能否比较。

十、最后的实验检查清单

  1. 是否明确反应式、测量终点和初始速率区间?
  2. 是否证明底物耗竭、产物积累和酶失活可忽略?
  3. 底物浓度是否覆盖低、中、高三个区域?
  4. 是否用原始速率做非线性拟合并检查残差?
  5. 是否报告置信区间、重复性和对照?
  6. 抑制剂实验是否同时改变底物和抑制剂浓度并做全局拟合?
  7. 偏离双曲线时是否先验证底物抑制、协同性、多底物或可逆反应?

结语

Michaelis–Menten 方程的价值不只是记住一条公式,而是把“结合、催化、测量和机制判断”连成一条可检查的链条。最稳妥的顺序是先验证假设,再测初始速率;先用原始数据做非线性拟合,再用诊断图解释偏差;最后才把参数外推到药物发现和体内模型。

整理来源:Henderson E. Mastering the Michaelis-Menten Equation: A Comprehensive Guide to Derivation, Assumptions, and Modern Applications in Drug Discovery. Enzyme Explorer, 2026-01-12。

标签: 暂无
最后更新:2026年9月25日

MSxiaoming

得失从缘,心无增减

点赞

MSxiaoming

得失从缘,心无增减

标签聚合
NONMEN POPPK QSP 群体药代动力学 R语言 PBPK模型 R语言,数据管理 NONMEM

COPYRIGHT © 2023 个人学习笔记. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

皖ICP备2023004398号-1