吳小翠 楊智華 謝毅強(qiáng) 李 瑞 程亞偉

1.海南省中醫(yī)院內(nèi)分泌科，海南海口 570203；2.廣州中醫(yī)藥大學(xué)，廣東廣州 510006；3.海南醫(yī)學(xué)院中醫(yī)學(xué)院，海南海口 571199；4.海南省中醫(yī)院眼科，海南?？?570203；5.海南省中醫(yī)院治未病中心，海南海口 570203

糖尿病及其并發(fā)癥是全球第九大死亡原因，糖尿病視網(wǎng)膜病變（diabetic retinopathy，DR）就是其中最常見(jiàn)的一種[1]。DR 的發(fā)生被認(rèn)為與低度炎癥[2]、氧化應(yīng)激[3]等多種分子途徑機(jī)制有關(guān)，其發(fā)生與發(fā)展與糖尿病血糖控制息息相關(guān)[1]。非增殖期DR（non-proliferative diabetic retinopathy，NPDR）是DR 的早期階段，其特征為視網(wǎng)膜出血、視網(wǎng)膜內(nèi)微血管異常、微動(dòng)脈瘤和棉絮斑等[4]，西醫(yī)治療常在控制血糖的基礎(chǔ)上給予營(yíng)養(yǎng)神經(jīng)、抑制血管生成、促進(jìn)微循環(huán)的藥物，并配合激光治療[5-7]。但常規(guī)的治療手段如糖皮質(zhì)激素往往會(huì)引起眼內(nèi)壓增高、惡心、嘔吐等不良反應(yīng)，而中醫(yī)藥治療NPDR 則著眼于辨證施治、調(diào)理臟腑[8-9]，治療疾病的同時(shí)個(gè)性化地改善患者的不良反應(yīng)，已逐漸成為NPDR治療中不可或缺的一環(huán)。本研究采用數(shù)據(jù)挖掘結(jié)合網(wǎng)絡(luò)藥理學(xué)、分子對(duì)接技術(shù)對(duì)中醫(yī)藥治療NPDR 的中藥方劑進(jìn)行系統(tǒng)分析，挖掘相關(guān)的分子機(jī)制，為后續(xù)的機(jī)制研究打下基礎(chǔ)。

1 資料與方法

1.1 中藥治療NPDR 處方數(shù)據(jù)挖掘

1.1.1 文獻(xiàn)來(lái)源 檢索中國(guó)知網(wǎng)、萬(wàn)方數(shù)據(jù)庫(kù)中有關(guān)中醫(yī)藥治療NPDR 的相關(guān)文獻(xiàn)，檢索期限為數(shù)據(jù)庫(kù)建庫(kù)至2020 年12 月。

1.1.2 納入標(biāo)準(zhǔn)①主題含“糖尿病視網(wǎng)膜病變”“糖網(wǎng)”或“消渴目盲”“消渴”；“中醫(yī)藥”“中草藥”或“中醫(yī)”或“中藥”。②使用口服單純中藥或以口服中藥為主結(jié)合其他方法治療NPDR 的臨床研究。

1.1.3 排除標(biāo)準(zhǔn)①處方組成有缺失；②綜述、單純動(dòng)物實(shí)驗(yàn)的文獻(xiàn)；③個(gè)人經(jīng)驗(yàn)、個(gè)案報(bào)道。

1.1.4 中藥規(guī)范化處理 參考《中藥學(xué)》[10]、《中華人民共和國(guó)藥典》[11]，對(duì)納入文獻(xiàn)中的藥物名稱及性味、歸經(jīng)進(jìn)行規(guī)范化處理。如“法半夏”“姜半夏”“制半夏”統(tǒng)一為“半夏”，“玄胡”“元胡”統(tǒng)一為“延胡索”等。

1.1.5 方劑錄入與核對(duì) 依照上述規(guī)范化處理原則，將整理好的處方一一錄入計(jì)算平臺(tái)的“方劑管理”模塊，一位研究人員錄入，另一位研究人員進(jìn)行二次審核校對(duì)，以確保數(shù)據(jù)資料準(zhǔn)確無(wú)誤。

