響應曲面法優化蟲草素和蟲草多糖的綜合提取工藝

鄧 黎，王曉虹，韓 濤，皮 立，韓 發

中國科學院西北高原生物研究所，西寧 810008

現代科學論證冬蟲夏草和蛹蟲草不僅具有特殊的營養價值，而且有明顯的藥用價值。其中尤以蟲草素、蟲草多糖、蟲草酸等多種生物活性物質的藥用價值最為顯著［1-4］。因此，冬蟲夏草和蛹蟲草受到了廣泛關注，已有大量藥品、保健食品上市，但這些產品多數停留在粗加工水平，缺乏市場競爭力，因此，要保證蟲草產業能夠持續穩定的發展，就必須要提高其加工產品的科技含量。因此，研究冬蟲夏草和蛹蟲草中主要生物活性物質的綜合提取工藝和提高其提取率具有重要的經濟價值和廣闊的市場應用潛力［5，6］。

近年來，已有研究人員對蟲草素和蟲草多糖的綜合提取進行了工藝改進和優化，但是所得到的優化工藝均是在單因素或正交試驗條件下獲得的［7，8］。響應面分析法(response surface methodology，RSM)是一種優化工藝條件的有效方法，可檢查一個或多個響應變量與一系列試驗變量之間的關系，確定試驗因素及其交互作用在工藝過程中對指標響應值的影響，并精確地表述因素和響應值之間的關系［9］。因此，本研究中將響應面分析的方法應用于蟲草中蟲草素和蟲草多糖的綜合提取工藝研究中，建立得率與考察因素之間關系的數學模型，為工業化生產提供技術參數和理論依據。

在比較水熱回流法和超聲波水提法對綜合提取蟲草素和蟲草多糖的影響基礎上，采用超聲波提取法進行初提，經醇沉、Sevage法除蛋白、乙醇分級沉淀后獲得蟲草多糖［3］，醇沉后獲得的上清液用大孔樹脂和反相硅膠分離純化，制備型液相色譜分離制備獲得蟲草素［10］。最終利用響應曲面法考察超聲提取時間、溫度和液固比對蟲草素和蟲草多糖綜合提取率的影響，得到最優提取工藝并予以驗證。

1 材料與方法

1.1 儀器、材料與試劑

Agilent 1200高效液相色譜儀，色譜柱Agilent Eslipse-XDB C18(5 μm，250 mm × 4.6 mm)，美國安捷倫公司;N-1001型旋轉蒸發儀(上海愛朗儀器有限公司);KQ 5200 DE型數控超聲波清洗器(昆山市超聲儀器有限公司);METTLER TOLEDO PL 203和 XS 204電子天平(瑞士梅特勒-托利多公司);P70D17TL-D 5微波爐(格蘭仕微波爐電器有限公司);制備型高效液相色譜儀(LC-8A，日本);DHG-9070 A型電熱恒溫故風干燥箱(上海精宏實驗設備有限公司)。

野生冬蟲夏草購于青海玉樹，人工培養蛹蟲草CSM-8菌株由西寧世峰生物科技有限公司提供，樣品經微波干燥后粉碎過60目篩備用;蟲草素標準品、D-葡萄糖標準品由國家標準物質研究中心提供;甲醇為色譜純，其余化學試劑均為分析純;實驗用水均為蒸餾水。

1.2 實驗方法

1.2.1 蟲草素、蟲草多糖提取率計算公式及含量測定方法

蟲草素提取率計算公式:

式中C1為依據標準曲線計算出的被測液中蟲草素濃度(mg/mL);C2為樣品溶液中蟲草素濃度(mg/mL)。

蟲草多糖提取率計算公式:

式中C1為依據標準曲線計算出的被測液中蟲草多糖濃度(mg/mL);C2為樣品溶液中蟲草多糖濃度(mg/mL)。

采用高效液相色譜法(HPLC)測定蟲草素含量，色譜柱為Agilent Eslipse-XDB C18柱;流動相為甲醇-水(12∶88);流速0.8 mL/min;紫外檢測波長260 nm;進樣量 10 μL;室溫。

苯酚-硫酸法測定蟲草多糖含量，樣品溶液多糖含量的測定在標準曲線項下操作，測定OD值并得出樣品中蟲草多糖含量。

