基于網絡藥理學的沙棘總黃酮治療心肌缺血的作用機制研究

魏志成　童東　楊娟　趙可惠　孟憲麗　張藝

[摘要] 該研究采用網絡藥理學虛擬篩選方法，旨在進一步研究沙棘總黃酮（total flavonoids of Hippophae rhamnoides，TFH）治療心肌缺血的活性成分和作用機制。首先通過TCMSP數據庫和PubChem數據庫，結合口服利用度和類藥性分析，篩選出沙棘化合物中的黃酮類化合物。篩選得到的7個化合物再通過ChemMapper服務器進行靶點的預測分析。將得到的潛在靶點導入MAS 3.0數據庫，結合KEGG數據庫進行靶點分析和通路分析。最后使用Cystoscope 3.3.0軟件繪制“化合物靶點作用通路”網絡圖。虛擬實驗預測得到了68個潛在靶點和60條信號通路，其中有31個靶點和23條通路直接或間接與心肌缺血有關。結果表明，TFH主要通過對鈣離子信號通路（calcium signaling pathway）、血管內皮生長因子信號通路（VEGF signaling pathway）和縫隙連接信號通路（gap junction）等信號通路的調節發揮協同作用，且與文獻報道吻合；提示其可能通過參與調節血小板聚集、脂質代謝、炎癥反應等過程，增強心功能和對血管內皮細胞保護作用，從而發揮抗心肌缺血的作用。

[關鍵詞] 沙棘； 沙棘總黃酮； 網絡藥理學； 心肌缺血； 機制研究

Action mechanism of total flavonoids of Hippophae rhamnoides in

treatment of myocardial ischemia based on network pharmacology

WEI Zhicheng， TONG Dong， YANG Juan， ZHAO Kehui， MENG Xianli， ZHANG Yi*

（Chengdu University of Traditional Chinese Medicine， Chengdu 611137， China）

[Abstract] In this study， a network pharmacological screening method was adopted to further study the active ingredients and action mechanism of total flavonoids of Hippophae rhamnoides（TFH） for the treatment of myocardial ischemia. Firstly TCMSP database and PubChem database were searched， and then the data were combined with oral bioavailability and drug analysis to screen flavonoids of H.rhamnoides compounds. Then predictive analysis was conducted for the 7 screened compounds by ChemMapper server.The obtained potential targets were imported into MAS 3.0. Database， and KEGG database was also used for targets analysis and pathway analysis. Finally Cystoscope 3.3.0 software was used to draw "compoundstargetspathway" network diagram. Virtual experiments predicted 68 potential targets and 60 signaling pathways， and 31 targets and 23 pathways of them were directly or indirectly associated with myocardial ischemia. The results showed that TFH played a synergistic rolemainly through the regulation of calcium signaling pathway， VEGF signaling pathway and gap junction signaling pathway， which was consistent with literature reports. These results indicated that it can enhance heart function， protect vascular endothelial cells， and fight against myocardial ischemia probably by regulating platelet aggregation， lipid metabolism， inflammation and other processes.

[Key words] Hippophae rhamnoides； TFH； network pharmacology； myocardial ischemia； mechanism

心肌缺血（myocardial ischemia），是指心臟的血液灌注減少，導致心臟的供氧減少，心肌能量代謝不正常，不能維持心臟正常工作的一種病理狀態[1]。冠狀動脈阻塞、主動脈供血減少和血壓降低，可直接導致心臟供血減少；心瓣膜病、血黏度變化和心肌本身病變同樣也會使心臟供血減少。其中冠狀動脈阻塞是造成心肌缺血的最主要也是最常見的原因，而使冠狀動脈阻塞的主要病因為冠狀動脈粥樣硬化。隨著人民生活水平的提高，不合理的飲食習慣和生活規律使得心肌缺血在我國的患病率呈逐年上升的趨勢，已成為中老年人的常見病和多發病之一，且逐年呈現年輕化的趨勢。