1.1.6 數(shù)據(jù)分析 采用中醫(yī)傳承計(jì)算平臺(tái)V3.0 的“藥物統(tǒng)計(jì)”對(duì)藥物的使用、性味、歸經(jīng)的頻率、頻次及分布進(jìn)行統(tǒng)計(jì)分析；運(yùn)用“數(shù)據(jù)分析系統(tǒng)”進(jìn)行組方規(guī)律分析及聚類分析，設(shè)置支持度為15%，置信度為0.5。

1.2 高頻藥對(duì)網(wǎng)絡(luò)藥理學(xué)分析

1.2.1 高頻藥對(duì)有關(guān)化學(xué)成分和基因的獲取 將數(shù)據(jù)挖掘所得到使用頻次最高的藥對(duì)組合輸入中醫(yī)藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫(kù)分析平臺(tái)、中醫(yī)藥綜合數(shù)據(jù)庫(kù)，并按照口服利用度≥30%，成藥相似性≥0.18 的篩選條件，在兩個(gè)數(shù)據(jù)庫(kù)中獲取該藥對(duì)組合的主要化學(xué)成分以及可能作用的靶點(diǎn)。將所得的靶點(diǎn)名稱輸入U(xiǎn)niprot數(shù)據(jù)庫(kù)以查找有關(guān)靶點(diǎn)對(duì)應(yīng)的基因名以及編號(hào)。最后，將靶基因和中藥主要化學(xué)成分作映射，構(gòu)建“成分-靶點(diǎn)”圖。

1.2.2 糖尿病視網(wǎng)膜病變的潛在靶點(diǎn)收集 在GeneCards數(shù)據(jù)庫(kù)以“non-proliferative diabetic retinopathy”“diabetic retinopathy”為主要檢索詞搜索與糖尿病視網(wǎng)膜病變有關(guān)的疾病作用靶點(diǎn)，選擇Score 值≥10 的疾病靶點(diǎn)與高頻藥對(duì)中的相關(guān)基因進(jìn)行交叉，得到高頻藥對(duì)治療糖尿病視網(wǎng)膜病變的作用基因。

1.2.3 蛋白質(zhì)-蛋白質(zhì)相互作用（protein-protein interaction，PPI）網(wǎng)絡(luò)的構(gòu)建 將高頻藥對(duì)治療糖尿病視網(wǎng)膜病變的作用基因輸入STRING 數(shù)據(jù)庫(kù)（https://string-db.org/），選擇“Homo sapiens”類別后，設(shè)置置信度為0.4。將所獲得的分析結(jié)果以tsv 文件形式導(dǎo)入Cytoscape 3.7.2 中，根據(jù)Cytoscape 軟件的程序包進(jìn)行計(jì)算，篩選出節(jié)點(diǎn)度值較大的靶點(diǎn)作為PPI 核心網(wǎng)絡(luò)。

1.2.4 KEGG 和GO 富集分析 將高頻藥對(duì)治療糖尿病視網(wǎng)膜病變的作用基因輸入R 軟件（R 4.0.2），調(diào)用clusterProfiler 包，對(duì)這些基因從細(xì)胞組分、分子功能、生物過(guò)程3 個(gè)方面進(jìn)行注釋。GO 分析是對(duì)基因本身功能的注釋，而KEGG 分析則是對(duì)于基因參與人體活動(dòng)的相關(guān)通路注釋。利用R 語(yǔ)言程序中的cluster-Profiler 程序包對(duì)關(guān)鍵靶點(diǎn)基因進(jìn)行GO 富集分析。GO 結(jié)果均以條形圖格式輸出。將關(guān)鍵靶點(diǎn)基因用R語(yǔ)言程序中的pathview 程序包進(jìn)行KEGG 富集通路分析，以基因比率為橫軸、通路名稱為縱軸、氣泡顏色代表P 值，繪制氣泡圖。

1.2.5 分子對(duì)接技術(shù) 在RSCB-PDB 數(shù)據(jù)庫(kù)（http://www.rcsb.org/）查找篩選出的排名最高的晶體結(jié)構(gòu)，利用ZINC 數(shù)據(jù)庫(kù)（http://zinc.docking.org/）查找篩選出的口服利用度及成藥相似度最高的化學(xué)成分的Mol2 結(jié)構(gòu)式，將活性成分與靶標(biāo)一并導(dǎo)入AutoDock 軟件，對(duì)活性成分和靶標(biāo)進(jìn)行虛擬分子對(duì)接分析。一般認(rèn)為，結(jié)合能＜0 的分子即可自行結(jié)合，而數(shù)值越小結(jié)合力越強(qiáng)，既往同行研究提示，結(jié)合能≤-5 kJ/mol 可視為該化學(xué)成分可與作用靶標(biāo)有效結(jié)合[12]。

