響應面法優(yōu)化離子液體雙水相萃取木瓜蛋白酶

曾穎，朱新儒，余壘，李雄，張海德

（海南大學食品學院，海南?？?70228）

番木瓜（Carica papaya L.）原產(chǎn)于墨西哥南部以及鄰近的美洲中部地區(qū)。在我國主要分布在廣東、海南、廣西、云南、福建、臺灣等省（區(qū)）。木瓜蛋白酶（papain）是存在于番木瓜的根、莖、葉和果實內(nèi)的一種低特異性蛋白水解酶，其中在未成熟的乳汁中含量最豐富。木瓜蛋白酶的活性中心含半胱氨酸，屬于含巰基（-SH）肽鏈內(nèi)切酶，它具有酶活高、熱穩(wěn)定性好、天然衛(wèi)生安全等特點，因此在食品[1]、醫(yī)藥[2]、飼料[3]、日化[4]、皮革[5]及紡織[6]等行業(yè)得到廣泛應用。木瓜蛋白酶是100%的天然物質(zhì)，因此在食品及飲料工業(yè)上有獨特的廣闊發(fā)展前途[7]。

離子液體，又稱低溫熔融鹽，是指在室溫或接近室溫下呈現(xiàn)液態(tài)的、完全由陰陽離子所組成的鹽。由于離子液體具有非常獨特的理化性能，非常適合于用作分離提純的溶劑。2003年，Gutowski等采用親水性離子液體1-丁基-3-甲基咪唑鹽酸鹽（[Bmim]Cl）和磷酸鉀（K3PO4）首次提出了離子液體雙水相體系這一概念。離子液體雙水相相較于傳統(tǒng)雙水相既為蛋白質(zhì)的萃取營造了溫和的環(huán)境，又避免了萃取過程中有機溶劑的使用，具備萃取率高、兩相界面清晰、黏度低、綠色環(huán)保、操作簡便等優(yōu)點。因此近年來在分離生物活性物質(zhì)方面迅速發(fā)展并在分離生物活性物質(zhì)[8-9]方面被廣泛研究。

本試驗利用離子液體雙水相對木瓜中的木瓜蛋白酶進行分離純化，旨在從方法上尋求新的突破，為木瓜蛋白酶的純化提供更多的理論依據(jù)和技術(shù)基礎。

1 材料與方法

1.1 材料與試劑

木瓜蛋白酶（BR，>2 000 U/mg）：生工生物（上海）股份有限公司；[C2mim]Ac、[C4mim]Ac、[C6mim]Ac、[C8mim]Ac：蘭州雨陸精細化工有限公司；K2HPO4（分析純）：廣州化學試劑廠。

1.2 儀器與設備

TU1901紫外可見分光光度計：北京普析通用有限責任公司；DK-98-1型恒溫水浴鍋 天津泰斯特儀器有限公司；HJ-5型多功能恒溫攪拌器：常州華奧儀器制造有限公司。

1.3 方法

1.3.1 蛋白質(zhì)測定方法

用Bradford法[10]來測定蛋白質(zhì)的含量。以牛血清白蛋白為標準品配制系列濃度的標準溶液，1.0 mL樣液加入5 mL考馬斯亮藍G-250染色液5.0 mL，反應5 min左右在595 nm波長下測定吸光度。以蛋白質(zhì)濃度為橫坐標，吸光值為縱坐標繪制標準曲線。標準曲線的回歸方程為：A595=0.007x+0.000 5，R2=0.999 2。

1.3.2 雙節(jié)線的測定

在295.15 K的條件下采用濁點法制作雙節(jié)線[11]。分別配制一定濃度（質(zhì)量分數(shù)為60%左右）的離子液體溶液于大試管中，逐滴加入一定濃度的鹽溶液，直到溶液剛好變渾濁，記下消耗的鹽溶液體積。再加入重蒸水使溶液澄清，重復上述操作。計算出在溶液渾濁時離子液體和鹽的質(zhì)量分數(shù)，以鹽的質(zhì)量分數(shù)為橫坐標，離子液體質(zhì)量分數(shù)為縱坐標繪制相圖。

