来源：Unraveling the tunnel dynamics of enzymes: a rational design strategy for robust biocatalysis in organic solvents（Food Chemistry）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
有机溶剂中酶的活性与稳定性往往此消彼长，成为工业生物催化的核心瓶颈。本研究提出"隧道动力学引导"的工程策略，通过分子动力学模拟揭示溶剂诱导脂肪酶产生"防御性压缩"、形成病态狭窄隧道的失活机制，并精准锁定八个非保守门控位点进行改造。所得突变体在66.7% DMSO和10%吡啶中的耐受性分别提升7.53倍和5.56倍，并在公斤级蔗糖月桂酸酯合成中实现96.5%转化率，为下一代工业生物催化剂的理性设计提供了新范式。
研究背景
生物催化凭借环境友好、可持续、安全和高特异性等优势，已成为传统有机合成的重要替代方案。然而，其工业应用常受限于反应环境中酶的催化效率不足和稳定性有限。核心瓶颈在于"溶解度-活性悖论"：极性底物在非极性有机溶剂中溶解性差，而酶在极性介质中活性又往往下降。有机溶剂通常通过多种机制损害酶功能：诱导构象变化、竞争酶表面必需结合水、竞争性抑制底物与活性位点结合、改变底物溶解度以及 destabilizing 酶-底物过渡态。尽管吡啶和DMSO等强极性溶剂能提供优异的底物溶解度、理论上可提高反应效率，但它们常是一把双刃剑——在许多情况下会剥夺酶表面的必需水分。
研究方法
2.1 材料与试剂
质粒提取试剂盒购自生工生物（上海，中国）；凝胶回收试剂盒购自Omega Bio-tek（Norcross, GA, USA）。氨苄青霉素（分析纯）、氯化钙（分析纯）、异丙基β-D-1-硫代半乳糖苷（IPTG，98%）、三(羟甲基)氨基甲烷（98%）和月桂酸乙烯酯（98%）均购自麦克林（上海，中国）。酵母提取物（分析纯）购自诺维信（北京，中国）。胰蛋白胨（分析纯）和蔗糖（98%）购自国药集团化学试剂有限公司（上海，中国）。氯化钠（分析纯）和琼脂粉（分析纯）购自Sigma-Aldrich（St. Louis, MO, USA）。氢氧化钠（分析纯）、磷酸二氢钠（98%）、磷酸氢二钠（98%）和琼脂糖（98%）由南京万青化工玻璃仪器有限公司提供（南京，中国）。甘油（98%）、乙二胺四乙酸（EDTA，98%）和乙酸（98%）购自安徽泽升科技有限公司（安庆，中国）。
【数据】试剂来源：生工生物、Omega Bio-tek、麦克林、诺维信、国药集团、Sigma-Aldrich、南京万青、安徽泽升；纯度：IPTG 98%、月桂酸乙烯酯 98%、蔗糖 98%、EDTA 98%等
2.2 实验与计算方法
#### 2.2.1 E. coli BL21(DE3)感受态细胞制备
取3 μL商业化E. coli BL21(DE3)感受态细胞划线接种于Luria-Bertani（LB）琼脂平板，37°C过夜培养。挑取单菌落接种至LB液体培养基中，于37°C、210 rpm摇床过夜培养。次日，将1 mL过夜培养物转移至含100 mL新鲜LB培养基的250 mL锥形瓶中，于37°C、210 rpm振荡培养，直至600 nm处光密度达到0.6。随后将锥形瓶置于冰上冷却30 min，细胞于4°C、12,000 rpm离心10 min收集。弃去上清，细胞沉淀用10 mL冰预冷的0.1 M CaCl₂溶液重悬，冰上孵育10 min。再次于4°C、8,000 rpm离心10 min收集细胞，重复该洗涤步骤一次。最后，细胞沉淀用4 mL含15%（v/v）甘油的冰预冷0.1 M CaCl₂溶液重悬。感受态细胞悬液分装为100 μL等份，于−80°C保存备用。
【数据】培养温度：37°C；摇床转速：210 rpm；OD600：0.6；离心条件：12,000 rpm/10 min/4°C、8,000 rpm/10 min/4°C；CaCl₂浓度：0.1 M；甘油浓度：15%（v/v）；分装体积：100 μL；保存温度：−80°C
#### 2.2.2 琼脂糖凝胶电泳
PCR扩增产物（通常取5 μL）与1 μL 6× DNA上样缓冲液混合，在含1× TAE（Tris-乙酸-EDTA）缓冲液的0.8%（w/v）琼脂糖凝胶上进行电泳分离。电泳在90 V恒定电压下运行50 min。分离完成后，DNA条带用4S Green核酸染料染色，并在蓝光透射仪下用凝胶成像系统观察。通过将相应条带与标准DNA ladder比对，确认特异性扩增片段大小。
【数据】琼脂糖浓度：0.8%（w/v）；电泳电压：90 V；电泳时间：50 min；缓冲液：1× TAE；上样量：5 μL PCR产物 + 1 μL 6× loading buffer
#### 2.2.3 重组质粒转化
经凝胶验证和纯化后，将突变重组质粒转化至E. coli BL21(DE3)感受态细胞中。凝胶纯化后的DNA用去离子水稀释至80 ng/μL浓度。在超净工作台中，取10 μL DNA加入100 μL刚刚解冻的感受态细胞中，轻轻混匀，冰上孵育25 min。随后在恒温混匀仪上42°C热激90 s，立即冰浴5 min。接着加入500 μL不含抗生素的LB培养基，于37°C、350 rpm振荡培养60 min。孵育结束后，3,000 rpm离心30 s收集细胞。弃去500 μL上清，用剩余培养基重悬细胞沉淀。将悬液均匀涂布于含氨苄青霉素（100 μg/mL）的LB平板上，倒置于37°C恒温培养箱中过夜培养。
【数据】DNA浓度：80 ng/μL；DNA用量：10 μL；感受态细胞用量：100 μL；冰浴：25 min；热激：42°C/90 s；恢复培养：37°C/350 rpm/60 min；离心：3,000 rpm/30 s；氨苄青霉素浓度：100 μg/mL；培养温度：37°C
#### 2.2.4 突变体表达
突变体库的表达和裂解在96深孔板中进行。首先，将预先用75%乙醇浸泡的96深孔板在超净工作台中晾干。每孔加入400 μL含氨苄青霉素（100 mg/L）的LB液体培养基。用无菌牙签从转化平板上挑取单菌落（包括野生型WT对照），插入相应孔中；对于饱和突变文库，接种全部96孔以确保完全覆盖。牙签放置1 min后取出。该板于37°C、300 rpm振荡培养12–14 h。过夜培养后，在新的无菌96深孔板中加入500 μL TB液体培养基（含100 mg/L氨苄青霉素）。用多通道移液器将20 μL过夜LB培养物转移至新鲜TB培养基中。这些板于37°C、300 rpm培养3 h。随后向每孔加入20 μL IPTG溶液（2.65 mmol/10 mL）诱导蛋白表达。摇床条件调整为18°C、300 rpm，诱导24 h。为裂解细胞，将板从摇床取出，每孔加入20 μL溶菌酶溶液（10 mg/mL），终浓度为0.4 mg/mL。涡旋确保细胞完全重悬，冰上孵育30 min。随后在37°C水浴中热激5 min，再冰浴30 min以释放胞内粗酶。最后，板于4,000 rpm离心30 min。收集含有粗酶的上清液，保存备用。
【数据】LB培养基体积：400 μL/孔；氨苄青霉素浓度：100 mg/L；培养条件：37°C/300 rpm/12–14 h；TB培养基体积：500 μL/孔；接种量：20 μL；诱导前培养：37°C/300 rpm/3 h；IPTG用量：20 μL（2.65 mmol/10 mL）；诱导条件：18°C/300 rpm/24 h；溶菌酶终浓度：0.4 mg/mL；冰浴：30 min；水浴热激：37°C/5 min；离心：4,000 rpm/30 min
#### 2.2.5 溶剂耐受性测定
分别配制10%（v/v）吡啶、50%（v/v）甲醇溶液、66.7%（v/v）DMSO溶液和50%（v/v）异丙醇溶液，用于评估突变体的溶剂耐受性。根据先前研究（Chong et al., 2024），采用pNPP法测定CALB变体的性能。首先，将40 μL酶溶液与20 μL配制好的溶液混合，在96孔板摇床中于40°C、500 rpm孵育18 min。孵育结束后立即加入180 μL pNPP A/B混合液，于40°C、500 rpm反应10 min。在OD 410处测定吸光度。该数据集代表脂肪酶的残余活性。初始活性测定方法为：用等体积TB培养基替代吡啶水溶液，并遵循相同的后续孵育和反应程序（表S1）。脂肪酶的溶剂耐受性可用上述测得的pNPP水解酶活性表示，计算公式如下：
残余活性% =（酶变体或WT在有机溶剂中的活性）/（酶变体或WT在缓冲液中的活性）× 100%
最后，选择溶剂耐受性显著提高的突变体进行复筛以进一步确认其性能。溶剂耐受性经再次确认后，突变体在LB培养基中培养并测序。菌株与50%甘油混合后，长期保存于−80°C冰箱中供后续研究使用。
【数据】溶剂体系：10%吡啶、50%甲醇、66.7% DMSO、50%异丙醇；孵育条件：40°C/500 rpm/18 min；反应条件：40°C/500 rpm/10 min；检测波长：OD 410；酶液用量：40 μL；溶剂用量：20 μL；pNPP底物用量：180 μL
#### 2.2.6 变体验证
为验证工程变体在不同溶剂体系中的催化性能，在15 mL反应小瓶中进行了三种不同的酯化反应。首先，蔗糖月桂酸酯的合成：将0.087 mmol蔗糖和0.348 mmol月桂酸乙烯酯与400 μL酶溶液混合于2 mL吡啶中。葡萄糖月桂酸酯的合成：将0.056 mmol葡萄糖和0.166 mmol月桂酸乙烯酯与400 μL酶溶液混合于2 mL DMSO/2M2B共溶剂混合物（4:1，v/v）中。最后，维生素E酯化反应：将0.39 mmol维生素E和1.19 mmol琥珀酸酐与400 μL酶溶液混合于2 mL DMSO中。所有反应均在磁力搅拌器上于40°C、400 rpm振荡24 h。
【数据】蔗糖用量：0.087 mmol；月桂酸乙烯酯用量：0.348 mmol（蔗糖酯化）/0.166 mmol（葡萄糖酯化）；葡萄糖用量：0.056 mmol；维生素E用量：0.39 mmol；琥珀酸酐用量：1.19 mmol；酶液用量：400 μL；溶剂体系：2 mL吡啶、2 mL DMSO/2M2B（4:1）、2 mL DMSO；反应条件：40°C/400 rpm/24 h
研究结果
3.1 分子动力学模拟揭示隧道动态模式
为阐明CALB在有机溶剂环境中的失活机制并为后续理性设计奠定基础，我们选择吡啶作为代表性有机溶剂，以水为对照体系，进行了100 ns分子动力学（MD）模拟，重建CALB的动态构象景观。尽管酶活性丧失传统上归因于整体蛋白质变性，但我们的均方根偏差（RMSD）分析显示，CALB在吡啶中虽比在水中表现出更高的结构波动（图S1a），但其整体折叠并未完全解体。相反，溶剂可及表面积（SASA）和回转半径（Rg）的降低揭示了一种更隐蔽的"防御机制"（图S1c、d）——酶通过采取更紧密的整体压缩构象来最小化与非极性溶剂的表面接触。
然而，这种溶剂诱导的"防御机制"（Segura et al., 2012）以关键功能区域的可及性为代价。最显著的结构重排发生在两个关键区域：底物结合口袋和通道入口。Fpocket分析（Schmidtke, Le Guilloux, Maupetit, & Tuffïry, 2010）显示，CALB的活性口袋体积从水相中的559.0 ų压缩至吡啶中的530.2 ų（表S3、S4）。这种物理阻塞直接阻碍了底物的有效锚定，从而抑制酶活性。
关键的是，使用Caver Analyst 2.0（Jurcik et al., 2018）对隧道动力学进行的定量分析阐明了底物转运受阻的微观机制。我们观察到有机溶剂显著改变了隧道的几何拓扑结构，平均瓶颈半径（Avg_BR）从1.665 Å显著收缩至1.411 Å，意味着隧道门被有效锁定在关闭状态。同时，一个极具揭示性的发现是隧道路径的显著延长——平均长度（Avg_L）从3.695 Å激增至6.033 Å。这种向"更窄更长"隧道构型的演变不仅增加了底物进入活性中心的空间位阻，还显著延长了扩散路径，最终导致理论底物通量的大幅下降（表1，图S2）。
均方根波动（RMSF）分析进一步定位了导致这种构象畸变的关键残基。特定区域（如残基42、47和144）在吡啶中表现出的异常柔性波动（图S1b）构成了观察到的瓶颈收缩和路径重建的结构基础。这些位于隧道入口和活性口袋外围的氨基酸残基充当着决定底物进入效率的"门控开关"。
总之，CALB在强极性有机溶剂中的失活不仅源于整体结构崩溃，更源于活性口袋的体积压缩和隧道动态构型的严重失衡。阐明这一机制为理性设计提供了清晰蓝图——在有机溶剂中恢复酶活性需要超越全局稳定性的工程策略，聚焦于精确重塑底物隧道和结合口袋的几何特征。
【数据】MD模拟时长：100 ns；活性口袋体积：559.0 ų（水相）→ 530.2 ų（吡啶）；Avg_BR：1.665 Å → 1.411 Å；Avg_L：3.695 Å → 6.033 Å；关键柔性残基：42、47、144
3.2 隧道动力学引导的突变体库构建
前期研究已证实底物通道是脂肪酶独特的结构特征，有助于底物和溶剂向活性位点的有效转运。在此转运过程中，氢键、静电力和疏水作用等特定分子间相互作用显著影响酶的溶剂耐受性（Koudelakova et al., 2013）。因此，修饰隧道结构并在活性位点内替换有利残基已成为增强溶剂抗性的有效策略。
基于此前提，我们首先使用Fpocket软件预测活性口袋景观，如图S3和表S4所示。该分析在CALB中鉴定出七个潜在底物结合口袋，其中Pocket 1为主要结合位点，包含G39、T40、T42、K47、W104、S105、D134、T138、L144、V154、Q157、E188、I189、V190、H224、A281、A282和I285等残基。这一鉴定结果与先前文献报道一致。
作为口袋分析的补充，我们利用Caver Analyst 2.0解析CALB隧道的动态构型，揭示了四条不同的底物通路，如图2所示。三条主要隧道的残基组成被详细绘制，其中第一条隧道由T40、D134、T138、L144、D145、V154、Q157、I189、H224和I285等残基定义。
然而，不加选择地诱变存在结构去稳定化的风险，尤其是在对蛋白质完整性至关重要的高度保守位点。为降低此风险，我们使用ConSurf软件（Ashkenazy et al., 2016）进行了进化保守性分析，筛选出保守性评分为6或更高的残基，从而排除可能导致酶失活的突变位点。通过将已鉴定的关键隧道和口袋残基与进化可变位点求交集，我们成功将选择范围缩小至八个具有高工程潜力的靶标残基，即T42、K47、L144、D145、V154、I189、A282和I285。这些位点策略性地位于口袋入口或隧道瓶颈处，且表现出足够的进化耐受性，可作为理想的功能开关。据此，我们构建了针对这些位置的饱和突变文库，目标是精确调节侧链体积、电荷和疏水性，以重新打开被溶剂压缩的隧道，实现活性和稳定性的双重提升。
【数据】保守性筛选阈值：评分≥6排除；靶标残基：T42、K47、L144、D145、V154、I189、A282、I285；预测口袋数量：7个；隧道数量：4条
3.3 具有优异溶剂耐受性的CALB工程改造
为建立稳健的高通量筛选方法，我们最初选择吡啶作为模型溶剂，因其对脂肪酶活性具有众所周知的损害作用（Cheng & Nian, 2023; Ramon Tomé Oliveira et al., 2025）。我们制备了浓度为5%、10%和50%的吡啶-水混合物，与野生型（WT）酶孵育以评估残余活性。结果表明，5%吡啶处理后WT的残余活性为46.2%，而10%吡啶溶液将活性维持在17.7%。在微量滴定板（MTP）上进行的多个平行实验证实了这些结果的高度可重复性，而50%吡啶则导致蛋白质变性沉淀。因此，我们选择10%吡啶水溶液作为溶剂耐受性筛选的标准条件。
筛选方法验证后，我们以WT为模板进行饱和突变。利用预实验中确定的吡啶浓度，测定了八个突变位点的溶剂耐受性。如图3a所示，对这些位点的筛选获得了20个溶剂耐受性显著提高的突变体，所有突变体在10%吡啶水溶液孵育后均表现出比WT更高的残余活性。尽管所有选定的突变体都表现出增强的耐受性，但不同突变体以及同一位置不同替换之间观察到相当大的溶剂抗性差异。值得注意的是，42位（包括T42K和T42Q）和285位（包括I285F、I285S、I285N、I285D等）的突变体表现尤为突出。
为进一步探索所选工程位点对CALB吡啶耐受性的影响，我们在双突变体NF、ND、TF和PF的基础上引入了额外突变。遗憾的是，第三个突变的引入通常导致耐受性普遍下降，只有KND保持了双突变体的效力，较WT提高5.4倍。进一步研究表明，这些突变体在各种酯类合成中表现出优异的性能。
【数据】5%吡啶残余活性：46.2%；10%吡啶残余活性：17.7%；筛选获得突变体数量：20个；KND较WT提高：5.4倍
讨论与解读
本研究针对脂肪酶在有机溶剂中面临的"活性-稳定性权衡"这一普遍科学挑战，成功建立并验证了一种名为"隧道动力学引导工程"的理性设计范式。通过整合MD模拟揭示隧道狭窄效应的机制见解与Fpocket、Caver和ConSurf等多维计算设计工具，从复杂蛋白质结构中精准定位了八个非保守门控位点，即残基42、47、144、145、154、189、282和285。实验结果表明，双突变体在吡啶耐受性方面表现出显著改善，其中145和285位的重组尤为突出。具体而言，突变体NF和ND在10%吡啶溶液中的耐受性分别比WT高5.46倍和5.56倍。双突变体在66.7% DMSO中耐受性提升7.53倍。机制分析鉴定出新形成的D285-K290盐桥，该盐桥刚性化了隧道入口并将底物路径缩短约40%，显著减轻了传质阻力。公斤级蔗糖月桂酸酯合成验证（96.5%转化率）确立了隧道动力学重塑作为设计下一代工业生物催化剂的稳健范式。
编译者解读：本研究巧妙地将"隧道动态"这一微观结构特征转化为可操作的工程靶点，跳出了传统"全局稳定性"改良的思维定式。从合成生物学应用角度看，该方法论的核心价值在于其可迁移性——隧道门控位点的识别逻辑可推广至其他具有底物通道的酶类。公斤级反应验证是亮点，直接回应了工业应用对规模化数据的需求。值得注意的是，三突变体普遍活性下降的现象提示多靶点改造存在非线性效应，未来或可结合机器学习预测组合突变的最优路径。D285-K290盐桥的机制解析也为定向进化提供了明确的选择压力设计依据。
参考来源
Cheng, W. (2026). Unraveling the tunnel dynamics of enzymes: a rational design strategy for robust biocatalysis in organic solvents. Food Chemistry, 2026. DOI: 10.1016/j.foodchem.2026.150333
DOI: 10.1016/j.foodchem.2026.150333