起始甲硫氨酸加工在核糖体组装与功能中的候选角色：一项跨物种普查

《Life》：Candidate Roles for Initiator Methionine Processing in Ribosome Assembly and Function: A Cross-Species Census

　　

起始甲硫氨酸加工如何与核糖体蛋白界面相关联仍不完全清楚。本研究整合了来自61个物种的68个核糖体结构（1873条家族–物种记录，95个家族），对末端及已发表的组装中间体进行了表征，以考察比较性关联并识别特定的候选界面。残基2的预测与12个物种中393个已表征末端中的388个一致。在58个可切换家族中，基于埋藏程度的RNA-only与暴露比值，预测保留的比值比（odds ratio）为0.31（95%置信区间0.18–0.53，按家族分层）、0.64（0.32–1.29，按家族×结构域分层）及0.76（0.29–1.99，按家族×进化支分层）；重建的家族×进化支估计值为1.21（0.48–3.07）。共同位置评分保留了合并的RNA关联，而匹配的系统发育估计仍不具结论性，蛋白质接触对比则强烈依赖于所评分的位置。一项针对21,485个参考蛋白质组的探索性分析在目（order）水平（所测试的五个层级之一）识别出细菌关联，该关联在结构域内多重校正后仍显著，但在跨三个结构域的校正后不再显著。保守的允许性家族在RNA接触组与蛋白质堆积组之间未检测到差异，且观察到的bS16保留与统一切割要求相悖。在单个界面上，uL15和eL24表现出替代性末端接触寄存器，而酵母uS12和eS27在pre-40S及成熟颗粒中分别占据受支持的RNA界面和蛋白质界面。这些比较通过识别特定的分子界面及其组装背景，扩展了与加工相关的功能研究。它们优先将uS12和eS27列为研究末端加工如何影响核糖体组织与功能的候选对象；各位点的功能后果仍有待确立。

论文解读：起始甲硫氨酸加工在核糖体组装与功能中的候选角色——一项跨物种普查

一、研究背景与问题

蛋白质合成通常以甲硫氨酸起始，甲硫氨酸氨肽酶（MetAP）可在第2位残基为小而中性氨基酸（Ala、Gly、Ser、Thr、Cys、Pro或Val）时共翻译切除起始甲硫氨酸。这一序列规则预测的是底物易感性而非完全加工。MetAP活性在大肠杆菌中是必需的，两种酵母酶联合必需。在真核生物中，切除可暴露底物供NatA依赖性乙酰化或Gly2供N-肉豆蔻酰化，乙酰化末端还可参与组装依赖性蛋白质质量控制。这些已确立的途径将加工与蛋白质功能联系起来，但未识别哪些核糖体界面可能对末端变化本身敏感。此前一项未经过同行评审的人类生物学组装普查发现，在1123个可分析蛋白质中，允许性残基2类别与预测保留甲硫氨酸的蛋白质相比，蛋白质界面参与频率较低（37.5%对47.5%；校正比值比0.65），但该比较涉及不同蛋白质且排除核酸接触。Nagai等人将酵母MetAP募集受损与选定核糖体蛋白上起始甲硫氨酸保留增加联系起来，并将翻译保真度表型定位到40S亚基，但负责的蛋白质或组合仍未解决。本研究通过整合序列衍生加工类别、末端观察、局部结构支持和已发表的组装状态，在跨物种同源核糖体家族中比较加工与末端环境之间的关联，并识别具体候选界面。

二、研究内容与结论

