趙元 李香斌 張尚斌 陳劍平*

1. 廣州中醫(yī)藥大學(xué)第四臨床醫(yī)學(xué)院，廣東 深圳 518033

2. 深圳市中醫(yī)院/深圳市醫(yī)院中藥制劑研究重點(diǎn)實(shí)驗(yàn)室，廣東 深圳 518033

3. 北京中醫(yī)藥大學(xué)深圳醫(yī)院，廣東 深圳 518172

2 型糖尿病性骨質(zhì)疏松(type 2 diabetic osteoporosis，T2DOP)是糖尿病常見(jiàn)的并發(fā)癥之一，是一種全身性代謝性骨病。T2DOP的發(fā)病率隨糖尿病患者群體的不斷擴(kuò)大而逐年升高，而中國(guó)糖尿病患者骨質(zhì)疏松的發(fā)病率為40%～66%[1]。T2DOP患者因骨脆性和微血管并發(fā)癥導(dǎo)致的微裂紋而極易引起骨折[2]，不僅嚴(yán)重影響患者生活質(zhì)量，而且給家庭和社會(huì)帶來(lái)沉重的經(jīng)濟(jì)負(fù)擔(dān)?，F(xiàn)代醫(yī)學(xué)的治療原則主要為控制原發(fā)病和補(bǔ)充鈣劑及維生素D[3]，但效果不顯著。中醫(yī)藥作用的多靶點(diǎn)性在臨床治療多因素影響的復(fù)雜的T2DOP有一定的優(yōu)勢(shì)。在傳統(tǒng)中醫(yī)學(xué)中，T2DOP屬于“消渴”合并“骨痿”范疇，腎虛和脾虛是其重要病機(jī)[1]。山藥和黃芪始載于《神農(nóng)本草經(jīng)》并列為上品，山藥，補(bǔ)脾養(yǎng)陰，生津益肺，強(qiáng)腎固精；黃芪，補(bǔ)氣升陽(yáng)，益水之源，生津養(yǎng)血。二藥相配補(bǔ)脾固腎，益氣生津，二者一陰一陽(yáng)，相互促進(jìn)，一則使脾氣升，散精達(dá)肺，輸布津液以止渴，二則使腎氣固，封藏精微以生髓，且山藥長(zhǎng)肌肉堅(jiān)強(qiáng)筋骨，治療“消渴”并“骨痿”有標(biāo)本兼治的效果。二者均是糖尿病合并骨質(zhì)疏松的常用藥物[4～6]，且現(xiàn)代藥理研究[7]顯示山藥能通過(guò)調(diào)控Wnt和p38信號(hào)通路而具有抑制骨丟失作用，作為PPARγ的激動(dòng)劑、AMPK激活劑而治療2 型糖尿病[8]。黃芪甲苷通過(guò)調(diào)節(jié)PI3K/Akt/FoxO1通路抑制糖尿病大鼠肝糖異生[9-10]，黃芪多糖可直接促進(jìn)成骨細(xì)胞質(zhì)增殖、提高成骨細(xì)胞活性[11]，但二者配伍治療T2DOP的作用機(jī)制卻未見(jiàn)研究。本研究通過(guò)網(wǎng)絡(luò)藥理學(xué)研究的方法探討黃芪和山藥配伍治療T2DOP的作用機(jī)制，為日后兩者配伍治療T2DOP研究提供方向。

1 材料與方法

1.1 材料

