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 变体催化羟基化的选择性和反应性。
核心内容
该研究通过分子对接、分子动力学模拟和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 |
关键发现
- FL#62 对青蒿素的对接表明 pro-C7(S) 结合能 -5.31 kcal/mol,优于 pro-C7(R) 的 -4.55 kcal/mol,与实验上 C7(S) 产物为主一致;X-E12 催化 C6a 羟化的氢提取平均活化势垒为
- 7 kcal/mol(最低
- 1 kcal/mol),def2-SVP 与 def2-TZVP 基组势垒差小于 2 kcal/mol,二重态和四重态势垒分别为
- 7 和
- 8 kcal/mol,高度简并;IV-H4 催化 C7(S) 羟化的最低活化势垒为
- 7 kcal/mol,生成 C7(R) 的势垒高 5-10 kcal/mol,与绝对 S 选择性一致;II-H10 催化 C7(R) 羟化时,Phe87-out 构象平均活化势垒约
- 9 kcal/mol,Phe87-in 构象约
- 0 kcal/mol,Phe87-in 与 out 的能量差为
- 7 kcal/mol;活化势垒与 QM/MM 优化反应物中 O-H 距离呈正相关,OH 距离最短的代表结构对应最低势垒。
研究结论
通过 MD 和 QM/MM 研究发现,P450BM3 变体对青蒿素羟化的选择性由 β' 螺旋及其相邻 loop 和 F 螺旋重塑底物结合口袋决定;活化势垒与 O-H 反应坐标线性相关,可作为预测酶反应活性的便捷指标,为 P450 选择性羟化复杂环状分子的理性设计提供指导。