傅馨瑩，孫正驥，吳露，張偉

基于整合藥理學與時空效應的防治動脈粥樣硬化中藥篩選研究

傅馨瑩1，孫正驥1，吳露2，張偉1

1.湖南中醫(yī)藥大學，湖南 長沙 410208； 2.湖南中醫(yī)藥大學附屬瀏陽醫(yī)院，湖南 瀏陽 410300

整合多數(shù)據(jù)庫資源，結(jié)合動脈粥樣硬化（AS）疾病時空效應的特點，針對性篩選有效中藥，以更好地指導臨床合理用藥。通過TCMIP、GeneCards、OMIM數(shù)據(jù)庫獲取AS及時空效應表型的相關(guān)靶點；基于上述靶點，通過TCMIP數(shù)據(jù)庫“整合藥理學分析工具集合”的“反向查找中藥”功能，篩選血管內(nèi)皮功能障礙（VED）、血管內(nèi)皮損傷（VEI）、粥樣斑塊形成（PF）3個時間階段及骨髓微環(huán)境（BMM）空間過程的防治中藥；通過String、WebGestalt數(shù)據(jù)庫對AS時空效應防治中藥進行蛋白互作及GO、KEGG富集分析；通過R語言構(gòu)建弦圖與熱圖，對有效中藥-時空效應進行關(guān)聯(lián)性、差異性及最佳用藥配伍聚焦。收集AS-VED、AS-VEI、AS-PF、AS-BMM共同靶點分別為924、935、885、1024個；反向查找中藥結(jié)果顯示，高良姜、漏蘆、葛根、降香等可預防早期VED，牡荊葉、茜草、高良姜、漏蘆等可治療中期VEI，吳茱萸、急性子、葛根、降香等可消退晚期PF，皂角刺、香加皮、益母草、三棱、瓜蔞等在BMM中發(fā)揮作用，其作用機制可能與激活Rap1、Hippo、MAPK、Toll-like、TNF、AMPK、PI3K-AKT、NF-κB等多信號通路，調(diào)控AKT1、IL-6、INS、TNF、NR3C1、PPARG、PTGS2等多靶點，以“多藥物-多靶點-多通路”模式在防治AS進程中發(fā)揮作用有關(guān)。本研究篩選出針對特異時間階段與空間過程的中藥，實現(xiàn)了微觀層次上的辨證論治，為中醫(yī)藥防治AS提供了針對性的中藥及靶點和通路的可能，為中醫(yī)藥輔助干細胞療法提供思路。

整合藥理學；動脈粥樣硬化；時空效應；血管內(nèi)皮細胞；中醫(yī)藥

《中國心血管病報告2018》顯示，心血管病病死率仍處在首位，占居民疾病死亡構(gòu)成的43.16%～45.50%，且發(fā)病率持續(xù)上升，在2.9億心血管病患者中有高血壓2.45億、腦卒中1300萬、冠心病1100萬[1]。造成腦卒中與冠心病的前期主要病理基礎(chǔ)為動脈粥樣硬化（atherosclerosis，AS），血管病變作為AS的始動因素被研究者廣泛關(guān)注。血管由血管內(nèi)皮細胞（VECs）、平滑肌細胞（SMCs）、成纖維細胞構(gòu)成，血管病變主要是各種致病因素作用下VECs功能障礙和結(jié)構(gòu)損傷，局部產(chǎn)生前列環(huán)素、一氧化氮、內(nèi)皮素-1等多種血管活性物質(zhì)，導致血管內(nèi)源性修復機制失常而出現(xiàn)的血管病理性重構(gòu)[2]。

本課題組前期研究顯示，黃芪、當歸及其組分配伍可減輕血管壁炎癥，抑制血管病理性重構(gòu)，能較好地防治心腦血管疾病[3-6]。血管內(nèi)皮病變貫穿AS疾病進程始終，故提出內(nèi)皮細胞修復存在早期血管內(nèi)皮功能障礙（VED）、中期血管內(nèi)皮損傷（VEI）、后期粥樣斑塊形成（PF）時間階段與骨髓微環(huán)境（BMM）空間過程。時空效應進程存在多環(huán)節(jié)、多靶點、多途徑共同作用的特點，在一定程度上與整合藥理學思想不謀而合。因此，本研究在中醫(yī)整體觀念與辨證論治理論指導下，結(jié)合現(xiàn)代整合藥理學思想[7]，分析TCMIP[8]等多數(shù)據(jù)庫資源，根據(jù)AS不同時間階段與空間過程特點，以靶點為橋梁，針對性篩選有效中藥，以更好地指導臨床合理用藥。

1 資料與方法

1.1 動脈粥樣硬化疾病靶點與時空效應表型靶點收集

利用TCMIP（http://www.tcmip.cn/TCMIP/index. php/Home/）、GeneCards（https://www.genecards.org/）[9]及OMIM（https://omim.org/）數(shù)據(jù)庫，以“Atherogenesis”“Atherosclerosis”為檢索詞[10-11]，整合3個數(shù)據(jù)庫結(jié)果以獲得AS疾病靶點。

利用GeneCards數(shù)據(jù)庫，分別以“Vascular Endothelial Dysfunction”“Vascular Endothelial Injury”、“Plaque Formation”為檢索詞[12-14]，以獲得VED、VEI、PF的表型靶點。