1.3.3 分配系數(shù)測定

配制一定量的的離子液體、無機鹽和木瓜蛋白酶溶液，震蕩混勻后以5 000 r/min離心10 min，讀取上下

式中：cT和cB分別為雙水相體系中上下相的蛋白酶的濃度，mg/mL。

1.3.4 單因素試驗

離子液體側(cè)鏈烷基對分配系數(shù)的影響：以相圖為指導，用不同的離子液體配制相同質(zhì)量組成的雙水相，固定離子液體和K2HPO4的質(zhì)量分數(shù)分別為25%、18%，加入等量適量的酶液，測定木瓜蛋白酶的分配系數(shù)。

離子液體濃度對分配系數(shù)的影響：配制離子液體濃度為20%、24%、28%、32%、36%，K2HPO4濃度固定為18%的雙水相，加入等量適量的酶液，測定木瓜蛋白酶的分配系數(shù)。

K2HPO4濃度對分配系數(shù)的影響：保持[C8mim]Ac濃度為24%不變，根據(jù)相圖，選取K2HPO4濃度為16%、18%、20%、22%、24%的雙水相體系。加入等量適量的酶液，測定木瓜蛋白酶的分配系數(shù)。

pH值對分配系數(shù)的影響：雙水相pH值分別調(diào)節(jié)為 5.0、6.0、7.0、8.0、9.0、10.0、11.0，K2HPO4濃度為 18%，[C8mim]Ac濃度為24%，加入等量適量的酶液，測定木瓜蛋白酶的分配系數(shù)。

1.3.5 響應面試驗

根據(jù)Box-Benhnken的中心組合試驗設計原理，根據(jù)單因素試驗結(jié)果，對離子液體濃度（A）、鹽濃度（B）、pH值（C）這3個因素進行優(yōu)化，采用三因素三水平的響應面分析方法進行試驗設計，以分配系數(shù)（K）作為考察指標，獲得木瓜蛋白酶的最佳條件萃取條件。試驗因素水平見表1。相的體積，測定上下相的蛋白濃度，分配系數(shù)（K）的計算公式如下：

表1 響應面試驗因素水平表Table 1 Factors and levels of RSM（response surface methodology）

2 結(jié)果與分析

2.1 雙節(jié)線的測定及關聯(lián)

本試驗研究了[Cnmim]Ac（n=2，4，6，8）和 K2HPO4的相形成能力，[Cnmim]Ac（n=4，6，8）與 K2HPO4的相圖見圖1。

圖1 [Cnmim]Ac（n=4，6，8）與 K2HPO4的相圖Fig.1 Phase diagram of the[Cnmim]Ac（n=4，6，8）with K2HPO4

雙節(jié)線越靠近原點，則成相能力越好，從圖1可以看出咪唑類醋酸鹽的成相能力為：[C8mim]Ac>[C4mim]Ac>[C6mim]Ac。一般情況下離子液體陽離子側(cè)鏈烷基越長，則離子液體疏水性越強，更容易形成雙水相。而[C2mim]Ac難以與K2HPO4形成雙水相。

用Merchuk[12]方程對雙節(jié)線數(shù)據(jù)進行關聯(lián)。

式中：WIL、WS分別表示離子液體和鹽的質(zhì)量分數(shù)，%；a、b、c為方程參數(shù)。雙節(jié)線數(shù)據(jù)關聯(lián)結(jié)果見表2。

表2 雙節(jié)線數(shù)據(jù)關聯(lián)結(jié)果Table 2 Binodal correlation results

由表2，可以看出對雙節(jié)線數(shù)據(jù)的關聯(lián)效果非常好。

2.2 單因素試驗