1.2.2 蟲草素和蟲草多糖綜合提取分離工藝研究

1.2.2.1 提取方法比較實驗

比較水熱回流法和超聲波水提法對綜合提取蟲草素和蟲草多糖的影響。

水熱回流提取法:精確稱取0.10 g冬蟲夏草和1.00 g預處理后的CSM-8子實體粉末，按液固比30 mL/g加入蒸餾水熱回流提取6h，HPLC測定蟲草素含量，苯酚-硫酸法測定蟲草多糖含量并計算提取率，實驗重復3次。

超聲波水提取法:精確稱取1.00 g預處理后的CSM-8和0.10 g冬蟲夏草子實體粉末，按液固比30 mL/g加入蒸餾水超聲波提取30 min，溫度55℃，HPLC測定蟲草素含量，苯酚-硫酸法測定蟲草多糖含量并計算提取率，實驗重復3次。

1.2.2.2 超聲波提取單因素實驗

分別考察超聲溫度(5、30、55、80、90 ℃)、超聲提取時間(20、40、60、80、100 min)及液固比(10、20、30、40、50 mL/g)、超聲提取次數(1 次、2 次、3 次、4次、5次)對人工蛹蟲草CSM-8菌株子實體干粉中蟲草素和蟲草多糖含量的影響。

1.2.2.3 響應曲面法實驗設計

在單因素實驗的基礎上選擇對響應值(蟲草素和蟲草多糖提取率)有顯著影響的因素，利用Design Expert軟件對實驗數據進行回歸分析。每一變量的低、中、高水平分別以-1、0、1進行編碼，模型選用二次方程:

其中Y為響應值(蟲草素和蟲草多糖提取率)，β0、β1、β2、β3、β12、β13、β23、β11、β22、β33為方程系數，x1、x2、x3為自變量編碼值。

通過Box-Behnken設計和響應面分析來進一步考察超聲提取時間、提取溫度和液固比之間的相互影響，并確定超聲波法綜合提取子實體中蟲草素和蟲草多糖的最佳工藝參數。

稱取1.00 g CSM-8菌株子實體粉末和0.10 g冬蟲夏草子實體粉末，分別按最佳提取工藝參數進行驗證實驗，HPLC測定蟲草素含量，苯酚-硫酸法測定蟲草多糖含量，并計算提取率，實驗重復3次。

2 結果與分析

2.1 蟲草素和蟲草多糖綜合提取分離工藝研究

2.1.1 提取方法比較實驗結果與分析

比較水熱回流法和超聲波水提法綜合提取蟲草素和蟲草多糖的提取效果，結果表明(表1):超聲波水提法優于水熱回流法，這是因為超聲波使菌絲細胞破碎，有利于細胞內物質的溶出，使提取效果更好。故選擇超聲波水提法為綜合提取蟲草素和蟲草多糖的方法。

表1 提取方法與蟲草素、多糖提取率的關系(n=3)Table 1 Comparison of extraction results using water reflux method and Ultrasonication method(n=3)

2.1.2 超聲波提取單因素實驗結果與分析

2.1.2.1 超聲溫度對蟲草多糖和蟲草素含量的影響

在超聲時間60 min、液固比7.5 mL/g、超聲提取次數2次的條件下考察微波功率對人工蛹蟲草子實體中蟲草多糖、蟲草素含量的影響，結果見表2。確定較優草上提取溫度范圍為30～80℃。

表2 超聲溫度對人工蛹蟲草CSM-8菌株子實體干粉中蟲草多糖、蟲草素含量的影響Table 2 Effect of extraction temperature on the content of polysaccharides，cordycepin in powdered fruit body of CSM-8

2.1.2.2 超聲時間對蟲草素和蟲草酸含量的影響

在超聲溫度為55℃、液固比7.5 mL/g、提取次數2次條件下考察超聲時間對人工蛹蟲草子實體中蟲草素、蟲草多糖含量的影響，結果見表3。確定較佳超聲時間范圍為40～80 min。

表3 超聲時間對人工蛹蟲草CSM-8菌株子實體干粉中蟲草多糖、蟲草素含量的影響Table 3 Effect of ultrasonic time on the content of polysaccharides，cordycepin in powdered fruit body of CSM-8

