基于指紋圖譜和化學(xué)模式識(shí)別分析的車(chē)前子與炒車(chē)前子質(zhì)量評(píng)價(jià)Δ

何子驥，伍斌璽，李雨昕，沈志濱，劉奇越，王秋紅，#（.廣東藥科大學(xué)中藥學(xué)院，廣州 50006；.黑龍江中醫(yī)藥大學(xué)教育部北藥基礎(chǔ)與應(yīng)用研究重點(diǎn)實(shí)驗(yàn)室/黑龍江省中藥及天然藥物藥效物質(zhì)基礎(chǔ)研究重點(diǎn)實(shí)驗(yàn)室，哈爾濱 50040；.廣州采芝林藥業(yè)有限公司，廣州 5045）

車(chē)前子為車(chē)前科植物車(chē)前Plantago asiaticaL.或平車(chē)前Plantago depressaWilld.的干燥成熟種子[1]，其性寒，味甘，具有清熱利尿通淋、滲濕止渴、明目、祛痰的功效，主要用于治療熱淋濕痛、水腫脹滿、暑濕泄瀉、目赤腫痛、痰熱咳嗽等[2]。車(chē)前子主要含有多糖類、苯乙醇苷類、環(huán)烯醚萜類、三萜類、黃酮類、甾醇及生物堿類等化學(xué)成分[3]，具有利尿、消炎、降血糖、降血壓、調(diào)血脂、抗氧化和調(diào)節(jié)免疫等藥理作用[4]。

炒車(chē)前子為車(chē)前子炒制后的炮制品，除用于便秘的療效顯著外[5]，還可增強(qiáng)清熱利尿、滲濕通淋的功效[6]。目前，關(guān)于炒車(chē)前子的研究主要為其炮制工藝優(yōu)化[7－9]，以及以京尼平苷酸、毛蕊花糖苷含量為指標(biāo)的質(zhì)量標(biāo)準(zhǔn)研究[10]，加之關(guān)于車(chē)前子炒制前后整體化學(xué)成分變化以及炒車(chē)前子指紋圖譜的研究較少，這使得難以綜合評(píng)價(jià)車(chē)前子炒制前后的整體質(zhì)量和特征成分。

中藥指紋圖譜技術(shù)結(jié)合聚類分析、主成分分析、正交偏最小二乘回歸分析等化學(xué)模式識(shí)別分析能更加客觀、全面地反映藥材質(zhì)量，已被廣泛用于中藥材的質(zhì)量控制研究[11]。本課題組前期以2020年版《中國(guó)藥典》（一部）車(chē)前子藥材項(xiàng)下的京尼平苷酸、毛蕊花糖苷含量為指標(biāo)，對(duì)車(chē)前子及炒車(chē)前子進(jìn)行含量測(cè)定，結(jié)果顯示，車(chē)前子中京尼平苷酸的含量為0.51%～1.38%，毛蕊花糖苷的含量為0.67%～1.06%；炒車(chē)前子中京尼平苷酸的含量為0.98%～1.47%，毛蕊花糖苷的含量為0.54%～0.94%，均符合藥典規(guī)定[1]。基于此，本研究在車(chē)前子及炒車(chē)前子含量符合《中國(guó)藥典》規(guī)定的基礎(chǔ)上，采用高效液相色譜（high performance liquid chromatography，HPLC）法建立了15批車(chē)前子與15批炒車(chē)前子的指紋圖譜，同時(shí)結(jié)合化學(xué)模式識(shí)別分析方法比較車(chē)前子與炒車(chē)前子的質(zhì)量差異，旨在為完善車(chē)前子與炒車(chē)前子的質(zhì)量標(biāo)準(zhǔn)提供參考。

1 材料

1.1 主要儀器

本研究所用主要儀器有Waters Alliance 型HPLC儀（美國(guó)Waters公司），MS105DU型電子天平（瑞士Mettler Toledo 公司），Milli-Q Advantage A10 型臺(tái)式純水系統(tǒng)（德國(guó)Millipore 公司），C21-WK2102 型多功能電磁爐（廣東美的生活電器制造有限公司）等。

1.2 主要藥品與試劑

京尼平苷酸對(duì)照品（批號(hào)MB6001-S，純度＞98%）、毛蕊花糖苷對(duì)照品（批號(hào)MB6742，純度＞98%）均購(gòu)自大連美侖生物技術(shù)有限公司；甲醇為色譜純，其余試劑均為分析純，水為超純水。

