潘多濤，史洪巖，袁德成，修志龍

?

釀酒酵母代謝過(guò)程的振蕩分析

潘多濤1,2，史洪巖2，袁德成2，修志龍1

（1大連理工大學(xué)生命科學(xué)與技術(shù)學(xué)院，遼寧大連116024；2沈陽(yáng)化工大學(xué)遼寧省化工過(guò)程控制技術(shù)重點(diǎn)實(shí)驗(yàn)室，遼寧沈陽(yáng) 110142)

振蕩現(xiàn)象是生物的固有特征，在不同生物動(dòng)態(tài)過(guò)程中都起著重要作用。為探究生物代謝過(guò)程中振蕩產(chǎn)生的條件，以典型的釀酒酵母發(fā)酵生產(chǎn)乙醇過(guò)程為研究對(duì)象，采用函數(shù)連續(xù)性分析方法對(duì)其數(shù)學(xué)模型進(jìn)行深入研究。首先對(duì)系統(tǒng)進(jìn)行仿真，結(jié)合相圖確認(rèn)該過(guò)程存在等幅持續(xù)振蕩現(xiàn)象（極限環(huán)振蕩）；以此為切入點(diǎn)，利用分叉分析方法考察模型參數(shù)對(duì)系統(tǒng)振蕩現(xiàn)象的影響；最后根據(jù)結(jié)果分析系統(tǒng)產(chǎn)生振蕩現(xiàn)象的條件。結(jié)果顯示，該代謝途徑中各步反應(yīng)均有參數(shù)會(huì)對(duì)系統(tǒng)的振蕩現(xiàn)象產(chǎn)生影響，得出產(chǎn)生極限環(huán)振蕩的參數(shù)范圍，同時(shí)根據(jù)結(jié)果分析，為今后抑制振蕩或利用振蕩有利特性提供指導(dǎo)。

系統(tǒng)生物學(xué)；糖酵解；振蕩；極限環(huán)；分支分析；非線性規(guī)劃

引 言

系統(tǒng)生物學(xué)是通過(guò)建模分析，結(jié)合理論實(shí)踐，系統(tǒng)研究生物過(guò)程的學(xué)科，其首要目標(biāo)是在現(xiàn)有實(shí)驗(yàn)數(shù)據(jù)的基礎(chǔ)上，建立一個(gè)能夠準(zhǔn)確預(yù)測(cè)生物過(guò)程的模型。隨著生物學(xué)研究進(jìn)入“后基因組”時(shí)代，其關(guān)鍵的問(wèn)題是如何解決生物過(guò)程網(wǎng)絡(luò)的復(fù)雜性，如生物系統(tǒng)經(jīng)過(guò)長(zhǎng)期進(jìn)化形成的魯棒性、高度非線性以及多回路反饋等特點(diǎn)[1-2]。

生物系統(tǒng)中普遍存在的振蕩現(xiàn)象是生物過(guò)程復(fù)雜性的另一種表現(xiàn)形式。生物振蕩是生物體不斷與外界發(fā)生物質(zhì)、能量交換，在特定階段會(huì)處于偏離穩(wěn)態(tài)，發(fā)生非線性、非平衡的狀態(tài)周期變化的現(xiàn)象[3]。由外部因素引發(fā)的振蕩稱之為強(qiáng)迫振蕩，而源于內(nèi)因的稱之為自發(fā)性振蕩。國(guó)內(nèi)外學(xué)者已對(duì)生物體系各方面的振蕩現(xiàn)象做了大量研究，如Longo等[4]考察了誘導(dǎo)因子信號(hào)轉(zhuǎn)導(dǎo)網(wǎng)絡(luò)的振蕩行為。而在代謝過(guò)程方面，特別是在連續(xù)發(fā)酵過(guò)程中的振蕩行為，越來(lái)越受到重視。Williamson等[5]研究了酵母的基因控制對(duì)代謝振蕩過(guò)程的影響；Gustavsson等[6]首次觀察到游離酵母細(xì)胞的振蕩現(xiàn)象；Preez等[7]則對(duì)已有的糖酵解振蕩模型的參數(shù)可適范圍進(jìn)行優(yōu)化，并校正模型。由于振蕩過(guò)程的復(fù)雜性，僅通過(guò)數(shù)學(xué)模型很難直觀地闡明其動(dòng)態(tài)特性，更為有效的方法是通過(guò)調(diào)整模型參數(shù)來(lái)分析系統(tǒng)的動(dòng)態(tài)行為。

