自選釀酒酵母（KDLYS9-3）產β-D-葡萄糖苷酶動力學研究

李艷

（1.河北科技大學生物科學與工程學院，河北石家莊050018；2.河北省發(fā)酵工程技術研究中心，河北石家莊050018）

自選釀酒酵母（KDLYS9-3）產β-D-葡萄糖苷酶動力學研究

李艷1，2

（1.河北科技大學生物科學與工程學院，河北石家莊050018；2.河北省發(fā)酵工程技術研究中心，河北石家莊050018）

自選高產β-D-葡萄糖苷酶的釀酒酵母（KDLYS9-3）在葡萄酒釀造過程中具有增強香氣的效果。依據菌體生長和產酶試驗，利用Logistic方程、DoseResp方程和Nelder方程建立了菌體生長和產酶，以及菌體生長速率與酶生成速率之間關系的動力學模型，通過Origin8.0軟件進行非線性擬合，并利用Lineweaver-Burk法作圖測定了該菌所產β-D-葡萄糖苷酶的動力學參數Km值和Vmax值。結果顯示：自選釀酒酵母KDLYS9-3的菌體生長與產酶的相關性為部分偶聯(lián)型，動力學模型與試驗值吻合度好，方程能夠反映菌體生長與產酶的變化規(guī)律。菌體生長8 h后開始產酶，菌體進入對數生長期時酶大量生成，到39 h菌體生長進入穩(wěn)定期，隨著菌體生長進入衰亡期后酶也隨之停止產生；利用Lineweaver-Burk法求得酶動力學參數Km=8.492 579 mmol/L，Vmax=1.030 715（μmol/L）/min。研究結果為該菌株的理論研究和實際應用奠定了基礎。

自選釀酒酵母（KDLYS9-3）；β-D-葡萄糖苷酶；動力學方程；米氏常數（Km）；最大反應速率（Vmax）

β-D-葡萄糖苷酶在葡萄酒生產過程中可增加葡萄酒的香氣[1-3]。目前，所用商品β-D-葡萄糖苷酶大多來自于黑曲霉菌發(fā)酵生產得到的酶制劑[4]，該酶的加入改善葡萄酒香氣的同時也給葡萄酒留下造成蛋白質不穩(wěn)定的隱患[5]。因此，篩選可以產生β-D-葡萄糖苷酶的釀酒酵母，在進行酒精發(fā)酵的同時提高香氣物質含量已成為葡萄酒行業(yè)關注的重要課題之一。

本研究是對本實驗室篩選得到的一株高產β-D-葡萄糖苷酶的釀酒酵母（KDLYS9-3）[6]的菌體生長、產酶以及彼此間的關系在實驗的基礎上，建立動力學模型，探索細胞生長對產酶的影響，并對該菌所產β-D-葡萄糖苷酶的動力學參數Km和Vmax進行測定，為理論研究和實際應用奠定基礎。

1 材料與方法

1.1 材料

1.1.1 釀酒酵母

自選釀酒酵母KDLYS9-3：本研究室于2010年采摘沙城產區(qū)夾河、沙營和東水泉3個鄉(xiāng)的龍眼葡萄，帶回實驗室后分別進行自然發(fā)酵，在自然發(fā)酵后期分離得到的一株酵母菌，對該菌株進行了WL培養(yǎng)基培養(yǎng)的菌落和顯微形態(tài)分類，結合5.8S rRNA-ITS區(qū)域及核糖體26S D1/D2區(qū)域RFLP分析的分子鑒定和基因測序證明是一株釀酒酵母[6]，后經王彩肖等[7]鑒定為高產β-D-葡萄糖苷酶的菌株。該菌株Genebank登錄號為：KX098507，核糖體26S D1/D2區(qū)域限制片長多態(tài)性（restriction fragment length polymorphism，RFLP）分析的基因序列為：

