热适应GH1 β-葡萄糖苷酶的βα5–βα7动态耦合核心：进化约束与温度依赖网络重构在(β/α)8 TIM桶中的汇聚

《PLOS Computational Biology》：Evolutionary and dynamic constraints converge on a conserved coupling sector in thermally adapted GH1 β-glucosidases

　　

蛋白质序列如何在热胁迫下演化以维持功能性构象动力学，仍是分子进化的重要问题。嗜热酶常被认为主要通过增强结构刚性实现稳定性，但热适应也可能涉及蛋白骨架中相关运动的保持或重组。在此，研究人员比较了三种糖苷水解酶家族1 β-葡萄糖苷酶（glycoside hydrolase family 1 β-glucosidases, GH1 β-glucosidases）——BsBgl1A、9UPT和LQBG8——它们共享保守的(β/α)₈ TIM桶折叠，但最适温度分别约为45、65和70 °C。通过整合比较结构分析、在300、325和350 K下三个独立的1 μs分子动力学模拟（molecular dynamics simulations, MD）、残基水平动态耦合的神经关系推断（neural relational inference, NRI）以及家族级共进化分析（coevolutionary analysis），研究人员识别了与其热表型相关的不同动态策略。所有三种酶在模拟温度下均保持整体紧凑性和主要二级结构框架，表明均方根偏差（root-mean-square deviation, RMSD）和残基柔性的差异并不反映大规模去折叠。BsBgl1A表现出最强的温度敏感性，在350 K时结构偏差和残基迁移率增加。其推断的相互作用网络呈非单调响应，特征为325 K时显著过度耦合，随后在350 K时连接性降低。相比之下，9UPT在不同温度下保持相对受限的残基动力学和高度可重复的耦合拓扑。LQBG8结合了显著的局部构象灵活性与持续的长程耦合，并在350 K时保持最高的有向网络密度（directed network density）。在三种酶中，βα5–βα7区域反复作为共享的结构和动态中心出现。共进化耦合也富集于TIM桶相关区域，这些区域在空间上汇聚于该区段。上述发现表明，GH1 β-葡萄糖苷酶的热适应不能仅由统一的结构刚性化充分解释。相反，具有相同保守折叠的酶可采用不同的相关运动组织方式，而更强的热适应与高温下以βα5–βα7为中心动态耦合模式的持续或重组相关。动态耦合与共进化富集之间的对应关系与该结构区段上的长期进化约束一致，尽管这并不能确立网络作为独立功能单元受到直接选择。因此，研究人员的上述结果提供了一个可检验的框架，其中热稳定性反映局部稳定化、区域柔性和高阶动态组织韧性之间的平衡。

研究背景与问题：蛋白质在高温下维持功能完整性与构象动力学是分子进化的重要问题。嗜热酶面临热波动增强、非共价相互作用减弱和构象熵增加等热力学挑战；经典模型强调疏水堆积、盐桥和环区柔性降低等静态因素，但过度刚性化可能在稳定折叠的同时损害催化效率，因此功能适应需要稳定与柔性之间的平衡。已有共进化分析和MD模拟分别揭示长期进化约束与温度依赖构象行为，但二者长期脱节，动态相互作用网络究竟只是折叠结构的物理涌现属性，还是进化保守的组织单元，仍不清楚。研究人员比较三种共享高度保守GH1 (β/α)₈ TIM桶折叠、最适温度约为45、65和70 °C的β-葡萄糖苷酶BsBgl1A、9UPT和LQBG8；其中BsBgl1A和LQBG8来自热泉宏基因组资源，9UPT为实验解析的中间温度代表结构，以此检验热适应是否与特定全局耦合架构的持续或重组相关，而非统一刚性化。

关键技术方法：研究人员采用比较结构分析；对三种酶在300、325、350 K各开展三个独立1 μs全原子MD模拟，分析RMSD、回转半径（radius of gyration, Rg）、均方根涨落（root-mean-square fluctuation, RMSF）、二级结构指派（DSSP）、接触存活概率及RMSD–Q自由能景观；用NRI从C_α轨迹推断残基水平和βα模块水平有向动态耦合网络，并评估保留测试集可重复性；从InterPro条目IPR001360收集GH-1结构，经400–500 aa、RMSD<3 ?和TM-score>0.80（相对BsBgl1A或LQBG8）筛选，最终178个代表性GH1结构用于esm-msa-1b家族级共进化分析。