2.1.2.3 液固比對蟲草素和蟲草多糖含量的影響在超聲溫度55℃、超聲時間60 min、超聲提取次數2次條件下考察液固比對人工蛹蟲草子實體中蟲草多糖、蟲草素含量的影響，結果見表4。確定較優液固比范圍為20～40 mL/g。

表4 液固比對人工蛹蟲草CSM-8菌株子實體干粉中蟲草多糖、蟲草素含量的影響Table 4 Effect of Water/material ratio on the content of polysaccharides，cordycepin in powdered fruit body of CSM-8

表5 提取次數對人工蛹蟲草CSM-8菌株子實體干粉中蟲草多糖、蟲草素含量的影響Table 5 Effect of times of extraction on the content of polysaccharides，cordycepin in powdered fruit body of CSM-8

2.1.2.4 超聲提取次數對蟲草素和蟲草酸含量的影響

在超聲溫度55℃、超聲時間60 min、液固比30 mL/g條件下考察超聲次數對人工蛹蟲草子實體中蟲草多糖、蟲草素含量的影響，結果見表5。確定較優提取次數為2次。

2.1.3 響應曲面實驗結果與分析

2.1.3.1 響應曲面法分析因素的選取及分析方案

在單因素實驗的基礎上固定提取次數為2次，采用Box-Behnken設計方案，分別以提取溫度、提取時間和液固比對應3個獨立變量X1、X2、X3，并以-1、0、1分別代表變量的水平。

對自變量進行編碼(表6)。式中xi為變量的編碼值，Xi為變量的真實值，X0為實驗中心點變量的真實值，ΔX為單變量增量，蟲草素和蟲草多糖提取率Y為響應值，實驗方案及結果見表6和7。

表6 實驗因素水平及編碼Table 6 Levels and codes of variable chose for Box-Behnken design

2.1.3.2 結果分析

Box-Behnken設計方案及蟲草素和蟲草多糖提取率的測定結果見表7。利用Design Expert軟件對表7數據進行二次多元回歸擬合，得到回歸方程為:

表7 Box-Behnken設計方案及蟲草素和蟲草多糖提取率的測定值Table 7 Box-Behnken design and the responses of the total polysaccharide and cordycepin yield

對回歸方程進行顯著性檢驗及方差分析，結果見表8。由表8得出:兩個模型的P均＜0.0001，表明兩個模型都極顯著且失擬項不顯著(模型1，P=0.9388＞0.05;模型2，P=0.3796＞0.05)，相關系數R2=0.9955、R2=0.9888，說明兩個模型擬合程度都很好。模型的調整確定系數 AdjR2=0.9897、AdjR2= 0.9745，說明模型 1能解釋98.97%的響應值的變化，模型2能解釋97.45%的響應值的變化，可以用模型1、2對超聲波水提法綜合提取蟲草素和蟲草多糖提取率進行分析及預測。

表8 回歸模型的方差分析Table 8 ANOVA analysis for response surface quadratic model

模型1回歸方程的系數顯著性檢驗如下:一次項A(P＜0.0001)、C(P＜0.0001)，二次項 A2(P＜0.0001)、B2(P＜0.0001)均達到極顯著水平，一次項B(P=0.0036＜0.01)達到顯著水平，表明提取溫度、超聲提取時間、液固比對蟲草素提取率有顯著影響;二次項AB(P=0.0014＜0.0001)、BC(P=0.0050＜0.001)，表明提取溫度與超聲提取時間、提取時間與液固比間的交互作用對蟲草素提取率的影響極顯著。

模型2回歸方程的系數顯著性檢驗如下:一次項A(P＜0.0001)、B(P=0.0003＜0.001)，二次項 A2(P ＜0.0001)、B2(P ＜0.0001)、C2(P ＜0.0001)均達到極顯著水平;一次項C(P=0.0020＜0.01)達到顯著水平，表明提取溫度、超聲提取時間、液固比對蟲草多糖提取率有顯著影響;二次項AB(P=0.0003＜0.001)、AC(P=0.0383＜0.05)、BC(P=0.0348 ＜0.05)，表明提取溫度與超聲提取時間、提取溫度與液固比以及提取時間與液固比的交互作用對蟲草多糖提取率的影響顯著。

