生物网络结构与网络模体¶

细胞中的分子很少孤立工作。转录因子结合调控区，受体把细胞外信号传入磷酸化级联，酶把一种代谢物转为另一种，蛋白质组装成复合物，两个基因的联合扰动又可能产生单独扰动无法预期的表型。把这些对象画成节点、把它们之间的关系画成边，便得到生物网络。网络的价值不只在于把许多名称放在一张图上，而在于保留关系的方向、符号和条件，从局部连接中寻找反复出现的回路，再检验这些回路怎样影响响应时间、噪声、记忆和资源分配。

网络结构是系统行为的骨架，具体行为还取决于结合强度、降解速率、逻辑门和初始状态；同一个三节点回路因而可以表现为延迟、脉冲或近似静态。同一条数据库边也可能代表直接接触、共同复合物或统计相关。节点—边语义、全局结构与局部网络模体构成这里的主要框架；速率方程、稳定点、双稳态和振荡见反馈、开关与振荡，网络数据的构建、数据库和算法分析见生信栏目的生物网络与系统生物学。

生物关系的图结构¶

一个网络可写作 \(G=(V,E)\)，其中 \(V\) 是节点集合，\(E\) 是边集合。节点可以是基因、转录本、蛋白质、蛋白质的活化状态、代谢物、反应、细胞或物种；边可以是结合、催化、激活、抑制、转运、遗传互作或信息传递。选择什么作为节点，实际是在规定系统的分辨率。把“未磷酸化的蛋白”和“已磷酸化的蛋白”压成同一节点会简化图形，也可能丢掉信号传播所依赖的状态转换；把代谢反应投影成代谢物两两相连，则会把三磷酸腺苷（adenosine triphosphate，ATP）、水和质子等高频参与者画成不具特异性的巨大枢纽。

无向边表达对称关系，有向边从调控者、反应物或发送者指向靶标、产物或接收者。边的正负号可区分激活与抑制、正遗传互作与负遗传互作；权重可表示速率、占有率、效应强度或证据置信度。多重边允许同一对节点同时具有物理结合、遗传互作和共同表达等不同关系，多层网络则把转录、信号、代谢和细胞间作用分别放置，再以共享分子或状态连接各层。节点名称相同并不保证各层中的生物学对象相同，例如一个基因节点、其蛋白产物和一个特定磷酸化蛋白节点属于不同层次。

网络还必须带有语境。细胞类型、发育阶段、亚细胞区室、营养条件、刺激强度和时间都会改变哪些节点存在、哪些边活跃。转录因子在某处有结合位点，说明一条可能的调控关系；在特定细胞中占据该位点、随后改变靶基因转录，才把这条边限定到当前条件。静态图通常是多次测量的汇总，它适合表示可能发生的关系，而实时发生的信号流需要由时间序列、分子状态和扰动实验进一步确定。

五类分子网络的边语义¶

不同网络可以共享相同的圆点和连线画法，边所陈述的生物学命题却很不一样。¹

转录调控网络最适合讨论有向、有符号的局部回路。转录因子可以调控自身，也可先调控第二个转录因子，再由二者共同调控靶基因。真核调控还包含增强子、染色质状态、共调节因子和远距离三维接触；把它们压成“转录因子—基因”一条边，可突出回路结构，却不会自动说明具体启动子逻辑。天然基因表达调控的分子机制见真核基因表达调控。

信号网络中的节点往往需要写成分子状态。激酶被磷酸化、GTP 酶结合三磷酸鸟苷（guanosine triphosphate，GTP）、转录因子进入细胞核，分别是同一分子的不同状态；激活和抑制也可能通过多步反应间接发生。受体向多条支路发散，使同一输入影响代谢、细胞骨架和转录；多条支路又可汇合到同一激酶或转录因子。支架蛋白和细胞器膜把反应物限制在局部，空间邻近因而与拓扑连通同样重要。

代谢网络天然包含反应方向和化学计量。一条反应同时消耗多个底物并产生多个产物，最忠实的图表示通常是“代谢物—反应”二部图或超图。糖、脂肪酸和氨基酸等多样输入通过较少的中心中间体汇合，再向能量、还原力和生物合成前体等多样输出发散，形成近似“蝴蝶结”结构；这种结构让许多输入共享核心反应，也使核心节点的容量与调控成为系统约束。² 具体反应与通路见代谢总览，通量平衡和代谢控制分析归入后续系统建模页面。

