Gα亚基与Gβγ结合界面中保守残基呈现Gα亚型特异性的Gβγ偶联作用

《Journal of Biological Chemistry》：Conserved Residues in the Gα interface show subtype specificity in Gβγ coupling

　　

尽管G蛋白(G protein)与G蛋白偶联受体(G-protein coupled receptor, GPCR)选择性相互作用的结构基础已被较好阐明，但不同Gα亚型(Gαs、Gαi、Gαq、Gα12等)与Gβγ选择性相互作用的分子机制仍不清楚。通常认为Gα亚型中的保守残基功能相似，但它们可能通过改变Gα:Gβγ界面接触的频率与稳定性来调控偶联选择性。研究人员采用全原子分子动力学(Molecular Dynamics, MD)模拟结合可解释的机器学习方法——贝叶斯网络模型(Bayesian Network Model, BNM)及蛋白-蛋白邻近(protein-protein proximity)生物发光共振能量转移(Bioluminescence Resonance Energy Transfer, BRET)实验，发现亲缘关系较近的Gαi/o与Gαq/11亚家族中的保守残基对Gβγ偶联产生差异性贡献。这些位于Gαi1和Gαq上的保守残基"热点(hotspot)"对Gβγ偶联产生不同的功能效应，表明序列保守并不等同于功能等效。研究结果表明，局部微环境和旁系同源亚型特异性别构偶联(allosteric coupling)共同塑造了保守界面残基对蛋白-蛋白偶联的贡献。该分析框架为解析亚型特异性机制提供了系统性方法，对药物设计及疾病相关变异的功能注释具有参考意义。本文所用计算方法可广泛适用于其他同源蛋白家族的分析。

论文解读：Gα界面保守残基对Gβγ偶联呈现Gα亚型特异性

研究背景与意义

异源三聚体G蛋白由Gα、Gβ、Gγ亚基组成，Gα激活后与Gβγ解离并启动下游信号。人类含有16种Gα亚型，分为Gα_s、Gα_i/o、Gα_q/11、Gα_12/13四个亚家族，各亚家族触发的效应通路不同。尽管GPCR与Gα选择性偶联的结构基础已有研究，但不同Gα亚型与Gβγ选择性偶联的机制尚不明晰，且普遍认为Gα亚型间保守残基功能相同。然而已有研究表明保守残基未必功能等效，其局部微环境和远程别构通讯可调节残基对信号输出的贡献。Gα突变与多种疾病相关（如McCune-Albright综合征、癫痫性脑病、葡萄膜黑色素瘤），明确保守残基在旁系同源Gα亚型中对Gβγ偶联的差异性影响，对理解G蛋白信号精细调控及疾病突变解读具有重要意义。本研究以序列和结构高度相似的Gα_i1（代表Gα_i/o亚家族）和Gα_q（代表Gα_q/11亚家族）为对象，探究Gα:Gβγ界面保守残基是否对Gβγ偶联具亚型特异性作用。研究结果发表于《Journal of Biological Chemistry》。

主要关键技术方法

研究人员以PDB 6CRK（Gα_i1:Gβ₁γ₂:GDP）和PDB 7F6G（apo Gα_q:Gβ₁γ₂）为起始结构，将GDP从6CRK移植至Gα_q构建Gα_q:Gβγ:GDP复合物模型，去除膜锚定所需的αN螺旋，进行体系构建与能量优化。采用CHARMM36m力场、GROMACS 2022.3执行全原子MD模拟，各复合物做5次独立1 μs模拟（合计各5 μs），取后600 ns轨迹分析。基于GetContacts计算Gα:Gβγ界面残基接触频率(contact frequency)与成对相互作用能（范德华+库仑），得每个Gα残基的平均接触频率和相互作用能。采用研究人员自研BaNDyT软件包构建贝叶斯网络模型(Bayesian Network Model, BNM)，以残基与周围12 ?内残基的相互作用能为输入，离散化后推断残基动态运动的协同依赖性(co-dependency)，计算节点加权度(weighted degree, WD)表征残基协同性。整合接触频率、相互作用能、BNM加权度三特征做分位数排序与距原点距离分析，预测亚型特异性热点。选取距原点较远的16个保守界面残基做丙氨酸扫描突变，利用TRUPATH BRET（基于神经降压素受体1型 NTSR1 激活诱导Gα–Gβγ解离致BRET比值下降）检测突变对Gα_i1和Gα_q与Gβγ解离的影响，验证计算预测。BNM进一步用于构建Gα内部区域间别构通讯的元网络(meta-network)及加权汉明距离(weighted Hamming distance)量化Gα_i1与Gα_q别构网络重排。

研究结果

接触频率和相互作用能显示Gα_i1与Gα_qα2螺旋处具亚型特异性热点

研究人员通过MD模拟计算Gα_i1与Gα_q各自Gα:Gβγ界面残基的接触频率与相互作用能。两亚型各有约59%–66%界面残基严格保守。对保守且持续接触（接触频率>20%）的16个残基分析发现，4个极性残基在两亚型间接触频率具显著差异（双尾Welch's t检验 p<0.05），差异残基集中于Ras样结构域的α2螺旋与β4-α3环——即Gα:Gβγ界面的外周区域。相互作用能分析揭示更高比例极性(40%)和非极性(33.3%)保守残基具显著亚型差异。接触频率和相互作用能均显示差异的共有残基为R205/R210（Gα_i1/Gα_q编号，下同）和H213/H218。结果表明仅靠热力学与静态接触不足以完全区分亚型特异性，需结合动态协同性分析。

贝叶斯网络分析揭示Gα_i1与Gα_qGα:Gβγ界面局部与别构协同依赖性

