灰色關(guān)聯(lián)分析法評價不同產(chǎn)地牡丹皮藥材質(zhì)量

王斌 梁偉龍 林欽賢

【摘 要】 目的：建立不同產(chǎn)地牡丹皮藥材的質(zhì)量評價方法。方法：以17個產(chǎn)地的牡丹皮藥材為研究對象，測定其沒食子酸、兒茶素、氧化芍藥苷、芍藥苷、苯甲酸、丹皮酚、苯甲酰芍藥苷、浸出物、總灰分的含量，結(jié)合灰色關(guān)聯(lián)分析法，建立牡丹皮藥材質(zhì)量的灰色關(guān)聯(lián)度評價模型。結(jié)果：17個產(chǎn)地牡丹皮藥材的相對關(guān)聯(lián)度在0.314～0.597之間，山西運城、安徽銅陵、安徽南陵產(chǎn)地牡丹皮樣品的相對關(guān)聯(lián)度排名靠前，藥材質(zhì)量較優(yōu)，湖北恩施、重慶墊江、安徽亳州、山東菏澤產(chǎn)地牡丹皮樣品的相對關(guān)聯(lián)度排名靠后，藥材質(zhì)量較差。結(jié)論：研究結(jié)果可為牡丹皮資源的開發(fā)與利用提供參考依據(jù)。

【關(guān)鍵詞】 牡丹皮;不同產(chǎn)地;灰色關(guān)聯(lián)分析法;相對關(guān)聯(lián)度;質(zhì)量評價

Abstract：Objective To establish the quality evaluation method of cortex moutan from different producing areas.Methods Seventeen samples of cortex moutan from differenct producing areas were collected，the contents of gallic acid，catechin，oxidation paeoniflorin，paeoniflorin，benzoic acid，paeonol，benzoyl-paeoniflorin，extract and total ash of cortex moutan were determined.The grey correlation degree model for evaluating the quality of cortex moutan was established by grey correlation analysis.Result The relative correlation degrees of all evaluation samples were in the range of 0.314～0.597.The samples of cortex moutan produced in Yuncheng City of Shanxi Province，Tongling City of Anhui Province and Nanling City of Anhui Province ranked the top few in the relative correlation degrees，which had better quality.The samples of cortex moutan from Enshi City of Hubei Province，Dianjiang County of Chongqing，Bozhou City of Anhui Province and Heze City of Shandong Province ranked the bottom few in the relative correlation degrees，indicating that the samples produced in these areas had poor quality.Conclusion The results of this study can provide a reference basis for the exploitation and utilization of this resources.

Keywords：Cortex Moutan;Different Producing Areas;Grey Correlation Degree Method;Relative Correlation Degrees;Quality Evaluation

牡丹皮為毛茛科植物牡丹Paeonia suffruticosa Andr.的干燥根皮，味苦、辛，微寒，歸心、肝、腎經(jīng)，具有清熱涼血、活血化瘀的功效[1]，主產(chǎn)于安徽、山東、河南、湖南、湖北、四川、重慶等地。牡丹皮中主含酚及其苷類、單萜及其苷類兩大類，亦含有機酸、三萜、黃酮、香豆素、揮發(fā)油、多糖及微量元素等多種成分[2-7]。

灰色關(guān)聯(lián)分析法為灰色系統(tǒng)理論的重要組成部分，主要針對信息不完全、機制不明確的灰色系統(tǒng)，根據(jù)數(shù)據(jù)序列曲線的相似程度分析因子間的量化和序化，實現(xiàn)對系統(tǒng)的預(yù)測、決策等功能[8-11]，目前已在許多中藥研究中應(yīng)用，如基于灰色關(guān)聯(lián)分析方法評價甘肅產(chǎn)紅芪質(zhì)量[12]、牡丹皮不同采收期質(zhì)量綜合評價研究[13]，采用灰色關(guān)聯(lián)度法綜合評價不同采收期的牡丹皮藥材質(zhì)量，取得了許多科技成果。為了客觀、全面地評價不同產(chǎn)地牡丹皮藥材的質(zhì)量，本實驗以17個不同產(chǎn)地牡丹皮藥材為研究對象，按《中華人民共和國藥典》2015版方法測定浸出物和總灰分，采用HPLC法同時測定其沒食子酸、兒茶素、氧化芍藥苷、芍藥苷、苯甲酸、丹皮酚和苯甲酰芍藥苷含量，結(jié)合灰色關(guān)聯(lián)分析法，建立灰色關(guān)聯(lián)度模型，為牡丹皮藥材質(zhì)量評價提供參考依據(jù)。

