反馈、开关与振荡¶

调控网络画出了谁影响谁，时间过程则取决于这些作用以多快的速度发生。一个转录因子可以抑制自身，但蛋白质若产生得很慢、降解得更慢，回路就会长期保留过去的输入；同一条正反馈在弱协同性下只会放大响应，在足够非线性时才可能把细胞分成低、高清晰的两个状态；负反馈通常稳定输出，若抑制来得足够迟，又可能使系统反复越过稳态而形成振荡。系统动力学的任务，就是把网络拓扑、反应速率、时间尺度和初始状态放进同一幅因果图景。

常微分方程和相图提供分析负反馈、适应、正反馈、双稳态、滞后和生化振荡的共同语言。节点与边的生物学语义见生物网络结构与网络模体，分子数涨落、随机转换和空间耦合见随机动力学、空间模式与群体系统，参数估计见系统建模、参数推断与实验检验。

状态变量与速率方程¶

动力学模型先规定系统边界和状态变量。状态变量可以是蛋白质浓度、磷酸化比例、受体活性、代谢物池、细胞数或某种表型的占比；输入可以是配体、营养、温度或实验者施加的脉冲。参数则描述当前模型中相对固定的规则，例如最大合成速率、结合常数和降解常数。区分二者很重要：蛋白质浓度会沿轨迹变化，降解常数通常在一次模拟中保持不变；若实验者逐步改变降解速率，它又成为控制系统行为的分岔参数。

对单一组分 \(x\)，最基本的质量平衡是

其中 \(v_{+}\) 汇总生成或流入，\(v_{-}\) 汇总降解、失活或流出。最简单的线性模型取 \(v_{+}=\alpha u\)、\(v_{-}=\beta x\)，恒定输入下的稳态为 \(x^{*}=\alpha u/\beta\)，时间过程为

因此 \(1/\beta\) 给出回到稳态的特征时间，半衰期为 \(\ln 2/\beta\)。稳态水平和响应速度由不同参数组合决定：提高生成率会抬高稳态，提高降解率既降低稳态又缩短记忆。把稳态表达量相同的两个系统放在一起，快速脉冲实验仍可显示完全不同的动力学。

真实生化网络通常含有多个状态变量。反应的化学计量决定每个变量得到或失去多少物质，速率律决定变化有多快。常微分方程默认变量连续、体系在所取区室内充分混合，适合分子数较多且空间梯度不占主导的情形；低拷贝随机性、扩散和细胞间耦合需要另加描述。简化模型的价值在于保留决定行为的回路和时间尺度，并不要求把每个中间体都写入方程。¹

稳态、扰动与局部稳定性¶

稳态是所有状态变量的净变化率同时为零。对一维系统，\(f(x^{*};u)=0\) 表示生成与清除相等。若在 \(x^{*}\) 左侧有 \(f>0\)、右侧有 \(f<0\)，状态被推回交点，这个稳态局部稳定；两侧箭头都指离交点时，微小偏离会被放大，它便不稳定。把 \(x=x^{*}+\delta x\) 代入并保留一阶项，可得

\(f'(x^{*})<0\) 时扰动衰减，\(f'(x^{*})>0\) 时扰动增长。这里的“稳定”只说明足够小的扰动会怎样变化，不等于系统不会被大扰动推到另一个吸引域，也不等于稳态对参数变化不敏感。

两个变量满足 \(\dot{x}=F(x,y)\)、\(\dot{y}=G(x,y)\) 时，\(F=0\) 和 \(G=0\) 分别形成两条零流线。交点是稳态，交点附近的速度方向由 Jacobian 矩阵

决定。对一般的二维连续系统，\(\operatorname{tr}J<0\) 且 \(\det J>0\) 给出局部稳定条件；实特征值常对应单调趋近，复特征值可对应衰减振荡。轨迹围绕稳态转几圈后停止，是阻尼振荡；收敛到闭合轨道并持续循环，才是极限环。相图由此把“最后到哪里”和“沿什么路径到达”同时表达出来。MIT 系统生物学课程也用生成—降解曲线、零流线和稳定性分析串联自调控、双稳态与振荡。²

非线性与阈值响应¶

分子结合、酶促修饰和转录调控很少在全部浓度范围内保持线性。常用 Hill 函数近似激活与抑制：

\(K\) 决定响应曲线的中点，\(n\) 描述曲线陡峭程度。拟合得到的 Hill 系数是有效协同性指标，不必等于复合物中配体数或结合位点数。多位点结合、零级超敏性、多步级联、底物竞争和分子滴定都能产生陡峭输入—输出关系。

