《Proteins: Structure, Function, and Bioinformatics》：Distinct Conformational Switches Control WPD-Loop Dynamics Across Protein Tyrosine Phosphatase Subfamilies
编辑推荐：
蛋白质酪氨酸磷酸酶（PTPs）通过一个保守的WPD-环在开放（非活性）与闭合（活性）构象之间切换，平衡磷酸酪氨酸水平，从而调控细胞信号传导。尽管这一构象变化对催化至关重要，但控制该运动的结构因素仍知之甚少。研究人员使用改良的局部空间模式（LSP）算法和XGBo
蛋白质酪氨酸磷酸酶（PTPs）通过一个保守的WPD-环在开放（非活性）与闭合（活性）构象之间切换，平衡磷酸酪氨酸水平，从而调控细胞信号传导。尽管这一构象变化对催化至关重要，但控制该运动的结构因素仍知之甚少。研究人员使用改良的局部空间模式（LSP）算法和XGBoost机器学习分析了551个PTP晶体结构，以识别其侧链构象影响WPD-环动力学的核心残基。分析揭示了两个不同的残基簇：一个以蛋白质核心中的F225（PTP1B编号）为中心，另一个位于距WPD-环约20 ?的α2螺旋上。值得注意的是，研究人员发现了涉及第225位保守苯丙氨酸的亚家族特异性调控机制。经典PTPs（如PTP1B）将该残基主要维持在“向下”构象，与WPD-环状态无关，而SHP2表现出独特的构象灵活性，其等效残基F469在“向上”和“向下”构象之间呈现双重占据。这表明不同的PTP亚家族在维持整体催化框架的同时，演化出了不同的变构控制机制。这些亚家族特异性的构象开关具有重要的治疗意义。PTP1B作为胰岛素信号传导的负调控因子，使其成为有吸引力的糖尿病治疗靶点，而SHP2在生长因子信号传导和肿瘤发生中的功能则提示了癌症治疗的机会。研究结果为理解亚家族特异性的WPD-环调控提供了结构框架，并确定了选择性治疗干预的位点，为开发规避传统活性位点导向方法局限性的高选择性抑制剂提供了策略。
**论文解读：蛋白质酪氨酸磷酸酶亚家族中WPD-环动力学的差异性构象开关调控机制**
**1. 研究背景与科学问题**
蛋白质酪氨酸磷酸酶（Protein Tyrosine Phosphatases, PTPs）通过催化底物蛋白上磷酸酪氨酸残基的去磷酸化，与蛋白酪氨酸激酶协同调控细胞信号转导网络，在细胞增殖、分化、代谢及免疫应答等众多生理过程中发挥关键作用。PTPs功能失调与多种人类疾病密切相关，如PTP1B作为胰岛素信号通路的负调控因子，其过度表达或活性增强与2型糖尿病和肥胖的发生发展紧密相关；而SHP2作为生长因子信号通路的正调控因子，其功能获得性突变则频繁见于多种白血病和实体瘤中。因此，PTPs已成为糖尿病和癌症等重大疾病药物研发的重要靶点。
PTP家族成员含有一个保守的催化结构域，其中一段柔性的WPD-环（包含催化天冬氨酸的环状结构）在催化过程中发生从开放（非活性）到闭合（活性）的构象转变，该转变过程是PTP催化反应的限速步骤。尽管WPD-环的构象变化对催化活性至关重要，但控制这一运动的精细结构因素仍不清楚。此外，由于PTPs活性位点在家族内高度保守，开发针对特定亚家族成员的高选择性抑制剂面临巨大挑战，传统靶向活性位点的抑制剂策略往往因选择性差和生物利用度低而受限。因此，深入理解PTPs中WPD-环运动的调控机制，对于阐明其功能多样性和开发新型选择性治疗策略具有重要意义。
**2. 研究方法概述**
研究人员从RCSB蛋白质数据银行（RCSB Protein Data Bank）中检索并筛选了551个I类经典PTP晶体结构（截至2023年11月14日），涵盖了不同家族成员及不同配体状态。研究采用了改良的局部空间模式（Local Spatial Patterns, LSP）算法来分析残基对在不同结构间的空间排列一致性。该算法以侧链重原子质心（Center of Mass, CoM）替代传统的Cβ原子，计算残基对间的Cα-Cα距离、Cα-CoM交叉距离及二面角，从而适应低序列同一性结构集的分析。随后，研究人员应用XGBoost梯度提升机器学习算法，以WPD-环的开放/闭合状态为分类目标，对LSP矩阵进行二分类建模，并通过五折交叉验证优化超参数。为验证家族水平分析的有效性，研究人员还进行了消融分析（Ablation Analysis），即仅使用PTP1B结构训练模型。此外，研究人员对PTP1B（PDB ID: 1ptv）和SHP2（PDB ID: 2SHP）分别开展了三次重复的分子动力学（Molecular Dynamics, MD）模拟（PTP1B为250 ns，SHP2为260 ns，前100 ns作为平衡舍弃），以直接观测关键残基F225和F469的构象动态行为。
**3. 研究结果**
**3.1 机器学习模型识别出的关键特征**
通过XGBoost模型的训练，研究人员获得了区分WPD-环开放与闭合状态的特征重要性排序。模型在测试集上取得了0.908的受试者工作特征曲线下面积（ROC-AUC），表明其分类性能优异。特征重要性分析显示，多个位于PTP结构域疏水核心的残基对（如F225-C226和F225-I261）对正确分类贡献显著。此外，特征重要性排名前10的残基对还包括F191-F194等。通过将每个残基参与的所有特征重要性求和，研究人员识别出了20个对预测WPD-环构象状态最重要的残基，其中F225得分最高（归一化得分为1.000）。
**3.2 高重要残基的结构定位**
将高重要性残基映射到PTP结构域的三维结构上，研究人员观察到两个明显的残基簇。第一个簇以F225为中心，包含C226、I261、R254、D265、Y271、V108、M109、F191、R199、L234、M235、S270和F194等残基，这些残基主要位于α3、α4、α6螺旋、β3折叠片和Q-环上。第二个簇包含位于α2螺旋上的残基，包括H94、W96、P89和N90，该螺旋距WPD-环超过20 ?。此外，还有少量残基（如I72和R254）显示高重要性但对间距离较远。总体上，两个高得分残基簇之间的连接较少，仅有通过β1、β2、β3折叠片的少数链接。使用Girvan-Newman社群算法对残基相互作用网络图进行分析，进一步证实了这些簇的相对独立性。
**3.3 F225的亚家族特异性构象行为**
针对重要性最高的F225残基，研究人员分析了其在所有551个PTP结构中的χ1二面角分布。结果显示，在PTP1B及大多数经典PTPs中，不管WPD-环处于开放还是闭合状态，F225都主要采用约180°的“向下”构象，朝向D229而非催化精氨酸R221。然而，一个显著的亚家族特异性差异被发现：F225（或等效位置）的“向上”构象（约290°）仅在SHP2的结构中被观察到。在SHP2的apo结构（PDB ID: 3ZM2）中，等效残基F469呈现出“向上”和“向下”两种构象各占0.5占据率的双重构象。这说明SHP2演化出了独特的调控机制，其F469能够在两种构象间切换，而这一动态灵活性在PTP1B等其他家族成员中是缺失的。唯一的例外是PTPN22的一个结合磷酸化肽段的结构（PDB ID: 3BRH），其同源位置为异亮氨酸而非苯丙氨酸。
**3.4 消融分析验证全家族分类方法的必要性**
为验证全家族XGBoost模型所识别的特征反映了真实的跨亚家族信号，研究人员仅用357个PTP1B结构重新训练模型。该PTP1B-only模型取得了0.877的平均交叉验证ROC-AUC。重要的是，残基重要性排序发生了根本性变化：在全家族模型中排名第一的F225在PTP1B-only模型中完全消失于前20名之外，取而代之的是F256、K255、L144、I145、V49、S50等反映PTP1B特有结构特征的残基。W96在两个模型中均出现。F225的消失可归因于：在PTP1B内部，F225在WPD-环开放和闭合状态下均采用“向下”构象，因此无法为分类提供判别信息。而在全家族模型中，F225的重要性源于其在不同亚家族间的构象变异性，特别是SHP2中独特的“向上”旋转异构体。
**3.5 分子动力学模拟确证F225/F469的动态构象开关**
MD模拟结果（图7）显示，PTP1B和SHP2两个体系在所有重复中均达到稳定平衡（RMSD集中在1.0-1.5 ?）。主链φ/ψ二面角分析表明，F225和F469位置的主链构象在两种蛋白质中是等效的，均形成一个紧密的单一簇。然而，侧链χ1/χ2二面角分布呈现显著差异：PTP1B的F225在三次重复模拟中始终占据单一的“向下”旋转异构体簇（χ1≈-100°，χ2≈80°），几乎未发生向其他构象的偏移，呈现构象锁定状态；而SHP2的F469在三次重复中均广泛采样不同的旋转异构体状态，同时涵盖与PTP1B F225等效的“向下”构象和另一种“向上”构象。该动态切换行为直接验证了SHP2 apo晶体结构（PDB: 3ZM2）中F469的双重占据现象，证明其并非晶体学伪影，而是SHP2内在的动态平衡特性。
**4. 讨论与结论**
本研究的核心发现是识别出PTP家族中一种此前未被认识的、由F225介导的亚家族特异性构象开关。结合已有的F225Y激活突变研究、共进化偶联分析以及MD模拟结果，充分确立了F225在PTP催化调控中的关键角色。在PTP1B中，F225的构象锁定可能通过维持刚性的疏水核心，牺牲动态切换能力以换取更稳定的能量景观，这与其作为组成型表达的胰岛素信号负调控因子的功能相适应。而SHP2中F469的动态切换则可能为其在生长因子信号级联中发挥正调控因子所需的精细调节提供了结构基础。值得注意的是，α2螺旋区域的残基簇（P89、N90、H94、W96）与研究团队此前对PTP1B闭合态模拟和TbPTP1/YopH模拟中识别的高中心性节点相吻合，提示该区域可能参与WPD-环闭合态的稳定。此外，图S1中识别出的已知变构残基（如F191、R199、K197）与本研究识别的簇分离，暗示存在一种独立的变构调控机制，且该机制在不同PTPs中并非完全保守。消融分析进一步证明了全家族分析对于揭示亚家族水平变构机制的不可或缺性。
总之，通过对551个PTP晶体结构的系统性分析，本研究揭示了PTP家族中WPD-环运动的两个关键调控簇：一个以F225为核心的疏水核心簇和一个位于α2螺旋的远端调节簇。F225/F469被证实是一种亚家族特异性的构象开关，在PTP1B中锁定于“向下”构象而在SHP2中呈现动态可切换性。这一发现为PTP亚家族特异性调控机制提供了结构基础，并为开发选择性变构抑制剂以治疗糖尿病（靶向PTP1B）和癌症（靶向SHP2）提供了新的药物设计位点和理论框架。该LSP-机器学习联合方法也为研究其他酶家族中类似的功能调控机制提供了有力的分析工具。
