Screening of A Neutral Protease-producing Strain and Optimization of Fermentation Conditions
-
摘要: 目前,低酶活及原料高成本严重制约了蛋白酶工业的发展,本研究采用透明圈法和摇瓶发酵法筛选出1株产中性蛋白酶的菌株,经形态学观察、生理生化实验以及16S rDNA序列分析鉴定为贝莱斯芽孢杆菌(Bacillus velezensis),并命名为ppr3,通过单因素实验和响应面试验对发酵培养基成分及发酵条件进行优化,结果表明:当发酵温度33 ℃,接种量8.5%(v/v),初始pH5.7,乳糖5.55 g/L、菜粕31.89 g/L、酵母浸粉7.09 g/L、吐温80 0.03%、MgSO4 0.04%、K2HPO4 0.04%时,蛋白酶活为500 U/mL,相较于未优化前提高了132%。将菜粕作为培养基的发酵氮源,在一定程度上降低了成本,为农副产品的再利用提供了一个新途径,也为后续工业化生产奠定了研究基础。Abstract: Currently, low enzyme activity and high raw material cost seriously restricted the development of the industrial production of protease. In this study, transparent circle and shaking flask fermentation methods were used to screen out a neutral protease-producing strain, which was identified as Bacillus velezensis by morphological observation, physiological and biochemical results and 16S rDNA analysis, and named ppr3. The medium composition and fermentation conditions for this strain to produce neutral protease were optimized by one-at-a-time strategy and response surface methodology. The results showed that: The fermentation temperature was 33 ℃, inoculation amount was 8.5% (v/v), initial pH was 5.7, lactose 5.55 g/L, rapeseed 31.89 g/L, yeast extract 7.09 g/L, tween-80 0.03%, MgSO4 0.04%, K2HPO4 0.04%, the protease activity was 500 U/mL. Compared with before optimization, it increased by 132%. The addition of rapeseed meal as nitrogen source to the medium significantly reduced the fermentation cost and provided a new way for the reuse of agricultural and sideline products, which had major implications in protease industrial fermentation.
-
蛋白质是体内重要的营养物质。摄入的蛋白质在胃肠道中被蛋白酶和肽酶降解,被吸收为小肽或游离氨基酸[1]。小肽比氨基酸更易被吸收,且具有增强免疫、抗氧化、抗菌、降血脂等多种生物学功能,被广泛应用于食品和饲料生产中[2-4]。含蛋白质的农业副产物对生物活性肽的生产越来越具有吸引力,这有助于废物的利用和价值的提高,如豆粕,菜粕,棉粕,花生粕等[5]。这些肽序列在含有完整的母蛋白序列时不显示任何生物活性,因此首先必须水解,一般是通过与一种或多种蛋白酶酶解、发酵或两者结合而催化水解。在发酵过程中,发酵底物中蛋白质的分子结构、空间结构效应以及理化性质会发生变化,产生独特的功能肽,如鲜味肽、抗菌肽[6-8]。但发酵效果易受微生物种类和所用发酵工艺的影响,乳酸菌、酵母菌在豆粕、菜粕脱毒和提升风味中效果显著,但由于其蛋白酶活低,使得豆粕、菜粕中的蛋白不能有效转化[9],菌酶协同是一种好的方式[10]。外源添加蛋白酶提高小肽的含量的同时,还可有效降低蛋白原料中的过敏原[11]。
中性蛋白酶具有催化速度快、催化过程温和、无污染等优点,常被用于食品、药品、饲料行业[12]。中性蛋白酶产生的游离氨基酸类型丰富,制备呈味肽的效果优于风味蛋白酶,其产生的不益物质以及水解度要比风味蛋白酶低,因此中性蛋白酶用于发酵工艺可以降低酶的用量并提高原料的回收率[13-14]。蛋白质是微生物产蛋白酶主要的底物来源,植物源蛋白可作为蛋白酶及多肽生产的潜在底物,如中性蛋白酶酶解红花籽粕制备抗氧化肽,酶解大豆粕蛋白制备小分子肽[15-16]。
植物源蛋白的应用逐渐成为当今研究的热点,我国非传统蛋白资源非常丰富,包括农产品加工副产物、糟粕类、植物及其加工副产物等,这些植物源蛋白的利用有助于减少酶或多肽的生产成本,同时也会获得多种活性功能的多肽。因此,适应菜粕等农业副产物的高性能、高产量的蛋白酶新型菌株和连续发酵工艺亟需推进。本研究从桑叶中筛选出一株高产中性蛋白酶的菌株,以菜粕作为氮源,通过优化培养基成分和发酵条件,使得酶活显著提高。该研究为植物源蛋白资源高值利用和植物蛋白向生物活性肽的转化提供了理论基础。
1. 材料与方法
1.1 材料与仪器
桑叶 −20 ℃下保藏,备用(品种为白桑叶,2021年5月采摘自某桑叶园);干酪素、三氯乙酸、硫酸铜、无水碳酸钠、葡萄糖、乳糖、曲拉通-100、吐温60、吐温80、尿素、氯化铵、硫酸镁、可溶性淀粉、氯化钠、氯化钾、硫酸锰、硫酸锌、硫酸铁、硫酸铜、麦芽糖 分析纯,均购于成都科隆化学品有限公司;福林酚试剂、蛋白胨、牛肉膏、酵母浸粉 生化级,均购于成都奥博星生物技术有限公司;蔗糖 化学纯,购于天津市科密欧化学试剂有限公司;司班60 化学纯,购于成都科隆试剂有限公司;菜粕、豆粕、玉米粉 食品级,市售;DNA试剂盒 购自天根生化科技有限公司。
酪素培养基(g/L):参考文献[17]并稍加修改,干酪素4 g,并滴入少量20 g/L NaOH溶液和蒸馏水,在沸水浴下溶解后调pH至7.5,再加入1.07 g Na2HPO4.7 H2O作为A液。0.36 g K2HPO4加入少量蒸馏水溶解作为B液。A液和B液混合加入20 g琼脂粉、0.5%酪素水解培养物2 mL,溶于1 L蒸馏水中,121 ℃灭菌15 min。
营养肉汤培养基(g/L):蛋白胨10 g/L、牛肉浸出粉3 g/L、氯化钠5 g/L、pH7.2±0.2,溶于1 L蒸馏水中,121 ℃灭菌15 min。
基础发酵培养基(g/L):葡萄糖2 g/L、蛋白胨10 g/L、氯化钠5 g/L、吐温80 0.4%,溶于1 L蒸馏水中,121 ℃灭菌15 min。
T-114分析天平 奥豪斯仪器有限公司;QYC-2102C恒温培养摇床 上海福马实验设备有限公司;DZKW-4电子恒温水浴锅、T6新世纪紫外可见分光光度计 北京普析通用仪器有限公司;HWS-260B恒温恒湿箱、LS-75HD高压蒸汽灭菌锅 江阴滨江医疗设备有限公司;BM1000显微镜 南京江南永新光学有限公司;AT-710 pH计 日本京都电子制造有限公司;SPX-100B-Z单人单面净化超净工作台 上海沪净医疗器械有限公司;Mini Amp PCR仪 上海艾研生物科技有限公司。
1.2 实验方法
1.2.1 产蛋白酶菌株的筛选
初筛:取20 g桑叶切碎加入150 mL生理盐水中,37 ℃、195 r/min摇床培养20 min,用梯度稀释法分别稀释10−1~10−5浓度,分别从10−1~10−5稀释液中取出100 μL涂布在酪素培养基平板上,于37 ℃培养24 h,将产生透明圈菌株点注在酪素培养基上,观察菌落与透明圈直径大小,直径较大的可以初步认为是产蛋白酶能强的菌株,并用接种环挑取产透明圈菌株进行纯种划线。
复筛:挑取产生透明圈的菌株接种于营养肉汤培养基中,37 ℃、195 r/min摇床培养35 h,测定中性蛋白酶活,结合初筛结果,挑选产中性蛋白酶活较高的菌株。
1.2.2 蛋白酶活力测定
根据GB/T 23527-2009《蛋白酶制剂》进行蛋白酶活测定,碱性蛋白酶pH10.5、中性蛋白酶pH7.5、酸性蛋白酶pH3.0条件下测定。以L-酪氨酸为标准曲线并建立回归方程,吸光度值为纵坐标,L-酪氨酸浓度为横坐标,回归方程为y=0.0106x+0.0021,R2=0.9998。
1.2.3 菌株鉴定
1.2.3.1 菌株形态学鉴定
所筛菌株在营养肉汤固体培养基上划线,37 ℃好氧培养24 h,观察菌落,通过革兰氏染色在显微镜下观察菌体形态。
1.2.3.2 生理生化特征分析
根据《伯杰氏细菌鉴定手册(第八版)》对该菌种进行糖发酵试验、VP试验、MR试验、吲哚试验、柠檬酸盐利用试验、苯丙氨酸脱羧酶测定、H2S试验、产气试验。
1.2.3.3 分子生物学鉴定
采用DNA试剂盒对筛选到的菌株进行全基因组提取,以27F(5’-AGAGTTTGATCCTGGCTCAG -3’)、1492R(5’-GGTTACCTTGTTACGACTT-3’)为引物对16S rDNA进行PCR扩增,经纯化后将样品送擎科生物公司测序,测序结果在NCBI数据库中进行BLAST分析,之后利用MEGA6进行基因排序匹配及系统发育树的构建,得到基因序列同源性较高的菌株。
1.2.4 菌株产酶曲线
对数生长期的菌液以3%(v/v)接种至营养肉汤培养基,37 ℃、195 r/min摇床发酵,每隔8 h取样,测定蛋白酶活,根据蛋白酶活的变化,确定最佳产酶时间。
1.2.5 不同培养基成分对菌株产酶的影响
1.2.5.1 不同碳源对菌株产酶的影响
对数生长期的菌液以3%(v/v)接种至基础发酵培养基,其中葡萄糖替换为2 g/L(蔗糖、乳糖、麦芽糖、果糖、可溶性淀粉、玉米粉),其余成分及含量不变,初始pH5.6、37 ℃、195 r/min摇床发酵40 h,测定蛋白酶活。最佳碳源以0.5、3.5、6.5、9.5、12.5 g/L添加进不同基础发酵培养基中发酵40 h后测定蛋白酶活,选取最优区间进行后续试验设计。
1.2.5.2 不同氮源对菌株产酶的影响
对数生长期的菌液以3%(v/v)接种至基础发酵培养基,其中蛋白胨替换为10 g/L(牛肉膏、菜粕、豆粕、尿素、硝酸钾、硫酸铵),其余成分及含量不变,初始pH5.6、37 ℃、195 r/min摇床发酵40 h,测定蛋白酶活。最佳氮源按1.5、11.5、19.5、27.5、35.5、43.5 g/L添加进不同基础发酵培养基中发酵40 h后测定蛋白酶活,选取最优区间进行后续试验设计。
1.2.5.3 不同金属离子对菌株产酶的影响
对数生长期的菌液以3%(v/v)接种至基础发酵培养基,其中NaCl替换为5 g/L(MgSO4、ZnSO4、CuSO4、FeSO4、K2SO4、MnSO4),其余成分及含量不变,初始pH5.6、37 ℃、195 r/min摇床发酵40 h,测定蛋白酶活。最佳金属离子按0.5、3、5.5、8、10.5 g/L添加进不同基础发酵培养基中发酵40 h后测定蛋白酶活,选取最优区间进行后续试验设计。
1.2.5.4 不同表面活性剂对菌株产酶的影响
对数生长期的菌液以3%(v/v)接种至基础发酵培养基,其中吐温80分别替换为0.4%(吐温60、Triton-100、司班60),其余成分及含量不变,初始pH5.6、37 ℃、195 r/min摇床发酵40 h,测定蛋白酶活。最佳表面活性剂按0.2、0.6、1、1.4、1.8 g/L添加进基础发酵培养基中发酵40 h后测定蛋白酶活,选取最优区间进行后续试验设计。
1.2.5.5 酵母浸粉和磷元素对菌株产酶的影响
对数生长期的菌液以3%(v/v)接种至基础发酵培养基,分别添加0、1、2、3、4 g/L的K2HPO4,0、3、6、9、12 g/L酵母浸粉,初始pH5.6、37 ℃、195 r/min摇床发酵40 h后测定蛋白酶活,探讨酵母浸粉和磷元素对蛋白酶活的影响。
1.2.6 发酵培养基优化试验
1.2.6.1 Plackett-Burman(P-B)试验
根据单因素实验得到的结果。采用P-B试验设计对6个培养基成分(乳糖、菜粕、MgSO4、酵母浸粉、吐温80、K2HPO4)进行方差分析,将每个因素分为两个水平,即低水平-1,高水平+1表示,试验设计如表1。
表 1 Plackett-Burman试验设计Table 1. Plackett-Burman experimental design因素 水平(g/L) −1 1 A 乳糖 0.5 4.5 B 菜粕 2 18 C MgSO4 0.4 3.5 D 酵母浸粉 0.25 3.75 E 吐温80 0.3 1.0 F K2HPO4 0.4 1.5 1.2.6.2 最陡爬坡试验
根据P-B试验设计得到的3个显著因素进行最陡爬坡试验,综合考虑各因素之间的正负影响并确定各因素步长和爬坡方向,其余成分按照低水平加入,试验设计如表2。
表 2 最陡爬坡试验设计Table 2. Design of steepest climbing test组数 成分 A 乳糖(g/L) B 菜粕(g/L) D 酵母浸粉(g/L) 1 0.5 2 0.25 2 2.5 10 2.00 3 4.5 18 3.75 4 6.5 26 5.50 5 8.5 34 7.25 6 10.5 42 9.00 7 12.5 50 10.75 1.2.6.3 Box-Behnken响应面优化设计
在上述试验基础上,利用Design-Expert 8.0软件进行响应面优化。以培养基成分中各显著因素为自变量,蛋白酶活为响应值(Y),设计3因素3水平试验,其余成分按照低水平加入。根据响应面试验模拟的最佳产酶的培养基成分及相应浓度进行验证。试验设计如表3。
表 3 培养基成分响应面优化设计Table 3. Response surface optimization design of medium components因素 水平(g/L) −1 0 1 A 乳糖 4.5 6.5 8.5 B 菜粕 18 26 34 D 酵母浸粉 3.75 5.50 7.25 1.2.7 发酵条件对菌株产酶的影响
除了培养基成分影响菌株产酶,发酵条件也是影响菌株产酶活力的重要因素。将菌体接入优化培养基中,过夜培养至对数生长期(OD≈0.5)作为种子液并进行如下实验。
1.2.7.1 温度对菌株ppr3产酶的影响
对数生长期的菌液以3%(v/v)接种至优化的培养基中,置于18、26、34、42、50 ℃,初始pH5.6、195 r/min摇床发酵40 h后测定蛋白酶活。
1.2.7.2 接种量对菌株ppr3产酶的影响
对数生长期的菌液分别以0.5%、3%、5.5%、8%、10.5%(v/v)接种至优化的培养基中,初始pH5.6、37 ℃、195 r/min摇床发酵40 h后测定蛋白酶活。
1.2.7.3 初始pH对菌株ppr3产酶的影响
对数生长期的菌液以3%(v/v)接种至优化的培养基中,初始pH分别调节为1、3、5、7、9,37 ℃、195 r/min摇床发酵40 h后测定蛋白酶活。
1.2.8 发酵条件的响应面优化设计
根据单因素实验结果,采用中心组合设计获得二次模型,中心点估计以蛋白酶活为指标,对3个因素进行响应面优化,根据响应面试验模拟最佳产酶的发酵条件进行验证。试验设计如表4。
表 4 发酵条件响应面优化设计Table 4. Response surface optimization design of fermentation conditions因素 水平 −1 0 1 a 发酵温度(℃) 26 34 42 b pH 3 5 7 c 接种量(%) 0.5 3 5.5 1.3 数据处理
运用SPSS 22进行单因素ANOVA分析,采用Desige-Expert12软件进行Plackett-Burman试验和Box-Behnken响应面优化试验设计,建立非线性回归数学建模,并对模型进行验证,该模型考虑了各因素之间的主次效应以及因素之间的二级交互作用。
2. 结果与分析
2.1 产蛋白酶菌株的筛选
从采集的桑叶中筛选得到7株产生透明圈的菌株,挑取7株菌株点注在酪素培养基,见图1,分别命名为ppr0、ppr2、ppr3、ppr4、ppr5、ppr6、ppr7,根据透明圈与菌株直径比例,初步筛选酶活较高的菌株(表5),比值最大为菌株ppr3。通过发酵实验对中性蛋白酶菌株复筛,结果显示,中性蛋白酶活最高的仍为ppr3菌株,达到215.4±2.26 U/mL。表6展示了该菌株在不同蛋白酶测定方法下的蛋白酶活,中性蛋白酶活最高,为220.4±3.47 U/mL,因此菌株ppr3产中性蛋白酶,后续试验以ppr3为试验菌株进行中性蛋白酶酶活的优化。
表 5 透明圈与菌株直径的比值及相应酶活Table 5. Ratio of transparent circle to strain diameter and corresponding enzyme activities表 6 不同pH下ppr3蛋白酶活Table 6. ppr3 Protease activity at different pH蛋白酶种类 pH ppr3蛋白酶活(U/mL) 酸性蛋白酶 3.0 138.17±6.77c 中性蛋白酶 7.5 220.40±3.47a 碱性蛋白酶 10.5 188.30±4.32b 表 7 Plackett-Burman试验设计及结果Table 7. Plackett-Burman experimental design and results试验号 A
乳糖
(g/L)B
菜粕
(g/L)C
MgSO4
(g/L)D
酵母浸粉
(g/L)E
吐温80
(g/L)F
K2HPO4
(g/L)Y 蛋白酶活
(U/mL)1 4.5 18 3.5 0.5 0.3 0.4 270±9.9d 2 0.5 18 0.4 3.5 1 0.4 234±4.24f 3 4.5 2 3.5 3.5 0.3 1.5 297±10.26b 4 4.5 18 0.4 0.5 0.3 1.5 287±3.45c 5 0.5 2 0.4 0.5 0.3 0.4 172±5.23h 6 4.5 18 0.4 3.5 1 1.5 322±2.55a 7 0.5 2 3.5 0.5 1 1.5 157±6.01i 8 0.5 18 3.5 3.5 0.3 0.4 264±4.66d 9 4.5 2 3.5 3.5 1 0.4 210±3.24g 10 0.5 18 3.5 0.5 1 1.5 254±5.77e 11 0.5 2 0.4 3.5 0.3 1.5 240±3.21f 12 4.5 2 0.4 0.5 1 0.4 232±2.28f 2.2 菌株鉴定
2.2.1 菌株形态学观察
由图2(左)可知,菌株近乳白色,边缘光滑规则,呈圆形,表面湿润。革兰氏染色结果阳性,菌体形态杆状(图2(右))。
2.2.2 生理生化特征
生理生化结果显示,菌株ppr3可分解葡萄糖、乳糖、蔗糖,不产气,VP试验、MR试验为阴性,不具有色氨酸酶、苯丙氨酸脱羧酶,不可利用柠檬酸盐作为碳源,不能分解含硫氨基酸,结合《伯杰氏细菌鉴定手册》,初步确定为芽孢杆菌属。
2.2.3 ppr3菌株的分子生物学鉴定
对菌株ppr3进行DNA提取,PCR扩增产物经琼脂糖凝胶电泳,如图3所示,ppr3菌株16S rDNA在1000~2000 bp之间,表明已成功扩增目的基因。建立系统发育树,比较菌株之间的同源性,如图4所示,ppr3与贝莱斯芽孢杆菌CBMB205同源性极高。结合菌株形态学观察、生理生化分析,初步认定ppr3菌株为贝莱斯芽孢杆菌,贝莱斯芽孢杆菌是一种新型益生菌,具有抗菌和改善肠道的能力[20],本研究筛选的贝莱斯芽孢杆菌可能在食品与农业生产中具有广泛潜力。
2.3 菌株的产酶曲线
测定ppr3不同发酵时间段的中性蛋白酶活,绘制产酶曲线。由图5(a)可知,菌株发酵过程中,蛋白酶活先升高后下降,40 h达到277.3±6.21 U/mL,最佳产酶时间为40 h,因此取40 h为发酵截止时间进行后续试验。
图 5 菌种的产酶曲线及发酵培养基成分对酶活的影响注:不同培养基成分的蛋白酶活(左轴),不同浓度的培养基成分蛋白酶活(右轴);b~e图中柱状图表示不同浓度的培养基成分蛋白酶活,折线图表示不同浓度的培养基成分蛋白酶活;图(f)左轴代表酵母浸粉蛋白酶活,右轴代表K2HPO4蛋白酶活;同一指标下不同字母表示差异显著(P<0.05),图7同。Figure 5. Enzyme production curve of strain and effect of medium composition on enzyme activity2.4 不同培养基成分对蛋白酶活的影响
添加不同碳源对菌株ppr3产酶的影响如图5(b)所示。结果表明,以乳糖为碳源时,蛋白酶活最高,酶活达297±6.66 U/mL,这与之前研究结果一致[21],乳糖作为缓速碳源,有助于蛋白酶等代谢产物的合成[22]。当乳糖浓度为3.5 g/L时,蛋白酶活最高达310±3.33 U/mL,随着乳糖浓度增加,酶活逐渐下降,这可能是微生物倾向于更多碳源代谢,影响了蛋白酶的产生,其他代谢产物的积累也抑制了菌体的生长。因此选择4.5和0.5 g/L作为P-B试验碳源的高低水平。
探究不同氮源对蛋白酶活的影响发现(图5(c)),以菜粕和豆粕为氮源时,蛋白酶活最高,分别达到了382.83±4.45、385.35±6.19 U/mL,二者差异不显著(P>0.05),相比于无机氮源,部分有机氮源更利于菌株产生蛋白酶。菜粕中含有丰富的蛋白质,可供蛋白酶的产生与利用[23-24],且菜粕价格优惠,故选择菜粕作为最佳氮源进行后续实验。当菜粕浓度为11.5 g/L时,蛋白酶活最高,为403±4.88 U/mL。继续增加菜粕浓度,酶活出现先平稳后下降的趋势。菜粕浓度的增加,发酵液呈粘性状态,当发酵液粘性过大时,可能不利于菌体产酶。因此选择18和2 g/L作为P-B试验氮源的高低水平。
不同金属离子对蛋白酶活的影响如图5(d)所示,添加Mg2+发酵组蛋白酶活最高,为305±8.67 U/mL,Na+、K+次之,添加Mn2+、Zn2+、Fe2+、Cu2+酶活最低,且低于空白组。Mg2+是蛋白酶的活性位点和活化剂,催化蛋白酶的产生,细胞中的某类蛋白系统对Mg2+浓度进行了储存,并运输到适当位置,组合成酶促系统,从而促进了蛋白酶的产生[25-26],因此选择Mg2+为最佳金属离子。Mg2+浓度为3 g/L时,蛋白酶活达到291±5.12 U/mL,与Mg2+浓度为5.5 g/L时的蛋白酶活无显著差异(P>0.05)。之后随着Mg2+浓度增加,酶活出现轻微下降,说明金属离子浓度过高抑制酶的产生,因此选择 3.5 和 0.4 g/L 作为 P-B 试验 Mg2+ 的高低水平。
表面活性剂对蛋白酶活的影响见图5(e),添加曲拉通-100对蛋白酶活有抑制作用,其他3种表面活性剂与对照组相比蛋白酶活无显著差异(P>0.05),非离子型表面活性剂可以增大细胞膜的渗透性,使得酶更容易分泌[27-28],选取吐温80进行后续实验。吐温80浓度为0.6 g/L时,蛋白酶活最高,达到了294±6.33 U/mL,适量的表面活性剂对蛋白酶活有所提高,因此选择1.0和0.3 g/L作为P-B试验吐温80的高低水平。
添加不同浓度的酵母浸粉、K2HPO4发现(图5(f)),随着浓度的增加,蛋白酶活呈现先增加后下降的趋势,当酵母浸粉、K2HPO4浓度分别为3 g/L和1 g/L时蛋白酶活最高,分别为387±5.88 U/mL和301±3.38 U/mL。即适量的酵母浸粉和K2HPO4可提高蛋白酶活。酵母浸粉中含有大量微生物可利用的生长因子[29],K2HPO4中的P元素是微生物代谢的重要组成部分,可为蛋白酶起到pH缓冲作用[30]。因此选择3.75和0.25 g/L作为P-B试验酵母浸粉高低水平,1.5和0.4 g/L作为P-B试验磷元素高低水平。
2.5 响应面优化
2.5.1 Plackett-Burman试验
P-B用来评估实验误差和一段模型的充分性[31],通过P-B试验设计可获得影响显著的变量,从而为后续变量优化提供理论基础。
依据单因素实验结果,确定因素与水平后,通过P-B试验(表1)筛选影响显著的重要因素。如表7所示以蛋白酶活(Y)为响应值,对数据进行回归分析,得到回归方程为Y=24.75A+26.92B−2.92C+16.25D−10.08E+14.58F。模型相关系数为R2=0.8854,P=0.0296<0.05,该模型可用。P-B试验的方差分析如表8所示,从各因素P值可知,乳糖、菜粕对蛋白酶活影响显著(P<0.05)。根据F值将显著因素由大到小排列为菜粕>乳糖>酵母浸粉>K2HPO4>吐温80>Mg2+,根据回归方程各变量系数可知,乳糖、菜粕、酵母浸粉对蛋白酶活均表现为正效应,选择F值前三的影响因素进行后续爬坡实验。
表 8 Plackett-Burman 设计方差分析Table 8. Plackett-Burman design variance analysis因素 F检验 P值 重要性 A 12.30 0.0172 ※※※※※ B 14.54 0.0125 ※※※※※※ C 0.1708 0.6966 ※ D 5.30 0.0696 ※※※※ E 2.04 0.2125 ※※ F 4.27 0.0937 ※※※ 注:※数量表示因素显著性,数量越多显著性越强。 2.5.2 最陡爬坡试验
根据已确定的因素正负效应,并依据前期试验确定爬坡方向及步长。由表9可知,第4组试验蛋白酶活最高,为483.79±6.65 U/mL,因此将第4组作为中心点进行响应面优化。
表 9 最陡爬坡试验Table 9. Steepest ascent experiment试验号 A 乳糖(g/L) B 菜粕(g/L) D 酵母浸粉(g/L) 蛋白酶活(U/mL) 1 0.5 2 0.25 221.18±4.47f 2 2.5 10 2 363.02±1.33d 3 4.5 18 3.75 415.84±6.66c 4 6.5 26 5.5 483.79±6.65a 5 8.5 34 7.25 409.25±8.68c 6 10.5 42 9 258.31±2.67e 7 12.5 50 10.75 437.49±4.78b 表 10 Box-Behnken试验设计及结果Table 10. Box-Behnken experimental design and results试验号 A 乳糖(g/L) B 菜粕(g/L) D 酵母浸粉(g/L) Y 蛋白酶活(U/mL) 1 4.5 34 5.5 451.23±3.21f 2 6.5 34 7.25 476.56±2.89de 3 6.5 18 7.25 378.12±4.33m 4 6.5 26 5.5 491.58±1.98a 5 8.5 26 7.25 413.02±2.44h 6 8.5 18 5.5 382.13±7.02lm 7 6.5 34 3.75 441.32±4.77g 8 6.5 26 5.5 486.53±5.03abc 9 4.5 26 3.75 388.49±3.01i 10 6.5 26 5.5 479.52±6.02cde 11 6.5 26 5.5 488.56±4.67ab 12 4.5 18 5.5 363.57±3.58n 13 6.5 18 3.75 349.43±3.28o 14 8.5 26 3.75 473.40±2.65e 15 6.5 26 5.5 482.00±4.33bcd 16 8.5 26 7.25 473.21±1.98e 17 8.5 34 5.5 441.23±5.67g 2.5.3 培养基成分的响应面优化设计
以乳糖(A)、菜粕(B)、酵母浸粉(D)为自变量,蛋白酶活为响应值(Y),通过软件Design-Expert 12对表10进行多元回归拟合,得自变量A、B、D与因变量Y间的编码回归方程Y=4.77A+42.14B+18.61D−7.14AB−35.06AD+1.64BD−17.64A2−58.46B2−15.82D2,实际回归方程Y=126.41A+55.02B+129.54D−0.45AB−10.02AD+0.12BD−4.41A2−0.9B2−5.17D2。表11方差分析可知,模型P值=0.0015<0.05,差异极显著,失拟值P=0.6681>0.05,不存在失拟现象,拟合度好。模型R2=0.9422,说明该模型能反映响应值的变化。
表 11 蛋白酶活回归分析Table 11. Regression analysis of protease activity来源 平方和 自由度 均方 F值 P值 模型 36533.35 9 4059.26 12.68 0.0015 A 127.19 1 127.19 0.3974 0.5485 B 14203.71 1 14203.71 44.38 0.0003 D 1939.09 1 1939.09 6.06 0.0434 AB 203.92 1 203.92 0.6371 0.4510 AD 2648.05 1 2648.05 8.27 0.0238 BD 10.73 1 10.73 0.0335 0.8599 A2 1068.84 1 1068.84 3.34 0.1104 B2 11741.24 1 11741.24 36.68 0.0005 D2 859.92 1 859.92 2.69 0.1452 残差 2240.54 7 320.08 失拟值 333.81 2 166.91 0.4377 0.6681 纯误差 1906.72 5 总和 38773.88 16 蛋白酶活响应面曲面图结果表明(图6),乳糖和酵母浸粉有交互影响,乳糖和菜粕无交互作用,菜粕和酵母浸粉无显著影响。利用软件预测响应值最大值,模型预测各因素组合为乳糖5.55 g/L、菜粕31.89 g/L、酵母浸粉7.09 g/L时,预测到最大蛋白酶活为501 U/mL,通过验证实验得蛋白酶活为497±4.56 U/mL,达到了预测值的99%。
2.6 发酵条件对菌株产酶的影响
除了发酵培养基会显著影响菌株发酵过程中的蛋白酶活,温度、初始pH等发酵条件也会对菌株的产酶起到关键作用。
如图7A所示,在发酵温度梯度变化下,蛋白酶活出现先升高后下降的趋势,34 ℃时蛋白酶活达到最高,为482.91±6.65 U/mL,50 ℃时达到最低,则此蛋白酶耐热性低。由于蛋白酶作为活性蛋白,对环境因素敏感,易失活,会使蛋白酶活降低,适宜的发酵温度可使酶变得活跃,从而增加蛋白酶活的产量[32]。
接种量对菌株产酶的影响见图7B,当接种量为3.0%时蛋白酶活达到468.68±1.64 U/mL,随着接种量不断增加,蛋白酶活和接种量3.0%时的蛋白酶活无显著差异(P>0.05)。可能是增加接种量导致产酶时间提前,蛋白质等营养成分也被更快消耗,没有足够的蛋白诱导,使得蛋白酶活无明显变化。
不同初始pH对菌株ppr3产酶的影响显示(图7C),当培养基初始pH5~7时蛋白酶活几乎相同,分别为493.83±10.28、494.42±3.18 U/mL,pH9时蛋白酶活有显著下降趋势(P<0.05),但仍然保持较高酶活,有研究表明贝莱斯芽孢杆菌生长最适pH5[33]。
2.6.1 发酵条件的响应面优化设计
响应面设计获得了自变量发酵温度(a)、pH(b)、接种量(c)与因变量Y间的编码回归方程(表12)Y=−13.69a+83.79b+12c+0.6475ab+2.56ac−6.13bc−60.67a2−113.49b2−7.98c2,实际回归方程Y=62.16a+327.9b+14.23c+0.04ab+0.13ac−1.23bc−0.95a2−28.37b2−1.28c2。通过方差分析(表13),该模型P值<0.01,表明差异极显著,失拟值>0.05,失拟项不显著,说明不存在失拟现象,拟合度好,R2=0.9987,因此可以用此模型进行验证分析。根据响应面曲线及椭圆形曲线探讨温度、pH、接种量之间的交互作用,图8显示了温度与pH之间无交互作用,接种量和pH之间有交互影响,温度和接种量之间无交互作用。
表 12 Box-Behnken试验设计及结果Table 12. Box-Behnken experimental design and results试验号 a 发酵温度(℃) b pH c 接种量(%) Y 蛋白酶活(U/mL) 1 26 5 0.5 421.65±0.98f 2 34 3 0.5 257.68±2.01k 3 34 7 0.5 441.98±1.67e 4 42 5 0.5 386.69±2.11i 5 26 3 3 238.45±2.09i 6 26 7 3 400.28±3.01h 7 34 5 3 478.15±1.87c 8 34 5 3 486.57±2.88a 9 34 5 3 483.54±2.56ab 10 34 5 3 486.56±3.11a 11 34 5 3 480.47±3.27bc 12 42 3 3 216.23±2.34m 13 42 7 3 380.65±1.01i 14 26 5 5.5 441.02±2.27e 15 34 3 5.5 293.46±3.17j 16 34 7 5.5 453.25±4.07d 17 42 5 5.5 412.31±3.22g 表 13 蛋白酶活回归分析Table 13. Regression analysis of protease activity来源 平方和 自由度 均方 F值 P值 模型 1.33E+05 9 14811.18 578.55 <0.0001 a 1499.33 1 1499.33 58.87 0.0001 b 56169.46 1 56169.46 2194.07 <0.0001 c 1152.96 1 1152.96 45.04 0.0003 ab 1.68 1 1.68 0.0655 0.8054 ac 26.27 1 26.27 1.03 0.3448 bc 150.19 1 150.19 5.87 0.0460 a2 15495.88 1 15495.88 605.29 <0.0001 b2 54231.73 1 54231.73 2118.38 <0.0001 c2 267.81 1 267.81 10.46 0.0144 残差 179.2 7 25.6 失拟值 123.59 3 41.2 2.96 0.1607 纯误差 55.62 4 13.9 总和 1.335E+05 16 该模型预测发酵条件为发酵温度33 ℃,接种量为8.5%,初始pH5.7,蛋白酶活达到最大值502 U/mL,通过3组平行验证实验,蛋白酶活的平均值为500±5.86 U/mL,达到了预测值的98%,说明实际值与预测值有很好的拟合性。
3. 结论
本研究筛选鉴定了一株高产中性蛋白酶的贝莱斯芽孢杆菌,蛋白酶活(215.4±2.26 U/mL)经统计学优化,提高到500 U/mL,卢超等[34]、许曌昕等[35]筛选到的产中性蛋白酶菌株,蛋白酶活分别为401.83、192.7 U/mL,均低于本研究的酶活。同时贝莱斯芽孢杆菌具有潜在的益生作用[21],这使得其在食品和农业的应用具有明显优势。通过方差分析和二阶模型回归方程,表明最优发酵条件为发酵温度33 ℃,接种量为8.5%,初始pH5.7,培养基成分为乳糖5.55 g/L、菜粕31.89 g/L、酵母浸粉7.09 g/L、吐温80 0.03%、MgSO4 0.04%,K2HPO4 0.04%。而菜粕作为替代蛋白,其成本相较于宋立立等[36]、郭艳霞等[37]显著降低。使农业副产物得到了有效利用,但其工业生产有待扩大研究,其产酶机理、酶学性质及益生潜力也需要进一步探索。
-
图 5 菌种的产酶曲线及发酵培养基成分对酶活的影响
注:不同培养基成分的蛋白酶活(左轴),不同浓度的培养基成分蛋白酶活(右轴);b~e图中柱状图表示不同浓度的培养基成分蛋白酶活,折线图表示不同浓度的培养基成分蛋白酶活;图(f)左轴代表酵母浸粉蛋白酶活,右轴代表K2HPO4蛋白酶活;同一指标下不同字母表示差异显著(P<0.05),图7同。
Figure 5. Enzyme production curve of strain and effect of medium composition on enzyme activity
表 1 Plackett-Burman试验设计
Table 1 Plackett-Burman experimental design
因素 水平(g/L) −1 1 A 乳糖 0.5 4.5 B 菜粕 2 18 C MgSO4 0.4 3.5 D 酵母浸粉 0.25 3.75 E 吐温80 0.3 1.0 F K2HPO4 0.4 1.5 表 2 最陡爬坡试验设计
Table 2 Design of steepest climbing test
组数 成分 A 乳糖(g/L) B 菜粕(g/L) D 酵母浸粉(g/L) 1 0.5 2 0.25 2 2.5 10 2.00 3 4.5 18 3.75 4 6.5 26 5.50 5 8.5 34 7.25 6 10.5 42 9.00 7 12.5 50 10.75 表 3 培养基成分响应面优化设计
Table 3 Response surface optimization design of medium components
因素 水平(g/L) −1 0 1 A 乳糖 4.5 6.5 8.5 B 菜粕 18 26 34 D 酵母浸粉 3.75 5.50 7.25 表 4 发酵条件响应面优化设计
Table 4 Response surface optimization design of fermentation conditions
因素 水平 −1 0 1 a 发酵温度(℃) 26 34 42 b pH 3 5 7 c 接种量(%) 0.5 3 5.5 表 5 透明圈与菌株直径的比值及相应酶活
Table 5 Ratio of transparent circle to strain diameter and corresponding enzyme activities
表 6 不同pH下ppr3蛋白酶活
Table 6 ppr3 Protease activity at different pH
蛋白酶种类 pH ppr3蛋白酶活(U/mL) 酸性蛋白酶 3.0 138.17±6.77c 中性蛋白酶 7.5 220.40±3.47a 碱性蛋白酶 10.5 188.30±4.32b 表 7 Plackett-Burman试验设计及结果
Table 7 Plackett-Burman experimental design and results
试验号 A
乳糖
(g/L)B
菜粕
(g/L)C
MgSO4
(g/L)D
酵母浸粉
(g/L)E
吐温80
(g/L)F
K2HPO4
(g/L)Y 蛋白酶活
(U/mL)1 4.5 18 3.5 0.5 0.3 0.4 270±9.9d 2 0.5 18 0.4 3.5 1 0.4 234±4.24f 3 4.5 2 3.5 3.5 0.3 1.5 297±10.26b 4 4.5 18 0.4 0.5 0.3 1.5 287±3.45c 5 0.5 2 0.4 0.5 0.3 0.4 172±5.23h 6 4.5 18 0.4 3.5 1 1.5 322±2.55a 7 0.5 2 3.5 0.5 1 1.5 157±6.01i 8 0.5 18 3.5 3.5 0.3 0.4 264±4.66d 9 4.5 2 3.5 3.5 1 0.4 210±3.24g 10 0.5 18 3.5 0.5 1 1.5 254±5.77e 11 0.5 2 0.4 3.5 0.3 1.5 240±3.21f 12 4.5 2 0.4 0.5 1 0.4 232±2.28f 表 8 Plackett-Burman 设计方差分析
Table 8 Plackett-Burman design variance analysis
因素 F检验 P值 重要性 A 12.30 0.0172 ※※※※※ B 14.54 0.0125 ※※※※※※ C 0.1708 0.6966 ※ D 5.30 0.0696 ※※※※ E 2.04 0.2125 ※※ F 4.27 0.0937 ※※※ 注:※数量表示因素显著性,数量越多显著性越强。 表 9 最陡爬坡试验
Table 9 Steepest ascent experiment
试验号 A 乳糖(g/L) B 菜粕(g/L) D 酵母浸粉(g/L) 蛋白酶活(U/mL) 1 0.5 2 0.25 221.18±4.47f 2 2.5 10 2 363.02±1.33d 3 4.5 18 3.75 415.84±6.66c 4 6.5 26 5.5 483.79±6.65a 5 8.5 34 7.25 409.25±8.68c 6 10.5 42 9 258.31±2.67e 7 12.5 50 10.75 437.49±4.78b 表 10 Box-Behnken试验设计及结果
Table 10 Box-Behnken experimental design and results
试验号 A 乳糖(g/L) B 菜粕(g/L) D 酵母浸粉(g/L) Y 蛋白酶活(U/mL) 1 4.5 34 5.5 451.23±3.21f 2 6.5 34 7.25 476.56±2.89de 3 6.5 18 7.25 378.12±4.33m 4 6.5 26 5.5 491.58±1.98a 5 8.5 26 7.25 413.02±2.44h 6 8.5 18 5.5 382.13±7.02lm 7 6.5 34 3.75 441.32±4.77g 8 6.5 26 5.5 486.53±5.03abc 9 4.5 26 3.75 388.49±3.01i 10 6.5 26 5.5 479.52±6.02cde 11 6.5 26 5.5 488.56±4.67ab 12 4.5 18 5.5 363.57±3.58n 13 6.5 18 3.75 349.43±3.28o 14 8.5 26 3.75 473.40±2.65e 15 6.5 26 5.5 482.00±4.33bcd 16 8.5 26 7.25 473.21±1.98e 17 8.5 34 5.5 441.23±5.67g 表 11 蛋白酶活回归分析
Table 11 Regression analysis of protease activity
来源 平方和 自由度 均方 F值 P值 模型 36533.35 9 4059.26 12.68 0.0015 A 127.19 1 127.19 0.3974 0.5485 B 14203.71 1 14203.71 44.38 0.0003 D 1939.09 1 1939.09 6.06 0.0434 AB 203.92 1 203.92 0.6371 0.4510 AD 2648.05 1 2648.05 8.27 0.0238 BD 10.73 1 10.73 0.0335 0.8599 A2 1068.84 1 1068.84 3.34 0.1104 B2 11741.24 1 11741.24 36.68 0.0005 D2 859.92 1 859.92 2.69 0.1452 残差 2240.54 7 320.08 失拟值 333.81 2 166.91 0.4377 0.6681 纯误差 1906.72 5 总和 38773.88 16 表 12 Box-Behnken试验设计及结果
Table 12 Box-Behnken experimental design and results
试验号 a 发酵温度(℃) b pH c 接种量(%) Y 蛋白酶活(U/mL) 1 26 5 0.5 421.65±0.98f 2 34 3 0.5 257.68±2.01k 3 34 7 0.5 441.98±1.67e 4 42 5 0.5 386.69±2.11i 5 26 3 3 238.45±2.09i 6 26 7 3 400.28±3.01h 7 34 5 3 478.15±1.87c 8 34 5 3 486.57±2.88a 9 34 5 3 483.54±2.56ab 10 34 5 3 486.56±3.11a 11 34 5 3 480.47±3.27bc 12 42 3 3 216.23±2.34m 13 42 7 3 380.65±1.01i 14 26 5 5.5 441.02±2.27e 15 34 3 5.5 293.46±3.17j 16 34 7 5.5 453.25±4.07d 17 42 5 5.5 412.31±3.22g 表 13 蛋白酶活回归分析
Table 13 Regression analysis of protease activity
来源 平方和 自由度 均方 F值 P值 模型 1.33E+05 9 14811.18 578.55 <0.0001 a 1499.33 1 1499.33 58.87 0.0001 b 56169.46 1 56169.46 2194.07 <0.0001 c 1152.96 1 1152.96 45.04 0.0003 ab 1.68 1 1.68 0.0655 0.8054 ac 26.27 1 26.27 1.03 0.3448 bc 150.19 1 150.19 5.87 0.0460 a2 15495.88 1 15495.88 605.29 <0.0001 b2 54231.73 1 54231.73 2118.38 <0.0001 c2 267.81 1 267.81 10.46 0.0144 残差 179.2 7 25.6 失拟值 123.59 3 41.2 2.96 0.1607 纯误差 55.62 4 13.9 总和 1.335E+05 16 -
[1] JAHAN-MIHAN A, LUHOVYY B L, KHOURY D E, et al. Dietary proteins as determinants of metabolic and physiologic functions of the gastrointestinal tract[J]. Nutrients,2011,3(5):574−603. doi: 10.3390/nu3050574
[2] MATTHEWS D M, CRAMPTON R F, LIS M T. Intestinal absorption of peptides[J]. Biochemical Society Transactions,1983,292(7568):639−640.
[3] CHEN P, OU Y, LI S. Effect of different protein sources on the concentrations of small peptides in the rumen of sheep[J]. Archives of Animal Nutrition,2007,61(4):301−307. doi: 10.1080/17450390701432696
[4] KUMAR M S. Peptides and peptidomimetics as potential antiobesity agents: Overview of current status[J]. Frontiers in Nutrition,2019,6:11. doi: 10.3389/fnut.2019.00011
[5] HAMEDY R J, FITZGERALD R J. Bioactive peptides from marine processing waste and shellfish: A review[J]. Journal of Functional Foods,2012,4(1):6−24. doi: 10.1016/j.jff.2011.09.001
[6] RUNNI M, RUNU C, DUTTA A. Role of fermentation in improving nutritional quality of soybean meal: A review[J]. Asian-Australasian Journal of Animal Sciences,2016,29(11):1523−1529.
[7] CASSONE M, JR L O. Synergy among antibacterial peptides and between peptides and small-molecule antibiotics[J]. Expert Rev Anti Infect Ther,2010,8(6):703−716. doi: 10.1586/eri.10.38
[8] ZHANG J N, ZHAO M M, SU G W, et al. Identification and taste characteristics of novel umami and umami-enhancing peptides separated from peanut protein isolate hydrolysate by consecutive chromatography and UPLC-ESI-QTOF-MS/MS[J]. Food Chemistry,2019,25(278):674−682.
[9] 王芳, 贾万利, 张浩男, 等. 混合菌发酵对豆粕品质的影响[J]. 甘肃农业大学学报,2017,52(4):45−51. [WANG F, JIA W L, ZHANG H N, et al. Effect of mixed bacteria fermentation on quality of soybean meal[J]. Journal of Gansu Agricultural University,2017,52(4):45−51. doi: 10.13432/j.cnki.jgsau.2017.04.007 [10] 毛银, 陆春波, 李国辉, 等. 菌酶协同发酵豆粕工艺的优化[J]. 食品与发酵工业,2019,45(14):108−114. [MAO Y, LU C B, LI G H, et al. Optimized fermentation process of soybean meal by bacteria and enzymes[J]. Food and Fermentation Industries,2019,45(14):108−114. doi: 10.13995/j.cnki.11-1802/ts.019892 [11] 汤恒, 张娟, 堵国成. 高效降解牛奶过敏原蛋白酶菌种筛选及其重组植物过氧化氢酶的功能验证[J]. 食品与生物技术学报,2021,40(3):28−34. [TANG H, ZHANG J, DU G C. Study on milk allergen protease screening and its recombinant vegetative catalase function verificition[J]. Journal of Food Science and Bictechnology,2021,40(3):28−34. doi: 10.3969/j.issn.1673-1689.2021.03.004 [12] WANG J, XU A, WAN Y, et al. Purification and characterization of a new metallo-neutral protease for beer brewing from Bacillus amyloliquefaciens SYB-001[J]. Applied Biochemistry & Biotechnology,2013,170(8):2021−2033.
[13] 唐霄, 孙杨赢, 江雪婷, 等. 不同蛋白酶制备鹅肉呈味肽的对比分析[J]. 食品科学,2019,40(22):141−146. [TANG X, SUN Y Y, JIANG X T, et al. Comparative analysis of flavor peptides prepared by enzymatic hydrolysis of goose meat with different proteases[J]. Food Science,2019,40(22):141−146. doi: 10.7506/spkx1002-6630-20181023-272 [14] 王亮亮, 李方方, 姜锡瑞. 我国蛋白酶的产业近况与展望[J]. 中国食品添加剂,2021,32(8):141−150. [WANG L L, LI F F, JIANG X R. Current situation and prospect of protease industry in China[J]. China Food Additives,2021,32(8):141−150. doi: 10.19804/j.issn1006-2513.2021.08.020 [15] 吕凯波, 朱文婷. 双酶分步酶解红花籽粕制备抗氧化肽[J]. 食品工业,2020,41(10):24−26. [LÜ K B, ZHU W T. The process of safflower seed meal and antioxidant activity by double enzymatic hydrolysis[J]. The Food Industry,2020,41(10):24−26. [16] 王玉敏, 齐红丽. 酶解大豆粕蛋白制备大豆低分子肽的研究[J]. 化学与生物工程,2010,27(4):64−66. [WANG Y M, QI H L. Study on preparation of low molecular weight peptides by enzymolysis of soybean protein[J]. Chemistry Bioengineering,2010,27(4):64−66. doi: 10.3969/j.issn.1672-5425.2010.04.018 [17] 王晓云. 高产蛋白酶枯草芽孢杆菌的筛选与诱变选育研究[D]. 泰安: 山东农业大学, 2015 WANG X Y. Screening and mutation breeding of Bacillus subtilis with high protease production[D]. Taian: Shandong Agricultural University, 2015.
[18] 周晶. 产低温蛋白酶动性球菌的筛选及其在低盐鱼露发酵中的应用[D]. 镇江: 江苏大学, 2020 ZHOU J. Screening of planococcus bacteria low-temperature proteases and its application in low-salt fish sauce fermentation[D]. Zhenjiang: Jiangsu University, 2020.
[19] 殷方荣. 贝莱斯芽孢杆菌角蛋白酶的基因克隆、表达、酶学性质及应用[D]. 无锡: 江南大学, 2021 YIN F R. Gene cloning, expression, characteristics and application of keratinase from Bacillus velezensis[D]. Wuxi: Jiangnan University, 2021.
[20] KHALID F, KHALID A, FU Y, et al. Potential of Bacillus velezensis as a probiotic in animal feed: A review[J]. The Journal of Microbiology,2021,59(7):627−633. doi: 10.1007/s12275-021-1161-1
[21] IBRAHIM, ABDELNASSER S, et al. Optimization of media and cultivation conditions for alkaline protease production by Alkaliphilic bacillus halodurans[J]. Current Research in Bacteriology,2009,4(7):251−259.
[22] 朱祥杰, 王震, 苑志欣, 等. 海洋芽孢杆菌N11-8产蛋白酶的发酵条件优化[J]. 渔业科学进展,2018,39(6):155−163. [ZHU X J, WANG Z, YUAN Z X, et al. Optimization of fermentation conditions of Bacillus sp N11-8 on the production of protease PBN11-8[J]. Progress in Fishery Sciences,2018,39(6):155−163. [23] 王永红, 李小斌, 徐磊, 等. 产蛋白酶菌株的筛选、鉴定及水解菜粕蛋白能力[J]. 生物资源,2018,40(2):135−140. [WANG Y H, LI X B, XU L, et al. Screening and identification of protease producing strains and their ability of hydrolysis of rapeseed protein[J]. Biological Resources,2018,40(2):135−140. [24] 杜静, 钮琰星, 周琦, 等. 固相微萃取条件优化及发酵菜粕风味物质分析[J]. 中国粮油学报,2017,32(7):114−120. [DU J, NIU Y X, ZHOU Q, et al. Optimization of solid phase microextraction (SPME) and analysis of volatile compounds in fermented rapeseed meal[J]. Chinese Journal of Cereals and Oils,2017,32(7):114−120. doi: 10.3969/j.issn.1003-0174.2017.07.019 [25] 马淑惠, 刘东彦. K+、Na+、Ca2+、Mg2+在生物过程中的作用[J]. 固原师专学报,2003(3):85−87. [MANG S H, LIU D Y. The roles of K+, Na+, Ca2+, Mg2+ in biological processes[J]. Journal of Guyuan Teachers College,2003(3):85−87. [26] ANANDHARAJ M, SIVASANKARI B, SIDDHARTHAN N, et al. Production, purification, and biochemical characterization of thermostable metallo-protease from novel Bacillus alkalitelluris TWI3 isolated from tannery waste[J]. Applied Biochemistry & Biotechnology,2016,178(8):1666−1686.
[27] DEIVE F J, CARVALHO E, PASTRANA L, et al. Strategies for improving extracellular lipolytic enzyme production by Thermus thermophilus HB27[J]. Bioresource Technology,2009,100(14):3630−3637. doi: 10.1016/j.biortech.2009.02.053
[28] SILVA W, MITIDIERI S, SCHRANK A, et al. Production and extraction of an extracellular lipase from the entomopathogenic fungus metarhiziumanisopliae[J]. Process Biochemistry,2005,40(1):321−326. doi: 10.1016/j.procbio.2004.01.005
[29] JACBO F F, STRIEGEL L, RVCHLIK M, et al. Spent yeast from brewing processes: A biodiverse starting material for yeast extract production[J]. Fermentation,2019,5(2):51. doi: 10.3390/fermentation5020051
[30] 贾天梅. 产耐热型碱性蛋白酶菌株筛选、发酵条件优化及其酶学性质研究[D]. 杭州: 浙江工商大学, 2018 JIA T M. Study on screening of thermostabilable alkaline protease strain, optimization of fermentation condition and its emzymatic properties[D]. Hangzhou: Zhejiang Gongshang University, 2018.
[31] LIU S, FANG Y, LYU M, et al. Optimization of the production of organic solvent-stable protease by Bacillus sphaericus DS11 with response surface methodology[J]. Bioresource Technology,2010,101(20):7924−7929. doi: 10.1016/j.biortech.2010.05.057
[32] LI X Q, CHAI X Q, LIU D Y, et al. Effects of temperature and feed processing on protease activity and dietary protease on growths of white shrimp, Litopenaeus vannamei, and tilapia, Oreochromis niloticus×O.aureus[J]. Aquaculture Nutrition,2016,22(6):1283−1292. doi: 10.1111/anu.12330
[33] 陈偿, 丁雄祺, 谢媚, 等. 一株具有高效降解水产养殖水体中饲料淀粉功能的贝莱斯芽孢杆菌D1及其应用: CN111235065A[P]. 2020-06-05 CHEN C, DING X Q, XIE M, et al. A Bacillus cereus D1 with efficient degradation of feed starch in aquaculture water and its application: CN111235065A[P]. 2020-06-05.
[34] 卢超, 陈景鲜, 王国霞, 等. 枯草芽孢杆菌L07产中性蛋白酶发酵条件优化[J]. 食品与发酵工业,2020,46(16):148−153. [LU C, CHEN J X, WANG G X. et al. Optimization of fermentation conditions for the production of neutral protease by Bacillus subtilis L07[J]. Food and Fermentation Industries,2020,46(16):148−153. doi: 10.13995/j.cnki.11-1802/ts.023566 [35] 许曌昕, 曾辉, 曾德勇, 等. 响应面法优化解淀粉芽孢杆菌产中性蛋白酶发酵条件[J]. 中国酿造,2017,36(4):78−82. [XU Z X, ZENG H, ZENG D Y. et al. Optimization of fermentation conditions for neutral protease-production Bacillus amyloliquefaciens by response surface methodology[J]. China Brewing,2017,36(4):78−82. doi: 10.11882/j.issn.0254-5071.2017.04.017 [36] 宋立立, 张泽生. 枯草芽孢杆菌产中性蛋白酶发酵条件的优化[J]. 饲料研究,2020,43(4):66−70. [SONG L L, ZHANG Z S. Optimization of fermentation conditions for the production of neutral protease by Bacillus subtilis[J]. Feed Research,2020,43(4):66−70. doi: 10.13557/j.cnki.issn1002-2813.2020.04.018 [37] 郭艳霞, 贾丽艳, 畅盼盼, 等. 响应面法优化贝莱斯芽孢杆菌YB19产中性蛋白酶发酵条件[J]. 中国调味品,2021,46(4):16−20, 37. [GUO Y X, JIA L Y, CHANG P P, et al. Optimization of fermentation conditions for neutral protease produced by Bacillus cereus YB19 by response surface methodology[J]. China Condiment,2021,46(4):16−20, 37. doi: 10.3969/j.issn.1000-9973.2021.04.004