3.1 进化约束定义GH1 β-葡萄糖苷酶中的保守核心：结构比较和结构导向序列比对显示，三种酶均具典型(β/α)₈ TIM桶；BsBgl1A与LQBG8主链RMSD仅0.693 ?，尽管序列一致性为40.59%；八个βα单元分别叠加显示β链和相邻α螺旋大部分空间重叠，差异主要在溶剂暴露环区和末端。结论：TIM桶中央骨架受强进化约束，外周环区具有较高结构可塑性。

3.2 高温适应与TIM桶核心的选择性稳定相关：温度分辨MD显示三者Rg均保持窄幅变化，无持续大规模膨胀；BsBgl1A在350 K时RMSD和RMSF显著增加、局部二级结构转换增多；9UPT残基波动最低、二级结构最连续；LQBG8本征柔性较高但保持紧凑折叠和主要二级结构；RMSD–Q景观显示BsBgl1A在350 K向高RMSD、低天然接触分数Q的盆地移动，9UPT最受限，LQBG8采样最广并保留较多天然接触；催化谷氨酸C_α–C_α距离总体稳定。结论：热适应并非统一刚性，而是局部稳定、区域柔性和天然接触保持的不同平衡。

3.3 热稳定性源于分布式动态相互作用网络的保持：NRI推断残基和模块有向耦合，阈值敏感性分析支持温度依赖趋势；BsBgl1A在325 K平均总度、加权节点强度和有向网络密度升高，350 K降低，呈非单调过度耦合后连接减少；9UPT网络组织较温度不敏感且较可重复；LQBG8保持广泛互联，350 K仍有较高有向网络密度和持续βα5–βα7耦合。三者中βα5–βα7反复为高权重耦合区；保留测试集模块级Spearman相关系数为0.548–0.952，模块级组织比单条边更可重复。结论：更强热适应与βα5–βα7中心耦合模式的持续或重组相关，但NRI概率是统计耦合而非物理因果通路。

3.4 进化约束锚定热应激下的动态网络核心：GH1数据集序列一致性矩阵显示广泛序列分化，而结构RMSD矩阵高度保守；共进化耦合非均匀，富集于β3–β4、β4–β5、α5、β6–α6、β6–β7、β7–β8等桶相关区，并在空间上汇聚于βα5–βα7；BsBgl1A和LQBG8均呈此分布，外周βα环和末端信号较弱。结论：动态通信核心与家族级共进化约束空间对应，提示该区段受长期进化约束，但相关不等于因果，不能证明网络作为独立选择单元。

讨论与结论：讨论指出，热稳定性应理解为稳定化相互作用与功能所需柔性之间的平衡，局部堆积、静电作用、区域刚性、溶剂组织、质子化状态和TIM桶内在力学均可能共同贡献三种酶的热表型。三种酶均保持TIM桶折叠；BsBgl1A温度敏感且网络响应非单调；9UPT构象受限且高阶耦合较可重复；LQBG8局部柔性与长程耦合并存。βα5–βα7的突出既可能源于TIM桶拓扑几何约束，也可能受序列适应调节。NRI边为统计耦合，需突变、扰动模拟、非平衡分析和实验验证；共进化与NRI重叠为相关性，需置换检验或结构感知零模型；模拟虽有三重复，网络收敛性仍需时间分辨NRI、分块比较、留一重复和重复训练评估；催化位点C_α距离不证明侧链几何或活性保持；同源模型在外周环区仍有不确定性；仅三种酶的结果应作为可检验假设。研究结论翻译如下：总之，BsBgl1A、9UPT和LQBG8在保守TIM桶骨架内以不同方式组织构象运动。BsBgl1A表现出显著温度敏感性，在350 K向高RMSD、低Q集合移动，并呈325 K中间温度过度耦合后连接减少的非单调网络响应。9UPT维持相对受限的构象景观和持续的高阶动态组织，而LQBG8结合显著构象可塑性与持续长程统计耦合及较高比例天然接触保持。在这些不同动态 regime 中，βα5–βα7区域一致作为突出的结构、动态和进化区段出现。上述发现支持如下假设：GH1 β-葡萄糖苷酶的热适应可能不仅涉及局部稳定化和刚性调节，还涉及保守桶核心内高阶相关动态模式的持续或重组。