跳转至

第7章：生物振子与开关（Biological Oscillators and Switches）

§1 作者

本章为 J.D. Murray《Mathematical Biology I》(Third Edition, 2002) 第 7 章,作者是 J.D. Murray 本人。承接 Ch 6 酶动力学的"非线性反应动力学", 本章进入生物振荡器与开关——心脏起搏器、果蝇羽化、糖酵解振荡、神经动作电位、睾酮分泌、白细胞 cAMP 振荡等。Murray 强调这些振荡器的共同特征:自主动力系统 (autonomous) 而非外部驱动, 周期从秒到小时到天, 是 ODE 系统极限环或混沌的体现。Murray 在 7.6 节给出的睾酮分泌模型是其与 W. R. Smith 1980s 早期合作工作的体现, 是"延迟微分方程" 在临床应用中的经典案例。

§2 内容概述

第 7.1 节给出生物振荡器综述: 心脏起搏器 (周期 ~1 秒), 神经动作电位 (Hodgkin-Huxley 1952, 获诺贝尔奖), 呼吸周期 (~5 秒), 糖酵解振荡 (Selkov 1968, 周期几分钟), 睾酮水平振荡 (2-3 小时), 黏菌 (Dictyostelium) cAMP 振荡 (几分钟), 细胞分裂周期 (小时)。 Murray 强调这些自主振荡器 (autonomous) 必须热力学开放系统——通过持续物质/能量输入维持周期行为。Lotka 1910 首次提出理论振荡反应, Bray 1921 首次实验观测, 但被怀疑违反热力学第二定律 (实则因系统开放, 振荡最终衰减)。

第 7.2 节给出反馈控制机制——Goodwin 1965 经典模型 (7.3):

\[ dM/dt = V/(D + P^m) - a M, \quad dE/dt = bM - cE, \quad dP/dt = dE - eP, \]

其中 \(M\) 是 mRNA, \(E\) 是酶, \(P\) 是产物; 产物 \(P\) 抑制 mRNA 转录 (负反馈)。 Griffith 1968 证明: Hill 系数 \(m > 8\) 时才有振荡——这在生物上不现实。 改进: 改 \(P\)-方程为 Michaelis-Menten 降解, 即使 \(m = 2\) 也可振荡。 推广到 \(n\) 物种负反馈 (7.4) \(du_r/dt = u_{r-1} - k_r u_r\), Tyson-Othmer 1978 证明 \(m > m_0(n) = \sec^n(\pi/n)\) 时失稳: \(n = 3\) 时 \(m_0 = 8\); \(n = 4\) 时 \(m_0 = 4\); \(n \to \infty\) 时 \(m_0 \to 1\)。 周期估计 (Hunding 1974, Rapp 1976): \(T \approx 2N/K\) (其中 \(K\) 是最小速率常数, \(N\) 是具有相近时间尺度的物种数)。

第 7.3 节给出2 物种反应动力学的定性结果——null cline 几何 + Poincaré-Bendixson 给出极限环存在条件。 关键观察: 当 \(\text{tr}\mathbf{A} > 0\) 且 \(|\mathbf{A}| > 0\) 时, 稳态是不稳定节点或螺旋, 配合 confined set 保证极限环。 稳态 sign 模式 (7.17) 4 种, 只有 (b) 和 (c) 可产生振荡。 "periodic bursting" (7.6 节) 描述多稳态区+慢变参数的复杂行为——与 Ch 1 Budworm 滞后突变 + Ch 6 自杀底物多重稳态同源。

第 7.4 节给出Schnackenberg 1979 最简振荡器 (7.19)-(7.20):

\[ du/dt = a - u + u^2 v, \quad dv/dt = b - u^2 v, \]