AGGGGCATTGCCTTAGTACGGCGAGTGAGCGGC AAAAGCTCAAATTTGAAATCTGGTACCTTCGGTGCCC GCGTTGTAATTTGGAGAGGGCAACTTTGGGGCCGTTC CTTGTCTATGTTCCTTGGAACAGGACGTCATAGAGGG TGAGAATCCCGTGTGGCGAGGAGTGCGGTTCTTTGTA AAGTGCCTTCGAAGAGTCGAGTTGTTTGGGAATGCA GCTCTAAGTGGGTGGTAAATTCCATCTAAAGCTAAAT ATTGGCGAGAGACCGATAGCGAACAAGTACAGTGAT GGAAAGATGAAAAGAACTTTGAAAAGAGAGTGAAAAAGTACGTGAAATTGTTGAAAGGGAAGGGCATTTGATCAGACATGGTGTTTTGTGCCCTCTGCTCCTTGTGGGTAGGGGAATCTCGACTTTCACTGGGCCAGCATCAGTTTTGGTGGCAGGATAAATCCATAGGAATGTAGCTTGCCTCGGTAAGTATTATAGCCTGTGGGAATACTGCCAGCTGGGACTGAGGACTGCGACGTAAGTCAAGGATGCTGGCATAATGGTTATATGCCGCCCGTCTTGAACCACGGACCA。

商品對照菌（安琪葡萄酒酵母）：中國湖北安琪活性干酵母股份有限公司。

培養(yǎng)基：斜面保藏培養(yǎng)基使用固態(tài)YEPD培養(yǎng)基；發(fā)酵培養(yǎng)基為液態(tài)YEPD培養(yǎng)基。

1.1.2 試劑

對硝基苯基-β-D-葡萄糖苷（pNPG）：上海寶曼生物科技有限公司；對硝基苯酚（pNP）：天津市博迪化工有限公司；磷酸二氫鈉、檸檬酸和無水碳酸鈉：天津永大化學試劑有限公司；蛋白胨、葡萄糖、酵母浸粉和瓊脂：北京澳博星生物技術有限責任公司。

1.1.3 儀器

Multiskan FC酶標儀：Thermo Fisher Scientific公司；DK-98-1型電熱恒溫水浴鍋：天津市泰斯特儀器有限公司；Jnoec-XS-212202型生物顯微鏡：南京雄姿光電儀器有限公司；高速離心機LG16-C：北京雷博爾離心機有限公司；ZQTY-70型恒溫振蕩培養(yǎng)箱：上海知楚儀器有限公司。

1.2 方法

1.2.1 酶活力的測定

粗酶液的獲得：挑取保藏斜面上的單菌落，接種到YEPD液態(tài)培養(yǎng)基中進行菌種的活化培養(yǎng)，當活菌數達到108～109個/mL時，以106個/mL的接種量接種到液態(tài)YEPD培養(yǎng)基中進行發(fā)酵產酶培養(yǎng)，培養(yǎng)條件為28℃、190 r/min振蕩培養(yǎng)48 h。培養(yǎng)結束后將菌懸液在4 000 r/min下離心15 min，取上清液視為粗酶液待用。

酶活力的測定參考王彩肖等[7]方法。

β-D-葡萄糖苷酶的酶活力單位定義：50℃，pH 5.0時，1 min內β-D-葡萄糖苷酶水解底物pNPG產生1 μmol的pNP所需的酶量。

1.2.2 酵母菌生長曲線的測定

測定酵母生長曲線是在發(fā)酵期間定時取樣，用無菌水適當稀釋，生物顯微鏡下血球計數板計數，得到菌體的濃度，計算公式為：x=50 000nb

式中：x為總菌數（個/mL）；50 000為生物顯微鏡的自身系數；n為顯微鏡中上下計數室菌數的平均值；b為稀釋倍數，試驗做3次平行。

1.2.3 動力學方程及參數的測定

