基于網絡藥理學探討隔山消積顆粒治療功能性消化不良的物質基礎及作用機制*

潘大波 ，姜德文 ，龍治峰 ，李小飄 ，王小霞 ，何文生

（1. 黔東南民族職業技術學院，貴州 凱里 556001； 2. 重慶市大足區人民醫院，重慶 402360）

功能性消化不良（FD）是常見的功能性胃腸道疾病，多表現為早飽、上腹痛等［1］，在普通人群中的發病率為5%～20%［2］。FD 發病機制復雜，病程較長，易反復發作，嚴重影響患者的生活質量［3］。中醫認為，FD 與“胃脘痛”“胃痞”“胃滿”的臨床癥狀相似。中醫治療FD 臨床療效顯著，中藥較緩和，長期服藥毒副作用小，患者臨床癥狀改善明顯［4］。目前，用于治療FD的中藥方劑有三仁湯、半夏瀉心湯、復方藏茴香和枳實消痞丸等。苗藥復方隔山消積顆粒（以下簡稱苗藥復方）由隔山消和馬蘭草2 味苗藥配伍。方中，隔山消具有補肝腎、強筋骨、益精血、健脾消食、解毒療瘡功效；馬蘭草又名魚鰍串、路邊菊等，有理氣、消食、清利濕熱之功效，用于治療胃脘脹痛、痢疾、水瀉、尿路感染等。本研究中采用網絡藥理學和分子對接技術，從苗藥復方多成分、多靶標的視角闡明其治療FD 的物質基礎與作用機制，為后續研究提供理論基礎。現報道如下。

1 資料與方法

1.1 活性成分篩選

通過SciFiner、中國知網、萬方、維普數據庫檢索獲取苗藥復方中隔山消、馬蘭草的化學成分。隔山消Cynanchum wilfordii是蘿藦科鵝絨藤屬植物，其基源植物有耳葉牛皮消Cynanchum auriculumRoyle ex Weight.、隔山 牛 皮 消Cynanchum wilfordii（Maxim.）Hemsl和戟葉牛皮消Cynanchum bungei，且隔山消與傳統中藥白首烏的主要來源一致，故以“隔山消”“Cynan?chum wilfordii”“ 耳 葉牛 皮消 ”“ 隔山 牛皮 消 ”“戟 葉 牛皮消”“白首烏”為關鍵詞檢索獲取隔山消的化學成分；馬蘭草為馬蘭Kalimeris indica（L.）Sch - Bip 的干燥全草，故以“馬蘭”“Kalimeris indica”“馬蘭草”為關鍵詞檢索獲取馬蘭草的化學成分。運用ADMETlab 2.0 軟件（https：// admetmesh. scbdd. com/）［5］在線預測其化學成分的性質，選取人體腸道吸收（HIA）在0～3 且符合Lipinski Rule（MW ≤ 500，logP ≤ 5，Hacc ≤ 10，Hdon ≤ 5，滿足3個及以上條件）的化學成分作為苗藥復方的活性成分。

1.2 活性成分靶標預測及FD 疾病靶標

采用基于配體的靶標預測在線服務器SwissTarget?Prediction［6］預測苗藥復方活性成分的作用靶標。以“Functional dyspepsia”作為關鍵詞，分別在GenCLip 3數據庫（http：//ci.smu.edu.cn/genclip3/analysis.php）［7］和 GeneCards 數據庫（https：// www. genecards. org/）［8］檢索收集靶標，合并去重后運用Uniprot數據庫（https：//www.uniprot.org/）校正整合，獲得FD 相關靶標。將FD疾病靶標與苗藥復方活性成分作用靶標進行疊合，獲得苗藥復方活性成分治療FD的靶標。

1.3 藥味-活性成分-治療FD 靶標相互作用網絡

