《Frontiers in Plant Science》：Artificial intelligence-based structural modeling of PSTOL1 receptor kinases in rice and sorghum
编辑推荐：
最初在水稻（Oryza sativa OsPSTOL1）中被鉴定出的磷胁迫耐受1（PHOSPHORUS-STARVATION TOLERANCE 1, PSTOL1）蛋白是多个植物物种（包括高粱 Sorghum bicolor，SbPSTOL1）中磷利用效率的
最初在水稻（Oryza sativa OsPSTOL1）中被鉴定出的磷胁迫耐受1（PHOSPHORUS-STARVATION TOLERANCE 1, PSTOL1）蛋白是多个植物物种（包括高粱 Sorghum bicolor，SbPSTOL1）中磷利用效率的重要决定因素。与OsPSTOL1不同，SbPSTOL1同源蛋白具有跨膜受体样激酶（receptor-like kinase, RLK）结构域架构，并带有具多态性的胞外结构域。研究人员采用人工智能驱动的计算流程来研究SbPSTOL1蛋白的结构与构象动力学。使用多个深度学习平台进行的单体预测所得结构与SbPSTOL1的RLK拓扑结构不一致，提示其固有单体不稳定性或自抑制状态。在验证步骤中，研究人员还对系统发育上相关的ZmWAKL RLK进行了建模，该激酶与其可能的共受体ZmWIK形成异源二聚体，并在计算机模拟中重现了二者的相互作用。寡聚体建模生成了在结构上合理的SbPSTOL1蛋白模型，克服了单体结构不稳定的问题。膜环境中的分子动力学模拟也支持SbPSTOL1寡聚体状态具有优先稳定性。区域分析表明，跨膜结构域和激酶结构域较为刚性，而胞外区域表现出较高的柔性，这与它们在环境感知中的功能一致。值得注意的是，寡聚化与催化位点的构象转变相关，可能使典型的αC-Glu?β3-Lys盐桥得以形成，该盐桥是激酶活化的标志。综上所述，这些结果支持这样一种模型：寡聚化可能是SbPSTOL1（但不一定是OsPSTOL1）获得催化能力的一个结构前提。这些基于人工智能的发现为未来实验验证提供了可检验的假说，并首次呈现了SbPSTOL1蛋白的计算机模拟结构表征，为阐明复杂的RLK介导的植物信号传导机制提供了见解。
**基于人工智能的结构建模揭示SbPSTOL1寡聚化对其催化活性构象的必要性：跨膜受体激酶结构与功能研究**
**一、研究背景与科学问题**
全球粮食安全正面临气候变化与人口增长的双重压力，联合国粮农组织估计到2050年粮食产量需增加70%方能满足需求。酸性土壤（pH<5.5）占全球可耕地面积高达50%，其上的作物生产受到铝毒和低磷（phosphorus, P）有效性等多重非生物胁迫的严重限制。磷是能量代谢和光合作用所必需的大量营养元素，但在酸性土壤中易与铁、铝氧化物结合而难以被植物吸收利用。因此，提高作物的磷获取效率（Phosphorus Acquisition Efficiency, PAE）成为可持续农业和粮食安全的关键研究方向。
高粱（Sorghum bicolor L. Moench）作为耐受多种环境胁迫的C4作物，是解析作物非生物胁迫耐受性的优良模式系统。研究表明，PAE是低磷土壤条件下高粱籽粒产量的主要决定因素。在水稻中鉴定的PHOSPHORUS-STARVATION TOLERANCE 1（OsPSTOL1）基因是理解PAE的里程碑式发现，该基因被预测编码胞质激酶，可促进根系生长和磷获取。在高粱中，关联分析表明SbPSTOL1同源基因与低磷条件下的根系形态改善和籽粒产量提高相关，后续研究在玉米、小麦等物种中也验证了PSTOL1基因的功能保守性。然而，与水稻OsPSTOL1的胞质激酶性质不同，高粱SbPSTOL1蛋白具有跨膜受体样激酶（Receptor-Like Kinase, RLK）拓扑结构，且其胞外结构域呈现多态性。
植物RLK采用典型的三模块架构，包括负责信号感知的胞外配体结合结构域、单次跨膜α螺旋和保守的胞质激酶结构域。经典激活模型认为，配体结合诱导受体形成同源或异源寡聚复合物，进而通过相互转磷酸化解除自抑制状态，建立催化活性构象。激酶结构域从非活性到活性状态的转变依赖于一个真核激酶共有的分子开关：αC螺旋由外旋位置（αC-out）向内旋位置（αC-in）移动，使得β3折叠链上的保守赖氨酸（β3-Lys）与αC螺旋上的谷氨酸（αC-Glu）形成稳定的盐桥，从而正确定位ATP结合口袋。然而，PSTOL1激酶属于非RD类激酶，缺乏催化天冬氨酸之前的保守精氨酸，其激活可能不严格遵循经典RD激酶范式。目前，PSTOL1蛋白从抑制态向活性态转变的具体结构机制尚不清楚。
膜蛋白结构建模面临固有权衡：虽然人工智能（Artificial Intelligence, AI）方法在蛋白质结构预测领域取得了重大突破，但自动置信度分数（如pLDDT）的解读存在关键障碍——高置信度模型可能呈现与膜环境不相容的“拓扑坍缩”构象。因此，预测模型的验证必须超越纯统计指标，整合细胞拓扑学知识和功能基序保守性信息。该研究正是基于这一背景，旨在利用AI驱动框架系统表征SbPSTOL1蛋白的结构特性，揭示可能对催化能力至关重要的结构模式。
**二、主要技术方法**
研究人员开发了一套整合性的AI驱动计算流程，专门针对大型受体样激酶的结构表征。（1）序列与功能注释：使用InterProScan对三个SbPSTOL1同源蛋白（Sb03g006765、Sb03g031690和Sb07g002840，对应高粱基因组v5.1版本编号为Sobic.003G080700、Sobic.003G263400和Sobic.007G032900）进行保守结构域和调控基序注释，使用Clustal Omega进行多序列比对，并以ZmWAKL、ZmWIK（GenBank编号分别为OQ421108、OQ421112）和OsPSTOL1（BAK26566）作为参照。ZmWAKL-ZmWIK异源二聚体验证体系和它们的玉米来源是验证结构预测流程的关键参照——其物理互作和相互转磷酸化已被实验证实。（2）多模态结构建模：整合AlphaFold2（v2.1、2.2、2.3）、AlphaFold3、RoseTTAFold（v1、v2）和RoseTTAFold-All-Atom多个深度学习平台进行单体和寡聚体构象预测。（3）膜环境嵌入与分子动力学模拟：使用CHARMM-GUI Membrane Builder将模型嵌入POPC脂双层，通过PPM 2.0服务器预测膜坐标并与InterProScan域注释交叉验证，使用NAMD3引擎和CHARMM36m力场在NPT系综、303.15 K条件下进行模拟。（4）界面分析与区域化稳定性评估：使用COCaDA软件量化分子间界面接触，通过自定义TCL脚本在VMD环境中对胞外、跨膜和胞内三个区域分别进行基于最小二乘拟合的均方根偏差（RMSD）和均方根波动（RMSF）计算。模拟时间由RMSD结构收敛决定：Sb03g006765×Sb03g031690异源二聚体为300纳秒，Sb03g031690×Sb07g002840为125纳秒，Sb03g006765×Sb07g002840为约100纳秒，异源三聚体为约60纳秒。
**三、研究结果**
**3.1 SbPSTOL1结构域分析**：三个SbPSTOL1蛋白均展现RLK典型的三模块结构，具有N端信号肽、单次跨膜螺旋和保守的胞内激酶结构域（激酶结构域大小仅略有差异，介于303至307个氨基酸之间）。与缺乏膜锚定基序的胞质OsPSTOL1形成鲜明对比。结构差异主要体现在胞外结构域：Sb03g006765和Sb03g031690拥有WAK家族结构域（GUB_WAK和WAK_Assoc），而Sb07g002840具有LysM结构域。在Andropogoneae族的代表性物种中（甘蔗属、须芒草属和芒属），蛋白结构域架构大体保守，表明SbPSTOL1的胞外结构域并非高粱特有，而是在该族中多态性保守。
**3.2 激酶核心的序列保守性**：激酶结构域的多序列比对显示高粱和水稻PSTOL1蛋白间有极高保守性。全长氨基酸序列同一性为44%-72%，而分离的激酶结构域同一性达61%-74%，覆盖率为100%。ATP结合位点同一性58%-71%，活性位点同一性达92%-100%。所有PSTOL1蛋白均属于非RD激酶，具有HFD型催化环基序，在催化天冬氨酸前携带非精氨酸残基。β3链基序分别为IAVK（Sb07g002840和OsPSTOL1）或VAVK（Sb03g006765和Sb03g031690）变体，其中的赖氨酸（K）残基构成形成经典盐桥的关键残基。
**3.3 PSTOL1蛋白的结构建模**：所有单体的AlphaFold和RoseTTAFold（除RoseTTAFold1外）预测均出现了系统性的结构坍缩——胞外结构域与胞内激酶结构域发生直接接触和拓扑重叠，跨膜螺旋严重弯曲变形，尽管局部置信度分数（pLDDT介于75-80）较高。RoseTTAFold1是唯一能够建立符合预期三模块架构的软件，其模型呈现严格的空间结构域分离。PPM 2.0预测的膜插入坐标与InterProScan结构域注释精确对应，表明跨膜螺旋正确穿越脂双层，使WAK和LysM结构域朝向质外体。
**3.4 膜环境中的催化口袋可及性**：在RoseTTAFold1模型中，Sb03g006765和Sb03g031690的激酶结构域出现旋转倾斜，ATP结合位点和活性位点朝向膜内侧，存在膜空间位阻。相比之下，Sb07g002840和OsPSTOL1中三个催化位点汇聚形成溶剂可及的催化口袋，激活环的空间取向不阻碍底物进入。
**3.5 ZmWAKL-ZmWIK复合物建模验证**：只有AlphaFold3能够生成无结构坍缩的ZmWAKL-ZmWIK异源二聚体模型，展现出严格的胞外-胞内空间分离。COCaDA界面分析鉴定出131个链间接触，分布于胞外（40个）、跨膜（10个）和胞内（81个）区域。激酶结构域呈“背对背”反向平行排列，使两个单体的催化裂缝位于相对面，保持溶剂可及。该几何构型的精确重现为后续SbPSTOL1寡聚体建模提供了方法论支持。
**3.6 SbPSTOL1寡聚复合物的结构预测与界面分析**：AlphaFold3预测的所有SbPSTOL1寡聚体均维持正确的膜拓扑，无球蛋白坍缩。三种异源二聚体的激酶结构域均呈反向平行（背对背）排列，催化口袋位于复合体相对面且无界面阻碍。界面主要由激酶结构域间的胞内相互作用稳定（接触数78-97个）。在异源三聚体中，Sb03g006765充当中心支架，Sb03g031690和Sb07g002840的激酶结构域呈面对面取向，导致两者的催化口袋可能被空间阻碍，而Sb03g006765的催化口袋完全暴露于胞质。值得注意的是，所有复合物的激活环均参与界面相互作用，主要由保守赖氨酸残基介导（Sb03g006765的469K、Sb03g031690的512K和Sb07g002840的470K）。
**3.7 αC-Glu?β3-Lys盐桥分析**：在所有SbPSTOL1单体中，αC-Glu与β3-Lys的距离为5.37-5.95 ?，超过盐桥形成阈值（4 ?），表明处于自抑制的非活性构象。在寡聚体中，该距离降至3.33-4.00 ?，符合盐桥形成标准。异源三聚体达到最短距离（3.33-3.79 ?），所有三个亚基同时达到活性状态。OsPSTOL1作为单体时即呈现3.57 ?的距离，表明其盐桥可固有形成，无需寡聚化。
**3.8 分子动力学模拟**：所有寡聚体的胞内区域RMSD均稳定收敛至平台期（异源二聚体低于3.0 ?，异源三聚体为3.25 ?），无去折叠迹象。区域化分析显示，跨膜结构域和激酶结构域高度刚性（平均RMSD分别为0.69-1.24 ?和2.77-3.09 ?），而胞外区域表现出高柔性（平均RMSD 4.22-7.56 ?），表明跨膜结构域可充当机械绝缘体。在Sb03g006765×Sb03g031690异源二聚体、Sb03g031690×Sb07g002840异源二聚体及异源三聚体中，盐桥在97%-100%的模拟时间内保持活性构象。而Sb03g006765×Sb07g002840异源二聚体表现出明显不稳定性，盐桥仅在28.1%（Sb07g002840）和61.7%（Sb03g006765）的时间内维持活性，这与其缺乏Sb03g006765亚基激活环的互惠接触有关。
**四、讨论与结论**
该研究首次建立了高粱SbPSTOL1蛋白的整合结构框架。核心发现之一是AI预测工具常用的置信度分数与SbPSTOL1蛋白三模块拓扑结构之间的明确脱节——AlphaFold模型尽管具有高pLDDT分数，仍系统性地产生拓扑不一致的坍缩模型。这凸显了膜蛋白结构建模中仅依赖统计指标的根本局限，预测验证必须整合蛋白质拓扑学知识和功能基序保守性信息。
关于激活机制，该研究提出了寡聚化依赖的SbPSTOL1激活模型：与可直接形成盐桥的胞质OsPSTOL1不同，SbPSTOL1单体处于自抑制的非活性状态，需要通过寡聚化触发构象转变，使αC螺旋正确取向并建立催化盐桥。这一模型与植物RLK的经典别构激活范式（如BRI1-SERK复合物）一致，但由于SbPSTOL1属于非RD激酶，其具体调控机制可能不依赖于经典的激活环磷酸化范式。激活环中保守赖氨酸残基（469K、512K、470K）介导的界面接触可能是稳定αC螺旋正确取向的关键结构触发器。分子动力学模拟证实了异源三聚体虽整体结构稳定，但仅Sb03g006765亚基的催化口袋完全可及，提示其催化能力可能具有功能不对称性。该研究还提出，SbPSTOL1与OsPSTOL1的结构差异可能反映了不同的信号传导策略：SbPSTOL1通过配体或接触诱导的寡聚化将质外体或细胞壁信号与胞内信号传导偶联，而OsPSTOL1可能以更组成性、细胞自主的方式发挥作用。
这项基于人工智能的发现代表了可检验的假说，为未来实验验证提供了明确路线图：通过共免疫沉淀、双分子荧光互补（BiFC）或分裂荧光素酶互补实验验证预测的蛋白质互作；通过定点突变结合体外激酶活性测定验证寡聚化是否确实是建立盐桥和获得催化能力的先决条件。此工作不仅首次提供了SbPSTOL1蛋白的计算机模拟结构表征，还建立了一个可迁移的膜蛋白建模计算框架，为解析复杂的RLK介导的植物信号传导机制提供了新的见解。