利用GeneCards數(shù)據(jù)庫，以“Mobilization of Bone AND Homing AND Proliferation AND Migration”為檢索式，整合檢索結(jié)果以獲得BMM空間過程[15-18]的表型靶點。

1.2 動脈粥樣硬化-時空效應共同靶點整合分析

分別取AS疾病靶點與時空效應表型靶點交集，以獲得AS疾病中VED、VEI、PF階段及BMM過程的“疾病-表型”共同靶點，利用Bioinformatics Gent（http://bioinformatics.psb.ugent.be/beg）中Van de Peer Lab繪制韋恩圖。

1.3 動脈粥樣硬化-時空效應共同靶點反向預測有效中藥及差異篩選

利用TCMIP數(shù)據(jù)庫“整合藥理學分析工具集合”的“反向查找中藥”功能，設定子數(shù)據(jù)庫為“中藥”，分別將AS-VED、AS-VEI、AS-PF時間表型靶點與AS-BMM空間表型靶點導入“自定義基因集”，以獲得干預AS時空效應的中藥，根據(jù)OR值與值選取前20味中藥為AS時空效應有效中藥。

取AS-VED、AS-VEI、AS-PF 3個階段有效中藥交集，以篩選AS時間效應的差異中藥，并利用Van de Peer Lab繪制韋恩圖。

1.4 動脈粥樣硬化-時間效應差異中藥作用靶點收集及蛋白互作網(wǎng)絡構(gòu)建

利用TCMSP（http://lsp.nwu.edu.cn/tcmsp.php）數(shù)據(jù)庫[19]，以人體藥代動力學ADME[20]參數(shù)中口服生物利用度（OB）≥30%、類藥性（DL）≥0.18為條件，獲取AS-時間效應差異中藥有效成分。通過PubChem（https://pubchem.ncbi.nlm.nih.gov/）數(shù)據(jù)庫[21]收集各有效成分SDF結(jié)構(gòu)式。通過Swiss Target Prediction（http://www.swisstargetprediction.ch/）數(shù)據(jù)庫[22]，將有效成分SDF結(jié)構(gòu)式導入，獲取差異中藥作用靶點。

利用String（https://string-db.org/）數(shù)據(jù)庫[23]，分別導入AS時間效應差異中藥的作用靶點，選擇Multiple Protein，限定物種為“Homo sapiens”，余為默認值設定，以獲得差異中藥作用靶點蛋白互作（PPI）信息。同時采用Cytoscape3.7.1軟件對PPI信息進行Degree分析并進行網(wǎng)絡可視化。

1.5 動脈粥樣硬化-時間效應差異中藥的作用靶點GO和KEGG富集分析

利用WebGestalt（http://www.webgestalt.org/#）數(shù)據(jù)庫[24-25]，初步設定基本參數(shù)，物種為“Homo sapiens”，采用ORA（over-representation analysis）算法分析，選用Funtional Database子數(shù)據(jù)庫，進行GO、KEGG富集分析。進階參數(shù)設定中，設定“Significance Level”＜0.05，余為默認設置。收集生物過程（BP）、細胞組分（CC）、分子功能（MF）及KEGG富集分析數(shù)據(jù)，進行加權(quán)集合覆蓋（Weighted set cover）分析，以富集的優(yōu)勢比（OR值）的log2對數(shù)為X軸，以偽發(fā)現(xiàn)率（FDR值）的負對數(shù)（-lg）為Y軸，繪制火山圖，OR值越大表明富集越顯著，F(xiàn)DR值越小表明富集錯誤發(fā)生率越低。

1.6 有效中藥-時空效應的關(guān)聯(lián)性與差異性分析及最佳用藥配伍聚焦

去除80味VED、VEI、PF階段及BMM過程有效中藥的重復值，以整合時空效應的相同/差異中藥。利用R語言Circlize包與Pheatmap包繪制有效中藥-時空效應弦圖與熱圖。從關(guān)聯(lián)性和差異性角度整合分析有效中藥-時空效應的作用關(guān)系，聚焦AS時空效應中各時間階段與空間過程的最佳用藥配伍。

2 結(jié)果

2.1 動脈粥樣硬化疾病靶點與時空效應表型靶點

通過TCMIP、GeneCards及OMIM數(shù)據(jù)庫，收集AS疾病靶點共1165個。通過GeneCards數(shù)據(jù)庫收集得到VED、VEI、PF、BMM時空效應表型靶點，分別為7002、5679、4024、17 606個。

2.2 動脈粥樣硬化-時空效應共同靶點整合分析

利用Bioinformatics Gent取交集分別獲得AS-VED、AS-VEI、AS-PF、AS-BMM共同靶點924、935、885、1024個，韋恩圖見圖1。

2.3 動脈粥樣硬化-時空效應共同靶點有效中藥

根據(jù)上述AS-時空效應共同靶點，利用TCMIP數(shù)據(jù)庫“反向查找中藥”功能，獲得AS時空效應80味中藥，前20味有效中藥見表1。