利用Design Expert 7.0.0軟件對表4數據進行二次多元回歸擬合，所得到的二次回歸方程Y1的響應面及其等高線見圖1、2和3。

從圖1可知，蟲草素提取率隨提取時間的延長先升高后下降，呈拋物線趨勢，當提取時間達到60 min左右時，提取率達到最大值，其原因可能是隨著時間的增大，可溶性成分溶出的越多，使提取液粘度增大，在后期操作中有所損耗;提取時間和提取溫度的交互作用對蟲草酸素提取的影響顯著，表現為曲面較陡，響應值隨提取溫度的變化率大于隨提取時間的變化率，說明提取溫度對提取率的影響大于超聲提取時間;蟲草素提取率隨提取溫度的升高而升高，當溫度達到67.5℃后，提取率隨溫度的升高變化不大。圖2可知，提取溫度和液固比的交互作用對蟲草素提取率的影響不顯著，表現為曲面平滑。從圖3可知，提取時間和液固比對提取率的影響均成拋物線形，即隨提取時間和液固比的增大，提取率呈先升高后降低的趨勢;提取時間和液固比的交互作用對蟲草素提取的影響顯著，表現為曲面較陡平滑;響應值隨液固比的變化率大于提取時間，說明液固比對提取率的影響大于提取時間。

圖1 提取溫度和提取時間影響蟲草素提取率的等值線與響應面Fig.1 Response surface plots(3-D)and contour plots(2-D)showing the effects of variables(X1:extraction temperature;X2:extraction time)on the response Y1

圖2 提取溫度和液固比影響蟲草素提取率的等值線與響應面Fig.2 Response surface plots(3-D)and contour plots(2-D)showing the effects of variables(X1:extraction temperature;X3:liquid-solid ratio)on the response Y1

圖3 提取時間和液固比影響蟲草素提取率的等值線與響應面Fig.3 Response surface plots(3-D)and contour plots(2-D)showing the effects of variables(X2:extraction time;X3:liquid-solid ratio)on the response Y1

圖4 超聲提取時間和提取溫度影響蟲草多糖提取率的等值線與響應面Fig.4 Response surface plots(3-D)and contour plots(2-D)showing the effects of variables(X1:extraction temperature;X2:extraction time)on the response Y2

圖5 提取溫度和液固比影響蟲草多糖提取率的等值線與響應面Fig.5 Response surface plots(3-D)and contour plots(2-D)showing the effects of variables(X1:extraction temperature;X3:liquid-solid ratio)on the response Y2

圖6 提取時間和液固比影響蟲草多糖提取率的等值線與響應面Fig.6 Response surface plots(3-D)and contour plots(2-D)showing the effects of variables(X2:extraction time;X3:liquid-solid ratio)on the response Y2

利用Design Expert 7.0.0軟件對表4數據進行二次多元回歸擬合，所得到的二次回歸方程Y2的響應面及其等高線見圖4、5和6。從圖4可知，提取時間和提取溫度的交互作用對蟲草多糖提取的影響較顯著，表現為曲面較陡，響應值隨提取溫度的變化率大于隨提取時間的變化率，說明提取溫度對提取率的影響大于超聲提取時間;蟲草多糖提取率隨提取時間的延長先升高后降低，成拋物線趨勢;蟲草多糖提取率隨提取溫度的升高而升高，隨后趨于平緩。從圖5可知，提取溫度和液固比的交互作用對蟲草素提取的影響顯著，響應值隨提取溫度的變化率大于隨液固比的變化率，說明提取溫度對提取率的影響大于液固比。從圖6可知，提取時間和液固比對蟲草多糖提取率的影響均成拋物線形，即隨提取時間和液固比的同時增大，提取率呈先升高后降低的趨勢;提取時間和液固比的交互作用對蟲草多糖提取的影響顯著，表現為曲面較陡;響應值隨液固比的變化率大于隨提取時間的變化率，說明液固比對蟲草多糖提取率的影響大于超聲提取時間。

利用Design Expert 7.0.0軟件對回歸模型進行進一步的典型性分析，得到最優提取工藝條件為:提取溫度72℃，液固比30 mL/g，提取時間77 min，提取2次，蟲草素提取率的理論值可達53.27%，蟲草多糖提取率的理論值可達41.92%。

