夏枯草中16种酚类化合物含量的测定及其与抗氧化活性的相关性分析

罗 敏，何 婷，龚 磊，陈荣祥

（1.遵义医科大学药学院，贵州遵义 563000；2.遵义医科大学附属医院药剂科，贵州遵义 563000；3.遵义医科大学生命科学研究院，贵州遵义 563000；4.遵义市理化分析测试工程技术研究中心，贵州遵义 563000）

夏枯草（Prunella vulgarisL.）是唇形科夏枯草属多年生草本植物，广泛分布于东亚和欧洲，在我国有悠久的食用和药用历史[1-2]。夏枯草的主要药用和食用部位为其干燥果穗，具有抗氧化、增强免疫、抗病毒、抗炎、降血压和抗肿瘤等作用[3-4]，可作为糖尿病的保健食品和药品[5]，也可用作汤料和凉茶[6]，其花蜜经充分酿造后具有很好的抗结肠炎和调节肠道菌群的作用[7]。夏枯草含有丰富的酚酸、黄酮、三萜酸等化学成分[2,8]，其中酚类化合物是夏枯草重要的活性成分[9]。酚类化合物可抵御紫外线或病原菌的侵袭，具有很好的抗氧化活性，可清除自由基对人体生命造成的损伤[10]，目前已经作为天然的抗氧化剂被广泛应用于食品方面[11]。夏枯草中的酚类化合物主要有迷迭香酸、丹参素、咖啡酸等，在《中国药典》2020年版中，迷迭香酸就被作为夏枯草含量测定的质量评价指标[12]。近年的研究也多采用高效液相色谱法测定酚类物质的含量，但是所能测定的指标多集中于迷迭香酸、丹参素、咖啡酸等少数几种化合物[13-14]。

相对传统的液相色谱-紫外检测法，液相色谱-串联质谱法（liquid chromatography-tandem mass spectrometry, LC-MS/MS）具有分析速度快、分离效能高、灵敏度高等优点，现被广泛应用于各种复杂化合物的分析，且测定低限可达纳克级[15-17]。如Ozdal等[18]采用LC-MS/MS法研究了土耳其蜂胶提取物中的多酚组成，并鉴定出32种酚类化合物；Schulz等[19]采用LC-MS/MS法定量了黄色番石榴果实酸水解后的23种酚类物质，不仅具有极高的灵敏度，而且还可提供常规方法难以获取的结构信息；罗弘杉等[20]采用液相色谱-高分辨质谱鉴定夏枯草中的酚类化合物，具有较高的分辨率，但未对多组分进行定量分析。

在色谱分离的基础上，将化学成分和生物活性成分相结合进行相关性分析，可明确不同成分对于生物活性的贡献。灰色关联分析（grey correlation analysis, GRA）和偏最小二乘回归（partial least squares regression, PLSR）分析是较常用的多因素统计方法，GRA将因素间发展的相似性程度作为衡量系统中各因素间的关联程度，即灰色关联度（gray relational degree, GRD），现已被广泛用于解决多因素、多变量之间复杂相互关系的问题；PLSR是一种对主成分分析和多元回归的特征进行泛化和组合的方法，能有效地解决变量间的多重共线性问题[17,21]。目前对于各种食品的品质调控多采用PLSR或GRA进行评价，并以此计算食品品质的影响程度，进而分析影响因素与综合评分的相关性[22-23]。

因此，本研究采用LC-MS/MS对夏枯草中16种酚类化合物进行定量分析，采用DPPH、ABTS自由基清除能力和铁离子还原能力（ferric reducing antioxidant power, FRAP）对夏枯草的抗氧化活性进行评价，并进一步采用GRA、PLSR法对24批夏枯草中16种酚类化合物的相关性进行分析，以此探讨夏枯草抗氧化的物质基础，为全面提高夏枯草的质量控制水平及其产品开发和应用提供参考。

1 材料与方法

1.1 材料与仪器

夏枯草样本分别采自广东、河南等9个省份，其详细信息见表1。样本经遵义医科大学顾丁副教授鉴定为唇形科植物夏枯草（Prunella vulgarisL）的干燥果穗。

表 1 24批夏枯草样品信息Table 1 Information of 24 batches of Prunella vulgaris

