枯草芽孢桿菌GXA-28發酵生產γ-聚谷氨酸動力學研究

王青龍，陳桂光，曾偉，鄭雙鳳，唐真，粱智群

（廣西大學生命科學與技術學院，廣西南寧530004）

枯草芽孢桿菌GXA-28發酵生產γ-聚谷氨酸動力學研究

王青龍，陳桂光，曾偉，鄭雙鳳，唐真，粱智群*

（廣西大學生命科學與技術學院，廣西南寧530004）

對耐熱高產菌株枯草芽孢桿菌GXA-28分批發酵生產γ-聚谷氨酸的動力學特征進行了研究，基于菌體生長特性，結合Logistic方程和Luedeking-Piret方程，提出了菌體生長、產物合成、葡萄糖消耗以及谷氨酸鈉消耗的動力學模型。應用Origin8.5對數據進行分析處理，得到了Bacillis subtilisGXA-28分批發酵合成γ-聚谷氨酸相應的動力學參數。將模型預測值和實驗值進行比較，結果表明，模型基本反映了枯草芽孢桿菌GXA-28分批發酵過程中的動力學特征。

枯草芽孢桿菌GXA-28；γ-聚谷氨酸；動力學模型

γ-聚谷氨酸（γ-polyglutamic acid，γ-PGA）主要是自然界中微生物來源的一種強水溶性的多聚氨基酸高分子化合物[1]。其具有良好的生物可降解性、可食用性、生物相容性和保濕性[2]，對人體和環境無毒無公害，可廣泛用于工業、農業、醫藥、食品、環境保護等[3]，還可應用于化妝品、食品用水凝膠[4]、黏稠劑[5]、土壤沙地的蓄水保水劑[6]、高吸水劑[7-8]、高強度纖維的制備以及生物醫學等領域[9]。近年來，隨著國內外對γ-聚谷氨酸研究的深入，γ-聚谷氨酸作為一種可降解高分子化合物，越來越顯現出其廣闊的科研價值和應用前景。目前由于受到發酵工藝、生產成本及生產菌株的限制，國內尚未見大規模生產的相關報道，因此對聚谷氨酸發酵過程的研究就有了很重要的現實意義。

國內目前主要是將菌種選育以及對發酵條件的優化作為研究重點，有關發酵動力學模型的報道相對較少[10]。本實驗考察了γ-PGA分批發酵過程中的菌體生長、產物合成、底物消耗情況。運用Origin 8.5數據處理軟件對枯草芽孢桿菌GXA-28發酵動力學模型進行非線性擬合，并獲得了相應的動力學模型及模型參數。利用擬合方程可以預測分批發酵過程中菌體量、γ-PGA量、葡萄糖量、谷氨酸鈉量的變化情況，為下一步的放大培養及補料分批發酵等提供了基礎數據。

1 材料與方法

1.1 材料與試劑

枯草芽孢桿菌GXA-28（B.subtilisGXA-28）：生命科學與技術學院自主篩選，菌株保藏編號CCTCC M2012347。

葡萄糖：天津市科密歐化學試劑有限公司；谷氨酸鈉（純度≥99%）：南寧市荷花味精有限公司；酵母膏：北京奧博星生物技術有限公司；KH2PO4：天津市紅巖化學試劑廠；MgSO4：天津博迪化工股份有限公司。

種子液培養基：葡萄糖10g/L、酵母膏5g/L、谷氨酸鈉5g/L、KH2PO40.5g/L、MgSO40.1g/L，pH 6.5±0.1。

發酵培養基：葡萄糖30g/L、酵母膏2.5g/L、谷氨酸鈉20g/L、KH2PO40.5g/L、MgSO40.1g/L，pH 6.5±0.1。

1.2 儀器與設備

YP1200型電子天平：上海精科實業有限公司；HVE-50自動滅菌鍋：日本HIRAYAMA公司；SPX-250型生化培養箱：上海躍進醫療儀器廠；EBA21型高速離心機：德祥科技有限公司；SBA-40D生物傳感分析儀：山東省科學院生物研究所；UV-1800紫外分光光度計：日本SHIMADIU公司；SKY-211B恒溫培養振蕩器：上海蘇坤實業有限公司。

1.3 方法

1.3.1 種子培養液的制備

取一管甘油保藏菌種接入裝有30mL種子培養基的250mL三角瓶中，于42℃、160r/min條件下培養16h，得到種子培養液。

