编辑推荐：
分子动力学模拟为研究木质纤维素中酶与底物的相互作用提供了高效方法，涵盖模型构建、软件与力场比较、动态机制解析及工业化应用优化。
王思佳|陈宏|陈子健|兰天青|傅向进|林秦璐
中南林业科技大学食品科学与工程学院，中国长沙天心区，410004
摘要
木质纤维素是地球上最丰富的可再生碳储存库，对可持续生物能源和生物材料的生产具有至关重要的意义。木质纤维素由三种聚合物组成——纤维素、半纤维素和木质素——它们的层次结构和分子间相互作用从根本上决定了木质纤维素生物质的酶法糖化效率。分子动力学（MD）模拟提供了一种有效的方法来探索酶与木质纤维素之间的相互作用，并且近年来已成为该领域广泛采用的研究方法。本综述系统地概述了当前的MD模拟方法，包括对常用软件平台和力场的比较分析，评估了它们在模拟木质纤维素系统方面的优缺点。随后介绍了纤维素、木质素和纤维素酶的基本结构及其建模策略。最后，总结了MD模拟在酶-木质纤维素相互作用中的应用，包括纤维素酶对纤维素的识别、吸附和水解作用，木质素对纤维素酶活性的影响，木质素与纤维素之间的相互作用，以及半纤维素与其他成分之间的相互作用。通过将方法论进展与代表性案例研究相结合，本综述为设计、验证和解释木质纤维素酶法水解的MD模拟提供了实用指导，最终目标是实现生物质分解过程的合理改进。
引言
木质纤维素是一种天然的大分子复合材料，主要由纤维素（占总质量的30–50%）、半纤维素（20–30%）和木质素（15–30%）组成。作为自然界中最丰富的可再生生物质资源，木质纤维素具有普遍可用性、低成本、可生物降解性和高可再生性等显著优势[1]，[2]，使其成为可持续生物能源生产的理想原料。
木质纤维素的酶法水解是将复杂的纤维素和半纤维素结构转化为纤维素乙醇的关键步骤[3]，这一过程的优化直接决定了木质纤维素资源转化的效率。因此，深入研究纤维素与纤维素酶之间的界面相互作用至关重要。该水解过程主要分为三个步骤：酶在底物上的吸附、酶的催化水解以及酶从底物上的脱附[4]。尽管大量的实验研究已经能够对酶促过程进行宏观表征，但目前的技术限制阻碍了对木质纤维素组分与纤维素酶之间原子级相互作用动态的全面阐明。一些关键问题尚未解决：首先，初始的酶-底物接触和吸附过程作为水解的启动步骤，在动力学参数和底物表面性质调控方面仍知之甚少。其次，活性位点与底物结合的分子识别机制涉及多种力的复杂相互作用，包括氢键网络、疏水相互作用、静电作用和范德华力，所有这些都需要在原子层面进行表征。此外，纤维素酶在催化过程中的动态构象变化，特别是那些调节底物可及性和催化三联体定位的变化，直接决定了水解效率，但在实验上仍难以捕捉。这些动态构象变化的分子机制需要进一步研究。最后，酶的方向、底物表面形态以及界面水如何共同影响结合亲和力、催化转化率和酶稳定性仍然没有明确答案。
近年来，随着木质纤维素分解研究的日益复杂化以及计算能力的提升，MD模拟技术已成为探索木质纤维素组分与酶之间相互作用的首选工具。MD模拟能够在原子尺度上对木质纤维素的酶法水解过程进行理论建模和计算，其结果与实验数据高度一致。这种计算方法在时间效率和成本效益方面明显优于传统的实验方法。具体来说，它有助于优化木质纤维素的分解效率，加深了对酶-底物识别模式的机制理解，并为合理设计纤维素酶以减少非生产性结合事件提供了依据。同时，MD模拟可以提供酶与底物之间相互作用的详细视图，包括原子间的接触、氢键的形成和断裂，以及疏水力和电荷相互作用的模式[5]，[6]。此外，MD模拟可以预测酶-底物相互作用的可能结果，为构建新的酶结构或添加新组分以提高水解效率的研究提供大量理论指导，这种预测能力显著减少了实验中的试错次数，提高了研究效率。
MD模拟的多尺度建模能力也是选择它的原因之一。结合量子力学计算、经典牛顿动力学和粗粒化模型，MD模拟可以在不同尺度上分析酶法水解过程，提供更全面的相互作用视角。在Kou等人的研究中[7]，系统探讨了木质素-碳水化合物复合物（LCCs）中组分大小对界面相互作用的影响。利用MD模拟，研究表明半纤维素的结构会影响木材的性质和组分之间的相互作用。虽然这项工作提供了对静态吸附机制的详细见解，但缺乏对LCCs长期动态演变的全面分析。文章侧重于静态吸附行为的分析，而对长期动态行为的研究不足。此外，对生物大分子相互作用的探索仍局限于特定类别，如纤维素-木质素界面，对更广泛的分子多样性的考虑有限。同时，Vermaas等人[8]利用MD模拟阐明了木质素的物理化学性质——特别是表面疏水性和官能团分布——在调控酶法水解效率中的作用。他们的发现表明，木质素的分子结构直接影响纤维素酶的可及性和非生产性结合倾向，从而调节糖化动力学。然而，该研究未能阐明木质素大分子与纤维素酶在酶法水解过程中氢键网络的动态形成-解离机制，也未阐明这些机制如何调节酶的构象转变和活性位点的可及性。在Zhao等人的综合综述[9]中，详细总结了木质素-酶相互作用研究的最新进展，特别强调了用于机制阐明的先进分析技术。然而，他们对MD模拟数据的分析仍然较为表面化。具体来说，综述忽略了木质素结构变化如何调控纤维素酶吸附动力学和催化能力的分子层面探究。此外，文章缺乏对实际酶法水解过程中木质素与酶之间相互作用动态变化的深入分析，例如木质素解聚引起的结合位点暴露或其对酶过程性的反馈效应。为了填补这些知识空白并充分发挥计算工具的潜力，本综述采用了以方法论为重点的视角，旨在为应用MD模拟解析木质纤维素酶法水解过程中的复杂相互作用提供一个系统框架。我们不仅列举了生物学见解，还强调了从模型构建和力场选择到模拟的整个流程。通过批判性地比较不同的MD方法并评估它们在木质纤维素-酶研究中的表现，本综述为研究人员选择、实施和推进该领域的模拟方法提供了实用指导。
基于上述讨论，本综述主要讨论了MD模拟在木质纤维素水解相互作用方面的一系列进展和未来发展。内容包括MD软件和力场的选择、纤维素、木质素和纤维素酶的分子建模，以及在探测酶-底物相互作用方面的实际应用（图1）。综述不仅总结了MD模拟在木质素、纤维素和纤维素酶相互作用研究中的当前应用现状，还思考并展望了未来如何更好地利用MD模拟方法来研究木质纤维素酶法水解的机制。
MD模拟概述
MD模拟基于牛顿力学，利用计算建模来追踪随时间变化的粒子轨迹，揭示微观动力学与宏观材料性质之间的基本关系。这种方法能够定量预测宏观物理性质，同时提供微观结构配置和粒子行为的原子级分辨率。该方法最初由Alder[10]提出，使用理想的“硬球”模型来研究方程
建立模型
建立准确的MD模型是成功模拟结果的基本前提，构建能够捕捉木质纤维素层次结构的真实模型仍是一个关键挑战。木质纤维素的复杂结构和酶的动态特性使得模型的构建和验证极为复杂。纤维素的结晶度、木质素的单体分布以及酶的灵活区域都影响着模型的构建
利用分子动力学模拟探索木质素、纤维素、半纤维素和纤维素酶之间的相互作用
木质纤维素酶法水解中的反应主要发生在四种组分之间（图6），影响这一过程的重要因素包括疏水相互作用、静电相互作用、氢键（包括π-π堆叠）和范德华力。疏水相互作用主导了纤维素酶在木质纤维素表面的初始吸附，而静电和氢键网络则调节结合的特异性和
结论与展望
MD模拟已成为解析木质纤维素-酶相互作用微观机制不可或缺的工具，提供了原子尺度上的独特时空分辨率，为生物质分解的物理化学性质提供了关键见解。本综述系统总结了该领域的关键进展，重点关注从模拟研究中获得的关于纤维素酶降解动力学和底物识别的深刻见解
作者贡献声明
王思佳：撰写——审稿与编辑，撰写——初稿。陈宏：撰写——审稿与编辑。陈子健：撰写——审稿与编辑。兰天青：撰写——审稿与编辑，概念化。傅向进：撰写——审稿与编辑，监督。林秦璐：监督。
利益冲突声明
作者声明他们没有已知的财务利益或个人关系可能影响本文报告的工作。
致谢
本工作得到了中国湖南省自然科学基金（资助编号：2023JJ30985）的支持。