物理互作网络中的边可以是直接接触，也可以是共同复合物或一定时空范围内的邻近。一个稠密子图常提示稳定复合物，但膜蛋白、短暂互作和低丰度蛋白可能因实验方法而缺边。遗传互作则记录双扰动对表型的非独立效应：负互作可以表现为比期望更严重的缺陷，正互作可表现为抑制或缓冲。遗传互作指向功能依赖，不能据此推断两个蛋白直接结合。两类网络叠合时，物理上接近且遗传上互作的节点可为同一功能模块提供互补证据。

稀疏性、异质连接与可达路径¶

含有 \(n\) 个节点的有向网络理论上可有近 \(n(n-1)\) 条边，真实细胞网络只实现其中一小部分。结合界面的化学特异性、基因调控位点、亚细胞区室、表达成本和进化历史共同造成稀疏连接。稀疏并不意味着彼此孤立：少数分支点、共享中间体和跨模块连接可以让扰动在较短路径上到达许多节点，同时保留大量局部专一关系。

节点度数表示直接相连的边数；有向图还区分入度和出度。高出度转录因子可协调许多靶基因，高入度启动子可整合多种信号；蛋白互作图中的高度节点可能是多亚基复合物的共同成员，代谢图中的高度代谢物也可能只是 ATP 或水一类广泛参与者。“枢纽”因此是相对于节点和边定义得到的结构属性，是否为功能瓶颈还要看边的方向、替代路径、分子丰度和扰动后果。

早期网络生物学常用幂律度分布和“无标度网络”描述少数高度节点与大量低度节点并存的情形，并由基因复制—分化或偏好连接解释其形成。有限、缺边且带测量偏倚的生物网络可由多种重尾分布产生相似图形。对近千个现实网络的系统比较显示，严格无标度结构远非普遍规律；生物网络可以具有显著度数异质性，却仍需比较幂律、对数正态等候选分布并报告取样范围。³ 度数分布也不足以单独推出网络的形成机制。

平均最短路径和聚类系数常用于描述“小世界”特征：全网节点之间的路径较短，而局部邻居又较容易彼此连接。投影方式会强烈改变这些量。把一个十亚基复合物投影成十个蛋白两两相连，会立即制造高聚类；把可逆代谢反应中的所有底物产物相连，会缩短许多并不存在直接物质流的路径。最短路径先是图上的可达关系，只有当沿途边的方向、符号、区室和时间能够连续组合时，才可成为候选生物通路。

网络模体与零模型¶

网络模体是相对于某个随机网络集合显著富集的小型连接模式。经典研究先在大肠杆菌转录调控网络中发现前馈环、单输入模块等重复模式，随后把三至四节点子图与保持网络基本特征的随机图比较，并将显著富集的模式称为模体。⁴⁵ 对第 \(i\) 种子图，可用

表示真实网络中的计数相对随机集合均值偏离了多少个标准差。\(Z_i\) 的意义完全依赖随机集合如何生成。只保持节点数和边数的均匀随机图，容易把由度数异质性造成的常见子图误作特殊设计；更严格的随机化通常保持各节点的入度和出度，有时还需保持层级、边符号、区室、时间或可检测边集合。算法怎样枚举和随机化子图属于生物网络与系统生物学的实现范围。

统计富集说明一种连接模式在所选背景下比预期常见，不会自行给出它的功能。一个前馈环可能来自共同祖先基因的复制，也可能由调控位点逐步获得；同样拓扑在不同参数下会有不同动态。把模体解释为功能单元，需要依次确定边的分子机制、测量各节点的时间响应、改变特定边或输入逻辑，并检验预言的延迟、脉冲或稳态是否随之改变。网络模体概念的强项正在于把大图缩成可实验操纵的局部假说。

自调控回路¶

转录因子调控编码自身的基因时形成自调控。负自调控中，蛋白积累后抑制自身进一步产生；正自调控中，蛋白促进自身表达。回路只有一条自环，却已经把产物浓度送回产生过程，是最小的反馈结构。

负自调控可缩短系统到达新稳态的时间：启动后先以较高速率产生调控蛋白，随着蛋白积累，反馈逐渐压低表达。Rosenfeld、Elowitz 和 Alon 在大肠杆菌合成转录单元中比较无自调控与负自调控回路，观察到后者的上升时间显著缩短，并与定量模型一致。⁶ 在反馈增益、分子数和噪声来源适当时，负自调控还可压缩个体间差异；强反馈也会改变平均表达与响应增益，并可能把某些高频噪声传入输出，其效果需结合具体时间尺度。

