计算设计提升卤代烷脱卤酶稳定性的研究

作者： Robert J. Floor、Hein J. Wijma、Dana I. Colpa 等
通讯作者： Dick B. Janssen，格罗宁根大学（University of Groningen）
期刊： *ChemBioChem*，2014年，第15卷，1660–1672页

研究背景与目的

酶作为高效的生物催化剂，在工业应用中常因热稳定性不足和有机溶剂耐受性差而受到限制。卤代烷脱卤酶 LinB（haloalkane dehalogenase LinB）来源于 Sphingomonas paucimobilis UT26，能催化多种氯代、溴代和碘代烷烃的水解脱卤反应，在环境污染物降解和立体选择性生物催化方面具有重要价值，但其中等的稳定性限制了其工程应用。本研究采用此前已建立的 FRESCO（Framework for Rapid Enzyme Stabilization by Computation，基于计算的酶快速稳定化框架）计算设计策略，目的是在不显著损失催化活性的前提下，系统提升 LinB 的热稳定性和有机共溶剂耐受性。此前该策略已在柠檬烯环氧化物水解酶（limonene epoxide hydrolase）上取得成功，将表观熔解温度（apparent melting temperature，T_m,app）提高了 32 °C。

研究流程

研究共分为五个主要阶段：(1) 远端点突变的设计与筛选；(2) 二硫键的设计与筛选；(3) 稳定突变的组合；(4) 活性位点附近突变的引入；(5) 组合突变体的功能验证与应用评估。

远端点突变设计。 首先排除距底物叔丁基-2-(2-溴丙酰胺基)乙酸酯（docked substrate tert-butyl-2-(2-bromopropanamido)acetate）9 Å 以内的 63 个残基，保留 233 个可替换残基以避免直接影响活性位点。利用 FoldX 和 Rosetta 两种程序预测所有单点突变的折叠自由能变化（ΔΔG_fold），其中 FoldX 预测出 150 个潜在稳定突变，Rosetta 预测出 314 个能量高于 3.5 kJ/mol 的突变，两者共同预测 57 个，共得到 407 个候选突变。经人工结构检查排除 219 个不可信突变，再进行分子动力学（molecular dynamics，MD）筛选，每个突变体运行五条 100 ps 的平行模拟，排除 79 个在动态模拟中表现出不利相互作用的突变，最终选定 109 个点突变进行实验验证。通过 QuikChange 定点突变在 96 孔板中构建突变体，利用亲和层析纯化，采用热荧光法（ThermoFluor）测定 T_m,app。以实验标准差的两倍（0.96 °C）作为稳定判据，从 99 个成功表达的可溶性突变体中发现 10 个稳定突变，成功率为 10%。其中 A197P 引入脯氨酸到环区，E15T、D166N 和 D255A 去除表面负电荷，A81K、D166K、E192K 和 E192R 引入表面正电荷，A53L 改善内部疏水堆积，G229Q 替换环区甘氨酸。活性测定显示几乎所有远端突变均保持与野生型相近的脱卤活性。

二硫键设计。 使用动态二硫键发现算法（dynamic disulfide discovery，DDD），同时结合 X 射线晶体结构和分子动力学模拟快照寻找适合引入二硫键的位置。基于晶体结构仅发现 7 对合适残基；利用 45 个 MD 快照进行构象采样后，额外发现 25 对，总数达 32 对。经人工检查和 MD 筛选排除 19 个可能破坏稳定的设计，剩余 13 个在 E. coli TOP10 和有利于胞质二硫键形成的 SHuffle 菌株中表达验证。仅 A5C/A185C 这一对二硫键变体在氧化态下稳定，提高 T_m,app 5 °C，活性与野生型相当。六种变体在氧化态下出现错误折叠，表明引入的半胱氨酸与内源半胱氨酸之间形成了非预期的二硫键。