对羟基苯甲酸、异槲皮苷、七叶内酯 Ark Pharm公司；对羟基肉桂酸、2,5-二羟基苯甲酸 Adamasbeta公司；迷迭香酸、绿原酸、咖啡酸、3,4-二羟基苯甲酸、金丝桃苷、阿魏酸、水杨酸、芦丁 上海阿拉丁生化科技股份有限公司；异迷迭香酸苷 成都麦德生科技有限公司；丹参素 辰光生物有限公司；3,4-二羟基苯甲醛 北京伊诺凯科技有限公司；所有对照品的纯度≥97%；乙腈 质谱纯，霍尼韦尔贸易有限公司；室验用水 屈臣氏；其余试剂为分析纯。

I-class+XevoTQ-S超高效液相色谱-三重四极杆质谱仪（配MassLynx V.4.1软件） 美国Waters；SB-5200DT型超声波清洗机 宁波新芝生物科技；ME104E电子天平 美国梅特勒-托利多。

1.2 实验方法

1.2.1 样品制备 将夏枯草粉碎，过60目筛。准确称取1 g样品粉末，加40 mL 80%的甲醇在室温条件下超声提取30 min，提取完成后将提取液置于-20 ℃条件下静置1 h，吸取上清液，用0.22 μm微孔滤膜过滤后进行色谱分析和抗氧化能力测定。

1.2.2 对照品配制 精密称取16种夏枯草对照品，分别加入甲醇溶解定容至10 mL，配制成浓度为1 mg/mL的标准品溶液，-20 ℃保存备用。使用时用80%甲醇将其稀释成所需的混合对照品工作溶液。

1.2.3 LC-MS/MS法测定16种酚类化合物

1.2.3.1 液相色谱条件 采用Waters ACQUITY UPLC BEH C18（2.1 mm×100 mm，1.7 μm）色谱柱，采用0.1%甲酸（A）-乙腈（B）作为流动相在体积流量为0.4 mL/min下进行梯度洗脱，洗脱程序为0～15 min，95%～70% A；15～19.5 min，70%～50% A；19.5～20.5 min，50%～5% A；20.5～21.5 min，5% A；21.5～22 min，5%～95% A；22～25 min，95% A。柱温为45 ℃；进样器温度为12 ℃；进样量为1 μL。

1.2.3.2 质谱条件 采用电喷雾离子源在负离子及多反应监测（multiple reaction monitoring, MRM）模式下检测16种酚类化合物，并采用两组离子对分别对16种酚类化合物进行定量和定性分析，采用外标法进行定量。毛细管电压为3 kV，蒸发温度为500 ℃、气流量为750 L/Hr。

1.2.4 抗氧化活性

1.2.4.1 铁离子还原能力 取醋酸盐缓冲溶液（pH3.6）、10 mmol/L TPTZ溶液、20 mmol/L FeCl3溶液加入到600 μL稀释后的样品溶液中混匀，暗反应30 min后于593 nm波长处测定样品吸光度。空白组中样品溶液的吸光度测定采用80%的甲醇代替。以Trolox为标准品绘制标准曲线，线性回归方程和相关系数分别为Y=0.0258X+0.0645和r=0.9997，其中，Y为吸光度（Au），X为Trolox浓度（mg/L）。所有样品测量结果均以Trolox当量表示（mg TE/g DW）。

1.2.4.2 ABTS自由基清除能力 ABTS工作液的配制：准确称取ABTS试剂0.0282 g，用蒸馏水溶解定容至7.5 mL，再加入7.5 mL 2.45 mmol/L过硫酸钾水溶液，在室温下避光放置8～12 h后用纯水稀释20倍后备用。取200 μL稀释100倍后的样品溶液加入到0.4 mL ABTS工作液中，暗反应30 min后于734 nm波长处测定样品吸光度。空白组中样品溶液的吸光度测定以80%甲醇代替，以Trolox为标准品绘制标准曲线，线性回归方程和相关系数分别为Y=-0.00901X+0.5311和r=0.9998，其中Y为吸光度（Au），X为Trolox浓度（mg/L）。所有样品测量结果的表示同“1.2.4.1”项。

