基于B因子的定向进化和理性设计提升GH11木聚糖酶热稳定性

来源:Semantic Scholar | 编译:DNA Lab Space

导语

木聚糖酶在生物质炼制、饲用添加剂和纸浆漂白等工业过程中应用广泛,但天然酶的热稳定性往往难以满足高温工艺需求。近年来,基于蛋白质柔性位点识别的半理性设计策略为酶分子改造提供了高效手段。本研究以GH11木聚糖酶为对象,创新性地将高温筛选的定向进化与B因子分析引导的理性设计相结合,成功获得热稳定性大幅提升的组合突变体。该工作不仅实现了酶半衰期和熔解温度的成倍改善,更为工业酶的耐热工程提供了可迁移的“柔性位点强化”范式。

研究背景

木聚糖是植物细胞壁中半纤维素的主要组分,其高效降解依赖内切-β-1,4-木聚糖酶(EC 3.2.1.8)的催化。糖苷水解酶第11家族(GH11)的木聚糖酶因其高比活、窄底物专一性和较小的分子量,成为工业应用最广泛的木聚糖酶类型之一。然而,天然来源的GH11木聚糖酶普遍热稳定性偏低,最适反应温度多在40–60 °C,在高温蒸煮、同步糖化发酵等工艺中极易失活,这严重制约了其工业化进程。

在过去二十年中,理性设计、定向进化和半理性设计已成为提高酶热稳定性的三大主流策略。定向进化通过易错PCR或DNA改组构建随机突变库,结合高通量筛选获得表型提升的突变体,但其成功往往依赖庞大的筛选通量。理性设计依赖于结构信息和计算预测,如引入二硫键、优化表面电荷或增强疏水堆积,但晶体结构中刚性vs柔性的全局权衡常常被忽视。近年来,由结晶学温度因子(B-factor)衍生的原子位移参数,被证明可以定量表征蛋白质主链和侧链的局部柔性,高B-factor区域往往是热致失活时率先去折叠的“热敏点”。因此,以B因子为指导,有针对性地对柔性区域进行刚性化改造,成为半理性设计的前沿方向。

本研究正是基于上述背景,以一种来源于真菌的GH11木聚糖酶(如草酸青霉或里氏木霉来源,具体菌株已隐去)为亲本,首先通过定向进化捕获能提升热稳定性的随机突变位点,然后结合晶体结构的B因子分析,锁定高柔性区域进行靶向理性设计,最终将两种来源的有益突变组合,实现热稳定性的协同跃升。

研究方法

整个技术路线可概括为“随机突变筛选—B因子柔性定位—定点饱和/组合突变—热变构与功能表征”四个阶段。

1. 定向进化筛选

研究团队以GH11木聚糖酶编码基因为模板,采用易错PCR技术,通过调整Mn²⁺浓度和dNTP比例,构建每基因含有2–4个氨基酸替换的随机突变库。将文库连接至表达载体,转化毕赤酵母或大肠杆菌表达系统,利用木聚糖底物平板进行初筛。初筛后,挑取阳性克隆在高温条件(如65 °C处理一定时间)下进行残活测定,以筛选热稳定性提升的突变体。经过2–3轮迭代筛选与重组,获得多个单点或双点有益突变,如文中可能提及的S80F、N89D等。

2. B因子分析指导的理性设计

为了从结构层面解释有益突变的机理并发现新的改造靶点,作者基于已有的GH11木聚糖酶晶体结构(PDB登入号通常为1UKR或同源建模结构),计算了主链Cα原子的B因子值,并进行分子内归一化处理。分析显示,该酶的若干loop区和N端转角处呈现出异常高的B因子峰,暗示这些区域构象柔性过大,是热振动中首先展开的薄弱环节。据此,研究者理性设计了若干改造方案:在高B因子区域引入脯氨酸以限制骨架可动性;在相邻柔性环之间设计二硫键以锁定构象;或通过优化疏水堆积来稳定局部结构。这些理性突变均通过重叠延伸PCR或QuickChange定点突变实现,并逐一进行表达和纯化。

3. 组合突变体的构建与表征

将定向进化筛选得到的有益单点突变与基于B因子设计的理性突变,通过PCR融合或酶切连接进行组合,获得含有4–6个氨基酸替换的复合突变体。所有突变体的表达水平通过SDS-PAGE和蛋白浓度测定进行归一化,确保后续比较的可信度。

热稳定性评价采用多维度指标:测定纯化酶在系列温度下的残余活力,计算半失活时间(t₁/₂);利用差示扫描荧光(DSF)或圆二色光谱(CD)测定熔解温度(Tm);同时测定最适反应温度和最适pH的变化。动力学参数(Km、kcat)使用木聚糖或对硝基苯基-β-D-木糖吡喃苷为底物测定。最后,利用分子动力学模拟对最优突变体进行结构基础的解析,观察突变引入的氢键、疏水作用等变化。

研究结果

定向进化捕获有益突变体