2 結(jié)果

2.1 中藥頻次

通過(guò)文獻(xiàn)篩選獲得107 首處方，共涉及151 味中藥，累計(jì)頻次達(dá)936 次。其中頻次≥10 次的有27 味，見(jiàn)表1。

2.2 中藥功效、歸經(jīng)及四氣五味

151 味中藥共涵蓋10 種功效，其中使用頻次較高的包括補(bǔ)虛類、清熱類、活血化瘀類。四氣以寒、平、溫性為主；五味主要以苦、甘味為主；歸經(jīng)主要?dú)w為肝、肺、腎經(jīng)。見(jiàn)圖1。

圖1 中藥功效、歸經(jīng)、四氣、五味統(tǒng)計(jì)圖

2.3 高頻中藥關(guān)聯(lián)規(guī)則分析

藥物核心組合的頻次及關(guān)聯(lián)程度見(jiàn)表2～3、圖2，其中中藥組合黃芪-生地出現(xiàn)頻次最高。

圖2 藥物關(guān)系網(wǎng)絡(luò)

表2 藥物核心組合頻次（前10 位）

表3 藥物核心組合關(guān)聯(lián)程度（前10 位）

2.4 高頻藥對(duì)黃芪-生地治療NPDR 的網(wǎng)絡(luò)藥理學(xué)分析

2.4.1 藥物成分及靶點(diǎn)篩選、NPDR 關(guān)鍵靶點(diǎn)篩選及網(wǎng)絡(luò)分析 得到黃芪、生地的活性成分共26 個(gè)，以及對(duì)應(yīng)的靶點(diǎn)1 280 個(gè)。剔除重復(fù)項(xiàng)和無(wú)對(duì)應(yīng)活性成分的靶點(diǎn)后，得到黃芪、生地的活性成分靶點(diǎn)共497 個(gè)。將NPDR 相關(guān)的疾病靶點(diǎn)與黃芪、生地的活性成分靶點(diǎn)匹配映射，得到共同靶點(diǎn)47 個(gè)，即黃芪-生地藥對(duì)治療NPDR 的關(guān)鍵靶點(diǎn)。將中藥活性成分、關(guān)鍵靶點(diǎn)導(dǎo)入Cytoscape 軟件，構(gòu)建中藥活性成分-關(guān)鍵靶點(diǎn)網(wǎng)絡(luò)，見(jiàn)圖3。網(wǎng)絡(luò)中有61 個(gè)節(jié)點(diǎn)（包括14 個(gè)活性成分及47 個(gè)關(guān)鍵靶點(diǎn)）。在該網(wǎng)絡(luò)中，藥物活性成分與關(guān)鍵靶點(diǎn)之間的連線代表該活性成分可能與這個(gè)關(guān)鍵靶點(diǎn)有聯(lián)系，連線越多提示該成分在治療NPDR 的過(guò)程中發(fā)揮的作用越大。連線較多的藥物活性成分有槲皮素、山柰酚、刺芒柄花素、異鼠李素、7-O-methylisomucronulatol 等。

圖3 活性成分-關(guān)鍵靶點(diǎn)網(wǎng)絡(luò)圖

2.4.2 關(guān)鍵靶點(diǎn)PPI 網(wǎng)絡(luò)分析 關(guān)鍵靶點(diǎn)的PPI 網(wǎng)絡(luò)圖見(jiàn)圖4。PPI 核心網(wǎng)絡(luò)圖中，節(jié)點(diǎn)度值越大提示該蛋白在PPI 網(wǎng)絡(luò)中的位置越重要，得到黃芪-生地藥對(duì)治療NPDR 的核心靶點(diǎn)有SPP1、TP53、SERPINE1、PTGS2、FOS、CCND1、VEGFA、TNF 等，見(jiàn)圖5。

