玫瑰茄干花萼，选用高原自然晒干熟透紫红黑亮饱满花萼，购自云南省武定县。纤维素酶，购自北京鸿润宝顺科技有限公司，酶活力为2000 U/g。盐酸、乙醇溶剂，购自国药集团化学试剂有限公司。本实验所用试剂均为分析纯。

本研究主要试验仪器与设备见下 表1 。

将适量的玫瑰茄干花萼平铺在干净的玻璃板上，(45 ± 1)℃下烘干至恒重(冷却10 min前后质量之差不大于2 mg)。取适量干燥后的花萼放入高速多功能粉碎机中打磨成粉，过80目筛密封避光保存备用。

称取1.0 g样品于50 mL离心管中，加入适量纤维素酶，按不同的料液比倒入配制好的提取剂(10% HCL (0.1 mol/L) + 90%乙醇(75%))。把离心管放入超声波清洗机中，设定超声时间为20 min、超声功率为200 W、超声温度为30℃。同时，为了避免样品在超声过程中温度升高，在超声波清洗机里适时加冷水，保持温度恒定。超声后在4℃下，4000 r/min离心10 min，上清液即为花青素提取液。

单因素试验因素、水平设置如下：纤维素酶浓度：0、0.2、0.4、0.6、0.8 U/mL；液料比：10:1、20:1、30:1、40:1、50:1 mL/g；超声温度：30、35、40、45、50℃；超声时间：5、10、15、20、25 min；超声功率：140、160、180、200、220 W。每个试验做三次平行，取平均值。

采用响应面中心组合试验设计(Box-Behnken)，设计自变量为液料比、超声温度、超声时间和超声功率的四因素三水平试验( 表2 )，以提取得率为响应值，确定各因素的最佳水平值。

采用pH示差法测定提取花青素的提取得率 [8] 。pH1.0与pH4.5缓冲液配制参照江岩的方法 [9] 。花青素在溶液介质中存在4种结构形式，在某一pH值下处于动态平衡 [10] 。pH值改变时动态平衡被打破，因此花青素溶液用缓冲液稀释后，必须静置一段时间，建立新的动态平衡后，才能测定其吸光度。本试验静置时间为2 h。

计算公式如下：

花青素得率(mg/g) $Y=\frac{\Delta A\cdot V\cdot n\cdot M_W}{{{\displaystyle \text{av}E}}_{1\text{cm}}^{1\%}\cdot \text{m}\cdot \text{L}}$

式中：△A = (A_{520nm, pH1.0} − A_{700nm, pH1.0}) − (A_{520nm, pH4.5} − A_{700nm, pH4.5})；V为待测花青素样品提取液总体积，mL；n为测定时提取液稀释倍数； $\text{av}{E}_{1\text{cm}}^{1\%}$为矢车菊色素-3-葡萄糖苷的消光系数，26,900；M_W为矢车菊色素-3-葡萄糖苷的相对分子质量，449.2；m为玫瑰茄干花萼粉末样品质量，g；L为比色光程长度，l cm。

以蒸馏水为对照，用A_700nm来消除样液浑浊的影响。

所有试验指标的测定均重复三次，结果取平均值并计算标准误差。数据采用Excel、Origin 2019b与Design Expert 12.0软件进行分析和绘制。用SPSS 17.0进行ANOVA单因素方差分析并检验数据的差异显著性。

吸取一定体积的玫瑰茄花青素提取液稀释后，将蒸馏水作为对照，波谱扫描400~600 nm的吸光度，由 图1 可知，提取液的特征吸收峰在520 nm左右处出现，这与师文添 [11] 、陆国权等 [12] 的研究结果一致，说明该提取液中的色素有花青素特征，为花青素类物质。

纤维素酶浓度对玫瑰茄花青素得率的影响如 图2 所示。添加纤维素酶处理的花青素提取得率显著高于未添加对照处理(P < 0.05)，最高达到9.58 mg/g，比对照提高了13.91%。随着纤维素酶浓度的升高，玫瑰茄花萼花青素提取得率不断提高，这可能是因为随着酶添加量的增大，增强了对花萼细胞壁的水解作用，底物与酶的接触更加充分，有利于花青素的释放，提高了提取得率，这与张辉青花等的研究一致 [6] 。当纤维素酶浓度达到0.4 U/mL后提取得率增加差异不显著(P > 0.05)，这可能是因为底物的量有限，此时反应已经达到饱和 [13] 。因此，从节约成本、提高生产经济效益考虑，选择纤维素酶浓度为0.4 U/mL。

液料比对花青素得率的影响如 图3 所示。在10:1~40:1之间时，花青素得率与液料比呈正相关性，在40:1时达到最大值，为9.62 mg/g。但方差分析结果表明，在20:1~40:1之间时，提取得率并无显著性增加(P > 0.05)。这可能是提取液体积较小时，整个提取反应体系较小，溶质浓度过高，流动性差，致使底物与酶的接触不充分、不彻底，提取得率较低。当提取液体积增加时，整个提取反应体系流动性增强，反应更为彻底，提取得率随之升高。但当提取液体积为样品质量数值的50倍时，提取得率反而降低，可能是因为反应体系过大，纤维素酶和底物的浓度降低，酶解效果并不理想，因此影响了提取得率 [14] 。考虑到便于后期的纯化和浓缩，选择液料比为20 mL/g。