1 儀器與材料

1.1 儀器 DIONEX Ultimate 3000高效液相色譜儀（賽默飛世爾科技有限公司）、MX5型、XS-204型電子天平（Mettler Toledo）、DHG-9245A電熱鼓風(fēng)干燥箱、HWS-26型電熱恒溫水浴鍋（上海一恒科學(xué)儀器有限公司）、SXL-1208型程控式箱式電爐（上海精宏實驗設(shè)備有限公司）

1.2 材料 沒食子酸（110831-201605，純度90.8%）、兒茶素（110877-201604，純度99.2%）、芍藥苷（110736-201539，純度96.4%）、丹皮酚（110708-201407，純度99.9%）和苯甲酸（100419-200301）購自中國食品藥品檢定研究院;氧化芍藥苷（MUST-16021505，純度99.93%）和苯甲酰芍藥苷（MUST-16061803，純度99.28%）購自成都曼斯特有限公司;HPLC級甲醇（Merck）和磷酸（阿拉丁）;其余試劑為分析純，超純水。牡丹皮藥材于2017年先后于安徽、山東、湖南、湖北、重慶、四川等牡丹皮藥材主產(chǎn)區(qū)收集，經(jīng)廣州市香雪制藥股份有限公司康志英高級工程師鑒定為毛茛科植物牡丹Paeonia suffruticosa Andr.的干燥根皮，樣品來源信息見表1。

2 方法與結(jié)果

2.1 浸出物 按《中華人民共和國藥典》2015年版四部醇溶性浸出物測定法（通則2201）項下的熱浸法測定，用乙醇作溶劑。

2.2 總灰分 按《中華人民共和國藥典》2015年版四部總灰分測定法（通則2302）進行測定。

2.3 多成分含量測定[14]

2.3.1 對照品溶液制備 取對照品沒食子酸、兒茶素、氧化芍藥苷、芍藥苷、苯甲酸、丹皮酚和苯甲酰芍藥苷適量，精密稱定，加50%甲醇水分別制成每1mL含沒食子酸0.5μg、0.5μg、15μg、20μg、20μg、50μg、50μg的混合對照品溶液，即得。

2.3.2 供試品溶液制備 取牡丹皮藥材粉末（過三號篩）約0.5g，精密稱定，置具塞錐形瓶中，精密加入甲醇50mL，稱定重量，加熱回流提取1h，放冷，再稱定重量，用甲醇補足減失的重量，搖勻，濾過。精密移取續(xù)濾液1mL，置10mL量瓶中，加50%甲醇至刻度，搖勻，過0.45μm的濾膜，即得。

2.3.3 色譜條件 Agilent ZORBAX Eclipse XDB-C18色譜柱（250mm × 4.6mm，5μm）;檢測器為DAD檢測器;檢測波長為230nm，258nm（氧化芍藥苷）;流動相為甲醇（A）-0.1%磷酸水溶液（B），梯度洗脫（0～15min，10%A→32%A;15～25min，32%A→48%A;25～40min，48%A→58%A）;流速為1mL·min-1;進樣量：10μL;柱溫：30℃。在此條件下，取對照品溶液、供試品溶液分別進樣10μL分析，HPLC色譜圖如圖1所示。

2.3.4 線性關(guān)系考察 精密吸取“2.3.1”項下的混合對照品溶液1、3、5、7、10、15μL，在“2.3.3”項色譜條件下分析，以對照品的量為橫坐標(biāo)（X），峰面積為縱坐標(biāo)（Y）進行線性回歸，結(jié)果見表2。

2.3.5 精密度試驗 精密吸取“2.3.1”項下的混合對照品溶液，連續(xù)進樣6次，記錄7個對照品峰面積，結(jié)果顯示7個成分的峰面積RSD均小于1.5%，表明儀器精密度良好。