沙棘Hippophae rhamnoides為胡頹子科酸棘屬植物，是藏族習用藥材，其果實可以除肺腫瘤，化血，治培根病[2]。據2015年版《中國藥典》記載，沙棘有健脾消食，止咳祛痰，活血散瘀之功效[3]，其主要活性成分為黃酮類化合物，國內外研究表明其具有抗血栓和體外抗血小板聚集活性[4]，能改善血液流變學[5]，改善心臟功能及血流動力學[67]等影響血液系統的藥理作用。以TFH為主要成分的已上市藥物心達康膠囊在臨床應用中也表明具有良好的抗心肌缺血和改善血流動力學作用[89]。雖然目前對TFH的藥理研究已有較多報道，但針對其多成分、多靶點及多途徑治療心肌缺血的作用機制研究仍不夠系統和全面。該文運用網絡藥理學對沙棘黃酮類成分抗心肌缺血作用機制進行系統研究，希望能為后續實驗提供一定的理論基礎。

1 材料與方法

1.1 材料 TCMSP數據庫（http：//lsp.nwsuaf.edu.cn/tcmsp.php）；PubChem（有機小分子生物活性數據庫，http：//pubchem.ncbi.nlm.nih.gov）；ChemBioOffice2010；ChemMapper數據庫（http：//lilab.ecust.edu.cn/chemmapper/index.html）；MAS 3.0（生物分子功能注釋系統，http：//bioinfo.capitalbio.com/mas3/）；KEGG數據庫（http：//www.genome.jp/kegg/）；Cystoscope 3.3.0。

1.2 沙棘化學成分分子庫的準備 通過文獻考察和TCMSP數據庫（http：//lsp.nwsuaf.edu.cn/tcmsp.php）查詢共獲得了沙棘223個化學成分的分子結構。將上述223個化合物代入軟件分析，通過口服利用度和類藥性這2個決定藥物分子是否具有較好的口服生物利用度以及成藥性優劣的關鍵參數進行篩選[1011]。設定口服利用度（OB）≥30%，和類藥性（DL）≥0.18 2個條件，得到同時符合條件個化合物有33個。這些化合物再通過PubChem（有機小分子生物活性數據庫，http：//pubchem.ncbi.nlm.nih.gov/）進行結構辨別，篩選出2個苯環（A環與B環）通過中央三碳相互聯接而成的這類化合物（黃酮類化合物），滿足黃酮類結構共7個化合物分別是槲皮素（quercetin，PubChem CID：5280343）、異鼠李素（isorhamnetin，PubChem CID：5281654）、山柰素（kaempferol，PubChem CID：5280863）、花葵素（pelargonidin，PubChem CID：5281605）、表兒茶素（epicatechin，PubChem CID：72322）、兒茶素（cianidanol，PubChem CID：2353）、AC1LDVZ[（2R）5，7dihydroxy2（3hydroxy4methoxyphenyl）2，3dihydrochromen4one，PubChem CID：3593]，化合物具體信息見表1。

1.3 沙棘黃酮類化合物潛在靶點預測與篩選 用ChemBioOffice2010繪制 7個沙棘黃酮類化合物的結構，將其存為MOL2格式，并導入 ChemMapper服務器進行潛在靶點的預測分析，參數設置如下，Options for ChemMapper：Target Navigator；3D similarity method：SHAFTS；Similarity Threshold：1.2；Bioactivity Database： DrugBank。ChemMapper[12]是基于分子相似性方法的靶效藥理學及化學關系發現的多功能平臺，本文通過ChemMapper整合Drugbank數據庫已知的包含生物活性及靶點注釋信息的藥物分子，然后應用三維分子相似性方法SHAFTS[13]來計算查詢分子與藥物分子的相似度，從而預測查詢分子的潛在靶點?；趯π募∪毖l病機制的總結，進一步通過KEGG數據庫靶點功能分析比對心肌缺血病理過程所對應靶點，篩選出相應的靶點進行下一步分析。