本文針對(duì)釀酒酵母在連續(xù)培養(yǎng)中的振蕩過(guò)程，在模型仿真分析的基礎(chǔ)上，利用非線性理論的方法對(duì)關(guān)鍵參數(shù)進(jìn)行分支分析，以期通過(guò)對(duì)系統(tǒng)等幅持續(xù)振蕩[8]的特性分析能夠進(jìn)一步為實(shí)驗(yàn)研究提供指導(dǎo)依據(jù)。

1 振蕩過(guò)程

振蕩是生物體固有的重要性質(zhì)之一，在體系內(nèi)的不同水平層次都有所體現(xiàn)，其宏觀表現(xiàn)形式一般被稱為節(jié)律[9]，如心肺張弛、生殖周期以及不同表現(xiàn)形式的生物鐘現(xiàn)象。而在生物細(xì)胞內(nèi)的代謝過(guò)程中出現(xiàn)的振蕩現(xiàn)象，通常是代謝物濃度出現(xiàn)周期性的波動(dòng)。

酵母發(fā)酵產(chǎn)生乙醇的過(guò)程是指以葡萄糖為代表的六碳糖共同經(jīng)歷的分解代謝生成丙酮酸并進(jìn)一步轉(zhuǎn)化為乙醇的過(guò)程。該過(guò)程中包含的糖酵解過(guò)程是生物界最原始獲取能量（ATP）的一種方式，是生物體共同經(jīng)歷的葡萄糖分解代謝的前期途徑，被認(rèn)為是大多數(shù)生物細(xì)胞代謝過(guò)程中最重要的途徑。在整個(gè)過(guò)程中，ATP-ADP以及NAD-NADH形成反饋回路，在一定條件下該過(guò)程就會(huì)表現(xiàn)出復(fù)雜的非線性動(dòng)力學(xué)特性[10]。在酵母發(fā)酵過(guò)程中，早在1957年Duysens等[11]就首次觀察到完整細(xì)胞中NADH的濃度波動(dòng)現(xiàn)象；之后陸續(xù)有人報(bào)道在無(wú)細(xì)胞的酵母提取物中存在糖酵解振蕩現(xiàn)象[12]；在糖酵解途徑中，磷酸果糖激酶催化的代謝物出現(xiàn)振蕩現(xiàn)象[13]。

2 方法

2.1 模型

生物代謝過(guò)程的數(shù)學(xué)模型一般采用微分方程（ordinary differential equation, ODE）來(lái)描述。以代謝物濃度為狀態(tài)變量的一組微分方程的向量形式如下

其中，=[12,…,x]T∈為系統(tǒng)狀態(tài)變量；=[1,2,…,k]∈為模型參數(shù)。式（1）描述了代謝物濃度隨時(shí)間的變化規(guī)律，一般以Monod或Michaelis-Menten方程描述其關(guān)系。為了盡可能提高模型的精確度，方程的形式越發(fā)復(fù)雜；代謝網(wǎng)絡(luò)的拓展也使需要考察的代謝物增多，對(duì)應(yīng)的系統(tǒng)模型維度也隨之不斷提高。這些因素導(dǎo)致模型整體的復(fù)雜度顯著增加，需要考慮更多的因素與非線性行為的關(guān)系。