1.1.1“山藥-黃芪”藥對(duì)化學(xué)成分：在中藥系統(tǒng)藥理數(shù)據(jù)庫(kù)和分析平臺(tái)(TCMSP, http://tcmspw.com/tcmsp.php)檢索山藥、黃芪對(duì)應(yīng)的化學(xué)成分，篩選條件為口服生物利用度(oral bioavailability, OB)>30%且類(lèi)藥性(drug likeness, DL)>0.18。

1.1.2T2DOP疾病相關(guān)基因：以“Type 2 diabetic osteoporosis”為關(guān)鍵詞，在Gene cards和 OMIM數(shù)據(jù)庫(kù)中檢索T2DOP相關(guān)基因，去除重復(fù)靶點(diǎn)，匯總得到與T2DOP相關(guān)的疾病靶點(diǎn)。

1.1.3“山藥-黃芪”藥對(duì)的活性成分靶標(biāo)：在TCMSP(http://tcmspw.com/tcmsp.php)及Swiss Target Prediction(http://www.swisstargetprediction.ch)數(shù)據(jù)庫(kù)中查詢(xún)1.1.1中所得化學(xué)成分的預(yù)測(cè)靶標(biāo)。并在Uniprot 數(shù)據(jù)庫(kù)(https://www.uniprot.org)中將預(yù)測(cè)靶標(biāo)轉(zhuǎn)換為基因名，獲得靶點(diǎn)的Uniprot ID。

1.2 方法

1.2.1構(gòu)建“山藥黃芪化學(xué)成分-山藥黃芪候選靶標(biāo)-T2DOP相關(guān)靶標(biāo)”網(wǎng)絡(luò)：利用Cytoscape 3.7.1軟件構(gòu)建“化學(xué)成分-預(yù)測(cè)靶標(biāo)-T2DOP相關(guān)靶標(biāo)”網(wǎng)絡(luò)，利用 Cytoscape軟件中Network Analyzer插件進(jìn)行網(wǎng)絡(luò)分析，并獲取網(wǎng)絡(luò)拓?fù)鋮?shù)，篩選大于度中心性、接近中心性、中介中心性/中間中心性三者均值的節(jié)點(diǎn)。

1.2.2基因功能估計(jì)分析：將篩選節(jié)點(diǎn)中的靶標(biāo)通過(guò)David 6.8數(shù)據(jù)庫(kù)(https://david.ncifcrf.gov/)進(jìn)行GO生物學(xué)富集分析，包括生物過(guò)程(biological process, BP)，分子功能(molecular function, MF)和細(xì) 胞 組 成(cellular component, CC)3 個(gè)部分，設(shè)定閾值P<0. 05，結(jié)果以氣泡圖形式進(jìn)行可視化展示。

1.2.3KEGG富集分析及可視化：運(yùn)用R語(yǔ)言軟件對(duì)關(guān)鍵靶點(diǎn)進(jìn)行信號(hào)通路富集分析，挑選與T2DOP相關(guān)的信號(hào)通路，最終設(shè)定閾值P<0. 05。

1.2.4“山藥黃芪化學(xué)成分-關(guān)鍵作用靶點(diǎn)-信號(hào)通路”網(wǎng)絡(luò)的構(gòu)建：利用Cytoscape 3.7.1軟件將1.2.1中得到的網(wǎng)絡(luò)節(jié)點(diǎn)及1.2.3中富集的信號(hào)通路構(gòu)成“山藥黃芪化學(xué)成分-關(guān)鍵作用靶點(diǎn)-信號(hào)通路”網(wǎng)絡(luò)，通過(guò)計(jì)算該網(wǎng)絡(luò)的拓?fù)涮卣髦?，判斷山藥黃芪治療T2DOP的關(guān)鍵化學(xué)成分、作用靶標(biāo)及信號(hào)通路。

2 結(jié)果

2.1 “山藥-黃芪”藥對(duì)中活性組分的篩選

在TCMSP中查詢(xún)到符合篩選條件的山藥活性成分16 個(gè)、黃芪20 個(gè)，最終共篩選出36個(gè)化合物，見(jiàn)表1。

表1 “山藥-黃芪”藥對(duì)的候選活性成分信息

2.2 T2DOP疾病靶基因的篩選結(jié)果

在Gene cards網(wǎng)站收集得到T2DOP疾病靶標(biāo)共3 180 個(gè)，OMIM網(wǎng)站收集得到183 個(gè)T2DOP靶標(biāo)，合并2 個(gè)數(shù)據(jù)庫(kù)數(shù)據(jù)，剔除重合基因，共檢索出3 305 個(gè)相關(guān)基因信息。

2.3 “山藥-黃芪”藥對(duì)中活性成分的靶標(biāo)預(yù)測(cè)

將篩選出的36 個(gè)成分分別在TCMSP平臺(tái)和Swiss Target Prediction網(wǎng)站中查詢(xún)預(yù)測(cè)的靶點(diǎn)，并通過(guò)Uniprot數(shù)據(jù)庫(kù)將蛋白靶點(diǎn)轉(zhuǎn)化為靶點(diǎn)基因，其中山藥16 個(gè)成分，370 個(gè)靶基因；黃芪20 個(gè)成分，547 個(gè)靶基因，最終共有36 個(gè)化合物，去除重復(fù)靶基因后共698條靶基因信息。