圖4 蛋白質(zhì)-蛋白質(zhì)相互作用網(wǎng)絡(luò)圖

圖5 蛋白質(zhì)-蛋白質(zhì)相互作用核心網(wǎng)絡(luò)圖

2.4.3 GO 與KEGG 富集分析 關(guān)鍵靶點(diǎn)GO 富集分析結(jié)果顯示，這些關(guān)鍵靶點(diǎn)基因涉及1 942 條通路，靶點(diǎn)基因數(shù)目前30 條繪制成柱狀圖見(jiàn)圖6，關(guān)鍵靶點(diǎn)KEGG 通路分析見(jiàn)圖7，可見(jiàn)黃芪-生地藥對(duì)治療NPDR 的主要通路包括有細(xì)胞對(duì)化學(xué)應(yīng)激的反應(yīng)、生殖結(jié)構(gòu)發(fā)育、生殖系統(tǒng)發(fā)育、上皮細(xì)胞增殖、細(xì)胞對(duì)脂多糖反應(yīng)等。

圖6 GO 富集通路柱狀圖

圖7 KEGG 富集通路氣泡圖

2.4.4 分子對(duì)接分析 根據(jù)“藥物-活性成分-疾病-靶點(diǎn)”網(wǎng)絡(luò)中交集靶點(diǎn)的運(yùn)算結(jié)果，篩選出黃芪-生地藥對(duì)中的核心活性成分有兩個(gè)，分別是槲皮素、山柰酚，其成分化學(xué)結(jié)構(gòu)見(jiàn)圖8。以網(wǎng)絡(luò)連線數(shù)（分值）為主要依據(jù)，篩選得到4 個(gè)核心靶點(diǎn)蛋白VEGFA、CASP3、ICAM1、SERPINE1，其關(guān)聯(lián)性均0.595 588 235，網(wǎng)絡(luò)分值及網(wǎng)絡(luò)連線數(shù)均為19，特征向量均為0.228 793 263。對(duì)核心活性成分與核心靶點(diǎn)蛋白進(jìn)行分子對(duì)接。核心靶點(diǎn)蛋白和核心活性成分對(duì)接的自由能計(jì)算結(jié)果顯示，槲皮素與VEGFA、ICAM1、SERPINE1的結(jié)合能分別為-8.1、-6.7、-6.7 kJ/mol，山柰酚與CASP3 的結(jié)合能為-7.3 kJ/mol，提示其與NPDR 中的關(guān)鍵核心靶點(diǎn)有良好的結(jié)合作用。見(jiàn)圖9～10。

圖8 黃芪-生地藥對(duì)核心活性成分化學(xué)結(jié)構(gòu)

圖9 槲皮素與 與VEGFA、ICAM1、SERPINE1 的對(duì)接示意圖

3 討論

NPDR 是由于血糖控制不佳引起的糖尿病常見(jiàn)微血管并發(fā)癥之一，屬于中醫(yī)范疇中的“消渴目病”范疇，多因?yàn)榕K腑功能失調(diào)，氣血痹阻、陰陽(yáng)平衡失調(diào)導(dǎo)致消渴病，而消渴病日久遷延至雙目所致，中醫(yī)藥對(duì)本病的治療效果確切[13]。為更好地總結(jié)近年來(lái)各家用中醫(yī)藥治療NPDR 的學(xué)術(shù)思想與用藥經(jīng)驗(yàn)，本研究同時(shí)采用數(shù)據(jù)挖掘與網(wǎng)絡(luò)藥理學(xué)分析方法，探討中藥治療NPDR 的一般規(guī)律和關(guān)鍵藥物的作用機(jī)制，可以為中醫(yī)藥治療NPDR 提供一個(gè)重要的參考依據(jù)。

圖10 山柰酚與CASP3 對(duì)接示意圖

