Enzymatic hydroxylation of an unactivated methylene C-H bond guided by molecular dynamics simulations
- 作者
- Alison R. H. Narayan, Gonzalo Jiménez-Osés, Peng Liu, Solymar Negretti, Wanxiang Zhao, Michael M. Gilbert, Raghunath O. Ramabhadran, Yun-Fang Yang, Lawrence R. Furan, Zhe Li, Larissa M. Podust, John Montgomery, K. N. Houk, David H. Sherman
- 单位
- University of Michigan, Ann Arbor, MI, USA; University of California, Los Angeles, CA, USA; University of California, San Francisco, CA, USA
- 杂志
- Nature Chemistry（2015）
- DOI
- 10.1038/nchem.2285
- 所属领域
- 合成生物学
- 关键词
- 阅读原文
- 查看 PDF 原文
解决的核心问题
实现未活化亚甲基C-H键的高位点和立体选择性氧化，特别是薄荷醇衍生物中远离导向基团的亚甲基，避免繁琐的蛋白质工程改造。
研究策略
利用分子动力学(MD)模拟和量子力学(QM)计算，结合PikC单加氧酶天然的盐桥锚定机制，通过底物工程（引入合成锚定基团）和蛋白质工程（理性突变）扩大底物范围并调控选择性。
核心内容
该研究以分子动力学模拟和量子力学计算为指导，将PikC单加氧酶从天然大环内酯底物扩展至薄荷醇等未活化亚甲基C-H键的氧化。通过底物工程引入合成锚定基团，并利用天然盐桥锚定机制，结合D50N/D176Q/E246A三突变消除不良相互作用，显著提升催化效率。计算预测与实验高度吻合，例如三突变体使(-)-28的总周转数从10提升至170，(-)-23从无产物提升至38。对(-)-薄荷醇，酶活性位点倾向氧化C4-H_eq位，生成单一反式二醇产物；对(+)-薄荷醇则检测到至少5种异构体。该平台成功氧化(+)-α-松油醇、(-)-异胡薄荷醇及稠环/桥环双环底物，均实现远端抗式羟基化且烯丙位C-H未被氧化。QM预测非酶模型轴向C4-H抽象能垒低0.7 kcal/mol，但MD和实验证明酶活性位点偏好赤道位，体现了酶环境对选择性的决定性影响。该计算驱动策略为C-H功能化生物催化提供了无需繁琐筛选的高效路径。
创新点
建立计算预测模型用于底物范围、位点选择性和立体选择性预测；将PikC底物从天然大环内酯扩展到薄荷醇等六元环和双环骨架；通过D50N/D176Q/E246A三突变消除不良相互作用，显著提高催化效率。
研究对象（5）
| 底盘细胞 | 代谢通路 | 基因 | 蛋白 | 功能 | 说明 |
|---|---|---|---|---|---|
| Streptomyces venezuelae ATCC 15439 | pikromycin生物合成通路 | pikC | PikC | 催化 | 细胞色素P450单加氧酶，催化底物羟基化，天然底物为YC-17和纳尔波霉素 |
| Streptomyces venezuelae ATCC 15439 | pikromycin生物合成通路 | pikC (D50N) | PikC_D50N | 催化 | 单突变体，提高催化效率，实验中使用RhFRED融合蛋白形式 |
| Streptomyces venezuelae ATCC 15439 | pikromycin生物合成通路 | pikC (D50N/D176Q/E246A) | PikC_D50N_D176Q_E246A | 催化 | 三突变体，去除D176和E246的羧酸基团，消除不良盐桥，使活性位点保持闭合构象，显著提高TTN |
| 无（体外酶反应） | 电子传递链 | RhFRED | RhFRED还原酶域 | 能量供给 | 来自Rhodococcus的NADPH依赖性还原酶，与PikC融合形成自给自足酶，提供还原力 |
| 无（体外酶反应） | NADPH再生 | zwf (来自Leuconostoc mesenteroides) | 葡萄糖-6-磷酸脱氢酶 | 催化 | 用于NADPH再生系统，偶联葡萄糖-6-磷酸 |
关键发现
- (-)-24和(-)-25的TTN分别为24和43
- (-)-26的TTN为96，优于间位(41)和邻位(10)类似物
- 三突变体PikC_D50N_D176Q_E246A使(-)-28的TTN从10提高至170，(-)-23从无产物提升至TTN 38。QM预测非酶模型轴向C4-H抽象能垒低0.7 kcal/mol（对应25°C下3:1动力学比），但MD模拟和实验证明酶活性位点倾向氧化(-)-薄荷醇的C4-H_eq，产物为diol 8。对(+)-薄荷醇底物，至少检测到5种异构体，TTN为2-28。该平台成功氧化(+)-α-松油醇、(-)-异胡薄荷醇及稠环/桥环双环底物，产物均为远端抗式羟基化，烯丙位C-H未被氧化。
研究结论
计算驱动的底物和蛋白质工程可以高效扩展P450单加氧酶的底物范围，实现对未活化亚甲基C-H键的位点和立体选择性羟基化，且预测与实验高度吻合，为C-H功能化生物催化提供新策略。