大麥苗總黃酮提取工藝的優化*

張輝，喬勇進，戚文元

1(上海市農業科學院作物育種栽培研究所，上海，201403)

2(上海理工大學醫療器械與食品學院，上海，200093)3(上海束能輻照技術有限公司，上海，201401)

大麥(Hordeum vulgare)，屬大麥屬(Hordeum)，早熟禾科(Poaceae)，其嫩苗富含葉綠素、類黃酮、維生素、抗氧化酶及蛋白質等多種功能營養成分［1－2］。據報道大麥苗黃酮類化合物主要有皂草苷(saponarin)和大麥黃苷(lutonarin)［3－4］，其自由基清除能力和抗氧化性能顯著高于其他植物黃酮和黃酮碳苷［5］，具有抗癌、抗炎及防治類風濕關節炎等生理功能［6－7］。目前，天然植物黃酮類化合物提取原料有銀杏葉［8］、杜仲葉［9］及玉米須［10］等。以大麥苗為原料進行總黃酮定向提取并對其提取工藝參數進行響應面法優化的研究鮮有報道。大麥苗因其原料具有種植周期短、成本低、產量大、易加工及副產物利用價值大等優點，可作為總黃酮提取的原料。已有研究表明超臨界流體萃取、超聲波及微波輔助提取可增加植物黃酮類化合物的提取量［11－13］，而本文以乙醇溶液為提取溶劑，對大麥苗總黃酮進行常溫浸漬提取，采用中心組合試驗設計(central composite design，CCD)，優化其最佳工藝條件。

1 材料與方法

1.1 材料與試劑

六葉期大麥苗(品種:大中88－91)，采收于上海市農業科學院試驗田;無水乙醇(分析純 )、NaNO2(BR)、AlCl3(BR)、NaOH(BR)、蘆丁(純度 ＞95%)，國藥集團化學試劑有限公司。

1.2 儀器與設備

1.3 試驗設計

1.3.1 葉綠素提取工藝流程

大麥苗鮮樣→挑選、清洗、瀝干→真空冷凍干燥(預凍溫度，－45℃;冷阱溫度，－37℃;升華溫度，25℃;真空度，170 Pa)→粉碎(100目)→浸提→過濾→濾液→減壓濃縮成膏狀(水浴溫度，30℃;真空度，110 Pa)

1.3.2 單因素試驗設計

根據預試驗結果，發現乙醇體積分數、液固比和提取時間等3個因素均對大麥苗總黃酮提取量影響較大。故分別取其不同水平進行單因素試驗，每組試驗均進行3次平行試驗，取其平均值進行研究。

1.3.2.1 不同乙醇體積分數對大麥苗總黃酮提取效果的影響

將粉碎干燥后的大麥苗過100目篩，稱取麥苗粉約2.50 g，置于250 mL錐形瓶中，按液固比為30∶1(mL∶g)加入乙醇體積分數分別為55%、65%、75%、85%、95%的乙醇溶液，用封口膜將錐形瓶密封并置于搖床中以150 r/min轉速搖晃浸提90 min，研究不同乙醇體積分數對大麥苗總黃酮提取效果的影響。

在經歷了這樣特殊的課程之后，從新加坡中學走出來的學生具有非常敏銳的“批判性思維”和“辯證思維”，他們的邏輯思維能力快速成長，也同樣對于他們的學術課程大有助益。

1.3.2.2 不同液固比對大麥苗總黃酮提取效果的影響

將粉碎干燥后的大麥苗過100目篩，稱取麥苗粉約2.50 g，置于250 mL錐形瓶中，按液固比分別為10∶1、20∶1、30∶1、40∶1、50∶1(mL∶g)加入乙醇體積分數為95%的乙醇溶液，用封口膜將錐形瓶密封并置于搖床中以150 r/min轉速搖晃浸提90 min，研究不同液固比對大麥苗總黃酮提取效果的影響。

1.3.2.3 不同提取時間對大麥苗總黃酮提取效果的影響

將粉碎干燥后的大麥苗過100目篩，稱取麥苗粉約2.50 g，置于250 mL錐形瓶中，按液固比為30∶1(mL∶g)加入乙醇體積分數為95%的乙醇溶液，用封口膜將錐形瓶密封并置于搖床中以150 r/min轉速搖晃分別浸提 30、60、90、120、150 min，研究不同提取時間對大麥苗總黃酮提取效果的影響。

1.3.3 中心組合試驗設計

采用中心組合試驗設計對大麥苗總黃酮提取工藝參數進行優化，根據單因素試驗結果，確定了CCD的試驗水平(表1)。根據試驗結果，可建立二次非線性回歸模型(式 1)［14］。