超敏性意味着输入在较窄区间变化时，输出跨过较大的动态范围；它本身仍可是一条单值曲线。阈值也不自动等于开关：一个单稳态系统可以响应很陡，却在输入降低时沿同一路径返回。双稳态还需要同一参数条件下存在两个稳定状态，滞后则要求上调和下调输入时跨越不同阈值。把这些现象分开，才能从剂量—响应曲线判断回路究竟是在放大、选择状态，还是保留历史。

负反馈的增益、速度与稳态¶

负反馈把输出的增加送回上游并抑制净生成。它常降低低频扰动对输出的影响，限制不同细胞或不同参数背景造成的稳态偏差，也可加快某些响应。负自调控就是最小例子：蛋白质尚低时生成较快，接近目标水平后生成被压低；相同最终表达量下，它可比无反馈系统更快进入稳态，相关合成回路实验已在自调控回路中介绍。

反馈效果取决于环路增益和各环节速度。增益过低时校正有限；增益提高可以减小慢扰动造成的偏差，却可能放大测量噪声、降低对真实输入的响应幅度，或在延迟存在时造成过冲。降解、转运、转录和翻译各有时间尺度，整个闭环的响应需要由各环节共同描述。负反馈是一种具有条件与时间尺度的动态结构。

比例反馈根据当前误差立即校正，通常仍留下与扰动大小有关的稳态误差。积分反馈则累积过去的误差。若控制变量 \(z\) 满足

只要闭环能到达稳定稳态，就有 \(\dot z=0\)，从而输出 \(y\) 回到目标值 \(y_{0}\)，即使恒定输入或部分参数已经改变。这种结构给“适应”增加了严格含义：系统在持续刺激下先产生瞬态响应，随后回到刺激前或规定的输出基线；它与受体彻底失去响应、输出耗竭或输入被清除需要实验区分。

细菌趋化的鲁棒精确适应¶

细菌趋化系统在吸引物阶跃增加后先降低翻滚概率，随后在吸引物仍存在时把受体—CheA 输出拉回接近原水平，使细胞重新对浓度变化敏感。受体的甲基化与去甲基化提供较慢的记忆变量；受体活性改变 CheB 等调节环节，形成与快速信号支路相反的反馈。Yi 等分析 Barkai–Leibler 模型时证明，其中具有积分反馈结构，并据此解释精确适应为何可在多种蛋白丰度和速率常数变化下保持。³

这里的鲁棒性质特指“稳态输出回到基线”。反馈增益、受体丰度和甲基化速度仍会改变峰值、适应时间及梯度追踪表现；输入若超出受体和修饰系统的动态范围，精确适应也会失效。积分反馈给出了稳态误差归零的结构条件，闭环稳定性和生化可实现范围仍需另行检验。

趋化还说明动态实验为何比终点测量信息更多。只比较刺激前后两个稳态，可能看不见中间的显著响应；只看到暂时回落，又不能证明对不同输入强度和参数扰动都保持精确。阶跃刺激、去除刺激、改变适应蛋白水平并连续测量通路输出，才能把感知、瞬态响应、反馈校正和鲁棒性连成一条证据链。

正反馈、双稳态与滞后¶

正反馈把输出增加送回生成过程。对一个自激活因子，可写成

左边的生成曲线若不够陡，只与线性降解曲线相交一次，系统仍是单稳态；反馈足够强且响应足够非线性时，可出现三个交点。低状态和高状态为稳定点，中间点不稳定，并成为两个吸引域的分界。小扰动在原吸引域内衰减，跨过分界的大扰动则把系统送入另一稳定状态。

随外部输入缓慢增加，低稳态可在一个鞍结分岔处消失，状态突然跳到高支；反向降低输入时，高稳态会在另一个阈值才消失。上、下阈值之间的同一输入支持两个稳定输出，当前状态取决于输入历史，这就是滞后。滞后环是双稳态的重要证据，但测量时必须慢到让系统接近准稳态，并排除蛋白长寿命、不可逆损伤、选择性死亡或群体成分变化造成的表观记忆。

大肠杆菌 lac 系统把通透酶摄取诱导物、诱导物解除 LacI 抑制以及更多通透酶产生连成正反馈。Ozbudak 等用单细胞报告和输入相图测量到低、高清晰状态及滞后，并通过改变碳源与反馈条件把滞后响应转成超敏但连续的响应，直接显示“有反馈”还不足以保证多稳态，参数区间决定实际行为。⁴ lac 操纵子的分子组成与碳源组合逻辑见原核与真核基因表达调控。

细胞命运开关的长期记忆¶