利用Van de Peer Lab獲取AS時間效應的差異中藥，共有中藥16味，AS-VED與AS-VEI差異中藥2味，AS-VED與AS-PF差異中藥2味，AS-VEI差異中藥2味，AS-PF差異中藥2味，韋恩圖見圖2。

圖1 AS-時空效應共同靶點韋恩圖

表1 AS-時空效應前20味有效中藥

序號AS-VED AS-VEI AS-PF AS-BMM 中藥OR值P值 中藥OR值P值 中藥OR值P值 中藥OR值P值 1皂角刺2.588.65E-14 香加皮3.231.58E-13 皂角刺2.618.65E-14 皂角刺3.408.79E-14 2香加皮3.149.11E-13 皂角刺2.552.80E-13 香加皮3.229.11E-13 香加皮2.542.90E-13 3益母草3.261.20E-11 瓜蔞2.142.09E-11 瓜蔞2.201.20E-11 益母草3.534.96E-12 4三棱2.501.32E-11 三棱2.464.08E-11 三棱2.501.32E-11 三棱2.212.96E-11 5瓜蔞2.107.52E-11 益母草3.185.47E-11 娑羅子5.357.52E-11 瓜蔞2.503.70E-11 6娑羅子5.111.01E-10 娑羅子5.139.51E-11 益母草3.201.01E-10 娑羅子5.001.58E-10 7黑豆2.237.12E-10 珠子參2.207.28E-10 藍布正2.007.12E-10 黑豆2.422.51E-10 8珠子參2.208.44E-10 補骨脂2.092.21E-9 黑豆2.248.44E-10 珠子參2.436.34E-10 9藍布正1.962.19E-9 苦玄參4.583.63E-9 馬兜鈴2.252.19E-9 藍布正2.511.79E-9 10竹根七2.282.23E-9 藍布正1.944.93E-9 苦玄參4.112.23E-9 竹根七2.142.56E-9 11補骨脂2.092.57E-9 黑豆2.175.65E-9 補骨脂2.172.57E-9 補骨脂1.862.63E-9 12苦玄參4.573.85E-9 竹根七2.236.09E-9 珠子參2.173.85E-9 苦玄參3.063.52E-9 13高良姜2.814.22E-9 馬兜鈴2.176.94E-9 三七2.114.22E-9 高良姜2.174.50E-9 14三七2.094.68E-9 三七2.071.15E-8 竹根七2.274.68E-9 三七2.184.61E-9 15葶藶子2.114.90E-9 葶藶子2.081.21E-8 葶藶子2.114.90E-9 葶藶子3.045.68E-9 16降香2.777.26E-9 漏蘆3.001.30E-8 穿山龍3.007.26E-9 降香4.935.78E-9 17馬兜鈴2.157.94E-9 穿山龍2.941.31E-8 急性子2.497.94E-9 馬兜鈴3.326.01E-9 18葛根2.001.20E-8 茜草2.161.48E-8 吳茱萸2.521.20E-8 葛根2.326.86E-9 19漏蘆3.001.41E-8 高良姜2.761.55E-8 葛根2.021.41E-8 漏蘆3.508.20E-9 20穿山龍2.921.43E-8 牡荊葉4.951.96E-8 降香2.791.43E-8 穿山龍1.941.24E-8

圖2 AS-時間效應異同中藥韋恩圖

2.4 動脈粥樣硬化-時間效應差異中藥的作用靶點蛋白互作網(wǎng)絡

利用String數(shù)據(jù)庫及Cytoscape3.7.1軟件對AS-時間效應差異中藥靶點進行蛋白共表達網(wǎng)絡分析，得出PPI網(wǎng)絡圖。設定Degree 1～19為粉色，20～29為藍色，≥30為紅色。

AS-VED與AS-VEI差異中藥作用靶點PPI網(wǎng)絡圖（圖3A）有48個節(jié)點，260條邊，＜1.0e-16，平均節(jié)點度10.8，聚類系數(shù)0.584，其中Degree≥20的靶點為AKT1（Degree＝30）、ALB（Degree＝24）、TNF（Degree＝21）、PPARG（Degree＝21）、NR3C1（Degree＝20）。

AS-VED與AS-PF差異中藥作用靶點PPI網(wǎng)絡圖（圖3B）有71個節(jié)點，503條邊，＜1.0e-16，平均節(jié)點度14.2，聚類系數(shù)0.632，其中Degree≥20的靶點為INS（Degree＝42）、AKT1（Degree＝38）、IL6（Degree＝38）、PPARG（Degree＝25）、NR3C1（Degree＝24）、PTGS2（Degree＝23）、ESR1（Degree＝22）、GASP3（Degree＝21）、IL1B（Degree＝21）、AR（Degree＝20）、TLR4（Degree＝20）。

AS-VEI差異中藥作用靶點PPI網(wǎng)絡圖（圖3C）有52個節(jié)點，253條邊，＜1.0e-16，平均節(jié)點度9.73，聚類系數(shù)0.521，其中Degree≥20的靶點為AKT1（Degree＝33）、TNF（Degree＝27）、PTGS2（Degree＝22）、NR3C1（Degree＝21）、CASP3（Degree＝20）、ESR1（Degree＝20）。