2.2.1 離子液體側(cè)鏈烷基對分配系數(shù)的影響

以相圖為指導，用3種離子液體配制相同質(zhì)量組成的雙水相，固定離子液體和K2HPO4的質(zhì)量分數(shù)分別為25%、18%，加入等量適量的酶液，結(jié)果如圖2所示。

木瓜蛋白酶的分配系數(shù)隨著離子液體陽離子側(cè)鏈烷基的增加而增大。咪唑類醋酸鹽離子液體的疏水性隨側(cè)烷基鏈增長緩慢加強，因此長鏈烷基可能更有利于木瓜蛋白酶被萃取到上相中。因此，選取萃取效果最好的[C8mim]Ac來進行下一步試驗。

2.2.2 離子液體濃度對分配系數(shù)的影響

離子液體質(zhì)量分數(shù)變化對木瓜蛋白酶分配行為的影響見圖3。

圖2 離子液體烷基側(cè)鏈變化對木瓜蛋白酶分配行為的影響Fig.2 Effects of the alkyl chain length of ionic liquids on the partition behavior

圖3 離子液體質(zhì)量分數(shù)變化對木瓜蛋白酶分配行為的影響Fig.3 Effects of the mass of ionic liquids on the partition behavior

固定K2HPO4濃度為18%，配制[C8mim]Ac濃度為20%、24%、28%、32%、36%的雙水相，由圖3可知K值先上升后下降，在[C8mim]Ac濃度為28%時取得最大值11.44。隨著[C8mim]Ac濃度增加（>28%），富含離子液體的上相體積也會增加，雖然可能對酶的萃取量變大了，但是上相酶濃度也會隨上相體積的變大而降低，所以總體的K值開始下降。

2.2.3 K2HPO4濃度對分配系數(shù)的影響

K2HPO4濃度變化對木瓜蛋白酶分配行為的影響見圖4。

圖4 K2HPO4濃度變化對木瓜蛋白酶分配行為的影響Fig.4 Effects of the mass of K2HPO4on the partition behavior

于是保持[C8mim]Ac濃度為28%不變，選取K2HPO4濃度為16%、18%、20%、22%、24%。圖4得出，K值先增大后減小，這可能是因為K2HPO4質(zhì)量分數(shù)的增加導致了鹽析作用增強，木瓜蛋白酶轉(zhuǎn)移到上相。但是同時上相的水分子慢慢減少，與蛋白質(zhì)表面的氨基酸殘基形成氫鍵能力減弱，當[C8mim]Ac濃度過大時上相的黏度增加，也會阻礙木瓜蛋白酶的轉(zhuǎn)移能力。

2.2.4 pH值對分配系數(shù)的影響

pH值變化對木瓜蛋白酶分配行為的影響見圖5。

圖5 pH值變化對木瓜蛋白酶分配行為的影響Fig.5 Effects of pH on the partition behavior

根據(jù)圖5得知：隨pH值的增大，K呈現(xiàn)先增大后下降的趨勢，并在pH=9時取得最大值。pH值對離子液體雙水相中酶分配系數(shù)的影響是通過酶與成相物質(zhì)之間的氫鍵效應和電荷效應來實現(xiàn)的。蛋白質(zhì)作為兩性分子，其帶電性質(zhì)受介質(zhì)值的影響[13]。當pH值大于等電點，蛋白質(zhì)帶負電荷，與的離子液體的咪唑陽離子之間發(fā)生靜電作用。咪唑類離子液雙水相體系中，蛋白質(zhì)表面的氨基酸殘基與離子液體陽離子之間的疏水相互作用是萃取分離蛋白質(zhì)的主要動力[14]。

2.3 響應面試驗設計及方差分析

2.3.1 響應面試驗結(jié)果