Logistic方程[8]是一種典型的“S”型的模型方程[9-10]，最初被用于人口增長趨勢的統(tǒng)計，該模型和菌體生長具有很大的相似性，所以利用Logistic方程建立菌體生長模型；利用DoseResp方程建立產酶的動力學模型；菌體生長與產酶之間的關系模型采用Nelder數學方程建立。通過Excel 2003和Origin 8.0軟件進行試驗值與方程計算值的線性和非線性擬合，并在模型的基礎上進行試驗驗證。

分別配制底物濃度為 5、10、20、30、35 mmol/L，在其他條件相同的情況下測定酶活力，利用Lineweaver-Burk法（簡稱L-B法，也稱雙倒數法）作圖求得自選釀酒酵母KDLYS9-3所產β-D-葡萄糖苷酶的動力學參數Km值和Vmax值。

2 結果與分析

2.1 自選釀酒酵母KDLYS9-3生長與產酶的測定

發(fā)酵過程中每間隔一定時間取樣分別測定菌體總數和酶活力，以時間為橫坐標，菌體總數和β-D-葡萄糖苷酶活力分別為縱坐標作圖得到圖1。

圖1 KDLYS9-3產酶與菌數曲線圖Fig.1 The curve of enzyme and colonies production for KDYLS9-3

可以看出KDLYS9-3存在4 h～8 h的延滯期，在延滯期內菌體適應新的生長環(huán)境，幾乎檢測不到β-D-葡萄糖苷酶的活力，說明菌種在這個時期產酶較少；8 h～39 h為菌體生長的對數期，菌體數量迅速增加。39 h～48 h為菌種生長的穩(wěn)定期。而菌體產酶活力在8 h～42 h快速提高，42 h～48 h期間維持較高活性，說明在菌體生長的對數期和穩(wěn)定期是菌體產酶的最佳時期；在48 h后菌體生長進入衰亡期，菌體產酶量也較低。從菌數和酶活力的關系可以得出結論，自選釀酒酵母KDLYS9-3產β-D-葡萄糖苷酶屬于對數和穩(wěn)定期型，即在菌體生長的初期和衰亡期不會產生大量的酶，對數期和穩(wěn)定期為產酶的最佳時期，說明菌體生長與產酶具有半偶聯(lián)的關系。

2.2 自選釀酒酵母KDLYS9-3菌體生長動力學模型的建立與擬合

由圖1可以看出菌體的生長符合“S”型曲線，Logistic方程[11]被廣泛地用于描述菌體生長的模型，通過不依賴于基質濃度的“S”形曲線，簡單計算得到具有生物學和幾何意義的發(fā)酵動力學參數[12-13]，說明本研究可以利用Logistic方程建立數學模型，基本公式為：

式中：X為菌體的總菌數，個/mL；Xm為最大細胞濃度，個/mL；μm為最大比生長速率，h-1。由于在 0～8 h時菌體數量基本不增加，所以對本方程進行改進為在t＞8 h進行積分，以 t=0 時 X=X0為初始條件，對（1）式積分得：

由式2可以整理為公式3：

由圖1中的生長曲線和試驗數據可知，Xm為最大的菌體濃度1.50×109個/mL，可將公式3看做為：

圖2 直線回歸方程Fig.2 Linear regression equation

由圖2中的線性回歸方程的直線斜率K=μm=0.241 2，又由截距可以求得 x0為 2.55×106，xm=1.50×109，μm帶入式2中最終得到KDLYS9-3的生長動力學模型為：

①當發(fā)酵時間在t=0時，方程為

②當發(fā)酵時間為t＞0 h

由模型計算所得曲線與試驗所得數據曲線的對比為圖3。

由圖3和圖4可以看出，所建模型與試驗所測得的數據具有很高的相似性，各數據點的殘差較小，說明模型和實測值擬合較好，所建模型可以很好的描述菌種細胞的增長趨勢，既菌體經過一段時間的延滯期，之后快速的增長，在39 h之后達到穩(wěn)定期，之后菌體受環(huán)境的影響，使之在一定的時間內保持穩(wěn)定。通過模型與試驗的擬合，可以說明所建模型可以應用到實際的發(fā)酵操作中。