式中:Y，總黃酮提取量預測值;Xi，Xj，試驗因素的編碼值;β0，βii，βij，截距及回歸常數(其中，Xi=，式中X為無量綱的編碼值，x為變量i的實ii際值，xio為變量i的零水平實際值。Δxi為變量i的步長實際值，即編碼值變化一個單位所對應的實際變化值)。

因子貢獻率采用式(2)計算［15］:

式中:SSj，因子平方和，SST，模型總平方和。

表1 中心組合試驗設計因子及編碼水平Table 1 Values factors and their coded levels in central composite experimental design

1.3.4 指標測定

浸提完成后將提取液置于10 000 r/min轉速下低溫(10℃)離心15 min，取1 mL上清液，采用硝酸鋁比色法［16］測定其總黃酮含量。總黃酮提取量采用式(3)進行計算。

式中:Y，大麥苗總黃酮提取量，μg/g;m，提取液中總黃酮含量，μg/mL;M，提取樣品質量，g。

1.3.5 數據處理

CCD試驗數據處理采用Design-Expert 8.0.6.1(State-Ease，Inc.，USA)、SAS 9.2(SAS Inc.，USA)，Origin 8.0(Electronic Arts Inc.，USA)繪制單因素圖及效應圖。

2 結果與分析

2.1 單因素試驗結果

圖1表明，大麥苗總黃酮提取量隨著乙醇體積分數的增大而增高。乙醇體積分數為95%時，總黃酮提取量最高。這可能是由于大麥苗中黃酮類化合物為弱極性物質，乙醇體積分數越大，其極性越接近大麥苗中黃酮類化合物的極性，導致總黃酮溶出量升高。

圖1 乙醇體積分數對總黃酮提取量的影響Fig.1 Effect of ethanol volume fraction on the flavonoids extraction yield

液固比對液相主體和固相主體之間的濃度差，即傳質推動力，影響較大。此外，液固比也會影響液固兩相的表觀傳質系數，從而影響傳質速率。圖2表明，隨著液固比的增加，大麥苗總黃酮提取量呈現先快速升高而后增加平緩的趨勢。這是因為增加液固比，傳質動力增加，黃酮就更容易溶出，損失減少;但液固比大于30 mL/g后增加速度趨于平緩。由于液固比較大會造成提取溶劑的浪費及后期濃縮的困難，從總黃酮提取效果和經濟邊際效應兩方面綜合考慮，認為進行大麥苗總黃酮提取時液固比不宜高于40 mL/g。

圖2 液固比對總黃酮提取量的影響Fig.2 Effect of liquid-solid ratio on the flavonoids extraction yield

由圖3可知，總黃酮提取量隨著提取時間的延長呈現先逐漸升高而后降低的趨勢。大麥苗總黃酮提取量在提取時間為120 min時達到最大值，而在提取時間為150 min時略微降低。這可能因為延長提取時間使得更多的黃酮類化合物溶出，但一些黃酮類化合物易于氧化，提取時間過長會導致總黃酮的溶出量小于其降解量，從而使得其提取量略微降低。

圖3 提取時間對總黃酮提取量的影響Fig.3 Effect of extraction time on the flavonoids extraction yield

2.2 CCD試驗結果

2.2.1 回歸模型的建立

試驗共設20次試驗，其中包括14次析因試驗及6次中心點試驗。所有試驗均進行3次平行試驗，取其平均值作為響應值，試驗按隨機順序進行，其結果見表2。

利用Design-Expert 8.0軟件對表2結果進行統計分析，并建立三元二次回歸方程(式3):

表2 CCD設計矩陣及響應值Table 2 Central composite design matrix and responses

2.2.2 回歸模型的統計檢驗

由表3可知，判定系數(R2)為0.994 3，說明模型相關性較好，試驗因子對響應值有較大影響。校正判定系數(adjR2)為0.989 1，表明98.91%的試驗數據的變異性可用此回歸模型來解釋。變異系數(C.V.)為1.66%，表示試驗的可信度和精確度較好。本試驗精密度(adeq precision)為40.536 ＞4.0，表明其為一個適宜的信號。模型 F值為192.61，表明其達到極顯著水平(P＜0.01)，即乙醇體積分數、液固比及提取時間等影響因素對總黃酮提取量的共同影響具有差異極顯著性，且模型一次項、二次項均達到顯著水平。此外，失擬F值為1.54(P=0.323 2＞0.05)，表明失擬值和純誤差沒有顯著性關系，即回歸模型在被研究的整個回歸區域不失擬。綜上可知，該所建立模型可用于指導試驗。此外，表3還表明，各一次項 、二次項和交互項X2X3、X1X3均對大麥苗總黃酮提取量影響極顯著(P＜0.01)，交互項X1X2對總黃酮提取影響不顯著 (P＞0.05)。