15 批車(chē)前子藥材（編號(hào)S1～S15）均于2019 年購(gòu)自廣州采芝林藥業(yè)有限公司，經(jīng)廣東藥科大學(xué)中藥學(xué)院中藥資源系劉基柱教授鑒定為車(chē)前科植物車(chē)前P.asiaticaL.的干燥成熟種子。取15 批干凈的車(chē)前子藥材（編號(hào)S1～S15）適量，置于預(yù)熱電磁爐中，炒至略有爆聲并有香氣逸出時(shí)，取出放涼，即得炒車(chē)前子（對(duì)應(yīng)編號(hào)為S16～S30）。15批車(chē)前子與15批炒車(chē)前子來(lái)源信息見(jiàn)表1。

表1 15批車(chē)前子與15批炒車(chē)前子來(lái)源信息

2 方法與結(jié)果

2.1 HPLC指紋圖譜的建立

2.1.1 色譜條件 以SunFire C18（250 mm×4.6 mm，5 μm）為色譜柱，以甲醇（A）-0.5%乙酸溶液（B）為流動(dòng)相進(jìn)行梯度洗脫（0～5 min，5%A；5～15 min，5%A→20%A；15～25 min，20%A→60%A；25～30 min，60%A；30～35 min，60%A→5%A）；流速為1.0 mL/min；檢測(cè)波長(zhǎng)為254 nm；柱溫為30 ℃；進(jìn)樣量為10 μL。

2.1.2 供試品溶液的制備 取車(chē)前子或炒車(chē)前子樣品粉末（過(guò)三號(hào)篩）1.0 g，置于具塞錐形瓶中，精密稱定，加60%甲醇50 mL，精密稱定；連接冷凝管加熱回流提取2 h，放冷，再次精密稱定，用60%甲醇補(bǔ)足減失的質(zhì)量，搖勻，濾過(guò)；取續(xù)濾液，于25 ℃以13 000 r/min離心10 min，取上清液，經(jīng)0.45 μm 微孔濾膜濾過(guò)，取續(xù)濾液，即得供試品溶液。

2.1.3 對(duì)照品溶液的制備 精密稱取京尼平苷酸、毛蕊花糖苷對(duì)照品各0.1 mg，分別加60%甲醇，制成京尼平苷酸、毛蕊花糖苷質(zhì)量濃度均為0.1 mg/mL 的單一對(duì)照品溶液，經(jīng)0.45 μm微孔濾膜濾過(guò)，取續(xù)濾液，即得。

2.1.4 精密度試驗(yàn) 取“2.1.2”項(xiàng)下供試品溶液（編號(hào)S11），按“2.1.1”項(xiàng)下色譜條件連續(xù)進(jìn)樣測(cè)定6 次，以京尼平苷酸為參照峰，計(jì)算各共有峰的相對(duì)保留時(shí)間和相對(duì)峰面積。結(jié)果顯示，各共有峰相對(duì)保留時(shí)間的RSD為0.03%～0.37%（n＝6），相對(duì)峰面積的RSD 為0.56%～1.62%（n＝6），表明儀器精密度良好。

2.1.5 穩(wěn)定性試驗(yàn) 取“2.1.2”項(xiàng)下供試品溶液（編號(hào)S11），分別于室溫下放置0、2、4、6、8、16、18、24 h 時(shí)按“2.1.1”項(xiàng)下色譜條件進(jìn)樣測(cè)定，以京尼平苷酸為參照峰，計(jì)算各共有峰的相對(duì)保留時(shí)間和相對(duì)峰面積。結(jié)果顯示，各共有峰相對(duì)保留時(shí)間的RSD 為0.09%～0.41%（n＝8），相對(duì)峰面積的RSD 為0.63%～1.95%（n＝8），表明供試品溶液于室溫下放置24 h 內(nèi)穩(wěn)定性良好。

2.1.6 重復(fù)性試驗(yàn) 取車(chē)前子（編號(hào)S11），共6 份，按“2.1.2”項(xiàng)下方法制備供試品溶液，再按“2.1.1”項(xiàng)下色譜條件進(jìn)樣測(cè)定，以京尼平苷酸為參照峰，計(jì)算各共有峰的相對(duì)保留時(shí)間和相對(duì)峰面積。結(jié)果顯示，各共有峰相對(duì)保留時(shí)間的RSD 為0.05%～0.41%（n＝6），相對(duì)峰面積的RSD 為0.27%～1.75%（n＝6），表明方法重復(fù)性良好。