圖3 菌體生長動力學模型預測值與實驗值對比圖Fig.3 Comparison of predicted values and experimental values of the growth kinetics model of the yeast

圖4 Origin8.0擬合各數據點殘差圖Fig.4 The fitting of residual for each data point in Origin8.0

2.3 產酶活力動力學模型的建立與擬合

由圖1可以看出，酶活曲線呈現類似“S”型，所以本研究使用DoseResp方程進行擬合。DoseResp方程是一種“S”型函數方程，被使用在類菌體生長、受菌體生長的產物生成、基質消耗等[14-15]，本研究對產酶活力的動力學建立模型。DoseResp方程公式為：

式中：P(t)為酶活，U/mL；A1、A2、LOGx0和 C 為待定常數，t為生長時間，h。

由實時酶活力測定值P和時間t的對應關系，通過軟件Origin 8.0軟件進行擬合，可得到常數A1、A2、LOGx0和C分別如表1。

表1 DoseResp方程擬合相關參數Table 1 Fitting parameters of DoseResp equation

由表1可以看出，常數的標準誤差都小于0.5，說明擬合效果較好，P＜0.05，P 值較小，R2=0.998 74，說明擬合模型極顯著，可以反映試驗結果。最終得到產酶酶活數學模型為：

由模型的預測值與試驗所得值進行比較，做圖5和圖6。

圖5 酶活模型預測值與實驗值對比圖Fig.5 Comparison of predictive value and experimental value of enzyme activity model

圖6 Origin8.0擬合各數據點殘差圖Fig.6 The fitting of residual for each data point in Origin8.0

由圖5和圖6可以看出模型預測值與實驗值吻合度很好，各個數據點與預測數據點之間的殘差都在-0.000 8～0.000 8之間，說明了所建的方程模型可以很好的反映酶活力與發(fā)酵時間之間的關系，可以為今后利用自選釀酒酵母KDLYS9-3進行釀酒試驗提供數據基礎。該模型顯示了很好的擬合效果，且具有快速和簡便的優(yōu)點，可以很準確的反應試驗結果。

2.4 菌體生長與產酶相關性模型的建立

由圖1可以看出，β-D-葡萄糖苷酶是在自選釀酒酵母KDLYS9-3的生長過程中產生的，菌體生長和產酶具有部分偶聯(lián)關系，本研究利用Nelder方程建立模型，并利用Origin軟件進行擬合。基礎方程式為公式9：

式中：P（x）為反應酶活（U/mL）；a、b0、b1和 b2為待定常數；x為菌體數量（106個/mLl）。

經過軟件Origin軟件擬合方程參數為表2。

表2 Nelder方程擬合參數Table 2 Fitting parameters of Nelder equation

由表2可以看出，常數的標準誤差都較小，說明擬合效果較好，P＜0.001，P 值較小，R2=0.984，說明擬合模型極顯著，模型預測結果可以很好的反映實驗結果。最終得到菌體生長與產酶的關系數學模型為：

由模型的預測值與試驗所得值進行比較，做圖7和圖8。

圖7 總菌數與酶活的關系Fig.7 The relationship between total colonies and enzyme activity

圖8 Origin8.0擬合各數據點殘差圖Fig.8 The fitting of residual for each data point in Origin8.0

由圖7和圖8可以看出，試驗值與模型預測值吻合較好，說明所建模型符合試驗的真實性，可以很好的對自選釀酒酵母KDLYS9-3的生長與β-D-葡萄糖苷酶的生成之間的關系進行描述，由殘差圖可以看出各數據點的殘差在-0.003～0.003之間，也說明了模型可以真實的反映試驗的真實值，可以用公式10描述菌體生長與產酶之間的半偶聯(lián)關系。

2.5 KDLYS9-3產β-D-葡萄糖苷酶的動力學參數Km和Vmax

利用米氏方程[16]和L-B作圖法對自選釀酒酵母KDLYS9-3產β-D-葡萄糖苷酶的動力學參數Km值和Vmax值進行測定，以商品安琪釀酒酵母為對照菌，分別以底物濃度為橫坐標和酶活力為縱坐標做圖9。