超声温度(酶解温度)对玫瑰茄花青素提取得率的影响如 图4 所示。在30℃~50℃范围内，玫瑰茄花萼花青素提取得率随着温度的升高显著地呈现先升高后降低的趋势，在35℃达到最大值，为9.41 mg/g。可能是因为随着温度的升高，纤维素酶活性增强，酶促反应加剧，有利于花青素的提取。但超声温度(酶解温度)持续升高，提取得率会有所下降，可能有两方面的原因：一方面，因温度过高，纤维素酶活性可能降低或失活，酶解速率下降，导致花青素提取得率下降；另一方面，超声温度过高，可能会引起部分花青素的分解。因此，选择超声温度(酶解温度)为35℃。

采用纤维素酶耦合超声波提取，超声时间与酶解时间相同，超声时间(酶解时间)对玫瑰茄花青素提取得率的影响如 图5 所示。在超声时间(酶解时间)为5~20 min时，随着反应时间的增加，花青素提取得率逐渐显著性升高(P < 0.05)，超声时间(酶解时间)在20 min时，提取得率达到最高，为9.92 mg/g。这可能是刚开始反应时，由于时间短，纤维素酶与超声波耦合作用还未充分发挥，花青素未能释放完全。当反应时间相应有所延长后，耦合作用效果进一步增强，反应更加完全彻底，因此提取得率显著性提高(P < 0.05)。但20 min后，如果继续超声(酶解)，花青素得率反而下降。很有可能是因为在反应已经达到饱和之后，继续超声，导致花青素的结构被影响了，或者是发生了降解，导致提取效果不佳。因此，选择超声时间(酶解时间)为20 min。

适当的超声处理有利于细胞内物质的浸出 [15] [16] 。超声功率对玫瑰茄花青素提取得率的影响如 图6 所示。玫瑰茄花萼花青素提取得率随着超声功率的升高，走向为先升高后降低，在160 W时达到最大值，为9.61 mg/g。这可能是因为超声可以增强细胞的空化作用并通过机械振动使分子运动加剧，细胞壁破碎程度加大，促进花青素溶出，同时由于颗粒尺寸减小和催化体系表面积增加，也可提高花青素的提取得率。当超声功率过小时，细胞壁的破碎程度较小，提取得率不高。当超声功率过大时，超强声作用导致花青素构象发生改变甚至破坏，从而降低了提取得率。因此，选择超声功率为160 W。

利用Design Expert软件对试验数据进行多元回归拟合，得到液料比(A)、超声温度/酶解温度(B)、超声时间/酶解时间(C)、超声功率(D)的多元回归模型为：Y = 11.38 + 0.4883A + 0.0508B + 0.1175C + 0.045D + 0.0475AB − 0.0675AC − 0.08AD + 0.0625BC + 0.0175BD − 0.0175CD − 0.3411A² − 0.4598B² − 0.5848C² − 0.3461D²。响应值结果见 表3 ，方差分析结果见 表4 。

续表

由 表4 可知，玫瑰茄花青素提取得率的整体模型方程的F值为68.11，具有极显著性(P < 0.01)，表明回归方程的因变量与自变量之间存在明显的线性关系，该试验方法可靠。失拟项F值为4.15，不具有显著性(P > 0.05)，同样证明了选择这个模型是比较合适的。方程的相关系数R²为0.9855，表示这个模型与试验的拟合程度较高。方程的校正决定系数R²_Adj为0.9711，预测R²_Pre为0.922，两者之差为0.0491，小于0.2，表明该模型准确解释数据的能力良好。因此，该模型方程可以用来分析与预测不同提取条件下玫瑰茄花青素得率的工艺条件。由F值分析结果判断各因素对玫瑰茄花青素得率的影响力为：A (液料比) > C (超声时间/酶解时间) > B (超声温度/酶解温度) > D (超声功率)。

注：P < 0.05，差异显著；P < 0.01，差异极显著。

利用Design Expert软件对试验数据进行分析，通过对6组响应面及等高线图( 图7 )的形状分析4个因素间的交互作用。

图7. 响应面及等高线图

进一步对拟合的回归方程进行分析可得，最优条件为：液料比27.19:1 mL/g，超声温度35℃，超声时间20.32 min，超声功率159.64 W，玫瑰茄花青素得率最高为11.56 mg/g。为了方便实际试验操作，将条件调整为：液料比27:1 mL/g，超声温度35℃，超声时间20 min，超声功率160 W。根据调整后的条件，花青素得率高达11.38 mg/g，把这个数据与模型预测值的11.56 mg/g相比，只有1.58%的相对误差。证明在此模型下得到的最优条件可靠性较高，可以作为纤维素酶耦合超声波提取玫瑰茄花青素工艺的回归分析和参数优化。