1.3.2 分批發酵培養

將培養得到的種子培養液取1mL接種于裝有50mL發酵培養基的三角瓶中，于42℃、160r/min條件下恒溫振蕩培養。考察不同時間段培養液中的菌體生長、產物合成、底物消耗情況。

1.3.3 殘糖、殘谷氨酸鈉含量測定

由上述培養條件獲得的發酵液采用生物傳感分析儀測定性能指標，溫度10～35℃，相對濕度＜90%，可同時測定β-D-葡萄糖、L-谷氨酸含量，測定范圍：0～100mg/100mL。

1.3.4 聚γ-谷氨酸產量的測定

按照參考文獻[11]的方法測定。

1.3.5 枯草芽孢桿菌GXA-28生物量的測定

發酵液稀釋10倍后，用可見分光光度計于波長660nm處測定其吸光度值。

2 結果與討論

2.1 動力學模型的建立

2.1.1 菌體生長動力學模型的建立

描述菌體生長的動力學模型最常用的有Monod方程和Logistic方程，Monod方程主要用來描述理想化的均衡生長動力學模型。Logistic方程[12]與Monod方程相比，其限制因子主要是菌體濃度或生長空間，而γ-PGA分批發酵過程后期發酵液黏度較大，因此菌體生長將受到生長空間及氧傳質等因素的限制，選用Monod方程進行建模會導致結果有較大偏差，因此本實驗選用Logistic方程對菌體生長規律進行描述。Logistic方程為：

式中：X為菌體生物量（用OD660nm值表示）；μm為最大比生長速率，h-1；Xmax為菌體最大生物量（用OD660nm值表示）。

將上式對時間t積分，得到Logistic方程的積分形式

式中：X0為菌體初始生物量（用OD660nm值表示）；μm為最大比生長速率，h-1；Xmax為菌體最大生物量（用OD660nm值表示）；t為時間，h。

2.1.2 產物合成動力學模型建立

通過跟蹤取樣分析發現，枯草芽孢桿菌GXA-28分批發酵生產γ-PGA分為2個階段：第一個階段主要是菌體快速的生長，產物合成量較少，此時基質的消耗主要用于菌體的生長；第二階段時γ-PGA開始快速合成，在此過程中菌體量已經接近穩定值，此時基質的消耗主要用于γ-PGA的合成。由此可以明顯看出，微生物發酵法生產γ-PGA屬于部分偶聯型[13]。因此本實驗選用由Luedeking和Piret提出的Luedeking-Piret方程對γ-PGA的合成規律進行描述如下：

將上式簡化得

式中：α，β為動力學模型參數；X為菌體生物量（用OD660nm值表示）；μm為最大比生長速率，h-1；Xmax為菌體最大生物量（用OD660nm值表示）；dP/dt為γ-聚谷氨酸產物合成速率，g（/L·h）。

將上式積分得

式中：P為γ-聚谷氨酸產量，g/L；α，β為動力學模型參數；X0為菌體初始生物量（用OD660nm值表示）；μm為最大比生長速率，h-1；Xmax為菌體最大生物量（用OD660nm值表示）；t為時間，h。

2.1.3 葡萄糖消耗動力學模型的建立

葡萄糖在聚谷氨酸發酵過程中是細胞生長和產物合成所需的重要營養成分，其消耗主要用于3個方面：①生成新的細胞；②用來維持細胞基本生命活動；③是γ-PGA的合成。因此葡萄糖消耗動力學可以用如下方程描述

式中：YX/S為葡萄糖用于菌體生長的得率常數，g/g；YP/S為葡萄糖用于產物積累的得率常數，g/g；mx為微生物的維持常數，g/（L·h）；X為菌體生物量（用OD660nm值表示）；dS1/dt為葡萄糖消耗速率，g/（L·h）。

在分批發酵過程中，菌體呼吸消耗的基質可以認為是菌體生長消耗的一部分[14-15]，因此葡萄糖的消耗就可以簡單的分為菌體生長的消耗和產物合成的消耗兩部分。故葡萄糖消耗的動力學模型可以簡化為

式中：k1為葡萄糖用于菌體生長的得率常數，g/g；k2為葡萄糖用于產物積累的得率常數，g/g。將（1）、（3）式子帶入得

式中：k1為葡萄糖用于菌體生長的得率常數，g/g；k2為葡萄糖