圖9 底物濃度與酶活的關系Fig.9 The relationship between substrate concentration and enzyme activity

由圖9可以看出，底物濃度與對酶活力的影響，在底物濃度5 mmol/L～35 mmol/L范圍內，酶活力隨底物濃度的升高而逐漸增大。以米氏方程11和12可建立1/S和1/V的一元線性關系；

以1/S為橫坐標和1/V為縱坐標做圖10。

圖10 兩種酵母菌所產酶1/V和1/S關系圖Fig.10 The relationship of 1/V and 1/S in the enzyme from the two yeasts

由圖10可以看出KDLYS9-3和安琪酵母所產酶的1/S和1/V的線性擬合方程分別為y=8.239 5x+0.970 2和y=17.671x+1.8 831，線性相關性系數R2分別為0.998 8和0.994 7，擬合直線與試驗值吻合度很高，說明擬合效果很好。由橫縱截距可以求出KDLYS9-3的酶動力學參數為Km=8.492 579 mmol/L，Vmax=1.030 715（μmol/L）/min；安琪酵母的 β-D-葡萄糖苷酶動力學參數為Km=9.383 994 mmol/L，Vmax=0.531 039（μmol/L）/min。

Km可以反映底物與酶的親和力大小，Km越小，說明酶與底物的親和力越強，反之則越弱[17-18]，由兩種菌的酶動力學參數Km和Vmax可以看出：KDLYS9-3的Km值小于安琪酵母的Km值，說明KDLYS9-3所產β-D-葡萄糖苷酶對pNPG底物的反應親和力要強于商品菌安琪釀酒酵母。

3 結論

本文應用Logistic方程、DoseResp方程和Nelder方程分別對自選高產β-D-葡萄糖苷酶的釀酒酵母KDLYS9-3生長、產酶、生長與產酶之間的半偶聯(lián)關系進行了數學描述，建立了數學模型，結果顯示出所建數學模型與實驗值吻合性很好，證明了生長和產酶存在的部分偶聯(lián)關系。對自選菌和商品菌所產β-D-葡萄糖苷酶的動力學參數Km和Vmax值進行測定，由此證明了KDLYS9-3所產β-D-葡萄糖苷酶對pNPG的底物親和力強于對照菌，該菌株具有很好的應用潛質。

[1]Gamero A,Hernández-Orte P,Querol A,et al.Effect of aromatic precursor addition to wine fermentations carried out with different Saccharomyces species and their hybrids[J].International journal of food microbiology,2011,147(1):33-44

[2]González-Pombo P,Fari?a L,Carrau F,et al.Aroma enhancement in wines using co-immobilized Aspergillus niger glycosidases[J].Food chemistry,2014,143:185-191

[3]ünal M ü,Aksoy V A,?ener A.Isolation,purification and determination of some biochemical properties of β-glucosidase from Muscat of Bornova grape[J].European Food Research and Technology,2014,238(1):9-15

[4]Rodríguez M E,Lopes C,Valles S,et al.Selection and preliminary characterization of β-glycosidases producer Patagonian wild yeasts[J].Enzyme and Microbial Technology,2007,41(6):812-820

[5]Ferreira R B,Pi?arra-Pereira M A,Monteiro S,et al.The wine proteins[J].Trends in food science&technology,2001,12(7):230-239

[6]趙靜靜,李艷.沙城產區(qū)釀酒酵母多樣性研究[J].食品科學,2012,33(5):224-228

[7]王彩肖,武偉偉,李艷.產β-D-葡萄糖苷酶酵母菌的篩選及產酶性質研究[J].釀酒科技,2014(10):14-18

[8]Komali A S,Zheng Z,Shetty K.A mathematical model for the growth kinetics and synthesis of phenolics in oregano(Origanum vulgare)shoot cultures inoculated with Pseudomonas species[J].Process Biochemistry,1999,35(3):227-235