長(zhǎng)期以來(lái)，人們?cè)噲D借助數(shù)學(xué)工具來(lái)研究發(fā)酵過(guò)程的振蕩機(jī)制，自20世紀(jì)70年代以來(lái)已有多個(gè)針對(duì)酵母振蕩過(guò)程的數(shù)學(xué)模型見(jiàn)諸報(bào)道[14-16]。為降低模型復(fù)雜度，本文采用Wolf等[14]給出的葡萄糖經(jīng)酵母代謝生成乙醇的途徑，描述了與實(shí)驗(yàn)觀測(cè)一致的酵母糖酵解途徑的最簡(jiǎn)核心模型，途徑中包含乙酸在細(xì)胞間交換的過(guò)程。如圖1所示。

為了便于研究，圖1給出的代謝過(guò)程是一種簡(jiǎn)化途徑，其中虛線框代表細(xì)胞膜，0和ex分別代表底物葡萄糖和代謝物乙酸在細(xì)胞間交換的流量。

Wolf等[14]依據(jù)圖1的代謝過(guò)程所建立的數(shù)學(xué)模型見(jiàn)式（2），其中相關(guān)參數(shù)值取自文獻(xiàn)[14]。這里，為了保證后續(xù)求解模型過(guò)程的順利，將NADH ? NAD、ATP ? ADP這樣的“二元環(huán)”方程用代數(shù)方程改寫(xiě)并求解。

為了數(shù)學(xué)表達(dá)式描述的方便，令x=[Glc, F16p, G3P, GP, Pyr, Acd, ExAcd, ATP, NAD]。

其中

2.2 仿真分析

通過(guò)動(dòng)態(tài)仿真分析，可初步了解系統(tǒng)的動(dòng)態(tài)行為，是接下來(lái)的參數(shù)分析的切入點(diǎn)。應(yīng)用前期工作中開(kāi)發(fā)的集成優(yōu)化工具GIEPT[17]，利用動(dòng)態(tài)優(yōu)化的方法進(jìn)行求解，方法見(jiàn)文獻(xiàn)[18-20]。

2.3 參數(shù)分支分析

穩(wěn)定的系統(tǒng)隨著一些參數(shù)連續(xù)變化，經(jīng)過(guò)某臨界值時(shí)，狀態(tài)會(huì)發(fā)生改變：一種情況是出現(xiàn)多穩(wěn)態(tài)現(xiàn)象；另一種情況是由穩(wěn)態(tài)變?yōu)檎袷?。這種由參數(shù)的變化導(dǎo)致的系統(tǒng)狀態(tài)躍遷的現(xiàn)象被稱為分支[21-22]。通過(guò)分支分析，能夠分析非線性系統(tǒng)中參數(shù)對(duì)狀態(tài)穩(wěn)定性的影響。

首先只考慮一個(gè)二階自治系統(tǒng)[23]，如式（3）所示。系統(tǒng)的平衡點(diǎn)可以由系統(tǒng)方程為零求得，在平衡點(diǎn)處，系統(tǒng)的動(dòng)態(tài)特性會(huì)有不同的表現(xiàn)。

根據(jù)系統(tǒng)連續(xù)性理論可知，系統(tǒng)平衡點(diǎn)處的Jacobi矩陣為式（4），可有助于分析其特性。

(4)

系統(tǒng)的特征值和特征向量由式（5）給出。

=v(5)

一般來(lái)說(shuō)，判斷特征值的屬性，可確定平衡點(diǎn)的特性。除此之外，還可通過(guò)的跡式（5）和行列式值式（6）來(lái)確定相關(guān)性質(zhì)。

(7)