AS-PF差異中藥作用靶點PPI網(wǎng)絡圖（圖3D）有54個節(jié)點，298條邊，＜1.0e-16，平均節(jié)點度11，聚類系數(shù)0.618，其中Degree≥20的靶點為INS（Degree＝37）、IL6（Degree＝31）、TNF（Degree＝29）、PPARG（Degree＝25）、TP53（Degree＝25）、PTGS2（Degree＝24）、IL1B（Degree＝23）、NR3C1（Degree＝21）、SREBF1（Degree＝21）、TLR4（Degree＝20）。

注：A. AS-VED與AS-VEI差異中藥作用靶點；B. AS-VED與AS-PF差異中藥作用靶點；C. AS-VEI差異中藥作用靶點；D. AS-PF差異中藥作用靶點

2.5 動脈粥樣硬化-時間效應差異中藥的作用靶點GO和KEGG富集分析

AS-VED與AS-VEI差異中藥作用靶點的BP富集分析了487個靶點，與198個生物過程相關(guān)，主要富集于肽絲氨酸修飾、肽基酪氨酸修飾、白細胞分化等過程。CC富集分析了487個靶點，與16個細胞組分相關(guān)，主要富集于神經(jīng)元-神經(jīng)元突觸、神經(jīng)元細胞體，突觸后修飾等組分。MF富集分析了95個靶點，與10個分子功能相關(guān)，主要富集于對生物的刺激反應、防御反應等功能。KEGG富集分析了487個靶點，與94條信號通路相關(guān)，主要富集于Rap1、Hippo、MAPK、Toll-like、cGMP-PKG等信號通路，結(jié)果見圖4A。

AS-VED與AS-PF差異中藥作用靶點的BP富集分析了71個靶點，與286個生物過程相關(guān)，分析主要富集于調(diào)節(jié)炎癥反應、脂質(zhì)定位、類固醇代謝等過程。CC富集分析了71個靶點，與3個細胞組分相關(guān)，主要富集于受體復合體、轉(zhuǎn)錄因子復合體、膜區(qū)等組分。MF富集分析了71個靶點，與32個分子功能相關(guān)，主要富集于轉(zhuǎn)錄因子結(jié)合DNA介導轉(zhuǎn)錄、類固醇結(jié)合、磷酸酶結(jié)合等功能。KEGG富集分析了71個靶點，與115條信號通路相關(guān)，主要富集于胰島素抵抗及TNF、AMPK、PI3K-AKT等信號通路，結(jié)果見圖4B。

AS-VEI差異中藥作用靶點的BP富集分析了52個靶點，與195個生物過程相關(guān)，主要富集于類固醇激素反應、細菌起源因子反應、調(diào)節(jié)炎癥反應等過程。CC富集分析了52個靶點，與2個細胞組分相關(guān)，主要富集于轉(zhuǎn)錄因子復合受體、受體復合物等組分。MF富集分析了52個靶點，與17個分子功能相關(guān)，主要富集于類固醇激素受體活性、轉(zhuǎn)錄因子結(jié)合DNA介導轉(zhuǎn)錄、類固醇結(jié)合等功能。KEGG富集分析了52個靶點，與75條信號通路相關(guān)，主要富集于胰島素抵抗、IL-17、MAPK、Toll-like、cGMP-PKG、自然殺傷細胞介導的免疫反應等信號通路，結(jié)果見圖4C。

AS-PF差異中藥作用靶點的BP富集分析了54個靶點，與235個生物過程相關(guān)，主要富集于調(diào)節(jié)免疫反應、調(diào)節(jié)脂質(zhì)代謝、脂質(zhì)定位等過程。CC富集分析了54個靶點，與3個細胞組分相關(guān)，主要富集于酶復合物-含磷基團轉(zhuǎn)移、受體復合體、膜區(qū)等組分。MF富集分析了54個靶點，與19個分子功能相關(guān)，主要富集于轉(zhuǎn)錄因子結(jié)合DNA介導轉(zhuǎn)錄、類固醇激素受體活性、類固醇結(jié)合等功能。KEGG富集分析了54個靶點，與106條信號通路相關(guān)，主要富集于胰島素抵抗、NF-κB、AMPK、PI3K-AKT、Th17細胞分化等信號通路，結(jié)果見圖4D。

注：A. AS-VED與AS-VEI差異中藥作用靶點；B. AS-VED與AS-PF差異中藥作用靶點；C. AS-VEI差異中藥作用靶點；D. AS-PF差異中藥作用靶點

2.6 有效中藥-時空效應的關(guān)聯(lián)性與差異性分析及最佳用藥配伍聚焦

TCMIP數(shù)據(jù)庫對AS時空效應反向預測80個有效中藥，主要聚焦于香加皮、皂角刺、益母草、瓜蔞、三棱、娑羅子、珠子參、黑豆、竹根七、補骨脂、藍布正、降香、三七、葶藶子、漏蘆、苦玄參、穿山龍、馬兜鈴、高良姜、葛根、急性子、吳茱萸、茜草及牡荊葉。將上述24味有效中藥與AS時空效應構(gòu)建“有效中藥-時空效應”弦圖（見圖5）與熱圖（見圖6）。

注：線條寬度默認值為1，中藥線條寬度默認值為4

注：黃色為差異中藥，藍色為共有中藥

