来源：Mechanistic insights into activity-enhancing mutations from comparative analysis of two Bacillus subtilis GH11 xylanases（International Journal of Biological Macromolecules）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
GH11木聚糖酶是生产低聚木糖的关键生物催化剂。本研究从枯草芽孢杆菌中筛选到两个仅存在八个氨基酸差异的GH11木聚糖酶BsXyn1和BsXyn2，通过定点突变鉴定出三个能显著提升酶活的关键位点（Y79I、V81C、S100G）。其中BsXyn1V81C突变体比活达14,940 U/mg，较野生型提升2.46倍，但热稳定性有所下降。分子动力学模拟揭示了thumb区域构象柔性增加与催化效率提升之间的关联，为木聚糖酶理性设计提供了新视角。
研究背景
低聚木糖（XOS）是一类具有调节肠道菌群、促进矿物质吸收等生物活性的功能性寡糖，在功能食品领域应用前景广阔。然而其天然含量低，促使研究者开发以木质纤维素生物质为原料的可持续生产平台。酶水解法因底物特异性强、副产物少、环境友好而成为首选策略。木聚糖酶（EC 3.2.1.8）是酶法生产XOS的关键酶，可将木聚糖水解为木二糖、木三糖等功能性组分。
糖苷水解酶家族11（GH11）木聚糖酶因底物特异性高、催化效率优异、分子量小（约20 kDa），被誉为"真正的木聚糖酶"。其结构采用特征性的β-果冻卷折叠，由两个反向平行β-折叠片构成深底物结合裂缝。GH11木聚糖酶含有一个独特的thumb环状区域，连接两个β-折叠片并含有保守的Pro-Ser-Ile基序。突变研究表明，增加thumb区域柔性可增强酶与底物的相互作用，但往往以牺牲稳定性为代价。
研究方法
2.1 化学试剂与材料
质粒构建使用大肠杆菌TOP10F'（Invitrogen, Carlsbad, CA），蛋白表达在大肠杆菌BL21(DE3)（Novagen, Madison, WI, USA）中进行。表达载体pF58来源于前期工作[25,26]。PCR扩增使用KOD FX DNA聚合酶（TOYOBO, Tokyo）。羧甲基纤维素钠（CMCNa, A501427-0250）和苯甲基磺酰氟（PMSF, A610425-0005）购自生工（上海，中国）。微晶纤维素（Avicel PH101, 11365-1 kg）购自Sigma-Aldrich（St. Louis, MO, USA）。海藻酸钠（S902506-500 g）和壳聚糖（C915936-100 g）购自Macklin（上海，中国）。玉米芯木聚糖、桦木木聚糖和稻壳木聚糖分别从玉米芯、桦木和稻壳中提取。其他化学试剂均为分析纯级，购自生工。所有寡核苷酸引物由通用生物（安徽，中国）商业合成，DNA测序服务由生工提供。
【数据】菌株：E. coli TOP10F'、E. coli BL21(DE3)；载体：pF58；试剂：CMCNa、PMSF、Avicel PH101、海藻酸钠、壳聚糖；底物：玉米芯木聚糖、桦木木聚糖、稻壳木聚糖
2.2 质粒构建
两个枯草芽孢杆菌木聚糖酶（BsXyn1和BsXyn2）的氨基酸序列从UniProt数据库获取（登录号分别为P18429和Q3HLJ4）。使用JCat软件（http://www.jcat.de/）对相应基因序列进行密码子优化后，由通用生物商业合成。采用OEPR克隆方法[27]，将BsXyn1和BsXyn2基因片段插入pF58载体，置于6×His标签序列下游，构建重组质粒pF58-BsXyn1和pF58-BsXyn2。以pF58-BsXyn1为模板进行定点突变，生成八个单点突变体：pF58-BsXyn1N63D、pF58-BsXyn1I77H、pF58-BsXyn1E78R、pF58-BsXyn1Y79I、pF58-BsXyn1Y80L、pF58-BsXyn1V81C、pF58-BsXyn1S100G和pF58-BsXyn1S134T。定点突变所用引物对列于补充表S1。
【数据】模板：pF58-BsXyn1；突变位点：N63D、I77H、E78R、Y79I、Y80L、V81C、S100G、S134T；数据库：UniProt P18429、Q3HLJ4
2.3 表达与纯化
将质粒转化至E. coli BL21(DE3)化学感受态细胞。转化细胞涂布于含100 μg/mL氨苄青霉素钠的LB琼脂培养基上，37°C过夜培养。次日挑取单菌落接种至4 mL含100 μg/mL氨苄青霉素钠的LB液体培养基中，37°C、220 rpm振荡活化。当培养物OD600达到0.6–0.8时，按1:1000（v/v）比例扩大培养至500 mL含100 μg/mL氨苄青霉素钠的LB液体培养基中，37°C、220 rpm培养3–4 h至OD600达0.6–0.8。最后加入IPTG至终浓度0.4 mmol/L诱导蛋白表达，16°C、220 rpm培养16 h。
IPTG诱导后，大肠杆菌细胞经5000 rpm离心15 min收集，弃上清。细胞沉淀按1:10（w/v）比例重悬于含0.5 mol/L NaCl和1 mmol/L PMSF的20 mmol/L Tris-HCl缓冲液（pH 8.0）中。使用超声波破碎仪（力辰，上海，中国）破碎细胞。裂解液于4°C、12000 rpm离心15 min，收集上清和沉淀组分。向上清中加入咪唑至终浓度10 mmol/L，用0.45-μm滤膜过滤。将澄清上清上样至预先用含10 mmol/L咪唑和0.5 mol/L NaCl的20 mmol/L Tris-HCl缓冲液（pH 8.0）平衡的Ni-NTA亲和层析柱。用3倍柱体积含30 mmol/L咪唑和0.5 mol/L NaCl的20 mmol/L Tris-HCl缓冲液（pH 8.0）洗涤后，用3倍柱体积含200 mmol/L咪唑和0.5 mol/L NaCl的20 mmol/L Tris-HCl缓冲液（pH 8.0）洗脱目标蛋白。洗脱蛋白用20 mmol/L PBS缓冲液（pH 7.5）透析脱盐，冻干后于−20°C保存。蛋白样品经SDS-PAGE电泳分析。
【数据】诱导剂：IPTG（终浓度0.4 mmol/L）；温度：37°C（活化）、16°C（诱导）；时间：3–4 h（扩大培养）、16 h（诱导）；转速：220 rpm；离心：5000 rpm/15 min、12000 rpm/15 min（4°C）；缓冲液：20 mmol/L Tris-HCl pH 8.0、0.5 mol/L NaCl、1 mmol/L PMSF；咪唑浓度：10/30/200 mmol/L；透析：20 mmol/L PBS pH 7.5；储存：−20°C
2.4 酶活测定
木聚糖酶活性采用3,5-二硝基水杨酸（DNS）法[28]测定，通过定量酶促反应中产生的还原糖（以木糖当量计）来计算。反应体系（总体积100 μL）包含90 μL溶于50 mmol/L磷酸-柠檬酸（McIlvaine）缓冲液（pH 6.0）的0.5%玉米芯木聚糖和10 μL适当稀释的酶液。反应混合物于50°C孵育10 min，然后加入100 μL DNS溶液终止反应。充分混匀后，沸水浴5 min，迅速冷却。随后加入300 μL去离子水混匀。取150 μL至微量板中，用酶标仪在540 nm处测定吸光度。基于木糖标准曲线定量木聚糖酶活性，以单位（U）表示，1 U定义为每分钟释放1 μmol还原糖所需的酶量。使用Bradford法测定蛋白浓度后计算比活。
【数据】底物：0.5%玉米芯木聚糖；缓冲液：50 mmol/L McIlvaine pH 6.0；温度：50°C；时间：10 min；检测波长：540 nm；酶活单位：1 U = 1 μmol还原糖/min
2.5 BsXyn1及其突变体底物特异性测定
底物特异性检测使用桦木木聚糖、玉米芯木聚糖、稻壳木聚糖、微晶纤维素、CMCNa、海藻酸钠、可溶性淀粉和壳聚糖。所有底物溶液配制成0.5%（w/v）浓度。分别测定BsXyn1、BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G对各底物的酶活。
【数据】底物浓度：0.5%（w/v）；底物种类：8种（桦木木聚糖、玉米芯木聚糖、稻壳木聚糖、微晶纤维素、CMCNa、海藻酸钠、可溶性淀粉、壳聚糖）
2.6 温度和pH对酶活的影响
木聚糖酶活性先在30°C至90°C温度范围内（间隔5°C）于50 mmol/L McIlvaine缓冲液（pH 6.0）中筛选，以确定最适温度。随后在50°C下用一系列pH值（4.0、5.0、5.5、6.0、6.0、7.0、8.0和9.0）确定最适pH，缓冲体系为50 mmol/L McIlvaine缓冲液（pH 4.0–8.0）和50 mmol/L甘氨酸-氢氧化钠缓冲液（pH 9.0）。两种测定中，观察到的最大木聚糖酶活性设为100%相对活性。所有测定重复三次。
【数据】温度范围：30–90°C（间隔5°C）；pH范围：4.0–9.0；缓冲液：50 mmol/L McIlvaine（pH 4.0–8.0）、50 mmol/L甘氨酸-NaOH（pH 9.0）；重复：三次
2.7 热稳定性测定
木聚糖酶热稳定性通过将酶在50 mmol/L McIlvaine缓冲液（pH 6.0）中于一系列温度（25°C、30°C、35°C、40°C、45°C、50°C、55°C和60°C）预孵育来测定。孵育0、0.5、1、2、4和6 h后，在50°C下测定酶活。相对活性以初始活性（0 h孵育）设为100%计算。所有测定重复三次。
【数据】预孵育温度：25–60°C（8个梯度）；孵育时间：0、0.5、1、2、4、6 h；测定温度：50°C；缓冲液：50 mmol/L McIlvaine pH 6.0；重复：三次
2.8 pH稳定性测定
为测定pH稳定性，将木聚糖酶稀释于50 mmol/L McIlvaine缓冲液中，pH范围5.0–8.0。在4°C孵育0、0.5、1、2、4和6 h后，于50°C测定酶活。相对活性以初始活性（0 h孵育）设为100%计算。所有测定重复三次。
【数据】pH范围：5.0、6.0、7.0、8.0；孵育温度：4°C；孵育时间：0、0.5、1、2、4、6 h；测定温度：50°C；重复：三次
2.9 EDTA和金属离子对酶活的影响
在终浓度2 mmol/L的EDTA或单一金属离子（Mg²⁺、Mn²⁺、Fe²⁺、Fe³⁺、Ni²⁺、Zn²⁺、Co²⁺、Cu²⁺、Ca²⁺和Cd²⁺）存在下评估木聚糖酶活性。相对活性以无添加剂时的酶活设为100%计算。所有测定重复三次。
【数据】添加剂浓度：2 mmol/L；金属离子：Mg²⁺、Mn²⁺、Fe²⁺、Fe³⁺、Ni²⁺、Zn²⁺、Co²⁺、Cu²⁺、Ca²⁺、Cd²⁺；重复：三次
2.10 有机溶剂对酶活的影响
木聚糖酶溶液与有机溶剂（甲醇、乙醇、正丁醇、异丙醇或乙酸乙酯）按50:50（v/v）体积比混合孵育。在25°C、200 rpm振荡孵育2 h。离心分离后，取水相于50°C、pH 6.0条件下测定木聚糖酶活性。相对活性以未处理酶液设为100%计算。
【数据】有机溶剂：甲醇、乙醇、正丁醇、异丙醇、乙酸乙酯；体积比：50:50（v/v）；温度：25°C；转速：200 rpm；时间：2 h；测定条件：50°C、pH 6.0
研究结果
3.1 两个枯草芽孢杆菌GH11木聚糖酶的比较分析
比较序列分析是定位控制木聚糖酶活性的关键残基、指导设计具有优良催化性能工程变体的有力工具[36]。通过系统筛选UniProt数据库，鉴定到两个来源于枯草芽孢杆菌的GH11木聚糖酶，命名为BsXyn1（UniProt ID: P18429）和BsXyn2（UniProt ID: Q3HLJ4）。BsXyn1最初克隆自枯草芽孢杆菌PAP115[39]，曾被用作定向进化以提高热稳定性的模板[37,38]。OmPraba等人在其引言中提出，来源于枯草芽孢杆菌AK1的BsXyn2在65°C和pH 8.5条件下表现出最适活性，并具有良好的热稳定性[40]。然而，除这一提及外，似乎缺乏系统研究和报道BsXyn2酶学特性的专门研究论文。
序列比对显示两个酶在以下位置存在八个氨基酸差异：63位（Asn到Asp）、77位（Ile到His）、78位（Glu到Arg）、79位（Tyr到Ile）、80位（Tyr到Leu）、81位（Val到Cys）、100位（Ser到Gly）和134位（Ser到Thr）（图1A）。值得注意的是，五个连续替换（残基77–81）发生在BsXyn1和其他芽孢杆菌GH11木聚糖酶高度保守的区域（补充图S1）。
因此，构建了BsXyn1和BsXyn2的重组表达载体，并在E. coli BL21(DE3)中表达蛋白。两种蛋白在相同条件下表达和纯化。SDS-PAGE分析显示BsXyn1以可溶形式高效表达（图1B），而BsXyn2主要以包涵体形式存在（图1C）。酶活测定表明，纯化的BsXyn1表现出高催化活性，比活为5969.33 ± 76.80 U/mg（图1D和表1），而BsXyn2因不可溶性未检测到活性。此外，研究者用BsXyn1对木四糖和木五糖进行了水解产物的TLC分析。如补充图S2所示，木四糖水解主要产生木二糖（X2），而木五糖同时产生木二糖（X2）和木三糖（X3），两种情况下均未检测到木糖。
为研究氨基酸变异的影响，通过定点突变构建了BsXyn1的八个单点突变体（BsXyn1N63D、BsXyn1I77H、BsXyn1E78R、BsXyn1Y79I、BsXyn1Y80L、BsXyn1V81C、BsXyn1S100G和BsXyn1S134T）。所有突变体在与BsXyn1和BsXyn2相同的条件下表达和纯化。它们均表现为可溶性表达，产量与野生型BsXyn1相当（补充图S3）。酶学分析显示，与野生型BsXyn1相比，三个突变体BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G对玉米芯木聚糖的催化效率显著提高，比活分别达7960.33 ± 128.29 U/mg（1.33倍）、14940.10 ± 939.43 U/mg（2.46倍）和8689.23 ± 446.66 U/mg（1.51倍）（图1D和表1）。SDS-PAGE分析确认纯化的BsXyn1及其突变体BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G纯度均大于95%（补充图S4）。圆二色光谱表明，BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G的二级结构相对于BsXyn1有轻微改变（补充图S5）。
值得注意的是，BsXyn1V81C突变体的活性增强最为显著。据研究者所知，BsXyn1V81C的比活远高于许多已表征的GH11木聚糖酶，包括来源于Bacillus pumilus HBP8（3986.7 U/mg）[9]、Neocallimastix patriciarum（3020.10 U/mg）[41]、Thermoanaerobacterium saccharolyticum（5622.0 U/mg）[42]和Paenibacillus campinasensis NTU-11（2750 U/mg）[9]的木聚糖酶。相反，其他突变体（BsXyn1N63D、BsXyn1I77H、BsXyn1E78R、BsXyn1Y80L和BsXyn1S134T）导致酶活显著降低（图1D）。值得注意的是，虽然BsXyn1E78R和BsXyn1Y80L表现出相对较高的表达产量（补充图S3），但其活性显著降低（图1D）。这种高表达与低活性之间的差异表明，这些突变可能影响折叠效率或内在催化能力。这些发现表明，79、81和100位残基对维持和增强GH11木聚糖酶的催化活性至关重要。
有益单点突变的组合诱变为进一步提高酶性能提供了途径。例如，来自Bacillus halodurans S7的木聚糖酶三突变体R69F/E137M/E145L比活较野生型提高了4.9倍，说明了多个突变的协同潜力[14]。据此，研究者构建了四个组合突变体（BsXyn1Y79I/V81C、BsXyn1Y79I/S100G、BsXyn1V81C/S100G和BsXyn1Y79I/V81C/S100G），并在与BsXyn1相同的条件下表达。与单点突变体的高溶解性和高产量相比，组合突变体表现出严重的表达受损（补充图S6）。只有BsXyn1V81C/S100G和BsXyn1Y79I/V81C/S100G以可溶形式表达（补充图S6）。虽然突变体BsXyn1V81C/S100G较野生型BsXyn1活性有所提高（6813.67 ± 223.48 U/mg，1.14倍）（补充图S7），但其活性仍显著低于任何一个单点突变体（BsXyn1Y79I、BsXyn1V81C或BsXyn1S100G）。因此，多个有益突变的组合可能导致潜在的拮抗相互作用，损害蛋白质生产和功能[43]。
GH11木聚糖酶以高底物特异性和可忽略的纤维素酶活性为特征，这与许多其他木聚糖酶不同。这一特性使其成为纸浆漂白应用的理想选择，可选择性地水解半纤维素而不损伤纤维素纤维[44]。因此，系统测定了BsXyn1及其有益突变体BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G对各种底物的催化活性（表1）。结果表明，所有突变体均保留了BsXyn1的木聚糖特异性催化能力，同时表现出显著增强的酶活。值得注意的是，突变体对不同来源木聚糖表现出不同的催化效率。BsXyn1V81C对玉米芯木聚糖（14,940.10 ± 939.43 U/mg）和桦木木聚糖（4504.90 ± 139.70 U/mg）的活性最高，较野生型分别提高150.23%和88.98%。相比之下，BsXyn1S100G在稻壳木聚糖上表现出最显著的活性（9312.47 ± 181.33 U/mg），较BsXyn1提高106.55%。底物特异性测定确认，所有测试酶严格靶向木聚糖底物，对羧甲基纤维素、微晶纤维素、可溶性淀粉、海藻酸钠或壳聚糖等非木聚糖底物均未检测到活性。这些发现表明，BsXyn1及其突变体具有高底物特异性，仅催化木聚糖水解，不表现纤维素酶、淀粉酶或其他糖苷水解酶活性。
【数据】BsXyn1比活：5969.33 ± 76.80 U/mg；BsXyn1Y79I比活：7960.33 ± 128.29 U/mg（1.33倍）；BsXyn1V81C比活：14940.10 ± 939.43 U/mg（2.46倍）；BsXyn1S100G比活：8689.23 ± 446.66 U/mg（1.51倍）；BsXyn1V81C/S100G比活：6813.67 ± 223.48 U/mg（1.14倍）；BsXyn1V81C对桦木木聚糖：4504.90 ± 139.70 U/mg（+88.98%）；BsXyn1S100G对稻壳木聚糖：9312.47 ± 181.33 U/mg（+106.55%）；纯度：>95%
3.2 温度和pH对BsXyn1及其有益突变体的影响
为表征野生型BsXyn1及其有益突变体BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G的酶学特性，首先测定了它们的最适反应温度（图2A和B）。四种酶均在50°C表现出最大活性，温度最适值相同（图2A）。然而，当温度超过55°C时酶活显著下降（图2A）。在60°C时，BsXyn1、BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G的残余活性分别降低了……
【数据】最适温度：50°C（四种酶一致）；活性显著下降温度：>55°C；60°C残余活性：原文未详述完整数值
（注：原文节选在此处中断，后续结果内容未包含在提供的原文节选中。根据图表清单，后续结果还包括热稳定性（图3）、pH稳定性（图4）、EDTA和金属离子影响（图5）、有机试剂影响（图6）、表面活性剂和还原剂影响（图7）、分子动力学模拟分析（图8）和残基能量分解分析（图9），但对应正文文本未在提供的节选范围内。）
讨论与解读
比较分析BsXyn1和BsXyn2鉴定出三个酶活显著增强的突变体BsXyn1Y79I、BsXyn1V81C和BsXyn1S100G。系统表征显示，尽管所有酶表现出相同的最适条件（50°C，pH 6.0），但所有突变体相对于BsXyn1热稳定性均有所降低。进一步的动力学和分子动力学模拟分析表明，thumb区域柔性的增加——通过升高的RMSF值和结合能分解来量化——是BsXyn1V81C催化效率提高的原因，其比活达14,940.10 ± 939.43 U/mg（较BsXyn1提高2.46倍），但伴随中等程度的稳定性折损。重要的是，所有突变体保持严格的木聚糖底物特异性和显著的pH稳定性。结合分子动力学模拟，结果表明增加thumb区域柔性以优化酶-底物相互作用往往伴随着活性-稳定性的权衡。更重要的是，本研究确定了thumb区域的动态特性是GH11木聚糖酶性能的关键决定因素。鉴于观察到的活性与稳定性之间的权衡，未来工作应优先考虑热稳定性优化，以使这些工程变体更适合工业应用。
编译者解读
本研究巧妙利用同源序列比较指导功能热点识别，避免了盲目随机突变，策略值得借鉴。V81C突变体14940 U/mg的比活数据在GH11家族中相当突出，但活性-稳定性权衡问题也提示：单一追求催化效率可能不适合高温工业场景。后续可考虑在V81C基础上引入稳定化突变（如二硫键或表面电荷优化）来补偿热稳定性损失。此外，组合突变体的表达受损现象提醒我们，多个"有益"突变叠加可能产生上位效应，理性设计需结合结构信息而非简单叠加。
DOI: 10.1016/j.ijbiomac.2026.151857