當(dāng)det()>0，且tr()<0時(shí)，平衡點(diǎn)是穩(wěn)定的；當(dāng)det()>0，且tr()>0時(shí)，平衡點(diǎn)是不穩(wěn)定的。當(dāng)特征值由全負(fù)實(shí)部變?yōu)榱銓?shí)部的情形時(shí)，通常是系統(tǒng)穩(wěn)定性發(fā)生變化的邊界點(diǎn)，在穩(wěn)態(tài)譜圖中即表現(xiàn)為分支現(xiàn)象。若系統(tǒng)具有一對(duì)純虛數(shù)特征值且其他特征值具有負(fù)實(shí)部，或det()=0，且tr(=0，則為Hopf分支，這一般表示系統(tǒng)由穩(wěn)定平衡點(diǎn)過(guò)度到極限環(huán)振蕩的臨界點(diǎn)[21]。

在二階系統(tǒng)基礎(chǔ)上，當(dāng)系統(tǒng)模型的維度較高、同時(shí)包含較多的參數(shù)時(shí)，分析過(guò)程計(jì)算較為復(fù)雜[24]，分支點(diǎn)的數(shù)量和類型也更多，甚至?xí)a(chǎn)生混沌現(xiàn)象[25]。

3 結(jié)果與分析

生化系統(tǒng)因其自身的復(fù)雜性蘊(yùn)含著大量的非線性特性，振蕩是非線性動(dòng)力系統(tǒng)和非線性理論研究的重要內(nèi)容，對(duì)微生物培養(yǎng)過(guò)程中出現(xiàn)的振蕩現(xiàn)象進(jìn)行研究有助于加深對(duì)生命活動(dòng)規(guī)律的進(jìn)一步認(rèn)識(shí)。

3.1 動(dòng)態(tài)過(guò)程仿真

使用GIEPT調(diào)用非線性規(guī)劃解題器進(jìn)行求解，結(jié)果如2所示。

圖中上段曲線為1,6-二磷酸果糖F16p的曲線，下段為葡萄糖Glu的變化情況（此處略去其他代謝物的計(jì)算結(jié)果）。需要說(shuō)明的是，為了使圖示效果更清晰可辨，圖2中的數(shù)據(jù)點(diǎn)每間隔20個(gè)繪制1次。

從圖2的動(dòng)態(tài)仿真結(jié)果可以看出，該系統(tǒng)具有明顯的等幅振蕩特性，除圖2中的兩種代謝物外，系統(tǒng)式（2）中其他代謝物均產(chǎn)生振幅不同的振蕩現(xiàn)象。

對(duì)于此類周期等幅振蕩現(xiàn)象，一般亦稱為極限環(huán)振蕩?？疾齑x物三磷酸腺苷（ATP）與輔酶I（NAD）的濃度變化關(guān)系，在二維平面中繪制兩者關(guān)系，結(jié)果如圖3所示。

結(jié)果顯示，糖酵解系統(tǒng)出現(xiàn)典型的極限環(huán)振蕩現(xiàn)象。由于NAD為糖酵解過(guò)程提供氧化還原當(dāng)量，而ATP則為代謝過(guò)程提供能量，因此這兩種物質(zhì)的變化情況與生物過(guò)程有著緊密的聯(lián)系。接下來(lái)著重分析NAD與ATP的模型參數(shù)與極限環(huán)振蕩的關(guān)系。

3.2 分支分析

首先對(duì)ATP降解為ADP的速率常數(shù)15進(jìn)行參數(shù)穩(wěn)態(tài)分析，使15參數(shù)值在15～45 min-1范圍內(nèi)變化同時(shí)計(jì)算系統(tǒng)特征值，在15分別為24.89 min-1及43.78 min-1時(shí)出現(xiàn)兩個(gè)系統(tǒng)特征值變化的臨界點(diǎn)，判斷為Hopf分支點(diǎn)。其中在15=24.89 min-1處，Hopf分支點(diǎn)的第一李雅普諾夫系數(shù)為正值，意味著此分支點(diǎn)為次臨界，同時(shí)此處存在一個(gè)不穩(wěn)定的極限環(huán)；而在15=43.78 min-1處分支點(diǎn)的第一李雅普諾夫系數(shù)為負(fù)值，這表示該分支點(diǎn)為超臨界且存在的極限環(huán)是穩(wěn)定的。結(jié)果如圖4所示。

