王猛 劉剛 路聊東 王東偉
摘要 目的:運(yùn)用網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接的方法探討“杜仲-牛膝”藥對(duì)治療骨質(zhì)疏松的分子機(jī)制。方法:通過(guò)中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫(kù)與分析平臺(tái)(TCMSP)篩選杜仲、牛膝含有的活性成分及作用靶點(diǎn);利用GeneCards、DisGeNET和CTD數(shù)據(jù)庫(kù)收集骨質(zhì)疏松(OP)相關(guān)靶點(diǎn)基因;采用Cytoscape 3.8.0軟件和STRING數(shù)據(jù)庫(kù)進(jìn)行“藥物-活性成分-靶點(diǎn)”網(wǎng)絡(luò)和蛋白質(zhì)-蛋白質(zhì)相互作用(PPI)網(wǎng)絡(luò)的構(gòu)建與分析;FunRich3.1.3軟件DAVID數(shù)據(jù)庫(kù)進(jìn)行基因本體(GO)富集分析和京都基因和基因組百科全書(shū)(KEGG) 富集分析(P<0.05);運(yùn)用Maestro version 11.5對(duì)主要活性成分及靶點(diǎn)蛋白進(jìn)行分子對(duì)接。結(jié)果:獲得“杜仲-牛膝”藥對(duì)治療OP的活性成分共42個(gè)和133個(gè)有效靶點(diǎn);獲得251條GO生物過(guò)程和90條相關(guān)信號(hào)通路,涉及PI3K-AKT、TNF、Wnt、HIF-1、MAPK信號(hào)通路等。結(jié)論:“杜仲-牛膝”藥對(duì)可能通過(guò)調(diào)控細(xì)胞增殖、分化與凋亡、骨代謝來(lái)治療OP,具有多成分、多靶點(diǎn)、多系統(tǒng)的特點(diǎn),為系統(tǒng)闡明“杜仲-牛膝”藥對(duì)主要成分治療OP的分子機(jī)制提供了參考。
關(guān)鍵詞 杜仲;牛膝;網(wǎng)絡(luò)藥理學(xué);分子對(duì)接;骨質(zhì)疏松;分子機(jī)制
Abstract Objective:To explore the molecular mechanism of Cortex Eucommiae-Achyranthes bidentata in the treatment of osteoporosis by network pharmacology and molecular docking.Methods:TCMSP database was used to screen the active components and action targets of Cortex Eucommiae-Achyranthes bidentata.We used GeneCards,DisGeNET and CTD database to collect osteoporosis(Osteoporosis,OP) related target genes; used Cytoscape 3.8.0 software and STRING database to construct and analyze the “drug-active ingredient-target” network and PPI network; FunRich3.1.3 Software DAVID database were used forGO analysis and KEGG enrichment analysis(P<0.05); we used Maestro version 11.5 to perform molecular docking of the main active ingredients and target proteins.Results:We obtained a total of 42 active ingredients and 133 effective targets of “Cortex Eucommiae-Achyranthes bidentata” drug for the treatment of OP; and obtained 251 GO biological processes and 90 related signal pathways,involving PI3K-AKT,TNF,Wnt,HIF-1,MAPK signal pathway etc.Conclusion:Cortex Eucommiae-Achyranthes bidentata drug pair may treat OP by regulating cell proliferation,differentiation and apoptosis,and bone metabolism.It has the characteristics of multi-component,multi-target and multi-system,which provides a reference for systematically elucidating the molecular mechanism of Cortex Eucommiae-Achyranthes bidentata drug pair in treating OP.
Keywords Cortex Eucommiae; Achyranthes bidentata; Network pharmacology; Molecular docking; Osteoporosis; Molecular mechanism
中圖分類號(hào):R285.6文獻(xiàn)標(biāo)識(shí)碼:Adoi:10.3969/j.issn.1673-7202.2021.20.002
骨質(zhì)疏松癥(Osteoporosis,OP)是最常見(jiàn)的骨骼疾病,是以全身性骨量減少,骨的微觀結(jié)構(gòu)退化為特征,表現(xiàn)為單位體積的骨量降低,礦鹽和骨基質(zhì)比例減少,生物脆性上升,導(dǎo)致骨的力學(xué)強(qiáng)度衰減,骨折風(fēng)險(xiǎn)上升的一種全身性骨代謝疾病[1]。2001年美國(guó)國(guó)立衛(wèi)生研究院(USA National Institutes of Health,NIH)將其定義為以骨折風(fēng)險(xiǎn)概率增加和骨的生物力學(xué)強(qiáng)度下降為特征的骨骼疾病,并提示骨量降低是骨質(zhì)疏松性骨折的主要危險(xiǎn)因素,但還存在其他危險(xiǎn)因素。我國(guó)OP的總患病率為6.6%~19.3%,平均為13%,其中約20%的患者會(huì)發(fā)生骨質(zhì)疏松性骨折[2]。許多臨床研究結(jié)果表明,臨床上常用的抗OP的西藥均具有一定的不良反應(yīng),如長(zhǎng)期使用激素治療的患者,其患骨壞死、冠心病、乳腺癌等疾病的概率將有不同程度的升高[3]。骨質(zhì)疏松的臨床癥狀與中醫(yī)疾病“骨痿”中的“不通則痛”“腰脊不舉”“兩足軟弱”“不任其身”等表現(xiàn)相符。中醫(yī)以“骨痿”論治骨質(zhì)疏松,重溫補(bǔ)腎陽(yáng)、強(qiáng)筋壯骨,兼以通經(jīng)活絡(luò)、祛風(fēng)除濕、補(bǔ)氣養(yǎng)血[4]。
杜仲味辛、平,性甘、溫,無(wú)毒,具有補(bǔ)肝腎,強(qiáng)筋骨,安胎等功效,《本草會(huì)言》曰:“腰膝之痛,非杜仲不除……補(bǔ)肝益腎,誠(chéng)為要藥?!迸Oノ犊唷⑺?,性甘平,主逐瘀通經(jīng),補(bǔ)肝腎,強(qiáng)筋骨,《醫(yī)學(xué)衷參西錄》認(rèn)為:“牛膝為補(bǔ)益之品而善引氣血下注,故善治腎虛腰痛,或腿痿不能任地?!監(jiān)P的主要病機(jī)在“腎虛”,腎為先天之本,人體的生長(zhǎng)壯老皆由腎氣的盛衰所主導(dǎo)[5]。而杜仲、牛膝補(bǔ)益肝腎、壯腰膝、強(qiáng)筋骨的功效顯著,其藥性溫和、不易產(chǎn)生不良反應(yīng),是治療OP的優(yōu)選中藥[6]。此兩藥最早配伍使用可追溯到《壽親養(yǎng)老新書(shū)》中“不老丸”一方?!岸胖?牛膝”藥對(duì)是治療OP的常用組方藥對(duì),但兩藥配對(duì)應(yīng)用后其多靶點(diǎn)作用機(jī)制尚不明確。
網(wǎng)絡(luò)藥理學(xué)是生物信息領(lǐng)域高速發(fā)展的一門新興學(xué)科,它集成了虛擬計(jì)算,高通量的基因組分析和網(wǎng)絡(luò)數(shù)據(jù)庫(kù)搜尋等方法,其基于系統(tǒng)生物學(xué)和藥物多靶標(biāo)多通路理論構(gòu)建“藥材-成分-靶標(biāo)”網(wǎng)絡(luò),分析通路機(jī)制,為藥物新方向預(yù)測(cè)和疾病的多向治療提供了導(dǎo)向[7]。分子對(duì)接是經(jīng)由受體的特性以及受體和藥物小分子之間的互相效用形式來(lái)施行藥物設(shè)計(jì)的門徑。主要探討分子間(如配體和受體)的相互作用,并預(yù)測(cè)其連接式樣和親和力的一種理論仿效方法[8]。本研究使用網(wǎng)絡(luò)藥理學(xué)及分子對(duì)接技術(shù),探討“杜仲-牛膝”藥對(duì)中有效成分與靶標(biāo)蛋白之間的互作關(guān)系,研究其藥對(duì)配伍后治療OP的可能作用機(jī)制,為后續(xù)的研究驗(yàn)證提供理論參考。
1 資料與方法
1.1 藥材中所含活性成分的篩選 在中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫(kù)與分析平臺(tái)(TCMSP)(http://tcmspw.com/tcmsp.php)以“杜仲”“牛膝”為關(guān)鍵詞,搜索杜仲、牛膝中含有的已知化學(xué)成分[9]。以藥物動(dòng)力學(xué)原理為基礎(chǔ),口服生物利用度(OB)是經(jīng)口服用的藥物中所含有的有效成分或活性物質(zhì)被人體吸收到達(dá)體循環(huán)的速度與程度,類藥性(DL)則是用以評(píng)價(jià)藥物分子成藥性高低的指標(biāo)[10]。結(jié)合文獻(xiàn)研究,將OB≥30%,DL≥0.18作為篩選標(biāo)準(zhǔn),篩選出杜仲、牛膝中的活性成分[11]。
1.2 藥物靶點(diǎn)的篩選及“藥物-活性成分-靶點(diǎn)”網(wǎng)絡(luò)的構(gòu)建 將1.1項(xiàng)下篩選出的活性成分導(dǎo)入TCMSP數(shù)據(jù)庫(kù)中進(jìn)行檢索,查找出藥物靶標(biāo)蛋白,通過(guò)Uniprot數(shù)據(jù)庫(kù)(https://www.uniprot.org/)進(jìn)行靶標(biāo)蛋白以及基因信息的校正。將藥物活性成分、靶點(diǎn)蛋白導(dǎo)入Cytoscape 3.8.0軟件,構(gòu)建“藥物-活性成分-靶點(diǎn)”網(wǎng)絡(luò)并進(jìn)行顯著性分析,篩選出關(guān)鍵活性成分。
1.3 “杜仲-牛膝”治療OP的潛在作用靶點(diǎn)獲取 ? 在GeneCards(https://www.genecards.org/)、DisGeNET(http://www.disgenet.org/)以及Comparative Toxicogenomics Database數(shù)據(jù)庫(kù)(http://ctdbase.org/)中以“Osteoporosis”或“OP”為關(guān)鍵詞,設(shè)定物種為“Homo Sapiens”檢索與OP相關(guān)的疾病靶點(diǎn),通過(guò)韋恩圖將“1.2”中篩選出的藥物靶點(diǎn)與相關(guān)疾病靶點(diǎn)映射篩選出共同靶點(diǎn),從而確定“杜仲-牛膝”藥對(duì)治療OP的潛在作用靶點(diǎn)。
1.4 蛋白質(zhì)-蛋白質(zhì)相互作用(PPI)網(wǎng)絡(luò)分析 將1.3項(xiàng)中篩選得到的治療OP潛在作用靶點(diǎn)導(dǎo)入STRING Version:11.0數(shù)據(jù)庫(kù)[12](https://string-db.org/),Organism(種屬)設(shè)為“Homo Sapiens(人類)”,得到PPI關(guān)系,保存數(shù)據(jù)為TSV格式。運(yùn)用Cytoscape3.8.0軟件構(gòu)建PPI網(wǎng)絡(luò)并分析篩選出核心靶點(diǎn)。
1.5 靶點(diǎn)富集分析與可視化 將“杜仲-牛膝”藥對(duì)治療OP的潛在作用靶點(diǎn)輸入到FunRich3.1.3軟件進(jìn)行基因本體(Gene Ontology,GO)富集分析,設(shè)定閾值為P<0.05,分別選取生物學(xué)過(guò)程(Biological Process,BP)、分子功能(Molecular Function,MF)、細(xì)胞組分(Cellular Component,CC)中顯著性前6的條目通過(guò)“Analysis”功能構(gòu)建圖形可視化。使用DAVID 6.8數(shù)據(jù)庫(kù)(https://david.ncifcrf.gov/)進(jìn)行京都基因和基因組百科全書(shū)(KEGG)富集分析[13],設(shè)置P<0.05,通路分析選取與OP關(guān)聯(lián)的且顯著性較高的信號(hào)通路借助Omicshare(http://omicshare.com/)進(jìn)行可視化。
1.6 分子對(duì)接 將1.2項(xiàng)下篩選出的主要活性成分輸入到ZINC(http://zinc.docking.org/)數(shù)據(jù)庫(kù)下載其3D結(jié)構(gòu),保存為“mol2”格式文件。將1.4項(xiàng)下篩選的核心蛋白靶點(diǎn)通過(guò)RCSB PDB數(shù)據(jù)庫(kù)(http://www.rcsb.org/)下載靶點(diǎn)蛋白的3D結(jié)構(gòu),保存為“pdb”格式文件,將靶點(diǎn)蛋白及活性成分的3D結(jié)構(gòu)導(dǎo)入到MAESTRO version 11.5中,靶點(diǎn)蛋白通過(guò)“Protein preparation”套件去水、加氫生成pH值為7.0±2的小分子構(gòu)象,活性成分通過(guò)“Ligprep”進(jìn)行預(yù)處理,“Receptor Grid Gener”套件確定活性位點(diǎn),運(yùn)用“Ligand Docking”面板進(jìn)行分子對(duì)接,對(duì)接方式為半柔性對(duì)接,以對(duì)接評(píng)分“Docking Score”評(píng)價(jià)“杜仲-牛膝”藥對(duì)主要活性成分與靶點(diǎn)蛋白的結(jié)合情況。
2 結(jié)果
2.1 藥物中主要化學(xué)成分的篩選 檢索得到“杜仲-牛膝”藥對(duì)中共290個(gè)化合物(5個(gè)為共有化合物),其中杜仲119個(gè),牛膝176個(gè)。篩選出42個(gè)活性化合物,其中25個(gè)來(lái)自杜仲,20個(gè)來(lái)自牛膝,3個(gè)為共有化合物。見(jiàn)圖1。“杜仲-牛膝”藥對(duì)中部分活性化合物基本信息見(jiàn)表1。
2.2 “藥物-活性成分-靶點(diǎn)”網(wǎng)絡(luò)的構(gòu)建與分析 ? 該網(wǎng)絡(luò)共有177個(gè)節(jié)點(diǎn),其中包含2個(gè)藥物節(jié)點(diǎn)、42個(gè)活性成分節(jié)點(diǎn)和133個(gè)靶點(diǎn)節(jié)點(diǎn),共有885條邊。其中藍(lán)色菱形代表杜仲、牛膝,黃色三角形代表活性成分,紅色箭頭代表靶點(diǎn),每條邊表示節(jié)點(diǎn)之間的相互作用關(guān)系。見(jiàn)圖2。借助“Network Analyze”功能對(duì)網(wǎng)絡(luò)進(jìn)行拓?fù)鋵傩苑治?。篩選度值(Degree)較大的節(jié)點(diǎn)進(jìn)行分析,這些節(jié)點(diǎn)在網(wǎng)絡(luò)中起到了樞紐作用,可能是關(guān)鍵化合物或靶點(diǎn)。該網(wǎng)絡(luò)中活性成分的平均度值為21.07。從活性成分方面分析,有51%的活性成分作用靶點(diǎn)≥21個(gè),其中槲皮素、β-谷甾醇、山柰酚、漢黃芩素、刺桐靈堿、豆甾醇分別能與20個(gè)以上的靶點(diǎn)蛋白發(fā)生作用,這些可能是“杜仲-牛膝”藥對(duì)中治療OP的主要物質(zhì)基礎(chǔ)。
2.3 “杜仲-牛膝”藥對(duì)治療OP的潛在作用靶點(diǎn) ? 基于GeneCards收集到26個(gè)OP靶點(diǎn),DisGeNET收集到63個(gè),CTD收集到820個(gè),三者合并去重后篩選共得到847個(gè)OP相關(guān)靶點(diǎn),將“杜仲-牛膝”藥對(duì)調(diào)控靶點(diǎn)和OP疾病靶點(diǎn)經(jīng)過(guò)Venny在線平臺(tái)(http://www.interactivenn.net/)分析后,獲得74個(gè)交集靶點(diǎn),即為藥物治療OP的潛在作用靶點(diǎn)。見(jiàn)圖3。
2.4 PPI網(wǎng)絡(luò)拓?fù)浞治?圖4顯示該網(wǎng)絡(luò)共有74個(gè)節(jié)點(diǎn)、856條邊,平均節(jié)點(diǎn)度值為23.14,圖中節(jié)點(diǎn)的大小、顏色和度值是正相關(guān)關(guān)系,節(jié)點(diǎn)的度值越大,顯示的節(jié)點(diǎn)越大且顏色越深。度值越大說(shuō)明與其他蛋白相互作用更強(qiáng)一些,可能在發(fā)揮生物學(xué)功能中起著重要的作用。度值排名前10的靶點(diǎn)為白細(xì)胞介素6(IL6)、腫瘤壞死蛋白P53(TP53)、血管內(nèi)皮生長(zhǎng)因子A(VEGFA)、腫瘤壞死因子(TNF)、JUN蛋白(JUN)、絲裂原活化蛋白激酶1(MAPK1)、表皮生長(zhǎng)因子(EGF)、過(guò)氧化氫合成酶2(PTGS2)、表皮生長(zhǎng)因子受體(EGFR)、絲裂原活化蛋白激酶8(MAPK8),這些靶點(diǎn)可能是“杜仲-牛膝”藥對(duì)發(fā)揮治療OP功效的關(guān)鍵靶點(diǎn)。
2.5 靶點(diǎn)富集分析 通過(guò)FunRich3.1.3軟件共篩選得到GO條目343個(gè)(P<0.05),其中BP條目251個(gè),MF條目64個(gè),CC條目28個(gè)。其中BP主要涉及序列特異性DNA結(jié)合轉(zhuǎn)錄因子活性的正調(diào)控、細(xì)胞增殖的負(fù)調(diào)控、凋亡過(guò)程的負(fù)調(diào)控、細(xì)胞增殖、凋亡過(guò)程、血管生成等;MF主要涉及蛋白質(zhì)結(jié)合、酶結(jié)合、蛋白質(zhì)均二聚活性、ATP結(jié)合、鋅離子結(jié)合等方面;CC方面主要涉及核、細(xì)胞質(zhì)、細(xì)胞外空間、胞質(zhì)溶膠、質(zhì)膜等方面。各類別前6的條目見(jiàn)圖5。
通過(guò)DAVID數(shù)據(jù)庫(kù)篩選得到90條(P<0.05)KEGG富集信號(hào)通路。主要涉及癌癥信號(hào)通路、鈣信號(hào)通路、PI3K-AKT、TNF、Wnt、cAMP、HIF-1等信號(hào)通路。富集顯著性前30的信號(hào)通路見(jiàn)圖6。縱軸代表通路名稱,橫軸代表富集因素的百分比;氣泡的大小和顏色分別代表富集靶點(diǎn)數(shù)和P值,顏色越紅P值越小。杜仲-牛膝”藥對(duì)和“PI3K-AKT信號(hào)通路的交集靶點(diǎn)見(jiàn)圖7。
2.6 分子對(duì)接 選取PPI網(wǎng)絡(luò)中度值排名前3的靶點(diǎn)蛋白,分別為IL6、TP53、VEGFA,與“藥物-活性成分-靶點(diǎn)”網(wǎng)絡(luò)中度值排名前3的活性成分,分別為槲皮素、β-谷甾醇、山柰酚,利用Maestro Version 11.5軟件進(jìn)行分子對(duì)接。見(jiàn)表2。Docking Score的絕對(duì)值越大,化合物與靶點(diǎn)結(jié)合力越強(qiáng),分子構(gòu)象越穩(wěn)定。根據(jù)對(duì)接分?jǐn)?shù)大小評(píng)判,絕對(duì)值>4.25代表結(jié)合活性一般,絕對(duì)值>5.0代表結(jié)合活性比較好,絕對(duì)值>7.0代表結(jié)合活性較強(qiáng)[14]。分子對(duì)接結(jié)果顯示IL6與槲皮素、山柰酚的親和力最強(qiáng),山柰酚與VEGFA的結(jié)合能力最強(qiáng)。圖8可以觀察到藥物活性成分與靶點(diǎn)蛋白相結(jié)合的具體位點(diǎn),黃色連線為二者之間的氫鍵相互作用力,為促使分子結(jié)合到活性位點(diǎn)的主要作用力。
3 討論
隨著人口老齡化進(jìn)程的加速,OP的發(fā)病率緊隨心血管病、糖尿病,已躍居至慢性疾病的第3位,成為世界的常見(jiàn)病和多發(fā)病[15]?,F(xiàn)代醫(yī)學(xué)對(duì)OP的治療主要為雌激素替代療法,但雌激素的長(zhǎng)期使用會(huì)增加乳腺癌、冠心病等疾病的患病概率,其他臨床上使用的化學(xué)合成藥物均存在一定的不良反應(yīng)[16]。而中藥治療OP具有不良反應(yīng)小、療效明確、經(jīng)濟(jì)負(fù)擔(dān)小的特點(diǎn)。中醫(yī)學(xué)認(rèn)為OP歸屬于“骨痿”“骨枯”“骨痹”等范疇。腎主骨生髓,腎精充盈,則骨髓生化有源,腎氣虧虛,則骨髓失養(yǎng)。由此,腎精虧損是該病的主要病機(jī)。而肝主筋,約束骨骼;脾主百骸,榮養(yǎng)四肢關(guān)節(jié);因此骨質(zhì)疏松的發(fā)病機(jī)制以腎虛為本,與肝虛、脾虛、血瘀密切相關(guān)[17]。杜仲、牛膝均可補(bǔ)肝腎,強(qiáng)筋骨。用于肝腎不足,腰膝酸痛,筋骨無(wú)力?,F(xiàn)代研究發(fā)現(xiàn),杜仲的化學(xué)成分主要有黃酮類、苯丙素類化合物以及多種氨基酸和微量元素,具有抗疲勞、抗骨質(zhì)疏松等多種藥理作用[18]。
根據(jù)網(wǎng)絡(luò)藥理學(xué)分析結(jié)果可以看出,槲皮素、山柰酚、β-谷甾醇、豆甾醇、黃芩素、刺桐靈堿、小檗堿是“杜仲-牛膝”藥對(duì)中主要的活性成分?;诜肿訉?duì)接模擬其結(jié)合活性,結(jié)果顯示槲皮素、山柰酚與IL6、TP53、VEGFA對(duì)接均具有較好的結(jié)合活性。槲皮素和山柰酚均為黃酮類化合物。鄭紅等[19]研究表明,槲皮素能夠有效促進(jìn)成骨分化,并通過(guò)上調(diào)BMP2和Smad4的表達(dá)改善OP。小鼠模型經(jīng)槲皮素處理后測(cè)定骨質(zhì),發(fā)現(xiàn)其骨鈣素含量上升,骨源性堿性磷酸酶(NBAP)及雌二醇的含量減少[20]。山柰酚能夠增加骨密度,改善骨微結(jié)構(gòu),實(shí)驗(yàn)研究表明Ca2+代謝可被高劑量的山柰酚所影響,通過(guò)調(diào)節(jié)代謝平衡來(lái)促進(jìn)構(gòu)造蛋白質(zhì)生成,改善骨小梁丟失以達(dá)到治療骨質(zhì)疏松的目的[21]。付方勝等[22]認(rèn)為黃芩素可以抑制破骨細(xì)胞的形成以及抑制破骨細(xì)胞的骨吸收功能,從而預(yù)防OP的發(fā)生。并且在一定程度上減緩骨流失并部分恢復(fù)骨量,同時(shí)在體外可以促進(jìn)骨髓間充質(zhì)干細(xì)胞(BMMSCs)的增殖和成骨分化[23]。小檗堿可以促進(jìn)p-AKT蛋白的表達(dá),并通過(guò)激活PI3K/AKT信號(hào)通路來(lái)促進(jìn)前成骨細(xì)胞的分化與礦化[24]。
本研究通過(guò)對(duì)“杜仲-牛膝”藥對(duì)治療OP的PPI網(wǎng)絡(luò)分析,發(fā)現(xiàn)IL6、TP53、VEGFA、TNF、JUN、MAPK1、EGF、PTGS2、EGFR、MAPK8等靶點(diǎn)是“杜仲-牛膝”藥對(duì)治療OP的主要作用靶點(diǎn)。其中IL6在成骨偶聯(lián)過(guò)程中發(fā)揮多重作用,成骨細(xì)胞與單核細(xì)胞均可合成IL6,IL6被合成后又直接作用于破骨細(xì)胞和成骨細(xì)胞,影響破骨細(xì)胞的活性,加速了OP的發(fā)生與發(fā)展[25]。VEGFA通過(guò)激活MAPK通路可以誘導(dǎo)BMSCs向成骨細(xì)胞分化,促進(jìn)了成骨。TNF屬于炎癥介質(zhì),它可以刺激破骨細(xì)胞分化形成,促進(jìn)分解代謝和骨細(xì)胞的破壞,是導(dǎo)致OP的重要環(huán)節(jié)[26]。TP53基因是一種腫瘤相關(guān)轉(zhuǎn)錄因子,在細(xì)胞應(yīng)激活動(dòng)中起重要作用。它通過(guò)激活多種細(xì)胞保護(hù)途徑,誘導(dǎo)細(xì)胞凋亡,調(diào)節(jié)細(xì)胞周期,參與成骨破骨過(guò)程[27]。EGF稱為表皮生長(zhǎng)因子,Vary等[28]研究發(fā)現(xiàn),成骨分化中的標(biāo)志基因血清堿性磷酸酶和Ocn的表達(dá)可被EGF抑制,說(shuō)明EGF可能影響成骨分化的負(fù)調(diào)控。JUN屬于JUN蛋白家族中的一種,而JUN蛋白是轉(zhuǎn)錄因子激活蛋白-1(Activatorprotein1,AP-1)家族中4個(gè)亞家族中的1個(gè),研究發(fā)現(xiàn)AP-1可通過(guò)調(diào)節(jié)相關(guān)基因表達(dá),來(lái)影響細(xì)胞的凋亡過(guò)程,促進(jìn)細(xì)胞凋亡,加快破骨前體細(xì)胞向破骨細(xì)胞分化[29-30]。MAPK1、MAPK8均屬于絲裂原活化蛋白激酶家族,其可通過(guò)對(duì)成骨細(xì)胞增殖和分化過(guò)程的促進(jìn),進(jìn)而誘導(dǎo)血管內(nèi)皮細(xì)胞的形成來(lái)發(fā)揮對(duì)OP的防治作用[31]。EGFR是ERK2的上游信號(hào)分子,具有酪氨酸激酶活性,可以通過(guò)與表皮生長(zhǎng)因子組合,啟動(dòng)細(xì)胞核內(nèi)基因,從而促進(jìn)細(xì)胞分裂[32]。
KEGG通路富集分析得到與OP較為相關(guān)的通路有PI3K-AKT、TNF、Wnt、HIF-1、MAPK信號(hào)通路。研究表明,PI3K-AKT信號(hào)通路參與骨質(zhì)疏松、骨性關(guān)節(jié)炎、骨壞死等病理性骨病,并參與調(diào)控破骨細(xì)胞和成骨細(xì)胞的增殖、分化及凋亡[33]。正常生理情況下,PI3K-AKT信號(hào)通路可以選擇性的影響成骨細(xì)胞和破骨細(xì)胞的生理功能,保護(hù)骨小梁結(jié)構(gòu),降低破骨細(xì)胞分化,通過(guò)降低Caspases-9活性,從而減輕細(xì)胞損傷所致的細(xì)胞凋亡,影響破骨細(xì)胞形成[34]。腫瘤壞死因子(TNF)信號(hào)通路具有調(diào)節(jié)破骨細(xì)胞活性及骨吸收的作用,其中分離出的糖蛋白骨保護(hù)素(Osteoprotegerin,OPG)對(duì)骨也具有積極的保護(hù)作用[35-36]。Wnt信號(hào)通路可以通過(guò)抑制MSCs向軟骨和成脂譜系的分化,來(lái)增強(qiáng)成骨細(xì)胞譜系分化[37-38]。成骨細(xì)胞和骨細(xì)胞的Wnt信號(hào)通過(guò)OPG的分泌從而間接抑制破骨細(xì)胞的分化和骨吸收[39]。Wnt信號(hào)通路直接或間接地對(duì)成骨細(xì)胞和破骨細(xì)胞進(jìn)行調(diào)控,在控制骨量與OP的發(fā)生發(fā)展方面具有重要作用。
綜上所述,本研究通過(guò)網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接相結(jié)合的方法,篩選出“杜仲-牛膝”藥對(duì)治療OP的活性成分,關(guān)鍵靶點(diǎn)和信號(hào)通路,發(fā)現(xiàn)“杜仲-牛膝”藥對(duì)可能通過(guò)調(diào)控細(xì)胞增殖、分化與凋亡、骨代謝來(lái)治療OP,具有多成分、多靶點(diǎn)、多系統(tǒng)的特點(diǎn),并基于分子對(duì)接初步模擬其可能的分子作用機(jī)制,為“杜仲-牛膝”藥對(duì)治療OP機(jī)制研究提供了理論依據(jù)和線索。鑒于網(wǎng)絡(luò)藥理學(xué)及分子對(duì)接技術(shù)的局限性,其預(yù)測(cè)結(jié)果還需要進(jìn)一步的研究驗(yàn)證。
參考文獻(xiàn)
[1]梁偉喬,鐘誠(chéng),李宇明.骨質(zhì)疏松癥的中醫(yī)病因病機(jī)認(rèn)識(shí)與治療進(jìn)展[J].中國(guó)骨質(zhì)疏松雜志,2020,26(1):135-139.
[2]白璧輝,謝興文,李鼎鵬,等.我國(guó)近5年來(lái)骨質(zhì)疏松癥流行病學(xué)研究現(xiàn)狀[J].中國(guó)骨質(zhì)疏松雜志,2018,24(2):253-258.
[3]單祎娜,王莉.骨質(zhì)疏松癥診斷和治療進(jìn)展[J].醫(yī)學(xué)綜述,2019,25(18):3652-3656,3661.
[4]向益,鄭烽,王顯.淺析骨質(zhì)疏松癥的中醫(yī)病名[J].內(nèi)蒙古中醫(yī)藥,2019,38(1):92.
[5]元唯安.中藥新藥研發(fā)關(guān)鍵問(wèn)題思考——以骨質(zhì)疏松癥為例[J].中國(guó)新藥雜志,2020,29(16):1836-1840.
[6]高衛(wèi)輝,吳芬芬,段小青,等.杜仲-牛膝藥對(duì)干預(yù)去卵巢骨質(zhì)疏松大鼠雌二醇和骨密度的影響實(shí)驗(yàn)研究[J].中南藥學(xué),2016,14(8):820-823.
[7]Zhang GB,Li QY,Chen QL,et al.Network pharmacology:a new approach for chinese herbal medicine research[J].Evid Based Complement Alternat Med,2013,2013:621423.
[8]Saikia S,Bordoloi M.Molecular docking:challenges,advances and its use in drug discovery perspective[J].Curr Drug Targets,2019,20(5):501-521.
[9]Hu AN,Liu JF,Zan JF,et al.Network pharmacology deciphering multiple mechanisms of volatiles of Wendan granule for treatment of senile dementia[J].Clin J Pharmacol Toxicol,2018,32(4):285-286.
[10]Xu X,Zhang C,Li P,et al.Drug-symptom networking:Linking drug-likeness screening to drug discovery[J].Pharmacol Res,2016,103:105-113.
[11]Wishart DS,Knox C,Guo AC,et al.DrugBank:a knowledgebase for drugs,drug actions and drug targets[J].Nucleic Acids Res,2008,36(Database issue):D901-906.
[12]von Mering C,Jensen LJ,Snel B,et al.STRING:known and predicted protein-protein associations,integrated and transferred across organisms[J].Nucleic Acids Res,2005,33(Database issue):D433-437.
[13]Dennis G Jr,Sherman BT,Hosack DA,et al.DAVID:database for annotation,visualization,and integrated discovery[J].Genome Biol,2003,4(5):P3.
[14]張英睿,干志強(qiáng),劉紫軒,等.基于網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接的清開(kāi)靈注射液治療新型冠狀病毒肺炎(COVID-19)的作用機(jī)制探尋[J].中草藥,2020,51(12):3201-3210.
[15]范勛奇,石磊.骨質(zhì)疏松癥相關(guān)研究進(jìn)展[J].廣東化工,2020,47(2):63-64,76.
[16]Zhao H,Zhao N,Zheng P,et al.Prevention and treatment of osteoporosis using chinese medicinal plants:special emphasis on mechanisms of immune modulation[J].J Immunol Res,2018,2018:6345857.
[17]張倩,王花欣,王媛,等.中藥防治絕經(jīng)后骨質(zhì)疏松癥的研究進(jìn)展[J].中國(guó)骨質(zhì)疏松雜志,2020,26(7):1083-1088.
[18]潘穎宜,邢蘊(yùn)蘊(yùn),王健英,等.杜仲骨保護(hù)的作用機(jī)制研究進(jìn)展[J].中華中醫(yī)藥學(xué)刊,2020,38(3):92-95.
[19]鄭紅,唐薇,角建林,等.槲皮素通過(guò)促進(jìn)成骨分化治療去勢(shì)骨質(zhì)疏松癥大鼠的分子機(jī)制[J].中藥藥理與臨床,2017,33(5):16-20.
[20]卞偉,孫宏,劉凱,等.槲皮素對(duì)骨質(zhì)疏松大鼠骨生物力學(xué)性能及骨代謝的影響[J].吉林中醫(yī)藥,2016,36(8):814-817.
[21]王雪峰,楊松,曹立新,等.山奈酚對(duì)卵巢去勢(shì)大鼠的促骨膠原生成和減少骨小梁丟失的保護(hù)作用研究[J].中國(guó)臨床藥理學(xué)雜志,2020,36(10):1302-1305.
[22]付方勝,丁佳昕,邵思遠(yuǎn),等.黃芩素抑制RANKL誘導(dǎo)的破骨細(xì)胞分化和功能[J].錦州醫(yī)科大學(xué)學(xué)報(bào),2019,40(2):18-20,118,119.
[23]桂林源.黃芩中小分子化合物對(duì)雌激素缺乏骨質(zhì)疏松的預(yù)防和治療作用的研究[D].西安:第四軍醫(yī)大學(xué),2016.
[24]江健,廖莉婭,容嬋,等.鹽酸小檗堿通過(guò)PI3K/Akt信號(hào)通路促進(jìn)小鼠前成骨細(xì)胞MC3T3-E1的分化[J].中國(guó)病理生理雜志,2018,34(10):1869-1875.
[25]許環(huán)順,候平,邢忠,等.鮭魚(yú)降鈣素對(duì)骨質(zhì)疏松癥患者血清胰島素樣生長(zhǎng)因子-1、白介素-6和腫瘤壞死因子-α的影響[J].疑難病雜志,2018,17(9):914-917,921.
[26]尹偉,楊斌輝,張波,等.骨質(zhì)疏松癥合并膝骨性關(guān)節(jié)炎患者骨密度及血清骨代謝指標(biāo)、炎癥因子水平研究[J].中國(guó)骨質(zhì)疏松雜志,2019,25(8):1121-1124.
[27]周旭艷,歐陽(yáng)桂芳,牧啟田,等.TP53突變骨髓增生異常綜合征患者的臨床特征及染色體核型分析[J].中國(guó)實(shí)驗(yàn)血液學(xué)雜志,2019,27(3):872-876.
[28]Vary CP,Li V,Raouf A,et al.Involvement of Ets transcription factors and targets in osteoblast differentiation and matrix mineralization[J].Exp Cell Res,2000,257(1):213-222.
[29]劉佳佳,李平,鄧琦,等.在Prrx1+骨髓間充質(zhì)干細(xì)胞中過(guò)表達(dá)肝素結(jié)合表皮生長(zhǎng)因子樣生長(zhǎng)因子抑制小鼠骨發(fā)育[J].中國(guó)組織工程研究,2019,23(29):4623-4628.
[30]劉夢(mèng)楠.CaMKⅡδ過(guò)表達(dá)對(duì)破骨細(xì)胞分化中c-fos、c-jun基因表達(dá)的影響[D].唐山:華北理工大學(xué),2019.
[31]陳鎮(zhèn)秋,洪郭駒,韓曉蕊等.骨科祛痰逐瘀方經(jīng)MAPK/ERK通路誘導(dǎo)血管內(nèi)皮細(xì)胞遷移與血管管腔形成[J].中國(guó)骨質(zhì)疏松雜志,2019,25(6):747-752.
[32]楊世春,溫本,黃小區(qū),等.探討ERα/EGFR在雙膦酸鹽對(duì)絕經(jīng)前后乳腺癌療效差異中的作用[J].中國(guó)處方藥,2019,17(8):142-143.
[33]Lin C,Shao Y,Zeng C,et al.Blocking PI3K/AKT signaling inhibits bone sclerosis in subchondral bone and attenuates post-traumatic osteoarthritis[J].J Cell Physiol,2018,233(8):6135-6147.
[34]強(qiáng)勝林,劉濤,賈永龍,等.補(bǔ)腎固本方對(duì)大鼠骨質(zhì)疏松模型中PI3K/AKT/mTOR信號(hào)通路表達(dá)的調(diào)控研究[J].中國(guó)骨質(zhì)疏松雜志,2020,26(7):983-987.
[35]田愛(ài)現(xiàn),馬劍雄,馬信龍.骨質(zhì)疏松關(guān)聯(lián)性分子信號(hào)通路及其臨床應(yīng)用[J].中華骨質(zhì)疏松和骨礦鹽疾病雜志,2019,12(3):282-291.
[36]Simonet WS,Lacey DL,Dunstan CR,et al.Osteoprotegerin:a novel secreted protein involved in the regulation of bone density[J].Cell,1997,89(2):309-319.
[37]Long H,Zhu Y,Lin Z,et al.miR-381 modulates human bone mesenchymal stromal cells(BMSCs) osteogenesis via suppressing Wnt signaling pathway during atrophic nonunion development[J].Cell Death Dis,2019,10(7):470.
[38]Yuan Z,Li Q,Luo S,et al.PPARγ and Wnt signaling in adipogenic and osteogenic differentiation of mesenchymal stem cells[J].Curr Stem Cell Res Ther,2016,11(3):216-225.
[39]王琰,王艷.Wnt信號(hào)通路在周圍神經(jīng)損傷后骨質(zhì)疏松中的研究進(jìn)展[J].中國(guó)康復(fù)理論與實(shí)踐,2020,26(8):930-935.
(2020-08-27收稿 責(zé)任編輯:魏慶雙,徐穎)