Catalytic Cycle of Haloalkane Dehalogenases towards Unnatural Substrates Explored by Computational Modeling
- 作者
- Sérgio M. Marques, Zuzana Dunajova, Zbynek Prokop, Radka Chaloupkova, Jan Brezovsky, Jiri Damborsky
- 单位
- Loschmidt Laboratories, Department of Experimental Biology and Research Centre for Toxic Compounds in the Environment RECETOX, Faculty of Science, Masaryk University; International Clinical Research Center, St. Anne's University Hospital Brno
- 杂志
- Journal of Chemical Information and Modeling（2017）
- DOI
- 10.1021/acs.jcim.7b00070
- 所属领域
- 酶工程
- 关键词
- 阅读原文
- 查看 PDF 原文
解决的核心问题
解析工程化脱卤酶DhaA31对非天然底物1,2,3-三氯丙烷(TCP)高活性的结构基础，阐明其催化循环中各个步骤（底物结合、反应复合物形成、化学步骤、产物释放）的分子机理及限速步骤的转变机制。
研究策略
采用分子动力学(MD)和加速分子动力学(aMD)模拟野生型DhaA和突变体DhaA31，结合分子对接、QM/MM绝热映射、MM/GBSA结合自由能计算等计算手段，系统解析TCP从结合到产物释放的完整催化循环；并构建中间突变体进行动力学实验验证理论预测。
核心内容
该研究通过分子动力学、加速分子动力学、QM/MM绝热映射及MM/GBSA结合自由能计算等手段，系统解构了工程化脱卤酶DhaA31对非天然底物1,2,3-三氯丙烷（TCP）的完整催化循环，涵盖底物结合、反应复合物形成、SN2化学步骤和产物释放各阶段。结果显示，DhaA31对TCP的kcat较野生型提高32倍，催化效率kcat/Km提高26倍，其中SN2化学步骤速率提升145倍。关键突变C176Y和V245F的贡献得到定量归因：C176Y主要增加SN2反应构象比例并降低能垒，但同时也使产物DCP的释放成为新的限速步骤；隧道动力学分析表明DhaA31的p1隧道开口率低于野生型，而活性位点水合程度无显著差异，推翻了此前相关假说。该工作首次完整解构了DhaA31的多步催化循环，明确了限速步骤从化学步骤向产物释放的转移机制，为平衡设计高效脱卤酶提供了分子水平指导，其多步骤计算方法亦可推广至其他酶体系。
创新点
首次对DhaA31的完整催化循环进行多步骤计算解构，推翻了此前关于活性位点水合减少的假说；明确C176Y和V245F是关键突变，前者主要提高SN2反应构象比例并降低能垒，但导致产物释放成为限速步骤；将计算模拟与突变体动力学分析相结合，实现了对催化机理的定量归因。
研究对象（6）
| 底盘细胞 | 代谢通路 | 基因 | 蛋白 | 功能 | 说明 |
|---|---|---|---|---|---|
| 无细胞体系（计算模拟） | 卤代烷烃降解途径 | dhaA | DhaA | 催化 | 野生型脱卤酶，来自Rhodococcus rhodochrous，对TCP活性低 |
| 无细胞体系（计算模拟） | 卤代烷烃降解途径 | dhaA31 | DhaA31 | 催化 | 含I135F、C176Y、V245F、L246I、Y273F五个突变，是目前对TCP活性最高的HLD变体 |
| 无细胞体系（计算模拟） | 卤代烷烃降解途径 | dhaA04 | DhaA04 (C176Y) | 催化 | 单突变体，用于实验验证C176Y对产物释放的影响 |
| 无细胞体系（计算模拟） | 卤代烷烃降解途径 | dhaA138 | DhaA138 (V245F) | 催化 | 单突变体，用于实验验证V245F对SN2步骤的影响 |
| 无细胞体系（计算模拟） | 卤代烷烃降解途径 | dhaA137 | DhaA137 (L246I) | 催化 | 单突变体，用于实验验证L246I对催化性能的影响 |
| 无细胞体系（计算模拟） | 卤代烷烃降解途径 | dhaA133 | DhaA133 (C176Y/V245F) | 催化 | 双突变体，用于验证关键突变的协同效应 |
关键发现
- DhaA31对TCP的kcat比野生型提高32倍(达32-fold)，催化效率kcat/Km提高26倍；SN2化学步骤的速率提高145倍；野生型DhaA的kcat=0.04 s^-1，kcat/Km=40 s^-1 M^-1；DhaA31中产物DCP的释放成为限速步骤，主要归因于C176Y突变；隧道动力学分析显示DhaA31的p1隧道开口率低于野生型（MD:
- 6%±1.0% vs
- 9%±1.7%；aMD:
- 4%±5.4% vs
- 6%±13.4%），但最大隧道半径相近（~3Å）；活性位点水合程度无显著差异。
研究结论
综合计算建模揭示了DhaA31催化性能增强的分子基础：SN2反应速率的提升源于反应构象比例增加和能垒降低，而C176Y突变同时导致产物释放变慢，使整个催化循环的限速步骤从化学步骤转移至产物释放。该研究为平衡设计更高效的脱卤酶提供了指导，所采用的多步骤催化循环计算方法可推广至其他酶系统。