將獲得的苗藥復方活性成分治療FD 的靶標映射到相應活性成分中，獲得活性成分與其對應治療FD 靶標信息。采用Cytoscape 3.8軟件［9］構建苗藥復方藥味-活性成分 - 靶標相互作用網絡，運用插件Centi?ScaPe2.2 分析，提取網絡中度值大于中介值的節點作為關鍵節點，即苗藥復方治療FD 的關鍵活性成分及關鍵作用靶標。

1.4 關鍵靶標-靶標相互作用網絡構建

將獲得的苗藥復方治療FD 的關鍵作用靶標導入String 11.5 數據庫（https：/ /string-db.org/）［10］，預測靶標與靶標間的相互作用，選取有證據的相互作用進行分析。運用Cytoscape 3.8軟件進行可視化顯示。

1.5 關鍵靶標GO 聚類分析和KEGG 通路富集分析

將獲得苗藥復方治療FD 的關鍵作用靶標導入David 6.8 數據庫（https：// david. ncifcrf. gov/）［11］進行GO 聚類分析和KEGG 通路富集分析。選取同時≥4 個基因且P<0.05的項目進行分析。

1.6 分子對接

選取度值排名靠前的苗藥復方治療FD 的關鍵作用靶標，運用Autodock vina 軟件［12］將苗藥復方活性成分對接到上述靶標中，考察活性成分與靶標的結合親和力及作用模式，對接打分熱圖用TBtools［13］顯示。

2 結果

2.1 活性成分及其作用靶標

通過文獻檢索，收集到隔山消和馬蘭草的化學成分分別為116個和114個。隔山消的化學成分主要為C21甾苷類、多糖類、苯酮類等，馬蘭草的化學成分主要為黃酮類、酚酸類、三萜類、揮發油等。采用Chemdraw 2015軟件繪制上述化學成分的結構并用ADMETlab 2.0軟件在線預測其性質，選取HIA 在0～3 且符合Lipinski Rule（MW ≤ 500，logP ≤ 5，Hacc ≤ 10，Hdon ≤ 5，滿足3個及以上條件）的化合物。符合條件的隔山消和馬蘭草的化學成分分別有35 個和71 個，其中1 個活性成分（β-谷甾醇）相同，即得到105個苗藥復方活性成分。

采用SwissTargetPrediction 數據庫預測其作用靶標，4 個化學成分（Diallyl carbonate，cynabungolide，coniferylalcohol和3，7，11-trimethyl-1，6-dodecadien-3，10，11 - triol）未預測到靶標，共作用于 821 個靶標，即101個化學成分作用于821個靶標。

2.2 FD 相關靶標

運用“Functional dyspepsia”作為關鍵詞在GenCLiP 3數據庫中檢索到127 個靶標、GeneCard 數據庫中大于Relevance score 中介值（2.87）的靶標 1 267 個，Uniprot數據庫標準化且去重后得到FD 相關靶標1 304 個。將FD 相關靶標與苗藥復方活性成分預測靶標疊合，得到共有靶標293 個（圖1），即293 個靶標可能是苗藥復方治療FD的作用靶標。

圖1 苗藥復方活性成分預測靶標與FD相關靶標韋恩圖Fig.1 Venn diagram of predicted targets of active ingredients of Miao medicine compound Geshan Xiaoji Granules and FD - related targets

2.3 藥味-活性成分-靶標相互作用網絡

將2.2 項下的共有靶標映射到100 個活性成分（其中1 個活性成分未找到），即苗藥復方的100 個活性成分通過作用于293 個靶標發揮治療FD 作用。獲得由2 442 條邊、395 個節點（其中，藥味 2 個、活性成分 100個、靶標293 個）組成的苗藥復方藥味-活性成分-靶標相互作用網絡。

度值不低于40的活性成分有20個（表1），其中來源于隔山消8個、馬蘭草12個。度值排名前10的化學成分有12-O-ikemaoyl-20-dihydro-lineolon、告達亭、芫花素、開德苷元、羅索他明、漢黃芩素、千層紙素A、山柰酚-7-葡萄糖苷、methyltrans-ferulate和加加米寧。