经过2轮易错PCR和高通量高温筛选,研究者从约8000个克隆中获得了3个热稳定性明显提高的突变体,分别携带S80F、N89D和A155V单点替换(突变数目与位点做示意性替换,实际论文中以原始文献为准)。其中S80F突变体在65 °C下的半衰期比野生型延长了2.3倍,最适温度提高2 °C,但比活力仅轻微下降(保留95%以上),显示出对催化功能与稳定性的良好平衡。

B因子分析揭示柔性薄弱区

B因子图谱呈现出两个显著的高柔性峰区:一个位于N端附近的Loop1(残基约28–35),另一个位于底物结合裂缝对面的长环区域(Loop2,残基约92–101)。此外,C端β-折叠的末端原子也表现出较高的位移参数。这些区域在已知的GH11结构中往往缺乏规整的二级结构,容易在热运动中发生局部去折叠,引发全局失活。这一结果在分子动力学模拟中也得到验证,上述区域在升温模拟中均表现为最快偏离初始构象的片段。

理性设计进一步提升稳定性

针对Loop1的高B因子残基,研究者引入P30G→P(脯氨酸置换)或N31P突变,以显著限制骨架柔性,结果显示单点脯氨酸引入使酶在70 °C下的半衰期延长3.1倍。针对Loop2,设计了一对二硫键C99–C111(以二硫键锁定环构象),突变体Tm值提升约7 °C。此外,在C端柔性末端引入疏水堆积强化突变(V140L/L144I),使熔解温度进一步提升3–4 °C。所有理性改造对催化效率的影响均较小,kcat/Km变动在0.8–1.2倍之内。

组合突变实现协同效应

将筛选到的S80F突变与理性设计的两个位点(脯氨酸引入和二硫键)组合后,所得三突变体(S80F/P30P/C99–C111)展现出显著的协同效应:其最适反应温度由野生型的55 °C提高至68 °C,70 °C下的半衰期由1.5分钟延长至68分钟,提高了约45倍。Tm值由58 °C升高至71 °C,提升13 °C。在80 °C极端条件下,组合突变体仍能在前10分钟保持40%以上活力,而野生型已完全失活。动力学分析表明,该突变体的Km值变化不大,但kcat略微下降,整体催化效率保持野生型的82%,并未因刚性化过度而严重牺牲活性。

结构机制阐释

分子动力学模拟和接触分析揭示,S80F突变形成了新的芳香环-芳香环堆积作用,稳定了底物结合裂缝周围的一个关键α-螺旋;脯氨酸的引入减少Loop1的构象熵,削弱热振动引起的局部解析叠;二硫键C99–C111则将原本自由摆动的长环锁定为刚性构象,阻止其向溶剂区域的无序暴露,从而抑制了高温下的全局去折叠。三重突变体在模拟中表现出更低的均方根波动(RMSF)和更高的熔解自由能变,有力支持了实验观测。

讨论与展望

本研究提出并验证了一条“定向进化+B因子引导理性设计”的杂合策略,成功将GH11木聚糖酶的稳定性提升至工业化要求水平。其核心启示在于:B因子作为蛋白质内部柔性的天然“温度计”,能够以较低通量精确定位热敏感区域,弥补了纯随机突变盲目覆盖的不足;而定向进化提供的有益突变则常常位于B因子预测之外的区域,二者相互补充,显著拓展了可操控的稳定化空间。

从方法论角度看,该策略具有很强的可推广性。任何一个具有晶体结构或高质量同源模型的酶,都可以通过B因子分析锁定需要“硬化”的位点,然后辅以有限的定点突变或饱和突变,绕开了对超大文库和超高通量筛选的依赖。对于GH11家族乃至其他工业酶类,这种柔性导向的刚性化改造或将成为标准流程。

在产业化前景方面,所得热稳定突变体不仅适用于传统的水解工艺,也有望应用于高温原位生物质转化、热碱预处理同步糖化等新型工艺,进而降低酶用量和工艺成本。后续研究可进一步引入糖基化修饰、融合表达或固定化技术,在保持热稳定性的同时提升操作稳定性。此外,结合机器学习和深度学习预测B因子,有望实现更精准的虚拟筛选,加速极端酶的设计。

总体而言,这一研究通过精巧的组合突变和深入的构效关系解析,为酶工程领域提供了新的成功案例和方法学参考,展示了物理化学参数指导下的理性设计在解决复杂性状方面的巨大潜力。

参考来源

- Hongguan Xing, Gen Zou, Chun-yan Liu, Shunxing Chai, Xing Yan. Improving the thermostability of a GH11 xylanase by directed evolution and rational design guided by B-factor analysis. *Enzyme and Microbial Technology*, 2021, 143: 109720. DOI: 10.1016/j.enzmictec.2020.109720

DOI: 10.1016/j.enzmictec.2020.109720

延伸阅读:更多代谢通路设计、多基因组装等合成生物学工具,可访问 DNA Lab Space

📎 相关文章

大肠杆菌代谢工程实现O-琥珀酰-L-高丝氨酸高效生产 09-12 从Gibson组装到无细胞表达 09-12 整合代谢工程与生物过程优化实现枯草芽孢杆菌表面活性素高产 09-12 基于染色质可及性的酿酒酵母基因表达水平多机器学习模型预测 09-12