研究人员整合了来自61个物种的68个核糖体结构，形成1873条家族–物种记录、95个家族。预测与12个物种中393个已表征末端中的388个一致（98.7%）。在58个可切换家族中，基于埋藏的RNA-only与暴露比值，预测保留的比值比为0.31（家族分层），但按家族×结构域分层后衰减至0.64，按家族×进化支分层为0.76；重建的家族×进化支估计为1.21。共同位置评分保留合并的RNA关联，但匹配的系统发育模型估计不具结论性（RNA-only对比暴露：成熟起始评分0.62，编码位置2评分0.58）。蛋白质接触对比强烈依赖于评分位置：蛋白质接触对比暴露从1.98变为0.53，蛋白质接触对比RNA-only从3.15变为0.90。对21,485个参考蛋白质组的探索性分析在细菌目水平识别出关联（位置特异性对比+22.0个百分点；置换p=0.0040/0.0019），在结构域内校正后显著（p=0.036/0.019）但跨三个结构域校正后不显著（p=0.117/0.057）。保守允许性家族在RNA接触组与蛋白质堆积组之间未检测到差异（基因组允许性分数91.3%对87.2%；Mann–Whitney p=0.32），且bS16在大肠杆菌、嗜热栖热菌和耐辐射奇球菌中观察到的Met-Val保留直接反驳统一切割要求。在单个界面上，uL15和eL24表现出替代性末端接触寄存器：uL15的Met1或后段亮氨酸占据同源疏水口袋，eL24的Ala/Pro末端接近酸性伴侣残基。酵母组装状态分析显示，uS12 Gly2在pre-40S和成熟颗粒中距RNA 2.38–3.24 ?而距蛋白质8.64–9.36 ?，eS27则距蛋白质2.85–3.91 ?而距RNA 4.78–6.77 ?，且eS27与eS21的受支持接触在6RBD、6WDR和成熟4V88中均存在。这些发现优先将uS12和eS27列为研究末端加工如何影响核糖体组织与功能的候选对象。

三、主要关键技术方法

本研究为观察性整合研究，未开展新生物学实验。主要方法包括：从RCSB蛋白质数据银行（PDB）检索并选定68个胞质核糖体结构（23细菌、10古菌、28真核生物），使用GEMMI读取坐标、SIFTS映射链、UniProt核对序列；按通用命名法分配家族并手工校正13条链；依据编码残基2分类MetAP允许性（A/G/S/T/C/P/V，预测切除）或甲硫氨酸保留（其他残基，预测保留）；使用FreeSASA测量埋藏程度，以≥25%蛋白质埋藏定义两态堆积标准，主三态分类器按4.0 ?距离阈值区分蛋白质接触、RNA-only和暴露；采用Mantel–Haenszel比值比及条件似然进行分层关联分析，使用PyMC构建带家族×结构域和系统发育相关物种截距的逻辑回归模型；对UniProt参考蛋白质组进行全基因组序列分析，采用整块置换检验评估五个分类层级；对七个酵母沉积结构进行全局比对和距离测量，以EM Q分数和包含度评估局部支持；对PRIDE PXD044467质谱数据库进行只读检查以审计末端证据。

四、研究结果

4.1 普查覆盖与观察到的末端加工：原始普查含1873条坐标合格记录，89个家族出现在至少三个物种中，其中58个切换预测加工类别、25个全程允许、6个全程预测保留。预测与388/393个独立表征末端一致，所有279个Ala/Gly/Pro/Ser末端均被加工，所有82个非允许性末端保留Met，Thr2加工15/16、Val2加工12/16。例外包括大肠杆菌、嗜热栖热菌和耐辐射奇球菌中的bS16 Met-Val，耐辐射奇球菌中的uS10 Met-Val，以及大肠杆菌中的bS1 Met-Thr。

4.2 家族内关联取决于末端环境：原始两态分析中，允许性记录12.5%蛋白质堆积，预测保留记录26.0%堆积，家族分层比值比0.34。纳入核酸埋藏逆转关联至1.73，促使分离蛋白质、RNA-only和暴露类别。在1351条切换家族记录中，预测保留比例在蛋白质接触（45%）、RNA-only（18%）和暴露（41%）间不同。家族分层保留比值比为蛋白质接触对比RNA-only 3.87、RNA-only对比暴露0.37、蛋白质接触对比暴露1.48。

4.3 分类学、评分位置与坐标合格性限定关联：原始RNA-only对比暴露估计从0.31（家族分层）衰减至0.64（家族×结构域）和0.76（家族×进化支）。共同编码位置2评分在1342对中给出0.33、0.55和0.33。骨干原子评分给出更强的家族×结构域关联（0.13），而末端氮评分（0.77）和前三残基评分（0.94）不具结论性。匹配系统发育模型给出RNA-only对比暴露0.62（成熟起始）和0.58（编码位置2），蛋白质接触对比暴露从1.98变为0.53。结构域条件估计显示细菌0.482、古菌1.689、真核生物0.146，异质性名义p=0.066，校正后不显著。

