基于網絡藥理學研究南五味子對肝損傷的保護作用機制分析

姜海慧，唐 凱，張化為，姜 祎，宋小妹，黃文麗，鄧 翀

陜西中醫藥大學 陜西省秦嶺中草藥應用開發工程技術研究中心 陜西省中藥基礎與新藥研究重點實驗室，咸陽 712046

肝臟是人體內新陳代謝的中心站，被喻為“人體最大的化工廠”。肝損傷包括物理性肝損傷、化學性肝損傷、酒精性肝損傷、環境性肝損傷等[1]。臨床上保護肝臟的特異性藥物數量少，且易產生毒副反應，而具有多成分、多靶點作用的中藥及中藥方劑對于肝損傷有一定的促進再生和修復作用，但大多作用機制尚不十分清楚。因此探討保肝中藥治療肝損傷的多靶點、多途徑作用機制是目前的研究熱點之一。

網絡藥理學是建立在計算機虛擬計算、高通量組學數據分析等技術基礎上，利用了生物信息學方法，研究多成分-多靶點-多疾病的整合機制，因而具有整體性和系統性的優勢，現廣泛應用于中藥復方及單味藥的作用機制研究中[2]。

南五味子為木蘭科植物華中五味子(Schisandrasphenanthera.Rehdet Wils.)的干燥成熟果實，具有收斂固澀、益氣生津、補腎寧心的功效[3]，含有有木脂素、多糖、三萜、揮發油等化學成分[4]。近年來，臨床和實驗研究表明南五味子具有保肝作用[5]，其中木脂素類具有降低谷丙轉氨酶、抗炎、抗氧化、抗腫瘤等藥理作用[6]。南五味子多糖能夠致小鼠急性肝損傷后的ALT和AST水平明顯降低，使小鼠肝臟損傷的到有效改善[7]。南五味子能夠促進膽汁分泌，使肝內有毒物質加速排泄，南五味子中含有多種活性成分能夠使SOD、GSH-Px等酶的生物活性增強，具有促進蛋白質生物合成、抗脂質過氧化、促進肝糖原生成、保護肝細胞膜等作用，同時南五味子能夠促進肝臟解毒過程，使肝臟收到良好的保護，并能夠使受肝損傷組織再生，能夠使外界因素導致肝組織中不斷升高的轉氨酶有效降低[8]。但對南五味子而言究竟哪些成分都發揮保肝功效，其治療肝損傷的作用機制均不清楚。基于此，本研究通過網絡藥理學，整理南五味子活性成分、關鍵靶點、作用通路之間的關系，結合分子對接驗證的方法，探討南五味子保肝的作用機理，從而為南五味子對治療肝損傷提供理論依據。

1 材料與方法

1.1 南五味子活性成分的篩選