表1 度值不低于40的活性成分Tab.1 Active ingredients with a degree ≥40

度值不低于22 的靶標有20 個（表2），度值排名前10 的靶標為細胞色素P450 19A1（CYP19A1）、雄激素受體（AR）、碳酸酐酶Ⅱ（CA2）、雌激素受體β（ESR2）、睪丸特異性雄激素結合蛋白（SHBG）、雌激素受體α（ESR1）、HMG-CoA 還原酶（HMGCR）、花生四烯酸5-脂肪氧化酶（ALOX5）、聚［ADP - 核糖］聚合酶 - 1（PARP1）和乙酰膽堿酯酶（ACHE）。選取度值大于中介值（12.364 5）的節點作為關鍵節點，提取關鍵節點間的相互作用網絡（圖2），共2 個藥味、70 個化學成分、61 個靶標和125 1條邊。

圖2 苗藥復方隔山消積顆粒藥味-活性成分-靶標相互作用網絡Fig.2 Drugs-active ingredients-targets interaction network of Miao medicine compound Geshan Xiaoji Granules

表2 度值不低于22的靶標Tab.2 Targets with a degree ≥22

2.4 靶標- 靶標相互作用網絡

將61個關鍵靶標導入Srting數據庫，用Centoscape 3.8軟件繪制靶標- 靶標相互作用網絡（圖3）。度值大于20 的靶標有17 個，排名靠前的有表皮生長因子受體（EGFR）、酪氨酸蛋白激酶（SRC）和環加氧酶2（PTGS2）等。

圖3 靶標-靶標相互作用網絡Fig.3 Targets-targets interaction network

2.5 靶標功能GO 聚類分析和KEGG 通路富集分析

將從苗藥復方藥味-化學成分-靶標相互作用網絡提取的61個關鍵靶標導入David數據庫進行GO功能聚類和KEGG通路富集分析，滿足≥4個基因且P<0.05的細胞組成（CC）、生物過程（BP）和分子功能（MF）的分別得到18，74，35條項目。CC 主要集中在nucleus（核）和nucleoplasm（核質）等，BP主要集中在signal transduction（信號轉導）和negative regulation of apoptotic process（凋亡過程的負調控）等，MF 主要集中在protein binding（蛋白質結合）和ATP binding（ATP 結合）等，苗藥復方關鍵作用靶標的GO 聚類分析見圖4。提示苗藥復方活性成分通過作用于細胞核、細胞核質及細胞膜等影響能量代謝、蛋白結合等，調控信號轉導和蛋白質磷酸化等，達到治療FD的目的。

圖4 苗藥復方關鍵作用靶標的GO聚類分析Fig.4 GO clustering analysis of key targets of Miao medicine compound Geshan Xiaoji Granules

滿足 ≥ 4 個基因且P<0.05 的 KEGG 通路有 85 條，其中與FD 相關的靶標調控的通路為Fc γ R 介導的吞噬作用（Fc gamma R-mediated phagocytosis）、逆行內源性大麻素信號（retrograde endocannabinoid signaling）、膽堿能突觸（cholinergic synapse）、Toll 樣受體信號通路（toll-like receptor signaling pathway）、自然殺傷細胞介導的細胞毒性（natural killer cell mediated cytotoxicity）、TRP 通道的炎性介質調節（inflammatory mediator regula?tion of TRP channels）和花生四烯酸代謝（arachidonic acid metabolism）等。苗藥治療FD 的KEGG 通路富集分析結果見圖5。

圖5 苗藥復方治療FD的KEGG 通路富集分析Fig.5 KEGG pathway enrichment analysis of Miao medicine compound Geshan Xiaoji Granules in the treatment of FD

2.6 分子對接分析