1.4 GO富集分析與通路注釋 將預測獲取的靶點蛋白信息導入MAS 3.0數據庫，即可得到靶點的相關通路信息。參數設置如下，Species： Homo sapiens （Human）；Molecule Type： Protein；Database Symbols： Gen Bank ACC；Functions： Gene，mRNA，Protein，Pathway；Regulation Threshold： Up 2.0，Down 0.5。另外，KEGG包含的通路數據庫儲存了基因功能的相關信息，通過圖形來表示細胞內的生物學過程。利用此數據庫可以對由MAS 3.0得到的通路進行注釋及GO富集分析，以及結合KEGG數據庫進行通路分析。

1.5 網絡構建 通過Cystoscope 3.3.0建立“化合物靶點通路”網絡模式，繪制沙棘黃酮類化合物分子與靶蛋白作用網絡圖。以選取的7個沙棘黃酮成分、預測出的靶點和通路及通路注釋為節點（node），建立相互關系，導入Cytoscape 3.3.0軟件，運用network analyzer插件分析藏藥沙棘黃酮類成分治療心肌缺血的“化合物靶點通路”網絡圖。在網絡中，以節點（node）表示活性成分、靶點和通路。以邊（edge）相連接表示某一靶點為某化合物的潛在作用靶點和某一靶點與某一通路相關。

2 結果

2.1 靶點預測與篩選 通過ChemMapper數據庫對TFH 7個黃酮類成分的分析，共預測得到了68個潛在靶點，經靶點功能分析后確認其中有31個靶點直接或間接與心肌缺血有關，相關靶點信息見表2。各化合物所對應的作用靶點數不盡相同，其中槲皮素14個、異鼠李素4個、山柰素14個、花葵素15個、表兒茶素24個、兒茶素10個，單從作用層次上，提示表兒茶素在網絡中貢獻最大。在預測的靶點中，PTGS1（前列腺素G/H合酶1）、PTGS2（環氧合酶2）是花生四烯酸代謝通路中重要節點；ADRB1（β1腎上腺素能受體），ADRB2（β2腎上腺素能受體），ADRB3（β3腎上腺素能受體），DRD1（多巴胺D1b受體），DRD5（多巴胺D1a受體）等與鈣離子通道有關；MAPK1（絲裂原活化蛋白激酶1），MAPK14（絲裂原活化蛋白激酶14），PIK3CG（磷脂酰肌醇激酶CG），PIK3R1（磷脂酰肌醇3激酶調節亞基1），PIK3R2（磷脂酰肌醇3激酶調節亞基2），CDK2（細胞分裂蛋白激酶2）等與細胞生存和凋亡有關；PIK3CG，PTGS1，PTGS2，MAPK1，IFNG（干擾素γ）等與炎癥反應有關；PIK3CG，PIK3R1，PIK3R2等與脂質代謝有關。由此推斷，TFH治療心肌缺血可能通過參與調節血小板聚集、脂質代謝、炎癥反應等過程，增強心功能和保護血管內皮細胞等作用有關。

2.2 GO富集分析及通路分析 通過GO富集分析發現，見圖1～3，預測靶點在生物過程中，蛋白偶聯受體蛋白信號通路和信號傳導排位靠前；在分子功能中，ATP結合、蛋白質結合、轉移酶活性和核苷酸結合相關性最大；在細胞組分中，質膜和細胞核所占比例最大。通過KEGG通路分析共得到60條信號通路，23條直接或間接與心肌缺血有關，見表3。其中有2條與血小板聚集有關，分別是Fc epsilon RI和arachidonic acid metabolism信號通路；有3條與脂質代謝相關，分別是Fc epsilon RI，mTOR和Insulin信號通路；有8條與炎癥相關，分別是Tolllike receptor，MAPK，JakSTAT，Arachidonic acid metabolism，Leukocyte transendothelial migration，TGFbeta，Proteasome和Focal adhesion信號通路；有6條與細胞凋亡相關，分別是VEGF，Apoptosis，Cell cycle，p53，Gap junction和ErbB信號通路；有3條與心臟功能和血管舒張有關，分別是Ca2+，Fc epsilon RI和Phosphatidylinositol signaling system信號通路；有3條與血管韌性和通透性有關，分別是Leukocyte transendothelial migration，Adherens junction和Tight junction信號通路。