2.3.6 穩(wěn)定性試驗 按“2.3.2”項下方法制備供試品溶液，在“2.3.3”項色譜條件下分別于0、3、6、9、12、18、24h進樣分析，記錄7個成分的峰面積，結(jié)果顯示7個成分峰面積RSD均小于2.0%，表明供試品溶液在24h內(nèi)基本穩(wěn)定。

2.3.7 重復(fù)性試驗 精密稱定同一樣品粉末6份，按“2.3.2”方法平行制備供試品溶液，進樣分析，記錄峰面積，計算含量。結(jié)果顯示7個成分的含量RSD均小于2.0%，表明該方法重復(fù)性良好。

2.3.8 加樣回收率試驗 精密稱取已知含量的牡丹皮粉末約0.25 g，共6份，加入與7個成分等量的對照品，按照“2.3.2”項下方法平行制備供試品溶液，計算加樣回收率，結(jié)果見表3。

2.4 各指標(biāo)成分含量測定 取不同產(chǎn)地牡丹皮藥材樣品，分別按“2.1”、“2.2”、“2.3”項下方法測定浸出物、總灰分、以及沒食子酸、兒茶素、氧化芍藥苷、芍藥苷、苯甲酸、丹皮酚和苯甲酰芍藥苷的含量，結(jié)果見表4。

2.5 灰色關(guān)聯(lián)度分析法的建立

2.5.1 參考序列選擇 設(shè)有n個研究樣品，每個樣品有m項評價指標(biāo)，組成評價單元序列，記為{Xik}（i=1，2，3……n;k=1，2，3……m;本研究中n=17，m=9）。為了統(tǒng)一藥材質(zhì)量評價標(biāo)準(zhǔn)，采用將低指標(biāo)高優(yōu)化的方法，本試驗將總灰分進行取倒數(shù)轉(zhuǎn)換，與沒食子酸、兒茶素、氧化芍藥苷、芍藥苷、苯甲酸、丹皮酚、苯甲酰芍藥苷和浸出物指標(biāo)共同組成指標(biāo)矩陣。藥材質(zhì)量評價指標(biāo)統(tǒng)一后，將最優(yōu)參考序列的各指標(biāo)對應(yīng)n個不同產(chǎn)地樣品對應(yīng)指標(biāo)的最大值，記為{Xsk}=max（1≤i≤n）{Xik}，最差參考序列的各指標(biāo)對應(yīng)n個不同產(chǎn)地樣品對應(yīng)指標(biāo)的最小值，記為{Xtk}=min（1≤i≤n）{Xik}。

2.5.2 原始數(shù)據(jù)規(guī)格化處理 以平均值變換最為常用，見公式（1）。

2.5.4 關(guān)聯(lián)度的計算 相對于最優(yōu)參考序列、最差參考序列的關(guān)聯(lián)度分別按公式（4）、（5）計算得到。

2.5.5 相對關(guān)聯(lián)度的計算 最佳評價單元是評價單元與最優(yōu)參考序列關(guān)聯(lián)程度最大，同時與最差參考序列的關(guān)聯(lián)程度最小值，所以將評價單元序列同時相對于最優(yōu)參考序列和最差參考序列的相對ri關(guān)聯(lián)度定義為公式（6）。

2.5.6 樣品數(shù)據(jù)集的建立 測定不同產(chǎn)地牡丹皮藥材中9種成分的含量，建立牡丹皮藥材質(zhì)量灰色模式識別數(shù)據(jù)集，見表4。

2.5.7 數(shù)據(jù)規(guī)格化處理 將原始測定數(shù)值按公式（1）進行數(shù)據(jù)規(guī)格化處理，結(jié)果見表5。

2.5.8 關(guān)聯(lián)系數(shù)的計算 分別按照公式（2）、（3）計算各評價單元相對于最優(yōu)、最差參考序列的關(guān)聯(lián)系數(shù)ζik（s）、ζik（t），結(jié)果見表6、表7。

2.5.9 關(guān)聯(lián)度與相對關(guān)聯(lián)度的計算 分別按公式（4）、（5）計算各評價單元相對于最優(yōu)、最差參考序列的關(guān)聯(lián)度ri（s）、ri（t），按公式（6）計算出17個不同產(chǎn)地牡丹皮藥材樣品的相對關(guān)聯(lián)度ri，并進行質(zhì)量排序，結(jié)果見表8。

