比较增强采样方法在探索β-Catenin17–48构象空间中的应用

《The Journal of Physical Chemistry B》：Comparing Enhanced Sampling Methods in Exploring the Conformational Space of β-Catenin17–48

　　

内在无序蛋白（IDPs）在细胞信号传导和调控中发挥关键作用，但其动态构象景观使其难以通过实验和计算方法进行表征。磷酸化是最常见的翻译后修饰之一，常发生在内在无序区域，并可调节蛋白质结构和功能。增强采样分子动力学方法为更高效地探索IDPs可及的多种构象提供了潜在途径，但对其在采样此类系统构象景观方面的性能进行系统比较的研究仍然有限。在本研究中，研究人员评估了内在无序β-catenin^17–48肽段在其非磷酸化和磷酸化状态下的构象采样，使用了三种增强采样方法：高斯加速分子动力学（GaMD）、元动力学（METAD）和加权系综模拟（WESTPA），并与常规分子动力学模拟进行比较。研究人员探索了两种集体变量（CVs）以引导采样：磷酸化位点Ser33和Ser37的?二面角以及肽段的末端距。研究人员发现，不同的增强采样方法探索了构象空间的不同区域，而非收敛于单一系综，其中GaMD与无偏模拟有较大重叠。值得注意的是，METAD和WESTPA更容易进入无偏模拟中未观察到的构象区域。对组合构象系综的分析识别出连接非磷酸化和磷酸化状态的中间构象，这些构象在使用自适应策略的模拟中被优先采样。此外，与末端距指标相比，Ser33/Ser37 ?角CV更能有效捕获磷酸化依赖的构象变化。总之，这些结果强调了增强采样策略的选择和集体变量的选择如何影响IDP构象景观的探索。

**论文解读：增强采样方法在β-Catenin^17–48构象空间探索中的比较研究**

**研究背景与问题**
内在无序蛋白（IDPs）在细胞信号传导和调控中扮演核心角色，但其缺乏稳定三级结构，以动态构象系综存在，这使其实验和计算表征极具挑战。磷酸化作为最常见的翻译后修饰之一，常发生在内在无序区域，可显著调节蛋白质构象与功能。然而，IDPs构象空间广阔且复杂，常规分子动力学（MD）模拟难以充分采样。增强采样方法如高斯加速分子动力学（GaMD）、元动力学（METAD）和加权系综模拟（WESTPA）有望提高采样效率，但针对IDPs的系统性性能比较仍十分有限。为此，研究人员以β-catenin N端内在无序区域肽段（残基17–48，β-catenin^17–48）为模型，评估不同增强采样策略在捕获磷酸化诱导构象变化方面的能力。该肽段具有明确实验基准：Ser33和Ser37磷酸化可诱导从无序系综到含N端α-螺旋片段的转变，且该区域磷酸化状态与β-catenin定位、活性及多种癌症相关。研究旨在回答：不同增强采样方法如何影响构象景观探索？能否识别磷酸化相关的中间态？集体变量（CV）选择如何影响采样效果？该研究发表于《The Journal of Physical Chemistry B》。

**主要关键技术方法**
研究人员采用三种增强采样方法：GaMD（通过非自适应提升势平滑能量面）、METAD（沿预定义CV自适应沉积偏置势）和WESTPA（通过自适应采样重分配计算资源），并与常规MD模拟对比。所有模拟均使用Drude-2019极化力场（Drude FF），以NMR实验结构为起点。定义了两类CV：局部CV（Ser33和Ser37的?二面角）和全局CV（末端距）。通过能量景观可视化方法（ELViM）进行降维投影，并采用k-means N-Ary Natural Initiation（NANI）聚类算法定量分析构象空间探索。模拟准确性通过计算系综平均NMR标量耦合（³J_HN-Hα）和化学位移（SPARTA+软件）并与实验数据比较来验证。

**研究结果**

**1. 构象系综的比较**
通过ELViM投影和聚类分析发现，不同增强采样方法探索了互补的构象区域，而非收敛于单一共享系综。GaMD采样与无偏MD高度重叠，尤其在磷酸化状态下，表明其主要富集已访问构象。METAD和WESTPA则更频繁地进入无偏模拟未观察到的区域，且METAD的探索能力受磷酸化状态影响：非磷酸化时与无偏模拟重叠少，磷酸化时重叠增加（尤其当使用?角CV时）。WESTPA在非磷酸化状态下桥接了无偏MD、GaMD和METAD的构象系综，但在磷酸化状态下与其他方法重叠最小。聚类分析进一步显示，METAD和WESTPA采样了无偏MD贡献极少的簇，体现了自适应策略对未探索区域的促进能力。

**2. 磷酸化对β-Catenin^17–48系综的影响**
ELViM投影表明，磷酸化使构象空间更紧凑且局部化，与非磷酸化状态明显分离。当使用末端距CV时，两状态完全分离；而使用?角CV时出现部分重叠，这些重叠区域主要由METAD和WESTPA采样贡献。聚类分析（11个簇）中，?角CV下出现了4个包含两状态至少5%贡献的混合簇（簇2、5、6、10），可能代表磷酸化诱导构象转变的中间态。进一步分析簇中心结构和螺旋度发现，簇5和6（主要来自磷酸化模拟）与含N端α-螺旋的磷酸化NMR结构相似；簇7和11（主要来自非磷酸化模拟）与紧凑无序的非磷酸化结构一致。簇10（主要由非磷酸化METAD采样）表现出两状态的螺旋度特征，提示非磷酸化系综可瞬时访问类似磷酸化的构象。这些中间态主要由METAD和WESTPA捕获，表明自适应方法能识别稀有但功能相关的构象。CV比较显示，?角CV比末端距更能灵敏反映磷酸化诱导的局部骨架变化，从而促进跨状态采样。

**3. 与NMR观测值的验证**
所有方法产生的构象系综与实验NMR数据（³J_HN-Hα耦合和Cα、Cβ化学位移）整体吻合良好。非磷酸化模拟与实验一致性较高；磷酸化模拟中N端残基的计算值略有偏差（耦合偏低、Cα化学位移偏高），但Cβ化学位移吻合较好。不同方法间RMSE值相近，其中WESTPA模拟略优（RMSE更低）。这表明尽管采样行为存在差异，各方法均能准确捕获IDP的系综性质。

**总结与结论**
本研究系统比较了GaMD、METAD和WESTPA在探索β-catenin^17–48构象景观中的表现。主要结论如下：增强采样方法探索了互补的构象空间，而非收敛于单一系综；GaMD主要富集无偏模拟已访问构象，而METAD和WESTPA更易进入未采样区域；自适应方法（METAD和WESTPA）能识别连接非磷酸化和磷酸化状态的中间构象；CV选择至关重要，局部?角CV比全局末端距更能捕获磷酸化依赖的构象变化；所有方法产生的系综均与实验NMR数据一致，验证了其准确性。这些结果不仅展示了增强采样方法在IDP构象研究中的能力，也为理解磷酸化调控IDP结构及功能的分子机制提供了见解。未来工作可探索更高维CV、并行偏置策略，并推广至更广泛的IDP体系。