将MD轨迹输入BaNDyT构建BNM，按加权度(WD)取前四分位"top residues"。Gα_i1独有的top residues含R144、V233、L234；Gα_q独有的含Y151、Q152、G207、E221，均为物种间高度保守残基。高WD残基空间聚集于Ras样结构域与α螺旋结构域(alpha-helical domain, AHD)交界的αD-αE环与β4-α3环——该区在Gα结构域开合运动中参与核苷酸交换。BNM识别出7个Gα:Gβγ界面保守残基（R144/R149、Y146/Y151、Q147/Q152、G202/G207、E216/E221、V233/V238、L234/L239）在Gα_i1与Gα_q间协同性不同，但在接触频率和相互作用能上无显著差异，说明协同运动信息是前两特征之外的重要补充维度。

基于距离法预测残基为亚型特异性热点

为综合三特征，研究人员将各界面残基在接触频率、相互作用能、WD三个维度的分位数排名在两亚型间求差值Δ，计算三维空间中距原点(Δ=0,0,0)的欧氏距离，假设距原点越远越可能是亚型特异性热点。选取距原点较远的16个保守界面残基做丙氨酸突变，通过TRUPATH BRET检测NTSR1激活下Gα–Gβγ解离。BRET剂量-响应曲线显示三种模式：(i)仅影响Gα_q（如R144A/R149A、E145A/E150A、Y146A/Y151A、V233A/V238A）；(ii)影响两亚型但程度不同（如G202A/G207A、H213A/H218A、L234A/L239A、E236A/E241A、N241A/N246A、F259A/F264A）；(iii)同等影响两亚型（如E245A/E250A）。16个测试残基中9个为Gα_i1热点，14个为Gα_q热点，即Gα_q对界面突变更敏感。10个被判定为亚型特异性热点，空间上分布于Ras样结构域与Gβγ接触区（G202/G207、H213/H218）、靠近GPCR界面的F259/F264，以及Ras-AHD"夹钳(clamp)"区（R144/R149、E145/E150、Y146/Y151、V233/V238）——该夹钳区在核苷酸交换时介导Ras样结构域与AHD的开合。累积分布函数(cumulative distribution function, CDF)分析证实亚型特异性热点距原点显著大于共同热点（CDF右移），其中接触频率单独区分能力最强。证明三特征整合可有效预测Gα:Gβγ界面保守残基的亚型特异性偶联作用。

两Gα亚型功能重要区域间的别构通讯

利用BNM边强度求和量化结构区域间协同运动。Gα_i1与Gα_q均显示Gα:Gβγ界面所在区（αD-αE环、α2螺旋、β4-α3环）与远处残基（如α1-β1环、αA螺旋、β5-αG环）具别构相关性，且GDP结合部位（β1-α1环）与Gα:Gβγ界面（β4-α3环）间存别构关联。值得注意的是Gα_q的α5螺旋除与Ras样结构域区域别构关联外还与AHD区域（α2、αC、αF）别构关联，而Gα_i1的α5螺旋仅关联Ras样结构域。构建以Gα与各互作蛋白（GPCR、Gβγ、RGS蛋白、RIC-8、GDP）界面为元节点的元网络，两亚型全局拓扑相似——GPCR:Gα界面与Gα:Gβγ界面协同性最强，其次为GPCR:Gα与RIC-8及RGS界面协同，反映GEF(GPCR/RIC-8)与GAP(RGS)调控环路的内在偶联。但加权汉明距离分析显示个体残基对间别构依赖显著重排(rewiring)：GDP结合位点?Gα:Gβγ界面边权重在Gα_q总体更强且分布更广，部分边在Gα_i1更强，说明虽整体别构通讯总量相当，但具体通路经亚型特异性重组，Gα_q更分散的连通性可能解释其对突变更高的敏感性。

讨论与结论

研究人员总结：Gα_i1与Gα_qGβγ结合界面上的保守残基并非功能等效，其通过残基接触的时间持续性、焓及协同运动三者共同作用产生亚型选择性偶联，局部微环境与旁系同源特异性别构网络塑造此差异。除经典的Switch II(α2螺旋)与Switch III区域外，AHD（特别是Ras-AHD夹钳区）保守残基突变可选择性影响Gα_q而非Gα_i1的Gβγ解联，提示AHD参与Gβγ解离调控。所整合的计算框架（接触频率＋相互作用能＋BNM加权度）成功预测并实验验证了亚型特异性热点，其中多个热点位点为已知疾病相关突变位点，说明这些位点在亚型中可能有不同致病机制。研究强调疾病相关Gα突变应在受体–效应器互作网络中评估，旁系同源信号蛋白即便序列结构高度保守也可具亚型特异性突变敏感性。该策略可推广至其他同源蛋白家族分析。

研究结论翻译

虽然蛋白质界面保守氨基酸通常被假定具相似效应，但本研究表明Gα_q与Gα_i1亚型上某些Gα:Gβγ结合界面的保守氨基酸对Gβγ亚基偶联产生差异性影响，这种差异主要源于Gα亚基的动态特性。通过整合源自分子动力学(MD)模拟的三项互补特征——(i)残基接触频率、(ii)残基相互作用能及(iii)贝叶斯网络分析得出的残基对协同依赖性——研究人员证实Gα_i1与Gα_q的亚型选择性偶联源自残基接触时间持续性及焓与协同效应的复杂汇合。此三重分析揭示了仅凭动态构象系综分析方显现的亚型特异性界面特征。经预测—实验验证迭代，证明旁系同源蛋白中保守残基位置可引发不同功能结果。