[论文解读] Robust gene regulation: Deterministic dynamics from asynchronous networks with delay

Konstantin Klemm, Stefan Bornholdt|ArXiv.org|Sep 25, 2003

Gene Regulatory Network Analysis参考文献 1被引用 10

一句话总结

该论文表明，异步布尔网络中的传输延迟和低通滤波可诱导出稳健的、确定性的动力学行为，其表现类似于同步更新，即使没有中心时钟。关键发现是：足够的延迟可稳定周期性吸引子，从而在异步更新下实现可重复的准确定性行为，即使在基因调控网络中亦如此。

ABSTRACT

We compare asynchronous vs. synchronous update of discrete dynamical networks and find that a simple time delay in the nodes may induce a reproducible deterministic dynamics even in the case of asynchronous update in random order. In particular we observe that the dynamics under synchronous parallel update can be reproduced accurately under random asynchronous serial update for a large class of networks. This mechanism points at a possible general principle of how computation in gene regulation networks can be kept in a quasi-deterministic ``clockwork mode'' in spite of the absence of a central clock. A delay similar to the one occurring in gene regulation causes synchronization in the model. Stability under asynchronous dynamics disfavors topologies containing loops, comparing well with the observed strong suppression of loops in biological regulatory networks.

研究动机与目标

• 理解在缺乏中心时钟的情况下，异步基因调控网络中如何涌现出稳健的、确定性的动力学行为。
• 研究在何种条件下，异步更新可再现离散动力学网络中同步更新的动力学行为。
• 探索为实现确定性行为而对异步系统施加的拓扑约束，特别是环路抑制。
• 使用离散时间、延迟、滤波的布尔网络对基因调控进行建模，采用阶跃函数激活和状态依赖更新。
• 确定传输延迟是否能稳定周期性吸引子并抑制异步更新序列中的随机性。

提出的方法

• 通过显式时间延迟 τ 和小时间步长 ε 的差分方程，形式化基因调控的离散时间动力学模型。
• 通过阈值浓度变量定义节点状态：x_i(t) = Θ(c_i(t) - 1/2)，其中 c_i 基于延迟输入信号演化。
• 通过在每个时间步随机选择节点并根据其调节子的延迟状态更新其 c_i，实现异步更新。
• 使用类似 sigmoid 的阶跃函数 f(s) ≈ Θ(s)，基于输入信号 s_i(t) = Σ_j w_ij x_j(t) - a_i 近似基因激活为二值。
• 通过 ε = αΔt 引入滤波以平滑过渡，通过 τ = ϑ/Δt 延迟以建模转录和翻译延迟。
• 比较随机异步更新与同步更新的动力学，通过时间序列间误差 E 衡量差异。

实验结果

研究问题

• RQ1 在随机顺序下，异步更新能否再现基因调控网络中同步更新的确定性动力学？
• RQ2 传输延迟在稳定异步更新网络中周期性吸引子方面起到什么作用？
• RQ3 过滤（ε）与延迟（τ）如何共同影响网络动力学的可重复性与周期性？
• RQ4 从异步动力学的角度看，为何环状基序（尤其是反馈环）在生物调控网络中受到抑制？
• RQ5 在何种拓扑条件下，异步动力学可保持同步与确定性，何时会失效？

主要发现

• 当传输延迟足够大（τ = 50,000）且滤波足够强（ε = 1/500）时，16节点网络在异步更新下遵循稳定的极限环，周期为28个宏观时间步，与同步更新的吸引子完全一致。
• 异步与同步动力学之间的误差 E 随延迟 τ 增大呈指数衰减，表明向确定性行为收敛。
• 当无延迟（τ = 0）时，即使有滤波，异步更新仍导致非周期性、不可重复的动力学，翻转序列不规则。
• 在无延迟情况下，具有相互激活的网络（如 N=2 个节点）表现出去同步：开启与关闭转换形成解耦链，最终坍缩为固定点。
• 具有可分离翻转序列的拓扑结构（如孤立环或竞争振荡器）无法维持同步，解释了此类基序在生物学中被抑制的原因。
• 延迟的存在使信号传播能够同步网络中各节点的更新，从而有效形成一种分布式、延迟诱导的钟控模式。