表3 回歸模型顯著性檢驗結果Table 3 Results of regression model significance test

圖4 因子貢獻率Fig.4 Percentage contributions of factors

2.2.3 模型效應分析

2.2.3.1 主效應分析

由于模型存在交互效應，故不宜采用回歸系數的絕對值大小來直接比較二次項、交互項作用的大小。因此，本試驗均采用因子貢獻率來衡量各因子對總黃酮提取量的影響大小。圖4－a表明，二次項效應＞一次項效應＞交互項效應。圖4－b表明，模型各項間乙醇體積分數的二次效應影響最顯著，其次是其液固比，而乙醇體積分數及液固比與乙醇體積分數的交互作用對黃酮提取量影響最小。3個試驗因素的效應大小依次為液固比 ＞提取時間 ＞乙醇體積分數。

2.2.3.2 兩因素效應分析

由圖5－a、5－c知，大麥苗總黃酮的提取量隨液固比的增加呈現先快速增長再緩慢增長的趨勢。液固比的增加可增加傳質動力，有助于總黃酮的溶出，但當液固比達到一定程度時，已溶出大部分總黃酮，若再增加液固比，則總黃酮提取量增加緩慢。由5－a、5－e可看出，總黃酮提取量隨著乙醇體積分數的增加呈現先增大而后減小的趨勢。這可能是由于在一定范圍內，乙醇體積分數的增大減小了提取液極性，使其更接近于總黃酮的極性，從而增加了總黃酮的溶出。但過度增加乙醇體積分數會導致提取液極性過小，從而偏離大麥苗總黃酮極性。5－c、5－e均表明，大麥苗總黃酮提取量隨著提取時間的延長曾先前期增加迅速，后期增加緩慢的趨勢。這可能是由于前期已溶出大部分黃酮，后期原料自身黃酮較少，故增加緩慢;也可能由于后期有部分黃酮發生降解，但降解速率仍小于溶出速率，從而呈現增長緩慢的現象。圖5－b為液固比與乙醇體積分數的等高線圖，其接近圓形，表明二者交互作用不顯著，圖5－d、5－e分別為液固比與提取時間的等高線圖及乙醇體積分數與提取時間的等高線圖，其均為橢圓形狀，表明二者交互效應顯著，這均與方差分析結果一致。

圖5 總黃酮提取量響應面圖及等高線圖Fig.5 Response surface plots and contour plots of total flavonoids extraction yield

2.3 反應條件的優化與模型驗證

規范分析［17］(canonical analysis)能求出響應面的拐點并判斷其類型(最大值點、最小值點、馬鞍形點)。對試驗結果進行規范分析(表4)，得到大麥苗總黃酮提取的一個穩定點，該點處理論響應值為950.36。由于變量X2的編碼值的特征值為0.02，其較接近零值，在試驗區域內可能會存在一個穩定嶺系統，故一般需進行嶺脊分析。

表4 規范分析結果Table 4 Results of canonical analysis

嶺脊分析(aidge analysis)是在以原始設計中心點為球心，以r為半徑的超球面與響應面的交點所形成的軌跡范圍內探測最佳響應值［17］。嶺脊分析結果(表 5)表明當超球面半徑(r)為 0.8，X1、X2、X3的特征值分別為 0.728 876、0.021 37、0.329 069 時，響應值達到最大(950.346 3 μg/g)，其與規范分析所得的穩定點較一致，故無穩定嶺系統的存在，穩定點為理論上的最優點。

表5 脊嶺分析結果Table 5 Results of ridge analysis

考慮到試驗操作的便利，將總黃酮提取工藝參數修正為液固比29、乙醇體積分數90% 、提取時間100 min。采用修正后的工藝參數進行3次平行驗證試驗，試驗結果測得大麥苗總黃酮提取量為946.63 μg/g，與模型預測值相差 3.73 μg/g，差異在誤差允許范圍內。

3 結論

