基于數(shù)據(jù)挖掘和網(wǎng)絡(luò)藥理學(xué)分析《中國藥典》治療心悸藥物配伍規(guī)律及其作用機制

徐寧陽，王 群，劉 悅，張 哲，楊關(guān)林*

1遼寧中醫(yī)藥大學(xué)，沈陽 110847；2遼寧中醫(yī)藥大學(xué)附屬醫(yī)院，沈陽 110032

心悸是中醫(yī)病證名，《傷寒雜病論》稱之為 “心動悸”“心下悸”“心中悸”“驚悸”?；颊咦愿袨榧聞硬话?、驚惕心慌甚至不能自主。情志不遂或者過度勞累時易發(fā)作，多伴有氣短、耳鳴、眩暈、胸悶、憋氣等癥狀。根據(jù)癥狀心悸在西醫(yī)中屬于心律失常范疇，是我國最為常見的心血管疾病類型，具有高發(fā)生率、死亡率，疾病臨床類型較多，常見有心動過速、心動過緩、心房顫動、傳導(dǎo)異常等[1,2]。如果不加以控制，不僅會對心肌細胞造成不可逆的損傷，還會引發(fā)心肌梗死等嚴重并發(fā)癥，對患者的生活有著難以改善的影響，降低患者的生活質(zhì)量的同時也會危急患者的生命。目前西醫(yī)治療多以抗心律失常藥物治療，但大量實踐表明仍有患者對西醫(yī)治療會產(chǎn)生復(fù)發(fā)或者藥物不能耐受等不良后果[3]。近幾年由于中醫(yī)藥在臨床治療中廣泛應(yīng)用，故有學(xué)者提出以中藥治療心律失常，心悸最早由張仲景提出，因此通過古代的中藥研究，如何能更好的醫(yī)治心悸，其具體分子靶點是怎樣的，本研究將從2015年《中國藥典》(一部)(以下簡稱《中國藥典》)入手，挖掘《中國藥典》中用于治療心悸中成藥的組方規(guī)律，分析出高頻藥物。

《中國藥典》是將藥品從研制到生產(chǎn)、從經(jīng)營到使用管理等相關(guān)從業(yè)者必須遵守的法定依據(jù)，是我國保證藥品質(zhì)量的法典。在前期文獻研究中發(fā)現(xiàn)有關(guān)學(xué)者已對藥典所收載的中藥飲片或毒性藥材展開歸類及配伍規(guī)律的研究[4]；對收載的成方制劑開展主治病證的中成藥組方規(guī)律分析[5,6]；對《中國藥典》收載的成方制劑和單味制劑進行系統(tǒng)數(shù)據(jù)挖掘，但均未發(fā)現(xiàn)對治療某一相關(guān)疾病尤其心血管疾病的中成藥進行數(shù)據(jù)挖掘及分析規(guī)律，故本研究以2015年版《中國藥典》中的治療心悸的中成藥的方劑名稱、中藥成分、中藥用量、功能與主治為研究對象，采用由中國中醫(yī)科學(xué)院中藥研究所研發(fā)的V2.5版中醫(yī)傳承輔助平臺分析用藥規(guī)律[7,8]，通過其主要功能關(guān)聯(lián)規(guī)則Apriori算法、復(fù)雜系統(tǒng)熵聚類等無監(jiān)督數(shù)據(jù)挖掘方法，分析中藥的使用頻次及中藥之間的關(guān)聯(lián)規(guī)則、處方規(guī)律以及新方的預(yù)測。

由上述軟件分析得出的用藥規(guī)律后，得出治療心悸的高頻藥物，通過由Hopkins等人提出的探索中藥與疾病直接關(guān)系的科研方法—網(wǎng)絡(luò)藥理學(xué)，系統(tǒng)的分析和挖掘，構(gòu)建“中藥-有效成分-疾病靶點-疾病”之間的關(guān)系，能夠為中藥機制分析及藥物研發(fā)提供新思路、新方法[9,11]。由于中藥具備多成分、多靶點的因素，所以網(wǎng)絡(luò)藥理學(xué)正是能夠更好的闡述中藥作用于心悸的作用機制。

1 資料和方法

1.1 治療心悸的中成藥收集

檢索2015年《中國藥典》(一部)中成方制劑與單位制劑部分的藥品說明書，在功能與主治中選取關(guān)鍵詞為：“心悸”“心律失?！钡茸謽佑糜谥委熜募碌闹谐伤?。在選取后的中成藥中排出劑型不同但藥方一致的方劑，如：歸脾丸、歸脾合劑、歸脾顆粒計為1個中成藥。