對(duì)處方進(jìn)行功效、歸經(jīng)分析，分析結(jié)果提示最常用的中藥為補(bǔ)虛藥，其次為清熱藥與活血化瘀藥，主要?dú)w肝、脾、腎經(jīng)。對(duì)納入的151 味中藥進(jìn)行關(guān)聯(lián)分析得到治療NPDR 的核心中藥，其中包括黃芪、生地、葛根、三七等。其中黃芪益氣固表，養(yǎng)血生肌；生地黃清熱生津、涼血止血；葛根清熱解毒，生肌止渴；三七活血化瘀。其中黃芪常與生地配伍共同發(fā)揮益氣養(yǎng)陰，清熱生津，涼血止血的功效。中醫(yī)認(rèn)為NPDR 的發(fā)病病機(jī)為氣陰兩虛、陰虛火旺或氣血瘀滯，提示治療應(yīng)當(dāng)以益氣養(yǎng)陰、補(bǔ)益肝腎、清熱涼血、化瘀通絡(luò)為主[14]，配以益氣健脾[15]。四氣五味統(tǒng)計(jì)結(jié)果提示治療NPDR的中藥多以寒、溫、平及甘、苦為主，甘味平性藥多為益氣健脾補(bǔ)虛中藥，性寒味苦的中藥多具有疏風(fēng)清熱涼血功效，辛溫的中藥多能疏肝健脾，因此在臨床上較為常用。目前中醫(yī)各家比較一致地認(rèn)為，消渴目病是由于先天稟賦不足或后天失養(yǎng)，導(dǎo)致肝腎陰虛，陰虛生火，虛火灼精，血液黏稠成瘀，精、氣、血不能上養(yǎng)于目所致[16]。因此NPDR 的中醫(yī)治療以標(biāo)本兼治為原則，治本以益氣養(yǎng)陰為主，治標(biāo)則當(dāng)活血化瘀。

治療NPDR 的核心中藥藥對(duì)黃芪-生地中的活性成分包括槲皮素、山柰酚、刺芒柄花素、異鼠李素、7-O-methylisomucronulatol 等。槲皮素能有效減輕高糖條件下人腎小球內(nèi)皮細(xì)胞損傷，從而抑制糖尿病進(jìn)展，這一作用與激活PI3K/Akt 通路調(diào)節(jié)細(xì)胞活力、細(xì)胞凋亡、凋亡基因mRNA 表達(dá)量有關(guān)[17-18]；亦可通過(guò)調(diào)控細(xì)胞自噬抑制高糖誘導(dǎo)的脈絡(luò)膜和視網(wǎng)膜血管生成[19]，對(duì)延緩、減輕DR 具有一定作用。山柰酚對(duì)高糖條件下人腎小球內(nèi)皮細(xì)胞氧化應(yīng)激及凋亡具有抑制作用，可控制糖尿病腎損傷[20]。徐曉紅等[21]的研究發(fā)現(xiàn)，山柰酚通過(guò)靶向VEGF 和PGF 基因抑制微血管發(fā)生，抑制Src-Akt1-Erk1/2 信號(hào)通路的激活，可能是治療DR 的潛在藥物。刺芒柄花素屬于黃芪總黃酮的一種，通過(guò)抑制缺氧ARPE-19 細(xì)胞分泌血管內(nèi)皮生長(zhǎng)因子，防止體內(nèi)缺氧誘導(dǎo)的視網(wǎng)膜血管新生[22]。異鼠李素能有效抑制氧化應(yīng)激損傷[23]，這可能對(duì)高糖狀態(tài)下的視網(wǎng)膜損害具有保護(hù)作用。PPI 網(wǎng)絡(luò)中黃芪-生地治療NPDR 的核心靶點(diǎn)包括SPP1、TP53、SERPINE1、PTGS2、FOS、MAPK1、VEGFA、TNF 等。其中TNF 基因編碼一種多功能促炎性細(xì)胞因子，屬于腫瘤壞死因子超家族，這種細(xì)胞因子主要由巨噬細(xì)胞分泌，參與調(diào)節(jié)胰島素抵抗，在糖尿病發(fā)展中起到不可忽視的作用[24]。SERPINE1 基因編碼絲氨酸蛋白酶抑制劑超家族成員，絲氨酸蛋白酶抑制劑具有增加胰島素敏感性、調(diào)節(jié)脂類代謝及調(diào)節(jié)血糖代謝等功能[25]，可通過(guò)調(diào)控TNF、SERPINE1 表達(dá)量，減輕胰島素抵抗，繼而延緩糖尿病進(jìn)展及并發(fā)癥的發(fā)生。VEGFA 基因編碼一種肝素結(jié)合蛋白，以二硫鍵連接的同二聚體存在，這種生長(zhǎng)因子誘導(dǎo)血管內(nèi)皮細(xì)胞增殖和遷移，參與視網(wǎng)膜病理血管生成[26]，對(duì)改善視網(wǎng)膜血管病變可能有直接作用。TP53 調(diào)節(jié)細(xì)胞代謝（包括細(xì)胞自噬及細(xì)胞成分的分解和重塑），通過(guò)這些代謝活動(dòng)維持細(xì)胞內(nèi)環(huán)境穩(wěn)定和支持細(xì)胞生長(zhǎng)[27]，其信號(hào)通路影響線粒體功能和視網(wǎng)膜色素上皮細(xì)胞增殖[28]。黃芪-生地藥對(duì)的核心有效成分與核心靶點(diǎn)蛋白進(jìn)行模擬對(duì)接的結(jié)合能＜-5 kJ/mol，提示有效成分和靶點(diǎn)蛋白之間的結(jié)合較為緊密，黃芪-生地藥對(duì)對(duì)NPDR 有潛在治療效果，結(jié)合網(wǎng)絡(luò)藥理學(xué)分析結(jié)果，后續(xù)可進(jìn)一步研究山柰酚、槲皮素、異鼠李素、刺芒柄花素等干預(yù)NPDR 的作用與機(jī)制，過(guò)表達(dá)或敲除相關(guān)核心靶點(diǎn)研究其上下游關(guān)系，探索中藥治療NDPR 新機(jī)制。