非洲爪蟾卵母细胞成熟把连续的孕酮信号转换为单细胞层面的全或无丝裂原活化蛋白激酶（mitogen-activated protein kinase，MAPK）活化和成熟响应。Ferrell 与 Machleder 的实验显示，中等激素剂量下，细胞群由未激活与充分激活的卵母细胞组成，群体平均曲线会掩盖这种单细胞开关。⁵

随后对 Mos–MEK–MAPK 回路的扰动和持续性测量表明，正反馈可形成双稳态记忆模块：足够强的启动信号把系统推入高 MAPK 状态，信号撤除后回路仍能维持成熟状态；破坏关键反馈会削弱这种记忆。⁶ 这类状态维持适合完成需要一致承诺的细胞命运转换，但“双稳态”仍不等于物理上永远不可逆。蛋白降解、细胞周期推进、反馈边改变或新的上游输入都可让参数越过分岔边界，重新安排可达状态。

细胞周期蛋白依赖性激酶 1（cyclin-dependent kinase 1，CDK1）与周期蛋白 B（cyclin B）的控制又把开关嵌入振荡器。活化 CDK1 促进其激活因子 CDC25，并抑制 Wee1/Myt1 类抑制因子，形成正反馈与双负反馈，使进入 M 期快速而协同。非洲爪蟾卵提取物实验观察到 CDK1 激活的滞后：启动有丝分裂所需的 cyclin B 水平高于维持已激活状态所需水平。⁷ 具体周期蛋白合成、后期促进复合体／周期体（anaphase-promoting complex/cyclosome，APC/C）降解和检查点机制见细胞周期与细胞分裂；这里关注这些分子组成的动力学开关。

延迟负反馈与持续振荡¶

一个没有显著延迟的简单负反馈常把偏离稳态的状态拉回去。若输出上升后，抑制信号要经过转录、翻译、逐级修饰或区室转运才返回，系统会在反馈到达前继续上升；抑制建立后又可能下降过头，反馈撤除同样迟到，于是反复越过稳态。理想化的延迟方程可写为

其中 \(\tau\) 是回路记忆。延迟必须与足够陡的非线性和合适的生成—清除时间尺度配合；否则轨迹只会单调回稳态或出现逐渐衰减的振荡。Novák 与 Tyson 将生化振荡的共同条件归纳为负反馈、足够的非线性、足够的记忆以及相配的反应时间尺度。⁸

参数变化可使稳定焦点失去稳定，并在附近产生极限环，这类转变常由 Hopf 分岔描述。极限环不是一串由外界重复触发的脉冲，而是恒定条件下系统经过暂态后自行保持的周期轨道。轨道上的状态在一次循环后返回相同位置，邻近轨迹则回到同一闭合轨道；初始状态主要改变最初经过的路径和相位，而长期振幅与周期由回路和参数决定。

Hes1 是延迟负反馈的经典实例。Hes1 蛋白抑制自身转录，Hes1 mRNA 和蛋白又都较不稳定，使旧信号能在下一轮开始前被清除。时间序列实验观察到约两小时尺度的 Hes1 mRNA 与蛋白振荡；阻断降解会延长组分寿命并破坏正常振荡，说明负反馈拓扑之外，转录—翻译延迟和清除速度同样是机制的一部分。⁹ 发育中的 HES/Notch 节律如何与细胞耦合并形成体节边界，见脊椎动物体节时钟。

嵌套反馈与细胞周期振荡器¶

早期胚胎细胞周期以 cyclin B 的合成与降解形成核心负反馈：cyclin B 积累并激活 CDK1，CDK1 促进 APC/C 激活，APC/C 随后降解 cyclin B，使 CDK1 活性下降并重置循环。若所有反应都缓慢、近线性地变化，这条负反馈可停在中间稳态。CDK1 开关的正反馈提供陡峭阈值和滞后，让 cyclin 积累先跨过较高的激活阈值，再由 APC/C 驱动下降并跨过较低的失活阈值，形成松弛振荡。

Tsai 等以非洲爪蟾胚胎细胞周期回路为原型，比较纯负反馈与正负反馈嵌套模型，并用随机参数采样检验不同结构的工作范围；嵌套结构扩大了支持振荡的参数区域，也使频率能够在较宽范围内变化而振幅相对稳定。¹⁰ 这不意味着每个生物钟都采用同一分子元件，而是说明两个回路可以分工：正反馈建立快速、可靠的状态转换，较慢负反馈负责重置。相似逻辑可由转录抑制、酶活化、蛋白降解或离子通量等不同分子机制实现。