1.2 數(shù)據(jù)的錄入、核對與處理

在《中國藥典》中，將符合“1.1”條件的中成藥方中的方劑名稱、中藥成分、中藥用量、功能與主治錄入Microsoft Excel for Mac；對所納入的中藥的名稱進行統(tǒng)一規(guī)范，如：沙棘鮮漿規(guī)范為沙棘、醋龜甲規(guī)范為龜甲。將所涉及的藏藥，蒙藥，苗藥或者中藥為同一藥材但名稱不同的，按照中國藥典及十三五規(guī)劃的中藥學(xué)教材標準統(tǒng)一錄入。規(guī)范后的數(shù)據(jù)錄入中醫(yī)傳承輔助平臺V2.5，采用兩人輸入，第三人校對審核，確保數(shù)據(jù)嚴謹、真實。

1.3 數(shù)據(jù)分析

1.3.1 中醫(yī)藥傳承輔助平臺對用藥規(guī)律的分析與預(yù)測

將“1.2”所錄入的數(shù)據(jù)在中醫(yī)藥傳承輔助平臺V2.5相應(yīng)模塊中進行“頻次統(tǒng)計”，“組方規(guī)律”，“新方分析”的用藥規(guī)律分析及預(yù)測。

1.3.2 網(wǎng)絡(luò)藥理學(xué)對藥物的深度分析

將“1.3.1”中所得到的最高頻藥物，在中醫(yī)藥系統(tǒng)藥理學(xué)平臺(The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform，TCMSP)[12]先后進行中藥有效成分篩選，同時在GeneCards數(shù)據(jù)庫以心律失常(arrhythmia)為關(guān)鍵詞，檢索心悸的相關(guān)疾病靶點，再使用Draw Venny Diagram在線程序構(gòu)建中藥-疾病交集靶點，在Cytoscape 3.7.1中構(gòu)建藥物-疾病-靶基因網(wǎng)絡(luò)圖，并將基因提交至String網(wǎng)站繪制蛋白質(zhì)相互作用關(guān)系圖(PPI)，再利用R3.6.1軟件進行GO富集分析和KEGG通路富集分析，來預(yù)測《中國藥典》中治療心悸的最高頻藥物的可能作用機制。

2 結(jié)果

2.1 藥頻統(tǒng)計--治療心悸的中成藥中各中藥的使用頻率

在《中國藥典》中共選取到93個中成藥方，其中14個為劑型不同但存在藥方一致，故剩余79個有效中成藥方。根據(jù)《中國藥典》，對79個有效中成藥中的176味中藥進行“頻次統(tǒng)計”，得到應(yīng)用頻次不少于8的中藥22個，其中用藥頻次最高的五位藥分別是：丹參31次；人參23次；麥冬18次；川芎和五味子各17次。詳見表1。

表1 《中國藥典》中治療心悸的中成藥的用藥頻次(≥8次)

2.2 關(guān)聯(lián)分析--治療心悸處方的組方規(guī)律分析

2.2.1 《中國藥典》中治療心悸常用藥對分析

根據(jù)Apriori算法，在設(shè)置支持度為7的情況下，得到25個高頻藥對。其中，出現(xiàn)頻度最高的四對配伍分別是：丹參-人參13次；川芎-丹參12次；麥冬-五味子和丹參-五味子各10次(詳見表2)。

表2 《中國藥典》中治療心悸處方中高頻藥對

續(xù)表2(Continued Tab.2)

2.2.2 《中國藥典》中治療心悸藥物應(yīng)用規(guī)則分析

同樣在Apriori算法中，設(shè)置置信度0.6，得到9種應(yīng)用規(guī)則，其中石菖蒲->丹參；紅花->川芎；當歸、遠志->酸棗仁的置信度最高，為0.875。規(guī)則詳見表3和圖1。

表3 《中國藥典》中治療心悸藥物規(guī)則分析

圖1 《中國藥典》中治療心悸藥物關(guān)聯(lián)規(guī)則可視化Fig.1 Visualization of association rules of drugs for treatment of arrhythmia in the Chinese Pharmacopoeia

2.3 聚類分析--基于熵聚類方法的新處方預(yù)測

在應(yīng)用“新方分析”模塊中，設(shè)置相關(guān)度8，懲罰度2，通過無監(jiān)督的熵層次聚類分析，得到核心組合16個，見表4。并進一步分析得到用于治療心悸新方組合10個(見表5，圖2)。

圖2 《中國藥典》中治療心悸藥物新處方網(wǎng)絡(luò)展示圖Fig.2 Network display diagram of new prescriptions for arrhythmia in the Chinese Pharmacopoeia