圖4中黑色實(shí)線為15變化的系統(tǒng)穩(wěn)態(tài)軌跡。15在兩處Hopf分支點(diǎn)間取值時(shí)系統(tǒng)出現(xiàn)極限環(huán)振蕩。圖中藍(lán)色區(qū)域?yàn)橄到y(tǒng)極限環(huán)隨著參數(shù)變化，其周期、軌跡形態(tài)的變化過(guò)程，特別是在15=23.61min-1處，環(huán)軌跡出現(xiàn)極限分支點(diǎn)，這表明環(huán)軌跡在此處出現(xiàn)折疊，在圖4中表現(xiàn)為自次臨界Hopf分支點(diǎn)（左）開(kāi)始產(chǎn)生極限環(huán)，在超臨界Hopf分支點(diǎn)（右）處結(jié)束。

由于系統(tǒng)涉及到乙酸在細(xì)胞間的交換過(guò)程，因此考慮乙酸降解速率常數(shù)16對(duì)NAD與ATP濃度波動(dòng)的影響。結(jié)果如圖5所示。使16在大于10 min-1的范圍變化，同時(shí)計(jì)算系統(tǒng)特征值，16在27.03 min-1處出現(xiàn)一個(gè)次臨界Hopf分支點(diǎn)（第一李雅普諾夫系數(shù)為正值），在16=10.93 min-1處出現(xiàn)環(huán)軌跡極限分支點(diǎn)。與15不同，16在經(jīng)過(guò)27.03 min-1處的Hopf分支點(diǎn)后，未檢測(cè)到新的分支點(diǎn)，而在環(huán)軌跡上表現(xiàn)為振幅變化不大、周期逐步放大的不穩(wěn)定狀態(tài)。

在圖5中，乙酸降解速率16決定了系統(tǒng)能否產(chǎn)生極限環(huán)振蕩。只要16超過(guò)27.03 min-1，系統(tǒng)就將產(chǎn)生振蕩，并且在這之后系統(tǒng)振蕩的振幅、周期與參數(shù)16的增長(zhǎng)呈弱相關(guān)性。由此可以推測(cè)，乙酸的降解量可能使得胞內(nèi)代謝物產(chǎn)生振蕩，這為通過(guò)產(chǎn)物控制代謝過(guò)程的振蕩提供了依據(jù)。

限于篇幅，其他參數(shù)的分支分析結(jié)果并未列出。通過(guò)對(duì)各個(gè)參數(shù)的分支分析，可以獲知能夠影響系統(tǒng)振蕩產(chǎn)生的參數(shù)取值范圍。同時(shí)，底物流加速率也是極限環(huán)振蕩能否產(chǎn)生的影響因素。

3.3 驗(yàn)證

經(jīng)過(guò)對(duì)參數(shù)的分支分析，能夠獲得系統(tǒng)振蕩的產(chǎn)生條件，通過(guò)調(diào)節(jié)參數(shù)，可獲得穩(wěn)定的代謝過(guò)程。然而，生物代謝過(guò)程中的振蕩現(xiàn)象分為兩類：自發(fā)振蕩和強(qiáng)制振蕩[26]。自發(fā)振蕩無(wú)須外部周期性變化因素來(lái)激勵(lì)，其維持動(dòng)力源于發(fā)酵體系內(nèi)部矛盾相互轉(zhuǎn)化；而強(qiáng)制振蕩一般受外部條件維持，如底物的消耗、產(chǎn)物的積累。強(qiáng)制振蕩較容易調(diào)控[27-28]。

在該模型的16個(gè)參數(shù)中，大多數(shù)參數(shù)為胞內(nèi)反應(yīng)速率常數(shù)，而16是胞外乙酸的降解速率常數(shù)，該參數(shù)值容易通過(guò)調(diào)節(jié)溫度等操作來(lái)改變。接下來(lái)，在分支分析的基礎(chǔ)上，調(diào)節(jié)16參數(shù)值，驗(yàn)證系統(tǒng)的穩(wěn)定性。由于16只有一個(gè)Hopf分支點(diǎn)，取其鄰近值，從仿真結(jié)果上可以看出，系統(tǒng)趨于穩(wěn)定（如圖6所示，這里只列出ATP、NAD的變化情況）；當(dāng)進(jìn)一步減小16時(shí)，結(jié)果會(huì)出現(xiàn)發(fā)散的現(xiàn)象（圖略）。

