基因调控网络的定性建模：从逻辑调控图到标准 Petri 网
Chaouiya, Remy, Ruet, Thieffry. ICATPN 2004, LNCS 3099, pp. 137–156. Springer.
核心贡献
提出一种系统性重写方法，把布尔逻辑基因调控图（Logical Regulatory Graph，R. Thomas 广义逻辑方法）机械转换为标准 Petri 网（称为布尔调控 Petri 网，BRPN），从而把逻辑方法的”反馈回路-动力学性质”直觉分析，与 Petri 网框架的代数/可达性分析工具结合起来。
背景与动机
- 调控网络（基因调控、代谢通路、信号传导）是全部生物功能的核心，其复杂性超出直觉，需要数学建模。
- 定量方法（ODE/偏微分/随机方程）需要精确参数，而生物系统普遍缺乏定量数据 → 定性方法更可行。
- 逻辑方法（布尔或多级离散）已有工具 GIN-sim，但状态数随基因数指数增长，需要解析方法处理大网络。
- 已有的 Petri 网基因网络应用多依赖仿真、且使用定制表示法，缺乏系统化的标准 PN 转换方法。本文填补这一空缺。
逻辑调控图（布尔情形）
一个调控图 R = (G, I, K)：
- G：节点集（基因），表达水平取值 {0,1}（产物缺失/达到足够水平）。
- I(j)：基因 g_j 的输入集（所有指向它的作用源基因）。
- K：逻辑函数集。K_j(X) 给出当 X（g_j 输入中”正在起作用”的交互子集）为真时，g_j 趋向的表达水平。
- 交互分激活/抑制；动力学图（dynamical graph）的顶点是系统状态（n 元组），边是状态迁移，通常按全异步或全同步假设生成。
BRPN 构造（核心定义）
初步重写
- 每个基因 g_j → 一个库所；每个参数 K_j(X) → 一个变迁 t_{g_j,X}。
- X 中的输入库所用普通弧连接，I(j)\X 中的用抑制弧连接；K_j(X)=1 时 g_j 是输出，=0 时是输入。
- 两个约束：① 调控产物不被消耗 → 用读弧（自环）；② 每库所托肯数限 1 → 容量限制。
补位变换（关键技巧）
- 对每个库所 g_i 定义互补库所 ¬g_i，两者托肯之和恒为 1。
- 变换后容量约束自动满足，抑制弧不再需要——得到纯标准 Petri 网（Pre/Post 取值 {0,1}）。
- 库所集 P = G ∪ ¬G（2n 个），变迁集 T = {t_{g_i,X}}。
- 注意（Remark 2）：BRPN 无法区分互补库所哪个代表”产物存在”，因此从 BRPN 不能唯一还原原调控图（Property 1 只对带抑制弧的初步网成立）。
结构性质
- Property 2（关联矩阵）：C = Post − Pre 是分块结构，每个基因对应一个 2×d_i 小块，元素取值 {-1,0,1}。由于读弧的存在，关联矩阵不反映底层图结构。
- Property 3（有效性守恒）：合法初始标识（M(g_i) = 1 − M(¬g_i)）下任何可达标识仍合法，BRPN 是 1-safe 的。
- Property 4（动力学等价）：布尔调控图的异步动力学图与 BRPN 的可达图双向等价——状态间存在迁移当且仅当对应变迁在标识下使能并可引发。
调控回路（Regulatory Circuits）分析
Thomas 法则的 PN 验证：
- 正回路（偶数个抑制）↔ 多稳态（multi-stationarity）；
- 负回路（奇数个抑制）↔ 持续周期行为。
关键结论：
- Property 5/6：n 基因孤立回路的 BRPN 有 2n 库所、2n 变迁，关联矩阵为 2n×2n 分块对角形，每块由激活（[[−1,1],[1,−1]] 型）或抑制（[[1,−1],[−1,1]] 型）块构成。
- Property 7（同构归约）：任意正回路 BRPN 同构于全正规范回路 C_n^+，任意负回路同构于 C_n^-（证明通过交换相邻负交互对应的变迁/库所做”翻转”）。→ 每种长度只需研究一个代表网，即可覆盖同类回路的全部动力学性质。
应用案例
论文给出两个真实生物学网络（从已发表的逻辑模型推导其 BRPN）：
- 果蝇（Drosophila melanogaster）早期细胞周期控制——合胞体分裂期的细胞周期转换调控。
- 拟南芥（Arabidopsis thaliana）花形态建成——ABE 模型（花器官发育的四类基因组合）。
可行动点 / 方法论启示
- 想把离散动力学系统（不限于基因网络，也适用逻辑控制电路、生化开关）纳入 Petri 网工具链时，“互补库所消抑制弧”是通用技巧：任何带抑制弧/读弧的布尔系统都能转成标准 P/T 网并用现成可达性/不变量工具分析。
- 回路符号（正/负）决定定性动力学类别这一结论，与硬件锁存器、正反馈双稳开关、负反馈振荡器（如环形振荡器）的工程直觉完全同构。
- 局限：仅布尔（二值）情形给出完整证明；多级（multi-level）推广是文中提出的后续方向。
关联
- 云盘
/田浩然上传的资料/Petri/目录收录 Application and Theory of Petri Nets 会议论文系列，本篇为该系列之一。 - 系统生物学 / 复杂网络稳定性分析可与控制理论笔记互参。