1.2.4.3 DPPH自由基清除能力 取200 μL稀释30倍后的样品溶液加入到400 μL 0.14 mol/L的DPPH溶液中，混匀，暗反应30 min后于517 nm波长处测定样品吸光度。空白组中样品溶液的吸光度测定以及所有样品测量结果的表示同“1.2.4.1”项。以Trolox为标准品绘制标准曲线，线性回归方程和相关系数分别为Y=-0.0155X+0.8361和r=0.9993。其中Y为吸光度（Au），X为Trolox浓度（mg/L）。

1.2.5 相关性分析

1.2.5.1 GRA 以16种酚类化合物含量和DPPH、ABTS自由基清除能力和FRAP效应值分别作为自变量和因变量。在判别系数设为0.5的情况下，计算GRD值，GRD越高，相应的酚类化合物抗氧化活性越大。

1.2.5.2 PLSR 采用16种酚类化合物的含量作为自变量，DPPH、ABTS和FRAP值作为因变量。将数据导入软件SIMCA 14.1中进行PLSR分析，得出各自的CoeffCS值和VIP值，以分析自变量与因变量的关系。

1.3 数据处理

采用Waters公司的MassLynx V4.1软件进行数据采集；采用Excel进行数据分析；采用Origin 2018版以及SIMCA 14.1版软件进行绘图。

2 结果与分析

2.1 夏枯草中16种酚类成分质谱条件优化及LCMS/MS色谱图

对16种化合物的标准品进行锥孔电压、碰撞能量等条件的优化，优化结果见表2。标准品图见图1A，除丹参素、迷迭香酸和异迷迭香酸苷浓度为5 mg/L外，其余化合物浓度均为0.5 mg/L。夏枯草样品图见图1B。

图 1 混合对照品（A）和样品（B）提取离子流色谱图Fig.1 The extract ion chromatograms of mixed reference (A)and testing sample (B)

2.2 方法学验证

2.2.1 线性关系 分别精密吸取“1.2.2”项下的对照品溶液配制成不同浓度的混标，按“1.2.3”项下的色谱条件进样分析，记录测定后16种酚类化合物的峰面积。以对照品的浓度（mg/L）为横坐标，各待测成分色谱峰峰面积为纵坐标绘制标准曲线，计算回归方程。结果显示，16种酚类成分在各自的线性范围内，线性关系良好（r＞0.99）。检出限（S/N=3）在2.80～9.47 ng/mL之间，其回归方程以及相关系数见表3。

表 2 16种主要成分的保留时间和质谱参数Table 2 Retention time and mass parameters of the 16 components

2.2.2 重复性 称取编号为S14的夏枯草粉末约1 g，按“1.2.1 样品制备”方法平行制备6份供试品溶液，按上述色谱和质谱条件，测定并记录夏枯草中16种化合物的峰面积，计算RSD值，结果显示RSD值≤5.71%。

2.2.3 精密度 精密量取标准品储备溶液并将其配制成0.5 mg/L的工作液（其中，异迷迭香酸苷、丹参素和迷迭香酸浓度为5 mg/L），按上述质谱和色谱条件，连续进样6次，进样量1 μL，测定并记录16种化合物的峰面积，计算RSD值，结果显示RSD值≤2.24%。

2.2.4 稳定性 称取编号为S14的夏枯草粉末约1 g，按“1.2.1 样品制备”方法制备供试品溶液，按上述测试条件在24 h（0、4、8 、12、16、24 h）内测定16种化合物的稳定性，记录夏枯草中16种化合物的峰面积变化并计算RSD值。结果显示24 h内RSD值≤5.06%。

2.2.5 加样回收实验 取同一批含量已知的夏枯草（S6）粉末3份，每份约1 g，按照含量1:1准确加入16种化合物，按“1.2.1 样品制备”方法制备供试品溶液，按上述测试条件进行测定，考察回收率，平均加标回收率在95.31%～103.74%之间，其RSD值≤6.24%。方法学考察结果见表3。

2.3 样品含量测定

分别称取夏枯草样品粉末，按“1.2.1 样品制备”方法制备供试品溶液，进样，按上述测试条件测定夏枯草中16种酚类化合物并记录各种化合物的峰面积，计算其含量，结果见表4。