2.4 “山藥黃芪化學(xué)成分-山藥黃芪候選靶基因-T2DOP相關(guān)靶基因”網(wǎng)絡(luò)的構(gòu)建

利用Cytoscape 3.7.1軟件構(gòu)建“山藥黃芪化學(xué)成分-山藥黃芪候選靶基因-T2DOP相關(guān)靶基因”網(wǎng)絡(luò)，見(jiàn)圖1。網(wǎng)絡(luò)中各個(gè)節(jié)點(diǎn)之間相互作用，根據(jù)網(wǎng)絡(luò)拓?fù)鋵W(xué)結(jié)構(gòu)分析結(jié)果，共512個(gè)節(jié)點(diǎn)，以大于該網(wǎng)絡(luò)的度中心性均值7.286 0、接近中心性0.316 8、中介中心性/中間中心性均值0.004 1為篩選條件篩選出關(guān)鍵節(jié)點(diǎn)67 個(gè)，其中含有靶標(biāo)40 個(gè)，MET，CYP19A1，AR，PTGS1，ADORA2A，PTGS2，CDK2，GSK3B，NCOA2，ESR2，KDR，ESR1，PTPN1，SRC，RXRA，EGFR，CASP3，GABRA1，CCNA2，PARP1，CA2，ALOX12，PGR，ADRB2，NOS2，PPARG，SYK，ABL1，MAPK14，CCNB1，BACE1共31 個(gè)靶標(biāo)為T(mén)2DOP疾病靶基因，此外尚有ACHE、ADORA1、ADRA1B、CA14、CA9、CDK1、CHRM1、HSD17B2、MAOB 9 個(gè)關(guān)鍵靶標(biāo)。

2.5 GO 功能分析結(jié)果

使用David 6.8數(shù)據(jù)庫(kù)對(duì)山藥黃芪治療T2DOP的40 個(gè)關(guān)鍵靶點(diǎn)進(jìn)行GO富集分析，篩選P<0.05的條目，共有87 個(gè)GO生物學(xué)條目，其中有56 個(gè)BP條目，21 個(gè)MF條目，10 個(gè)CC條目，生物過(guò)程包括肽基絲氨酸磷酸化、調(diào)節(jié)細(xì)胞增殖、跨膜受體蛋白酪氨酸激酶信號(hào)通路、視黃酸的反應(yīng)、轉(zhuǎn)錄、ERK1和ERK2級(jí)聯(lián)的正調(diào)控、子宮內(nèi)胚胎發(fā)育、RNA聚合酶II啟動(dòng)子轉(zhuǎn)錄的正調(diào)控、炎癥反應(yīng)、凋亡過(guò)程的負(fù)調(diào)控等；分子功能涉及ATP結(jié)合、鋅離子結(jié)合、類(lèi)固醇激素受體活性、序列特異性DNA結(jié)合、碳酸鹽脫水酶活性等。富集于細(xì)胞核、細(xì)胞膜、細(xì)胞溶質(zhì)、細(xì)胞質(zhì)核周區(qū)域、高爾基體、質(zhì)膜胞質(zhì)側(cè)的外源成分、受體復(fù)合體、核染色質(zhì)、細(xì)胞連接、RNA聚合酶II轉(zhuǎn)錄因子復(fù)合物等區(qū)域。見(jiàn)圖2～圖4。

圖2 “山藥-黃芪”抗T2DOP潛在靶點(diǎn)生物過(guò)程的GO分析前10 位

圖3 山藥-黃芪抗T2DOP潛在靶點(diǎn)分子功能的GO分析

圖4 山藥-黃芪抗T2DOP潛在靶點(diǎn)細(xì)胞組分的GO分析

2.6 KEGG信號(hào)通路富集分析與化合物-靶點(diǎn)-通路互相作用網(wǎng)絡(luò)構(gòu)建