圖3 丹參活性化合物靶點與心悸靶點的交集基因Fig.3 Intersection genes between the target of Salviae Miltiorrhizae Radix et Rhizoma active compound and the target of arrhythmia

圖4 丹參與心悸蛋白質(zhì)相互作用核心網(wǎng)絡(luò)圖Fig.4 The core network diagram of Salviae Miltiorrhizae Radix et Rhizoma involvement in arrhythmia protein interaction

表4 《中國藥典》中治療心悸藥物聚類核心組合

表5 《中國藥典》中治療心悸藥物新處方預(yù)測結(jié)果

2.4 網(wǎng)絡(luò)藥理學(xué)分析--高頻藥物丹參治療心悸的潛在作用機制分析

2.4.1 丹參的靶標預(yù)測結(jié)果

根據(jù)“2.1～2.3”所得到的最關(guān)鍵性藥物丹參進行網(wǎng)絡(luò)藥理學(xué)分析。在TCMSP中預(yù)測到丹參的65個有效成分(篩選要求OB大于30%，DL大于0.18[13,14])。具體信息見表6。

表6 丹參的有效成分表

2.4.2 丹參治療心悸的蛋白、基因分析

將TCMSP中得到的藥物靶點和在GeneCards數(shù)據(jù)庫中得到的疾病靶點進行交集基因預(yù)測，其中有79個交集靶點(見圖3)。在String網(wǎng)站中，生成PPI互作網(wǎng)絡(luò)圖(見圖4)。并根據(jù)其蛋白質(zhì)-蛋白質(zhì)作用數(shù)量，生成蛋白統(tǒng)計條圖(見圖5)。最后對65個中藥靶點(藍色)及79個疾病靶點(綠色)構(gòu)建丹參-有效成分-疾病靶點-疾病的網(wǎng)絡(luò)圖，交集線比較集中的部分可能是其關(guān)鍵靶點，整個網(wǎng)絡(luò)中包括138個節(jié)點，642條邊，degree排名前五名的化合物有木犀草素(luteolin)、丹參酮(tanshinone)、4-亞甲丹參新酮(4-methylenemiltirone)、隱丹參酮(cryptotanshinone)、鼠尾草酚酮(salviolone)，排名前五的基因有環(huán)加氧酶(PTGS2)、β2腎上腺素能受體(ADRB2)、人源全長重組蛋白(SCN5A)、阿片受體編碼基因(OPRM1)、毒蕈堿型膽堿受體M1(CHRM1)(見圖6)。

圖5 蛋白統(tǒng)計條圖Fig.5 Protein statistics bar graph

圖6 丹參-有效成分-靶點-疾病網(wǎng)絡(luò)圖Fig.6 Salviae Miltiorrhizae Radix et Rhizoma-active ingredients-targets-disease network diagram

續(xù)表6(Continued Tab.6)

續(xù)表6(Continued Tab.6)

2.4.3 丹參治療心悸的通路可視化富集分析

為了進一步探究丹參對于心悸的作用功能及通路，利用R3.6.1軟件進行GO富集分析和KEGG通路富集分析，其中GO生物學(xué)功能分析得到結(jié)果見圖7。KEGG通路富集分析得到PI3K-AKT信號通路、IL-17信號通路、TNF信號通路等通路(見圖8)。

圖7 GO生物學(xué)功能富集分析圖Fig.7 Analysis diagram of GO biological function enrichment

圖8 KEGG通路富集分析圖Fig.8 KEGG pathway enrichment analysis diagram

3 討論

心律失常作為比較常見的心系疾病，是臨床心血管疾病常見癥狀，不及時治療可致惡性心律失常，嚴重者危及生命安全。中醫(yī)在治療心悸有著很好的療效。因此筆者在《中國藥典》中選取心悸進行數(shù)據(jù)挖掘與分析，得到高頻藥物丹參、人參、麥冬、川芎、五味子等。其中丹參，人參，麥冬等藥在Ling和Hao等[15,16]的心悸用藥規(guī)律分析中也是高頻用藥。雖然丹參是本文挖掘出的最高頻藥物，但這些高頻藥物又是如何作用在心律失常疾病的機理上，故本文后續(xù)通過網(wǎng)絡(luò)藥理學(xué)平臺對丹參治療心律失常進行了靶點分析。

3.1 心律失常疾病靶點分析

