延迟遗传调控网络稳定性及分支.pdfVIP

◎ 生物数学学报2008，234()：647-655B—ioJoIYul壁rnamll!of 延迟遗传调控网络稳定性及分支 王晶囡 ，。 蒋卫华 。 f1哈尔滨理工大学 应用数学系，黑龙江 哈尔滨 150080；2哈尔滨工业大学 理学院数学系。黑龙江 哈尔滨 150001) 摘 要： 本文主要研究了延迟遗传调控网络的局部稳定性和该网络的Hopf分支存在条 件．延迟遗传调控网络是无穷维系统，此类系统在平衡点线性化后的特征方程为超越方程。通过 对此超越方程进行研究，得到 了系统系数不同时的系统稳定的条件及相关结论，又进一步说明了 此系统的Hopf分支存在条件．最后，举一个例子进行了数值仿真验证了所得到的结论． 关键词：遗传调控网络；稳定性；延迟；Hopf分支 中图分类号：0175．14 MR分类号t92B20；34K20；34K15 文献标识码：A 文章编号：1001—9626(2008)04-0647-09 0 引 言 近些年来，由于遗传调控网络被广泛应用于生物工程中，如基因序列分析和基因识别等， 所以有很多人研究遗传网络 [1-6l，文 1【】给出了此类 自调节编译转录因子的遗传回路模型， 文 2【】研究的是此网络的频率选择，多稳定和能量振荡的内容，文 [3】的作者研究了延迟对二 维此网络的周期振荡的影响，文 [4】涉及的是遗传调控网络细胞分裂动力系统方面的知识．文 5【】用微分方程建立了此网络的数学表达式，由于在转录函数和翻译函数中都存在延迟，文 6【】 在此基础上，给出了延迟遗传调控网络的数学模型，在特殊情况下 (系数都相同)，分析了此 网络的局部稳定性，简要地描述了分支．本文是在文 f6】的基础上，去掉系统系数相同这一限 制，在较一般情况下，对网络进行了稳定性分析，给出了网络稳定性判别定理并明确的给出了 延迟遗传调控网络Hopf分支存在的条件，并利用计算机仿真检验了所得到的结果． 数学表达式 6【】 J (￡)=一Vm(t)+c(p(t一 ))， ， l ()=一，p()+d(m(t一7m．))， 、 这里m={ml，m2，…，m )∈R ，P= {pl，P2，…，P)∈R ，且 和 是正实数对角矩阵， 和 叩分别表示mRNA和蛋白质延迟量．m(t一 )={m1 一7m1)，…，m (t—Tm))， p(t—Tm)= {pl(t一丁m1)，…，Pn(—Tm))和 d(m)= {d1(m)，…， (m))∈R“，其中c)= {c1(p)，…，cn ))∈R”是非线性向量值函数． 收稿 日期：2007．01—12 基金项且：黑龙江省自然科学基金 (A200806)；哈尔滨工业大学理学研究基金 (HITC200704) 作者简介：王品囡(1978一)，女，硕士． E-mail：jnwang2006@yahoo．corn．cn 生 物 数 学 学 报 第 23卷 1 稳定性研究 假设在转录函数c(p)中只有一个蛋白质浓度做变量，同样在翻译函数 d(m)中也只有一个 mRNA浓度做变量．可以将 (1)式写成具体形式 {】 』帆()=一{mt()+勺t(-))， (21 I (￡)：一k~pi(t)+dt(m{(￡一7_m{))， 这里‰ 和 为对正的实的角阵．即 ：diag(k~nl，…， )， =diag(kpl，…，七pn)， 其中，七 {0，km0， (i，J=l，2，…，n)，J=i一1，当J=0时，且记J=札． 方程 (2)的不动点为 (m ，P)T，假设cj(pj)和di(mi)是满足光滑的函数，将 (t)=mt—m： 和玑(t)=Pi—P：代入 (2)在零点线性化，此线性化系统为 一 )+ 0