有觀點(diǎn)[29]指出，產(chǎn)物乙醇的反饋抑制作用，是產(chǎn)生振蕩的一個(gè)內(nèi)因。通過(guò)以上分析，乙酸的降解速率很有可能是產(chǎn)生振蕩的另一因素。

綜上，借助數(shù)學(xué)工具對(duì)模型進(jìn)行參數(shù)分支分析，一方面，在連續(xù)培養(yǎng)過(guò)程中，可以利用這些信息為降低振蕩的不利影響提供指導(dǎo)；另一方面，可以利用振蕩產(chǎn)生后有利的影響，如利用產(chǎn)物濃度周期振蕩的特性，可以設(shè)計(jì)優(yōu)化的產(chǎn)品回收及分離工藝。

4 結(jié) 論

以酵母發(fā)酵生產(chǎn)乙醇為代表的生物過(guò)程是復(fù)雜并高度集成化的系統(tǒng)，難以借助傳統(tǒng)的實(shí)驗(yàn)手段高效準(zhǔn)確地加以分析。對(duì)這類復(fù)雜的生物系統(tǒng)過(guò)程的研究，客觀上需要系統(tǒng)生物學(xué)和生物信息學(xué)提供方法上的支持。

針對(duì)典型的釀酒酵母發(fā)酵生成乙醇的振蕩過(guò)程，對(duì)其數(shù)學(xué)模型進(jìn)行了深入分析，獲得有意義的結(jié)果如下。

（1）對(duì)過(guò)程模型進(jìn)行仿真，并通過(guò)狀態(tài)關(guān)系分析確定該過(guò)程能夠產(chǎn)生極限環(huán)振蕩。

（2）對(duì)模型參數(shù)進(jìn)行分支分析，其中15具有一個(gè)超臨界Hopf分支點(diǎn)和一個(gè)次臨界Hopf分支點(diǎn)，這是能夠引起極限環(huán)振蕩的范圍。

（3）根據(jù)分支點(diǎn)臨界值的分析結(jié)果，調(diào)整16參數(shù)，可使系統(tǒng)趨于穩(wěn)定，提出乙酸降解速率可能是引起振蕩行為的一個(gè)重要因素。

從數(shù)學(xué)模型參數(shù)出發(fā)，考察其如何影響振蕩的產(chǎn)生，是后續(xù)研究振蕩機(jī)制的基礎(chǔ)。

References

[1] HEINZLE E, BIWER A P, CHARLES L C. Development of Sustainable Bioprocesses: Modeling and Assessment[M]. Sussex: Wiley, 2007.

[2] KARSENTI E. Self-organization in cell biology: a brief history[J]. Nat. Rev. Mol. Cell Biol., 2008, 9(3): 255-262.

[3] KRUSE K, JULICHER F. Oscillations in cell biology[J]. Curr. Opin. Cell Biol., 2005, 17(1): 20-26.

[4] LONGO D M, SELIMKHANOV J, KEARNS J D,. Dual delayed feedback provides sensitivity and robustness to the NF-κB signaling module [J]. PLoS Comput. Biol., 2013, 9(6): e1003112.

[5] WILLIAMSON T, ADIAMAH D, SCHWARTZ J M,. Exploring the genetic control of glycolytic oscillations in[J]. BMC Syst. Biol., 2012, 6: 108.

[6] GUSTAVSSON A K, VAN NIEKERK D D, ADIELS C B,. Heterogeneity of glycolytic oscillatory behaviour in individual yeast cells[J]. FEBS Lett., 2014, 588(1): 3-7.

