[论文解读] On the quasi-steady-state approximation in an open Michaelis--Menten reaction mechanism

Justin Eilertsen, Marc R. Roussel+2 more|arXiv (Cornell University)|Jan 14, 2021

Nonlinear Dynamics and Pattern Formation参考文献 50被引用 18

一句话总结

本文通过奇异摄动理论与不变流形理论，严格分析了具有底物输入的开放Michaelis-Menten机制中的拟稳态近似（QSSA）。研究识别出特征慢流形，并在标准奇异摄动假设之外建立了QSSA有效性的条件，表明即使不存在临界流形，慢流形依然可以存在并吸引轨迹，从而在精确的时间尺度估计与高阶修正下，拓展了开放生化系统中QSSA的理论基础。

ABSTRACT

The conditions for the validity of the standard quasi-steady-state approximation in the Michaelis--Menten mechanism in a closed reaction vessel have been well studied, but much less so the conditions for the validity of this approximation for the system with substrate inflow. We analyze quasi-steady-state scenarios for the open system attributable to singular perturbations, as well as less restrictive conditions. For both settings we obtain distinguished invariant slow manifolds and time scale estimates, and we highlight the special role of singular perturbation parameters in higher order approximations of slow manifolds. We close the paper with a discussion of distinguished invariant manifolds in the global phase portrait.

研究动机与目标

• 将拟稳态近似（QSSA）的理论依据从封闭系统扩展至具有底物输入的开放反应机制。
• 研究在标准奇异摄动假设不成立时，开放Michaelis-Menten系统中不变慢流形的存在性与结构。
• 通过推导时间尺度估计与慢流形的高阶修正，量化QSSA约化的精度。
• 分析开放系统的全局相图，包括利用Poincaré球分析无穷远处的行为，并通过Poincaré变换对无穷远处的平衡点进行分类。

提出的方法

• 将奇异摄动理论应用于开放Michaelis-Menten机制，识别出与快速酶-底物结合和慢速产物形成相关的微小参数ε。
• 利用Tikhonov–Fenichel参数值（TFPV）确定系统存在非孤立平衡点、形成临界流形的参数区域。
• 采用正规型理论与中心流形约化方法，计算在临界流形不存在时的慢流形高阶近似。
• 利用Poincaré球分析全局动力学，包括无穷远处的行为，并通过Poincaré变换对无穷远处的平衡点进行分类。
• 应用Butler-McGehee定理与Poincaré-Bendixson理论，研究ω-极限集及平衡点的全局吸引性。
• 推导出至三阶与四阶的显式正规型展开，以表征不变流形的稳定性和几何结构。

实验结果

研究问题

• RQ1 在标准奇异摄动框架不适用时，具有底物输入的开放Michaelis-Menten机制中，拟稳态近似在何种条件下仍有效？
• RQ2 当临界流形不存在时，是否存在特征不变慢流形并能吸引附近轨迹？
• RQ3 当系统不适用于标准Tikhonov-Fenichel约化时，高阶修正在近似慢流形中起什么作用？
• RQ4 开放系统的全局相图，包括无穷远处的动力学，如何影响QSSA的有效性与精度？
• RQ5 慢流形的精确几何与动力学结构是什么？其与平衡点稳定性及ω-极限集的关系如何？

主要发现

• 当满足某些通用性条件时，即使标准奇异摄动框架失效，开放Michaelis-Menten系统仍存在特征不变慢流形。
• 当临界流形不存在时，慢流形依然存在并吸引邻近轨迹，表明QSSA可在经典Tikhonov-Fenichel理论之外成立。
• 在原点处的正规型展开至三阶为 ˙y3 = (k2eT − k0)y3³ + ...，表明原点的稳定性取决于 (k2eT − k0) 的符号，当 k2eT = k0 时为退化吸引节点。
• 在无穷远处，系统在P2处表现出鞍结点，上半球面具有排斥节点部分，且原点不稳定流形中轨迹的ω-极限集必须包含P1或P3的对径点，排除排斥节点。
• 当 k2eT = k0 时，四阶正规型为 ˙y3 = −(k−1 + k2)k2eT/k1 y3⁴ + ...，确认P1为退化吸引节点，并对第一象限内所有初值全局吸引。
• 全局相图被完全刻画：当 k2eT < k0 时，第一象限内所有轨迹收敛至原点P0；当 k2eT > k0 时，所有轨迹收敛至P1，且在后一条件下P1为全局吸引。