《International Journal of Molecular Sciences》：Evolutionary Adaptability Coupled with Computation-Driven Engineering for Thermostability Enhancement of a Deoxynivalenol-Detoxifying Fusion Enzyme
编辑推荐：
脱氧雪腐镰刀菌烯醇（DON）是一种常见于谷物及其制品中的单端孢霉烯类真菌毒素，对人类和动物健康构成显著风险。前期工作中，研究人员构建了一种由脱氢酶DADH与醛酮还原酶AKR13B3组成的融合酶，可一步将DON转化为无毒的3-表-DON（3-epi-DON）。然
脱氧雪腐镰刀菌烯醇（DON）是一种常见于谷物及其制品中的单端孢霉烯类真菌毒素，对人类和动物健康构成显著风险。前期工作中，研究人员构建了一种由脱氢酶DADH与醛酮还原酶AKR13B3组成的融合酶，可一步将DON转化为无毒的3-表-DON（3-epi-DON）。然而该融合酶热稳定性差，限制了其工业应用。本研究利用EVcouplings与GRAPE-WEB平台鉴定调控AKR13B3–DADH融合酶热稳定性的关键氨基酸残基。通过单位点突变筛选与有益位点组合，获得三突变体M361L/T508Y/Y603F（M1）。M1在50 °C下的半衰期约500 min，较野生型提升16.4倍，催化活性提升2.7倍，表观熔解温度（Tm）升高约4 °C。分子动力学模拟证实热稳定性提升源于关键区域构象柔性降低、结构堆积增强及氢键网络强化。结果表明研究人员成功开发了热稳定DON解毒融合酶，为其工业应用提供了实践基础。
研究背景与立项依据：脱氧雪腐镰刀菌烯醇（DON）作为单端孢霉烯类真菌毒素广泛污染谷物及衍生品，传统理化脱毒法存在去除不彻底、产生有害副产物及损伤感官营养等缺陷。生物脱毒凭借底物专一性强、无二次污染和高催化效率成为优选，其中C3差向异构化路径（DON→3-酮-DON→3-表-DON）毒性递减。研究人员前期从土壤Devosia菌株A6-243鉴定出DADH与AKR13B3，并构建AKR13B3–DADH融合酶实现一步转化，但该酶50 °C半衰期仅约30 min，而饲料制粒工艺需耐受65–80 °C，热稳定性短板严重制约产业化。既有文献未见针对DON解毒融合酶的蛋白工程改造，故本研究引入计算驱动半理性设计予以突破。
主要关键技术方法：研究人员以实验室保存的AKR13B3–DADH融合酶基因为对象，先后采用EVcouplings对DADH（比对8998条同源序列）与AKR13B3（比对89276条同源序列）域分别做共进化耦合分析，以独立分数阈值1筛选非绝对保守且强耦合残基；继而将全序列导入GRAPE-WEB，联合FoldX（截值?1.5 kcal/mol）、Rosetta（截值?1.0 REU）、ABACUS（截值?2.5 AEU）三类能量函数算法打分，聚类后构建74个单点突变库。实验验证依托大肠杆菌BL21(DE3)/pET28a表达系统，经Ni亲和层析获纯酶，以50 °C孵育后残活测定初筛，继以50 °C半衰期（t1/2）复筛；热力学参数用差示扫描荧光法（DSF）测表观熔解温度（Tm）；机制解析采用AMBER 24全原子100 ns分子动力学模拟（RMSD、RMSF、Rg、SASA、氢键统计）及AlphaFold3建模与AutoDock Vina对接。
研究结果
2.1 融合酶热稳定关键残基识别：经EVcouplings共进化分析剔除完全保守位，获DADH域696个、AKR13B3域203个候选；GRAPE-WEB三算法分别产出271、265、212个候选，整合聚类得74突变体库，兼顾通量与命中率。
2.2 突变对催化与热稳影响：74单点50 °C处理后，6株残活显著优于野生型，9株催化升而不损热稳，另5株（S174L等）失活揭示其为活性构象关键。二级筛选15株中M361L单点t1/2达97.6 min（3.2倍），T508Y催化达287.2%。六株优株（A160V、M361L、A387T、S418P、T508Y、Y603F）两两组合得15双点，其中T508Y/Y603F催化2.7倍、t1/2 247.55 min（8倍）；以此为基础嵌M361L成三点M1（M361L/T508Y/Y603F），50 °C t1/2约500 min（16.4倍）、催化2.71倍，420 min残活>50%；四点延伸未超越M1，提示M1为局部适应峰。
2.3 酶学性质表征：野生型与M1最适pH均为7，M1在pH 5.5–8.5残活>50%（野生型6–8.5），全pH 5–10活力更高；最适温由35 °C移至40 °C，45 °C时M1残活48%而野生型近失活；55 °C下M1 120 min残活约50%（野生型t1/2 14 min），60 °C 40 min残活>20%，65 °C 5 min残活>60%（野生型3 min内灭）。DSF测得野生型Tm 50.8 °C，M1为54.9 °C（升约4 °C）。
2.4 分子动力学模拟：100 ns内M1的RMSD更快 Plateau于0.5–0.6 nm，波动低于野生型；RMSF在50–300及400–800残基区柔性普遍降低，近300位野生型峰超1.0 nm而M1消减；Rg更低稳，SASA峰值345 nm2（野生型）降至更窄带，氢键数稳态1–4条（野生型更少）。结构显示三突变均在DADH域：M361L距底物29.11 ?优化疏水核堆积，T508Y位于活性袋外表面芳香簇（Y492/F509）加固刚性，Y603F距8.36 ?调活性袋边缘疏水。外周-远端稳定而非活性中心直改，规避稳定性-活性权衡。
讨论与结论翻译：本研究采用半理性设计策略，以EVcouplings生成同源比对、GRAPE-WEB预测热稳残基，经单点筛选与组合获得三重突变体M1，其50 °C半衰期约500 min，为野生型16.36倍，分子动力学模拟揭示M1具更稳构象、更紧堆积及缓和表面暴露。M1在热稳提升同时催化效率改善，成为DON工业酶解候选，亦可作为二轮蛋白工程起点。该工作发表于《International Journal of Molecular Sciences》。