2.5.10 藥材質(zhì)量評價 相對關(guān)聯(lián)度ri越大，表明該產(chǎn)地牡丹皮藥材的質(zhì)量評價越高，ri排序結(jié)果可作為藥材質(zhì)量高低評價的依據(jù)。由表8可見，17個產(chǎn)地牡丹皮藥材樣品的相對關(guān)聯(lián)度在0.314～0.597之間，其中6個產(chǎn)地樣品的相對關(guān)聯(lián)度>0.50，5個產(chǎn)地樣品的相對關(guān)聯(lián)度在0.45～0.50之間，6個產(chǎn)地樣品的相對關(guān)聯(lián)度<0.45，表明不同產(chǎn)地牡丹皮藥材質(zhì)量存在差異。山西運城絳縣（S07）、山西運城聞喜（S06）、安徽銅陵（S01）、安徽南陵（S09）和陜西商洛（S15）的相對關(guān)聯(lián)度分別為0.597、0.583、0.546、0.537、0.532，排名前5位，藥材質(zhì)量評價優(yōu)。湖北恩施利川（S10）、重慶墊江（S17）、亳州十九里（S03）、亳州五馬（S04）、亳州華佗（S02）和山東菏澤（S05）的相對關(guān)聯(lián)度分別為0.414、0.401、0.392、0.381、0.376、0.314，排名后6位，藥材質(zhì)量評價較差。

3 討論

3.1 流動相系統(tǒng)的選擇 在建立牡丹皮藥材多成分含量測定的過程中，曾對流動相洗脫系統(tǒng)進行優(yōu)選，選擇甲醇-0.1%磷酸水溶液、甲醇-0.1%甲酸水溶液、乙腈-0.1%甲酸水溶液、乙腈-0.1%醋酸水溶液等多種流動相洗脫系統(tǒng)進行比較，結(jié)果顯示，甲醇-0.1%磷酸水溶液流動相系統(tǒng)基線平穩(wěn)，干擾較小，各指標(biāo)成分峰分離度良好，故選擇甲醇-0.1%磷酸水溶液為牡丹皮的流動相系統(tǒng)。

3.2 提取方式的選擇 本試驗考察超聲提取、回流提取的提取方式，結(jié)果顯示，回流提取的芍藥苷和丹皮酚峰面積均高于超聲提取，且重復(fù)性結(jié)果優(yōu)于超聲提取。比較回流提取0.5h、1h、1.5h，回流1h和1.5h含量相近，而芍藥苷峰面積略高于0.5h，綜合考慮，確定牡丹皮供試液的提取方式為回流提取1h。

3.3 粉末粒徑的選擇 本試驗考察2號篩、3號篩、4號篩不同篩目的粉末粒徑，結(jié)果顯示，2號篩樣品丹皮酚的重復(fù)性較3號篩、4號篩的差。3號篩與4號篩樣品的重復(fù)性總體情況良好，兩者含量相近，其芍藥苷峰信號均高2號篩粉末樣品，綜合考慮，確定選擇牡丹皮粉末粒徑為3號篩。

3.4 定容溶劑的選擇 本試驗在稀釋定容過程中考察0%～100%的甲醇溶液，結(jié)果顯示，0%～40%的甲醇溶液稀釋的色譜圖中色譜峰數(shù)量較少，而60%～100%的甲醇溶液稀釋的部分色譜峰的峰型較差，分離度低，部分色譜峰出現(xiàn)裂峰，而50%的甲醇溶液稀釋較為良好，色譜峰數(shù)量較多，信息量較大。綜合考慮，確定選用50%的甲醇溶液進行稀釋定容。

3.5 檢測波長的選擇 本試驗考察230nm、258nm、265nm、274nm、292nm等檢測波長，結(jié)果顯示，檢測波長為258～292nm，色譜峰數(shù)量減少，包括芍藥苷在內(nèi)的一些特征信息未表達(dá)出來;檢測波長為335nm時，基本上只有丹皮酚一個色譜峰，而檢測波長為230nm時，色譜圖的峰數(shù)目、信息量較大，綜合考慮，確定選擇230 nm為牡丹皮的檢測波長。

4 結(jié)論