[7] DU PREEZ F B, VAN NIEKERK D D, KOOI B,. From steady-state to synchronized yeast glycolytic oscillations (Ⅰ): Model construction[J]. FEBS J., 2012, 279(16): 2810-2822.

[8] GUSTAVSSON A K, VAN NIEKERK D D, ADIELS C B,. Sustained glycolytic oscillations in individual isolated yeast cells[J]. FEBS J., 2012, 279(16): 2837-2847.

[9] GOLDBETER A, GERARD C, GONZE D,. Systems biology of cellular rhythms[J]. FEBS Lett., 2012, 586(18): 2955-2965.

[10] YAMAZAKI S, MIKI K, KANO K,. Mechanistic study on the role of the NAD+-NADH ratio in the glycolytic oscillation with a pyruvate sensor[J]. Journal of Electroanalytical Chemistry, 2001, 516(1/2): 59-65.

[11] DUYSENS L N, AMESZ J. Fluorescence spectrophotometry of reduced phosphopyridine nucleotide in intact cells in the near-ultraviolet and visible region[J]. Biochim. Biophys. Acta, 1957, 24(1): 19-26.

[12] CHANCE B, ESTABROOK R W, GHOSH A. Damped sinusoidal oscillations of cytoplasmic reduced pyridine nucleotide in yeast cells[J]. Proc. Natl. Acad. Sci. USA, 1964, 51: 1244-1251.

[13] RICHARD P. The rhythm of yeast[J]. Fems. Microbiol. Rev., 2003, 27(4): 547-557.

[14] WOLF J, PASSARGE J, SOMSEN O J,. Transduction of intracellular and intercellular dynamics in yeast glycolytic oscillations[J]. Biophys. J., 2000, 78(3): 1145-1153.

[15] WOLF J, HEINRICH R. Effect of cellular interaction on glycolytic oscillations in yeast: a theoretical investigation[J]. Biochem. J., 2000, 345: 321-334.

[16] DU PREEZ F B, VAN NIEKERK D D, SNOEP J L. From steady-state to synchronized yeast glycolytic oscillations(Ⅱ): Model validation[J]. FEBS J., 2012, 279(16): 2823-2836.

