《Biochemistry》：Allostery between Distant Structural Regions Dictates Selectivity in GPCR:G Protein Coupling
编辑推荐：
本文深入探究了G蛋白偶联受体（GPCR）与Gα蛋白选择性偶联的分子机制，其核心是远端结构区之间的变构协同性。研究发现，Gα蛋白核心内的远端残基群落（如N末端、h4s6环和H5螺旋）通过复杂的变构依赖关系，共同调控着GPCR:Gα界面处的协同性相互作用，从而决定了Gα亚型（Gαs、Gαq、Gαi）的偶联选择性。通过整合可解释机器学习（贝叶斯网络模型）、分子动力学模拟和功能实验，作者不仅揭示了过去未知的协同热点，还成功通过Gαs核心的定点Gαq样突变，工程化改造了其偶联偏好，为理性设计和调控GPCR信号通路提供了新策略。
引言
G蛋白偶联受体（GPCR）是最大的膜蛋白超家族，也是主要的药物靶点。激动剂结合GPCR后，会激活受体并与特定的Gα蛋白亚型（如Gαs、Gαq、Gαi）偶联，从而触发下游信号通路。这种Gα亚型特异性相互作用塑造了下游信号输出，影响着多种生理和病理过程的功能特异性。尽管已有大量高分辨率结构研究，但GPCR:Gα蛋白偶联选择性的分子机制仍未完全阐明。这主要是因为界面是动态的，仅凭静态结构分析不足以解释选择性的决定因素，这一点在丙氨酸扫描突变研究中得到证实：某些在三维结构中不直接接触的界面残基，突变后仍能显著影响受体活性。理解GPCR如何在与Gα蛋白界面特征重叠的情况下实现选择性偶联，是设计能够精确调控信号通路疗法的关键一步。
以往研究通常将Gα蛋白的C端最后10个氨基酸残基（H5螺旋的“尖端”）和其余部分（“核心”）分开讨论。虽然H5尖端在Gα蛋白偶联选择性中的重要性已被广泛认知，但Gα核心对选择性的影响同样重要。然而，核心中赋予选择性的具体残基和/或结构区域尚不清楚，且对核心内不同结构区域如何变构地影响GPCR:Gα蛋白偶联界面的理解严重缺乏。
研究方法
本研究采用了一种可解释的机器学习方法——贝叶斯网络（BN）建模，对广泛的分子动力学（MD）模拟轨迹结合实验数据进行完全数据驱动的分析，以揭示Gα蛋白核心中赋予G蛋白偶联选择性的特定结构区域和残基。
模拟体系构建：研究基于多个GPCR:G蛋白三聚体复合物的实验三维结构（如6GDG、3SN6等）作为MD模拟的起点，涵盖Gαs、Gαq、Gαi偶联的受体。将复合物嵌入明确的POPC双层膜中，并使用CHARMM36m力场进行参数化。
分子动力学模拟：使用GROMACS2019对每个GPCR:Gα蛋白复合物进行了5微秒（5 x 1 μs）的全原子MD模拟。系统经过最小化、加热和平衡后，进行无约束的生产模拟，并监测骨架原子的均方根偏差以确保稳定性。
界面接触分析与贝叶斯网络建模：使用Python工具“GetContacts”分析GPCR与Gα亚基之间的分子间接触，包括盐桥、氢键、范德华力等。为每一帧生成二元相互作用指纹（存在接触为1，否则为0）。这些二元矩阵被输入到MD专用的BN建模软件BaNDyT中。BN模型将界面中的残基对表示为节点，将它们的直接概率依赖关系表示为有向边。通过最小不确定性评分函数和随机重启确保模型收敛。节点的加权度（所有边强度之和）被用来衡量其协同性，高加权度的接触可被视为协同接触。
功能实验验证：通过cAMP积累实验在HEK293ΔGsix细胞（缺乏六种功能性G蛋白的细胞）中验证预测的突变对G蛋白偶联功能的影响。利用SPASM传感器（受体与Gα蛋白通过柔性连接子串联）测试Gαs-to-Gαq突变体与Gq偶联受体（如血管加压素V1AR）的偶联能力。
结果
协同接触的鉴定与Gα亚型特异性
比较不同Gα亚型（s, q, i）的GPCR:Gα界面接触发现，超过50%的界面接触在不同亚型间是共享的。然而，高协同性的接触（加权度排名前25%）则具有高度的Gα亚型特异性：Gαs偶联复合物中88.6%的协同接触是Gαs特异性的，Gαq和Gαi中该比例分别为66.7%和64.7%。这表明Gα特异性相互作用驱动了界面协同性，可能对偶联选择性有重要贡献。
协同接触分布与选择性调控
通过分析协同接触在GPCR和Gα蛋白各结构区域的分布（以热图形式呈现），发现了Gα亚型特异的协同模式。Gαi和Gαq偶联复合物中，高协同性接触主要位于Gα的H5螺旋与GPCR的TM5、TM6之间。此外，Gαi复合物中，hns1和s2s3环与GPCR胞内环2（ICL2）之间存在协同接触。而在Gαs偶联复合物中，协同接触更均匀地分布在H5螺旋、hgh4区域、H4螺旋和h4s6环。这些差异化的协同性模式强化了假说：这些区域可能通过贡献于受体与Gα蛋白之间的选择性相互作用来允许偶联选择性。
基于协同性热点的功能转换
为了探究特定协同结构区域如何贡献于Gα蛋白偶联选择性，研究设计了一种将Gαs的偶联特征转换为Gαq的策略。选择突变位点的标准是：（1）在Gαs和Gαq中都具有高协同性评分（加权度）；（2）在这些蛋白质中氨基酸类型不同。据此预测了Gαs核心中多个结构片段的突变位点（，包括N末端、s2s3、H4、h4s6、S6和H5（C端）区域。
功能实验表明，将Gαs的H5螺旋突变为Gαq样序列，能显著降低由异丙肾上腺素（iso）诱导的、通过内源性β2AR的cAMP积累。更重要的是，将核心区域（N末端、h4s6等）突变为Gαq样序列，也能破坏Gαs信号，并且当与H5突变结合时（全核心突变体），效果更显著，类似于未转染细胞。这提示核心和H5螺旋区域的协同相互作用对于有效的Gαs信号传导都至关重要。
反过来，测试这些突变体与Gq偶联受体V1AR的偶联能力发现，Gαq样的H5突变体、全核心突变体以及核心+H5突变体，都能在血管加压素（AVP）刺激下引起V1AR的cAMP积累增加，获得了偶联至Gq通路的功能。特别是N末端+h4s6组合突变体引起的cAMP积累甚至比单独的Gαq样H5突变更显著，表明核心区域内的协同相互作用可能在调节偶联选择性中起关键作用。
H5螺旋受核心结构区域的变构调控
贝叶斯网络模型可量化不同结构区域之间的变构依赖关系。通过计算每个Gα核心结构区域与H5螺旋接触节点之间的边权重总和，可以衡量该核心区域对H5螺旋的协同影响强度。分析发现，h4s6环和N末端区域对H5螺旋具有很强的协同影响（）。这指导了选择N末端和h4s6区域进行旨在实现功能转换的突变分析。
对Gαs蛋白BN子网络的深入分析显示，核心/N末端与H5螺旋之间存在大量变构依赖关系（接触间距大于10?），表明一个区域的变化可以传播并影响另一个区域。关键接触对的相互作用能量时间序列分析进一步支持了这种变构关联。例如，N末端接触与H5螺旋接触的相互作用能量呈负相关，提示一种相互作用可以维持与受体的偶联。当在HN.52引入Gαs-to-Gαq转换突变后，使得Gαq样Gαs蛋白能够采用有利于成功Gαq偶联的构象。这些分析共同证明，核心区域、N末端和H5螺旋之间的协同依赖关系对于介导Gα蛋白与GPCR偶联的功能转换至关重要。
嵌合Gα蛋白的动态与协同性分析
Gα嵌合体（如Gαs核心-Gαq H5）已被广泛用于稳定和解析GPCR:Gα蛋白复合物的三维结构。通过MD模拟比较此类嵌合体与野生型Gα蛋白的动态差异发现，嵌合体的N末端在模拟中会从更接近Gαs WT的构象转变为更接近Gαq WT的构象。而嵌合体的H5螺旋在三维结构中与Gαq WT高度一致，但在模拟中会与Gαq WT构象发生轻微偏离，同时保持相对于Gαs WT的稳定构象，这很可能源于H5螺旋基部与保留Gαs样特征的核心之间的相互作用。
对5HT2A受体与Gαq嵌合体（Gαs核心-Gαq H5）复合物进行BN分析发现，其协同相互作用主要集中在H5介导的接触中，核心和N末端的贡献极少，这与野生型Gαq或Gαs偶联复合物形成对比。此外，嵌合体中核心与H5螺旋之间的变构依赖数量减少，且N末端+h4s6区域对H5螺旋的总协同性评分低于野生型Gαq蛋白。这表明嵌合体表现出Gα结构区域间变构影响的减弱，可能导致Gα与GPCR偶联的差异。这些发现与此前研究一致：Gαq蛋白的偶联更依赖于H5螺旋，而Gαs的偶联则很大程度上由Gα核心主导。
结论
本研究表明，GPCR:G蛋白偶联选择性由Gα亚基内的协同和变构相互作用所主导，其范围超出了经典的H5螺旋，延伸至更远端的核心区域。通过结合广泛的MD模拟与可解释的贝叶斯网络模型，研究鉴定出了G蛋白亚型特异性的协同残基群落，特别是涉及Gα蛋白N末端和h4s6环的区域，它们能变构调控H5螺旋并塑造选择性受体偶联。在这些预测的协同性热点指导下，对Gαs核心的定点突变成功改变了其偶联特征，将受体信号转向Gαq样信号。这些发现确立了Gα蛋白核心内的协同变构是GPCR信号特异性的一个重要机制，并为理性调控蛋白质-蛋白质相互作用提供了一个可推广的框架。