利用R語(yǔ)言對(duì)40 個(gè)關(guān)鍵靶標(biāo)基因進(jìn)行KEGG 通路分析。在通路富集分析中，得到102條信號(hào)通路，其中與T2DOP相關(guān)的17 條通路，分別為VEGF信號(hào)通路(VEGF signaling pathway)、甲狀腺素信號(hào)通路(thyroid hormone signaling pathway)、松弛素信號(hào)通路(relaxin signaling pathway)、雌激素信號(hào)通路(estrogen signaling pathway)、內(nèi)分泌抵抗(endocrine resistance)、催乳素信號(hào)通路(prolactin signaling pathway)、細(xì)胞衰老(cellular senescence)、卵巢類(lèi)固醇生成(ovarian steroidogenesis)、血小板活化(platelet activation)、黏著斑(focal adhesion)、FoxO信號(hào)通路(FoxO signaling pathway)、脂肪細(xì)胞中脂解的調(diào)節(jié)(regulation of lipolysis in adipocytes)、p53信號(hào)通路(p53 signaling pathway)、PI3K-Akt信號(hào)通路(PI3K-Akt signaling pathway)、鈣信號(hào)通路(calcium signaling pathway)、IL-17信號(hào)通路(IL-17 signaling pathway)、花生四烯酸代謝(arachidonic acid metabolism)。見(jiàn)圖5。

將篩選出的68 個(gè)節(jié)點(diǎn)及17 條T2DOP相關(guān)的KEGG通路導(dǎo)入Cytoseape軟件構(gòu)建“化合物-靶基因-信號(hào)通路”互相作用網(wǎng)絡(luò)。利用Cytoscape 3.7.1軟件構(gòu)建“山藥-黃芪”藥對(duì)中化學(xué)成分所對(duì)應(yīng)的靶點(diǎn)與通路，見(jiàn)圖5。網(wǎng)絡(luò)中各個(gè)節(jié)點(diǎn)之間相互作用，根據(jù)網(wǎng)絡(luò)拓?fù)鋵W(xué)結(jié)構(gòu)分析結(jié)果，共84 個(gè)節(jié)點(diǎn)，該網(wǎng)絡(luò)的度中心性均值26.928 6、接近中心性0.442 4、中介中心性/中間中心性均值0.015 8的2倍為篩選條件篩選出關(guān)鍵節(jié)點(diǎn)28 個(gè)，其中含有靶標(biāo)15 個(gè)，均為T(mén)2DOP疾病靶基因，按度值排序?yàn)镃DK2、PTGS2、PTGS1、GSK3B、SRC、CYP19A1、EGFR、ESR1、ESR2、AR、MAPK14、ADORA2A、NCOA2、KDR、NOS2；13 個(gè)節(jié)點(diǎn)為化合物，按度值排序依次為山奈酚、異鼠李素、槲皮素、(6aR,11aR)-3-羥基-9,10-二甲氧基紫檀烷、7-O-甲基異戊醇、3,9-di-O-美迪紫檀素、海風(fēng)藤烯酮、芒柄花黃素、熊竹素、毛蕊異黃酮、異黃烷酮、蓽茇明寧堿、(3R)-3-(2-hydroxy-3,4-dimethoxyphenyl)chroman-7-ol。文獻(xiàn)研究表明，山奈酚[12-13 ]對(duì)2 型糖尿病及骨質(zhì)疏松均有調(diào)節(jié)作用，槲皮素[14]具有降糖作用，異鼠李素[15]可能抑制經(jīng)典N(xiāo)F-κB通路對(duì)破骨細(xì)胞的分化具有抑制作用, 毛蕊異黃酮、芒柄花黃素可通過(guò)調(diào)節(jié)骨形態(tài)發(fā)生蛋白2的表達(dá)促進(jìn)成骨細(xì)胞的成骨作用[16]，說(shuō)明山奈酚、異鼠李素、槲皮素等13 個(gè)化合物可能是山藥黃芪藥對(duì)發(fā)揮抗T2DOP的活性化合物。度值排前三位的信號(hào)通路分別為PI3K-Akt信號(hào)通路、鈣信號(hào)通路及雌激素信號(hào)通路，說(shuō)明這些信號(hào)通路在調(diào)節(jié)T2DOP中起關(guān)鍵作用。

圖5 “山藥-黃芪”藥對(duì)治療T2DOP的 KEGG 信號(hào)通路分析

圖6 “山藥-黃芪”藥對(duì)中化學(xué)成分所對(duì)應(yīng)的靶點(diǎn)與通路

3 討論