稳定突变的组合。 将三个远端稳定点突变（G229Q、D166K、E192K）和一个二硫键（A5C/A185C）组合构建变体 G1，其 T_m,app 比野生型高 11 °C（62 °C）。动力学参数显示 G1 对 1-溴丙烷的 k_cat 与野生型相似，K_m 略有升高，总体催化效率基本保持。对 11 种不同卤代烷的活性测定均与野生型相当。热失活实验显示 G1 在 50 °C 的半衰期延长 40 倍，在 54 °C 延长 20 倍；最适反应温度由野生型的 40 °C 提升至 50 °C，且在最适温度下 G1 的活性为野生型的 1.5 倍。

活性位点附近突变。 由于活性位点覆盖区域（α4 螺旋第 144–148 位残基和 α5 螺旋第 169–173 位残基）具有较高的 B 因子（B-factor）和均方根波动（root-mean-square fluctuation，RMSF），属于柔性区域，研究进一步在 G1 背景中引入距结合氯离子 9 Å 以内的突变。经计算流程（能量预测、结构检查和 MD 筛选）获得 27 个候选突变，24 个成功表达，其中 7 个稳定，成功率为 29%，远高于远端突变的 10%。五种突变改善了活性位点附近疏水内部的堆积（D147M、D147L、F169V、A247F、T249L），其中 D147 位点的四种疏水替换（H、Y、M、L）使 T_m,app 提高 3.0–5.5 °C，但均显著降低催化活性；仅 F169V 和 T249L 在提高稳定性的同时保持了活性。

最终稳定变体 G3。 将远端突变 E15T/A53L/A81K/A197P/D255A 和近活性位点突变 F169V/A247F 进一步引入 G1，得到含 10 个点突变和 1 个二硫键的变体 G3。差示扫描量热法（differential scanning calorimetry，DSC）和热荧光法均测得 T_m,app 为 74 °C，比野生型高 23 °C。G3 在 60 °C 的半衰期为 200 min，而野生型在 55 °C 半衰期仅为 1.0 min。在最适温度 55 °C 时 G3 的活性与野生型在其最适温度 40 °C 时的活性相当。尽管 G3 对部分底物的活性较野生型降低 2–10 倍，但其获得的高温下催化能力超过了野生型。

共溶剂耐受性与应用。 G1 和 G3 在 50% DMSO 和 25% 1,4-二氧六环中的稳定性均显著优于野生型。在含 25% DMSO 的体系中，G1 于 45 °C 下 8 h 内完全转化 180 μM 的 β-六氯环己烷（β-HCH），而野生型仅转化 26% 并在 30 min 内失活。在含 10% THF 的条件下，G1 能将底物浓度提高三倍（5 mM）应用于叔丁基-2-(2-溴丙酰胺基)乙酸酯的动力学拆分，持续工作超过 72 h，转化率达 53%，对映选择性 E 值为 28–30（野生型无共溶剂时 E 值为 67），而野生型在 5 h 内转化不足 3%。

结论与价值

本研究成功将 FRESCO 计算策略应用于较大的单体酶 LinB，仅实验测试约 150 个突变体，即获得 18 个稳定突变和一个热稳定性提高 23 °C 的最终变体 G3。研究发现：提高柔性区域的局部折叠稳定性对整体稳定性贡献最大；引入二硫键可有效阻止 N 端的解析折叠；共溶剂耐受性与热稳定性之间存在正相关关系。应用价值在于，稳定变体 G1 和 G3 显著提升了对 DMSO、THF 等有机共溶剂的耐受性，使其能够高效降解水中溶解度极低的环境污染物 β-HCH，并实现以较高底物浓度进行的对立选择性生物催化。

研究亮点

本研究的突出之处在于：一是仅用计算筛选取代了绝大部分实验筛选，实验测试数量不到计算评价突变数的 3%，显著降低了成本；二是系统地比较了远端突变与活性位点附近突变对稳定性和活性的影响，明确了柔性区域（尤其是活性位点覆盖螺旋）的折叠稳定性是整体稳定性的限速环节；三是证实了两个计算程序（FoldX 和 Rosetta）联合使用并结合 MD 筛选可有效识别稳定突变，同时指出表面暴露的疏水残基错误预测是计算方法的主要局限之一。