3 討論

整合藥理學是基于網(wǎng)絡大數(shù)據(jù)，整合中藥學、藥理學、藥代動力學、計算機學等多學科知識，主要采用多維度、多層次的研究方法，具有復雜的網(wǎng)絡特性和整體、動態(tài)的特點[26-27]，實現(xiàn)了“中藥成分-體內(nèi)代謝-分子網(wǎng)絡-藥理作用-中醫(yī)理論”多層次可視化的網(wǎng)絡分析，對中醫(yī)藥與人體相互作用規(guī)律的研究具有重要作用。中藥復方在人體產(chǎn)生作用的機制及代謝過程，體現(xiàn)了復雜的動態(tài)網(wǎng)絡特性。因此，基于整合藥理學研究中醫(yī)理論，可揭示中藥復方多成分、多靶點、多途徑作用于人體的整體效應。

VECs是依附于血管壁緊密排列的單層扁平細胞，具有屏障、分泌和血管生成的功能，在調(diào)節(jié)血管功能中發(fā)揮重要作用[28-29]。VECs功能障礙與損傷是AS疾病的始動因素與病理基礎(chǔ)[30]。在異常血流剪切力、高血脂、氧自由基等因素作用下，內(nèi)皮細胞周圍微環(huán)境改變，分泌功能紊亂，一氧化氮、內(nèi)皮素-1、白細胞介素（IL）等血管活性物質(zhì)分泌異常[31]，VECs出現(xiàn)衰老、凋亡、剝脫的結(jié)構(gòu)損傷，使緊密和黏附連接結(jié)構(gòu)為基礎(chǔ)的屏障功能減弱或喪失，血管壁通透性增加，內(nèi)膜完整性受到破壞，從而使炎性因子、血小板黏附因子、單核細胞及脂質(zhì)等向VECs下遷移與沉積，最終導致PF。因此，VED、VEI、PF是AS病理過程中的3個時間階段。生理條件下，由內(nèi)皮祖細胞（EPCs）介導的修復代償機制與VECs介導的病理損傷過程存在著動態(tài)平衡[32-33]，EPCs存在于骨髓、胚胎細胞臍帶血中，VECs損傷后，EPCs經(jīng)歷動員（募集到外周血）、歸巢（黏附整合到損傷部位）、增殖及分化（增殖、分化為VECs）等環(huán)節(jié)，從空間層次上對AS病理過程進行自我修復[34]。因此，尋找特異性修復3個時間階段及促進EPCs修復血管的空間過程的藥物，可達到更高效修復血管病變的目的，為中醫(yī)藥輔助干細胞療法提供思路。

本研究通過TCMIP數(shù)據(jù)庫，以疾病和表型相關(guān)靶點為橋梁，反向預測具有特異修復AS時空效應的中藥共80味，其中時間效應60味，整合后為24味，空間效應20味。結(jié)果顯示高良姜、漏蘆、葛根、降香等中藥可預防早期VED，牡荊葉、茜草、高良姜、漏蘆等中藥可治療中期VEI，吳茱萸、急性子、葛根、降香等可消退晚期PF，皂角刺、香加皮、益母草、三棱、瓜蔞等在BMM中發(fā)揮作用。研究表明，高良姜的主要有效成分高良姜素能通過調(diào)控血管內(nèi)皮細胞生長因子（VEGF），在血管新生中起作用；也能激活NF-κB/ERK通路調(diào)控TNF-α、IL-1β、IL-6因子，介導血管氧化應激和炎癥反應[35-36]。葛根的主要有效成分葛根素與梓醇配伍可通過激活ERK/PI3K-AKT/mTOR信號通路，上調(diào)HIF-1α，抗VECs凋亡，促進血管新生及修復和EPCs歸巢[37]。茜草炭可增強VECs活性，降低ox-LDL對VECs的損傷[38]。吳茱萸有效成分吳茱萸次堿通過激活辣椒素受體促進降鈣素基因相關(guān)肽（CGRP）的合成與釋放，與VECs中的CGRP受體結(jié)合以促進VECs增殖，抑制其凋亡，修復VEI[39]。因此，在AS時空效應轉(zhuǎn)化進程中，針對各時間階段及空間過程合理選用中藥，可更大程度改善各病理階段的血管病變，促進EPCs遷移過程，更好地防治AS。

本研究基于時空效應差異中藥，進一步探討中醫(yī)藥防治AS的可能作用機制，PPI網(wǎng)絡及GO、KEGG富集分析結(jié)果表明，AS時空效應過程受到多靶點相互作用、多信號通路共同介導，具有多靶點-多通路特點，但是存在偏重效應。其中，多靶點主要涉及AKT、TNF、ALB、NR3C1、IL1B、GASP3等，多通路主要與TNF、PI3K-AKT、VEGF信號通路的多步驟連續(xù)調(diào)控有關(guān)，在時空效應層面上具有偏重效應，如AKT1在VED-VEI轉(zhuǎn)化階段Degree為30，在VED-PF轉(zhuǎn)化階段為38，在VEI階段為33，在PF階段未見AKT1分子相互作用，因此AKT1可能主要參與VEI、VED階段，在轉(zhuǎn)化階段可能起到重要調(diào)控作用，在PF階段作用可能不顯著。此預測結(jié)果可作為中醫(yī)藥防治AS的新思路，為后續(xù)研究提供參考。

