《Process Biochemistry》：Charge-Dependent Interfacial Mechanism and Conformational Dynamics of Lipase Immobilization on ZIF-8: Insights from MD Simulation and Experiments
编辑推荐：
脂酶表面电荷差异调控金属有机框架（ZIF-8）吸附机制及催化性能研究：通过分子动力学模拟与实验验证，发现负电性ROL与ZIF-8形成疏水-静电协同作用，活性提升131%；正电性RNL因静电排斥导致构象扭曲及活性下降58%。揭示了表面电荷影响酶-载体界面吸附及构象稳定性的分子机制，为工业生物催化固定化系统设计提供理论依据。
梁欣|黄一萌|朱旭成|Prasanna J Patil|李伟伟|李秀婷|孙宝国
教育部北京工商大学老年营养与健康重点实验室，北京工商大学，北京100048，中国
摘要
在金属有机框架（MOFs）上固定化过程中保持脂肪酶活性仍然是一个主要挑战。本研究系统地探讨了疏水性且易于合成的沸石咪唑框架-8（ZIF-8）如何通过原子级别的界面工程影响具有不同表面电荷的脂肪酶。分子动力学模拟显示，带负电荷的米曲霉脂肪酶（ROL）与ZIF-8的甲基咪唑基团形成协同的疏水性和静电相互作用，并构建稳定的氢键网络，使酶保持开放构象，从而提高催化效率（达到游离酶的131%）。相反，带正电荷的米平曲霉脂肪酶（RNL）与ZIF-8发生静电排斥，导致构象重排、盖子闭合以及催化位点变形，仅保留了58%的活性。结构分析、ζ电位测定和吸附动力学研究证实了这些机制。这项工作为脂肪酶与MOFs之间的相互作用提供了分子层面的见解，并为工业生物催化中高性能固定化酶系统的合理设计提供了指导。
引言
脂肪酶是工业生物催化领域的重要工具酶，应用于食品、乳化剂、生物柴油、香料、药品、化妆品、废水处理等领域，因为它们能够催化多种反应，如酯水解、酯化、酯交换、酸解和醇解[1]、[2]、[3]。然而，脂肪酶在催化过程中的构象适应性伴随着内在的不稳定性和对环境扰动的敏感性增加，这可能会破坏内部氢键并改变活性位点的可及性，从而降低催化效率[4]。作为关键的生物催化工程策略，酶的固定化可以通过空间限制来调节酶的微环境，从而提高结构稳定性、拓宽工作范围、增强对环境变化的适应性并延长保质期[2]。
金属有机框架（MOFs）是由金属节点与有机配体自组装形成的多孔晶体材料，由于其高比表面积、孔径可调性和表面功能化能力，在极端操作条件下构建酶保护系统方面表现出独特优势[5]、[6]。在各种MOFs中，沸石咪唑框架（ZIFs）因其优异的化学稳定性、温和的合成条件和良好的生物相容性而特别受到青睐。在过去十年中，ZIF-8已成为研究最广泛的用于脂肪酶固定化的材料之一，具有优异的疏水性、低毒性和易于合成等优点[7]、[8]、[9]、[10]。尽管通过载体修饰和优化固定化策略，MOF工程显著提高了脂肪酶的性能[2]、[6]、[11]、[12]、[13]，但仍经常报告活性损失，表明通过MOF工程提高性能的机制尚未完全理解（见表S1）[14]、[15]、[16]。核心挑战在于酶-载体界面的机制尚未完全阐明。大多数研究主要集中在载体的制备、固定化酶的催化性能和应用等方面[17]、[18]、[19]。然而，活性增强或失活的微观起源，特别是界面微环境调节对脂肪酶构象动态的影响，尚未得到深入探索[20]。
对于脂肪酶而言，催化效率与控制底物通过界面激活机制进入活性位点的盖子区域的动态开闭密切相关[21]。MOF表面性质（如电荷分布和疏水性）如何调节盖子动态和催化中心几何结构仍知之甚少，这成为合理设计高性能固定化脂肪酶系统的关键障碍。分子动力学（MD）模拟作为研究生物分子和纳米材料界面行为的核心计算工具，能够提供关于结合方式、相互作用力和构象转变的原子级别见解，这些在实验中难以捕捉[22]、[23]、[24]、[25]。近年来，MD模拟为ZIF-8-脂肪酶相互作用提供了重要见解。Xia等人通过MD模拟表明，黑曲霉脂肪酶（AOL）与ZIF-8的结合不仅依赖于范德华力（vdW），还涉及静电相互作用的协同效应[26]。He等人进一步发现，ZIF-8的表面电荷分布通过静电-疏水性协同作用调节了光滑假丝酵母脂肪酶（CRL）的吸附取向和构象稳定性，从而影响催化性能[20]。然而，这些研究主要集中在吸附能量和整体结构稳定性上，而对盖子动态、催化中心几何结构和脂肪酶表面电荷异质性的系统研究仍然有限。值得注意的是，脂肪酶家族成员之间的内在结构异质性，如活性位点拓扑结构、表面电荷分布和灵活区域的差异，可能导致它们与MOFs的特定结合模式，进而产生不同的催化性能[27]。目前仍缺乏对这种异质性如何影响脂肪酶-MOF兼容性的全面理解。
基于此，本研究通过结合MD模拟和实验验证，系统地揭示了具有显著不同电荷特性的米曲霉脂肪酶（ROL，pI 6.40）和米平曲霉脂肪酶（RNL，pI 8.14）在ZIF-8上的固定化系统的结构-活性关系和分子机制。该研究不仅阐明了脂肪酶表面电荷分布对ZIF-8结合模式的调节作用，还建立了固定化过程中不同pI值脂肪酶的构象动态与催化活性之间的关系。这为未来通过计算机模拟筛选合适的脂肪酶生物相容性载体提供了理论指导，并为高性能固定化脂肪酶系统的合理设计奠定了科学基础。
部分摘录
化学物质和脂肪酶
ROL（产品代码80612）和RNL（产品代码62310）购自日本Sigma-Aldrich公司。2-甲基咪唑（纯度≥99%）和乙酸锌（纯度≥99.99%）购自美国Merck公司。4-硝基苯基丁酸（p-NPB，纯度≥98%）购自Sigma-Aldrich公司。用于蛋白质检测的双缩肌酸（BCA）试剂购自中国Solarbio公司。本研究中的所有其他化学物质均为分析级。MD模拟
ZIF-8的结构（CCDC ID：1579412）数据来源于剑桥数据库
用于MD模拟的固定化脂肪酶的整体结构分析
蛋白质主链的均方根偏差（RMSD）在最初的0-100 ns内达到稳定平台，并在此后保持稳定，表明系统已达到动态平衡。关键结构参数，包括回转半径（Rg）和溶剂可及表面积（SASA），在模拟的最后100 ns内也没有系统性变化（图1b,c,e,f）。此外，ROL-ZIF-8和RNL-ZIF-8都在50 ns内快速完成了结合过程结论
在本研究中，我们通过分子动力学模拟结合实验验证，揭示了米曲霉脂肪酶（ROL）和米平曲霉脂肪酶（RNL）与ZIF-8之间的差异化界面相互作用机制及其构象关系。结果表明，表面富含负电荷残基的ROL通过与ZIF-8的疏水-静电协同作用形成了多个稳定的键合
CRediT作者贡献声明
朱旭成：撰写、审稿和编辑。李伟伟：项目管理。Prasanna J Patil：撰写、审稿和编辑。李秀婷：撰写、审稿和编辑、监督、资源获取、资金申请。孙宝国：监督。梁欣：撰写、初稿撰写、可视化、软件应用、方法学设计、实验研究、数据分析、概念构建。黄一萌：撰写、初稿撰写和验证。致谢
作者们对中国国家自然科学基金（项目编号32372280）和中国国家重点研发计划（项目编号2022YFD2101401）的支持表示衷心的感谢。同时，我们也非常感谢Lu Hongyun在分子动力学模拟方面提供的帮助。
