Page 50 - 《中国药房》2025年18期

1 200 1 2 3 4 6 7 8 9 10 11 的离散程度较大，表明不同批次草豆蔻药材的质量差异
1 100 5 R
1 000 S18 较大。但样品检测数据均在95%置信区间内，表明样品
S17
900 S16
S15
响应值/mAU 700 S13 检测数据无异常。18批草豆蔻药材大致可分为3组，即
800
S14
S12
600
S11
S10
500
S1～S10为一组，S11～S13为一组，S14～S18为一组，这
S9
400
S7
300 S8
S6
S5
200 S4 一结果与CA结果一致。
100 S3
S2
0 S1
0 10 20 30 40 50 2 S5
时间/min S2 S6
6：山姜素；7：乔松素；8：小豆蔻明；11：桤木酮。 1
S13 S14 S15
图1 18批草豆蔻药材的HPLC指纹图谱 t[2] 0 S7 S4 S11 S12 S16 S18
S3
－1 S9 S8 S17
150 S1 S10
－2
125 －3
－10 －8 －6 －4 －2 0 2 4 6 8
t[1]
100
响应值/mAU 75 图4 18批草豆蔻药材的PCA得分图
50 山姜素 桤木酮 2.4.3 正交偏最小二乘法判别分析
乔松素 为明确引起不同产地药材质量差异的因素，本研究
小豆蔻明
25
运用SIMCA 14.1软件对18批草豆蔻药材共有峰峰面积
0 数据进一步分析，构建正交偏最小二乘法判别分析（or‐
0 10 20 30 40 50
时间/min thogonal partial least squares-discriminant analysis，OPLS-
图2 混合对照品溶液的HPLC图 DA）模型（图 5）。该模型的累积解释率参数 R x＝
2
2
0.865、R y＝0.704，累积预测能力参数 Q ＝0.556，均大
2
2.4 指纹图谱的化学模式识别
2.4.1 聚类分析 于 0.5，表明所构建的模型稳定、可靠、预测能力强。本
研究对所构建的 OPLS-DA 模型随机排列置换检验 200
采用SPSS 23.0软件对18批草豆蔻药材进行聚类分
2
2
次，可得 R 截距值为 0.022（＜0.3），Q 截距值为－0.486
析（cluster analysis，CA）。以 11 个共有峰的峰面积为变
量，选用质心聚类法进行分析，结果见图 3。如图 3 所 （＜0），表明所构建的模型结果可靠。
示，当欧氏距离为 5 时，18 批草豆蔻药材可聚为 3 类：我 2.0
1.5
国广东、广西产药材（即S1～S10）为第一类，我国海南产 1.0 S17 S13 S11
药材（即S14～S18）为第二类，印度尼西亚和我国云南产 0.5 0 S15 S12 S4 S7 S5 S3 S1
药材（即S11～S13）为第三类。 1.032 59＊t[2] －0.5 S16 S18 S14 S6 S2 S9
欧氏距离 －1.0 S8 S10
－1.5
0 5 10 15 20 25 －2.0
S3 －2.5
S7 －10 －8 －6 －4 －2 0 2 4 6 8
S1 1.001 18＊t[1]
S5
S8 图5 18批草豆蔻药材的OPLS-DA得分图
S10
S9
S4 变量重要性投影（variable importance in projection，
S2
S6 VIP）是筛选质量差异标志物的重要指标，可直观反映共
S15
S17 有峰对应成分间差异的权重大小 ，VIP值越大，表明该
[11]
S16
S18
S14 变量（成分）对样品质量差异的贡献度越大，一般以VIP
S12 [12]
S13 值＞1.0 为筛选重要变量的标准 。本研究共筛选出 9
S11 个 VIP 值＞1.0 的成分，依次为 4、9、3、2、7（乔松素）、8
图3 18批草豆蔻药材的CA结果
（小豆蔻明）、6（山姜素）、10、11（桤木酮）号峰对应的成分，
2.4.2 主成分分析 其 VIP 值分别为 1.190 22、1.042 30、1.040 34、1.038 90、
为进一步验证不同产地草豆蔻药材的差异性，本研 1.020 08、1.018 70、1.015 34、1.014 93、1.002 01，详见图6。
究利用SIMCA 14.1软件对18批草豆蔻药材的共有峰峰 2.4.4 综合评分
面积数据进行主成分分析（principal component analysis， 根据上述OPLS-DA结果可知，有9个成分在草豆蔻
PCA）。以特征值＞1为标准进行主成分提取，共提取到 药材的质量差异中起主要作用，故将18批草豆蔻药材中
2 个主成分，累计方差贡献率达 89.798%，表明这 2 个主 上述 9 个成分共有峰的峰面积分别按 F＝（Fx－Fmin ）/
成分能代表11个指标在样品中的大部分信息。PCA得 （Fmax－Fmin ）进行归一化处理（式中，F表示归一化处理后
2
分图见图 4。由图 4 可知，模型累积解释率参数 R x＝ 的色谱峰峰面积，Fx表示色谱峰原始峰面积，Fmin表示色
0.898（＞0.5），表明该模型拟合性较好。18 批药材之间 谱峰最小峰面积，Fmax表示色谱峰最大峰面积）；并按照
· 2264 · China Pharmacy 2025 Vol. 36 No. 18 中国药房 2025年第36卷第18期