×

生物分子反应网络跳跃扩散近似的混合主方程。 (英语) Zbl 1443.60073号

总结:细胞反应具有多尺度性质,即分子物种的丰度和反应速率的大小可以在数量级上发生变化。这种多样性自然会导致混合模型,将连续和离散建模模式结合起来。为了捕捉这种多尺度特性,我们在之前的研究中提出了跳跃-扩散近似。其关键思想是将反应划分为快组和慢组,然后将慢组的马尔可夫跳跃更新方案与快组的扩散(Langevin)更新方案相结合。在本研究中,我们证明了反应计数过程中跳跃-扩散近似的联合概率密度函数满足混合主方程,该混合主方程结合了化学主方程和福克-普朗克方程的项。受条件矩方法的启发,我们提出了一种有效的方法来求解这个主方程,即在给定慢反应的反应计数器的情况下,利用快速反应的反应计数矩来求解主方程。然后,对于每个感兴趣的时间点,我们求解一组最大熵问题,以便从其矩中恢复条件概率密度。这最终使我们能够重建所有反应计数器上的完整联合概率密度,从而获得混合主方程的近似解。最后,我们展示了该方法应用于一个简单的多尺度转换过程的准确性。

MSC公司:

60J60型 扩散过程
92C45型 生化问题中的动力学(药代动力学、酶动力学等)
65C20个 概率模型,概率统计中的通用数值方法

软件:

CVX公司
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Abramov,RV,《多维最大熵矩问题:数值方法综述》,Commun。数学。科学。,8, 2, 377-392 (2010) ·Zbl 1193.49034号 ·doi:10.4310/CMS.2010.v8.n2.a5
[2] Anderson,D.F.,Kurtz,T.G.:化学反应网络的连续时间马尔可夫链模型。收录:Koepll,H.,Setti,G.,Bernardo,Md,Densmore,D.(编辑)生物分子电路的设计与分析。柏林施普林格出版社(2011)
[3] 安德烈琴科,A。;米基耶夫,L。;Wolf,V.,基于动量的随机化学动力学模型重建,ACM Trans。模型。计算。模拟。,2015年2月25日·Zbl 1384.92066号 ·doi:10.1145/2699712
[4] 博伊德,S。;Vandenberghe,L.,凸优化(2004),剑桥:剑桥大学出版社,剑桥·Zbl 1058.90049号
[5] Bretthorst,GL,矩的最大熵方法和贝叶斯概率理论,AIP Conf.Proc。,1553, 3-15 (2013) ·数字对象标识代码:10.1063/1.4819977
[6] Bronstein,L。;Koeppl,H.,生物分子反应网络动力学的矩-闭合近似变分方法,J.Chem。物理。,148, 1, 014105 (2018) ·数字对象标识代码:10.1063/1.5003892
[7] 谢瓦利埃,A。;Engblom,S.,空间随机动力学多尺度变量分裂方法的路径误差界,SIAM J.Numer。分析。,58, 1, 469-498 (2018) ·Zbl 1382.65013号 ·doi:10.1137/16M1083086
[8] 克鲁杜,A。;德彪西,A。;Radulescu,O.,多尺度基因网络的混合随机简化,BMC系统。生物学,389(2009)·doi:10.186/1752-0509-3-89
[9] Cseke,B。;施诺尔,D。;Opper,M。;Sanguinetti,G.,连续时间随机过程的期望传播,J.Phys。数学。理论。,49, 49, 494002 (2016) ·Zbl 1360.60150号 ·doi:10.1088/1751-8113/49/49/494002
[10] Doob,JL,Markoff链-可数外壳,Trans。美国数学。《社会学杂志》,58,3455-473(1945)·Zbl 0063.01146号 ·doi:10.2307/1990339
[11] 邓肯,A。;Erban,R。;Zygalakis,K.,《随机化学动力学模拟的混合框架》,J.Compute。物理。,326, 398-419 (2016) ·Zbl 1373.92041号 ·doi:10.1016/j.jp.2016.08.034
[12] Eldar,A。;Elowitz,MB,《遗传电路中噪声的功能作用》,《自然》,467167-173(2010)·doi:10.1038/nature09326
[13] Engblom,S.,计算主方程高维解的矩,应用。数学。计算。,180498-515(2006年)·Zbl 1103.65011号
[14] Engblom,S.、Hellander,A.、Lötstedt,P.:随机反应扩散网络的多尺度模拟。计算细胞生物学的随机过程、多尺度建模和数值方法,第55-79页。柏林施普林格出版社(2017)·Zbl 1401.92091号
[15] Fedoroff,N。;Fontana,W.,《大分子的小数目》,《科学》,2971129-1131(2002)·doi:10.1126/science.1075988
[16] 弗里德曼,N。;蔡,L。;Xie,X.,通过活细胞中的单分子实验观察到的基因表达的随机性,Isr。化学杂志。,49, 333-342 (2010) ·doi:10.1560/IJC.49.3-4.333
[17] Ganguly,A。;Altñtan,D。;Koeppl,H.,随机反应动力学的跳跃-扩散近似:误差界和算法,多尺度模型。模拟。,13, 4, 1390-1419 (2015) ·Zbl 1339.60088号 ·doi:10.1137/140983471
[18] 加德纳,CW,《随机方法手册》(1985),柏林:施普林格,柏林
[19] M.吉布森。;Bruck,J.,《具有多物种和多通道的化学系统的高效精确随机模拟》,J.Phys。化学。A、 104、9、1876-1889(2000)·doi:10.1021/jp993732q
[20] Gillespie,D.,《化学主方程的严格推导》,《物理学A》,188,404-425(1992)·doi:10.1016/0378-4371(92)90283-V
[21] Gillespie,DT,《化学朗之万方程》,J.Chem。物理。,113, 1, 297-306 (2000) ·doi:10.1063/1.481811
[22] Gillespie,DT,化学动力学的随机模拟,年。物理版。化学。,58, 35-55 (2007) ·doi:10.1146/annurev.physchem.58.032806.104637
[23] Grant,M.,Boyd,S.:CVX:Matlab软件,用于严格的凸规划,2.1版(2014)
[24] Grima,R。;托马斯,P。;Straube,AV,非线性化学Fokker-Planck方程和化学Langevin方程的精度如何?,化学杂志。物理。,135, 8, 084103 (2011) ·数字对象标识代码:10.1063/1.3625958
[25] Hasenauer,J。;Wolf,V.公司。;Kazeroonian,A。;Theis,FJ,《化学主方程的条件矩法(MCM):矩法和混合随机确定性模型的统一框架》,J.Math。生物学,69,3,687-735(2014)·Zbl 1302.92070号 ·doi:10.1007/s00285-013-0711-5
[26] Hellander,A。;Lötstedt,P.,化学主方程的混合方法,J.Compute。物理。,227, 1, 100-122 (2007) ·Zbl 1126.80010号 ·doi:10.1016/j.jcp.2007.07.020
[27] Jahnke,T.,关于化学主方程的简化模型,多尺度模型。模拟。,9, 4, 1646-1676 (2011) ·Zbl 1244.65005号 ·doi:10.1137/10821500
[28] Jaynes,ET,信息理论和统计力学,物理。修订版,106、4、620(1957年)·Zbl 0084.43701号 ·doi:10.1103/PhysRev.106.620
[29] Kampen,N.G.v.:马尔可夫过程的扩散近似。摘自:《生物过程的热力学和动力学》,第185-195页。Walter de Gruyter公司(1982年)
[30] Kurtz,TG,密度相关Markov cahins的强逼近定理,Stochast。过程。申请。,6, 3, 177-191 (1978) ·Zbl 0373.60085号 ·doi:10.1016/0304-4149(78)90020-0
[31] 李,CH;Kim,KH;Kim,P.,随机反应网络的矩闭包方法,J.Chem。物理。,130, 134107 (2009) ·数字对象标识代码:10.1063/1.3103264
[32] Menz,S。;拉托雷,JC;舒特,C。;Huisinga,W.,化学主方程的混合随机确定性解,多尺度模型。模拟。,10, 4, 1232-1262 (2012) ·兹比尔1302.92047 ·doi:10.1137/10825716
[33] 奥顿,DL;Vedula,P.,非线性Fokker-Planck方程基于矩的求积方法,J.Stat.Mech。理论实验,2011,9,09031(2011)·doi:10.1088/1742-5468/2011/09/P09031
[34] Pawula,RF,Fokker-Planck-Kolmogorov方程的推广和扩展,IEEE Trans。Inf.理论,13,1,33-41(1967)·Zbl 0154.42401号 ·doi:10.1109/TIT.1967.1053955
[35] 皮克勒,L。;马苏德,A。;洛杉矶伯格曼;Papadrakakis,M。;Stefanou,G。;Papadopoulos,V.,《用有限差分法和有限元法对福克-普朗克方程进行数值求解——一项比较研究》,《随机动力学中的计算方法》,《应用科学中的计算法》,69-85(2013),多德雷赫特:施普林格·Zbl 1304.65222号
[36] Risken,H。;Haken,H.,《福克-普朗克方程:求解方法和应用》(1989),柏林:施普林格出版社,柏林·Zbl 0665.60084号
[37] Salis,H。;Kaznessis,Y.,耦合化学或生物化学反应系统的精确混合随机模拟,J.Chem。物理。,122, 054103 (2005) ·数字对象标识代码:10.1063/1.1835951
[38] 辛格,A。;Hespanha,JP,化学反应系统的近似力矩动力学,IEEE Trans。自动化。对照,56,2144-418(2011)·Zbl 1368.80011号 ·doi:10.1109/TAC.2010.2088631
[39] Whittle,P.,《关于在随机过程处理中使用正态近似》,J.R.Stat.Soc.Ser。B(Methodol.),19,2,268-281(1957)·Zbl 0091.30502号
[40] 沃尔夫,V。;Goel,R。;Mateescu,M。;TA Henzinger,使用滑动窗口求解化学主方程,BMC Syst Biol,4,42(2010)·doi:10.1186/1752-0509-4-42
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。