4.4 八个额外物种的回顾性比较：额外样本提供192条记录，两态物种分层比值比0.25，距离分类RNA-only对比暴露0.25（去除微孢子虫后0.17），蛋白质接触对比暴露不具结论性（1.36）。

4.5 基因组转移关联具有层级特异性和多重性敏感性：细菌目水平严格分析含73,364条记录、6924个基因组、18个家族、201个块。RNA标记进化支在残基2允许性上超出其他标记进化支19.9个百分点，位置特异性对比+22.0点。整块置换p=0.0040/0.0019，十检验校正后p=0.036/0.019，但三十检验跨结构域校正后p=0.117/0.057。家族水平效应未通过校正，古菌和真核生物结果稀疏。

4.6 同源蛋白质界面容纳替代性末端寄存器：4038个末端伴侣接触中，uL15–uL4、uL6–eL20、eL24–uL3等伴侣身份跨加工类别维持。uL15在保留细菌中将Met1置于同源疏水口袋，而鲍曼不动杆菌和柑橘黄单胞菌将后段亮氨酸置于该处；eL24在古菌Ala/Pro末端接近人类Met1接触的酸性伴侣残基。

4.7 保守类别与已发表加工响应优先化特定候选：25个原始保守允许性家族跨RNA-only、蛋白质接触和暴露环境。uS12在41个物种中使用六种允许性残基，扩展后49/52物种允许；eS27为29/31。Nagai等人突出显示的九个40S蛋白质中，uS12、eS17、eS27和uS4具有建模成熟界面。未修饰保留uS12肽段在ES27L和map1比较中均增加，而eS27同一保留肽段在ES27L中比值为46.37/43.179但在map1中为0.945。RNA接触组与蛋白质堆积组的基因组允许性分数无显著差异。

4.8 uS12和eS27的不同界面在成熟前后均有代表：七个酵母条目中，uS12 Gly2距RNA 2.38–3.24 ?而距蛋白质8.64–9.36 ?，局部支持在6FAI、6RBD、6WDR和两个成熟拷贝中通过。eS27距蛋白质2.85–3.91 ?而距RNA 4.78–6.77 ?，受支持接触涉及eS21。eS17早期模型Q分数低，仅6WDR和一个成熟拷贝通过主要筛选。eS24提供暴露对照，模型化起始距RNA 9.59–12.22 ?、距蛋白质8.01–13.57 ?。

五、讨论与结论

本研究最清晰的贡献是识别出赋予核糖体加工分子背景的特定末端界面和替代性结构排列。uS12和eS27将已发表的加工响应与pre-40S和成熟颗粒中受支持的RNA和蛋白质接触联系起来；uL15和eL24证明同源蛋白质口袋容纳不同末端残基和寄存器。比较普查暴露测量依赖性：合并RNA-only关联在共同位置评分下持续，但匹配家族×结构域和系统发育区间不具结论性；蛋白质接触对比在相同编码位置评分下显著改变。一般性陈述某一环境独立决定加工将超出这些结果。uL15和eL24的替代性寄存器显示局部序列和伴侣化学的重要性，但自然发生的排列支持跨序列背景的结构容纳，既不显示保留Met在每个允许性同源物中都会破坏性，也不确立容纳会保持功能。Nagai等人的对照排除简单NatA介导解释，但肽段丰度、末端占据和功能因果仍是不同测量。已确立的生物发生途径提供具体背景：人类耗竭研究将uS12和eS24置于早期pre-rRNA加工类别，裂殖酵母中uS12 N端区域参与组装前与Ofd1的关联。多项观察限定范围：加工基准覆盖12个物种且含五个例外，直接bS16保留与组装细菌核糖体相容，RNA接触对比蛋白质堆积的保守家族比较不具结论性，基因组转移具有层级特异性和多重性敏感性。MetAP必需性涵盖多个底物和途径，NatA依赖性质量控制和N-肉豆蔻酰化提供从末端加工到蛋白质可用性的已确立联系。疾病和生物体表型确立蛋白质的重要性而非加工特异性疾病机制。一个聚焦的后续问题是改变末端加工是否改变已识别的uS12–RNA或eS27–eS21相互作用，进而影响亚基组织或性能。本研究所组装的结构和加工证据使这些蛋白质特异性问题变得可处理。