2.2 最佳工藝驗證實驗

最優工藝條件的驗證實驗結果如表9所示，3次重復實驗的蟲草素提取率的平均值為50.90%，RSD=1.09%，蟲草多糖提取率的平均值為37.28%，RSD=1.42%，表明二次多項式建立的數學模型具有良好的預測性，所選工藝條件重現性好。

表9 驗證實驗結果(n=3)Table 9 Verification of optimal extraction conditions for the total polysaccharide yield(n=3)

2.3 超聲綜合提取冬蟲夏草子實體中蟲草素和蟲草多糖

表10 超聲提取冬蟲夏草子實體中蟲草素和蟲草多糖實驗結果(n=3)Table 10 Extraction yields of cordycepin and polysaccharide of C.sinensis using the optimized ultrasonic extraction conditions(n=3)

按照響應曲面法確定的最佳提取工藝參數，稱取0.10 g冬蟲夏草子實體粉末進行超聲提取，實驗結果見表10。設定提取溫度73℃，液固比35 mL/g，提取時間65 min，提取2次，超聲綜合提取冬蟲夏草子實體中蟲草素的提取率平均值為50.20%，RSD=0.76%，蟲草多糖的提取率平均值為36.92%，RSD=0.64%，表明所選工藝條件重現性好。

3 結論

綜合提取蟲草素和蟲草多糖工藝條件和路線:采用超聲波提取法進行初提，提取液濃縮、醇沉、離心，沉淀用Savage法除蛋白，乙醇分級沉淀獲得蟲草多糖，上清液用大孔樹脂和反相硅膠分離純化，制備型液相色譜分離制備獲得蟲草素。經響應曲面法確定的最優提取工藝條件為:提取溫度72℃，液固比30 mL/g，提取時間77 min提取2次。以冬蟲夏草和CSM-8菌株為原料，按照最佳提取工藝條件進行超聲提取，測得蟲草素提取率平均值分別為50.90%和50.20%，蟲草多糖提取率平均值分別為37.28%和36.92%。

1 Liu DZ(劉東澤)，Chen W(陳偉)，Gao XH(高新華)，et al.Progressofresearch on cordycepin (3.-Deoxyadenosine).Acta Agric Shanghai(上海農業學報)，2004，20(2):89-93.

2 Kim JR，Yeon SH，Kim HS，et al.Larvicidal activity against Plutellaxylostella of cordycepin from the fruiting body of Cordyceps militaris.Pest Manage Sci，2002，58:713-717.

3 Yu RM，Wang L，Zhang H.Isolation，purification and identification of polysaccharides from cultured Cordyceps militaris.Fitoterapia，2004，75:662-666.

4 Cunningham KG，Manson W，Sping FS，et al．Cordycepin，a metabolic product isolated from cultures of Cordyceps militaris.Nature，1950，66:949.

5 Zou XY(鄒向英).Comparative analysis of effective compositions in three Cordyceps.J Northeast Agric Univ(東北農業大學學報)，2010，41:139-142.

6 Chen W(陳偉)，Wu W(吳畏)，Gao XH(高新華)，et al.Study on continuous separation technology of active components from Cordyceps militaris.Acta Agric Shanghai(上海農業學報)，2007，23(3):62-65.

7 Wang YJ(王英娟)，Li DW(李多偉)，Wang YC(王義潮)，et al.Integrated extracting technology of cordycepin and polysaccharides in Cordyceps militaris.Acta Bot Bor-occid Sin(西北植物學報)，2005，25:1863-1867.

8 Qi B(祁波).Integrated extracting technology of cordycepin and polysaccharides in Cordyceps militaris.Sci Techol Inf(科技信息)，2009:341-343.

9 Wang YF(王永菲)，Wang CG(王成國).The application of response surface methodology.J Central Univ Nation，Nat Sci(中央民族大學學報，自科版)，2005，14:236-240.

10 Chen W(陳偉)，Wu W(吳畏)，Gao XH(高新華)，et al.Research on separation and enrichment of cordycepin by column chromatography.Acta Agric Shanghai(上海農業學報)，2006，22(4):69-71.