[9]鄭毅,張志國,關雄.蘇云金芽孢桿菌蛋白酶發(fā)酵動力學模型的構建[J].生物數學學報,2008,23(4):727-734

[10]劉建忠,翁麗萍.黑曲霉過氧化氫酶發(fā)酵過程的數學模型[J].工業(yè)微生物,2002,32(3):6-9

[11]Gabriel J P,Saucy F,Bersier L F.Paradoxes in the logistic equation[J].Ecological Modelling,2005,185(1):147-151

[12]張凱麗,鄭晗青,牛啟啟,等.釀酒酵母工程菌分批發(fā)酵產UMP動力學模型[J].食品工業(yè)科技,2015,36(9):158-161

[13]Mahajan P M,Gokhale S V,Lele S S.Production of nattokinase using Bacillus natto NRRL 3666:media optimization,scale up,and kinetic modeling[J].Food Science and Biotechnology,2010,19(6):1593-1603

[14]李俠,彭丹,張春暉,等.紅棗酒發(fā)酵動力學研究[J].中國釀造,2012,31(1):190-193

[15]王虎虎,董洋,徐幸蓮,等.酸化亞氯酸鈉對雞肉致病菌的致死模型研究[J].中國食品學報,2013(9):33-38

[16]Epstein E,Hagen C E.A kinetic study of the absorption of alkali cations by barley roots[J].Plant physiology,1952,27(3):457

[17]Jacobo-Velázquez D A,Hernández-Brenes C,Cisneros-Zevallos L,et al.Partial purification and enzymatic characterization of avocado(Persea americana Mill,cv.Hass)lipoxygenase[J].Food research international,2010,43(4):1079-1085

[18]Anese M,Sovrano S.Kinetics of thermal inactivation of tomato lipoxygenase[J].Food chemistry,2006,95(1):131-137

Research on the Enzyme Kinetics of Poducing Beta-D-Glucosidase from Indigenous Saccharomyces cerevisiae KDLYS9-3

LI Yan1，2
（1.College of Bioscience and Bioengineering，Hebei University of Science and Technology，Shijiazhuang 050018，Hebei，China；2.Hebei R&D Center for Fermentation Engineering，Shijiazhuang 050018，Hebei，China）

Indigenous Saccharomyces cerevisiae KDLYS9-3 produces high activity beta-D-glucosidase，and it has the effect of enhancing aroma during the wine making process.Based on the experiment of bacteria growth and enzyme production，the Logistic equation，DoseResp equation and Nelder equation were used to build a dynamic model for the yeast growth and enzyme production，and the relationship of yeast growth rate and enzyme producing rate in this study.After the nonlinear fitting by Origin 8.0 software，and graphed by the method of Lineweaver-Burk，the kinetic parameters of Kmand Vmaxvalues for the enzymes of this yeast producing beta-D-glycosidase were detected.The result shows that the bacteria growth and enzyme production dynamic model of KDLYS9-3 were correlation with the test results，the dynamic model equation can reflect the variety regulation of bacteria growth and enzyme production，the correlation was part of the coupling relationship.The yeast began to produce enzymes after 8 h growth，and it produced a large number of enzymes when entered the logarithm of growth.The growth of yeast was stable after 39 h.With the decline phase of the yeast，the enzyme also stopped producing.The kinetic parameters of the KDLYS 9-3 enzyme were obtained by the double reciprocal method.The Kmwas 8.492 579 mmol/L，and the Vmaxwas 1.030 715 （mu/L）/min.The results helped set the basement of the strains'theoretical foundation and practical application.

Indigenous Saccharomyces cerevisiae KDLYS9-3；beta-D-glucosidase；kinetic equation；Km；Vmax

10.3969/j.issn.1005-6521.2017.23.032

河北省自然科學基金項目（C2011208028）；河北省科技支撐計劃項目（16227106D）

李艷（1958—），女（漢），教授，本科，研究方向：傳統(tǒng)發(fā)酵工程創(chuàng)新技術、葡萄酒和果酒釀造。

2017-08-14