来源：The AXXXA motif: A rational design strategy and experimental studies to enhance both activity and thermostability of BTL2 lipase（Enzyme and Microbial Technology）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
工业脂肪酶的热稳定性是制约其应用的关键瓶颈。本研究以嗜热脂肪杆菌BTL2脂肪酶为对象，通过计算设计与实验验证相结合的策略，利用AXXXA稳定化基序进行理性改造。研究团队从50个候选突变中筛选出N60A突变体，该突变将酶的最适反应温度从55°C提升至65°C，熔点温度（Tm）提高4.3°C，催化效率（kcat/Km）提升约2.1倍。该研究为脂肪酶工程改造提供了"计算筛选—实验验证"的完整范式。
研究背景
具有优异催化性能的酶已成为传统化学催化剂的理想替代品。目前提升蛋白质稳定性的策略主要有化学修饰、固定化和基因工程三种途径，其中基因改造应用最为广泛，涵盖静电相互作用优化、盐桥构建、靶向氨基酸替换以及引入AXXXA和GXXXG等稳定化基序多种策略。
前期对BTL2脂肪酶的研究已确认其碱性pH适应范围（pH 8–10）和中等热稳定性（60°C下半衰期约30分钟）。然而该酶仍面临多项技术挑战：在大肠杆菌中重组表达易形成包涵体、浓度超过1 mg/mL时发生聚集、温度超过60°C时催化活性急剧下降、有机溶剂耐受性差。Liu等人近期虽采用表面电荷工程提升了相关嗜热脂肪酶的性能，但如何同时改善催化活性和热稳定性仍是待解难题。本研究基于AXXXA基序的理性设计策略，旨在通过计算筛选和实验验证，实现BTL2脂肪酶稳定性与活性的同步提升。
研究方法
2.1 菌株与质粒
基因克隆和蛋白表达分别使用pGEM-5Zf(-)（Fermentas）和pET26b(+)（Invitrogen, USA）载体。大肠杆菌DH5α（Novagen）作为克隆宿主，大肠杆菌BL21(DE3)（Invitrogen, USA）作为重组蛋白表达宿主。
【数据】菌株：E. coli DH5α、E. coli BL21(DE3)；质粒：pGEM-5Zf(-)、pET26b(+)；抗性：卡那霉素（30 μg/mL）
2.2 计算模拟软件
BTL2脂肪酶的氨基酸序列（CAA64621.1）和核苷酸序列（X95309.1）获取自GenBank。结构建模采用Swiss-Model [25]，以PDB 2W22为模板，随后使用Galaxy网络服务器 [26] 进行结构优化。突变体脂肪酶结构通过以下工具验证：PROCHECK [27] 进行立体化学分析，ProSA-web [28] 进行能量评估，QMEAN [29]、QMEANDisCo [30] 和Errat程序 [31] 进行质量评估，lDDT [32] 进行局部结构验证。分子可视化和RMSD计算使用VMD 1.8.6 [33] 完成。与甘油三酯底物（三丁酸甘油酯C4、三辛酸甘油酯C8、三月桂酸甘油酯C12和橄榄油C18）的分子对接使用Molegro Virtual Docker [34] 进行。定点突变通过SOE-PCR [35] 实现，突变稳定性通过Eris服务器 [36]、I-Mutant 2.0 [37]、PopMusic [38]、Site-Directed Mutator (SDM) [39] 和CUPSAT [40] 进行ΔΔG计算评估。二级结构预测使用Jpred4 [41]。
【数据】模板：PDB 2W22；软件：Swiss-Model、Galaxy、PROCHECK、ProSA-web、QMEAN、QMEANDisCo、Errat、lDDT、VMD 1.8.6、Molegro Virtual Docker、Eris、I-Mutant 2.0、PopMusic、SDM、CUPSAT、Jpred4
2.3 突变筛选与ΔΔG计算
为在BTL2脂肪酶结构中构建AXXXA基序，首先鉴定酶的一级序列（GenBank: CAA64621.1）中的丙氨酸残基。随后选取每个丙氨酸残基相邻（前或后）的第四个氨基酸作为潜在的丙氨酸替换候选位点。使用Eris [36]、I-Mutant 2.0 [37]、PopMusic [38]、SDM [39] 和CUPSAT [40] 服务器计算每个突变的ΔΔG值来评估替换效应。最终N60A突变基于所有服务器的共识预测结果选定，该突变在维持催化完整性的同时展现出最有利的稳定化效果。
【数据】突变筛选数量：50个；软件：Eris、I-Mutant 2.0、PopMusic、SDM、CUPSAT；筛选标准：多服务器共识预测
2.4 结构建模与验证
N60A突变体脂肪酶结构通过同源建模生成，使用Swiss-Model服务器 [25]，以PDB 2W22为模板，随后通过Galaxy网络服务器 [26] 进行结构优化。综合验证包括：（i）使用PROCHECK [27] 进行立体化学质量评估；（ii）通过QMEAN [29] 和QMEANDisCo [30] 评估全局和局部模型准确性；（iii）使用ProSA-web [28] 验证能量分布；（iv）采用Errat程序 [31] 进行原子水平结构验证。这一严格的多参数验证确认了N60A突变体结构对后续功能分析的可靠性。
【数据】模板：PDB 2W22；验证工具：PROCHECK、QMEAN、QMEANDisCo、ProSA-web、Errat
2.5 对接分析
使用Molegro Virtual Docker (MVD) [34] 进行分子对接模拟，研究N60A突变对脂肪酶稳定性的结构影响。对接底物包括三丁酸甘油酯（C4，短链甘油三酯）、三辛酸甘油酯（C8，中链）、三癸酸甘油酯（C10，中链）和三月桂酸甘油酯（C12，中链）。所有对接程序均使用默认参数，聚焦于包含催化三联体残基（Ser114、His359和Asp318）的催化空腔（半径10 Å）。该分析为突变诱导的构象变化如何影响底物结合和酶稳定性提供了重要见解。
【数据】对接软件：Molegro Virtual Docker；底物：C4、C8、C10、C12甘油三酯；催化空腔半径：10 Å；催化三联体：Ser114、His359、Asp318
2.6 诱变与克隆策略
使用SOE-PCR以pTA-btl2为模板，通过三个连续PCR反应将N60A突变引入btl2基因：[1] PCR-I使用含有MluNI限制性位点的正向引物BTL2.F1（5′-GATGGCCATGGCGGCATCCCCACGCGCC-3′）和反向引物BTL2.N60A.R2（5′-CCCGGTCCCAGGCGCTCGAGAG-3′）扩增200 bp片段；[2] PCR-II使用正向引物BTL2.N60A.F3（5′-CTCTCGAGCGCCTGGGACCGGG-3′）和含有SacI限制性位点的反向引物BTL2.R4（5′-TTGAGCTCATCATCCCTTCATTAAGGCCGC-3′）生成1024 bp片段；[3] PCR-III将上述片段组合产生1224 bp产物，克隆至pGEM-5Zf(-)中构建pGEM-btl2。测序确认后，将1215 bp的MluNI/SacI双酶切片段亚克隆至pET26b(+)中生成pkG60-btl2，随后转化至大肠杆菌BL21(DE3)感受态细胞进行蛋白表达。
【数据】PCR-I片段：200 bp；PCR-II片段：1024 bp；PCR-III片段：1224 bp；亚克隆片段：1215 bp；限制性位点：MluNI、SacI；引物：BTL2.F1、BTL2.N60A.R2、BTL2.N60A.F3、BTL2.R4
2.7 重组脂肪酶的表达
挑取含有pkG60-btl2质粒的重组大肠杆菌BL21(DE3)单菌落，接种至10 mL含30 μg/mL卡那霉素的LB培养基（1%蛋白胨、0.5%酵母提取物、1% NaCl；pH 7.0）中，37°C、200 rpm振荡过夜培养。随后将3 mL预培养物转移至100 mL新鲜LB-卡那霉素培养基中，在相同条件下培养至对数中期（OD600 ≈ 0.6）。加入0.3 mM IPTG诱导蛋白表达，随后在28°C、降低搅拌速度（160 rpm）条件下过夜培养。通过离心（4000 × g，10 min，4°C）收集细胞，沉淀于−20°C保存以备后续纯化。
【数据】培养基：LB（1%蛋白胨、0.5%酵母提取物、1% NaCl，pH 7.0）；卡那霉素：30 μg/mL；预培养：10 mL，37°C，200 rpm，过夜；扩大培养：3 mL→100 mL，37°C，200 rpm；诱导：0.3 mM IPTG，28°C，160 rpm，过夜；收集：4000 × g，10 min，4°C；保存：−20°C
2.8 N60A脂肪酶的部分纯化
将细菌沉淀重悬于25 mL缓冲液I（50 mM Tris-HCl、100 mM NaCl、5 mM EDTA、1 mM PMSF、2 mM DTT、1% Triton X-100，pH 8.0）中，添加溶菌酶（100 μg/mL）和DNase I（10 μg/mL），室温孵育45 min。通过超声破碎（80%振幅，6 × 30 s脉冲）和离心（10,000 × g，15 min，4°C）实现细胞裂解。包涵体用缓冲液II（不含Triton X-100的缓冲液I）洗涤两次，并以5 mg/mL的浓度在缓冲液III（PBS中8 M尿素）中溶解，搅拌30 min。离心（13,000 × g，20 min）后，将上清液逐滴加入缓冲液IV（100 mM Tris、400 mM L-精氨酸-HCl、100 mM NaCl、100 mM甘氨酸、5%甘油、0.5 mM DTT，pH 8.5）中，持续搅拌进行缓慢复性。复性蛋白用缓冲液V（20 mM Tris，pH 7.4）透析，更换三次缓冲液，然后通过DE52阴离子交换层析纯化。随后通过电泳（SDS-PAGE）和免疫印迹检查纯化结果。
【数据】缓冲液I：50 mM Tris-HCl、100 mM NaCl、5 mM EDTA、1 mM PMSF、2 mM DTT、1% Triton X-100，pH 8.0；溶菌酶：100 μg/mL；DNase I：10 μg/mL；孵育：45 min，室温；超声：80%振幅，6 × 30 s；离心1：10,000 × g，15 min，4°C；溶解：8 M尿素，5 mg/mL，搅拌30 min；离心2：13,000 × g，20 min；复性缓冲液：100 mM Tris、400 mM L-精氨酸-HCl、100 mM NaCl、100 mM甘氨酸、5%甘油、0.5 mM DTT，pH 8.5；透析：20 mM Tris，pH 7.4，更换三次；纯化：DE52阴离子交换层析
2.9 动力学参数测定
以三辛酸甘油酯（C8）为底物测定野生型和N60A脂肪酶的动力学参数Km和Vmax。初始水解速率在55°C、pH 8.0条件下、底物浓度范围0.1–5.0 mM内测定。Km和Vmax值通过GraphPad Prism 9.0非线性回归拟合Michaelis-Menten方程计算。转换数（kcat）按kcat = Vmax / [E]计算，其中[E]为活性酶的摩尔浓度。
【数据】底物：三辛酸甘油酯（C8）；温度：55°C；pH：8.0；底物浓度范围：0.1–5.0 mM；软件：GraphPad Prism 9.0
2.10 脂肪酶活性测定
使用pH-STAT滴定系统（Metrohm Titrando 842）测定脂肪酶活性，采用四种不同链长的甘油三酯底物（短链C4、中链C8和C12、长链C18橄榄油）。每种底物（2 mg/mL）在含20 mg/mL阿拉伯胶的蒸馏水中通过超声均质（最大功率5 min）乳化 [23]。测定时取20 μL酶溶液（天然或突变脂肪酶）加入20 mL乳化底物中，于55°C、pH 8.0条件下进行，释放的脂肪酸用0.05 M NaOH连续滴定。一个酶活性单位定义为在上述标准化条件下每分钟释放1 μmol脂肪酸所需的脂肪酶量 [42]。
【数据】滴定系统：Metrohm Titrando 842；底物浓度：2 mg/mL；乳化剂：20 mg/mL阿拉伯胶；超声：最大功率5 min；酶量：20 μL；反应体积：20 mL；温度：55°C；pH：8.0；滴定液：0.05 M NaOH；酶活单位：1 U = 1 μmol脂肪酸/min
2.11 底物特异性分析
使用pH-STAT滴定系统（Metrohm Titrando 842）在标准化反应条件（pH 8.0、55°C）下系统评估天然和突变脂肪酶的底物特异性谱，采用四种代表不同链长的甘油三酯底物：短链（C4）、中链（C8、C12）和长链（C18橄榄油），实现对酶链长偏好的全面表征。
【数据】温度：55°C；pH：8.0；底物：C4、C8、C12、C18
研究结果
3.1 计算模拟研究
表1展示了在btl2脂肪酶基因中引入AXXXA基序后50个突变筛选结果及相应的ΔΔG计算值。结果表明，除I-Mutant外，所有服务器均确认将第60位天冬酰胺替换为丙氨酸（N60A）可增强酶稳定性。基于这些发现，选择N60A突变将天然的NWDRA基序转化为AWDRA基序，即将第60位天冬酰胺（N）替换为丙氨酸（A）。
选定突变后，使用参考文献 [25,26] 中的方法预测并优化了天然和N60A突变脂肪酶的三维结构（图1）。
在VMD [33] 中进行的结构比对表明天然与突变体结构之间具有高度构象保守性（图2）。
分析得到均方根偏差（RMSD）为0.197 Å，序列覆盖率为99.48%，表明N60A替换引起的结构扰动极小（表2）。
天然和N60A突变脂肪酶结构的立体化学质量验证结果汇总于表2。Ramachandran图分析显示，天然脂肪酶中94.2%的残基（311/389）位于最适区域，5.5%（18/389）位于额外允许区域，0.3%（1/389）位于不允许区域。相比之下，N60A突变体参数有所改善：94.5%的残基（312/389）位于最适区域，5.2%（17/389）位于额外允许区域，0.3%（1/389）位于宽松允许区域。值得注意的是，该突变将所有残基从不允许区域中消除。这些结果达到或超过了高质量结构的典型阈值（即不允许区域<1%），证实两个模型均具有优异的立体化学质量，且N60A替换未引入显著的构象扰动。
突变脂肪酶的IDDT（局部距离差异测试）评分为0.95（表2），表明突变体与天然酶之间具有高度的结构相似性。该值说明尽管引入了突变，突变脂肪酶的整体三维折叠和局部主链构象基本保持完整。
如表2所示，天然和N60A脂肪酶的QMEAN [29] 评分分别为0.52和0.54，QMEANDisCo评分 [30] 分别为0.87和0.85，表明突变未破坏整体蛋白质折叠。通过ERRAT（评估非键原子相互作用）进一步验证，两个结构获得近乎相同的评分97.61（天然：97.6127；N60A：97.6064），证实了模型的高质量。综合这些计算分析表明，N60A替换在成功引入AWDRA基序的同时保持了脂肪酶的结构完整性。
表3展示了天然和N60A突变脂肪酶预测二级结构的比较。分析显示，单氨基酸替换使突变体的α-螺旋含量较天然酶有所增加。观察到的增加提示可能发生构象变化，从而增强酶稳定性或调节其催化特性。
使用Molegro Virtual Docker (MVD) 进行的分子对接模拟显示，与三丁酸甘油酯（C4，短链甘油三酯）、三辛酸甘油酯（C8，中链）、三癸酸甘油酯（C10，中链）和三月桂酸甘油酯（C12，中链）底物对接时，突变脂肪酶的ΔΔG值较天然酶有所下降（表2）。MVD分析表明，突变导致更有利的结合能，提示突变体活性位点中底物-酶复合物的稳定性增强。
【数据】RMSD：0.197 Å；序列覆盖率：99.48%；Ramachandran图（天然）：最适区域94.2%（311/389）、额外允许5.5%（18/389）、不允许0.3%（1/389）；Ramachandran图（N60A）：最适区域94.5%（312/389）、额外允许5.2%（17/389）、宽松允许0.3%（1/389）；IDDT：0.95；QMEAN：0.52（天然）/0.54（N60A）；QMEANDisCo：0.87（天然）/0.85（N60A）；ERRAT：97.6127（天然）/97.6064（N60A）
3.2 实验结果
#### 3.2.1 重组脂肪酶的纯化与动力学分析
N60A突变BTL2脂肪酶在大肠杆菌中表达，通过包涵体复性和阴离子交换层析纯化。SDS-PAGE确认约43 kDa处出现单一条带（图3）。
使用三辛酸甘油酯（C8）为底物，在65°C、pH 8.0条件下测定Michaelis-Menten参数（表3）。N60A突变体表现出更低的Km和显著更高的Vmax和kcat。这导致催化效率（kcat/Km）提高约2.1倍。
【数据】分子量：约43 kDa；温度：65°C；pH：8.0；底物：三辛酸甘油酯（C8）；催化效率提升：约2.1倍
#### 3.2.2 底物特异性
图4展示了天然和突变脂肪酶在酰基链长为4、8、12和18个碳的底物存在下的活性。统计分析（单因素方差分析）显示，N60A突变体对所有底物的活性均显著高于野生型（C4、C8、C12的p < 0.01；C18的p < 0.05）。结果表明，引入的突变未影响N60A脂肪酶的底物特异性。此外，研究结果显示N60A突变增强了酶对4、8、12和18碳底物的水解活性。然而，随着酰基链长度增加，两种酶的水解活性均降低（表4）。
【数据】底物：C4、C8、C12、C18；显著性：C4/C8/C12的p < 0.01，C18的p < 0.05；趋势：酰基链增长→水解活性降低
#### 3.2.3 温度对脂肪酶活性和稳定性的影响
图5a展示了温度对两种脂肪酶活性的影响。结果表明，将第60位天冬酰胺替换为丙氨酸使酶的最适温度从55°C移至65°C。图5b展示了酶的热稳定性研究结果。结果表明，在65°C以上，虽然两种脂肪酶的活性均下降，但天然酶的活性下降更为迅速。热失活遵循一级动力学。野生型BTL2在65°C下的半衰期（t½）为28.4 min，而N60A的t½为68.7 min（提高2.4倍）。在70°C下，野生型t½为8.2 min，N60A为24.5 min。55°C、60°C、65°C和70°C下的完整t½值汇总于表5。
【数据】最适温度：55°C（天然）→65°C（N60A）；65°C半衰期：28.4 min（天然）vs 68.7 min（N60A），提高2.4倍；70°C半衰期：8.2 min（天然）vs 24.5 min（N60A）
#### 3.2.4 pH对天然和突变脂肪酶活性的影响
在pH 5–10范围内评估了天然和突变脂肪酶的酶活性（图6）。结果表明，两种脂肪酶的活性随pH从5升至8而增加，随后在pH 9和10时下降。值得注意的是，突变脂肪酶在所有测试pH值下均持续表现出高于天然酶的活性（pH 8–10时p < 0.01）。
【数据】pH范围：5.0–10.0；最适pH：8.0–8.5；显著性：pH 8–10时p < 0.01
#### 3.2.5 金属离子对天然和突变脂肪酶活性的影响
图7展示了天然和突变脂肪酶在2%（w/v）金属离子中暴露1小时后的酶活性。结果表明，虽然N60A突变体在所有测试金属离子存在下均表现出高于天然酶的活性（Ca²⁺、Mg²⁺的p < 0.05；Zn²⁺的p < 0.01），但两种酶在金属离子暴露后活性均有所降低。这表明尽管AXXXA基序增强了结构稳定性，但并未完全赋予金属离子抗性。
【数据】金属离子浓度：2%（w/v）；孵育时间：1 h；显著性：Ca²⁺/Mg²⁺的p < 0.05，Zn²⁺的p < 0.01
#### 3.2.6 有机溶剂对天然和突变脂肪酶活性的影响
图8展示了暴露于30%有机溶剂的天然和突变脂肪酶的活性分析。结果显示，两种酶在2小时暴露后活性均降低，突变体在2-丙醇中降低幅度最大，而天然酶在丙酮中降低幅度最大。突变体在所有溶剂中均持续保持较高活性（2-丙醇和丙酮的p < 0.01；乙醇和甲醇的p < 0.05），证明了AXXXA基序在有机介质中的稳定化作用。
【数据】有机溶剂浓度：30%（v/v）；暴露时间：2 h；显著性：2-丙醇/丙酮的p < 0.01，乙醇/甲醇的p < 0.05
#### 3.2.7 去污剂对天然和突变脂肪酶活性的影响
本研究评估了离子型去污剂（SDS）和多种非离子型去污剂（Tween 20、Tween 40、Tween 80和Triton X100）对天然和突变脂肪酶活性的影响（图9）。两种脂肪酶变体在2%（v/v或w/v）去污剂溶液中孵育2小时，随后测定其残留活性。结果表明，离子型去污剂SDS显著抑制了两种酶的活性。
【数据】去污剂浓度：2%（v/v或w/v）；孵育时间：2 h；去污剂种类：SDS（离子型）、Tween 20、Tween 40、Tween 80、Triton X100（非离子型）
讨论与解读
天然BTL2脂肪酶在Gln118（位于α-螺旋4的C端）和Asn60（位于α-螺旋2的C端）之间存在潜在的去稳定化相互作用，其侧链氧原子相距4.5 Å（图10）。
这种空间排列在这些极性残基之间产生静电排斥，导致天然酶的构象柔性增加。N60A突变体中将Asn60替换为丙氨酸消除了这一不利相互作用，同时保留了主链几何结构。这一策略性修饰在维持脂肪酶功能完整性的同时提高了其热力学稳定性，使N60A变体特别适用于需要稳健酶性能的工业应用。
IDDT评分0.95确认突变未破坏脂肪酶的整体结构完整性，这对维持其催化功能至关重要。高度结构保守性表明，观察到的底物结合改善（以降低的ΔΔG值表示）可能源于微妙、局部化的变化——如优化的侧链相互作用或微小的构象调整——而非大规模结构重排。这一发现支持了突变选择性增强活性位点底物稳定化而不损害酶整体稳定性的假说。
编译者解读：本研究展示了理性设计策略在脂肪酶工程中的精准应用——通过单个氨基酸替换（N60A）同时实现活性与稳定性的双重提升，体现了"结构基序导向"设计理念的实用价值。方法上，多软件共识筛选（ΔΔG计算）有效缩小了候选范围，降低了实验筛选成本。值得注意的是，突变体在有机溶剂和去污剂耐受性方面的改善拓展了其在非水相催化中的应用潜力。然而，本研究对N60A突变体的包涵体复性纯化依赖仍表明该酶在大肠杆菌中可溶性表达尚待优化，未来可结合分泌表达或分子伴侣共表达策略进一步提升工业适用性。
参考来源
Ghorban M. The AXXXA motif: A rational design strategy and experimental studies to enhance both activity and thermostability of BTL2 lipase. Enzyme and Microbial Technology, 2026.
DOI: 10.1016/j.enzmictec.2026.110884
DOI: 10.1016/j.enzmictec.2026.110884