2.3 網絡構建 利用Cytoscape 3.3.0軟件建立了TFH的“化合物靶點通路”網絡圖，見圖4。直觀地呈現了6個化合物（其中化合物AC1LDVZP無對應靶點）、31個靶蛋白、23條作用通路與心肌缺血疾病之間的關系。分析結果表明， Calcium signaling pathway，VEGF signaling pathway和Gap junction 3條信號通路在網絡中網絡節點度排位靠前（作用于通路的靶點數均≥8），為網絡中的樞紐節點，提示它們可能是TFH抗心肌缺血核心作用途徑；化合物與靶點間存在著“一對多、多對一”的關系，體現了TFH其多成分、多靶點的抗心肌缺血作用機制；沙棘黃酮類成分在血小板聚集、脂質代謝、炎癥反應、細胞凋亡和心血管功能等相關通路上具有協同作用。

3 討論

網絡藥理學是在基于高通量組學數據分析、計算機虛擬計算及網絡數據庫檢索構建 “藥物基因靶點疾病”相互作用網絡，進一步利用網絡平衡的整體觀來分析藥物與機體相互作用的一門學科[14]。網絡藥理學的發展打破了“一個藥物、一個靶標、一種疾病”的傳統理念，將高選擇性單一靶點的新藥研究模式轉化為針對機體網絡調控的組合模式[15]。網絡藥理學這一研究策略為中藥和民族藥藥理學研究提供了嶄新的思路與方法，為中藥和民族藥多靶點、多途徑效應提供了可量化數據，實現了對單味藥多成分或者復方的作用機制研究的整體化、系統化[16]。

心肌缺血主要由于血液所含物質的變化和血管腔內產生病理改變的結果[1]。大多起因于膳食不合理、運動過少等生活因素導致的高脂血癥、高膽固醇血癥。久而久之導致動脈血管形成粥樣斑塊，而冠狀動脈粥樣硬化直接造成冠狀動脈血管通道變得狹窄，甚至阻塞，這也是心肌缺血最主要和常見的原因之一[17]。藏醫理論認為，心是臟腑之君，《四部醫典》[18]對心臟病的病因描述為“愁悲煩亂心不安，腹肌失眠瞋粗暴生”。心肌缺血的癥狀屬于藏醫“心悸”“心刺痛”“心熱癥”等病的范疇[19]，而冠狀動脈粥樣硬化所導致的心肌缺血，藏醫認為，一些病因一方面會使隆行太過發生偏盛的病變，或隆亂血塞；另一方面作用于血使赤巴受害而打破三因平衡[20]。

本研究發現大量預測靶點與血管內皮細胞保護和心肌細胞保護有關，且各個靶點、通路相互協同作發揮作用，提示這可能是TFH抗心肌缺血的關鍵。心肌缺血/再灌注損傷引起炎癥反應，大量炎癥介質的表達，造成心肌細胞的進一步凋亡，從而加劇心肌缺血/再灌注損傷[21]。另一方面，血小板聚集、炎癥反應、氧化應激反應會進一步損傷內皮細胞，促進血栓形成和加劇炎癥反應，進而加速冠狀動脈硬化和心肌缺血病理過程?；ㄉ南┧岽x信號通路在整個抗心肌缺血過程扮演著重要的角色，其代謝產物參與膽固醇代謝、內皮細胞單核細胞黏附、血小板聚集，調節血液黏度和血細胞功能等過程[22]；而且Src酪氨酸激酶的激活能增加細胞內鈣離子水平和激活磷脂A2（PLA2），進而激活花生四烯酸在質膜磷脂中的合成[4]。TFH對活性氧自由基（ROS）的清除作用[23]能減少心肌缺血造成的損傷。心肌缺血會導致心肌細胞產生大量ROS，再灌注也會進一步促使心肌細胞產生ROS[24]，加劇心肌細胞和內皮細胞的凋亡。而TFH可能通過MAPK14，PIK3CG，PTK2B，PIK3R1和PIK3R2等靶點調控白細胞跨內皮遷移信號通路（Leukocyte transendothelial migration），抑制NADPH氧化酶表達，進而減少ROS的分泌，保護心肌細胞和內皮細胞受到氧化應激損傷。

