全 文 :第 13卷第 6期
2015年 11月
生 物 加 工 过 程
Chinese Journal of Bioprocess Engineering
Vol 13 No 6
Nov 2015
doi:10 3969 / j issn 1672-3678 2015 06 012
收稿日期:2015-03-11
基金项目:国家自然科学基金(21306084)
作者简介:章人川(1990—),男,福建福州人,硕士研究生,研究方向:微生物次级代谢产物;黄 和(联系人),教授,E⁃mail:biotech@ njtech.edu.cn
基于渗透压调控的纽莫康定 B0发酵动力学
章人川,郑嘉琪,肖 潇,冯昆达,秦婷婷,宋 萍,黄 和
(南京工业大学 生物与制药工程学院,江苏 南京,211800)
摘 要:以海洋丝状真菌 Glarea lozoyensis Q45为研究对象,研究不同渗透压条件下 Glarea lozoyensis Q45的生长、底
物甘露醇消耗以及纽莫康定 B0的生产特性。 基于 Logistics方程、Luedeking Piret方程和 Luedeking Piret Like 方
程,得到了描述纽莫康定 B0发酵模型的动力学参数。 基于发酵动力学模型,提出了在菌体生长稳定期添加 1 61
g / L NaCl的调控策略。 在该策略下,纽莫康定 B0的产量达到 1 75 g / L,比初始条件的结果提高 24 6%。
关键词:纽莫康定 B0;渗透压; 发酵动力学;Glarea lozoyensis
中图分类号:Q93 文献标志码:A 文章编号:1672-3678(2015)06-0065-05
Fermentation kinetic model of pneumocandin B0 based on
osmotic pressure stage⁃controlled
ZHANG Renchuan,ZHENG Jiaqi,XIAO Xiao,FENG Kunda,QIN Tingting,SONG Ping,HUANG He
(College of Biotechnology and Pharmaceutical Engineering,Nanjing Tech University,Nanjing 211800,China)
Abstract:The fermentation kinetic models for Glarea lozoyensis Q45 were established to describe the effect of
different osmotic pressure conditions for cell growth,consumption of mannitol and production of pneumocandin
B0.Based on Logistics equation,Luedeking⁃Piret equation and Luedeking⁃Piret Like equation, the kinetic
parameters of different osmotic pressure were analyzed.An osmotic pressure stage⁃controlled strategy,with 1 61
g / L NaCl added during the stationary phase,was suggested and experimentally verified.Following this strategy,
the final production of pneumocandin B0 was 1 75 g / L,increased by 24 6%
Keywords:pneumocandin B0;osmotic pressure;fermentation kinetic;Glarea lozoyensis
海洋丝状真菌 Glarea lozoyensis 的次级代谢产物
纽莫康定 B0是抗真菌药物卡泊芬净的前体,由一个
环状的六肽母核以及一个脂肪酸链组成。 卡泊芬净
作为一种新型的广谱抗真菌药物———棘白菌素类抗
真菌药物,能选择性地抑制真菌细胞壁 β 1,3 葡聚
糖的合成,在抗真菌药物领域受到人们广泛的
关注[1-3]。
渗透压对大部分微生物的生长具有显著影响,细
胞会通过合成用于调节胞内渗透压的物质来适应环
境的变化[4]。 王玉磊等[5]在研究产朊假丝酵母发酵
的盐胁迫时发现,添加适当的 Na+和 K+能提高谷胱甘
肽和 S 腺苷蛋氨酸的联产产量。 曾艳等[6]研究裂殖
壶菌时发现,一定浓度的渗透压胁迫使二十二碳六烯
酸(DHA)产量提升 46 84%。 姜玮等[7]在研究游动
放线菌合成阿卡波糖的过程中发现,在前期维持较低
的渗透压 772 5 kPa,在中后期维持 1 030~1 288 kPa,
最终 50 L发酵罐的效价提升 50%。
目前,相关学者已就渗透压对 G. lozoyensis 发酵
的影响[8]以及 G.lozoyensis发酵产纽莫康定 B0的工业
放大[9 11]进行过研究,但对 G.lozoyensis发酵动力学进
行研究的很少。 发酵动力学能通过数学模型描绘微
生物代谢来减少传统经验型发酵调控的盲目性。 因
此,本研究中,笔者构建基于渗透压调控的G.lozoyensis
产纽莫康定 B0的分批发酵动力学模型,研究在最优
初始渗透压发酵条件下 G.lozoyensis 发酵过程中的菌
体生长,底物消耗以及产物合成的潜在规律,并初步
探讨渗透压对菌体生理代谢的影响方式。
1 材料与方法
1 1 实验材料
1 1 1 菌种
菌株 Glarea lozoyensis Q45,中国典型培养物保
藏中心(CCTCC)保藏,保藏日期 2014年 9月 14日,
保藏编号为 CCTCC M 2014416。
1 1 2 培养基
PDA固体培养基(g / L):去皮马铃薯 200 0,蔗
糖 20 0,琼脂 20 0。
种子培养基(g / L):葡萄糖 40 0,黄豆粉 20 0,
玉米浆 10 0,KH2PO4 10 0,微量元素 10 0。
微量元素(g / L):FeSO4·7H2O 1 0,MnSO4·H2O
1 0,ZnSO4·7H2O 0 2,CaCl2·2H2O 0 1,HBO3 0 056,
CuCl2·2H2O 0 025,(NH4)5MO7O24·4H2O 0 003。
发酵培养基(g / L):甘露醇 80,玉米蛋白粉 10,
K2HPO4 2 5。
1 1 3 实验材料
甘露醇标准品购自西格玛奥德里奇(上海)贸
易有限公司,其余试剂均购自国药集团试剂公司的
分析纯试剂。
LC 20AB型高效液相色谱购自岛津(中国)有
限公司,U3000型高效液相色谱购自赛默飞世尔科
技(中国)有限公司。
1 2 实验方法
1 2 1 培养方法
发酵培养:自斜面取 2 cm2菌苔至 50 mL种子培
养基中,26 ℃、220 r / min 培养 5 d。 以 10%的接种
量接入装有 50 mL发酵培养基的 250 mL锥形瓶中,
26 ℃、220 r / min培养 15 d。
1 2 2 检测方法
甘露醇检测采用岛津高效液相色谱 LC Solution工
作站测定。 色谱柱 Aminex HPX 87P柱(300 mm×7 8
mm),流动相为超纯水,流速 1 mL / min,检测器为示差
折光检测器(RID 10A),柱温 80 ℃,进样量 20 μL。
生物量检测采用菌体干质量浓度法,80 ℃烘干
至恒质量。
纽莫康定 B0 检测采用戴安高效液相色谱
U3000测定。 色谱柱 Venusil MP C18 柱(250 mm×
4 6 mm),流动相 A 为乙腈,流动相 B 为 0 1 %的
H3PO4。 洗脱程序 :0 ~ 20 min,V(A) ∶ V(B)= 40 ∶
60;20~40 min,V(A) ∶ V(B) = 50 ∶ 50,流速 1 mL /
min。 检测器为紫外检测器,检测波长 210 nm,柱温
35 ℃,进样量 20 μL。
1 2 3 渗透压计算方法
培养基中总渗透压为电解质和非电解质之和,
计算方法参照文献[12]。
目前,用于模拟菌体生长的非结构模型主要为
Monod方程和 Logistic 方程[13]。 G lozoyensis Q45 发
酵生产纽莫康定 B0的过程中,菌体生长符合 Logistic
方程,Logistic菌体生长动力学模型见式(1)。
dX
dt
= μm 1 -
X
Xm
æ
è
ç
ö
ø
÷ X (1)
式中:X为菌体质量浓度(g / L);Xm为最大菌体质量
浓度(g / L);μm 为最大比生长速率(d
-1); t 为发酵
时间(d)。
产物生成动力学模型:纽莫康定 B0的合成与菌
体的生长应归于部分生长偶联型。 因此,选用
Luedekin Piret方程来描述 G.lozoyensis Q45菌合成
纽莫康定 B0的速率[14],见式(2)。
dP
dt
= α dX
dt
+ βX (2)
式中:P表示纽莫康定 B0质量浓度(g / L),α表示生长
偶联的产物合成系数,β表示非生长偶联的产物合成系
数。 当 α≠0、β=0时,为产物形成与细胞生长偶联型;
α≠0、β≠0时,为产物形成与细胞生长部分偶联型;α=
0、β≠0时,为产物形成与细胞生长非偶联型。
底物消耗动力学模型:纽莫康定 B0分批发酵过程
中,甘露醇作为碳源和能源用于菌体生长,维持菌体新
陈代谢以及次级代谢产物的合成。 底物消耗动力学可
采用 Luedekin Piret Like方程来表示[15],见式(3)。
- ds
dt
= 1
YX/ S
dX
dt
+ mX + 1
YP / S
dP
dt
(3)
式中:ds
dt
表示甘露醇消耗速率(g·L-1·d-1),YX/ S表示
菌体相关系数,m 为菌体维持自身代谢与底物消耗
相关系数,YP / S表示产物相关系数。
66 生 物 加 工 过 程 第 13卷
发酵动力学参数计算:用 Matlab7 软件工具包
中 lsqcurvefit参数对上述公式进行非线性最小二乘
法拟合。
2 结果与讨论
图 1 渗透压对纽莫康定 B0发酵的影响
Fig 1 Effects of the osmotic pressure on
pneumocandin B0 fermentation
2 1 渗透压对 G. lozoyensis Q45 产纽莫康定 B0的
影响
考察渗透压对菌体底物消耗与产物合成的相
互影响,结果如图 1所示。 由图 1可知:高浓度甘露
醇 80 g / L(1 132 31 kPa)条件下纽莫康定 B0的产量
是低浓度 40 g / L(566 16 kPa)的 4 96 倍。 在低浓
度甘露醇条件下,添加终质量浓度为 6 44 g / L
(566 16 kPa)的 NaCl 使渗透压回到高浓度甘露醇
的渗透压,纽莫康定 B0 产量较低浓度时提高
191 4%。 底物消耗方面,NaCl 的添加使底物甘露
醇的消耗速度减缓。 因此,渗透压对 G. lozoyensis
Q45的发酵产纽莫康定过程和菌体代谢有着重要影
响。 但是纽莫康定 B0产量并没有回复到高浓度的
水平,可能是初始甘露醇浓度不够,没有足够的底
物维持细胞和产量的合成。 因此,在下阶段实验
中,采用 80 g / L甘露醇作为底物质量浓度。
图 2 NaCl质量浓度对纽莫康定 B0发酵过程的影响
Fig 2 Effects of NaCl concentration on pneumocand
in B0 fermentation
2 2 NaCl 浓度对 G. lozoyensis Q45 产纽莫康定
B0的影响
考虑到 6 42 g / L NaCl所提供的渗透压(相当于
40 g / L甘露醇)已经较大,所以在 80 g / L甘露醇基础
上添加 NaCl 的时候,以添加 6 42 g / L NaCl 为最高
值,以倍数关系做递减,考察不同初始 NaCl溶液条件
下纽莫康定 B0的发酵情况,结果如图 2所示。 由图 2
可知:G.lozoyensis Q45生长过程可分为 3个阶段。 在
迟滞期(Ⅰ阶段,0~2 d),碳源消耗速度缓慢,菌体生长
缓慢,同时产物少量合成;在对数期(Ⅱ阶段,2~6 d),
菌体快速生长,产物合成以及碳源消耗速度略微加
快;在稳定期(Ⅲ阶段,6~15 d),菌体生长趋于稳定,
此时,在碳源大量消耗的基础上,产物开始大量合成。
同时,随着 NaCl 浓度的增加,甘露醇消耗速度减缓,
生物量的积累减少,纽莫康定 B0的合成速度和最大
产量先增加后减少。 在迟滞期添加 1 61 g / L NaCl,纽
莫康定 B0产量相比添加 3 22 g / L 和 6 24 g / L NaCl
76 第 6期 章人川等:基于渗透压调控的纽莫康定 B0发酵动力学
的实验组分别提高 15 4%和 21 5%,因此,确定 1 61
g / L NaCl为最优添加量。
2 3 发酵动力学模型的建立与分析
2 3 1 发酵动力学模型的建立
建立 80 g / L 甘露醇的发酵动力学模型。 菌体
生长模型 X = e
0 73t
0 34+0 024e0 73t
,底物消耗模型 S =
80-0 43
e0 73t
1-0 067(1-e0 73t)
-1æ
è
ç
ö
ø
÷ -8 85ln(1-0 067·
(1-e0 73t)),产物生成模型 P = 167 79ln(1-0 067·
(1-e0 73t )) - 0 87
e0 73t
1-0 067(1-e0 73t)
-1æ
è
ç
ö
ø
÷ 。 模型相
关系数 R2分别为 0 997 8、0 994 8和 0 927 1。
为了验证模型的适用范围,分别以 8% ~ 12%初
始菌体浓度(即接种量)以及 60 ~ 100 g / L 初始底物
质量浓度作为考察范围,研究不同接种量计算值与实
验值的相对误差,结果如表 1 所示。 由表 1 可知:当
接种量在±20%之内时,菌体生长、底物消耗的实验值
与计算值的平均相对误差均在 5%之内,产物合成的
平均相对误差在 10%以内。 当初始底物质量浓度在
70~90 g / L范围内时,菌体生长、底物消耗和产物合
成模型的相对误差均控制在 5%以内,当初始底物质
量浓度在 60和 100 g / L时,菌体生长和产物合成的平
均相对误差均超过 10%,说明底物浓度相比发酵接种
量对发酵阶段菌体生长的影响更大(表 2)。
表 1 接种量计算值与实验值的相对误差
Table 1 Comparison between experimental dates and
predictive values of different inoculum sizes
接种
量 / %
菌体生长平均
相对误差 / %
底物消耗平均
相对误差 / %
产物合成平均
相对误差 / %
8 3 36 1 21 7 45
9 1 32 0 68 3 98
11 1 52 0 96 2 01
12 2 69 1 04 8 65
表 2 底物浓度计算值与实验值的相对误差
Table 2 Comparison between experimental dates and
predictive values of different substrate
concentrations
ρ(初始底物) /
(g·L-1)
菌体生长平均
相对误差 / %
底物消耗平均
相对误差 / %
产物合成平均
相对误差 / %
60 15 51 12 21 11 69
70 3 21 2 61 4 25
90 3 25 3 51 2 86
100 20 32 10 73 13 29
该模型在初始菌体质量浓度(2 76±0 552) g / L、
初始底物质量浓度(80±10) g / L 范围内能得到很好
的应用。
2 3 2 渗透压对发酵动力学参数影响
添加 NaCl后,菌体生长方程、底物消耗方程以及
产物合成方程分别为 X = e
0 67t
0 33 + 0 031e0 67t
, S =
80 -0 52
e0 67t
1 - 0 086(1 - e0 67t)
- 1æ
è
ç
ö
ø
÷ - 6 39ln(1 -
0 086(1 - e0 67t)) ,P = 17 46 ( e
0 67t
1 - 0 086(1 - e0 67t)
-
1 ) + 174 7ln(1 - 0 086(1 - e0 67t)) ,NaCl的添加使最
大生物量 Xm降低 21 68%,最大菌体比生长速率 μm降
低 7 72%。 由此可知 NaCl添加所产生的渗透压对菌
体产生一定的毒性作用,降低了菌体的生长速度。 比
较最大菌体相关系数 YX/ S与最大产物相关系数 YP / S
可以发现, NaCl的添加使最大菌体相关系数降低,用
于产物合成的相关系数提高,说明此时消耗的底物更
多地向产物合成方向分配。 维持菌体生长的底物消
耗相关系数可以进一步证明,NaCl 的添加使底物消
耗更多,用于产物的合成而不是菌体的生长。
2 4 基于发酵动力学研究的渗透压调控策略
基于以上的动力学参数分析,NaCl 对纽莫康定
B0发酵的影响如下:①降低菌体的生长速度;②促
使底物的消耗流向纽莫康定 B0合成方向。 考虑到
纽莫康定 B0是次级代谢产物,产物合成需要生物量
的积累[16],因此,研究不同阶段 NaCl的添加对产物
合成的影响,结果如图 3所示。 由图 3可知:稳定期
(12 d)添加 1 61 g / L NaCl使纽莫康定 B0的产量提
升 24 6%,与在迟滞期(0 天)和对数期(4 d)添加
1 61 g / L NaCl 相比,分别提高 21 28%和 15 47%。
因此,此结果验证了基于发酵动力学参数分析提出
的渗透压调控策略是可行的。
可以从以下两个方面尝试解释 NaCl 对纽莫康
定 B0合成存在的双面影响:①由于纽莫康定 B0合成
菌株 G.lozoyensis Q45属于海洋丝状真菌,其天然的
生长环境本身就存在一定程度的渗透压,因此 NaCl
的添加一定程度上模拟了菌体的最适生长环境,导
致产物合成的增加;②Na+的添加会改变微生物细
胞膜上 Ca2+通道的流量,使得细胞膜上的磷脂酶信
号被激活,从而影响钙调磷激酶响应的锌指蛋白,
最终提高胞内次级代谢产物的合成[17], 因此,推测
86 生 物 加 工 过 程 第 13卷
在纽莫康定 B0的合成中,Na
+也起到同样效果。
图 3 不同生长阶段 NaCl添加对产量的影响
Fig 3 Effects of NaCl added in different growth
phase for production
3 结论
通过添加不同浓度 NaCl,研究不同渗透压下
G. lozoyensis产纽莫康定 B0的分批发酵动力学模型,该
菌株生长符合 Logistic 方程,底物消耗符合 Luedekin
Piret Like方程,产物合成符合 Luedekin Piret方程,
动力学方程分别为 X = e
0 728
0 337 9 + 0 024 3e0 728t
,
S = 80 - 0 429 2
e0 728t
1 - 0 067 2(1 - e0 728t)
- 1æ
è
ç
ö
ø
÷ -
8 847 1ln(1-0 067 2(1-e0 728t)) , P = - 0 873 5·
e0 728t
1 - 0 067 2(1 - e0 728t)
- 1æ
è
ç
ö
ø
÷ + 167 79ln(1 - 0 067 2
(1 -e0 728t)) ,NaCl溶液添加会降低菌体的最大比
生长速率、最大生物量以及菌体用于自身代谢的相
关系数,同时会提高产物合成相关系数。
根据不同渗透压条件下动力学参数的分析,提
出在生长稳定期添加 NaCl 的调控策略。 在前期维
持较高的菌体比生长速率,使得菌体达到更高的生
物量,在后期添加 NaCl 溶液,降低菌体用于维持自
身代谢的相关系数,提高产物合成相关系数,使菌
体积累足够生物量的同时实现产物的高浓度积累。
对所提出的发酵控制策略进行验证,通过在稳定期
中期流加 NaCl的策略,产量与 80 g / L条件相比,提
高 24 6%,达到 1 75 g / L。
参考文献:
[ 1 ] Wichmann C F, Liesch J M, Schwartz R E. L⁃671 329, a new
antifungal agent:Ⅱ.structure determination[J] .J Antibiot,1989,
42(2):168⁃173.
[ 2 ] Fromtling R A,Abruzzo G K.L⁃671 329,a new antifungal agent:
Ⅲ.invitro activity,toxicity and efficacy in comparison to aculeacin
[J] .J Antibiot,1989,42(2):174⁃178.
[ 3 ] Sawistowskaschroder E T, Kerridge D, Perry H. Echinocandin
inhibition of 1,3⁃beta⁃d⁃glucan synthase from Candida albicans
[J] .FEBS Lett,1984,173(1):134⁃138.
[ 4 ] Yancey P H,Clark M E,Hand S C,et al.Living with water⁃stress⁃
evolution of osmolyte systems[J] .Science,1982,217:1214⁃1222.
[ 5 ] 王玉磊,叶子优,贺秋萍,等.S 腺苷蛋氨酸和谷胱甘肽联产
发酵中盐胁迫的作用[J] .生物加工过程,2012,10(3):33⁃38.
[ 6 ] 曾艳,冯云,仝倩倩,等.渗透压调控对裂殖壶菌发酵产 DHA
的影响[J] .生物加工过程,2012,10(6):24⁃28.
[ 7 ] 姜玮,生英涛,蔡月明,等.麦芽糖浓度和渗透压对游动放线
菌生长及阿卡波糖生物合成的综合影响[ J] .中国医药工业
杂志,2010,41(3):178⁃182.
[ 8 ] Connors N, Petersen L,Hughes R, et al. Residual fructose and
osmolality affect the levels of pneumocandins B0 and C0 produced
by Glarea lozoyensis [ J] . Appl Microbiol Biotechnol, 2000, 54
(6):814⁃818.
[ 9 ] Pollard D J, Buccino R, Connors N C, et al. Real⁃time analyte
monitoring of a fungal fermentation, at pilot scale, using in situ
mid⁃infrared spectroscopy[J].Bioproc Biosyst Eng,2001,24(1):
13⁃24.
[10] Pollard D J, Hunt G R, Kirschner T F, et al. Rheological
characterization of a fungal fermentation for the production of
pneumocandin[J] .Bioproc Biosyst Eng,2002,24(6):373⁃383.
[11] Pollard D J,Kirschner T F,Hunt G R,et al.Scale up of a viscous
fungal fermentation: application of scale⁃up criteria with regime
analysis and operating boundary conditions [ J ] . Biotechnol
Bioeng,2007,96(2):307⁃317.
[12] Wright P J,Reed R H.Method for determination of turgot pressure
in macroalgae,with particular reference to the phaeophyta [ J] .
Marine Biology,1988,99:473⁃480.
[13] Feng J,Zhang J S,Jia W,et al.An unstructured kinetic model for
the improvement of triterpenes production by Ganoderma lucidum
G0119 based on nitrogen source effect [ J] . Bioprocess Biosyst
Eng,2014,19(4):727⁃732.
[14] Luedeking R, Piret E L. A kinetic study of the lactic acid
fermentation batch process at controlled pH [ J ] . Biotechnol
Bioeng,2000,67(6):636⁃644.
[15] 张扬,冯小海,李莎,等.Kitasatospora sp.MY5 36 产 ε 聚赖
氨酸分批补料发酵动力学[ J] .生物加工过程,2012,10(3):
23⁃27.
[16] Papagianni M. Fungal morphology and metabolite production in
submerged mycelial processes[J] .Biotechnol Adv,2004,22(3):
189⁃259.
[17] Xu Y N, Xia X X, Zhong J J. Induced effect of Na plus on
ganoderic acid biosynthesis in static liquid culture of Ganoderma
lucidum via calcineurin signal transduction [ J ] . Biotechnol
Bioeng,2013,110(7):1913⁃1923.
(责任编辑 管 珺)
96 第 6期 章人川等:基于渗透压调控的纽莫康定 B0发酵动力学