用于產物積累的得率常數，g/g；X為菌體生物量（用OD660nm值表示）；Xmax為菌體最大生物量（用OD660nm值表示）；α，β為動力學模型參數；μm為最大比生長速率，h-1。

將（8）式子積分得

式中：S1為葡萄糖的質量濃度，g/L；S0為葡萄糖初始質量濃度，g/L；X0為菌體初始生物量（用OD660nm值表示）；k1為葡萄糖用于菌體生長的得率常數，g/g；k2為葡萄糖用于產物積累的得率常數，g/g；Xmax為菌體最大生物量（用OD660nm值表示）；α，β為動力學模型參數；μm為最大比生長速率，h-1。

2.1.4 谷氨酸鈉消耗動力學模型的建立

枯草芽孢桿菌GXA-28屬于底物依賴型菌株，因此谷氨酸納的消耗主要用于γ-PGA的合成，谷氨酸鈉消耗的動力學模型可以簡化為：

式中：dS2/dt為谷氨酸鈉消耗速率，g/（L·h）；k3為谷氨酸鈉用于產物合成的得率常數，g/g。

對上式積分得

式中：S2為谷氨酸鈉質量濃度，g/L；S0為葡萄糖初始質量濃度，g/L；X0為菌體初始生物量（用OD660nm值表示）；k3為谷氨酸鈉用于產物合成的得率常數，g/g；Xmax為菌體最大生物量（用OD660nm值表示）；α為動力學模型參數；μm為最大比生長速率，h-1。

2.2 動力學模型參數求解

根據實驗數據和上述的模型方程應用Origin 8.5軟件對實驗值進行非線性擬合，獲得最優化參數估計值。并得出實驗值與預測值之間的方差，分別為菌體生長動力學模型R2=0.997，產物合成動力學模型R2=0.993，葡萄糖消耗動力學模型R2=0.991，底物谷氨酸鈉模型R2=0.989。參數估計值結果見表1。

表1 分批發酵生產γ-PGA的動力學參數估計值Table 1 Estimated kinetic parameter of γ-PGA batch fermentation

將擬合所得模型參數代入方程可以得出γ-PGA分批發酵的動力學模型分別為：

菌體生長的動力學模型

產物γ-PGA合成的動力學模型

葡萄糖消耗的動力學模型

谷氨酸鈉消耗的動力學模型

2.3 曲線擬合分析

為了驗證上述動力學模型的可靠性，本實驗運用Origin 8.5數據處理軟件對γ-PGA分批發酵過程進行非線性擬合分析，得到了γ-PGA分批發酵過程中菌體生長、產物合成、葡萄糖消耗以及谷氨酸鈉消耗的擬合曲線分別如圖1～圖4所示。從圖1～圖4可以看出，除了在發酵開始階段和結束階段的個別點的相對誤差較大外，在發酵中期，所建動力學模型都能對真實的發酵情況進行很好的擬合。可以認為擬合曲線基本反映了γ-PGA真實的發酵過程。

圖1 菌體生物量實驗值與擬合曲線的比較Fig.1 Comparison of fitted curve and measured biomass

圖2 γ-PGA產量實驗值與擬合曲線的比較Fig.2 Comparison of fitted curve and measured γ-PGA yield

圖3 剩余葡萄糖量實驗值與擬合曲線比較Fig.3 Comparison of fitted curve and measured remaining glucose content

圖4 谷氨酸鈉剩余量實驗值與擬合曲線比較Fig.4 Comparison of fitted curve and measured sodium glutamate content

3 結論

通過對枯草芽孢桿菌GXA-28發酵合成γ-PGA特征分析，運用Origin 8.5數據處理軟件進行非線性擬合，得出其分批發酵過程中菌體生長、γ-PGA合成、葡萄糖消耗、底物消耗的動力學模型及相關參數，通過擬合曲線圖可以看出，實驗值與模型曲線非常接近。將預測值與實驗值進行對比可以看出，除發酵初期和發酵末期的個別點相對誤差較大外，在發酵中期動力學模型能很好的對實驗數據進行擬合。因此可以用該模型來預測發酵液中菌體量、產物合成量、葡萄糖以及谷氨酸鈉在發酵過程的變化情況，為以后工業化生產提供理論依據及補料分批發酵提供參考。

[1]趙紫華，馬霞，劉蕊，等.高產γ-聚谷氨酸菌株的選育與鑒定[J].中國釀造，2007，26（7）：32-34.