本文通過預測得到的Calcium signaling pathway，VEGF signaling pathway和Gap junction 3個通路為網絡中樞紐節點，與已有文獻報道吻合，證實了預測實驗的準確性。TFH對鈣離子轉運的影響，能增強心肌細胞的收縮力[25]和減少鈣離子超載[26]造成的細胞損傷。TFH通過多個腎上腺素能受體和多巴胺受體對鈣離子信號通路進行調控，提示這可能是總黃酮保護內皮細胞的機制之一。心肌和血管在缺血性缺氧環境下，機體會上調低氧誘導因子（HIF1α）的表達，誘導體內一些血管活性物質，如血管內皮生長因子（VEGF）和內皮素1（endothelin1，ET1）的表達[27]。VEGF在正常血管中無表達，但在動脈粥樣硬化斑塊中表達顯著增加，這也是斑塊形成的重要特征之一。而沙棘提取物能顯著減少由缺氧誘導的促炎細胞因子和VEGF水平的增加[28]，提示TFH能通過調節VEGF信號通路上游靶點（SRC，PIK3CG，PTGS2，MAPK14等）改善機體缺血缺氧狀態，從而下調VEGF的表達?？p隙連接和心律失常及細胞損傷等有著密切的關系，在心肌缺血早期給予阻斷縫隙連接細胞間通訊，發現可以明顯降低隨后發生的心肌壞死范圍[29]。已有文獻報道表明，槲皮素能顯著抑制縫隙連接細胞間通訊，從而達到心肌保護的作用[3031]。

TFH其主要有效成分為槲皮素、異鼠李素、山柰素等。目前TFH運用于人體主要是通過口服吸收，故在分子庫的準備階段對沙棘所含化合物進行了口服生物利用度和類藥性篩選，然后共得到7個滿足條件的化合物，這些化合物良好的口服生物利用度和類藥性決定了它們具有很好的成藥性，從藥物代謝角度上看可以很好代表了TFH的藥效。最后以上述7個化合物為研究對象，通過網絡藥理學的手段分析了上述成分可能的作用靶點及信號通路，初步揭示了其多成分、多靶點、多途徑的協同效應的作用特點。本文為下一步深入研究TFH抗心肌缺血藥效基礎研究及機制研究提供了一定的理論依據。

[參考文獻]

[1] Sheps.Lifethreatening ventricular arrhythmias in patients with silent myocardial ischemia due to coronary artery spasmNEJM[J].New Engl J Med，1992，326（22）：1451.

[2] 帝瑪爾·丹增彭措.晶珠本草[M].上海：上?？萍汲霭嫔?，2012.

[3] 中國藥典.一部[S].2015：184.

[4] Cheng J， Kondo K， Suzuki Y， et al. Inhibitory effects of total flavones of Hippophae rhamnoides L on thrombosis in mouse femoral artery and in vitro platelet aggregation[J]. Life Sci，2003，72（20）：2263.

[5] 王秉文，成華，袁秉祥，等.沙棘總黃酮對正常大鼠血液流變學的影響[J].沙棘， 2004，17（1）：28.

[6] 王秉文，馮養正，于佑民，等.沙棘總黃酮對正常人心功能及血流動力學的影響[J].西安交通大學學報：醫學版，1996（4）：37.

[7] 吳英，王毅.沙棘總黃酮對急性心衰犬心功能和血流動力學的影響[J].中國中藥雜志，1997（7）：429.

[8] 向楠，胡有志，石杰，等.心達康膠囊治療冠心病心絞痛（胸痹心血瘀阻證）臨床研究[J].中國臨床藥理學與治療學，2004，9（9）：1033.

[9] 陳光宇，賈秀蘭.心達康膠囊與片劑治療冠心?。ㄐ慕g痛型）181例臨床療效對比[J].中國中醫急癥，2000，9（6）：249.

[10] 張建永.丹參山楂組分配伍抗動脈粥樣硬化及作用機制研究[D].北京：中國中醫科學院，2013.