2 型糖尿病(T2DM)是一種以缺乏胰島素和高糖為特征的慢性代謝性疾病，占糖尿病患者的90%以上，可引起骨質(zhì)疏松，導(dǎo)致T2DOP。T2DOP在T2DM并發(fā)癥中病發(fā)率最高可達(dá)60%[17-18]，西醫(yī)主要以對(duì)癥治療為主，然而糖尿病合并骨質(zhì)疏松的發(fā)病機(jī)制比較復(fù)雜，早期癥狀不明顯，部分降糖西藥亦會(huì)影響到骨密度[19-20]，使其在防治中加大了難度。中醫(yī)藥作用的多靶點(diǎn)性在臨床治療復(fù)雜發(fā)病機(jī)制的T2DOP方面有一定的優(yōu)勢(shì)，網(wǎng)絡(luò)藥理學(xué)以生物學(xué)網(wǎng)絡(luò)節(jié)點(diǎn)的連接和關(guān)系分析為基礎(chǔ)構(gòu)建“化學(xué)成分-預(yù)測(cè)靶標(biāo)-信號(hào)通路”的多層次網(wǎng)絡(luò)，然后綜合分析中藥復(fù)方的作用，既反映了中藥復(fù)方多成分、多靶點(diǎn)、多通路的特點(diǎn)，又與疾病結(jié)合起來(lái)，與中醫(yī)整體觀、辨證論治的基本原則相吻合。

通過(guò)檢索TCMSP 數(shù)據(jù)庫(kù)，山藥-黃芪藥對(duì)中含有OB>30 %且DL>0.18 %的化學(xué)成分36個(gè)，經(jīng)過(guò)進(jìn)一步檢索得到預(yù)測(cè)作用靶點(diǎn)698 個(gè)，T2DOP共得到3 305 個(gè)相關(guān)靶點(diǎn)，山藥黃芪化學(xué)成分-山藥黃芪候選靶基因與T2DOP相關(guān)靶標(biāo)作用網(wǎng)絡(luò)圖中以大于該網(wǎng)絡(luò)的度中心性、接近中心性、中介中心性/中間中心性均值為篩選條件篩選出關(guān)鍵節(jié)點(diǎn)67 個(gè)，其中含有靶標(biāo)40 個(gè)。以40 個(gè)靶標(biāo)進(jìn)行GO分析及KEGG 通路分析，結(jié)果表明，GO分析中有56個(gè)BP條目，10 個(gè)CC條目，21 個(gè)MF條目，結(jié)果表明山藥-黃芪藥對(duì)主要調(diào)節(jié)肽基絲氨酸磷酸化、調(diào)節(jié)細(xì)胞增殖、跨膜受體蛋白酪氨酸激酶信號(hào)通路、視黃酸的反應(yīng)、轉(zhuǎn)錄、ERK1和ERK2級(jí)聯(lián)的正調(diào)控、子宮內(nèi)胚胎發(fā)育、RNA聚合酶II啟動(dòng)子轉(zhuǎn)錄的正調(diào)控、炎癥反應(yīng)、凋亡過(guò)程的負(fù)調(diào)控等生物過(guò)程。通路富集分析表明，主要影響對(duì)糖尿病和骨質(zhì)疏松均有調(diào)節(jié)作用的卵巢類(lèi)固醇生成及雌激素信號(hào)通路、糖尿病相關(guān)信號(hào)通路(松弛素信號(hào)通路、黏著斑、催乳素信號(hào)通路、甲狀腺素信號(hào)通路、內(nèi)分泌抵抗、血小板活化[21]、FoxO信號(hào)通路[22])、與骨質(zhì)疏松成骨作用相關(guān)的作用通路(脂肪細(xì)胞中脂解的調(diào)節(jié)、VEGF信號(hào)通路、p53信號(hào)通路、IL-17信號(hào)通路)、與骨質(zhì)疏松破骨作用相關(guān)的信號(hào)通路(鈣信號(hào)通路、細(xì)胞衰老、花生四烯酸代謝)、與骨質(zhì)疏松骨成骨作用破骨作用均具有調(diào)節(jié)作用的PI3K-Akt信號(hào)通路[23-24]等。