本實驗結(jié)果表明，不同產(chǎn)地牡丹皮藥材質(zhì)量存在差異。山西運城產(chǎn)地牡丹皮藥材中兒茶素、氧化芍藥苷、芍藥苷和浸出物含量明顯高于其他產(chǎn)地，說明其藥材質(zhì)量評價較高;安徽銅陵、南陵所產(chǎn)的牡丹皮藥材質(zhì)量排序第3、第4名，說明其藥材品質(zhì)較優(yōu)，這也與安徽銅陵和南陵為牡丹皮的道地產(chǎn)區(qū)相吻合;湖北恩施、重慶墊江、亳州、山東菏澤產(chǎn)地牡丹皮樣品的相對關(guān)聯(lián)度排名靠后，說明其藥材質(zhì)量較差。

牡丹皮化學(xué)成分眾多，成分的多樣性致使其質(zhì)量很難控制，單一的指標(biāo)成分不能準(zhǔn)確反映牡丹皮的藥效和質(zhì)量，使其在分析評價藥材的質(zhì)量時損失信息量太多。本研究以相對關(guān)聯(lián)度為測度，結(jié)合灰色關(guān)聯(lián)分析法，建立一種牡丹皮藥材質(zhì)量綜合評價的模型，客觀整體地反映牡丹皮藥材的內(nèi)在質(zhì)量，為牡丹皮藥材的綜合質(zhì)量評價提供了新思路，也為牡丹皮資源的開發(fā)利用提供參考依據(jù)。

參考文獻(xiàn)

[1]國家藥典委員會.中華人民共和國藥典（一部）[S].北京：中國醫(yī)藥科技出版社，2015：172.

[2]張艷，范俊安.中藥材牡丹皮研究概況III 丹皮化學(xué)成分研究概況[J].重慶中草藥研究，2008（2）：24-31.

[3]DING Y，WU E Q，CHEN J B，et al.Quality evaluation of moutan cortex radicis using multiple component analysis by high performance liquid chromatography [J].Bull.Korean Chem.Soc.，2009，30（10）：2240-2244.

[4]LIN H C，DING H Y，WU T S，et al.Monoterpene glycosides from Paeonia Suffruticosa[J].Phy-tochemistry，1996，41（1）：237-242.

[5]許舜軍，楊柳，張勉，等.牡丹皮化學(xué)成分的液相色譜-飛行時間串聯(lián)質(zhì)譜分析[J].藥學(xué)學(xué)報，2006，41（9）：852-856.

[6]王子堯，陳彥，孫文秀，等.超濾法分離丹皮多糖的研究[J].膜科學(xué)與技術(shù)，2009，29（4）：98-100.

[7]鄧斌，鄧勝軍，章愛華.原子吸收光譜法測定牡丹皮中的微量金屬元素鐵、鋅、銅和錳[J].光譜實驗室，2008，25（4）：630-632.

[8]鄧聚龍.灰理論基礎(chǔ)[M].武漢：華中理工大學(xué)出版社，2002：5-10.

[9]劉思峰，蔡華，楊英杰，等.灰色關(guān)聯(lián)分析模型研究進展[J].系統(tǒng)工程理論與實踐，2013，33（8）：2041-2046.

[10]劉思峰，楊英杰，吳利豐.灰色系統(tǒng)理論及其應(yīng)用[M].7版.北京：科學(xué)出版社，2014：63-83.

[11]劉思峰，楊英杰.灰色系統(tǒng)研究進展（2004-2014）[J].南京航空航天大學(xué)學(xué)報，2015，47（1）：1-6.

[12]李成義，強正澤，王燕，等.基于灰色關(guān)聯(lián)分析方法評價甘肅產(chǎn)紅芪質(zhì)量[J].中國實驗方劑學(xué)雜志，2016，22（12）：60-65.

[13]張洪坤，路麗，黃玉瑤，等.牡丹皮不同采收期質(zhì)量綜合評價研究[J].世界科學(xué)技術(shù)-中醫(yī)藥現(xiàn)代化，2019，21（2）：240-247.

[14]張洪坤，王其豐，郭長達(dá)，等.不同加工方法牡丹皮中7種指標(biāo)性成分的含量測定及質(zhì)量評價[J].中國藥房，2018，29（22）：3063-3068.