正自调控使少量产物进一步促进自身产生，可放大弱输入并延长开启过程。若启动子响应具有足够协同性和非线性，低表达与高表达状态可以同时稳定，产生阈值、滞后和记忆；若非线性不足，则可能只是较慢、较强的单稳态响应。拓扑提示可能存在自维持，真正的开关仍由结合曲线、合成与降解速率决定。

前馈环¶

前馈环由三个节点组成：上游调控者 \(X\) 调控中间节点 \(Y\)，\(X\) 和 \(Y\) 又共同调控输出 \(Z\)。三条边各可为激活或抑制，因而共有八种有符号类型。若 \(X\) 对 \(Z\) 的直接作用与经 \(Y\) 的间接作用最终同号，称相干前馈环；二者最终异号，称非相干前馈环。两条路径把同一输入以不同长度送到输出端，由此产生时间比较。⁷

经典的相干 1 型前馈环中，\(X\) 激活 \(Y\) 和 \(Z\)，\(Y\) 也激活 \(Z\)。若 \(Z\) 需要 \(X\)、\(Y\) 同时存在，即采用近似 AND 逻辑，输入开启后必须等待 \(Y\) 积累，输出才启动；输入短暂出现时，\(Y\) 尚未越过阈值，\(Z\) 不响应。输入关闭后 \(X\) 立即消失，输出可迅速终止。这构成对开启方向的延迟和持续性检测。若启动子采用 OR 逻辑，\(X\) 或 \(Y\) 任一存在便可维持输出，延迟更可能出现在关闭方向。于是，“相干前馈环造成延迟”还必须说明输入方向和启动子逻辑。

非相干 1 型前馈环中，\(X\) 同时激活 \(Z\) 和抑制者 \(Y\)，\(Y\) 再抑制 \(Z\)。输入突然开启时，直接支路先使 \(Z\) 上升，较慢的间接抑制随后把输出拉回，形成脉冲或适应；在另一些参数区间，它能加快响应并扩大对相对变化的敏感性。前馈环的八种符号组合还可对输入增加与减少表现出不同方向的加速或延迟。Mangan 和 Alon 的分析及后续实验把这些动力学与大肠杆菌、酵母的具体转录回路联系起来，显示模体功能来自拓扑、逻辑门和时间尺度的共同作用。⁷⁸

多输出前馈环把同一 \(X\)—\(Y\) 组合接到一组靶基因。各启动子的阈值、\(X\) 与 \(Y\) 的组合逻辑和产物稳定性不同，可使靶基因按顺序开启或关闭。大肠杆菌鞭毛基因网络的定量研究测量了调控系数和启动子输入函数，并通过改变调控者水平或结合位点重新安排时间表达程序，说明一张有数字和逻辑的回路图能够提出可检验预测。⁹

单输入模块与组合调控¶

单输入模块由一个主要调控者控制一组功能相关的靶基因，调控者常同时调控自身。若各靶启动子具有不同结合阈值，调控者浓度逐渐上升时，靶基因可按固定顺序开启；浓度下降时，关闭顺序可能相反。共同调控使同一过程的多个组分协调表达，阈值差异又让组装、运输或代谢步骤获得时间次序。单输入模块因此比“一个因子调许多基因”的星状图更具体：靶基因共享调控输入，但各自的启动子函数保留差异。

稠密重叠调控子由多个调控者以不同组合连接一组输出基因。每个启动子相当于一个组合逻辑接口，可响应环境信号、代谢状态和细胞命运因子的不同搭配。细菌中多个全局与局部转录因子可共同决定碳源利用或应激程序；真核增强子则常通过多个因子、共调节复合物和染色质环境实现更丰富的组合调控。网络图显示“哪些输入可能汇合”，具体输出仍需由结合位点、协同性和染色质状态确定。

级联把调控者按层次依次连接，适合分阶段推进缓慢过程；双节点或多节点正反馈可保持细胞状态，负反馈环可稳定输出或在延迟下产生振荡。它们与前馈环经常嵌套：一个上游主调控因子可通过前馈环启动多个模块，各模块内部再由反馈确定稳态。真实网络中的局部模体彼此嵌套，一条边常同时属于多个模体。

模块、层级与跨模块连接¶

功能模块是一组共同完成相对独立任务的分子，例如核糖体装配、脱氧核糖核酸（deoxyribonucleic acid，DNA）复制起始、糖酵解或趋化感知。Hartwell 等提出的模块化细胞观强调，理解模块需要把整体行为的实验描述与组成分子的机制研究结合起来。¹⁰ 图论中的社区则是内部边相对密集、与外部连接相对稀疏的节点集合。功能模块可能在网络上形成社区，但两者并不总是一一对应：线性代谢通路内部未必稠密，一个蛋白也可同时属于复合物、信号支路和应激过程。