在中藥系統藥理學數據庫與分析平臺TCMSP(http://lsp.nwu.edu.cn/tcmsp.php)和一些已發表的文獻中檢索南五味子的所有化學成分，通過OB(口服生物利用度)≥20%，DL(類藥性)≥0.1及Lipinski原則[9]篩選出活性成分。

1.2 南五味子靶點的整理

在PubChem(https://pubchem.ncbi.nlm.nih.gov/)檢索活性成分對應的2D結構式，通過Swiss數據庫(http://www.swisstargetprediction.ch/)進行靶點的預測富集，將TCMSP數據庫檢測到的靶點與其整合，使用Uniprot數據庫(https://www.uniprot.org)，搜索靶點的基因名，去重整合。

1.3 收集疾病相關靶點

在Genecard與OMIM中搜索關鍵詞“hepatic injury”，整合肝損傷疾病相關靶點的基因名。然后通過Venny2.1(http://bioinfogp.cnb.csic.es/tools/venny/index.html)把南五味子的活性成分的預測靶點與肝損傷疾病的靶點繪制韋恩圖并取交集，整理共同的靶基因。

1.4 構建關鍵靶點的蛋白互作網絡

STRING數據庫(https://string-db.org/)收錄了計算、實驗、文獻等方法的蛋白相互作用關系，是目前常用構建蛋白-蛋白相互作用(PPI)網絡的數據庫。將上述分析取到的關鍵靶基因導入STRING數據庫，選擇物種為人，得到南五味子治療肝損傷的PPI網絡，在構建的PPI網絡中，節點的大小和顏色深淺反映自由度(degree)值大小，邊粗細反映組合分數(combine score)大小。分析PPI網絡，計算網絡中節點度值，度值越大代表節點在網絡中越處于核心地位。用其得到的“TSV”文件通過Cytoscape(版本3.2.1)中NetworkAnalyzer工具進行拓撲分析。

1.5 構建“藥物成分-靶點”網絡并分析

Cytoscape3.2.1(可視化網絡構建與分析軟件)是構建與研究中藥網絡的常用軟件。導入南五味子活性成分與肝損傷相關疾病的靶點，構建“藥物成分-靶點”網絡。

1.6 關鍵靶點的GO分析與KEGG富集分析

使用半編程性質的腳本軟件R project，可以用于實現多種生物信息分析與結果可視化。使用其生物學分析功能，對南五味子治療肝損傷的關鍵靶點進行GO與KEGG分析，設置P<0.05，Q<0.05，其中GO分析包括生物過程(BP)、細胞成分(CC)和分子功能(MF)分析，分析結果以柱形圖的形式予以展示。

1.7 活性成分與關鍵靶點的分子對接驗證

將上述篩選出的南五味子的活性成分在PubChem平臺中查詢其結構式，并通過SMILES導出，在DISCOVERY STUDIO軟件中打開并在Macromolecules模塊下對每個小分子進行優化處理，PPI網絡中“degree”值前20的靶點作為受體，在PDB數據庫(http://www.rcsb.org/)中下載其靶點的蛋白結構，然后導入DISCOVERY STUDIO軟件進行分子結構常規預處理，接著選擇Libdock對接模塊中設定docking preference為user specified，設定minimum LibDockscore為105分，其他參數取默認值，進行分子對接。打分函數LibDock Score分數越高，其對接的結合的活性越好。

2 結果

2.1 南五味子活性成分的篩選

一共整理了南五味子的化學成分共149個，經篩選，共有31個活性成分納入研究(見表1)。

表1 南五味子的主要活性成分Table 1 Main active ingredients of Schisandra sphenanthera

2.2 收集南五味子有效活性成分的作用靶點

通過TCMSP和Swiss檢索的靶點，收錄的南五味子中成分的作用靶點共1 639個。使用Uniprot數據庫收集作用靶點的基因名，刪除無效與重復靶點，共獲取南五味子有效活性成分作用靶點537個。

2.3 收集肝損傷疾病的相關基因

檢索OMIM和Genecard數據庫，整理并去重，共收集到肝損傷相關疾病靶點3 493個。將上述南五味子有效成分的作用靶點映射到其中，共獲得交集358個，見圖1，將其定義為南五味子治療肝損傷的“關鍵靶點”。

圖1 南五味子作用靶點與肝損傷相關靶點交集Fig.1 Intersection of target of S.sphenantheraand related target of hepatic injury

2.4 構建并分析關鍵靶點的PPI網絡

將上述358個關鍵靶點輸入STRING數據庫，構建成分與疾病靶點的PPI網絡，PPI網絡中包含節點308個，邊2 036條。將其“TSV”文件導入Cytoscape，構建南五味子與肝損傷靶點的PPI網絡，并通過R語言制作柱形圖，其中展示了PPI網絡中Degree值較高的前20位，見圖2，將其定義為“核心靶點”。選取 degree值排名前六的核心靶點：PIK3CA、APP、SRC、MAPK1、MAPK3、STAT3。其中PIK3CA(磷脂酰肌醇-3-激酶)在細胞的生長、增殖、凋亡、血管生成、自吞噬等過程中發揮著極其重要的生物學功能[10]，MAPK1與MAPK3(絲裂原活化蛋白激酶)調節著細胞的生長、分化、對環境的應激適應、炎癥反應等多種重要的細胞生理或病理過程[11]，STAT3(信號傳導及轉錄激活因子3)具有顯著的抗凋亡和促有絲分裂活性,能夠上調一系列細胞存活及增殖相關基因[12]。

圖2 PPI網絡中degree值較高的靶點Fig.2 Targets with higher degree in PPI networks

2.5 構建“藥物成分-靶點”網絡并分析

將南五味子活性成分、關鍵靶點輸入Cytoscape軟件，構建“藥物成分-靶點”網絡圖(見圖3)，在網絡圖中，藍色圓形代表活性成分所作用的靶點，綠色棱形代表南五味子中所含有的活性成分。分析網絡圖，可知網絡圖中共包含節點(node)388個，節點與節點之間的作用關系，即邊(edge)1 267條，展示了南五味子治療肝損傷多成分、多靶點的作用。

2.6 GO生物學過程富集與KEGG通路富集分析

分析358關鍵靶點主要富集的細胞學組分、分子功能與生物學過程。主要進行靶點在生物學過程的分析，結果顯示南五味子治療肝損傷的關鍵靶點共在2 794個GO term上富集。其中所涉及的生物信息學功能與過程包括positive regulation of MAPK cascade、response to toxic substance、peptidyl-serine phosphorylation等生物學過程密切相關。圖4展示了前20位GO分析BP的結果。

對南五味子治療肝損傷的358個關鍵靶點進行KEGG富集分析，分析結果顯示共富集在168條通路上，其中包括Neuroactive ligand-receptor interaction、Proteoglycans in cancer、Prostate cancer、EGFR tyrosine kinase inhibitor resistance、Hepatitis B、ErbB signaling pathway等信號通路。見圖5，展示了前20位KEGG通路富集分析結果。其中參與EGFR tyrosine kinase inhibitor resistance通路的核心靶點有9個，分別為EGFR、PIK3CA、MAPK1、JAK2、AKT1、STAT3、MAPK3、JAK1、SRC(圖6)，在排名前五的信號通路里核心靶點存在數量最高，因此，將此信號軸作為南五味子對肝損傷治療發揮藥效的通路之一。并且PIK3CA、MAPK1、STAT3三個靶點在前10的信號通路里涉及到4個以上的信號通路，這些通路主要涉及到炎癥的發生、氧化應激、癌細胞的抑制等生物學過程。

圖3 南五味子成分與靶點網絡Fig.3 Compounds-targets interaction network of S.sphenanthera

圖4 關鍵靶點的GO富集分析Fig.4 GO enrichment analysis for key targets

圖5 關鍵靶點的KEGG通路富集分析Fig.5 KEGG pathway enrichment analysis of key targets

圖6 EGFR酪氨酸激酶抑制劑耐藥Fig.6 EGFR tyrosine kinase inhibitor resistance

2.7 分子對接結果

將南五味子31個活性成分與導致肝損傷疾病重點相關的靶蛋白PIK3CA、MAPK1、STAT3做分子對接。在DISCOVERY STUDIO中對接并計算RMSD(均方根偏差)，普遍認為RMSD≤0.2 nm時能說明此方法能較好的重現結合模式，對接的結果比較可靠[13]。當對接大于105或者高于原配體的打分值時，可認為對接有較好的結合活性，結果顯示出南五味子活性成分中有14個成分與PIK3CA有較好的的結合活性，有15個成分與MAPK1有較好的結合活性，有12個成分與STAT3有較好的結合活性。見表2。其中新南五味子酸C(neokadsuranic acid C)與PIK3CA和STAT3結合(圖7和8)，和STAT3結合五味子酯乙(interiotherin B)與MAPK1結合(圖9)其對接效果較好。

表2 分子對接結果Table 2 Molecular docking results

續表2(Continued Tab.2)

靶點名稱Target namePBD IDMOLID化合物Composition對接分數LibDock scoreMOL009226Schisantherin J107.734MOL009228Schisantherin M116.231MOL009235Angusifolin B107.446MOL000492(+)-Catechin107.924MOL000432Linolenic acid115.775MOL000749Linoleic119.706

圖7 PIK3CA與neokadsuranic acid C對接Fig.7 The molecular docking of PIK3CA and neokadsuranic acid C

圖8 STAT3與neokadsuranic acid C對接Fig.8 The molecular docking of STAT3 and neokadsuranic acid C

3 討論

經本次實驗文獻挖掘，發現安五脂素((+)-anwulignan)、長南酸(changnanic acid)、胡蘿卜苷(daucosterol)、亞麻酸(linolenic acid)新南五味子酸C(neokadsuranic acid C)、五味子酯乙(interiotherin B)等31種具有潛在抗肝損傷成分，其涉及到治療肝損傷的機制包含抑制氧化應激、炎癥反應、細胞纖維化、細胞凋亡等[14]。

根據KEGG通路分析表明了肝損傷的相關通路可能通過及神經信號傳遞通路(Neuroactive ligand-receptor interaction)、癌癥蛋白多糖信號通路(Proteoglycans in cancer)、前列腺癌信號通路(Prostate cancer)、EGFR酪氨酸激酶抑制和耐藥信號通路(EGFR tyrosine kinase inhibitor resistance)、乙型肝炎病毒信號通路(Hepatitis B)、ErbB信號通路(ErbB signaling pathway)等信號通路起治療作用。其中富集到EGFR tyrosine kinase inhibitor resistance通路上的關鍵靶點有EGFR、PIK3CA、MAPK1、JAK2、AKT1、STAT3、MAPK3、JAK1、SRC，在排名前五的信號通路里存在數量最高，EGFR(表皮生長因子受體)屬于是受體型酪氨酸激酶(RTKs)家族的成員之一[15]，是一個巨大的跨膜糖蛋白，EGF可以刺激細胞分裂，與組織的修復和再生有關。EGF與EGFR結合，激活酪氨酸激酶區，使其磷酸化，促進細胞的分裂增殖及正常細胞的惡性突變，它在細胞的凋亡、增殖、分化、遷移和細胞周期循環發揮著十分重要的作用[16,17]。有研究表明，在肝損傷和修復過程中，ErbB1的表達升高[18]，從早期炎癥和肝細胞增殖到纖維化發生和腫瘤轉化，EGFR信號通路在肝損傷和炎癥中都擔任著重要的角色[19]，故以此通路為切入點進一步尋找南五味子治療肝損傷的作用靶標。

圖9 MAPK1與interiotherin B對接Fig.9 The molecular docking of MAPK1 and interiotherin B

挖掘到潛在抗肝損傷成分的核心靶點為EGFR、PIK3CA、MAPK1、JAK2、AKT1、STAT3、MAPK3、JAK1、SRC等，其中PIK3CA、MAPK1、STAT3三個靶點在前10的信號通路里涉及到4個以上的信號通路，這些通路主要涉及到炎癥的發生、氧化應激、癌細胞的抑制等生物學過程。基于此本研究以PIK3CA、MAPK1、STAT3為靶點，與31個潛在抗肝損傷成分進行分子對接分析，對其生物活性進行預測。

機體遭受有害刺激引起活性氧簇(ROS)生成增加或從而導致氧化-抗氧化失衡，引起細胞氧化應激狀態(OS)[20]。氧化應激可激活MAPK/ERK、PI3K/Akt、JAK/STAT、Nrf2/ARE、Wnt/β-catenin 等一系列信號通路及其相關響應因子，造成細胞外基質(ECM)成分代謝及相關因子表達異常，從而介入骨節、心血管、肝、腎等組織器官ECM代謝異常疾病的發生發展[21,22]。另外MAPK、STAT激活的信號通路對于維持細胞的正常生理功能以及在炎癥反應中占有重要地位，如LPS可以誘導炎癥細胞產生致炎因子TNF-α、IL-1、IL-6等均與這些信號通路緊密相關[23]。表明PIK3CA、MAPK1、STAT3三個靶點與肝損傷密切相關。

通過相關網絡圖的構建可直觀的顯示出南五味子治療肝損傷的多成分、多靶點、多途徑等共同作用的特點，符合中醫藥治療疾病的特點，相比于西藥單基因-單靶點的作用機制，優勢明顯。本研究通過分子對接驗證了網絡藥理學靶點預測的可靠性，為后續進一步驗證奠定了基礎，所得結果也與現有文獻報道高度符合，但本研究結果只是理論上的預測，藥物在體內作用機制的確認仍需要進一步通過整體動物實驗的驗證。