来源：Semantic Scholar | 编译：DNA Lab Space
导语
酶设计被誉为合成生物学的“圣杯”，但真正能将非天然反应催化得像天然酶一样高效的案例依然凤毛麟角。Morita–Baylis–Hillman（MBH）反应作为构建碳-碳键的有力工具，其人工酶的设计尝试一再受挫。近期，一项基于过渡路径抽样的分子动力学研究深入原子尺度，揭示出人工MBH酶催化效率低下的根本原因——不仅在于活性位点的静态几何匹配，更关键的是酶骨架动态耦合网络的缺失。这一发现为未来理性酶设计提供了全新视角，指出了从“结构模拟”走向“动力学模拟”的必要路径。
研究背景
酶催化以其惊人的速率加速和专一性，一直是合成生物学和生物催化领域追求的黄金标准。然而，自然界的催化工具箱并未覆盖所有化学家渴望的反应类型，Morita–Baylis–Hillman反应便是其中之一。MBH反应通过偶联α,β-不饱和羰基化合物与亲电试剂，生成高功能化的烯丙醇产物，在精细化学品和药物中间体合成中占据独特地位。由于自然界尚未发现能高效催化该反应的天然酶，科学家转而通过计算设计人工酶来填补这一空白。
过去二十年间，计算酶设计取得了长足进步，成功案例包括Kemp消除酶、Diels–Alder反应酶等。然而，针对MBH反应的多轮设计虽能获得具有微弱催化活性的蛋白骨架，催化效率普遍仅相当于背景反应的数十倍，远不及天然酶动辄10⁶–10⁸倍的速率提升。这一现象引发了一个根本性的追问：为何MBH酶的设计如此艰难？传统观点多聚焦于活性位点残基的几何排列和过渡态稳定化能，却忽视了蛋白质是处于永恒热运动中的动态实体。越来越多的证据表明，酶的催化能力不仅取决于静态活性口袋形状，更依赖于整个蛋白质骨架的构象涨落和动力学耦合网络。因此，要破解MBH酶设计困境，必须从蛋白质动力学维度重新审视催化循环中的决速步骤。
研究方法
为深入探测人工MBH酶的催化动态过程，研究团队采用了过渡路径抽样这一稀有事件模拟技术。与常规分子动力学模拟只能观察稳定态附近涨落不同，过渡路径抽样专门“捕捞”反应物穿越过渡态区域的成功轨迹，能够在不预先假设反应坐标的前提下，揭示真实的反应机理和决速步。
研究以先前通过Rosetta等计算平台设计获得的若干MBH人工酶为模型体系，这些酶骨架分别源自不同天然蛋白，如酮类固醇异构酶、丙糖磷酸异构酶等，活性位点引入关键催化残基（半胱氨酸、组氨酸、天冬氨酸等）以模拟小分子催化的酸碱双功能机制。作者首先为每个设计变体构建了酶-底物复合物的全原子模型，包含显式水分子和生理盐浓度离子。力场参数针对底物的非标准原子类型进行专门拟合，确保过渡态附近区域能量面的可靠性。
过渡路径抽样的实施分为多步：先通过伞形采样或约束动力学获得初始反应通道，再采用射箭算法从近过渡态区域发射大量无偏轨迹，收集反应性路径的系综。为确保统计收敛，每个体系产出了上百条独立过渡路径。分析阶段则不仅提取平均力势和自由能垒，更重点计算了反应坐标与蛋白质骨架运动的交叉相关函数、动力学群落分析以及过渡路径时间分布。这一套分析框架巧妙地将酶催化问题从“热力学稳定化”拓展至“动力学协同”层面。
研究结果
模拟结果揭示出令人深思的图景。在所测试的多种MBH设计酶中，即便活性位点残基的几何排列与小分子催化模型高度相似，其反应自由能垒仍比理论预期高出5–12 kcal/mol。更关键的是，通过对过渡路径的分解，研究者发现催化效率低下的根源并非热力学因素——底物结合和过渡态结合能基本合理——而是动力学因素：成功穿越过渡态的路径极度稀少，且路径所需的蛋白质骨架重排时间远长于天然酶的微秒级构象交换。
具体而言，过渡路径抽样显示，在活性较高的设计变体中，反应性路径仅占全部尝试的不足0.1%。这与天然酶通常近乎100%的反应性路径形成鲜明对比。进一步的分析发现，人工酶活性位点附近约5–8 Å范围内的残基构象涨落与反应坐标之间缺乏紧密耦合。在天然酶中，底物进入活性位点时引发的微小构象变化会像多米诺骨牌一样通过氢键网络和疏水堆积传递至全酶，从而有效降低活化熵；而设计酶的这些通讯通路是断裂的，导致局部重构难以被“感知”和“响应”，使得大量动能耗散在非生产性运动中。
一项引人注目的量化指标是反应协同因子：天然MBH反应的碳-碳键形成与质子转移具有高度协同性，但人工酶中这两个化学步骤的去耦合概率显著上升，瞬态中间体寿命延长，为副反应和溶剂猝灭提供了窗口。此外，水分子在活性位点的滞留行为也耐人寻味——人工酶的疏水口袋构型不够精确，无法像天然酶那样在过渡态时刻彻底排尽活性位点水分子，导致过渡态极性环境被扰动，进一步抬高了有效活化能。
研究者还通过群落分析将蛋白质划分为若干动力学单元，发现设计酶中连接活性位点与远距离骨架的“变构高速公路”未形成。即便有少量反应性轨迹成功抵达产物态，轨迹回溯显示这些稀有事件往往对应着非典型的、瞬时的残基偶联模式，类似于“阵发性”的网络接通。这提示设计酶的瓶颈不在于能否实现一次催化转化，而在于其无法像天然酶那样稳定地、可重复地维持催化活性构象的系综。
讨论与展望
这项研究从根本上重新定义了MBH酶设计的困难：它不仅是“钥匙与锁”的形状适配问题，更是一道构建分子内动力学通讯网络的难题。如果蛋白质骨架无法在皮秒至纳秒的时间尺度上与化学反应坐标起舞，那么即便活性位点照片与设计模型完美重合，催化效率仍将远低于预期。这为“催化抗体”与“天然酶”在效率上的差距提供了一种机制性解释——催化抗体往往仅优化了静态结合，而缺乏进化的动力学优选过程。
从方法学视角，过渡路径抽样的引入为合成生物学提供了超越自由能微扰的全新诊断工具。未来人工酶设计不应止步于过渡态静态类似物的几何优化，而必须将动力学匹配纳入打分函数，例如通过设计阶段引入弹性网络模型筛选动力学耦合强度，或在分子动力学后处理中评估变构通路完整性。此外，定向进化不仅是热力学结合力的筛选，其实质很可能正是逐步修复动力学缺陷、重建骨架通讯网络的随机行走过程。将过渡路径分析的前沿计算与实验定向进化深度整合，有望最终打破MBH反应乃至更广泛非天然反应酶设计的困局。
参考来源
- Balasubramani S G, Korchagina K, Schwartz S D. Transition Path Sampling Study of Engineered Enzymes That Catalyze the Morita–Baylis–Hillman Reaction: Why Is Enzyme Design so Difficult? *Journal of Chemical Information and Modeling*, 2024. DOI: 10.1021/acs.jcim.4c00045
DOI: 10.1021/acs.jcim.4c00045