生物医学网络科学Biomedical Network Science
经典论文难度:入门

蛋白质网络中的致死性与中心性

Lethality and centrality in protein networks

Jeong H, Mason SP, Barabási A-L, Oltvai ZN · Nature 2001 · DOI: 10.1038/35075138

蛋白质相互作用网络酵母必需基因枢纽致死性–中心性规则

一句话

在酵母蛋白质相互作用网络中发现度分布重尾,且度越高的蛋白敲除后越可能致死,提出“致死性–中心性”关系。

相关专题:生物分子网络与拓扑特征、小世界与无标度网络模型、网络可控性与稳定性

一句话看懂

把酵母蛋白质相互作用网络与基因敲除表型对照:连接很少的蛋白绝大多数可以删除而不致死,连接很多的“枢纽”蛋白被删除后致死的比例高得多。一个蛋白在网络中的位置,能部分预测它对细胞存活的重要性。

研究背景与问题

传统细胞生物学按“单个分子的功能”(催化、信号、结构)来理解蛋白。后基因组时代的高通量酵母双杂交实验第一次给出了全蛋白组规模的相互作用图,引出一个新视角:蛋白的作用还取决于它在相互作用网络中的“语境”。与此同时,系统性的敲除实验表明酵母能耐受大量单基因缺失。

作者把两件事放在一起问:酵母蛋白质网络的整体结构是均匀的还是高度不均匀的?如果是后者,Albert 等(2000)在模型网络中发现的“耐随机故障、怕攻击枢纽”是否也体现在基因敲除的表型上?

核心思想与方法

这是一篇只有两页的短文,方法很简单:

  1. 构建网络:整合以酵母双杂交为主的系统实验数据和数据库数据,得到 1870 个蛋白、2240 条直接物理相互作用;没有方向信息,所有边按无向处理。最大连通分支约含 78% 的蛋白。
  2. 度分布:拟合带指数截断的幂律,截断位置约在 kc≈20k_c\approx20,即度超过 20 的蛋白略少于纯幂律的预期。
  3. 模拟扰动:随机删除蛋白几乎不改变网络整体拓扑;按度从大到小删除时,网络直径迅速增大。
  4. 与表型对照:从 Proteome 数据库取 1572 个已知表型的突变体,按度分组,计算每组中“单基因缺失即致死”(必需基因)的比例。

必需基因(essential gene):单独敲除即导致细胞死亡(在标准培养条件下)的基因。

致死性–中心性规则(centrality–lethality rule):蛋白在相互作用网络中的连接数越多,其编码基因越可能是必需的。

主要发现

  1. 度分布高度不均匀:酵母 PPI 网络的度分布是带指数截断的幂律,与幽门螺杆菌的蛋白网络和多种生物的代谢网络类似;少数高度连接的蛋白居于网络中心。
  2. 耐随机删除、怕删除枢纽:随机删除不影响整体拓扑;删除连接最多的蛋白使网络直径迅速增大。
  3. 度与致死性正相关:度不超过 5 的蛋白约占 93%,其中只有约 21% 是必需的;度超过 15 的蛋白只占约 0.7%,但约 62% 删除后致死。按作者的说法,高连接蛋白成为必需基因的可能性约为低连接蛋白的三倍。图 1c 的说明中报告必需蛋白比例与连接数的皮尔逊相关系数为 0.75。
  4. 作者据此推测,代谢网络和蛋白网络共享的重尾结构可能是进化选择的结果,酵母对突变的鲁棒性部分来自相互作用的组织方式。

怎样读这篇论文

全文就是一个长段落加一张三联图。重点看图 1:(a) 网络全图,节点按敲除表型着色(红色致死、绿色非致死),可以直观看到红色节点多集中在高连接区域;(b) 度分布,注意纵轴做了截断修正;(c) 必需蛋白比例随度的变化,这是全文的核心证据。读完后建议马上读 He 与 Zhang(2006)或 Zotenko 等(2008),体会同一相关性如何被不同机制解释。

局限与后续工作

  • 数据质量:早期酵母双杂交数据假阳性、假阴性都很高,不同来源的数据集重叠很少;被研究得多的蛋白测到的相互作用也更多,度本身可能带有研究偏倚。
  • 机制解释之争:作者倾向于“删除枢纽破坏网络连通性”的解释。后续提出了多种竞争解释:He 与 Zhang(2006)认为少数相互作用本身是必需的,度越大越可能包含这样的相互作用;Zotenko 等(2008)认为必需性主要来自必需蛋白复合物(如核糖体),复合物成员天然度高。还有研究在部分高质量数据中只观察到较弱的相关。
  • “枢纽即靶点”的推论需谨慎:Goh 等(2007)发现大多数人类遗传病基因并不是枢纽,枢纽的必需性与“致病性”是两回事。
  • 方法学启示:稳健的相关可以有多个机制解释,区分它们需要设计让不同解释给出不同预测的分析,例如控制复合物成员身份后再看度的预测力。

关联知识

延伸思考

  1. 如果相互作用数据存在“被研究得越多、测到的边越多”的偏倚,而被研究得多的蛋白又往往是重要蛋白,致死性–中心性相关会被如何放大?你会设计什么分析来检验这一点?
  2. 请为“必需相互作用”解释和“必需复合物”解释各提出一个可检验的预测,并说明需要什么数据。
  3. 度只是一种中心性。如果改用介数或 k-core 指数,你预期与致死性的相关会更强还是更弱?为什么?