5.集群
我们和其他小组最近的一些研究，强调了在与情况的因果互动和内部信息处置惩罚中适应性和鲁棒性之间的平衡（Walker et al. 2016），已知信息处置惩罚特性体现出从稳定阶段（以信息存储为主）到混沌阶段（以信息通报为主）之间临界转变距离的强烈依赖性。众所周知，临界系统具有生物系统的两个基本特征：对种种差异情况条件的适应性和鲁棒性。我们期望基因调控中的临界性普遍存在，这将为生物网络奇特的信息处置惩罚和控制特性研究提供参考。
Control kernel 控制核
庞大系统研究的主旨之一是理解并理预测出系统涌现的全局动态，例如动态轨迹或最终状(Barabási 2016）。这些代表系统集群行为的全局状态可以通过施加外部因果干预来改变。引导整个系统从初始状态到达期望的最终状态。这为许多领域应用都提供了强大的工具，包罗从舆论动态相关的问题到流行病，从细胞的分化到社会性昆虫的巢穴选择问题等等。
最近的研究生长了庞大网络控制相关的数学基础，试图通过理解因果结构与各组成部门动态的更新规则来寻找严格的控制机制。Liu 等人（2011）生长了他们的前瞻性框架，将具有线性动态历程的有向庞大网络可控性问题简化为一个图论问题。利用他们的理论可以确定一组能够完全控制整个网络动态的驱动节点（driver nodes）。他们还发现，所确定的驱动节点往往是低度节点，而且总数主要由网络的度漫衍决定。
生物系统模型大多是非线性动力学系统，完全可控性既不行行也没有须要。相反，找到一种特定的控制机制，能够将一个给定的系统从特定的初始状态或轨迹引导到另一组期望状态或轨迹荟萃，是更为现实的，往往也足够的（Cornelius et al. 2013）。例如，控制基因调控网络的主要目的之一是将系统引导到特定的细胞状态，而非任意的表达水平。这种表型控制通常被定位在整个网络的一小部门中，以能驱动系统到期望吸引子的控制核来界说。不外，控制核的巨细通常比一般随机网络还要更大（Liu et al. 2011; Kim et al. 2013; Gates &Rocha 2016; Choo et al. 2018）。这讲明生物系统的可控性水平可能存在一定的规模，可以通过控制核的巨细来量化。因此我们会将可控性水平视为生物系统信息架构的基本特征（Walker et al. 2016; Kim et al. 2015），以区别于所有非生物集群。
Causal emergence 因果涌现
识别庞大系统中的因果关系是理解其涌现行为最基本的要领之一。然而这却是一项困难的任务：不仅因为庞大系统中的因果结构很庞大（坦率地讲它们确实很庞大），而且因为因果结构可以在多个时空尺度上进行分析，往往让我们搞不清楚哪种尺度是分析的首选（如果有的话）。恒久以来，人们一直认为主因结构只能在最低（微观）尺度上界说，更高尺度上的因果结构是有用的，但只是对主因结构的简陋描述。即认为微观层面的因果关系确定了所有更高条理的因果结构，更高（宏观）条理的因果并没有实际上的因果孝敬（Kim 1993，2000）。
这种看法下的还原论要领在科学各领域取得了广泛的乐成，包罗物理学和化学的重大生长。尽管如此，因果涌现的另一种可能——宏观尺度是因果结构的重要驱动，被一再提出来，因为它能对种种类型的生物集群行为提供更简朴的解释：从表观遗传学到蚁群的群体决策，再到大脑中的心理状态等等。以往对这种因果涌现的研究大多局限于定性论证，而最近 Hoel 等人（2013）和 Hoel（2017）所做的定量研究讲明，在一些系统中，因果结构可能无法完全被微观尺度结构捕捉，相反宏观层面的因果结构会更有效、具有更富厚的信息。
如上所述，我们认为寻找一种庞大系统的控制机制，可以成为研究差异条理因果结构的有用工具，它能从局部（微观）条理因果深入洞悉驱动节点对全局（宏观）系统动态转向的影响。
生物集群系统信息架构与因果结构之间的耦合关系
如上节所讨论，在以往研究中已经体现了生命系统信息架构与因果结构之间的关系。然而它们之间是否存在有意义的关系、若存在是否能够区分差异生物集群行为，这些问题却仍然是模糊、没有定论的。在这里我们提供了一个定量的例子，基于经验数据集建设了一个生物网络的动力模型，其规模较小故允许进行统计分析。由于系统生物学中高通量技术和数据驱动要领的泛起，细胞生物通路的布尔网络模型是满足这一须要条件的最佳候选之一。
基于最近 Daniels 等人（2018）提供的迄今为止最全面的临界性视察，对从细胞集群数据库（Helikar et al. 2012）中获得的67个细胞调控布尔网络进行了研究，显示所有网络都接近临界状态。另一方面，Kim 等（2013）分析了种种生物分子的调控布尔网络，并确定了控制核——控制细胞网络到达预期状态所必须的最小节点集。研究讲明，对于大多数网络而言，控制核只是网络组件的一小部门。因此，利用同类生物分子的调控布尔网络来研究生物系统、并将信息架构与其因果机制结合起来，是我们自然的选择。
在这里，我们以粟酒裂殖酵母（S.pombe，Fission Yeast）的细胞周期历程为案例，展示了如何在生物系统中研究信息架构与因果机制之间的关系（Walker al. 2016；Kim al. 2015），并作为生物中心功效体现的一个简朴实例。
裂殖酵母细胞周期的布尔模型描述
以布尔网络对裂殖酵母的细胞周期进行建模（图5），共有9个要害卵白组成（Davidich & Bornholdt，2008），通过被称为 G1-S-G2-M 四个步骤的细胞功效进行控制：从细胞生长到DNA复制，再到细胞破裂成两个子细胞。在网络中，每个代表卵白质的节点都有一个0或1的布尔值，划分体现其缺失或存在，每个节点状态的组合集被称为网络状态或网络配置。单个节点的状态，是与其它节点生化作用状态的函数，类型包罗交互、抑制和激活。
将局部动力规则应用于每个节点，以模拟系统一个时间步长的演化。在对应细胞周期起点的特定网络配置上重复进行时间演化，从而以准确的顺序再现生物历程各个阶段序列。从512种可能的网络状态中的每一种开始生成网络状态空间，即所有可能时间演化的动力学轨迹，其景观就代表了裂变酵母细胞周期布尔网络的全局动态。网络状态空间包罗13个吸引子、一个或一组多个网络状态，其中任意时间演化轨迹都市收敛。由于动力特性是确定性的，凭据每个网络状态收敛到哪个吸引子，整个网络状态空间可以划分为13个不相交的子群。其中，约74%网络状态结束所在的主吸引子是唯一具有生物功效的吸引子（Davidich & Bornholdt 2008）。
Kim 等人（2013）最近关于控制这些生物调控布尔网络的时间演化轨迹的研究讲明，存在一个控制核——能通过改变所有吸引子状态景观收敛到主吸引子，从而决定整个系统命运节点的最小子集。对于裂殖酵母细胞周期网络，研究讲明整个网络中有4个卵白是控制核——即当 Rum1、 we1、 Ste9 和 Cdc25 的状态牢固为1、1、1、0 ，与它们所在生物功效主吸引子中的状态相同时，从每种可能网络配置开始的所有时间进化轨迹都市趋于生物功效吸引子。
信息通报与控制核的关系——以裂变酵母细胞周期网络为例
在最近的事情中（Kim et al. 2015），我们通过盘算网络中每一对节点之间的通报熵来量化裂变细胞周期网络的信息架构，从结果中识别出生物学的特有的模式，并通过控制核证明它与控制机制之间的密切关系。为了区分生物系统表征的信息模式，我们将通报熵的分析结果与对随机网络集成进行的相同分析进行了比力，这些网络被构建为与每个酵母网络共享某些相同拓扑特征，但不执行相同的生物功效。
为了进行比力，我们利用了 Erds-Rényi （ER）和无标度（SF）两种类型的随机网络集成。ER随机网络与裂变酵母网络，除了网络巨细外没有配合的结构偏向。相比之下，SF随机网络集成与细胞周期网络一样，每个节点保持相同数量的激活和抑制连接，这使得SF与生物网络有相同的拓扑度漫衍。我们将单个网络处置惩罚的总信息界说为其中所有节点对之间通报熵的总和。结果显示，酵母细胞周期网络处置惩罚的信息量比任何荟萃中大多数随机网络都多（图6a）。这是一个直接丈量结果，支持了与信息处置惩罚和临界点相关的假设，即在生物系统中临界点四周存在着信息传输的优化，不像以前没有基于经验数据集的定性或一般模型的研究，在这里这点被定量地清楚展示了出来。
通过凭据每对节点通报熵的丈量值分析缩放标度，我们还研究了处置惩罚后的信息在各个网络中的漫衍。标度关系讲明，与随机网络相比，生物网络具有明显异常的值（图6b），生物治理信息最显著的特点是，生物网络中信息被处置惩罚得越发均匀。标度模式还体现出了统计学上奇特的高-中品级节点对规模，与随机网络中相应尺度的节点相比，生物网络的通报熵明显更大。而且我们发现，这种生物学上奇特的规模是由控制核节点和非控制核节点之间的信息处置惩罚主导的（图6c）。
细胞周期网络在信息结构上体现出与随机网络相区此外特征模式，控制核节点通过控制有关结构的全局因果（状态转换图中的因果流），在网络信息通报中起主导作用，调治着系统全局状态转换。
图5：裂殖酵母细胞周期调治的布尔网络模型。节点代表调控卵白，边缘体现节点之间的两种生化作用：激活（以箭头结尾）和抑制（以条形结尾）。其中红色节点是控制核，当牢固在特定值时，它可以调治网络的全局行为。该图凭据 Davidich 和 Bornholdt（2008）和 Kim 等人（2013）的数据重新生成。
图6：裂变酵母细胞周期网络中的信息架构和因果结构。a ER（绿色）和 SF（蓝色）网络集成处置惩罚总信息的漫衍。每个数据点体现在x轴上，传输的给定信息总量的各个荟萃中的单个网络的数量。b 裂变酵母细胞周期中成对节点间通报熵的标度，ER和SF网络划分用红色、绿色和蓝色体现。对于1000个网络的样本，盘算每个随机网络荟萃的平均值和尺度差。虚线之间的区域体现裂变酵母的信息标度在生物学上奇特的规模。c 裂变酵母细胞周期网络信息通报的标度关系。CK 体现控制核节点，而 NCK 则是一组不包罗在控制核内的节点。所示数据与中间面板中裂变酵母通报熵的标度模式相同。凭据一对中每个节点是在CK照旧在NCK中，标度模式被分为四类信息传输。这些数字来自 Kim 等（2015）（在线彩色图）