[11] Giulio Vistoli， Alessandro Pedretti， Bernard Testa. Assessing druglikeness——what are we missing？[J].Drug Discov Today，2008， 13（7/8）：285.

[12] Gong J， Cai C， Liu X， et al. ChemMapper： a versatile web server for exploring pharmacology and chemical structure association based on molecular 3D similarity method[J]. Bioinformatics， 2013， 29（14）：1827.

[13] Liu X， Jiang H， Li H. SHAFTS： a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening[J]. J Chem Inf Model， 2011， 51（9）：2372.

[14] Hopkins A L. Network pharmacology[J]. Nat Biotechn， 2007， 25（10）：1110.

[15] Zhang Yanqiong， Mao Xia， Guo Qiuyan， et al. Network pharmacologybased approaches capture essence of Chinese herbal medicines[J]. Chin Herb Med， 2016， 8（2）： 107.

[16] Li S， Fan T P， Jia W， et al. Network pharmacology in traditional Chinese medicine[J]. Evid Based Complement Alternat Med（eCAM）， 2014， 2014（1）：57.

[17] Cocco G， Gasparyan A Y. Myocardial ischemia in Wegener′s granulomatosis： coronary atherosclerosis versus vasculitis[J]. Open Cardiovas Med J， 2009， 4（1）：57.

[18] 宇妥·元丹貢布.四部醫典（藏文）[M].拉薩： 西藏人民出版社，1982.

[19] 第司·桑杰嘉措.藍琉璃[M].上海：上?？茖W技術出版社，2012.

[20] 趙可惠，王靜，仁青東主，等.藏醫娘塞病與中醫胸痹心痛辨證治療對比[J].時珍國醫國藥，2015，26（12）：2967.

[21] 董六一，陳志武.心肌缺血/再灌注損傷與炎癥反應[J].中國臨床藥理學與治療學，2008，13（5）：582.

[22] 郭歌，祖凌云，高煒.花生四烯酸代謝產物與動脈粥樣硬化[J].心血管病學進展，2014（5）：526.

[23] 李小潔，趙保路，侯京武，等.沙棘總黃酮對活性氧自由基的清除作用[J]. 中國藥理學通報，1990（2）：97.

[24] 吳蕊，凃玲.氧化應激與心血管疾病[J].心血管病學進展，2007， 28（1）：110.

[25] 王正榮，王玲，尹華虎，等.沙棘總黃酮對牽張所致心衰心肌細胞的收縮力學和鈣離子轉運的影響[J].航天醫學與醫學工程，2000， 13（1）：6.

[26] 程嘉藝. 沙棘總黃酮抗栓作用機制及有效成分初探[D]. 沈陽：遼寧中醫藥大學， 2008.

[27] 吳冬嗚， 唐燕， 韓藝，等. 大鼠心肌缺血再灌注后不同時間點心臟及肝臟中HIF1、VEGF的表達[J]. 南京醫科大學學報：自然科學版， 2012（8）：1047.

[28] Purushothaman J， Suryakumar G， Shukla D， et al. Modulation of hypoxiainduced pulmonary vascular leakage in rats by seabuckthorn （Hippophae rhamnoides L.）.[J]. Evid Based Complement Alternat Med： eCAM， 2011， 2011（5）： 574524.

[29] Tetsuji M， Yoshito O， Atsushi K， et al. Protective role of gap junctions in preconditioning against myocardial infarction[J]. Ajp Heart Circulat Physiol， 2004， 286（1）：214.

[30] Kim Y J， Sang G S， Choi K， et al. Recovery effect of onion peel extract against H2O2induced inhibition of gapjunctional intercellular communication is mediated through quercetin[J]. J Food Sci， 2014， 79（5）：H1011.

[31] Dongeun L， Bongjik S， Haengjeon H， et al. Quercetin， the active phenolic component in kiwifruit， prevents hydrogen peroxideinduced inhibition of gapjunction intercellular communication[J]. Bri J Nutr， 2010， 104（2）：164.

[責任編輯 張寧寧]