24批样品中化合物的含量差异很大。其中，迷迭香酸、异迷迭香酸苷含量最高，丹参素次之，其含量分别为205.97～1036.73 mg/100 g、44.73～120.67 mg/100 g和16.24～59.24 mg/100 g。刘月新等[24]的研究表明，夏枯草在不同生长阶段其化学成分含量不同，迷迭香酸和异迷迭香酸苷在枯萎期含量明显增加。刘光敏等[25]的研究指出同一品种夏枯草药材因产地不同其酚类成分含量也有较大差异。因此，样品间的差异不仅与产地有关，而且受采收时间影响。

2.4 抗氧化活性研究

DPPH、ABTS自由基清除能力和铁离子还原能力结果见表5。在24批夏枯草样品中，DPPH、ABTS自由基清除能力和铁离子还原能力的范围分别为7.53～43.83 mg TE/g DW、10.75～34.93 mg TE/g DW和8.53～47.76 mg TE/g DW。其中，抗氧化能力最好的夏枯草样品为S6，较差的样品为S7和S24。

2.5 夏枯草中16种酚类化合物与抗氧化能力的关联性分析

2.5.1 GRA分析 利用GRA分析夏枯草中16种酚类化合物含量与抗氧化活性之间的相关性，以评价这些成分的活性贡献。GRD值用来描述成分与抗氧化活性之间的关系，GRD大于0.6说明化合物与抗氧化活性之间具有相关性，GRD大于等于0.8说明相关性较强[26-27]。16种酚类化合物与抗氧化活性相关性的GRD值如表6所示，16种化合物的GRD值均大于0.6，可见这些化合物均和抗氧化活性相关。其中，迷迭香酸与DPPH、ABTS、FRAP的关联度均大于0.8，芦丁和DPPH的关联度大于0.8，咖啡酸、芦丁和FRAP的关联度大于0.8，它们具有相对较高的GRD值可以认为是夏枯草中关键的抗氧化成分。

表 3 方法学考察Table 3 Method validation of the proposed method

表 4 24批夏枯草样品中16种酚类化合物含量测定结果（mg/100 g）Table 4 The contents of 16 phenolic compounds in 24 batches of Prunella vulgaris (mg/100 g)

2.5.2 PLSR分析 采用PLSR分析DPPH、ABTS、FRAP法与夏枯草中16种酚类化合物的相关性，以16种酚类化合物含量为X变量，以DPPH、ABTS、FRAP值为Y变量进行相关性分析，结果见图2。CoeffCS系数的正负表示不同化合物同抗氧化活性正相关或负相关，VIP值代表相关性强弱。从结果可以看出，对3种抗氧化活性均呈正相关且VIP值大于1.0的化合物有迷迭香酸、芦丁、金丝桃苷、异槲皮苷、咖啡酸、水杨酸、2,5-二羟基苯甲酸。

图 2 16种化合物的偏最小二乘回归相关性分析CoeffCS （A）和VIP （B）图Fig.2 Correlation coeffCS (A) and VIP (B) values of the 16 compounds analyzed by PLSR

大量研究表明，酚类化合物是食用/药用植物抗氧化的主要来源[28-29]。本研究结合GRA与PLSR相关性分析可知，夏枯草中16种酚类化合物均具有抗氧化活性，GRD值均大于0.6，且GRA与PLSR相关性分析均表明迷迭香酸抗氧化能力最好。

表 5 24批夏枯草抗氧化能力结果Table 5 The antioxidant capacity in 24 batches of Prunella vulgaris

表 6 GRA结果Table 6 The results of GRA

3 结论

综上所述，本研究成功运用LC-MS/MS法同时测定夏枯草中16种酚类化合物的含量，该方法线性关系良好（r＞0.99），灵敏度较高。24批夏枯草中不同酚类化合物含量在0.02～1036.73 mg/100 g之间，其中，迷迭香酸含量最高，为205.97～1036.73 mg/100 g。进一步采用GRA和PLSR分析了化学成分与DPPH、ABTS自由基清除能力和铁离子还原能力的相关性。GRA结果表明16种酚类化合物与抗氧化活性的关联系数均大于0.6，均与抗氧化活性相关。PLSR的结果表明迷迭香酸、芦丁、金丝桃苷、异槲皮苷、咖啡酸等化合物对抗氧化活性的贡献较大。本研究不仅测定了夏枯草中主要化合物的含量，而且明确了主要化学成分对抗氧化的贡献，可为夏枯草质量控制与开发利用提供参考。