[2]SHIH I L,VAN Y T.The production of poly-(γ-glutamic acid)from microorganisms and its various applications[J].Bioresource Technol, 2001,79(3):207-225.

[3]ASHIUCHI M.Occurrence and biosynthetic mechanism of poly-gamma-glutamic acid[M]//HAMANO Y.Amino-acid homopolymers occurring in nature.Berlin,Springer Berlin Heidelberg,2010:77-93.

[4]KONNO A,TAGUCHI T,YAMAGUCHI T.Bakery products an noodles containing poly glutamic acid:US,4888193[P].1989-11-19.

[5]YAHATA K,SADANOBU J,ENDO T.Preparation of poly-benz-γ-glutamic acid fiber[J].Polym Prepr Jpn,1992,42:1077-1082.

[6]CHOI H J,KUNIOKA M.Preparation conditions and swelling equilibria of hydrogel prepared byγ-irradiation from microbial poly(γ-glutamic acid)[J].Radiat Phys Chem,1995,46(2):175-179.

[7]BORBELY M,NAGASAKI Y,BORBELY J,et al.Biosynthesis and chemical modification of poly(γ-glutamic acid)[J].Polym Bull,1994, 32(2):127-132.

[8]KUNIOKA M,FURUSAWA K.Poly(γ-glutamic acid)hydrog prepared from microbial poly(γ-glutamic acid)acid alkaced amine with watersoluble carnodiimide[J].J Appl Polym Sci,1997,65(10):1889-1896.

[9]YASUYOSHI A,YUJI F,SHUICHI K.Health food with poly(γ-glutamic acid)as the chief ingredient:Japan,5095767A2[P].1993-04-20.

[10]任尚美，王曉麗，馬霞，等.枯草芽孢桿菌發酵生產γ-聚谷氨酸的動力學研究[J].中國釀造，2008，27（21）：59-62.

[11]ZENG W,CHEN G G,ZHANG Y K,et al.Studies on the UV spectrum of poly（γ-glutamic acid）based on development of a simple quantitative method[J].Int J Biol Macromol,2012,51(1):83-90.

[12]王濟川，郭志剛.Logistic回歸模型—方法與應用[M].北京：高等教育出版社，2001.

[13]聶光宇.枯草芽孢桿菌生產γ-聚谷氨酸的培養條件及動力學研究[D].杭州：浙江大學碩士論文，2006.

[14]STEFAN S,KHOLODENKO B N,WESTERHOFF H V.Cellular information transfer regarded from a stoichiometry and control analysis perspective[J].Biosystems,2000,55(1-3):73-78.

[15]GIUSEPPIN M L,VAN RIEL N A.Metabolic modelin ofSaccharomyces cerevisiaeusing the optimal control of homestasis:a cybernetic model definition[J].Metab Eng,2000,2(1):14-21.

Fermentation kinetics model of γ-polyglutamic acid production byBacillis subtilisGXA-28

WANG Qinglong,CHEN Guiguang,ZENG Wei,ZHENG Shuangfeng,TANG Zhen,LIANG Zhiqun*
(College of Life Science and Technology,Guangxi University,Nanning 530004,China)

The kinetics feature of a heat-resistant and high-yieldγ-polyglutamic acid strainBacillis subtilisGXA-28 was investigated.Based on bacterial growth characteristics,combining with Logistic and Luedeking-Piret equation,the kinetic model of microbial growth,γ-polyglutamic acid formation and glucose consumption during fermentation were proposed.The data was analyzed with origin 8.5 to generate the kinetics parameters ofγ-ploy glutamic acid production byB.subtilisGXA-28 batch fermentation.The prediction value was compared with the experimental value;it showed that the model could provide reasonable description for the batch fermentation process ofB.subtilisGXA-28.

Bacillis subtilisGXA-28;γ-polyglutamic acid;kinetics model

Q936

A

0254-5071（2014）04-0075-04

10.3969/j.issn.0254-5071.2014.04.018

2014-03-13

廣西研究生教育創新計劃資助項目（YCBZ2012004）；2013年廣西大學亞熱帶農業生物資源保護與利用國家重點實驗室研究課題；2013年廣西大學微生物及植物遺傳工程教育部重點實驗室主任基金課題

王青龍（1988-），男，碩士研究生，研究方向為發酵工程及生物制藥。

*通訊作者：梁智群（1959-），男，教授，博士，研究方向為生物工程學。