AKT分子是PI3K-AKT信號通路中重要的下游靶激酶，磷酸化的AKT通過介導SMCs遷移、血小板活化、白細胞黏附等，促進AS的發(fā)生發(fā)展。有研究顯示，在AS疾病進程中，氧化應激、炎癥因子及ox-LDL可導致PI3K和磷酸化的AKT失活，通過上調(diào)凋亡基因Bax、Bak表達，使線粒體膜上形成通道并與之結(jié)合，促進細胞色素穿過線粒體膜進入胞質(zhì)，導致VECs凋亡、脫落，增加血管內(nèi)膜通透性[40-41]。TNF-α可直接損傷血管內(nèi)膜，促進VECs表面黏附分子異常表達，促進穿孔蛋白在VECs表面形成膜間通道，使B粒酶進入胞內(nèi)，激活Caspases相關(guān)蛋白活性，破壞VECs結(jié)構(gòu)及DNA斷裂，引起VECs死亡、脫落，血管內(nèi)膜通透性增加，進而促進白細胞、血小板黏附于內(nèi)皮表面，SMCs增殖，脂質(zhì)向內(nèi)膜下沉積，巨噬細胞吞噬脂質(zhì)分化成泡沫細胞，內(nèi)膜增厚，斑塊形成，從而導致AS的發(fā)生發(fā)展。研究表明，EPCs能對AS損傷VECs進行自我修復，主要通過激活VEGF信號通路，促進EPCs動員、歸巢、增殖、分化為VECs。VEGF是骨髓動員的重要因子，VEGF可上調(diào)SDF-1與CXCR4，激活PI3K和AKT，促進eNOS磷酸化，釋放一氧化氮，同時VEGF與EPCs表面的VEGFR2受體結(jié)合，促進EPCs動員進入外周血[42]。內(nèi)膜受損部位VECs釋放VEGF等趨化因子與進入外周血的EPCs結(jié)合，促進EPCs歸巢，增殖分化為VECs參與血管內(nèi)膜修復。因此，多靶點、多通路貫穿AS時空效應始終，綜合發(fā)揮調(diào)控作用，為靶向防治策略提供了理論依據(jù)，為中醫(yī)藥多靶點防治提供了可能。

本研究進一步對有效中藥-時空效應進行關(guān)聯(lián)與差異性聚焦分析，得出24味中藥與時空效應各階段的對應關(guān)系，同時整合中藥作用靶點PPI網(wǎng)絡和KEGG富集分析可得出高良姜、葛根、茜草、吳茱萸等各時空效應相關(guān)中藥可能是通過激活Rap1、Hippo、MAPK、Toll-like、TNF、AMPK、PI3K-AKT、NF-κB等多信號通路，調(diào)控AKT1、IL-6、INS、TNF、NR3C1、PPARG、PTGS2等多靶點，以“多藥物-多靶點-多通路”模式在防治AS進程中發(fā)揮重要作用。

綜上所述，本研究基于整合藥理學分析AS時間效應3個階段與空間效應中BMM的靶點，通過TCMIP數(shù)據(jù)庫反向預測修復3個階段及促進EPCs修復血管的空間過程的藥物，與益氣活血中藥黃芪-當歸配伍，根據(jù)疾病不同階段特點實現(xiàn)辨證論治，為中醫(yī)藥防治AS時空效應的靶點和通路提供依據(jù)。

[1] 胡盛壽,高潤霖,劉力生,等.《中國心血管病報告2018》概要[J].中國循環(huán)雜志,2019,34(3)：209-220.

[2] WONG B W, MARSCH E, TREPS L, et al. Endothelial cell metabolism in health and disease：impact of hypoxia[J]. EMBO J,2017,36(15)：2187-2203.

[3] 張偉,鄧常清.從氣血理論探討病理性血管重構(gòu)的研究思路[J].中華中醫(yī)藥雜志,2019,34(10)：4713-4717.

[4] 彭熙煒,閻卉芳,黃娟,等.黃芪-當歸不同配伍比例對大鼠血管內(nèi)膜增生的影響[J].中國中醫(yī)藥信息雜志,2017,24(6)：56-59.

[5] 閻卉芳,徐昊,彭熙煒,等.黃芪和當歸配伍對大鼠血管內(nèi)膜增生血管平滑肌細胞增殖的影響[J].中國臨床藥理學與治療學,2018,23(4)：361-369.

[6] YAN H, PENG X, XU H, et al. Inhibition of aortic intimal hyperplasia and vascular smooth muscle proliferation and extracellular matrix protein expressions by Astragalus-Angelica combination[J]. Evid Based Complement Alternat Med,2018,2018：1508637.

[7] 許海玉,楊洪軍.整合藥理學：中藥現(xiàn)代研究新模式[J].中國中藥雜志,2014,39(3)：357-362.

[8] XU H Y, ZHANG Y Q, LIU Z M, et al. ETCM：an encyclopaedia of traditional Chinese medicine[J]. Nucleic Acids Res,2019,47(D1)：D976-D982.

