Page 61 - 《中国药房》2022年10期
P. 61

250                                            以组间连接法进行聚类分析。结果显示,当距离为 20
            200                                            时,13 批 QTXT 可聚为 3 类,其中 S2 为一类,S1、S9、S10
                                                           为一类,S3~S8、S11~S13 为一类,表明不同批次样品
            150
           mAU             6                               的化学成分存在差异。结果见图4。
            100          4
                                   10
                                                           2.5 主成分分析
             50                8
                          5        9                           以 17 个共有峰峰面积为变量,采用 SPSS 22.0 软件
                                              14
             0
                                                           进行主成分分析。结果显示,前3个成分的特征值均大
              0   5  10  15  20  25  30  35  40  45  50  55  60
                                 t/min                     于1,累计方差贡献率为85.120%,可涵盖QTXT共有峰
                             A.混合对照品溶液                     中 85.120%的信息,具有一定的代表性。碎石图显示,
            650
                           6                               当特征值<1 时,曲线趋于平缓,表明提取特征值>1 的
            520
                                                           成分合理。主成分因子载荷矩阵表明,4、7、10~13号峰
            390          4                                 在主成分1中具有较高的载荷值;1、2、5、6、15号峰在主
           mAU  260                 10                     成分 2 中具有较高的载荷值;7~9、14、16、17 号峰在主
            130                                            成分3中具有较高的载荷值。结果见表1、图5、表2。
                          5     8  9                           根据下式计算 3 个主成分的主成分得分(F1、F2、
             0                                14
                                                           F3 ):F=ZX×A(式中,F表示主成分得分,ZX表示原始数
              0      10    20    30     40    50     60
                                 t/min                     据标准化后得到的标准化矩阵,A表示标准化的正交特
                           B.供试品溶液(编号S5)
                                                           征向量矩阵) ;以 3 个主成分对应的方差贡献率为权
                                                                       [19]
            650
                                                           重,得到综合得分(F 综 ):F 综=0.344 83×F1+0.279 00×F2+
            520
                                                                                                 [19]
                                                           0.227 37×F3,综合得分越高,表明质量越好 。结果显
            390                                            示,不同批次 QTXT 样品质量存在差异,S2、S3、S5 综合
           mAU  260                                        评分较高,表明这3批样品的质量较好。结果见表3。
            130                                                将17个共有峰峰面积导入SIMCA 14.1软件得到主
                                                                                                      2
             0                                             成分分析得分图。结果显示,模型解释率参数(R X)为
                                                                                    2
              0   5  10  15  20  25  30  35  40  45  50  55  60  0.983,模型预测能力参数(Q )为0.621,均大于0.5,表明
                                 t/min                                     [21]
                               C.空白溶液                      该模型稳定、可靠 。13 批样品可分为 4 类,其中 S2 为
           4:柚皮苷;5:橙皮苷;6:新橙皮苷;8:盐酸黄连碱;9:盐酸巴马汀;             一类,S1为一类,S9、S10为一类,S3~S8、S11~S13为一
        10:盐酸小檗碱;14:洛伐他汀                                   类。与聚类分析结果相比,主成分分析进一步将 S1 与
        图2   柚皮苷等成分的混合对照品溶液、供试品溶液、空                        S9、S10样品区分开。结果见图6。
             白溶液的HPLC图                                     2.6  正交偏最小二乘法-判别分析
        2.3.6  相似度评价      采用《中药色谱指纹图谱相似度评                      以 17 个共有峰峰面积为变量,采用 SIMCA 14.1 软
        价系统(2012版)》,以对照指纹图谱为参照,进行整体相                       件进行有监督模式识别的正交偏最小二乘法-判别分
                                                                                         2
                                                                                             2
        似度评价。结果显示,13批样品与对照指纹图谱的相似                          析。结果显示,模型解释率参数(R X、R Y)分别为0.892、
                                                                                    2
        度 分 别 0.989、0.990、0.997、0.999、0.992、0.993、0.992、   0.962,模型预测能力参数(Q )为0.719,均大于0.5,表明
                                                                           [22]
        0.992、0.988、0.987、0.997、0.997、0.998,表明不同批次         该模型稳定、可靠 。13 批 QTXT 可分为 4 类,S2为一
        QTXT中的化学成分种类相似度较高。                                 类,位于得分图右上方;S1 为一类,位于得分图右下方;
        2.3.7  共有峰归属      取各单味饮片样品溶液和各阴性                   S9、S10为一类,位于得分图右侧中间;S3~S8、S11~S13
        样品溶液,按“2.1”项下色谱条件进样测定,记录色谱图,                       为一类,位于得分图左侧。与主成分分析模型相比,不同
        将上述溶液色谱图与对照指纹图谱进行对比。结果显                            批次样品间的差异更加明显。结果见图7。
        示,7、8、9、10 号峰为黄连的专属峰;3、4、5、6、11、12、13                  变量重要性投影(variable importance in projection,
        号峰为麸炒枳实的专属峰;14号峰为红曲的专属峰;1号                         VIP)为筛选差异性标志物的重要指标,其值越大,表明
                                                                                       [21]
        峰为黄连和红曲的共有峰;2、15号峰为麸炒枳实和红曲                         对组间质量差异的影响越明显 。以 VIP 值>1 为标
                                                             [19]
        的共有峰;16、17 号峰为 6 味饮片的共有峰。结果见                       准 ,筛选影响QTXT质量的标志性成分。结果显示,有7
        图3。                                                个共有峰的VIP值>1,从大到小依次为10号峰(盐酸小
        2.4  聚类分析                                          檗碱)、9号峰(盐酸巴马汀)、5号峰(橙皮苷)、11号峰、8号
            以 17 个共有峰峰面积为变量,采用 SPSS 22.0 软件                峰(盐酸黄连碱)、12号峰、6号峰(新橙皮苷),表明这7个
        进行 Z-score 标准化处理后,利用平方欧氏距离为测度,                     成分可能是影响QTXT质量的标志性成分。结果见图8。


        中国药房    2022年第33卷第10期                                            China Pharmacy 2022 Vol. 33 No. 10  ·1207 ·
   56   57   58   59   60   61   62   63   64   65   66