模块边界通常具有通透性。共享辅因子、能量与还原力把许多代谢模块耦合起来，支架和磷酸酶可连接不同信号支路，全局转录因子又跨越多个局部调控子。少量跨模块节点能够协调资源和状态，也可能成为传播扰动的瓶颈。评价模块时需综合内部连接、边的语义、亚细胞位置和共同输入输出，社区检测算法的分区只是其中一项证据。

层级组织把较小模块嵌入较大的功能单元。转录网络可由环境感知因子进入全局调控因子、局部调控子和效应基因；代谢网络可由反应簇组成通路，再由通路汇入核心代谢与生物合成系统。早期对 43 个物种代谢网络的比较报告了小而高聚类的拓扑模块嵌套成较大单元，并在大肠杆菌中与部分已知代谢功能对应。¹¹ 这种结果依赖网络投影和当时的数据库覆盖，但“多尺度嵌套、边界可重叠”仍是理解生物层级的重要框架。

生物网络的层级同时包含分支和闭环。代谢物反馈抑制上游酶，末端细胞状态改变转录因子活性，细胞间信号又把一个细胞的输出接回另一个细胞的输入。上层输入可以协调多个模块，下层反馈则改变上层决策；这些反馈关系是理解系统稳态和状态转换的关键结构。

网络结构的形成与演化¶

基因复制可同时复制蛋白结构、表达模式和部分相互作用，随后突变使两个拷贝丢失、保留或获得不同边；启动子与增强子中的结合位点周转会重接转录网络，结构域重组改变蛋白互作，水平基因转移则可一次引入带内部连接的代谢或调控模块。保守的生化约束、细胞资源与发育过程限制可行连接，自然选择又可保留在特定环境中改善响应、稳态或繁殖的网络变体。

模体的重复出现可由多种演化过程产生。度数约束、基因复制方式、突变靶大小、重组和遗传漂变都能改变某些三节点图样的出现概率；网络取样也会放大研究充分的节点。群体遗传学分析显示，部分调控网络的定性特征可由非适应过程形成。¹² 较强的功能论证需要重建边的演化变化，并测量改变这条边对回路动力学和适合度的影响。

同一模体在不同物种中出现，也可能属于趋同、共同祖先保留或统计上的独立富集。比较时应先确保节点粒度、边语义和网络覆盖相当，再追踪直系同源基因、调控位点和生化机制。拓扑保守可提示共同设计问题，参数与表达环境的改变仍能让动态功能发生分化。

拓扑候选的实验检验¶

一张网络图要转化为系统生物学解释，至少经历四层证据。第一层明确节点和边：结合实验、催化反应、遗传效应或调控扰动各自支持什么命题。第二层测量状态与时间：节点的浓度、活化形式、区室和变化顺序是否符合所画方向。第三层建立定量模型，把阈值、协同性、合成、降解和传输延迟加入局部回路。第四层实施干预，删除一条边、改变结合位点、替换启动子逻辑或控制输入脉冲，再比较输出是否按模型改变。

边的扰动往往比整节点删除更有辨别力。删除转录因子会同时破坏它的全部靶边，并可能引起广泛补偿；只突变某一启动子的结合位点，可以检验特定前馈支路。蛋白互作研究可用界面突变或结构域替换分离直接结合与共同复合物效应；代谢网络可通过同位素示踪、酶活扰动和通量测量检查静态反应图中哪些路径真正承载物质流。救援实验和正交测量帮助确认表型来自预期边，而非表达负担或非特异损伤。

网络结构还应在不同条件下复验。一个前馈环可能只在某种碳源、细胞阶段或输入持续时间下发挥延迟作用；模块边界可在应激时重排；遗传互作的符号也会随背景和表型尺度改变。可靠的结论会同时写明网络边界、细胞条件、时间尺度和读出，并把没有测到的边保留为未知，而非自动视为不存在。

从这里继续学习时，可以沿两条路线展开：一条在反馈、开关与振荡中把回路写成速率方程和相图，比较稳定性、双稳态、滞后、振荡与鲁棒性；另一条在系统建模、参数推断与实验检验中讨论模型边界、参数可辨识性、组学约束和扰动预测。网络拓扑提出可检验的局部结构，动力学与实验决定这些结构在真实系统中做了什么。

参考资料与延伸阅读¶