[9] STELZER G, ROSEN N, PLASCHKES I, et al. The GeneCards suite：From gene data mining to disease genome sequence analyses[J]. Curr Protoc Bioinformatics,2016,54(6)：1.30.1-1.30.33.

[10] BROWN R A, SHANTSILA E, VARMA C, et al. Current understanding of atherogenesis[J]. Am J Med,2017,130(3)：268-282.

[11] EMINI VESELI B, PERROTTA P, De Meyer G, et al. Animal models of atherosclerosis[J]. Eur J Pharmacol,2017,816(12)：3-13.

[12] MADKHALI H A. Morin attenuates high-fat diet induced-obesity related vascular endothelial dysfunction in Wistar albino rats[J]. Saudi Pharm J,2020,28(3)：300-307.

[13] LIN F, YANG Y, WEI S, et al. Hydrogen sulfide protects against high glucose-induced human umbilical vein endothelial cell injury through activating PI3K/Akt/eNOS pathway[J]. Drug Des Devel Ther, 2020,14：621-633.

[14] ZHAO J, MI X, ZHAO L, et al. Validation of PM2.5 model particle through physicochemical evaluation and atherosclerotic plaque formation in ApoE-/-mice[J]. Ecotoxicol Environ Saf,2020,192(4)：110308.

[15] LENKIEWICZ A M, ADAMIAK M, THAPA A, et al. The Nlrp3 inflammasome orchestrates mobilization of bone marrow-residing stem cells into peripheral blood[J]. Stem Cell Rev Rep,2019, 15(3)：391-403.

[16] CHAVAKIS E, DIMMELER S. Homing of progenitor cells to ischemic tissues[J]. Antioxid Redox Signal,2011,15(4)：967-980.

[17] DEL PAPA N, PIGNATARO F. The role of endothelial progenitors in the repair of vascular damage in systemic sclerosis[J]. Front Immunol,2018,9(6)：1383.

[18] POTENTE M, CARMELIET P. The link between angiogenesis and endothelial metabolism[J]. Annu Rev Physiol,2017,79(2)：43-66.

[19] RU J, LI P, WANG J, et al. TCMSP：a database of systems pharmacology for drug discovery from herbal medicines[J]. J Cheminform,2014,6(4)：13.

[20] ZHANG W, CHEN Y, JIANG H, et al. Integrated strategy for accurately screening biomarkers based on metabolomics coupled with network pharmacology[J]. Talanta,2020,211(5)：120710.

[21] KIM S, CHEN J, CHENG T, et al. PubChem 2019 update：improved access to chemical data[J]. Nucleic Acids Res,2019,47(D1)：D1102- D1109.

[22] GFELLER D, MICHIELIN O, ZOETE V. Shaping the interaction landscape of bioactive molecules[J]. Bioinformatics,2013, 29(23)：3073-3079.

[23] SZKLARCZYK D, GABLE A L, LYON D, et al. STRING v11：protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets[J]. Nucleic Acids Res,2019,47(D1)：D607-D613.

[24] WANG J, VASAIKAR S, SHI Z, et al. WebGestalt 2017：a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit[J]. Nucleic Acids Res,2017,45(W1)：W130-W137.

[25] LIAO Y, WANG J, JAEHNIG E J, et al. WebGestalt 2019：gene set analysis toolkit with revamped UIs and APIs[J]. Nucleic Acids Res, 2019,47(W1)：W199-W205.

[26] 許海玉,劉振明,付巖,等.中藥整合藥理學計算平臺的開發(fā)與應用[J].中國中藥雜志,2017,42(18)：3633-3638.

[27] 劉志華,孫曉波.網(wǎng)絡藥理學：中醫(yī)藥現(xiàn)代化的新機遇[J].藥學學報, 2012,47(6)：696-703.

[28] CHIU J J, CHIEN S. Effects of disturbed flow on vascular endothelium：pathophysiological basis and clinical perspectives[J]. Physiol Rev,2011,91(1)：327-387.

[29] DAIBER A, STEVEN S, WEBER A, et al. Targeting vascular (endothelial) dysfunction[J]. Br J Pharmacol,2017,174(12)：1591- 1619.

[30] 胡月,李潞,趙紅麗.血管內(nèi)皮功能與冠狀動脈粥樣硬化發(fā)生發(fā)展的研究進展[J].中西醫(yī)結(jié)合心腦血管病雜志,2018,16(11)：1525-1528.

[31] 饒璇,李元建.內(nèi)皮細胞損傷與修復的研究進展[J].中國動脈硬化雜志,2017,25(5)：531-535.

[32] MILUTINOVI? A, ?UPUT D, ZORC-PLESKOVI? R. Pathogenesis of atherosclerosis in the tunica intima, media, and adventitia of coronary arteries：an updated review[J]. Bosn J Basic Med Sci, 2020,20(1)：21-30.

[33] HAYBAR H, SHAHRABI S, REZAEEYAN H, et al. Endothelial cells：from dysfunction mechanism to pharmacological effect in cardiovascular disease[J]. Cardiovasc Toxicol,2019,19(1)：13-22.

[34] YU B, CHEN Q, LE BRAS A, et al. Vascular stem/progenitor cell migration and differentiation in atherosclerosis[J]. Antioxid Redox Signal,2018,29(2)：219-235.

