基于标准人类基因组图谱的蛋白质互动网络：HuRI的整体结构与功能模块
针对人类蛋白质二元相互作用网络（Interactome）全局拓扑特征与功能模块映射关系尚缺乏系统量化评估的问题，本研究基于人类参考相互作用组（HuRI）图谱，在纯拓扑分析视角下对包含8,245个蛋白质及52,068对二元相互作用的网络进行了图论与复杂网络分析。研究证实HuRI网络具有重尾无标度分布、小世界特性及高模块化结构，并成功在无先验功能标签的前提下恢复出包括细胞骨架、跨膜运输、转录调控与剪接体组装在内的10大核心细胞功能模块，为理解人类细胞分子拓扑逻辑提供了基准图谱。
研究问题
蛋白质很少孤立发挥作用，绝大多数细胞功能由物理相互作用的蛋白质复合物及网络协同完成。解析蛋白质相互作用网络（PPI）的全局拓扑结构与功能模块，是网络生物学的核心目标。尽管已有多种局部或混合数据源的相互作用图谱，但在消除实验偏好与文献偏见的前提下，针对高质量、无偏的人类二元相互作用参考图谱（HuRI）开展系统性的网络图论解析仍有待深入。具体而言，该网络是否同时具备无标度（Scale-free）、小世界（Small-world）与高度模块化（Modular）的经典复杂网络特征？在排除节点度数分布的影响后，其聚类与路径特征是否依然显著？仅凭纯粹的网络拓扑结构，能否准确还原出人类细胞的主要生物学功能过程与核心控制节点？
研究路径
为回答上述问题，研究团队将HuRI数据库中的原始边列表构建为无向简单图（剔除480个自环后，包含8,245个蛋白节点和52,068对二元相互作用），并重点针对包含98.84%节点的最大连通分量（LCC，8,149个节点、52,016条边）展开计算。
研究路径分为四个步骤： 1. 全局拓扑与度分布分析：计算平均度、网路密度、全局群聚系数（Transitivity）及局部群聚系数，通过幂律拟合评估度分布的重尾特征。 2. 严格的小世界网络检验：构建100个保持节点度数不变的随机重构零模型（Degree-preserving rewired null models），采用Humphries-Gurney公式计算小世界系数 S，并评估真实网络与零模型在群聚系数与平均最短路径上的统计显著性（Z-score）。 3. 功能模块检测与一致性校验：分别运用Louvain与Leiden社区发现算法对LCC进行模块划分，计算模块化 Q 值及算法间一致性（调整兰德指数 ARI）。 4. 功能富集与枢纽节点识别：基于超几何分布检验，对前10大社区开展GO生物过程（Biological Process）富集分析（FDR <0.05）；同时利用Brandes算法计算介数中心性（Betweenness Centrality），识别调控网络流量的关键枢纽蛋白（Hubs）。
核心发现
研究结果显示，人类二元蛋白质相互作用网络在全局拓扑上呈现出高度协同的复杂网络规律，且拓扑模块与细胞生理功能高度重合。
- 网络高度连通且呈现重尾无标度分布。HuRI网络由单一巨大连通分量主导（覆盖98.84%的节点与99.90%的边），平均度为12.63，最大节点度达498。整体度分布表现出明显的重尾幂律特征（幂律指数 \gamma \approx 1.13），证实网络由少数高连接度的枢纽节点所统领。
- 零模型检验确立了显著的小世界拓扑架构。相比于100个保持度数不变的随机零模型，HuRI真实网络的全局群聚系数显著高于随机期望（Z = 64.4），而平均最短路径（3.839）仅比零模型略长9.5%（Z = 19.6），计算得出小世界系数 S = 1.60（显著大于1），严格证实了其小世界网络特性。
- 模块化分割高度清晰并精准映射核心细胞功能。Louvain与Leiden算法分别划分为65和58个社区，模块化 Q 值均超过0.51（显著高于随机网络的0.3），且两算法一致性良好（ARI = 0.62）。GO富集分析表明，前10大社区（覆盖77.9%的蛋白质）精准对应了细胞骨架组织（C0）、跨膜运输（C1）、转录调控（C2）、细胞外基质（C3）、代谢过程（C4）、体液免疫/泛素化降解（C5/C6）、Ras信号转导（C7）、纤毛运输（C8）及剪接体组装（C9，富集倍数达23.9倍）等核心生理过程。
- 介数中心性识别出掌控网络交通的关键瓶颈枢纽。介数中心性分析表明，MEOX2（最高介数0.094，度数498）、LNX1、CYSRT1、UBQLN2及TLE5等蛋白占据了网络最短路径的交汇点，作为模块间或子网间的桥梁节点（Bottlenecks），在细胞信息流动与稳态维持中发挥核心控制作用。
总结而言，本论文证明了仅依靠高质量二元相互作用图谱的纯拓扑结构分析，即可无偏地还原出人类细胞的各大核心功能分子机器与网络控制枢纽。
适用边界
本研究结论的推导与应用受限于特定条件。首先，HuRI数据集基于酵母双杂交（Y2H）技术构建，天然偏向于可溶性、核/质定位的二元直接相互作用，对于膜蛋白、瞬时相互作用及依赖多蛋白复合物的间接接触覆盖相对不足。其次，平均最短路径计算采用了1,000个随机节点对的采样估算，虽满足大型图计算需求并具备极高Z值，但仍存在轻微采样波动。此外，GO富集分析包含了电子推断（IEA）注释，且社区发现算法在微观边界划分上存在一定的随机性与算法敏感性。
基于现有证据，该研究能够支持在整体系统生物学视角下将HuRI作为人类二元相互作用的基准图谱，用于网络拓扑特征分析、模块功能预测及关键控制节点的优先筛选；但尚不足以证明特定拓扑社区或枢纽蛋白在所有动态生理条件、特定组织器官或疾病状态下均保持完全相同的物理结合与功能表达。