昼夜钟、细胞周期、钙振荡和发育分节都可产生节律，却读取不同的时间尺度与输出。昼夜系统还要接受光暗、进食和温度等周期输入，把内源周期锁定到环境；其生理输出见松果体褪黑素与昼夜时相。把外部周期与内源振荡器区分开来很重要：周期输入可以迫使一个本不自主振荡的系统产生周期响应，也可使已有振荡器发生同步。

周期、振幅、相位与同步¶

振荡至少要描述四个量。周期 \(T\) 是相邻同相位事件之间的时间，频率为 \(1/T\)；振幅描述一个周期内偏离基线的范围；相位给出系统在周期中的位置；占空比则表示输出处于高状态的时间比例。不同回路可以保持周期而改变振幅，也可以以近似固定振幅改变频率。只报告“有振荡”会掩盖这些可独立调节的输出维度。

短暂扰动可提前或推迟下一次峰值，由扰动施加时的相位决定。把相位移动对初始相位作图得到相位响应关系，它能够定位振荡器对输入最敏感的窗口。周期性外部信号接近内源频率并具有足够强度时，振荡器可按固定相位差与外界同步；光暗同步昼夜钟、相邻细胞通过 Notch 信号同步分节时钟，都属于相位动力学问题，具体耦合和空间波将在后续随机—空间页面展开。

群体平均尤其容易误判节律。单细胞若各自持续振荡但相位逐渐分散，群体信号会表现为振幅衰减；反过来，细胞群中不同亚群交替占优也可能制造周期性平均值。长时间单细胞追踪、谱系信息和相位分布能够区分单细胞阻尼、相位去同步与真正的极限环。荧光报告蛋白的成熟和降解时间还会滤去快速变化，因此观测到的周期与相位要和报告系统本身的动力学一起解释。

鲁棒性的扰动对象与保持性质¶

鲁棒性描述某项输出对某类扰动的敏感程度。趋化可对蛋白丰度变化保持适应终点，却改变响应速度；细胞周期振荡器可在参数变化下保持振幅，却改变周期；双稳态开关可保持两个状态，却让切换阈值移动。分析时要写明扰动的是输入、初始状态、参数还是网络结构，以及保持的是稳态值、峰值、周期、相位还是状态数目。

反馈会带来取舍。高增益负反馈能压低慢扰动造成的偏差，却可能增加能量与分子周转、传递高频测量噪声，并在延迟下损害稳定性；陡峭正反馈能提高状态分离和记忆，却可能放大细胞间参数差异，或让错误越阈值后的恢复更困难；长延迟有利于振荡，却也增加周期对生成和运输时间的依赖。多重反馈、冗余通路和模块耦合可以拓宽可工作参数区间，同时引入新的串扰与资源约束。

几何分析能直观看见这类取舍。参数变化移动或弯曲零流线，稳态随交点移动；两个交点在鞍结分岔处相遇并消失，产生开关；稳定点在 Hopf 分岔附近失稳，产生振荡。分岔图表达的是模型在参数空间中的行为区域，不是对每个细胞必然命运的宣告。细胞间参数分布、分子噪声和缓慢环境变化会让群体跨越不同边界，确定性相图需要与单细胞数据共同解释。

动力学的输入波形与边级扰动¶

动力学模型应接受能够区分机制的输入。阶跃输入适合测量响应时间、过冲和适应终点；短脉冲检验系统是否具有持续性检测或记忆；缓慢上、下扫描输入可寻找不同切换阈值；周期输入则测量频率响应、相位移动和同步区间。终点剂量—响应只能看到轨迹落点，不能替代这些时间信息。

对开关，需要在单细胞层面同时测量状态分布和历史。先从低态与高态出发，在相同输入下观察是否长期维持两个群体，再以反馈边突变、诱导物转运改变或蛋白降解调节移动阈值，可区分真正双稳态、缓慢弛豫与群体异质性。对振荡器，则应在恒定条件下证明多轮持续循环，改变延迟、非线性或负反馈强度并检验周期、振幅和稳定区是否按模型移动。

边级扰动通常比删除整个节点更有解释力。改变一个结合位点可削弱特定反馈而保留蛋白其他靶标；替换降解标签可选择性改变时间尺度；快速可逆的光控或化学扰动能在不同相位施加同样输入，测量相位响应。模型若能预言此前未测条件下的切换阈值、适应时间或振荡周期，并由独立实验验证，才从事后拟合走向机制检验。模型选择、参数可辨识性、灵敏度和扰动设计将在系统建模、参数推断与实验检验中继续。

参考资料与延伸阅读¶