由 trimolecular reaction \(2X + Y \to 3X\) 推导。 稳态 \(u_0 = a + b\), \(v_0 = b/(a + b)^2\)。 参数域 \(b - a > (a + b)^3\) 决定不稳定 (7.23)。 参数化方法 (Murray 1982, 7.27): 边界曲线 \(a = u_0(1 - u_0^2)/2, b = u_0(1 + u_0^2)/2\)。

\(\lambda\)-\(\omega\) 系统 (7.30) 给出精确极限环解——但难以从化学反应推导, 主要是数学原型。 在极坐标下 \(dr/dt = r\lambda(r), d\theta/dt = \omega(r)\), 极限环 \(r = r_0\)。

第 7.5 节给出FitzHugh-Nagumo 模型 (7.39) 作为 Hodgkin-Huxley 1952 (获诺贝尔奖) 的简化:

\[ dv/dt = v(a - v)(v - 1) - w + I_a, \quad dw/dt = bv - \gamma w, \quad 0 < a < 1, \tag{7.39} \]

\(v\) 模拟膜电位, \(w\) 替代 Hodgkin-Huxley 三个门变量 \(m, n, h\)。 \(I_a\) 是注入电流。 \(I_a = 0\) 时系统可兴奋 (excitability, 类似 Ch 6 阈值现象); \(I_1 < I_a < I_2\) 窗口内出现 Hopf 分岔→极限环→周期性神经发放 (periodic firing)。 几何条件 (7.42): 极限环存在要求 \(v\) null cline 斜率 \(\tan\theta < b/\gamma\)。 该模型仍是当代神经动力学的核心——22 年后 (2026) 仍是计算神经科学的入门模型。

第 7.6 节给出睾酮分泌的延迟模型——Smith-Murray 1980s 早期工作 (7.44):

\[ dR/dt = f(T) - b_1 R, \quad dL/dt = g_1 R - b_2 L, \quad dT/dt = g_2 L(t - \tau) - b_3 T, \]

其中 \(R\) 是 LHRH, \(L\) 是 LH, \(T\) 是睾酮, \(\tau\) 是血液循环时间。 特征方程:

\[ \lambda^3 + a\lambda^2 + b\lambda + c + d e^{-\lambda \tau} = 0, \quad a = b_1 + b_2 + b_3, \quad d = -f'(T_0) g_1 g_2 > 0. \tag{7.49} \]

Murray 用 \(\lambda \to w\) 平面映射证明存在临界 \(\tau_c\), 当 \(\tau > \tau_c\) 时稳态失稳→睾酮周期振荡, 周期 \(2\pi/\nu\)。 临床应用: 化学去势 (chemical castration) 治疗前列腺癌与性犯罪者——欧洲使用后, 儿童性侵者再犯率从 50% (美国, 不用药) 降至 2% (欧洲, 用 Depo-provera)。

§3 核心方程与概念

3.1 通用振荡器 \(n\) 物种负反馈

(7.4) \(du_r/dt = u_{r-1} - k_r u_r\) (慢变量) + \(du_1/dt = f(u_n) - k_1 u_1\), 反馈函数 \(f(u) = V/(D + u^m)\) (Hill)。 振荡存在条件 (Tyson-Othmer 1978): \(m > m_0(n) = \sec^n(\pi/n)\)。 周期 \(T \approx 2N/K\)。

3.2 Poincaré-Bendixson 与 2 维极限环

2 维 ODE (7.10) \(du/dt = f, dv/dt = g\)。 稳态 sign 模式 (7.17) 4 种, 极限环存在当 (b) 或 (c) + \(\text{tr}\mathbf{A} > 0\) + confined set。

3.3 Schnackenberg 振荡器 (7.20)

\[ du/dt = a - u + u^2 v, \quad dv/dt = b - u^2 v. \]

稳态 \((u_0, v_0) = (a + b, b/(a + b)^2)\)。 \(\text{tr}\mathbf{A} = b - a - (a + b)^3\), \(|\mathbf{A}| = (a + b)^2\)。 失稳条件 \(b - a > (a + b)^3\) (7.23)。 参数化边界 (7.27):

\[ a = u_0(1 - u_0^2)/2, \quad b = u_0(1 + u_0^2)/2, \quad u_0 > 0. \]

3.4 FitzHugh-Nagumo 神经模型 (7.39)

\[ dv/dt = v(a - v)(v - 1) - w + I_a, \quad dw/dt = bv - \gamma w, \quad 0 < a < 1. \]

\(v\)-null cline 是 \(w = v(a - v)(v - 1) + I_a\) (三次), \(w\)-null cline 是 \(w = bv/\gamma\) (直线)。 极限环存在当 \(I_1 < I_a < I_2\) (窗口), 几何条件 (7.42) \(\tan\theta < b/\gamma\)。 Hopf 分岔: \(I_a = I_1\) 时, 稳态从稳定焦点变不稳定螺旋, 产生小振幅极限环, 周期 \(T = 2\pi/|\text{Im}\lambda_c|\)。

3.5 睾酮分泌延迟模型 (7.44)

\[ dR/dt = f(T) - b_1 R, \quad dL/dt = g_1 R - b_2 L, \quad dT/dt = g_2 L(t - \tau) - b_3 T. \]

特征方程 \(\lambda^3 + a\lambda^2 + b\lambda + c + d e^{-\lambda \tau} = 0\), 其中 \(a = b_1 + b_2 + b_3, b = b_1 b_2 + b_2 b_3 + b_3 b_1, c = b_1 b_2 b_3, d = -f'(T_0) g_1 g_2\)。 \(\tau = 0\) 时 Routh-Hurwitz 条件保证稳定; 临界 \(\tau_c\) 由 \(\text{Re}\lambda = 0\) 给出, 可用 \(\lambda \to w\) 平面映射估计。

3.6 \(\lambda\)-\(\omega\) 系统 (7.30)

\[ du/dt = \lambda(r) u - \omega(r) v, \quad dv/dt = \omega(r) u + \lambda(r) v, \quad r = (u^2 + v^2)^{1/2}. \]

极坐标 \(dr/dt = r \lambda(r), d\theta/dt = \omega(r)\)。 极限环 \(r = r_0\) 当 \(\lambda(r) > 0\) (\(r < r_0\)) 且 \(\lambda(r) < 0\) (\(r > r_0\))。

§4 关键结论

  1. Goodwin-Hill 模型显示 Hill 系数 \(m > 8\) 才能振荡——这在生物上不现实, 提示经典负反馈不是细胞周期振荡的主要机制。 现代分子生物学发现两个负反馈 (例如双负调节, Tsai 2008) 或正-负反馈组合 (双稳态+振荡, Tyson-Novák 2010) 是真实细胞周期振荡的核心。

  2. Tyson-Othmer 1978 的 \(m_0(n) = \sec^n(\pi/n)\) 是负反馈振荡的"普适阈值"——增加反馈环长度可大幅降低 \(m_0\)。 这一观察 22 年后 (2026) 仍是合成生物学"工程化振荡器" 的设计原则 (Elowitz 2000 repressilator 用了 \(n = 3\) + 高 Hill 系数, 后续设计改用更长环路 + 较低 Hill)。

  3. Schnackenberg 1979 振荡器的参数域 \(b - a > (a + b)^3\) 是最简单的非线性极限环存在条件——可在课堂上 5 分钟讲完, 但内涵涵盖了所有 2 维 ODE 振荡器理论。 Murray 1982 的"参数化" 方法 (用 \(u_0\) 作参数) 避免解三次方程——这是 22 年后 (2026) 仍在使用的教学技巧。

  4. FitzHugh-Nagumo (7.39) 捕获 Hodgkin-Huxley 4 维模型的关键现象——\(I_a = 0\) 时的可兴奋性 + \(I_1 < I_a < I_2\) 的极限环发放。 22 年后 (2026) 该模型仍是计算神经科学的"教科书"——但实际的神经元有 10+ 离子通道, 简化模型只保留钠/钾的定性特征。

  5. 睾酮分泌延迟模型预测 \(\tau > \tau_c\) 时周期振荡——这是临床上化学去势的科学基础。 LHRH 类似物 (Goserelin) 通过负反馈压低 LH 与睾酮, 周期振荡消失。 Murray 1980s 早期工作的临床影响在 22 年后(2026) 仍深远, 是"数学→临床" 成功案例。

  6. 2 物种 ODE 的极限环存在性需要 trimolecular 反应 (Hanusse 1972)——单/双分子反应在 2 维 ODE 框架内不能产生极限环。 这是 Murray 7.4 节强调的关键: "为了得到振荡, 你需要至少三阶反应"——这与 Ch 6 自催化 \(x^2 y\) 同源。

  7. periodic bursting (7.7 节) 是"快-慢"系统的两尺度耦合——快子系统是振荡器, 慢子系统是调制参数, 产生"一连串脉冲" 行为。 这在神经元的 burst firing、胰脏 \(\beta\) 细胞的胰岛素分泌、心脏的早搏中都观察到, 22 年后(2026) 仍是计算神经科学的活跃课题。

  8. Lotka 1910 与 Bray 1921 是"化学振荡" 概念的开端——Murray 7.1 节强调它们的跨时代重要性: 50 年的"违反热力学" 怀疑 → Belousov 1951 重新发现 → Belousov-Zhabotinskii 反应成为 Ch 8 的核心。

§5 挑战和开放性问题

  1. 细胞周期振荡的完整机制: Goodwin-Hill 模型是简化, 真实细胞周期涉及 cyclin-CDK 的多反馈、磷酸化-去磷酸化级联、检查点 (checkpoints)。 Tyson-Novák 2001/2010 的"翻转开关 + 振荡器" 模型远超 Murray 1977 的 Goodwin-Hill。 22 年后 (2026) "数学细胞周期" 已成熟为独立子领域。

  2. 神经发放的离子通道细节: FitzHugh-Nagumo 保留 Hodgkin-Huxley 4 维模型的定性特征, 但钙通道、钾通道亚型 (Kv1, Kv2 等) 数十种 在真实神经元中起关键作用。 22 年后 (2026) 神经动力学的多尺度模型 (Hodgkin-Huxley 类型 + 突触 + 网络) 仍是热点。

  3. 周期倍化与混沌: Murray 在 7.3 节末提到 Decroly-Goldbeter 1987 在 Dictyostelium cAMP 模型中观察到 period doubling→chaos, 但没深入。 22 年后(2026) 这是细胞信号转导的成熟课题 (Goldbeter 1996, 2002)。

  4. 生物开关的多稳态: Murray 7.3 节用 null cline 几何给出开关的定性条件, 但真实的生物开关 (如 \(\lambda\) 噬菌体溶源/裂解、细胞命运决定) 涉及基因调控网络, 通常需要 2+ 维相空间。 Gardner-Cantor 2000 等的"synthetic toggle switch" 是工程化实例。

  5. 可兴奋性与 Hopf 分岔的关系: Murray 7.5 节强调 FitzHugh-Nagumo 的可兴奋性是"线性稳定, 大扰动大响应", 7.5 节后半段提示"在 \(I_1 < I_a < I_2\) 窗口内是 Hopf 失稳→极限环"。 22 年后(2026) 这仍是神经动力学的标准框架, 但可兴奋性 vs 振荡的过渡细节仍在研究 (Klaus-Dilão 2012)。

  6. 睾酮模型对个体差异的预测: Murray 1980s 模型假设所有男性的激素参数相同, 但实际个体差异 (年龄、遗传、疾病) 影响响应。 现代精准医学需要参数化的个体模型, 远超 Murray 1977 的简化。

  7. 化学去势的伦理: Murray 7.6 节提到"用 Depo-provera 治疗性犯罪者", 这是"科学-政策-伦理" 的复杂议题。 22 年后 (2026), 这仍是争议——美国部分州 (California, Florida, Texas) 立法化学去势, 但伦理问题 (自主性、报复性 vs 治疗性) 仍在持续。

  8. 周期性爆发的几何分析: Murray 7.3 节末用图 7.6-7.7 描述 periodic bursting, 但解析理论 (不变流形分析、慢-快分离) 直到 2010s 才有系统化 (Desroches 2012 等)。 Murray 2002 的定性讨论是这一方向的开端。

§6 个人反思与批判性分析

Murray 的"反馈控制机制"教学: 7.2 节的 Goodwin-Hill 模型 + Tyson-Othmer \(m_0(n)\) 是该领域最核心的解析结果——它给出"为什么真实细胞需要多个反馈环" 的数学答案。 Murray 没明确指出双负反馈是真实细胞周期 (如 CDK-cyclin) 的核心, 但 22 年后(2026) 这已是被反复验证的共识。

对 FitzHugh-Nagumo 的反思: 该模型是 22 年后(2026) 计算神经科学的"基础", 但 Murray 2002 没提它的两个关键扩展: (a) 噪声注入 (stochastic FHN) 用于"神经发放的统计性", (b) 空间耦合 (FHN + 扩散) 用于"神经传播"。 这两个方向 Murray 留给 Ch 12 (相互作用振荡器) 和 Vol II Ch 1 (反应扩散)。

对睾酮模型的时间敏感性: Murray 7.6 节的延迟模型是 1980s 早期工作, 22 年后(2026) 临床已用GnRH 类似物 (Goserelin, Leuprolide) 治疗前列腺癌 30+ 年, 但化学去势的伦理 仍争议。 Murray 提到"欧洲降低儿童性侵者再犯率到 2%" 的数据出自 1990s, 22 年后(2026) 该数据被复现/质疑。 模型的临床意义稳固, 但应用伦理仍在演化。

Murray 的"参数化方法" (7.4 节末 (7.27)) 是一个未充分发扬的技巧——把"在 \((a, b)\) 空间画稳定域" 转化为"用 \(u_0\) 作参数画曲线", 避免了"解三次方程" 的代数困难。 这种"聪明的参数化"在 22 年后(2026) 的计算物理中仍是核心技巧, 但 Murray 没强调它的普遍性。

对"反应+扩散" 的伏笔: 7.1 节末 Murray 提示"下一章和后续章节讨论波的传播", 这是 Ch 8 (BZ 反应+反应扩散) 的桥梁。 Murray 把"化学反应的空间耦合" 留到 Vol II, 这是一个合理的教学安排——避免在该章同时讨论时间+空间复杂性。

对 Goodwin 模型的批评: Goodwin 1965 提出 mRNA→酶→产物→抑制 mRNA 的 3 维负反馈, Hill 系数 \(m > 8\) 才能振荡。 这是过度简化的抽象——真实细胞中, mRNA 转录调控远比 Hill 方程复杂 (转录因子结合、染色质修饰、非编码 RNA 等)。 Murray 没明确指出 Goodwin 模型的"经验机制薄弱"——只把它作为"教学原型"。 22 年后(2026) Goodwin 模型在合成生物学中被工程化使用 (Elowitz 2000 repressilator), 但真实细胞使用更复杂的调控网络。

关于"period doubling → chaos" 的现代修正: Murray 7.3 节末提到 Decroly-Goldbeter 1987 观察到 Dictyostelium 模型的 period doubling, 但真实 cAMP 振荡是单周期 (无 chaos)。 这提示"模型→真实" 的差距——模型的复杂数学行为不一定对应生物学现实。 Murray 没强调这一警示, 但 22 年后(2026) 这是计算生物学的共识: "模型行为的丰富 ≠ 真实生物的丰富"。

与自身研究的相关性: 我做血管生物力学, 关注的"血管平滑肌细胞 (SMC) 周期" 与"内皮细胞周期" 是与细胞周期相关的, 但更受关注的是血管壁的应力-生长律——血管重塑涉及持续的细胞增殖/凋亡/迁移, 这些过程在某些病理条件下会出现振荡 (如动脉粥样硬化的斑块形成-破裂-修复循环)。 Murray 7.6 节的延迟模型可以借鉴到"血管壁的反馈-延迟" 建模 (WSS → 内皮 → SMC → 血管直径变化 → WSS 变化), 周期约为数月到数年。 这是一个未充分开发的应用方向。

关于"周期 vs 可兴奋性" 的神经科学反思: Murray 7.5 节强调 FitzHugh-Nagumo 在 \(I_a = 0\) 时是可兴奋 (Ch 6 阈值现象的神经版), \(I_a > I_1\) 时是振荡 (极限环发放)。 这一连续过渡 在真实的神经元中也观察到——刺激电流增大时, 神经元从"单脉冲响应" 过渡到"周期发放"。 这种过渡的精确参数依赖 (Klaus-Dilão 2012) 是 22 年后(2026) 计算神经科学的标准知识, Murray 2002 的定性描述是起点。

对"反馈环长度" 的哲学反思: Tyson-Othmer \(m_0(n) = \sec^n(\pi/n)\) 提示长反馈环降低 Hill 阈值。 这与生物界真实存在的级联调控一致 (如 MAPK 级联、细胞周期级联)——演化"选择"了长反馈环来降低对单步调节精度的要求。 Murray 没在哲学层面展开这一观察, 但 22 年后(2026) "演化-工程" 类比 (Tsai 2008) 是热门方向。

§7 重要参考文献

[X1] Hodgkin, A.L., Huxley, A.F. (1952). "A quantitative description of membrane current and its application to conduction and excitation in nerve." Journal of Physiology 117: 500-544. (本章 7.5 节的核心原始文献, 4 维神经膜模型, 1952 年, 获 1963 年诺贝尔生理学/医学奖。)

[X2] FitzHugh, R. (1961); Nagumo, J., Arimoto, S., Yoshizawa, S. (1962). "An active pulse transmission line simulating nerve axon." Proceedings of the IRE 50: 2061-2070. (本章 7.5 节的简化模型 (7.39), 2 维 ODE 保留 Hodgkin-Huxley 的关键定性特征。)

[X3] Goodwin, B.C. (1965). "Oscillatory behavior in enzymatic control processes." In Advances in Enzyme Regulation, Pergamon. (本章 7.2 节 mRNA→酶→产物负反馈模型的原始来源, Griffith 1968 证明 \(m > 8\) 才有振荡。)

[X4] Tyson, J.J., Othmer, H.G. (1978). "The dynamics of feedback control circuits in biochemical pathways." Progress in Theoretical Biology 5: 1-62. (本章 7.2 节 \(m_0(n) = \sec^n(\pi/n)\) 的原始推导, 是负反馈振荡理论的奠基性工作。)

[X5] Tyson, J.J. (1979). "Periodic enzyme synthesis and the repressilator: a comment on the design of a biological clock." Journal of Theoretical Biology 80: 27-38. (本章 7.2 节末"周期估计 \(T \approx 2N/K\)" 的原始来源, 是 Goodwin 模型定量化的关键。)

[X6] Hastings, S., Tyson, J., Webster, D. (1977). "Existence of periodic solutions for negative feedback cellular control systems." Journal of Differential Equations 25: 39-64. (本章 7.2 节"周期解存在的一般结果" 的数学严格化, 是 Goodwin-Tyson 框架的解析基础。)

[X7] Schnackenberg, J. (1979). "Simple chemical reaction systems with limit cycle behaviour." Journal of Theoretical Biology 81: 389-400. (本章 7.4 节最简振荡器 (7.20) 的原始来源, trimolecular reaction \(2X + Y \to 3X\)。)

[X8] Murray, J.D. (1982). "Parameter space for limit cycle solutions of a simple chemical reaction model." IMA Journal of Applied Mathematics 29: 221-224. (本章 7.4 节末"参数化方法" (7.27) 的原始论文, 避免解三次方程。)

[X9] Goldbeter, A. (1996). Biochemical Oscillations and Cellular Rhythms: The Molecular Bases of Periodic and Chaotic Behaviour. Cambridge University Press. (本章 7.1 节引用的"生化振荡综述" 经典专著, 涵盖糖酵解、cAMP、细胞周期等。)

[X10] Lotka, A.J. (1910). "Contribution to the theory of periodic reaction." Journal of Physical Chemistry 14: 271-274. (本章 7.1 节末"化学振荡" 的理论开山, 早于 Bray 1921 的实验。)

[X11] Bray, W.C. (1921). "A periodic reaction in homogeneous solution and its relation to catalysis." Journal of the American Chemical Society 43: 1262-1267. (本章 7.1 节首次实验观察到的化学振荡——H2O2-碘酸盐, 50 年后才被广泛接受。)

[X12] Belousov, B.P. (1951, 1959). (俄文文献, Murray 引用) (本章 7.1 节末引用的"化学振荡" 实验里程碑, 是 Ch 8 Belousov-Zhabotinskii 反应的前身。)

[X13] Selkov, E.E. (1968). "Self-oscillations in glycolysis." European Journal of Biochemistry 4: 79-86. (本章 6.6 节 / 7.1 节糖酵解振荡的原始工作, 2 维 ODE 包含自催化 \(x^2 y\)。)

[X14] Smith, W.R. (1980). "A model for the control of testosterone secretion in the male." Journal of Theoretical Biology 87: 401-422. (本章 7.6 节睾酮反馈模型的原始来源, Smith-Murray 1980s 早期合作。)

[X15] Cartwright, M., Husain, H. (1986). "A model for the control of testosterone secretion in the male." Journal of Theoretical Biology 123: 239-250. (本章 7.6 节末"扩展延迟模型" 的来源, 包含更复杂的生理学。)

[X16] Decroly, O., Goldbeter, A. (1987). "From simple to complex oscillatory behaviour: analysis of a model of Dictyostelium discoideum." Proceedings of the National Academy of Sciences 84: 4143-4147. (本章 7.3 节末"period doubling→chaos" 的来源, 黏菌 cAMP 模型的复杂动力学。)

[X17] Elowitz, M.B., Leibler, S. (2000). "A synthetic oscillatory network of transcriptional regulators." Nature 403: 335-338. (Murray 没引用, 但 22 年后 (2026) 的"repressilator" 是 Goodwin 模型的合成生物学实现, 是 Goodwin-Tyson 框架的工程验证。)

[X18] Hanusse, P. (1972). "De l'existence d'un cycle limite dans un système de réactions chimiques." Comptes Rendus de l'Académie des Sciences 274: 1245-1247. (本章 7.4 节首"2 物种 ODE 极限环需要 trimolecular" 的原始证明, 是 Schnackenberg 模型的理论基础。)

[X19] Gardner, T.S., Cantor, C.R., Collins, J.J. (2000). "Construction of a genetic toggle switch in Escherichia coli." Nature 403: 339-342. (Murray 没引用, 但 22 年后 (2026) 的"synthetic toggle switch" 是 7.3 节"生物开关" 概念的工程实现, 是 Goodwin-Hill 框架的扩展。)

[X20] Murray, J.D. (1977). Lectures on Nonlinear-Differential-Equation Models in Biology. Oxford University Press. (Murray 7.2 节末引用的"延迟控制模型" 早期专著, 是 Ch 7 的前身, 也包含本章睾酮模型的思想萌芽。)