這與文獻(xiàn)調(diào)研[25]不謀而合，雌激素可通過(guò)胰島β細(xì)胞中的雌激素受體促進(jìn)胰島細(xì)胞存活，脂質(zhì)穩(wěn)態(tài)和胰島素的生物合成和分泌；雌激素同時(shí)是骨骼代謝的主要激素調(diào)節(jié)劑，雌激素缺乏會(huì)增強(qiáng)骨吸收，導(dǎo)致骨重建失衡，還可促進(jìn)骨形成，提高骨密度，對(duì)糖尿病及骨質(zhì)疏松均具有較好的調(diào)節(jié)作用。糖尿病從多方面影響骨的形成以及骨吸收過(guò)程，誘發(fā)骨質(zhì)疏松，因此，干預(yù)糖尿病的發(fā)展可減緩骨質(zhì)疏松的發(fā)展。研究報(bào)道，松弛素信號(hào)通路[26]、黏著斑[27]、催乳素信號(hào)通路[28]、甲狀腺素信號(hào)通路[29]對(duì)糖尿病均具有調(diào)節(jié)作用。骨代謝的穩(wěn)態(tài)取決于成骨細(xì)胞誘導(dǎo)的骨形成和破骨細(xì)胞誘導(dǎo)的骨吸收之間的平衡。糖尿病性骨質(zhì)疏松主要是由于成骨細(xì)胞和破骨細(xì)胞之間的不平衡導(dǎo)致骨形成減少，骨密度和骨強(qiáng)度受損，以及骨微結(jié)構(gòu)的破壞。強(qiáng)化和恢復(fù)成骨細(xì)胞的成骨分化，重新平衡成骨細(xì)胞和破骨細(xì)胞的數(shù)量，可能有助于治療糖尿病性骨質(zhì)疏松癥，如雙歧曲霉多糖引起抗骨質(zhì)疏松作用的機(jī)制是由VEGF和OPG的上調(diào)介導(dǎo)的[30]；p53抑制成骨細(xì)胞增殖并增強(qiáng)成骨細(xì)胞凋亡，從而減少骨形成[31]；IL-17在骨質(zhì)流失中起重要作用，IL-17促進(jìn)了骨形成負(fù)調(diào)節(jié)劑硬化蛋白基因誘導(dǎo)的脂肪形成，并增強(qiáng)了對(duì)成骨的抑制作用[32]；NFATc1是破骨細(xì)胞分化的主要轉(zhuǎn)錄因子，可調(diào)控許多破骨細(xì)胞的特異功能基因，RANKL激活鈣信號(hào)通路，是激活NFATc1所必需的[33]；說(shuō)明VEGF信號(hào)通路、p53信號(hào)通路、IL-17信號(hào)通路等主要影響成骨作用，鈣信號(hào)通路等主要調(diào)節(jié)破骨細(xì)胞的功能。

綜上所述，本文采用網(wǎng)絡(luò)藥理學(xué)的分析方法探討山藥-黃芪藥對(duì)治療T2DOP的可能作用靶標(biāo)及途徑，并建立“山藥黃芪藥對(duì)化學(xué)成分-靶標(biāo)基因-信號(hào)通路”相互作用網(wǎng)絡(luò)圖，揭示了山藥黃芪藥對(duì)治療T2DOP的多維調(diào)控網(wǎng)絡(luò)，其中起主要作用的化學(xué)成分可能為山奈酚、異鼠李素、槲皮素、(6aR,11aR)-3-羥基-9,10-二甲氧基紫檀烷、7-O-甲基異戊醇、3,9-di-O-美迪紫檀素、海風(fēng)藤烯酮、芒柄花黃素、熊竹素、毛蕊異黃酮、異黃烷酮、蓽茇明寧堿、(3R)-3-(2-hydroxy-3,4-dimethoxyphenyl)chroman-7-ol 共13 個(gè)化合物，山藥黃芪藥對(duì)治療T2DOP的靶標(biāo)可能主要為CDK2、PTGS2、PTGS1、GSK3B、SRC、CYP19A1、EGFR、ESR1、ESR2、AR、MAPK14、ADORA2A、NCOA2、KDR、NOS2 共15 個(gè)靶標(biāo)，度值相對(duì)較高的信號(hào)通路為PI3K-Akt信號(hào)通路、鈣信號(hào)通路及雌激素信號(hào)通路，最終達(dá)到調(diào)節(jié)T2DOP的作用。