2.1.7 HPLC 指紋圖譜的建立 分別取15 批車(chē)前子及15 批炒車(chē)前子樣品，按“2.1.2”項(xiàng)下方法制備供試品溶液，再按“2.1.1”項(xiàng)下色譜條件進(jìn)樣測(cè)定，將得到的色譜數(shù)據(jù)導(dǎo)入《中藥色譜指紋圖譜相似度評(píng)價(jià)系統(tǒng)（2012版）》，采用中位數(shù)法，車(chē)前子以S11樣品為對(duì)照?qǐng)D譜，炒車(chē)前子以S26 樣品為對(duì)照?qǐng)D譜（因S11、S26 樣品的色譜峰峰形較好），時(shí)間窗寬度設(shè)置為0.2，得到15 批車(chē)前子及15 批炒車(chē)前子的HPLC 疊加指紋圖譜和對(duì)照?qǐng)D譜。結(jié)果顯示，15批車(chē)前子有18個(gè)共有峰，15批炒車(chē)前子有13個(gè)共有峰。經(jīng)對(duì)比發(fā)現(xiàn)，車(chē)前子與炒車(chē)前子共有8個(gè)共有峰，即車(chē)前子指紋圖譜中的峰1～峰3、峰6、峰13～峰16，炒車(chē)前子指紋圖譜中的峰1～峰3、峰7、峰10～峰13（車(chē)前子中的峰6、峰13～峰16 分別與炒車(chē)前子中的峰7、峰10～峰13對(duì)應(yīng)，為方便后文描述，將對(duì)應(yīng)的共有峰依次標(biāo)記為峰A～峰H）。15批車(chē)前子與15批炒車(chē)前子的HPLC疊加指紋圖譜見(jiàn)圖1，對(duì)照?qǐng)D譜見(jiàn)圖2。

圖1 15批車(chē)前子與15批炒車(chē)前子的HPLC疊加指紋圖譜

圖2 車(chē)前子與炒車(chē)前子的對(duì)照?qǐng)D譜

2.1.8 相似度評(píng)價(jià) 采用《中藥色譜指紋圖譜相似度評(píng)價(jià)系統(tǒng)（2012 版）》對(duì)15 批車(chē)前子及15 批炒車(chē)前子進(jìn)行相似度評(píng)價(jià)。結(jié)果顯示，30批樣品與對(duì)照?qǐng)D譜的相似度均大于0.920，表明車(chē)前子炒制前后的化學(xué)成分種類相似。結(jié)果見(jiàn)表2。

表2 15批車(chē)前子與15批炒車(chē)前子的相似度評(píng)價(jià)結(jié)果

2.1.9 共有峰的指認(rèn) 將車(chē)前子及炒車(chē)前子的對(duì)照?qǐng)D譜（圖2）與對(duì)照品的色譜圖（圖3）進(jìn)行比較，共指認(rèn)了2個(gè)共有峰，分別為車(chē)前子中的峰6（京尼平苷酸）、峰13（毛蕊花糖苷），炒車(chē)前子中的峰7（京尼平苷酸）、峰10（毛蕊花糖苷）。由于京尼平苷酸的分離度較好，含量較高，鄰近峰較少，易于區(qū)分，故選擇京尼平苷酸為參照峰，分別計(jì)算各共有峰的相對(duì)保留時(shí)間和相對(duì)峰面積。結(jié)果顯示，車(chē)前子中各共有峰相對(duì)保留時(shí)間的RSD 為0.08%～0.55%，相對(duì)峰面積的RSD為3.86%～44.09%；炒車(chē)前子中各共有峰相對(duì)保留時(shí)間的RSD 為0.08%～0.64%，相對(duì)峰面積的RSD 為12.51%～58.61%。各樣品相對(duì)保留時(shí)間的RSD 差異較小，相對(duì)峰面積的RSD差異較大，表明同一品種、不同產(chǎn)地、不同批次車(chē)前子炒制前后所含化學(xué)成分的種類具有較高的相似度，但含量存在較大差異。

圖3 京尼平苷酸與毛蕊花糖苷對(duì)照品的HPLC圖

2.2 聚類分析