通過研究心律失常疾病的靶點中，發(fā)現(xiàn)心律失常與離子通道基因表達異常緊密相關(guān)，多個離子通道基因的突變可引起各種心律失常。目前，已知絕大多數(shù)的原發(fā)性心電異常都是由編碼各主要離子通道亞單位的基因突變引起的，統(tǒng)稱“離子通道病”。Nav1.5通道是人類主要的心臟鈉離子通道類型，由SCN5A基因編碼。近年來發(fā)現(xiàn)SCN5A基因突變與房性心律失常，室性心律失常，病態(tài)竇房結(jié)綜合征密切相關(guān)[17]。不僅如此，研究表明像ESR等靶點也能誘發(fā)心律失常[18]。除此之外，微循環(huán)障礙會對心血管造成影響，引發(fā)相關(guān)心血管疾病，因此有研究表明，丹參治療微循環(huán)障礙也主要體現(xiàn)在PTGS2上等[19]，不僅如此，也有相關(guān)研究表明，在丹參治療冠心病相關(guān)靶點中也與PTGS2、PTGS1密切相關(guān)[20]。并且丹參的活性化合物分子與相關(guān)心血管疾病治療密切相關(guān)，丹參通過作用于多個靶點來發(fā)揮藥理學(xué)作用，主要體現(xiàn)在PTGS1、PTGS2、ADRA1A、ADRA1B、CHRM2中[21]。

3.2 丹參治療心律失常的活性化合物預(yù)測

在我國丹參作為活血祛瘀之要藥，在治療心血管疾病已有兩千多年的歷史，丹參活性成分及丹參制劑具有一定的抗動脈粥樣硬化、改善心肌重構(gòu)、抗心肌缺血、抗心律失常的心血管保護作用。但對心律失常的報道還是鮮見，通過本研究中發(fā)現(xiàn)，在預(yù)測的活性化合物中，丹參酮ⅡA通過減少體內(nèi)MDA含量、加強SOD活力從而提高HUVEC細胞抗氧化應(yīng)激損傷能力，降低細胞自噬水平達到保護血管的作用[22]，同時Sun等[23]研究發(fā)現(xiàn)，丹參酮IIA是一種特異性激活劑，可用于預(yù)防和治療由心肌細胞動作電位復(fù)極延遲和IKs通道功能降低引起的心律失常。不僅是丹參酮IIA，在丹參的活性化合物中丹酚酸B同樣重要，研究發(fā)現(xiàn)，丹酚酸B具有抗氧化、抗炎、抗腫瘤、抑制細胞凋亡等多種藥理活性，在臨床可廣泛用于治療心腦血管疾病[24]。因此通過藥理學(xué)發(fā)現(xiàn)丹參的許多化合物都可以預(yù)防甚至治療心律失常。

3.3 丹參治療心律失常的信號通路研究

影響心律失常的原因有很多，其中心肌缺血再灌注后發(fā)生的心律失常是重要的臨床表現(xiàn)形式，因此提高心肌細胞的抗氧化能力，抑制細胞的過度自噬顯得非常重要，在這個過程中，藥理學(xué)研究發(fā)現(xiàn)可以上調(diào)PI3K信號通路的表達，進一步上調(diào)Akt和mTOR的表達，也就是通過PI3K/Akt/mTOR信號通路起到抗心肌缺血再灌注損傷致心律失常的作用[25]。不僅如此，在免疫應(yīng)答方面，T輔助型細胞為T細胞的一種，主要有Th1、Th2及Th17三個亞型。其中Th17主要分泌IL-17、IL-17F等。大量實驗證明IL-17信號通路參與了心臟纖維化的進程[26]，而相關(guān)基礎(chǔ)研究也證實了IL-17參與了AngII誘導(dǎo)的心房纖維化，對心律失常的產(chǎn)生有促進作用。同時對于TNF通路而言，對心律失常相關(guān)研究報道鮮有發(fā)生，但發(fā)現(xiàn)其能夠通過誘導(dǎo)誘導(dǎo)型一氧化氮合酶的表達，從而導(dǎo)致微血管內(nèi)皮細胞凋亡[27]，只能說TNF通路與微循環(huán)障礙的發(fā)生、發(fā)展密切相關(guān)，又因微循環(huán)障礙可能與心血管疾病有關(guān)，故可能對心律失常有影響。

4 結(jié)論

本文通過對《中國藥典》的數(shù)據(jù)挖掘，對治療心悸的中成藥有了系統(tǒng)的分析，并在結(jié)果上通過藥理學(xué)分析，得出丹參治療心悸的具體分子靶點機制，為后續(xù)的臨床選方用藥科研試驗提供參照。但因為篇幅有限，對人參，麥冬等高頻藥物亦可以進行后續(xù)相關(guān)研究，為臨床提供更多證據(jù)。