Survey of nucleotide-specific Rab GTPase interactions reveals multiple Rab effectors

技术摘要：核苷酸特异性 Rab GTPase 相互作用调查

问题陈述
Rab GTPase 是作为核苷酸依赖性开关发挥作用的，它们调节特定亚细胞区室内的蛋白质-蛋白质相互作用。它们的活性受结合 GTP（活性状态）与结合 GDP（失活状态）时“开关”区域构象变化的影响，这种变化决定了它们对下游效应因子的招募。虽然 Rab 蛋白的一般功能已得到理解，但其核苷酸特异性相互作用组的全面图谱仍不完整。以往的高通量筛选往往无法区分 GDP 结合态与 GTP 结合态，且许多现有的酵母双杂交（Y2H）文库存在偏向于 3' 非翻译区或缺乏足够的蛋白质编码结构域覆盖的问题。此外，在相互作用筛选中使用全长蛋白可能会掩盖由分子内遮蔽介导的离散结构域间的相互作用。因此，完整的 Rab 效应因子库，特别是那些需要特定核苷酸构象进行结合的效应因子，尚未被完全阐明。

方法论
作者采用了一种改进的大规模 Y2H 方法，称为 DEEPN（用于评估蛋白质网络的动态富集），并结合了定量高通量测序。

• 文库构建： 为了克服传统 cDNA 文库中的偏差，作者利用 Human ORFeome v8.1 构建了一个高密度的 Y2H 捕获（prey）文库。他们扩增了开放阅读框（ORF），通过酶促片段化处理以避免剪切伪影，并将它们连接到 Gal4 激活结构域下游的捕获载体中。这 resulting 的文库包含约 103 万个独特元件，其中约 19% 编码框内编码序列（CDS）片段，相比之下，商业 cDNA 文库的框内 CDS 覆盖率仅为约 6.7%，这是一个显著的提升。
• 诱饵（Bait）设计： 本研究利用 36 种人类 Rab GTPase 作为诱饵。每种 Rab 都被工程化为两种不同的构象：GTP 锁定状态（通过 Q>L 或 S/T>N 突变）和 GDP 锁定状态。C 末端 CAAX 异戊烯化基序被突变为丙氨酸尾部，以防止膜结合，从而确保其在酵母中的胞质表达。
• 筛选与分析： 将 Rab 诱饵在批次液体培养中针对 Human ORFeome 文库进行筛选。在选择前（缺乏组氨酸的培养基中生长）和选择后回收捕获质粒，并进行双端 Illumina 测序。
• 数据处理： 作者利用更新后的 DEEPN 软件套件（已移植至 Python 3/PyQt5 并兼容 Apple Silicon）来分析测序数据。工作流程包括：
  • 将读取序列（reads）比对至 hg38 基因组。
  • 识别捕获基因片段并利用接头分析确定其阅读框。
  • 使用结合贝叶斯后验富集概率（>0.8）和基于 DESeq2 的 log2 倍数变化（>3 倍富集，p ≤ 0.01）的统计学命中调用（hit-calling）。
  • 进行严格过滤，以去除“仅含载体”的背景相互作用物以及比对至不可及区域（如膜蛋白的胞外结构域）的片段。
• 验证： 通过二元 Y2H 点样实验、GST 拉下实验（pull-down）以及核磁共振波谱（HSQC）验证候选相互作用，以确认直接的物理结合。

关键结果
本研究生成了一个包含 527 个高置信度 Rab-捕获物相互作用的数据集，涵盖了所有 36 种 Rab GTPase 的 337 个独特捕获基因。

• 核苷酸特异性： 大多数识别出的相互作用特异于 GTP 结合构象，证实了 Rab 与效应因子结合的核苷酸依赖性质。
• 多 Rab 效应因子： 筛选识别出了几种能与多种 Rab GTPase 相互作用的蛋白质：
  • MICAL-L1： 通过其 bMERB 结构域，证实了与 Rab8、Rab10、Rab13、Rab15 以及新发现的 Rab12 的相互作用。
  • Optineurin (OPTN)： 验证了已知的 Rab8 相互作用，并发现了与 Rab1、Rab10、Rab35、Rab38 和 Rab43 的新型 GTP 依赖性相互作用。片段映射表明，其 N 端结构域内存在针对不同 Rab 伙伴的截然不同的结合表面。
  • RABEP1 (Rabaptin-5)： 证实了 Rab4 和 Rab5 具有不同的结合区域。值得注意的是，Rab4 结合区域被发现能与更广泛的 GTP 结合态 Rabs（包括 Rab8、Rab10、Rab13、Rab33A、Rab35 和 Rab43）发生相互作用，这暗示了一种内体成熟的反馈机制。
• 新型 Rab 结合结构域： 一个重要的发现是鉴定出 SNX13 和 SNX14 的 RGS（G 蛋白信号调节因子）结构域为直接的 Rab 结合模块。
  • SNX13 RGS： 以 GTP 依赖方式结合 Rab5，并以 GDP 和 GTP 两种构象结合 Rab9。
  • SNX14 RGS： GTP 依赖性地结合 Rab11、Rab32、Rab9 和 Rab38。
  • 直接结合通过 GST 拉下实验和 NMR 化学位移扰动得到了确认。
• 调节与疾病相关性：
  • 磷酸化： SNX14 中的磷酸化模拟突变（S382D/S388D）选择性地破坏了与 Rab32 的结合，同时保留了与 Rab11 和 Rab9 的相互作用，表明磷酸化可作为伙伴选择的调节开关。
  • 疾病突变： 与 SCAR20 相关的错义变体 E370K 在 SNX14 中被证明会消除 Rab11 结合，但保留 Rab9 相互作用。尽管该突变降低了蛋白质稳定性，但这种选择性的丢失将特定的 Rab11-SNX14 相互作用与脊髓小脑共济失调的病理生理学联系起来。
  • 定位： SNX13 的 RGS 结构域被证明对于其向 Rab5 阳性内体的招募至关重要，从而将 RGS-Rab 相互作用与亚细胞定位联系起来。

意义与主张
本文声称，该 DEEPN 数据集提供了一个高分辨率、具有核苷酸特异性的 Rab 相互作用组图谱，其在特异性和结构域水平的分辨率上都超越了以往的研究。通过使用基于片段的文库和严格的统计过滤，作者成功地区分了真实的、具有诱饵特异性的相互作用与背景噪声。

研究强调，许多 Rab 效应因子拥有比此前认为的更广泛的伙伴库，通常通过特定结构域与多种 Rab 同工型发生作用。识别出 SNX13 和 SNX14 的 RGS 结构域作为 Rab 结合模块，表明存在一种趋同进化机制，即通常与 Gα 亚基相关的 RGS 结构域也调节着 Rab GTPase。这一发现架起了 ER 定位的分选连接蛋白（sorting nexins）与内体区室之间的桥梁。

作者强调，尽管某些相互作用仍需进一步的生物学验证，但该数据集为定义 Rab 介导的细胞过程的机械机制提供了强大的资源。这项工作表明，理解 Rab 生物的全貌仍然是不完整的，且核苷酸状态特异性的筛选对于揭示这些 GTPases 的功能作用至关重要。更新后的 DEEPN 软件以及人类 ORFeome 片段文库的可用性，被作为促进蛋白质相互作用网络进一步探索的易用工具进行展示。