根據(jù)Box-Benhnken的中心組合試驗設計原理，根據(jù)單因素試驗結(jié)果，對離子液體濃度（A）、鹽濃度（B）、pH 值（C）這 3個因素進行優(yōu)化，以分配系數(shù)（K）作為考察指標。根據(jù)響應而分析試驗設計原理進行響應面試驗，響應面設計及結(jié)果如表3所示。

表3 響應面試驗設計方案及結(jié)果Table 3 Experimental design and results of RSM

續(xù)表3 響應面試驗設計方案及結(jié)果Continue table 3 Experimental design and results of RSM

2.3.2 回歸模型建立及顯著性分析

利用軟件進行多元回歸擬合，得到分配系數(shù)K的二次多項式回歸模型如下：

對回歸模型進行方差分析和回歸系數(shù)的顯著性檢驗，結(jié)果如表4和表5所示。

表4 回歸模型的方差分析Table 4 Variance analysis for the regression model

表5 回歸模型的顯著性檢驗Table 5 Significance test for the regression model

模型的失擬項均不顯著（P=0.479 9>0.05），說明方程對試驗擬合度很好。回歸方程與實測值之間的擬合度可以通過相關系數(shù)（R2）和校正決定系數(shù)（Adj.R2）來解釋。方程3的相關系數(shù)為99.7%，說明模型實測值和預測值之的擬合度都非常好。校正決定系數(shù)為99.31%，說明該回歸模型可以解釋99.31%的數(shù)據(jù)的變異性，因此可利用該模型對分配系數(shù)進行分析和預測。

從回歸方程一次項系數(shù)的絕對值大小可評價該因素對試驗指標影響程度[15]。方程2顯示，各因素對分配系數(shù)K的影響順序為pH值>[C8mim]Ac的質(zhì)量分數(shù)>K2HPO4的質(zhì)量分數(shù)在該水平范圍內(nèi)，AB、BC、AC、A2、B2、C2對K的影響極顯著，A、B對K的影響顯著。

2.3.3 響應曲面分析

在回歸分析的基礎上，根據(jù)回歸方程擬合繪制響應面圖，如圖6所示。

圖6 分配系數(shù)K的響應面和等高線圖Fig.6 Responsive surface and contour plot of partition coefficient K

兩因素間交互作用的強弱可以通過等高線圖來判斷，若等高線為圓形，則表示兩者交互作用不顯著，等高線圖呈橢圓形，說明交互作用顯著[16-17]。AC的交互作用較AB、BC而言對K的影響最為顯著。

2.3.4 最優(yōu)萃取條件的確定及模型驗證

通過Design-Expert 8.0.5軟件進行系統(tǒng)分析，得出 [C8mim]Ac-K2HPO4系統(tǒng)萃取木瓜蛋白酶的最佳萃取條件為：離子液體[C8mim]Ac濃度為28.03%，K2HPO4濃度為20.08%，pH值為8.55，此時分配系數(shù)K的預測值為11.65。按照優(yōu)化條件進行3次平行驗證試驗以考察模型的可靠性，得分配系數(shù)的平均值為11.85，與預測值11.65非常接近，相差僅為1.69%。

3 結(jié)論

在同樣的條件下，分別測定 [Cnmim]Ac（n=2、4、6、8）離子液體和K2HPO4的成相能力，[C8mim]Ac的成相效果比[C4mim]Ac、[C6mim]Ac都要好。利用Merchuk方程對雙節(jié)線數(shù)據(jù)進行關聯(lián)，相關系數(shù)都大于0.99，說明關聯(lián)效果非常好，這為進一步試驗奠定了基礎。通過對各單因素的考察可知，[C8mim]Ac的萃取效果優(yōu)于[C6mim]Ac和[C4mim]Ac，然后釆用響應面法優(yōu)化離子液體濃度、鹽濃度、pH值這3個因素，最佳條件下分配系數(shù)可以達到11.65。與預測值11.85非常接近，說明采用響應面法優(yōu)化的提取條件可靠。