黃芪-生地參與治療NPDR 的主要通路可能包括細(xì)胞對(duì)化學(xué)應(yīng)激的反應(yīng)、生殖結(jié)構(gòu)發(fā)育、細(xì)胞對(duì)脂多糖反應(yīng)、上皮細(xì)胞增殖等。氧化應(yīng)激損傷胰島β 細(xì)胞和降低外周組織對(duì)胰島素的敏感性，導(dǎo)致糖尿病的發(fā)生發(fā)展[3]；生物活性氧包括O2、OH-、H2O2、NO 等，是線粒體氧化代謝的產(chǎn)物，為氧化應(yīng)激的主要標(biāo)志物?；钚匝踉谔悄虿〖安l(fā)癥機(jī)制中的作用主要是對(duì)組織和細(xì)胞的毒性損傷，還是重要的細(xì)胞內(nèi)信使，能夠活化許多信號(hào)傳導(dǎo)通路，間接損傷組織和細(xì)胞[29-30]。上皮細(xì)胞增殖同樣屬于糖尿病視網(wǎng)膜病損的重要因素，在高糖狀態(tài)下，包括視網(wǎng)膜色素上皮細(xì)胞在內(nèi)的細(xì)胞更新（凋亡、增殖、肥大）模式會(huì)隨著線粒體活性氧和氧化應(yīng)激的增加，最終導(dǎo)致視網(wǎng)膜血管生成和組織重塑[3]。

綜上所述，中醫(yī)藥治療NPDR 以益氣養(yǎng)陰、補(bǔ)益肝腎、清熱涼血、化瘀通絡(luò)為常規(guī)治則，用藥以補(bǔ)虛最為多見(jiàn)，次而配伍清熱涼血。文獻(xiàn)處方中的核心組合黃芪-生地，可能主要通過(guò)槲皮素、山柰酚、刺芒柄花素等活性成分，作用于SPP1、TP53、SERPINE1、PTGS2、FOS、CCND1、VEGFA、TNF 等相關(guān)疾病靶點(diǎn)，調(diào)控對(duì)化學(xué)應(yīng)激的反應(yīng)、生殖結(jié)構(gòu)發(fā)育、生殖系統(tǒng)發(fā)育、上皮細(xì)胞增殖等相關(guān)細(xì)胞信號(hào)通路，共同發(fā)揮NPDR 的治療作用，體現(xiàn)了中醫(yī)藥治療針對(duì)多靶點(diǎn)、多通路的特點(diǎn)。但是本研究?jī)H在理論上研究中醫(yī)藥治療NPDR 的一般規(guī)律和相關(guān)機(jī)制，后續(xù)仍需要更多臨床實(shí)踐和進(jìn)行相關(guān)實(shí)驗(yàn)予以驗(yàn)證。