《Biochemical Journal》：Docking in the
Abies grandis abietaenol synthase to investigate the mechanism underlying conifer resin acid diterpene synthase evolution
编辑推荐：
本研究聚焦松科植物二萜树脂酸生物合成，针对其关键酶——密切相关的(类别I)二萜合酶的功能分化机制尚不明确这一问题，通过分析巨冷杉（Abies grandis）的冷杉萜醇合酶（AgAS），探究其向(异)海松二烯合酶演化的关键残基。研究人员发现，将AgAS第723位丙氨酸（A）替换为苏氨酸（T）的A723T突变体能几乎完全将其产物由冷杉萜醇转变为异海松-7,15-二烯。通过应用TerDockin计算方法，揭示该苏氨酸残基很可能直接作为催化碱基，通过去质子化初始形成的异海松-13E-烯-8-基碳正离子中间体，从而“短路”了天然反应路径。这为理解二萜合酶催化的复杂碳正离子级联反应机制提供了新见解，并对更广泛的萜类合酶反应建模具有重要意义。
松科植物（如冷杉、松树、云杉）是地球上重要的针叶树种，它们分泌的油性树脂是天然的防御屏障。这些树脂中的关键成分是一类被称为二萜树脂酸的化合物，它们具有(异)海松烷型或重排的冷杉烷型三环骨架。在自然界中，这些骨架各异但化学结构相近的化合物，却是由一群关系很近的二萜合酶“催化工厂”生产出来的。这些酶如何通过微小的改变，就实现了从生产A化合物到生产B化合物的“业务转型”？这背后复杂的碳正离子级联反应机制是怎样的？理解这一点，不仅有助于我们窥探植物次生代谢进化的奥秘，也对酶工程和合成生物学具有重要价值。
长期以来，来自巨冷杉（Abies grandis）的冷杉萜醇合酶（AgAS）一直是研究这类酶的模型。它原本的工作是将底物香叶基香叶基焦磷酸（GGPP）经过一系列复杂的环化、重排和质子转移，最终主要转化为带有羟基的冷杉萜醇。然而，在松科植物的进化历程中，从类似AgAS的“祖先”酶独立演化出了至少三次，产生出了主要制造异海松二烯这类不含羟基的烯烃产物的新酶。这种功能的“分道扬镳”究竟是如何在分子水平上实现的？
为了解答这个问题，发表在《Biochemical Journal》上的这项研究，以AgAS为研究对象，展开了一场从序列比对到计算模拟的深入探索。研究人员首先通过序列比对，在与AgAS亲缘关系更近的香脂冷杉（Abies balsamea）的AS和PS配对中，发现了一个关键差异：在PS中，对应AgAS第723位的丙氨酸（A）被替换成了苏氨酸（T）。受此启发，他们在AgAS中构建了A723T等一系列定点突变体。实验结果令人惊讶：A723T突变体几乎完全（97%）将产物从冷杉萜醇转换成了异海松-7,15-二烯，效果比之前报道的A723S突变（产生58%的异海松-7,15-二烯）更为彻底。这一发现强烈暗示，第723位引入的羟基（来自丝氨酸S或苏氨酸T）直接参与了催化。
为了深入探究其背后的机理，研究人员采用了名为TerDockin的计算方法。这是一种化学信息引导的、基于分子力学的对接方法，能够将反应中间体准确地对接进酶的活性位点。他们首先对AgAS:A723T突变体进行了模拟，尝试解释苏氨酸如何促使反应“短路”。然而，最初的模拟结果与实际观察不符，预测产物主要是另一个异构体。研究团队意识到，他们忽略了一个关键“演员”：在天然反应中，一个水分子会加成到最终的碳正离子上形成冷杉萜醇的羟基。当他们在模拟中引入并合理约束这个“反应水”的位置后，计算结果发生了根本性改变，成功预测了异海松-7,15-二烯作为主要产物。模拟结构显示，A723T突变引入的苏氨酸羟基，在G1螺旋偶极矩的激活下，其位置恰好可以直接作为催化碱基，去夺取异海松-15-烯-8-基碳正离子中间体（记作A）C7位上的一个质子，从而直接生成异海松-7,15-二烯，阻止了反应继续进行后续的重排和水合步骤。进一步对野生型AgAS的模拟则显示，在没有突变的情况下，碳正离子中间体A会经过一系列复杂的重排（包括环翻转、乙烯基旋转、分子内1,4-质子转移和1,2-甲基迁移）转化为另一个碳正离子中间体B，随后被附近的水分子进攻，最终生成冷杉萜醇。
本研究主要应用了以下关键技术方法：1) 蛋白序列比对与系统发育分析：用于识别不同松科植物二萜合酶中VSIAL基序的保守性与变异，指导关键残基的选取。2) 定点诱变与体外产物分析：在大肠杆菌（Escherichia coli）中共表达GGPP合酶与AgAS突变体，通过有机溶剂萃取和GC-MS（气相色谱-质谱联用）分析产物组成变化。3) 核磁共振（NMR）分析：用于精确鉴定野生型AgAS主要产物的化学结构（包括C13差向异构体）。4) TerDockin计算建模：这是一种专门用于研究萜类合酶反应机理的计算方法，通过对接优化的反应中间体构象库并施加化学约束，来模拟活性位点内的反应过程。
研究结果部分：
关键基序的保守性模式识别
通过比对松科植物TPS-d3亚家族的酶序列，研究人员确认了VSIAL基序在产冷杉烷型二萜的合酶中高度保守。然而，在产(异)海松二烯的合酶中，该基序发生了变异。特别是在云杉属（Picea）中，丙氨酸（A）被丝氨酸（S）取代；而在香脂冷杉（Abies balsamea）的PS中，则对应为苏氨酸（T）取代。这些发现为后续的突变实验提供了进化依据。
产物分析
对野生型AgAS的产物进行GC-MS和NMR分析证实，其主要产物是一对差向异构的冷杉萜醇（6a和6b），其中13α-羟基异构体（6a）占约70%。这表明AgAS本质上是一个冷杉萜醇合酶。
突变分析
在AgAS中引入基于序列比对发现的突变。关键发现包括：A723S突变体主要产生异海松-7,15-二烯（7, 58%）；而A723T突变体几乎完全（97%）转向生产7。I722T和L724G突变也显著改变了产物谱，但效果不同，这强调了精确的底物定位对产物特异性的关键作用。
TerDockin分析
计算模拟表明，在AgAS:A723T突变体中，只有当模拟中包含了用于生成冷杉萜醇的“反应水”分子并加以合理约束时，预测的产物分布才与实验观测（主要生成7）相符。模拟结构支持苏氨酸723的侧链羟基直接作为碱基，去质子化碳正离子中间体A的C7位，从而生成7。对野生型AgAS的模拟则揭示了从中间体A到B的复杂重排路径，并表明水分子在反应后期才参与，其位置可能靠近焦磷酸盐副产物。
结论与讨论：
本研究强调了VSIAL基序在决定松科植物二萜合酶产物特异性中的核心作用。特别值得注意的是，将第723位丙氨酸替换为苏氨酸（A723T）能近乎完美地将AgAS从冷杉萜醇合酶转变为异海松二烯合酶。TerDockin计算模拟为此提供了强有力的机理解释：引入的苏氨酸羟基很可能直接充当催化碱基，通过去质子化关键的碳正离子中间体A，从而“短路”了原本漫长的天然反应途径。丝氨酸（S）和苏氨酸（T）仅一个甲基之差，却对催化效率产生显著影响，这凸显了萜类合酶活性位点内底物定位的极端精确性。
这项工作的意义在于双方面：首先，它深化了我们对TPS-d3亚家族二萜合酶活性位点相互作用及其进化路径的理解，揭示了自然界如何通过单个氨基酸的替换实现酶功能的重大转变。其次，它验证并拓展了TerDockin这一计算工具的应用。研究表明，要准确模拟此类复杂反应，必须考虑所有反应组分（如本研究中的水分子）的约束。尽管该方法在完全复现实验观测的高特异性方面仍有改进空间，但其相对低廉的计算成本允许进行迭代式探索，为理解其他萜类合酶的催化机制提供了有价值的范式和见解。总之，这项研究在酶学机制与计算生物学交叉领域做出了重要贡献。