[35] CHEN D, LI D, XU X B, et al. Galangin inhibits epithelial- mesenchymal transition and angiogenesis by downregulating CD44 in glioma[J]. J Cancer,2019,10(19)：4499-4508.

[36] HUANG Y C, TSAI M S, HSIEH P C, et al. Galangin ameliorates cisplatin-induced nephrotoxicity by attenuating oxidative stress, inflammation and cell death in mice through inhibition of ERK and NF-kappaB signaling[J]. Toxicol Appl Pharmacol,2017, 329(8)：128-139.

[37] 劉楊.基于HIF-1α梓葛凍干粉針抗腦缺血血管內(nèi)皮細胞凋亡和促新血管形成作用及其機制研究[D].重慶：西南大學,2018.

[38] 肖林焱,王海麗,陳毅,等.茜草炭不同極性部位化學成分組成及保護氧化損傷活性比較[J].中國中藥雜志,2017,42(23)：4604-4610.

[39] 邱模昌,余艷榮,張瑜,等.吳茱萸次堿對內(nèi)皮細胞損傷的保護效應及機制[J].時珍國醫(yī)國藥,2013,24(3)：580-582.

[40] HANSSON G K, ROBERTSON A K, S?DERBERG-NAUCLéR C. Inflammation and atherosclerosis[J]. Herz,2006,29(8)：297-329.

[41] F?RSTERMANN U, XIA N, LI H. Roles of vascular oxidative stress and nitric oxide in the pathogenesis of atherosclerosis[J]. Circ Res,2017,120(4)：713-735.

[42] DU F, ZHOU J, GONG R, et al. Endothelial progenitor cells in atherosclerosis[J]. Front Biosci (Landmark Ed),2012,17(6)：2327- 2349.

Study on Screening of TCM Prevention and Treatment of Atherosclerosis Based on Integrative Pharmacology and Spatiotemporal Effects

FU Xinying1, SUN Zhengji1, WU Lu2, ZHANG Wei1

To integrate multi-database resources to select effective Chinese materia medica combined with the characteristics of spatiotemporal effect of atherosclerosis (AS); To better guide the rational medication in clinic.The related targets of atherosclerosis and “spatiotemporal effect” phenotype were obtained by TCMIP, GeneCards and OMIM database, and TCM prevention and treatment of “time effect” of vascular endothelial dysfunction (VED), vascular endothelial injury (VEI), plaque formation (PF) and “space effect” of bone marrow microenvironment (BMM) were screened through the function of “reverse search for Chinese materia medica” in the collection of integrated pharmacological analysis tools in TCMIP database. The enrichment analysis of PPI, GO and KEGG was carried out on TCM prevention and treatment for spatiotemporal effect of AS through String and WebGestalt database; chord map and heat map were constructed by R language to focus on the correlation and difference of “effective Chinese materia medica spatiotemporal” and the best compatibility of medicine.The common targets of AS-VED, AS-VEI, AS-PF and AS-BMM were 924, 935, 885 and 1024 respectively. The results of reverse search of Chinese materia medica showed that Alpiniae Officinarum Rhizoma, Rhapontici Radix, Puerariae Lobatue Radix, and Dalbergiae Odoriferae Lignum could prevent vascular endothelial dysfunction in the early stage; Viticis Negundo Folium, Rubiae Radix et Rhizoma, Alpiniae Officinarum Rhizoma, Rhapontici Radix could treat vascular endothelial injury in the middle stage; Euodiae Fructus, Impatientis Semen, Puerariae Lobatue Radix, and Dalbergiae Odoriferae Lignum could eliminate atherosclerotic plaque in the advanced stage, and Chinese materia medica such as Gleditsiae Spina, Periplocae Cortex, Leonuri Herba, Sparganii Rhizoma, and Trichosanthis Fructus could play a role in bone marrow microenvironment, and their mechanism may be related to activating Rap1, Hippo, MAPK, Toll-like, TNF, AMPK, PI3K-AKT, NF-κB signal pathway and regulating multiple targets such as AKT1, IL-6, INS, TNF, NR3C1, PPARG, PTGS2, etc. The “multi-drug-multi-target-multi-pathway” model could play a more important role in the prevention and treatment of AS.The Chinese materia medica with specific time stages and spatial processes screened out in this study realize the differentiation and treatment at the micro level, provide targeted Chinese materia medica and the possibility of targets and pathways for TCM prevention and treatment of AS, and provide ideas for TCM to assist stem cell therapy.

integrated pharmacology; atherosclerosis; spatiotemporal effect; vascular endothelial cells; TCM

R271.917.41；R285.5

A

1005-5304(2021)04-0032-09

10.19879/j.cnki.1005-5304.202003451

國家自然科學基金（81774032）；湖南省自然科學基金杰出青年基金（2020JJ2024）；湖南省自然科學基金青年基金（2018JJ3391）；湖南省教育廳科學研究重點項目（19A374）；湖南中醫(yī)藥大學研究生培養(yǎng)質(zhì)量工程專項（2020CX62）

張偉，E-mail：zeen6463@sina.com

（收稿日期：2020-03-16）

（修回日期：2020-04-18；編輯：季巍?。?/p>