陳朝喜,賀冬梅,文 娟,湯 承
(西南民族大學(xué)生命科學(xué)與技術(shù)學(xué)院,四川成都610041)
藏獸醫(yī)學(xué)是祖國(guó)獸醫(yī)學(xué)的重要組成部分,也是藏族勞動(dòng)人民同家畜疾病做斗爭(zhēng)的經(jīng)驗(yàn)總結(jié),是藏族勞動(dòng)人民在長(zhǎng)期的醫(yī)療實(shí)踐中,吸取中醫(yī)理論和外來醫(yī)學(xué)的精華,再結(jié)合本地區(qū)的特色建立起來的具有獨(dú)特理論體系的民族獸醫(yī)學(xué),對(duì)于保證廣大牧區(qū)畜牧業(yè)的發(fā)展有著重要的作用和貢獻(xiàn)。然而,長(zhǎng)期以來,同時(shí)由于歷史條件的限制,對(duì)藏獸醫(yī)驗(yàn)方的收集和整理只限于部分畜牧業(yè)發(fā)達(dá)的地區(qū),使得藏獸醫(yī)驗(yàn)方資源的發(fā)展和傳承受到一定的限制。因此,本研究利用中醫(yī)傳承輔助系統(tǒng),對(duì)收集到的藏獸醫(yī)驗(yàn)方和文獻(xiàn)資料記載的方劑進(jìn)行數(shù)據(jù)挖掘,對(duì)收集到的630個(gè)藏獸醫(yī)驗(yàn)方和26個(gè)含有纖毛婆婆納的驗(yàn)方進(jìn)行組方規(guī)律分析,以期將寶貴的藏獸醫(yī)驗(yàn)方得到系統(tǒng)的整理,同時(shí)為后期系統(tǒng)研究收集到藏獸醫(yī)驗(yàn)方的現(xiàn)代藥理學(xué)機(jī)理奠定基礎(chǔ)。
1.1.1 處方來源 630個(gè)藏獸醫(yī)方劑主要來源于《藏獸醫(yī)經(jīng)驗(yàn)選編》[1]、《常用藏藥志》[2]、《現(xiàn)代中獸醫(yī)大全》[3]、《藏獸醫(yī)驗(yàn)方選》[4]、《高原中草藥治療手冊(cè)》[5]及中國(guó)知網(wǎng)(CNKI)檢索“藏獸醫(yī)”和“藏獸藥”關(guān)鍵詞篩選到的相關(guān)文獻(xiàn)(1975年-2013年)[6-13]。
1.1.2 分析軟件 “中醫(yī)傳承輔助系統(tǒng)(V2.0)”軟件,由中國(guó)中醫(yī)科學(xué)院中藥研究所提供。
1.2.1 處方的輸入和核對(duì) 將1.1文獻(xiàn)資料中篩選的藏獸醫(yī)驗(yàn)方方劑錄入中醫(yī)傳承輔助系統(tǒng)。為確保數(shù)據(jù)源和篩選方劑的數(shù)據(jù)挖掘結(jié)果的準(zhǔn)確性,避免錄入過程可能出現(xiàn)遺漏或重復(fù)輸入等錯(cuò)誤,在完成篩選方劑的錄入后,再由錄入人員以外人員負(fù)責(zé)數(shù)據(jù)源的審核和修正。
1.2.2 數(shù)據(jù)分析 用中醫(yī)傳承輔助系統(tǒng)“數(shù)據(jù)報(bào)表系統(tǒng)”模塊的“方劑分析”功能,對(duì)630個(gè)藏獸醫(yī)方劑進(jìn)行“常用藥物頻次分析”和“方劑藥性分析”。同時(shí),基于關(guān)聯(lián)規(guī)則和熵聚變算法對(duì)含有藏藥關(guān)鍵詞“纖毛婆婆納”的26個(gè)方劑進(jìn)行組方規(guī)律分析。
利用中醫(yī)傳承輔助系統(tǒng)軟件對(duì)收集到的630個(gè)藏獸醫(yī)方劑進(jìn)行藥物頻次、歸類統(tǒng)計(jì)、四氣統(tǒng)計(jì)、五味統(tǒng)計(jì)和歸經(jīng)統(tǒng)計(jì),分析其組方規(guī)律。
2.1.1 藥物頻次 從參考文獻(xiàn)中共檢索到630個(gè)藏獸醫(yī)驗(yàn)方(其中四川省409 個(gè),西藏自治區(qū)119個(gè),甘肅省59個(gè),云南省3個(gè),青海省15個(gè),未提供驗(yàn)方來源單位的驗(yàn)方25 個(gè)),共有1 143 味藥物。對(duì)630個(gè)處方中的藥物頻次進(jìn)行統(tǒng)計(jì),使用頻次高于20的藥物共有32味,其中植物藥物25味,動(dòng)物藥物2味,礦物藥物5味(表1)。
2.1.2 藥物歸類 對(duì)使用頻次高于100的科屬進(jìn)行歸類分析,有8個(gè)種屬的藥物使用頻次較高,菊科植物的使用頻次最高,其次為毛茛科、豆科、傘形科、唇形科、姜科、使君子科和玄參科(圖1)。值得一提的是盡管其他科屬藥物使用頻次在100次以下,但其使用頻次遠(yuǎn)遠(yuǎn)高于傳統(tǒng)中獸醫(yī)臨床用藥頻次(如鹿科,數(shù)據(jù)未列出),這可能與藏區(qū)藥物資源優(yōu)勢(shì)和藏獸藥的使用習(xí)慣有一定的關(guān)系。
表1 630個(gè)方劑中頻次20以上的藥物Table 1 Drugs of frequency over 20times in 630prescriptions
圖1 630個(gè)方劑藥物科統(tǒng)計(jì)分析Fig.1 Statistical analysis of family in 630prescriptions
2.1.3 四氣統(tǒng)計(jì)分析 630個(gè)方劑中藥物的四氣統(tǒng)計(jì)分析結(jié)果見圖2,溫性藥物的使用頻次最高,共使用1 371次;熱性藥物使用頻次最低,共使用217次。
2.1.4 五味統(tǒng)計(jì)分析 630個(gè)方劑中藥物的五味統(tǒng)計(jì)分析結(jié)果見圖3,苦味藥物的使用頻次最高,共使用2 110次;淡味藥物使用頻次最低,共使用33次。其他性味使用頻次順序?yàn)樾粒靖剩緷舅幔鞠獭?/p>
2.1.5 歸經(jīng)統(tǒng)計(jì)分析 630個(gè)方劑中藥物的歸經(jīng)統(tǒng)計(jì)分析結(jié)果見圖4,歸經(jīng)頻次最高的是肝臟,共1 406次;歸經(jīng)頻次最低的是三焦,共74次。其他歸經(jīng)的統(tǒng)計(jì)結(jié)果除了歸經(jīng)為膽、膀胱、小腸和心包的藥物使用頻次低于260次外,其他歸經(jīng)的藥物使用頻次都在600次以上。
圖2 630個(gè)方劑藥物四氣統(tǒng)計(jì)分析Fig.2 Statistical analysis of four properties in 630prescriptions
圖3 630個(gè)方劑藥物五味統(tǒng)計(jì)分析Fig.3 Statistical analysis of five tastes in 630prescriptions
圖4 630個(gè)方劑藥物歸經(jīng)統(tǒng)計(jì)分析Fig.4 Statistical analysis of channel-tropism in 630prescriptions
從630個(gè)藏獸醫(yī)方劑中共篩選到26個(gè)含有“纖毛婆婆納”的方劑,共166味藥物。利用中醫(yī)傳承輔助系統(tǒng)“方劑分析系統(tǒng)”模塊,對(duì)其組方規(guī)律進(jìn)行分析,從而篩選出新的組方。
2.2.1 基于關(guān)聯(lián)規(guī)則分析的組方規(guī)律研究 ①用藥模式分析。在設(shè)置支持度=20%和置信度=90%的條件下對(duì)26個(gè)含有“纖毛婆婆納”的藏獸醫(yī)方劑進(jìn)行組方規(guī)律分析,26個(gè)方劑中使用頻次在5以上的組合見表2。使用頻次在5以上的用藥組合共35個(gè),其中排在前5位的組合分別是“五靈脂,纖毛婆婆納”、“訶子,纖毛婆婆納”、“廣木香,纖毛婆婆納”、“波棱瓜子,纖毛婆婆納”和“五脈綠絨蒿,纖毛婆婆納”。②規(guī)則分析。26個(gè)方劑藥用藥規(guī)則見表3,共有12條規(guī)則。從表中看出,在2味藥模式中,牛黃、波棱瓜子、五靈脂、五脈綠絨蒿、訶子和廣木香與纖毛婆婆納合用的幾率較高;在3味藥模式中,“五靈脂,五脈綠絨蒿”和“紅花,五靈脂”組合與纖毛婆婆納合用的幾率較高。③網(wǎng)絡(luò)展示:對(duì)26個(gè)方劑的12條規(guī)則進(jìn)行網(wǎng)絡(luò)顯示,直觀展示藥物間相互關(guān)系(圖5)。
續(xù)表2
表3 26個(gè)方劑中藥物組合的用藥規(guī)則(置信度>0.9)Table 3 Rules of herbal combinations in 26prescriptions(confidence level>0.9)
2.2.2 基于熵聚變的方劑組方規(guī)律研究 ①基于改進(jìn)的互信息法的藥物間關(guān)聯(lián)度分析:根據(jù)方劑數(shù)量的多少,同時(shí)結(jié)合經(jīng)驗(yàn)判斷,在設(shè)置相關(guān)度為8,懲罰度為2的情況下進(jìn)行聚類分析,得到26個(gè)方劑中兩兩藥物間的相關(guān)度,共得到149條藥物關(guān)聯(lián)信息,將關(guān)聯(lián)系數(shù)0.08以上的藥對(duì)列表(表4)。②基于復(fù)雜系統(tǒng)熵聚類的藥物核心組合分析:以改進(jìn)的互信息法的藥物間關(guān)聯(lián)度分析結(jié)果為基礎(chǔ),按照相關(guān)度為8和懲罰度為2的約束規(guī)則,基于復(fù)雜系統(tǒng)熵聚類,得到4個(gè)藥物核心組合,分別是“翼首草,訶子,廣木香”、“余甘子,白檀香,高山龍膽”、“波棱瓜子,牛黃,藏馬兜鈴”和“余甘子,白檀香,甘青青蘭”。③基于無(wú)監(jiān)督熵層次聚類的新處方分析:在提取核心組合基礎(chǔ)上,運(yùn)用無(wú)監(jiān)督熵層次聚類算法,得到2個(gè)新處方,分別是“翼首草,訶子,廣木香,波棱瓜子,牛黃,藏馬兜鈴”和“余甘子,白檀香,高山龍膽,甘青青蘭”,候選新方的網(wǎng)絡(luò)展示見圖6。
藏獸醫(yī)藥學(xué)是我國(guó)古老的科學(xué)文化遺產(chǎn)之一,也是我國(guó)傳統(tǒng)獸醫(yī)藥學(xué)的重要分支。藏獸醫(yī)藥學(xué)隨著藏區(qū)科學(xué)文化和畜牧業(yè)生產(chǎn)技術(shù)的進(jìn)步而產(chǎn)生和發(fā)展。解放后,藏區(qū)獸醫(yī)教育迅速發(fā)展,改變了解放前由于歷史條件的限制使得很多古方和驗(yàn)方僅能采用口傳心授方式傳承下來的局面,出版了大量的獸醫(yī)學(xué)著作并培養(yǎng)了大批獸醫(yī)人員,使得藏獸醫(yī)資源得到記錄和傳承。
中醫(yī)傳承輔助系統(tǒng)軟件能夠?qū)︶t(yī)案進(jìn)行數(shù)據(jù)挖掘和發(fā)現(xiàn)用藥規(guī)律的總結(jié),為中藥的二次開發(fā)及古方和驗(yàn)方的開發(fā)研究尤其是中藥新藥的創(chuàng)制也具有重要的價(jià)值[13]。利用中醫(yī)傳承系統(tǒng)軟件網(wǎng)絡(luò)資源對(duì)藏獸醫(yī)經(jīng)驗(yàn)方進(jìn)行統(tǒng)計(jì)和處理,同時(shí)對(duì)相關(guān)數(shù)據(jù)進(jìn)行分析和挖掘,總結(jié)藏獸藥的用藥規(guī)律,對(duì)藏獸藥的開發(fā)和傳承都有一定的指導(dǎo)意義。同時(shí),能夠?yàn)椴孬F藥的現(xiàn)代藥理學(xué)研究提供新的思路和方法,為篩選更好的新型方劑和開拓新的應(yīng)用領(lǐng)域奠定基礎(chǔ)。
圖5 支持度為5,置信度為0.9的藥物網(wǎng)絡(luò)展示Fig.5 Network display of new prescriptions on condition that the relevancy is 5and confidence level is 0.9
表4 基于改進(jìn)的互信息法的藥物間關(guān)聯(lián)度分析Table 4 Correlation coefficient between herbs based on revised mutual information method
本研究對(duì)檢索到的630個(gè)藏獸醫(yī)驗(yàn)方組方藥物進(jìn)行四氣五味和歸經(jīng)等分析,同時(shí)以纖毛婆婆納為例進(jìn)行新藥方篩選,在結(jié)果處理過程中發(fā)現(xiàn)收集到的驗(yàn)方組方有以下特征:
3.2.1 藥物組成數(shù)量 收集到的630個(gè)藏獸醫(yī)驗(yàn)方中,藥物組成數(shù)量差別很大:最少的只使用1味藥物,如沙棘膏用沙棘果實(shí)制得;最多的方劑使用25味藥物,如方劑“五味子25味”和“獐芽菜25味”,分別使用五味子和獐芽菜等25味藥物組方。
3.2.2 藥物使用劑量 在藥物劑量方面,很多藥物包括毒性較大的藥物使用量比較大,如瑞香狼毒在治療牛羊疥癬時(shí)使用劑量為15g~30g[1],考慮到個(gè)體差異因素很可能會(huì)引起動(dòng)物中毒事件,在使用中必須引起重視。
圖6 相關(guān)度為8,懲罰度為2的新處方藥物網(wǎng)絡(luò)展示Fig.6 Network display of new prescriptions drugs on condition that the relevancy is 8and penalty index is 2
3.2.3 藥物使用種類 藏區(qū)中草藥資源非常豐富,在藏獸藥的配方和使用中,其組方藥材的種類也較多使用了雪蓮花、紅景天和藏木香等大批高原地區(qū)稀有的藥材資源。此外,除了植物藥物外,也較多使用了動(dòng)物藥物和礦物藥物(如五靈脂和硇砂等)。王賽等[14]也通過對(duì)藏區(qū)特色植物綠蘿花進(jìn)行降血糖的研究,結(jié)果表明其降血糖效果顯著,并且能顯著提升造模小鼠血清胰島素的水平。馮慎遠(yuǎn)等[15]也利用藏藥五根散對(duì)缺血性神經(jīng)損傷大鼠模型進(jìn)行研究,結(jié)果表明五根散能夠降低缺血性大鼠新皮質(zhì)和海馬中內(nèi)皮素的含量從而保護(hù)腦組織。
本研究中收集的630個(gè)方劑,雖然都來自書面文件,但較少來自專題研究資料,很多驗(yàn)方未能在實(shí)踐中通過大量的動(dòng)物實(shí)驗(yàn)驗(yàn)證,而且部分藥物的學(xué)名、科屬和性味功能等信息還需要進(jìn)一步考證。因此,在后續(xù)的研究中,我們有必要對(duì)收集的驗(yàn)方方劑進(jìn)行進(jìn)一步考究,同時(shí)根據(jù)具體情況和處方來源,對(duì)組方藥材的配伍和用量等信息進(jìn)行核對(duì),并通過現(xiàn)代藥理學(xué)方法進(jìn)一步研究其藥效學(xué)和藥動(dòng)學(xué),以使其在藏區(qū)養(yǎng)殖業(yè)中更好地發(fā)揮作用,達(dá)到真正的防治疾病的目的。
本研究利用中醫(yī)傳承系統(tǒng)軟件篩選出2個(gè)新處方,分別與尕爾戈報(bào)道的“感冒散”和“龍膽十五味”組方有一定相似性[10],為研究篩選的新方劑的藥理學(xué)作用提供了參考依據(jù)和研究思路。然而,新篩選藥方能否和報(bào)道的驗(yàn)方一樣發(fā)揮藥理學(xué)作用,還需要通過與相似度高的驗(yàn)方進(jìn)行藥理學(xué)比較分析,同時(shí)結(jié)合理論分析和臨床驗(yàn)證才能做出全面的評(píng)價(jià),從而為藏獸醫(yī)的豐富和發(fā)展提供新的材料和思路。
致謝:軟件的安裝和數(shù)據(jù)處理過程得到中國(guó)中醫(yī)科學(xué)院中藥研究所研究生申丹的幫助,在此表示感謝!
[1] 《藏獸醫(yī)經(jīng)驗(yàn)選編》編寫組.藏獸醫(yī)經(jīng)驗(yàn)選編[M].北京:農(nóng)業(yè)出版社,1979:34-174.
[2] 田淑琴,鄧孝廷,張 琦.常用藏藥志[M].四川成都:四川科學(xué)技術(shù)出版社,1997:225-227.
[3] 于 船,陳子斌.現(xiàn)代中獸醫(yī)大全[M].廣西南寧:廣西科學(xué)技術(shù)出版社,2000:695-725.
[4] 若爾蓋縣紅星鄉(xiāng)獸防站.藏獸醫(yī)驗(yàn)方選[M].四川成都:四川民族出版社,2011:1-180.
[5] 姚海潮,色 珠,拉巴次旦,等.一復(fù)方藏獸藥體外抗菌活性及藥效學(xué)實(shí)驗(yàn)研究[J].廣東畜牧獸醫(yī)科技,2010,35(6):33-35.
[6] 何 文,安海宏,張富亨,等.藏獸藥治療羔羊拉?。跩].畜牧獸醫(yī)雜志,2009,28(6):118.
[7] 魯國(guó)義.藏獸醫(yī)對(duì)家畜肺炎和結(jié)癥的治療驗(yàn)方[J].中獸醫(yī)學(xué)雜志,1988(1):40-41.
[8] 魯國(guó)義.藏獸醫(yī)治療家畜疥癬病驗(yàn)方[J].中國(guó)獸醫(yī)科技,1986(12):56.
[9] 陸宏開.藏獸醫(yī)治療幾種牦牛病的經(jīng)驗(yàn)[J].西南民族大學(xué)學(xué)報(bào):自然科學(xué)版,1984(4):66-67.
[10] 尕爾戈.藏獸醫(yī)方藥選輯(一)[J].西南民族大學(xué)學(xué)報(bào):自然科學(xué)版,1993,19(3):229-237.
[11] 索南木.藏獸醫(yī)驗(yàn)方、單方和偏方選[J].西南民族大學(xué)學(xué)報(bào):自然科學(xué)版,1983,19(3):59-61.
[12] 江士辰.藏獸醫(yī)驗(yàn)方選(五)——色達(dá)縣大則公社藏獸醫(yī)驗(yàn)方(續(xù))[J].西南民族學(xué)院學(xué)報(bào):畜牧獸醫(yī)版,1980(1):36-41.
[13] Yang H J,Chen J X,Tang S H,et al.New drug R&D of traditional Chinese medicine-role of data mining approaches[J].J Biol Syst,2009,17(3):329-347.
[14] 王 賽,劉 群,韓金潭,等.藏藥綠蘿花對(duì)高血糖小鼠模型的治療試驗(yàn)[J].動(dòng)物醫(yī)學(xué)進(jìn)展,2015,39(4):67-72.
[15] 馮慎遠(yuǎn),陳廣迪.藏藥五根散對(duì)缺血性腦損傷大鼠新皮質(zhì)與海馬中內(nèi)皮素含量變化的影響[J].動(dòng)物醫(yī)學(xué)進(jìn)展,2003,24(2):83-84.