本試驗通過CCD試驗設計建立了影響因子與大麥苗總黃酮提取量之間的三元二次回歸模型，確定了大麥苗總黃酮最佳提取工藝參數為液固比29 mL/g、乙醇體積分數90% 、提取時間100 min，在此工藝條件下，大麥苗總黃酮提取量為946.63 μg/g，其與預測值950.36 μg/g較一致。故該模型能反映大麥苗總黃酮提取過程的內部規律，可較好地預測大麥苗總黃酮的提取情況，并能為黃酮類化合物提取原料的來源、進一步提高大麥嫩苗資源利用率及大麥苗黃酮類化合物工業化生產提供理論參考。

［1］ Singh N，Verma P，Pandey B R.Therapeutic potential of organicTriticum aestivumLinn.(wheat grass)in prevention and treatment of chronic diseases:an overview ［J］.International Journal of Pharmaceutical Sciences and Drug Research，2012，4(1):10－14.

［2］ Peeva V N，Tóth S Z，Cornic G，et al.Thermoluminescence and P700 redox kinetics as complementary tools to investigate the cyclic/chlororespiratory electron pathways in stress conditions in barley leaves［J］.Physiologia Plantarum，2012，144(1):83－97.

［3］ Benedet J A，Umeda H，Shibamoto T.Antioxidant activity of flavonoids isolated from young green barley leaves toward biological lipid samples［J］.Journal of Agricultural and Food Chemistry，2007，55(14):5 499－5 504.

［4］ Osawa T，Katsuzaki H，Hagiwara Y，et al.A novel antioxidant isolated from young green barley leaves［J］.Journal of Agricultural and Food Chemistry，1992，40(7):1 135－1 138.

［5］ Kamiyama M，Shibamoto T.Flavonoids with potent antioxidant activity found in young green barley leaves［J］.Journal of Agricultural and Food Chemistry，2012，60(25):6 260－6 267.

［6］ Ben-Arye E，Goldin E，Wengrower D，et al.Wheat grass juice in the treatment of active distal ulcerative colitis:a randomized double-blind placebo-controlled trial［J］.Scandinavian Journal of Gastroenterology，2002，37(4):444－449.

［7］ Dey S，Sarkar R，Ghosh P，et al.Effect of Wheat grass Juice in supportive care of terminally ill cancer patients-a tertiary cancer center experience from India［J］.J Clin Oncol，2006，24:8 634.

［8］ Chen S，Xing X，Huang J，et al.Enzyme-assisted extraction of flavonoids from Ginkgo biloba leaves:Improvement effect of flavonol transglycosylation catalyzed byPenicilliumdecumbens cellulose［J］.Enzyme and Microbial Technology，2011，48(1):100－105.

［9］ Huang W，Xue A，Niu H，et al.Optimised ultrasonic-assisted extraction of flavonoids from Folium eucommiae and evaluation of antioxidant activity in multi-test systems in vitro［J］.Food Chemistry，2009，114(3):1 147 － 1 154.

［10］ Liu J，Lin S，Wang Z，et al.Supercritical fluid extraction of flavonoids from Maydis stigma and its nitrite-scavenging ability［J］.Food and Bioproducts Processing，2011，89(4):333－339.

［11］ Qu Y，Tong G，Wang B，et al.Optimized microwave-assisted extraction of flavonoids from Tamarix chinensis and evaluation of antioxidant activityin vitro［C］ //Remote Sensing，Environment and Transportation Engineering(RSETE)，2011 International Conference on.IEEE，2011:7 697－7 700.

［12］ Zhang Q A，Fan X H，Li T，et al.Optimisation of ultrasound extraction for flavonoids from semen astragali complanati and its identification by HPLC‐DAD‐MS/MS［J］.International Journal of Food Science＆ Technology，2013.doi:10.1111/ijfs.12178.

［13］ Shan B，Xie J，Zhu J，et al.Ethanol modified supercritical carbon dioxide extraction of flavonoids fromMomordicacharantiaL.and its antioxidant activity［J］.Food and Bioproducts Processing，2012，90(3):579－587.

［14］ Maran J P，Manikandan S.Response surface modeling and optimization of process parameters for aqueous extraction of pigments from prickly pear(Opuntia ficus-indica)fruit［J］.Dyes and Pigments，2012，95(3):465 －472.

［15］ Maran J P，Manikandan S，Thirugnanasambandham K，et al.Box-behnken design based statistical modeling for ultrasound-assisted extraction of corn silk polysaccharide［J］.Carbohydrate Polymers，2013，92(1):604－611.

［16］ Jia Z S，Tang M C，Wu J M.The determination of flavonoid contents in mulberry and their scavenging effects on superoxide radicals［J］.Food Chemistry，1999，64(4):555－559.

［17］ Draper N R.Ridge analysis of response surfaces ［J］.Technometrics，1963，5(4):469－479.