基于Markov链的重组细胞恒化培养的随机建模分析.pdf

第34卷 第2期 工 程 数 学 学 报 v01．34No．2 2017年04月 CHINESEJOURNALOFENGINEERINGMATHEMATICS Apr．2017 doi：10．3969／j．issn．1005—3085．2017．02．001 文章编号：1005—3085(2017)02—0111—13 基于Markov链的重组细胞恒化培养的随机建模分析木 李小月，姬雪晖，原三领十 (上海理工大学理学院，上海 200093) 摘 要：本文研究了一类重组细胞恒化培养的连续时间Markov链模型．首先利用累积母函数 表示出数字特征所满足的矩方程，然后通过对数正态分布近似的矩封闭技术得到了封 闭后的矩方程，最后运用Euler—Maruyama方法构建了时间和状态都是连续的It6随机 微分方程．为了验证矩封闭的合理性，利用数值模拟给出了确定模型、随机模型和矩 封闭后的方程的比较，并分析了重组细胞的变化趋势，结果表明其随机游走趋势与相 应确定性模型是一致的． 关键词：连续时间Markov链；累积母函数；对数正态分布；矩封闭；It5随机微分方程 分类号：AMS(2000)60H10；60J05；60J25 中图分类号：O211．6 文献标识码：A 1 引言 恒化器是实验室中设法使培养液流速保持不变，且微生物始终在低于其最高生长速率 条件下生长繁殖的一种连续培养装置，最初 由Monod创制[1j2】．在对微生物进行连续培养 的过程 中，实验人员可 以通过控制营养物质的输入率和稀释率来调节微生物的增长速度 ， 从而提高微生物培养的效率，其在生物发酵工程方面有着非常重要的作用．因此，以恒化 器为研究背景的动力学建模研究吸引着国内外许多生物学者、实验室技术人员和数学工作 者的广泛关注[3—61． 基因重组是改变生物原有遗传特性，获得新品种的生物技术 ．应用基因工程技术可 以 进行大规模 的生产，但是在生产过程 中却存在一个重大的问题 ：提取出来的DNA大分子 质粒不稳定．这种不稳定性会使携带质粒的微生物变成不携带质粒的微生物，而携带质粒 微生物拥有我们需要的优 良特性．因此，在现实培养过程中需要创造有利于携带质粒微生 物的竞争环境，从而获得更多携带质粒的微生物 ．恒化器中携带质粒与不携带质粒微生物 的竞争可 以由下面的模型来刻画 7【18] (1) 收稿日期：2015—06—23．作者简介：李小月(1990年5月生)，女，硕士生．研究方向：生物数学． 基金项目：国家自然科~ 11671260)；中国沪江基金(B14005)；上海市教委重点创新项目(13ZZl16) 十通讯作者：原三领 E-mail：sanling@usst．edu．ca l12 工 程 数 学 学 报 第34卷 其中s(t)表示t时刻营养液的浓度， 1(t)，X2(t)分别表示t时刻恒化器 中携带质粒和 不携带质粒的微生物的浓度，D为稀释率，S0为输入营养 的浓度 ，q(0q 1)为 携带质粒细胞在 自我复制过程中因为失去质粒而变成不携带质粒细胞的比率，即丢 失率，fl(8)和 。(s)分别为两种微生物的功能性反应函数．Hsu等[718]详细地研究了模 型f11具有抑制或非抑制功能性反应函数情况下的全局动力学行为．在模型 (1)的基础 上，Yuan和Zhang[0l考虑了周期营养输入和营养循环时滞并研究了模型的动力学行为． 尽管模型f1)能够很好地反映重组细胞培养的增长过程 ，但在微生物的培养过程中存 在许多不确定因素如温度、光照、测量仪器的精确度等，这些因素会对微生物培养过程产 生一定的影响[10】．因此，利用随机模型来研究携带质粒与不携带质粒的重组质粒微生物培 养显得尤为必要．本文基于确定性模型f1)，考虑随机因素对微生物培养过程的影响，建 立并研究了连续时间Markov链模型，即生灭