[17] 潘多濤, 黃明忠, 張學(xué)軍, 等. 工程規(guī)劃建模研究及流程優(yōu)化的實(shí)現(xiàn)[J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2012, 38(10): 1486-1490. PAN D T, HUANG M Z, ZHANG X J,. Engineering modeling programming and implementation of process optimization[J]. Journal of Beijing University of Technology, 2012, 38(10): 1486-1490.

[18] SHIVAKUMAR K, BIEGLER L T. Simultaneous dynamic optimization strategies: recent advances and challenges [J]. Computers and Chemical Engineering, 2006, 30(10/11/12): 1560-1575.

[19] 潘多濤, 史洪巖, 黃明忠, 等. 復(fù)雜生物過(guò)程的仿真分析[J]. 系統(tǒng)仿真學(xué)報(bào), 2014, 26(3): 670-674. PAN D T, SHI H Y, HUANG M Z,. Simulation and analysis of complex biological processes [J]. Journal of System Simulation, 2014, 26(3): 670-674.

[20] 潘多濤, 史洪巖, 黃明忠, 等. 非線性規(guī)劃在生物代謝仿真過(guò)程中的應(yīng)用[J]. 控制工程, 2014, (6): 896-899. PAN D T, SHI H Y, HUANG M Z,. Application of nonlinear programming for simulation of biological metabolic process [J]. Control Engineering of China, 2014, (6): 896-899.

[21] CRAWFORD J D. Introduction to bifurcation theory[J]. Reviews of Modern Physics, 1991, 63(4): 991-1037.

[22] KUZNETSOV Y A. Elements of Applied Bifurcation Theory[M]. New York: Springer-Verlag, 2004.

[23] GOVAERTS W. Numerical bifurcation analysis for ODEs[J]. Journal of Computational and Applied Mathematics, 2000, 125(1): 57-68.

[24] DHOOGE A, GOVAERTS W, KUZNETSOV Y A. MATCONT: a MATLAB package for numerical bifurcation analysis of ODEs[J]. ACM Transactions on Mathematical Software, 2003, 29(2): 141-164.

[25] DENG B. Food chain chaos due to junction-fold point[J]. Chaos, 2001, 11(3): 514-525.

[26] 申渝. 酵母細(xì)胞超高濃度乙醇連續(xù)發(fā)酵振蕩行為的研究[D]. 大連: 大連理工大學(xué), 2009. SHEN Y. Exploration for oscillation in continuous VHG ethanol fermentation with[D]. Dalian: Dalian University of Technology, 2009.

[27] 申渝, 葛旭萌, 李寧, 等. 高濃度乙醇連續(xù)發(fā)酵振蕩過(guò)程中代謝通量分析及誘發(fā)機(jī)理[J]. 化工學(xué)報(bào), 2009, 60(6): 1519-1528. SHEN Y, GE X M, LI N,. Metabolic flux analysis and mechanistic study of process oscillation in continuous VHG ethanol fermentation with[J]. CIESC Journal, 2009, 60(6): 1519-1528.

[28] 楊蕾, 陳麗杰, 白鳳武. 高濃度酒精連續(xù)發(fā)酵過(guò)程中振蕩行為的模擬及填料弱化振蕩的機(jī)理[J]. 化工學(xué)報(bào), 2007, 58(3): 715-721. YANG L, CHEN L J, BAI F W. Dynamic models of VHG continuous ethanol fermentation and mechanisms of oscillation attenuation by packing[J]. Journal of Chemical Industry and Engineering(China), 2007, 58(3): 715-721.

[29] PORRO D, MARTEGANI E, RANZI B M,. Oscillations in continuous cultures of budding yeast: a segregated parameter analysis[J]. Biotechnol. Bioeng., 1988, 32(4): 411-417.

Analysis of metabolic oscillation processes in

PAN Duotao1,2, SHI Hongyan2, YUAN Decheng2,XIU Zhilong1

(1School of Life Science and Biotechnology, Dalian University of Technology, Dalian 116024, Liaoning, China;2Chemical Control Technology Key Laboratory of Liaoning Province, Shenyang University of Chemical Technology, Shenyang 110142, Liaoning, China)

Oscillation phenomenon is an inherent characteristic in biological systems and plays an important role in many dynamic bioprocesses. In order to explore the certain conditons that could possibly boost a oscillation, the metabolic pathway of the, glycolysis were researched, and the parameters of mathmatical model was analyzed. Firstly, the simulation results associated the phase diagrams indicated that the model exists sustained oscillations with constant amplitude (limit cycle oscillations). Next, bifurcation analysis approach was used to investigate the infulence of parameters for the system producing oscillations. The results showed that several parameters of the model could lead to oscillations and the range of parameters’ value was obtained, which could be applied to direct the manipulation of metabolic oscillations.

systems biology; glycolysis; oscillation; limit cycle; bifurcation analysis; nonlinear programming

10.11949/j.issn.0438-1157.20161628

TQ 920. 1

A

0438—1157（2017）03—0964—06

國(guó)家自然科學(xué)基金項(xiàng)目（21476042）；遼寧省教育廳科研一般項(xiàng)目（L2014168）；遼寧省博士科研啟動(dòng)基金項(xiàng)目（201501072）。

2016-11-16收到初稿，2016-11-26收到修改稿。

聯(lián)系人：修志龍。第一作者：潘多濤（1979—），男，博士，講師。

2016-11-16.

Prof.XIU Zhilong, Zhlxiu@dlut.edu.cn

supported by the National Natural Science Foundation of China (21476042), the Educational Commission of Liaoning Province (L2014168) and the Doctoral Research Fund of Liaoning Province(201501072).