来源：Semantic Scholar | 编译：DNA Lab Space
# 最小化实验筛选：半理性策略提升色氨酸2-单加氧酶热稳定性
导语
色氨酸2-单加氧酶（TMO）能够催化L-色氨酸转化为重要的植物生长调节剂吲哚-3-乙酰胺，但天然酶的热稳定性极差，严重制约其工业应用。传统酶工程改造常陷于规模庞大的实验筛选泥潭，而对TMO的稳定性设计更因结构-功能关系的复杂性而充满挑战。Kongjaroon等在《RSC Chemical Biology》发表的最新研究，巧妙运用以结构计算为先导的半理性设计策略，通过仅有数个变体的实验验证，即成功获得熔解温度跃升超过14°C且催化活力几乎不受损的耐热突变体，真正实现了“少实验、高命中”的酶改造目标。
研究背景
色氨酸2-单加氧酶（EC 1.13.12.3）是一种黄素依赖型氧化酶，催化L-色氨酸的氧化脱羧，生成吲哚-3-乙酰胺和二氧化碳。该产物是合成植物激素吲哚-3-乙酸的关键前体，在农业和园艺领域具有巨大潜力。与此同吋，TMO所参与的代谢网络也被认为与某些病原菌的致病机制以及抗肿瘤活性分子的生物合成存在关联，这使得TMO成为合成生物学与化学生物学交叉研究中的热点酶种。然而，绝大多数天然来源的TMO热稳定性极低，在稍高于室温的条件下即迅速失活。这一缺陷直接限制了其在体外催化、酶固定化以及连续流反应器中的实际应用，使得提升TMO热稳定性成为该领域亟待攻克的核心难题。
针对酶热稳定性的理性改造，研究者已发展出多种策略。定向进化虽然不依赖结构信息，但其成功往往依赖于海量突变文库的构建和高通量筛选手段，时间与资源成本高昂。纯粹的理性设计，如基于单一能量函数进行点突变推荐，则因对蛋白质构象动态和远端效应的预测能力不足，命中率长期偏低。半理性设计恰恰弥合了二者的间隙：它利用序列保守性分析、晶体结构、分子动力学模拟以及各种计算预测工具，将突变搜索空间大幅压缩至若干极具潜力的候选位点，随后仅需构建一个小型高质量的突变库进行实验验证。近年来，通过结合B因子分析、折叠自由能变化（ΔΔG）计算、二硫键工程以及共识序列设计等方法，多个酶的热稳定性已获得显著突破。然而，以最小化实验投入为前提，针对氧化酶类——尤其是需要维持精细活性构象和辅因子结合能力的TMO——开展系统半理性设计，此前尚未见报道。本研究正是以此切入点，意图建立一套高效快捷且普适性强的TMO稳定性工程改造流程。
研究方法
研究团队建立了一条从计算筛选到实验验证的闭环式半理性路线。首先从蛋白质数据库获取目标TMO的高分辨率晶体结构，若缺乏实验结构则采用AlphaFold2构建可靠模型。在开始突变设计之前，开展全原子分子动力学模拟，分别在生理温度（300 K）和升高温度（如320-340 K）下采集纳秒级轨迹，计算每个残基的均方根涨落（RMSF）及模拟B因子，锁定构象高度柔性的环区、loop区域以及热振动剧烈的二级结构末梢，这些区域被普遍认为是引发热失活的结构薄弱点。
随后，团队并行使用了三套预测策略以交叉验证，避免单一算法偏差。第一，采用FoldX的PositionScan模块对所有位点进行虚拟饱和突变扫描，计算突变体的折叠自由能差值（ΔΔG），保留ΔΔG < -1.5 kcal/mol的稳定化突变。第二，利用Rosetta的ddg_monomer模块和笛卡尔空间弛豫协议，进一步评估上述突变的能量效益，仅将通过双平台检验的候选列入初筛清单。第三，基于同源家族的超过150条序列的多序列比对，计算每个位置的氨基酸相对熵值和保守性评分，并通过共识序列分析，提取可能提高结构刚性的天然高频突变类型。此外，还特异性地调用Disulfide by Design 2.0工具判断引入二硫键的几何可能性，并对满足立体化学标准的位置对计算虚拟突变后的稳定化能。
所得所有候选突变按下述规则进行严格过滤：剔除位于活性中心5 Å范围内的残基，以避免对催化造成直接干扰；排除可能扰乱FAD辅因子结合界面的突变；删除过度改变表面静电势而不利于胶体稳定性的替换。最终，仅有8个位点的单点突变和2对协同突变组合脱颖而出，形成仅含10个成员的极小突变库。全基因合成这些变体，于大肠杆菌BL21(DE3)中异源表达，镍亲和层析一步
DOI: 10.1039/d4cb00102h