選取苗藥復方藥味-活性成分-靶標相互作用網絡中靶標度值較大，且與目前研究發現FD 的發病因素關系較大的靶標CA2，HMGCR，ALOX5，PTGS2，ACHE作為受體（表3），將苗藥復方活性成分作為配體對接到受體中，分子對接打分熱圖見圖6。

圖6 苗藥復方活性成分分子對接打分熱圖Fig.6 Heat map of the molecular docking scoring of active ingredients of Miao medicine compound Geshan Xiaoji Granules

表3 分子對接的FD相關的靶標蛋白Tab.3 FD-related target proteins in the molecular docking

對接打分低于-8.0 kcal/mol的活性成分-靶標作用有119 對，其中打分最好的是活性成分cynandione E（白首烏酮E）與ACHE 間的作用，為 - 11.4 kcal/ mol。為進一步探討靶標與活性成分的相互作用模式，分別選取每個靶標與活性成分打分最優的組合進行分析（圖7）。

與CA2 結合最好的活性成分為加加米寧，對接打分為 - 9.3 kcal/ mol。加加米寧（gagamine）與 CA2 和THR200 形成 1 個氫鍵、ASN62 側鏈氨基形成 2 個氫鍵、ASN67側鏈羰基形成2個氫鍵，共5個氫鍵，與HIS64和HIS94 形成π - π 堆積作用（圖7 A）。在HMG - CoA 還原酶體系中，與青陽參苷元（qingyangshengenin）對接打分為- 8.8 kcal/ mol。由于青陽參苷元含有羥基較多，與HMG-CoA 形成很強的極性相互作用，共存在11 個氫鍵（圖7 B）。α-香樹脂醇（α-amyrin）嵌入ALOX5口袋中，與ASN554形成2個氫鍵（圖7 C）。cynandione E 與ACHE 的對接打分為 - 11.4 kcal/ mol。cynandione E 占據ACHE的2個活性空腔，且在2個空腔均與ACHE有氫鍵作用，即與SER293，PHE295，TYR341，THR75形成5個氫鍵（圖7 D）。大黃酚（chrysophanol）很好嵌入PTGS2 活性口袋中，與SER530和TYR355形成2個氫鍵（圖7 E）。

圖7 苗藥復方活性成分與靶標作用分子對接模式Fig.7 Molecular docking mode of the active ingredients of Miao medicine compound Geshan Xiaoji Granules and targets

3 討論

研究表明，FD的發病原因包括胃腸動力障礙、內臟高敏感性、腦- 腸軸功能紊亂、精神心理因素、十二指腸黏膜屏障受損、低度炎癥和腸道菌群失調等。苗藥復方由苗族常用藥隔山消和馬蘭草配伍，采用水煎煮制備而成，具有消積行氣功效，用于脾胃氣滯所致的食積內停、脘腹脹痛等。

本研究結果顯示，苗藥復方隔山消積顆粒中有100個活性成分與FD 相關的293 個靶標作用。上述靶標主要 集 中 在 AR，ESR2，ESR1，CA2，HMGCR，ALOX5，ACHE 等。其中，AR 能通過調控 Bmp 分泌，調控 Wnt 信號而控制腸干細胞增殖［14］；ESR1 和 ESR2 與胃食管反流病、食管癌、消化性潰瘍、胃癌、炎性腸病、腸易激綜合征和結腸癌的發生關系密切［15］；CA2 主要分布在胃腸道、腎遠曲小管等，與消化系統炎癥相關［16］；HMGCR抑制劑對幽門螺桿菌誘導的胃炎模型小鼠具有抗炎作用［17］；ALOX5是花生四烯酸合成白三烯的關鍵酶，產物白三烯與炎性疾病密切相關，包括胃炎、哮喘等［18］；ACHE 抑制劑阿考替胺能通過抑制ACHE 的活性來增強胃肌肉收縮，提高胃動力［19］。選擇性PTGS2抑制劑不影響胃十二指腸黏膜，同時具有與非選擇性非甾體抗炎藥（NSAID）相似的抗炎和鎮痛功效［20］。上述靶標主要調控幽門螺桿菌感染中的上皮細胞信號傳導、膽堿能突觸、腫瘤壞死因子（TNF）信號通路和花生四烯酸代謝等通路，涉及炎癥、腦- 腸軸功能紊亂、精神心理因素等，與FD的發生、發展密切相關。