以15批車(chē)前子與15批炒車(chē)前子共有峰的峰面積為變量進(jìn)行標(biāo)準(zhǔn)化處理，使用SPSS 25.0軟件，采用質(zhì)心聚類法，以平方歐氏距離為測(cè)度進(jìn)行聚類分析。結(jié)果顯示，當(dāng)歐氏距離為5時(shí)，30批樣品可聚為6類，其中S1為一類，S2～S5 為一類，S6～S15 為一類，S16～S17 為一類，S18～S20 為一類，S21～S30 為一類；當(dāng)歐氏距離為10 時(shí)，30 批樣品可聚為3 類，其中S1～S5 為一類，S6～S15、S21～S30為一類，S16～S20為一類；當(dāng)歐氏距離為16 時(shí)，30 批樣品可聚為2 類，S1～S5、S16～S20 為一類，S6～S15、S21～S30為一類。結(jié)果見(jiàn)圖4。

圖4 15批車(chē)前子與15批炒車(chē)前子的聚類分析樹(shù)狀圖

2.3 主成分分析

以15批車(chē)前子與15批炒車(chē)前子共有峰的峰面積為變量，使用SPSS 25.0軟件進(jìn)行主成分分析。結(jié)果顯示，前2 個(gè)主成分的累計(jì)方差貢獻(xiàn)率為82.575%，表明前2個(gè)主成分可以反映樣品的質(zhì)量。通過(guò)碎石圖可知，前2個(gè)主成分的曲線陡峭，后面6個(gè)主成分的曲線較平緩，表明前2個(gè)主成分在車(chē)前子及炒車(chē)前子中具有重要作用。結(jié)果見(jiàn)表3、圖5。

圖5 車(chē)前子與炒車(chē)前子的主成分分析碎石圖

表3 8個(gè)共有峰的特征值及方差貢獻(xiàn)率

主成分載荷的絕對(duì)值越大，表示主成分的變量代表性越大[12]。采用SPSS 25.0 軟件進(jìn)行主成分因子載荷矩陣。結(jié)果顯示，峰A、峰C、峰E、峰F、峰H在主成分1上具有較高的載荷值，表明這5 個(gè)峰對(duì)主成分1 的影響較大；峰D、峰G在主成分2上具有較高的載荷值，表明這2個(gè)峰對(duì)主成分2 的影響較大；峰B 在主成分1、主成分2上的載荷值均不高，表明其對(duì)主成分1、主成分2的影響較小。結(jié)果見(jiàn)表4。

表4 車(chē)前子與炒車(chē)前子的主成分因子載荷矩陣

綜合評(píng)分能夠反映樣品的整體質(zhì)量，評(píng)分越高表明樣品的整體質(zhì)量越好[13]。采用方差貢獻(xiàn)率計(jì)算15 批車(chē)前子與15批炒車(chē)前子的綜合評(píng)分。使用SPSS 25.0軟件得到主成分1（Z1）和主成分2（Z2）的評(píng)分分別為Z1＝0.445×峰A－0.326×峰B+0.347×峰C－0.165×峰D－0.433×峰E+0.418×峰F+0.266×峰G－0.339×峰H；Z2＝0.081×峰A－0.378×峰B+0.062×峰C+0.631×峰D+0.199×峰E－0.029×峰F+0.532×峰G+0.356×峰H。其中，峰A～峰H 為共有峰的峰面積，Z1、Z2的系數(shù)為主成分1、2的載荷值與特征值開(kāi)平方的比值，按公式計(jì)算綜合評(píng)分（Z），Z＝Z1×主成分1的方差貢獻(xiàn)率+Z2×主成分2的方差貢獻(xiàn)率。結(jié)果顯示，車(chē)前子的綜合評(píng)分均小于0，炒車(chē)前子的綜合評(píng)分均大于0（S24樣品除外），表明綜合評(píng)分可以較好地區(qū)分車(chē)前子與炒車(chē)前子。結(jié)果見(jiàn)表5。

表5 15 批車(chē)前子與15 批炒車(chē)前子的綜合評(píng)分結(jié)果及排序

2.4 正交偏最小二乘回歸分析

以15批車(chē)前子與15批炒車(chē)前子共有峰的峰面積為變量，采用SIMCA-P 14.1軟件進(jìn)行正交偏最小二乘回歸分析。結(jié)果顯示，車(chē)前子均位于得分圖的左側(cè)，炒車(chē)前子均位于得分圖的右側(cè)，二者能明顯分開(kāi)。模型累計(jì)解釋能力參數(shù)R2X（cum）為0.991，R2Y（cum）為0.965，模型預(yù)測(cè)度Q2（cum）為0.942，均大于0.500，表明模型構(gòu)建成功，預(yù)測(cè)能力較強(qiáng)，擬合度好[14]。結(jié)果見(jiàn)圖6。

