[论文解读] A dynamic p53-mdm2 model with delay kernel

Florin George Horhat, Mihaela Neamţu+1 more|ArXiv.org|Aug 22, 2006

Cancer-related Molecular Pathways参考文献 11被引用 3

一句话总结

本文提出了一种具有分布时滞和弱核函数的p53-Mdm2反馈环延迟微分方程模型，用于研究DNA损伤诱导的振荡动力学。通过霍普夫分支分析，表明该系统在时滞参数变化下表现出亚临界霍普夫分支，导致不稳定或稳定的周期振荡，为实验中观察到的哺乳动物细胞内p53振荡提供了更精细的动力学解释。

ABSTRACT

Specific activator and repressor transcription factors which bind to specific regulator DNA sequences, play an important role in gene activity control. Interactions between genes coding such transcripion factors should explain the different stable or sometimes oscillatory gene activities characteristic for different tissues. In this paper, the dynamic P53-Mdm2 interaction model with distributed delays and weak kernel, is investigated. Choosing the delay or the kernel's coefficient as a bifurcation parameter, we study the direction and stability of the bifurcating periodic solutions. Some numerical examples are finally given for justifying the theoretical results.

研究动机与目标

• 为实现更符合生物实际的动力学，采用分布时滞和弱核函数对p53-Mdm2反馈环进行建模。
• 通过霍普夫分支分析，研究系统表现出振荡行为的条件。
• 通过改变时滞或核函数参数，确定分岔周期解的方向、稳定性和周期。
• 利用合理的生物参数值进行数值模拟，验证理论结果。
• 相较于先前的离散时滞模型，通过引入分布时滞和更平滑的核函数加以改进。

提出的方法

• 为p53和Mdm2的mRNA及蛋白浓度构建包含分布时滞的延迟微分方程（DDE）系统。
• 采用弱核函数 $ k_i(s) = \frac{q^{p+1}_i}{p!} s^p e^{-q_i s} $，其中 $ p=0 $ 表示弱核，允许分布时滞。
• 将时滞 $ \tau $ 或核系数 $ q $ 视为分岔参数，分析平衡态的局部稳定性，并确定霍普夫分岔的条件。
• 应用中心流形理论和规范型方法（Hassard方法），计算分岔周期解的方向、稳定性和周期。
• 基于系统参数和时滞，推导出关键分岔量的显式表达式：$ \mu_2 $、$ \beta_2 $ 和 $ T_2 $。
• 使用Maple 9.5进行数值模拟，采用符合生物实际的参数值，可视化波形和相图。

实验结果

研究问题

• RQ1 在具有分布时滞的p53-Mdm2系统中，何种条件下会发生霍普夫分岔？
• RQ2 时滞核的选择（弱核与离散核）如何影响振荡解的稳定性和性质？
• RQ3 什么决定了分岔周期解的方向（超临界/亚临界）和稳定性（稳定/不稳定）？
• RQ4 振荡周期如何随时滞参数变化而改变？
• RQ5 该模型是否能比先前的离散时滞模型更准确地再现实验中观察到的p53振荡？

主要发现

• 对于离散时滞模型（$ \tau_2 = 3 $），系统在 $ \tau_0 = 9.54 $ 处发生亚临界霍普夫分岔，$ \mu_2 = -0.420 $，表明当 $ \tau > \tau_0 $ 时存在分岔解。
• 周期解为轨道不稳定（$ \beta_2 = 0.280 $），且随着时滞超过临界值，周期增加（$ T_2 = 0.0005 $）。
• 对于弱核模型（$ q_2 = 0.5 $），霍普夫分岔为亚临界（$ \mu_2 = -0.599 $），周期解为轨道稳定（$ \beta_2 = -0.748 $），且当 $ \tau_1 > 32.37 $ 时周期增加（$ T_2 = 0.180 $）。
• 该模型成功再现了与DNA损伤后p53动力学实验观察一致的振荡行为。
• 采用分布时滞与弱核函数可实现更平滑的动力学，产生比离散时滞模型更符合生物实际的振荡。
• 数值模拟证实了理论预测，显示p53和Mdm2蛋白水平呈现持续振荡，具有特征性波形和相平面图样。