← 返回列表

Regio- and Stereoselectivity of CYP450BM3-Catalyzed Hydroxylation of Complex Terpenoids: A QM/MM Study

作者
Chenggong Hui, Warispreet Singh, Derek Quinn, Chun Li, Thomas S. Moody, Meilan Huang
单位
Department of Chemistry & Chemical Engineering, Queen's University, Belfast, BT9 5AG Northern Ireland, United Kingdom; Almac Sciences, Department of Biocatalysis and Isotope Chemistry, Almac House, 20 Seagoe Industrial Estate, Craigavon BT63 5QD, Northern Ireland, United Kingdom; Institute for Synthetic Biosystem/Department of Biochemical Engineering, School of Chemistry and Chemical Engineering, Beijing Institute of Technology, Beijing 100081, P.R. China; Key Lab of Industrial Biocatalysis Ministry of Education, Institute of Biochemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, P.R. China; Arran Chemical Company Limited, Unit 1 Monksland Industrial Estate, Athlone, Co. Roscommon, Ireland
杂志
Physical Chemistry Chemical Physics(2020)
DOI
10.1039/D0CP03083J
所属领域
合成生物学
阅读原文
查看 PDF 原文

解决的核心问题

P450BM3 对复杂萜类底物(如青蒿素)羟基化的区域选择性和立体选择性机制尚不明确,限制了其作为生物催化剂用于萜类生物合成的理性工程改造。

研究策略

结合分子对接、分子动力学模拟和量子力学/分子力学(QM/MM)计算,以青蒿素为模型底物,研究 P450BM3 变体催化羟基化的选择性和反应性。

研究路线图

100%
研究问题:P450BM3对青蒿素羟基化的区域/立体选择性机制不明
策略:分子对接+MD模拟+QM/MM计算
关键改造1:β1折叠片与β'螺旋构象变化调控底物定位
关键改造2:FL#62模板结合能差异(pro-C7S优于R)
关键改造3:X-E12催化C6a羟化势垒22.7 kcal/mol
关键改造4:IV-H4绝对S选择性(势垒18.7 vs 高5-10)
关键改造5:II-H10 Phe87-in/out构象决定R选择性
结果:活化势垒与O-H反应坐标正相关,作为活性预测描述符
结论:β'螺旋重塑口袋决定选择性,势垒-O-H距离线性关系指导P450理性设计
Text is not SVG - cannot display

下载 .drawio 源文件

核心内容

该研究通过分子对接、分子动力学模拟和QM/MM计算,系统解析了P450BM3变体催化青蒿素羟基化的区域与立体选择性机制。以FL#62、X-E12、IV-H4和II-H10四种变体为对象,计算结果显示FL#62对pro-C7(S)结合能(-5.31 kcal/mol)优于pro-C7(R)(-4.55 kcal/mol),与实验上C7(S)产物为主一致;IV-H4催化C7(S)羟化最低活化势垒为18.7 kcal/mol,而生成C7(R)的势垒高出5-10 kcal/mol,解释了绝对S选择性。研究还发现底物入口处β1折叠片构象变化及β'螺旋位移对底物定位至关重要,其中II-H10的Phe87-out构象平均活化势垒约18.9 kcal/mol,显著低于Phe87-in的约25.0 kcal/mol。关键的是,活化势垒与QM/MM优化反应物中O-H距离呈线性正相关,该描述符可用于便捷预测P450BM3反应活性。这项工作首次报道了P450对非天然环状萜类的QM/MM机理研究,为P450选择性羟化复杂环状分子的理性设计提供了理论指导。

创新点

1. 首次报道 P450 对非天然环状萜类底物羟基化的 QM/MM 机理研究;2. 揭示底物入口处 β1 折叠片的构象变化和 β' 螺旋位移对底物进入和定位的关键作用;3. 发现活化势垒与反应坐标(O-H 距离)之间的线性相关性,可作为预测 P450BM3 反应活性的简便描述符。

研究对象(4)

底盘细胞代谢通路基因蛋白功能说明
Bacillus megaterium 青蒿素羟化(萜类代谢) CYP102A1 P450BM3(FL#62 模板) 催化 模板变体,催化青蒿素 C7 位羟化,以 C7(S) 产物为主
Bacillus megaterium 青蒿素羟化(萜类代谢) CYP102A1 P450BM3 X-E12 催化 C6a 选择性 94%,QM/MM 平均活化势垒 22.7 kcal/mol
Bacillus megaterium 青蒿素羟化(萜类代谢) CYP102A1 P450BM3 IV-H4 催化 C7 位绝对 S 选择性,最低活化势垒 18.7 kcal/mol
Bacillus megaterium 青蒿素羟化(萜类代谢) CYP102A1 P450BM3 II-H10 催化 C7 位 100% R 选择性,Phe87-out 构象活化势垒约 18.9 kcal/mol,Phe87-in 约 25.0 kcal/mol

关键发现

  1. FL#62 对青蒿素的对接表明 pro-C7(S) 结合能 -5.31 kcal/mol,优于 pro-C7(R) 的 -4.55 kcal/mol,与实验上 C7(S) 产物为主一致;X-E12 催化 C6a 羟化的氢提取平均活化势垒为
  2. 7 kcal/mol(最低
  3. 1 kcal/mol),def2-SVP 与 def2-TZVP 基组势垒差小于 2 kcal/mol,二重态和四重态势垒分别为
  4. 7 和
  5. 8 kcal/mol,高度简并;IV-H4 催化 C7(S) 羟化的最低活化势垒为
  6. 7 kcal/mol,生成 C7(R) 的势垒高 5-10 kcal/mol,与绝对 S 选择性一致;II-H10 催化 C7(R) 羟化时,Phe87-out 构象平均活化势垒约
  7. 9 kcal/mol,Phe87-in 构象约
  8. 0 kcal/mol,Phe87-in 与 out 的能量差为
  9. 7 kcal/mol;活化势垒与 QM/MM 优化反应物中 O-H 距离呈正相关,OH 距离最短的代表结构对应最低势垒。

研究结论

通过 MD 和 QM/MM 研究发现,P450BM3 变体对青蒿素羟化的选择性由 β' 螺旋及其相邻 loop 和 F 螺旋重塑底物结合口袋决定;活化势垒与 O-H 反应坐标线性相关,可作为预测酶反应活性的便捷指标,为 P450 选择性羟化复杂环状分子的理性设计提供指导。