圖6 15 批車(chē)前子與15 批炒車(chē)前子的正交偏最小二乘回歸分析得分圖

當(dāng)變量重要性投影（variable importance in the projection，VIP）值大于1 時(shí)，表示變量對(duì)整體的貢獻(xiàn)較大[15]，故以VIP 值大于1 為標(biāo)準(zhǔn)，得到貢獻(xiàn)相對(duì)較大的3個(gè)峰，分別為峰E（毛蕊花糖苷）、峰D（京尼平苷酸）、峰G，表明這3個(gè)峰對(duì)車(chē)前子與炒車(chē)前子的貢獻(xiàn)較大，推測(cè)其可能是影響車(chē)前子與炒車(chē)前子質(zhì)量差異的標(biāo)志性成分。變量的點(diǎn)距離原點(diǎn)的距離越遠(yuǎn)，表示變量對(duì)模型的影響越大[16]。結(jié)果顯示，峰E、峰D、峰G 對(duì)車(chē)前子與炒車(chē)前子的影響較大。結(jié)果見(jiàn)圖7、圖8。

圖7 8個(gè)共有峰的VIP值

圖8 15 批車(chē)前子與15 批炒車(chē)前子的正交偏最小二乘回歸分析載荷圖

3 討論

本課題組前期分別對(duì)不同流動(dòng)相（甲醇-0.1%甲酸溶液、甲醇-0.5%乙酸溶液）、不同流速（0.5、1.0 mL/min）、洗脫時(shí)間（35、45、55 min）、柱溫（30、35 ℃）等進(jìn)行考察。結(jié)果顯示，采用本研究所用的色譜條件時(shí)，各色譜峰均能較好地分離，基線相對(duì)平穩(wěn)。

指紋圖譜結(jié)果顯示，15批車(chē)前子有18個(gè)共有峰，炒車(chē)前子有13個(gè)共有峰，車(chē)前子與炒車(chē)前子共有8個(gè)共有峰，共指認(rèn)了2 個(gè)共有峰，分別為車(chē)前子中的峰6（京尼平苷酸）、峰13（毛蕊花糖苷），炒車(chē)前子中的峰7（京尼平苷酸）、峰10（毛蕊花糖苷）。15批車(chē)前子與對(duì)照?qǐng)D譜的相似度均大于0.980，15 批炒車(chē)前子與對(duì)照?qǐng)D譜的相似度均大于0.920，相似度較高，表明不同產(chǎn)地車(chē)前子、炒車(chē)前子的化學(xué)成分種類相似，具有較好的一致性。

聚類分析結(jié)果顯示，當(dāng)歐氏距離為10時(shí)，30批樣品可聚為3類，其中S1～S5為一類，S6～S15、S21～S30為一類，S16～S20 為一類；當(dāng)歐氏距離為16 時(shí)，30 批樣品可聚為2 類，S1～S5、S16～S20 為一類，S6～S15、S21～S30 為一類，表明江西省和山東省產(chǎn)車(chē)前子有一定的相似性，能與湖北省產(chǎn)車(chē)前子區(qū)分開(kāi)。主成分分析結(jié)果顯示，前2個(gè)主成分的累計(jì)方差貢獻(xiàn)率為82.575%，車(chē)前子的綜合評(píng)分均小于0，炒車(chē)前子的綜合評(píng)分均大于0（S24 樣品除外），表明炒車(chē)前子的整體質(zhì)量較好。正交偏最小二乘回歸分析結(jié)果顯示，峰E（毛蕊花糖苷）、峰D（京尼平苷酸）和峰G 可能是影響車(chē)前子與炒車(chē)前子質(zhì)量差異的標(biāo)志性成分，其可能與地理環(huán)境、氣候條件、種植模式、采收期等有關(guān)。

綜上所述，本研究所建指紋圖譜穩(wěn)定、簡(jiǎn)便快速，結(jié)合化學(xué)模式識(shí)別分析可用于評(píng)價(jià)車(chē)前子與炒車(chē)前子的質(zhì)量。毛蕊花糖苷、京尼平苷酸和峰G所代表的成分可能是影響車(chē)前子與炒車(chē)前子質(zhì)量差異的標(biāo)志性成分。