由表1可知，苗藥復方活性成分-靶標相互作用網絡獲得度值大于40的20個活性成分有告達亭、開德苷元、加加米寧、漢黃芩素、neoechinulin A和大黃素甲醚。其中，告達亭在隔山消的含量為（1.354 ± 0.609）μg/g［21］，從5.7 kg 耳葉牛皮消中可提取80 mg 告達亭、32 mg 開德苷元和 50 mg 加加米寧［22］；從 20 kg 馬蘭全草中可提取漢黃芩素22 mg、千層紙素A 6 mg、7，4'-二羥基異黃酮10 mg［23］；從 15 kg 馬蘭全草中可提取芫花素 22 mg、aurantiamide金色酰胺醇10.5 mg、neoechinulin A 16.3 mg、鄰苯二甲酸二丁酯 5.6 mg、大黃素甲醚 53.4 mg［24-26］；從11 kg 馬蘭中可提取山柰酚-7-葡萄糖苷3.31 ng、methyltrans- ferulate 3.8 mg、槲皮素 6.3 mg、山柰酚5.7 mg［27］。在隔山消或馬蘭草中的含量較大，均超過1 mg/kg，可能是苗藥復方的特異性成分。

隔山消中的C21甾苷類活性成分在治療FD中起著關鍵作用。白首烏提取物能調控Toll 樣受體4/骨髓分化因子88/核轉錄因子κB（TLR4/MyD88/NF- κB）信號通路，抑制炎性反應［28］。C21 甾苷類成分能上調核因子E2 相關因子2（Nrf2）和血紅素氧合酶 - 1（HO - 1），抑制H2O2誘導的氧化損傷和炎性反應，其中告達亭能減輕肝中的炎性癥狀，下調白細胞介素6（IL-6）、白細胞介素1β（IL-1β）和單核細胞趨化因子的表達［29］。加加米寧抑制脂質過氧化物的形成，具抗氧化性［30］，且抑制5 - 脂肪氧化酶的半數抑制率（IC50）為26 μmol/ L。千層紙素A 通過保護腸黏膜和改變腸道微生物群來緩解結腸炎［31］。馬蘭草中的漢黃芩素激活Nrf2，促進抗氧化酶包括醌氧化還原酶-1（NQO-1）、谷胱甘肽硫轉移酶（GST）、HO-1、超氧化物歧化酶1（SOD1）和超氧化物歧化酶2（SOD2）在肝細胞中的表達，并抑制NF-κB 調節的促炎信號傳導，即通過激活Nrf2 的同時改善氧化應激和炎性反應［32］。槲皮素可抑制ERK和p38MAPK 的磷酸化，穩定 NF - κB/IκB 而阻止IκB 的脫落，抑制下游IL - 1β，IL - 6，TNF - α，NO/ iNOS 的表達，降低炎性因子的水平，同時能抑制COX - 2 的表達［33］。Cynan?dione A 通過調控 α7 - nAChR 依賴性 IL - 10/ β - 內啡肽信號通路抑制神經性疼痛［34］。隔山消積顆粒的活性成分通過作用于胃腸道炎性相關靶點IL-6，IL-1β，TNF - α，5 - 脂肪氧化酶，ESR2，ESR1，CA2，AR，HMGCR和ALOX5等，促進胃動力靶點和神經調控靶點ACHE等，通過調控炎性信號通道TLR4/MyD88/NF-κB、腦- 腸軸功能紊亂、胃腸動力障礙和精神心理因素等相關通路達到治療FD的作用。