《Protein Science》：Distinct roles of monosaccharides in the inhibition of the carbohydrate-binding module and catalytic domain of Cel7A
编辑推荐：
木质纤维素生物质的酶促降解是生物化学品可持续生产的关键步骤。纤维二糖水解酶Cel7A通过水解纤维素在该过程中发挥核心作用。其主要反应产物（纤维二糖）和次级单糖（葡萄糖、甘露糖和木糖）通过与Cel7A的碳水化合物结合模块（CBM）和催化域（CD）相互作用抑制其活
木质纤维素生物质的酶促降解是生物化学品可持续生产的关键步骤。纤维二糖水解酶Cel7A通过水解纤维素在该过程中发挥核心作用。其主要反应产物（纤维二糖）和次级单糖（葡萄糖、甘露糖和木糖）通过与Cel7A的碳水化合物结合模块（CBM）和催化域（CD）相互作用抑制其活性。前期工作已证实，葡萄糖直接降低酶的催化速率（kcat），而甘露糖和木糖主要降低其结合速率（kon）和整体持续加工性。为揭示这些观察结果背后的微观机制，研究人员采用全原子分子动力学模拟结合Kirkwood–Buff溶液理论。结果表明，两个域的溶剂化偏好存在意想不到的差异：糖类优先从CD中被排除，而CBM则表现出优先的糖结合。研究人员观察到沿CBM平面结合面（残基Y5、Y31和Y32）以及I11、S14、S21、V27、L28和P30附近存在高密度的糖。木糖作为实验中最有效的抑制剂，在构成纤维素结合面的CBM残基附近显示出最高密度，这可能干扰酶对纤维素的锚定。此外，研究人员观察到糖在CD隧道入口处及隧道结合环（B2、B3和A3）内积累，为酶持续加工循环受阻提供了结构解释。通过协调更广泛的溶剂化偏好与高度特异的位点相互作用，研究结果为葡萄糖和半纤维素来源糖类如何抑制Cel7A提供了图景。
**研究背景与意义**
木质纤维素生物质是可持续生产生物燃料和工业生物产品的丰富可再生资源，主要由纤维素（30%–45%）、半纤维素（20%–40%）和木质素（10%–25%）组成。纤维素是由β-1,4-连接的D-葡萄糖单元构成的线性多糖，是发酵糖的主要来源。将结晶纤维素转化为可发酵糖依赖于物理化学预处理与多组分酶系统的双重策略。然而，生物质降解过程中产生的产物——包括纤维二糖、葡萄糖、甘露糖、木糖及各类酚类化合物——是已知的酶抑制剂。即使经过严格预处理，残留组分仍常吸附于纤维素表面，共同构成复杂的抑制环境，严重制约酶系统效率。
在生物质降解酶中，Cel7A是丝状真菌Trichoderma reesei（或Hypocrea jecorina）分泌的主要纤维二糖水解酶之一，通常占总蛋白分泌谱的一半以上。Cel7A具有模块化结构，由碳水化合物结合模块（CBM）通过富含脯氨酸、丝氨酸和苏氨酸的内在无序糖基化连接肽与催化域（CD）相连。CD具有长隧道状底物结合裂隙和八个主要环（A1–A4和B1–B4），覆盖由β-折叠三明治形成的裂缝。TrCel7A的CD具有约50 ?长的底物结合位点，结合纤维素链并从还原端水解糖苷键产生纤维二糖。其CBM主要通过三个酪氨酸（Y5、Y31、Y32）促进酶在纤维素表面的吸附，这些酪氨酸通过疏水相互作用、CH–π和π–π堆积排列在结晶纤维素的平坦表面上。
尽管已有多种实验技术报道了糖类对Cel7A抑制的宏观影响，但CBM和CD差异化作用的分子层面解释仍难以获得。Nong等近期提出纤维二糖通过高亲和力结合CD隧道中的"前门"（W40）和"后门"（W376）位点，降低酶的着陆速率和持续移动。除纤维二糖外，葡萄糖、甘露糖和木糖等单糖也具有抑制作用。值得注意的是，这些单糖并非以相同方式破坏酶的功能：葡萄糖直接损害酶的核心催化活性，而甘露糖和木糖则主要降低酶在纤维素表面有效结合和移动的能力。Jiang等近期的动力学模型强有力地支持了这一机制区分，其工作表明甘露糖和木糖并非影响实际催化速率（kcat），而是通过有效降低酶的结合速率（kon）和整体持续加工性来驱动抑制。
为揭示这些抑制背后的微观机制，研究人员对酶（CBM和CD两个域）在不同浓度的葡萄糖、甘露糖和木糖溶液中进行了分子动力学（MD）模拟，采用最小距离分布函数（MDDFs）和Kirkwood–Buff（KB）溶剂化理论表征溶剂化结构。研究结果揭示了一个显著偏好：糖类优先从CD中被排除，但相对于水优先结合于CBM。这种偏好尤其在构成酶平坦纤维素结合面的特定保守酪氨酸残基（Y5、Y31和Y32）周围最为显著。历史上，CBM对抑制的敏感性因其对结晶纤维素的天然高亲和力而被低估。然而，模拟表明这些关键酪氨酸与周围糖类——尤其是木糖——之间的相互作用可能改变控制初始底物接触的局部热力学，为先前观察到的酶着陆速率（kon）降低提供了分子基础。同时，研究人员在CD处观察到相反的模式：虽然CD的广泛表面通常排除这些糖类，但在关键结构位点仍发生特异性相互作用。具体而言，研究人员注意到隧道入口处及隧道包围环（B2、B3和A3）周围存在特异性和非特异性相互作用。这些局部糖相互作用揭示了催化隧道内可能与底物结合竞争的高亲和力位点，并可能阻碍沿纤维素链的持续移动。综合而言，研究结果提供了一个将域特异性溶剂化偏好与残基级结构相互作用联系起来的坚实热力学框架，为这些糖在持续水解过程中如何抑制Cel7A提供了分子层面的解释。
**关键技术方法**
研究人员采用全原子分子动力学（MD）模拟结合Kirkwood–Buff（KB）溶液理论，对TrCel7A的CD（PDB ID: 8CEL）和CBM（PDB ID: 2CBH）在葡萄糖、甘露糖和木糖的α-和β-异头物溶液中（0.1–2.0 mol L?1，八个浓度）进行模拟。CBM每个系统运行20个独立副本×80 ns，CD每个系统5个副本×1000 ns。采用CHARMM36m力场和TIP3P水模型，糖分子使用ADD力场参数化。通过最小距离分布函数（MDDFs）和Kirkwood–Buff积分（KBIs）量化局部和全局溶剂化模式，计算优先相互作用参数（Γ）。此外，利用STRIDE算法评估二级结构，RMSF评估局部柔性，RMSD表征B3环构象变化，WASA量化溶剂可及表面积。
**研究结果**
**3.1 Cel7A CBM在糖溶液中的溶剂化性质**
研究人员首先聚焦糖类及其浓度对Cel7A CBM溶剂化性质的影响。MDDFs显示两个明确的溶剂化壳层，分别对应氢键相互作用和非特异性及水介导的相互作用。随葡萄糖浓度增加，水的MDDF峰略升高，而葡萄糖的MDDFs降低。KBIs表明CBM的排除体积被水在2.0 ? < r < 4.0 ?处的局部密度增强部分补偿。在较大距离（r > 7.0 ?）处，水的KBI与葡萄糖浓度无关，收敛于约?3.0 L mol?1。相反，葡萄糖的KBIs表现出对共溶剂浓度的显著依赖性，最低浓度（0.1 M）下对CBM的亲和力最高。优先相互作用参数（Γglucose）揭示CBM在高达1.5 mol L?1时被葡萄糖优先溶剂化，在0.5 mol L?1处达到最大值，在2.0 mol L?1处发生向优先水合的转变。葡萄糖积累的详细图谱显示，葡萄糖主要与纤维素结合面上的芳香族残基（Y5、Y31和Y32）以及疏水和极性残基（包括I11、S14、S21、V27、L28和P30）相互作用，这些相互作用主要为非特异性的。
甘露糖和木糖对CBM的溶剂化遵循与葡萄糖相似的趋势，但存在特殊性。甘露糖和木糖的MDDFs也显示两个溶剂化壳层，但与葡萄糖相比，其第一溶剂化壳层密度较小，而第二溶剂化壳层密度较大。木糖在第二溶剂化壳层的增加最为显著。将总MDDF分解为极性和非极性残基贡献表明，极性残基对所有三种糖的总CBM溶剂化提供最大绝对贡献，这主要来自糖羟基与亲水蛋白表面之间广泛的氢键网络。然而，糖与极性和非极性残基之间的非特异性相互作用从葡萄糖到木糖增加，在酪氨酸MDDFs的r ~ 2.5 ?处尤为明显，表明木糖增强的溶剂化由其与芳香三联体进行疏水相互作用的更高倾向驱动。差异密度图揭示甘露糖尤其是木糖在芳香三联体（Y5、Y31和Y32）及相邻疏水残基（I11、V27、L28和P30）的第二溶剂化壳层内积累增加，木糖显示三种共溶剂中最高的局部密度。这一趋势与三种糖的相对疏水性平行（木糖 > 甘露糖 > 葡萄糖），结构上与木糖缺乏葡萄糖和甘露糖中存在的C6羟甲基一致。木糖的KBIs始终保持较高，证实其对CBM表面的更强亲和力。值得注意的是，虽然Γxylose在2.0 mol L?1处降低，但仍为正值，表明木糖抵抗葡萄糖和甘露糖所观察到的向优先水合的转变。
α-异头物显示更高的整体优先相互作用参数，尽管异头物形式在各浓度下表现出相似趋势。基于水溶液平衡丰度的加权平均优先相互作用参数严格保持相对顺序（木糖 > 甘露糖 > 葡萄糖）。残基逐一分析Γ贡献表明，β-D-木糖不仅与更多残基优先相互作用，而且这些相互作用也明显更强。三种糖的主要个体贡献者始终是纤维素结合面上的芳香族残基（Y5、Y31或Y32）。β-D-木糖的独特之处在于更广泛的正面贡献分布于显示增强相互作用的侧链附近残基，这种更广泛的优先相互作用范围维持了整体正Γ值。在低浓度下，这些芳香位点被有效占据，更有效地贡献正优先相互作用。然而在2.0 mol L?1处，虽然相互作用位点保持一致，但这些残基周围糖的相对密度降低，表明位点部分饱和。值得注意的是，Y13（不属于保守结合面的残基）周围的糖密度保持较低。由于这三个核心酪氨酸对介导与结晶纤维素表面的CH–π堆积相互作用至关重要，其被游离糖分子的持续占据可能遮蔽CBM的主要结合界面，这为酶抑制在结晶纤维素上比无定形底物上更显著提供了合理解释。
**3.2 Cel7A CD在水和糖溶液中的结构性质**
与CBM相对平坦的结合表面不同，CD具有复杂的隧道状结构，由一系列柔性环（A1–A4和B1–B4）包围，共同形成催化隧道的顶部。在无底物状态（apo-CD）下，A和B环及N端和C端区域表现出天然柔性，A1、A2、A4和B4环的高RMSFs尤为显著。模拟揭示糖的存在改变环动力学。A2环（残基Q98至K102）表现出局部α-螺旋含量，通常被糖溶液稳定。B3环（G244–G253）在所有浓度下糖存在时RMSF系统性降低，其RMSD分布向较低值偏移。先前的主成分分析（PCA）已识别apo-CD执行协同"呼吸"运动，其特征在于隧道包围环的开合波动对底物加工至关重要，而糖存在下环的较小移动性可抑制这种运动。B3环对维持TrCel7A在结晶纤维素上的高持续加工性很重要，通过催化后稳定纤维素链还原端于产物结合亚位点（+1/+2）。糖诱导的B3柔性限制可能调节产物释放或纤维素链的前向推进。尽管存在这些局部动态变化，CD的整体折叠不受糖存在影响。
在TrCel7A CD的底物结合隧道内，四个保守色氨酸残基（W40、W38、W367和W376）通过在?7（W40）、?4（W38）、?2（W367）和+1（W376）亚位点提供芳香族堆积相互作用促进持续水解。研究人员观察到糖添加后几乎整个酶表面的水可及性（WASA）系统性降低，对这些特定芳香族残基的影响尤为显著。即使在低糖浓度下，这些色氨酸的WASA降低超过50%，尽管其结构移动性基本不变。这种显著降低表明糖有效遮蔽芳香族残基的水合，与纤维素链疏水堆积能力降低一致，为底物结合受损提供了分子解释。
**3.3 Cel7A CD在糖溶液中的溶剂化性质**
CD中糖和水的MDDFs与CBM观察到的相似。在共溶剂中，葡萄糖表现出最高的与蛋白表面氢键结合倾向，其次是甘露糖和木糖。木糖MDDF峰变化适度。研究人员绘制了隧道残基（距底物路径5 ?以内的残基）对糖分布的特定贡献。糖优先在隧道入口处（涉及残基Q7、W40、S47、N49、D52、G53、T81、A100、Q101、K102和V104）及包围隧道的柔性环周围——特别是B2（N200）、B3（R251、Y252）和A3（Y370）——积累。相反，糖密度在核心催化三联体（E212、D214和E217）和内部芳香族残基（W38、W367和W376）周围保持显著较低。B3环——特别是R251和Y252残基周围——的独特溶剂化是一个值得注意的观察，该区域对产物协调和持续移动至关重要。糖在此区域的局部积累为上一节讨论的结构稳定化或RMSF降低提供了清晰的分子解释。通过与被认为固有柔性的残基相互作用，糖约束了酶沿纤维素链推进所需的动态波动。差异密度图进一步表明葡萄糖对催化酸/碱残基（E217）的偏好显著高于甘露糖或木糖，这支持了Zhai等提出的行为区分，为葡萄糖主要影响催化活性而甘露糖和木糖通过降低酶持续加工性抑制整体水解提供了结构解释。此外，甘露糖和木糖在W40附近的积累倾向表明隧道入口处有利的CH–π和疏水相互作用。由于W40在底物识别中发挥核心作用，这些单糖对其占据可能阻碍纤维素链进入催化隧道的初始结合和后续通过。木糖在T81周围的更高密度也值得注意，T81位于催化域入口处并通过氢键协调底物，木糖对其持续占据可能破坏这些初始极性相互作用，为持续加工性降低提供进一步潜在机制。
**3.4 溶剂化的域和糖依赖性**
CBM的水KBIs在不同糖类型和浓度下基本不变，而CD在相同条件下水KBIs显示轻微但渐进的增加。两个域的糖KBIs随糖浓度升高而降低，木糖始终对两个域保持最高亲和力。CBM在大多数浓度下表现出糖的优先溶剂化（Γ > 0），木糖在1.5 mol L?1处达到峰值0.96 ± 0.1。相反，CD始终产生负Γ值，在最高浓度下降至?27.6 ± 0.8。这种CD的一般排除与CBM的积累并存，表明两个域通过不同机制介导抑制。在CBM处，糖主要通过沿纤维素结合面的酪氨酸残基的芳香族–碳水化合物相互作用积累。同时，在CD处，酶效率可能由隧道入口附近及B2、B3和A3环周围的局部糖相互作用调节。模拟中观察到的CBM酪氨酸残基处更高的优先溶剂化与现有实验数据相关，后者显示木糖比甘露糖更有效地抑制Cel7A的生产性结合和持续加工性。通过更有效地占据这些关键芳香位点，木糖积累与CBM锚定结晶纤维素表面的可用性降低一致，为其相对于甘露糖对生产性结合更大影响提供了结构解释。
**讨论与结论**
开发更高效的酶混合物以减轻木质纤维素解聚过程中的产物抑制，需要对抑制剂相互作用的精确分子层面理解。模拟揭示了葡萄糖和半纤维素糖（甘露糖和木糖）与Cel7A两个主要功能域相互作用的显著差异。具体而言，CBM被这些糖优先溶剂化——主要由非特异性相互作用和纤维素结合面芳香族酪氨酸处的局部积累驱动——而CD在所研究浓度范围内保持优先水合。这种CBM结合界面的竞争性遮蔽与纤维素结合面锚定可用性降低一致，为先前报道的结合速率kon降低提供了结构解释。此外，尽管CD的广泛表面通常排除糖类，研究人员在隧道入口和隧道包围环（特别是B2、B3和A3）内识别了关键的局部相互作用。由此产生的结构稳定化和糖积累——特别是在催化出口附近的B3环内——为观察到的Cel7A整体持续加工性降低提供了分子机制。
