通过进化引导的残基谱分析解码G蛋白偶联受体（GPCR）的功能特化

《British Journal of Pharmacology》：Decoding functional specialization in G protein-coupled receptors (GPCRs) through evolution-guided residue profiling

　　

背景与目的：G蛋白偶联受体（GPCR）是整合膜蛋白，通过使细胞能够检测并响应多种刺激来介导生理过程。尽管许多亚家族特异性功能热点已被描述，但家族范围内共同与亚家族特异性功能的决定因素仍不完整。实验方法：本研究开发了一种进化框架，利用直系同源物内的保守性和旁系同源物间的变异，将位置分类为常见残基（CR）或选择性残基（SR）。关键结果：常见残基（CR）聚集在与结构稳定性和激活相关的位点，而SR集中在参与配体和转导子结合的选择性界面。不同家族的SR分布表明，某些类别主要通过配体识别位点的改变实现多样化，而其他类别则通过转导子结合界面的改变实现多样化。研究人员还发现了参与家族特异性和跨家族基序的CR，包括保守的二硫键和胆固醇接触位点。结论与意义：这些发现共同提供了家族范围特征的进化蓝图，强化了已知关联，并提供了可检验的假设，尤其对研究不足的家族具有重要价值。

# 论文解读：通过进化引导的残基谱分析解码GPCR功能特化

## 研究背景与问题

G蛋白偶联受体（GPCR）是最大的跨膜（TM）蛋白超家族之一，通过将外部刺激转化为细胞反应，在细胞通讯中起关键作用。人类GPCR传统上分为七大家族（A、B1、B2、C、F、味觉2和嗅觉受体），各家族具有独特的功能和进化轨迹。尽管已有研究描述了亚家族特异性功能热点，但家族范围内决定共同与特化功能的机制仍不完整。例如，Flock等人（2017）揭示了G蛋白偶联选择性的进化基础，但主要强调G蛋白侧，受体侧的选择性驱动机制未充分探索。此前的研究多集中于特定受体亚群或区域，缺乏对GPCR功能进化全局原则的理解。为此，研究人员开发了一种系统进化框架，通过分析直系同源物内的保守性和旁系同源物间的变异，鉴定家族水平的功能热点，并揭示进化如何塑造不同家族的共同与特化功能。该研究发表在《British Journal of Pharmacology》。

## 主要技术方法概述

研究人员首先利用隐马尔可夫模型（HMM）搜索从UniProtKB人类蛋白质组中检索GPCR序列，共获得837个人类受体序列。随后，通过BLAST搜索鉴定每个人类受体的直系同源物，并构建多序列比对（MSA）和最大似然（ML）系统发育树。采用自定义算法剔除高度分化的序列，保留真实直系同源物。基于直系同源MSA计算每个位点的保守性百分比，并设定分层阈值（≥90% for >50 sequences, 95% for 25–50 sequences，排除<25个直系同源物的受体）。对于每个家族，通过CD-HIT聚类（0.65同一性阈值）生成家族MSA，计算家族水平保守性（平行同源物中保守位点的比例）和香农熵（Shannon entropy）。根据熵值，将家族内高度保守的位点分为常见残基（CR，低熵）和选择性残基（SR，高熵）。阈值通过手动优化，结合已知功能基序、结构聚类和最小残基集标准。此外，利用临床变异（ClinVar）和六个深度突变扫描（DMS）数据集进行基准验证，并计算Pearson相关系数。

## 研究结果

### 3.1 A类受体
A类GPCR（293个受体）作为基准验证。CRs富集于受体核心，包括关键脯氨酸、保守二硫键、CWxP、DRY、NPxxY、PIF/PIW基序和Na⁺口袋等，与结构稳定性和激活相关。SRs主要位于配体结合口袋和转导子结合界面，其中配体结合位点熵值显著高于转导子结合位点，表明配体识别驱动了更大程度的多样化。临床变异和DMS数据验证显示，家族水平保守性与突变不耐受性正相关（r=0.815），熵值与突变不耐受性负相关，且常见激活网络残基（Zhou et al., 2019）中CRs占主导。与先前进化追踪方法比较，本方法更好地区分了CRs和SRs。

### 3.2 嗅觉受体
嗅觉受体CRs聚集于三个主要区域：受体核心与转导子界面（共享Golf激活机制）、N端尾与ECL1/ECL2交叉处（独特胞外结构）、以及一对相互作用残基（1x43和2x58）。保留了A类部分基序（如NPxxY、DRY、TM3-ECL2二硫键），但缺失PIW和CWxP基序，且Na⁺口袋中S3x39被天冬氨酸取代。SRs主要位于TM3-4-5和ECL2区域，表明配体结合口袋形状和位置与A类不同，与首个嗅觉受体实验结构一致。

### 3.3 味觉2受体（T类）
T类受体CRs聚集于转导子结合界面、受体核心和胞内环2（ICL2），ICL2形状独特，与G_gust偶联相关。SRs分布于配体结合口袋、胞内界面周边及TM7与螺旋8之间。保留了PIF基序中的P5x50和F6x44、NPxxY中的脯氨酸等，但Na⁺口袋2x50位点带正电，抑制钠离子结合。

### 3.4 B类受体
B1受体CRs富集于PxxG、NPGQ、HETX基序、二硫键、配体结合口袋、受体核心和转导子界面，包括四个与A类同源的胆固醇结合位点（3x38、4x50、4x49¹、4x53）。SRs主要位于配体结合口袋和转导子界面周边，与B1主要偶联G_s一致。B2（粘附受体）CRs位于GAIN结构域、二硫键、胆固醇结合位点、受体核心和束缚激动剂结合口袋，而SRs集中于转导子结合界面和GAIN结构域周边，与其偶联多种G蛋白的特性相符。此外，识别了B1和B2中与RAMP1结合的SRs，提示其在RAMP选择性中的作用。

### 3.5 F类受体
F类受体包含SMO和10个FZD。CRD结构域中CRs包括八个形成二硫键的半胱氨酸和两个脯氨酸，SRs富集于CRD-WNT界面，尤其指状界面比拇指界面有更多SRs，决定WNT选择性。TM结构域中SRs位于胞外环，CRs包括YHW、WGM基序和R/K6x32–W7x55开关。胆固醇结合位点（W4x50及其它四个芳香族残基）为CRs，与近期实验验证一致。胞内区域CRs位于转导子结合界面中心，SRs位于周边，包括螺旋8。

### 3.6 C类受体
C类受体（含VFT结构域）CRs在VFT结构域核心和配体结合位点富集，SRs集中于配体结合口袋和上叶相互作用网络。CRD中八个半胱氨酸形成四个二硫键（零熵），SRs位于二聚化界面。TM结构域中CRs包括ECL2-TM3二硫键、变构调节位点（6x50、6x53）和胞内盐桥网络（3x50、6x35、7x51、12x53）。SRs富集于变构调节位点、TM3-TM5界面、胆固醇结合位点（与其它家族不同）以及TM7和C端区域，可能参与转导子选择性。

## 讨论与结论

讨论部分总结了进化原则：CRs富集于共同功能位点，SRs定位于配体结合和转导子偶联的选择性界面。这一模式源于基因复制后，新旁系同源物承受强纯化选择以保留基本功能（CRs），而选择性界面经历放松和谱系特异性选择（SRs），导致配体和偶联偏好多样化。A、B1、F、C、T和嗅觉受体在配体结合界面有大量SRs，而B2受体因束缚激动剂变化小，其配体口袋富集CRs。C类受体在TM变构调节位点意外发现SRs，提示可能存在内源性变构调节剂。转导子结合界面方面，A和B2类受体SRs富集程度更高，如β-肾上腺素受体与α2-肾上腺素受体具有不同的G蛋白相互作用SRs。

跨家族比较显示，除C类外所有家族共享TM3-ECL2二硫键；胆固醇结合位点W4x50在A、B、F和嗅觉受体中部分或完全保守；脯氨酸（P5x50、P6x50、P7x50）和芳香族开关W/Y6x48在多个家族中保守，支持所有GPCR的共同进化起源。

研究局限性包括：自动直系同源鉴定可能包含未有效过滤的差异序列；F类受体数量少限制了分辨率；CRs和SRs的阈值具有主观性，不应视为严格生物学类别。

结论：这一进化蓝图系统揭示了家族范围内的功能特征，强化了已知关联，并提供了可检验的假设，尤其对研究不足的家族具有重要价值。通过靶向SRs设计药物可提高特异性并减少脱靶效应，而SRs也可作为别构调节剂的新结合位点。数据公开于gpcrevolution.org，以供社区使用。