点击上方蓝字,关注我们~

来源：《世界中医药》，网络首发日期：2021-12-03

基于网络药理学探讨黄芪治疗慢传输型便秘的作用机制

施敏 刘富林 夏旭婷 张婷

（湖南中医药大学，长沙，410208）

慢传输型便秘（Slow Transit Constipation，STC）是指因为结肠传输功能障碍而导致的肠内容物滞留或排出不顺利的状态。PrestonDM在1986年首次提出STC的概念，认为临床主要以便意减退或消失、便次减少，粪质干结、排便不畅并伴随腹部不适等症状为主要表现。相关流行病学研究显示，STC在慢性便秘中发病率占首位，其发病机制尚不明确，现有研究认为可能与肠神经递质、Cajal间质细胞、肠道菌群等多种因素有关。STC不仅严重影响患者的生活质量，且易导致肛肠、心脑血管等危险疾病发生。目前西医主要有药物和手术治疗,而西药有一定的依赖性和毒副作用,手术治疗存在风险性，传统中医治疗STC有其独特的优势作用。STC多以虚证为主,少有实证。《丹溪心法·脾约丸论》云:“胃强脾弱,约束津液不得四布……故曰脾约……脾弱不得通也。”故STC治疗重点在于健脾益气。黄芪有健脾益气的功效，现代药理学研究证明,黄芪有调节机体免疫功能的作用,尤其适用于体虚的老年便秘患者。单崇武等用黄芪汤加味治疗将50例肺脾气虚型功能性便秘患者，治疗组总有效率为88.0%,高于对照组50例的70.0%。

本研究基于系统生物学相关理论，借助网络药理学方法和分子对接技术，探讨黄芪的活性化合物与靶点之间的关系，探究黄芪多成分、多靶点、多通路治疗STC的作用机制，以期为黄芪治疗STC提供科学依据与临床应用参考。

1.1黄芪活性成分和靶蛋白筛选

应用中药系统药理分析平台（TCMSP，http://lsp.nwu.edu.cn/tcmsp.php）检索黄芪的活性成分，以类药性（DL）≥0.18和口服生物利用度（OB）≥30%进行筛选，得到黄芪的有效成分，并下载其分子结构，同时检索活性成分对应的靶蛋白。筛选完成后统一在Uniprot（https://www.uniprot.org/）蛋白质数据库、NCBI（https://www.ncbi.nlm.nih.gov/）等数据库输入靶点蛋白名称，限定物种为“Homo sapiens”将化合物作用的蛋白质靶点进行规范。

1.2STC疾病靶点筛选

在OMIM数据库（http://omim.org/）和Genecards数据库（https://www.genecards.org）中输入关键词“slow transit contipation”检索疾病靶点。去除重复靶点并经Uniprot校正整合得到STC相关靶点。运用R语言(VennDiagram包)将黄芪有效成分的靶蛋白与STC疾病靶点取交集，得到黄芪治疗STC的靶点。

1.3构建黄芪中药成分－疾病靶基因网络图

借助Cytoscape3.7.2软件构建中药－活性成分－疾病基因网络，并计算其相关网络特征值。

1.4蛋白质相互作用网络的构建及拓扑分析

将黄芪与STC交集靶基因导入String（https://string-db.org/）数据库，将物种设定为“Homo sapiens”，互相作用阈值设定为“median confidence”（0.4），其余均为默认设置，获取交互作用，构建蛋白互作（PPI）网络并计算其网络中每个靶基因degree,采用R语言绘制degree最大的前20个靶基因的柱状图。将得到的PPI网络导入CytoScape3.7.2，借助CytoNCA插件分析靶点的连接度（Degree）、介度（Betweenness）等网络拓扑参数，根据Degree和Betweenness≥2倍中位数筛选核心靶点。

1.5GO富集分析、KEGG通路分析

应用R语言（dose、cluster Profiler、enrichplot包）从生物过程(BP)、细胞成分(CC)、分子功能(MF)进行GO通路富集分析，同时进行KEGG通路富集分析，设定阈值P<0.05，绘制条形图和气泡图以进一步揭示黄芪治疗 STC 的作用机制。

1.6分子对接

在PDB数据库（http://www.rcsb.org/）下载核心靶蛋白的3D结构pdb格式文件，用Pymol软件移去水分子和配体分子，在TCMSP数据库中下载化合物3D结构的mol格式。导入SwissDock在线分子对接工具（http://www.swissdock.ch/）将靶点与黄芪活性成分进行对接，用以评价活性成分与关键靶点之间的对接效果。

2.1筛选黄芪的活性成分和靶蛋白

通过TCMSP检索并结合文献筛选到黄芪所含丁子香萜、熊竹素、常春藤皂苷元、β-谷甾醇、异鼠李素、3,9-二-O-甲基尼森香豌豆紫檀酚、5’-羟基异微凸剑叶莎醇-2’,5’-二-O-葡萄糖苷、7-O-甲基-异微凸剑叶莎醇tol、9,10-二甲氧基紫檀烷-3-O-β-D-葡萄糖苷、黄芪紫檀烷苷、联苯双酯、芒柄花素、二氢异黄酮、毛睿异黄酮、山柰酚、叶酸、7,2ʹ-二羟基-3ʹ,4ʹ-二甲氧基异黄烷、黄芪异黄烷苷、1,7-二羟基-3,9-二甲氧基紫檀烯、槲皮素共20个活性成分，筛选出作用靶基因953个，去重后获得相关靶基因189个。

2.2STC疾病靶点筛选

经Genecards和OMIM数据库检索，分别得到STC靶基因3278、143个，整合后共3817个，经过Uniprot校正整合后得到3174个，将其与2.1得到的药物作用靶基因189个借助R语言取交集得到PGR、NOS2、PTGS1、AR、SCN5A等共143个，交集韦恩图见图1。

2.3构建黄芪中药成分－疾病靶基因网络图

整合黄芪、活性成分和STC疾病靶基因之间的作用关系，采用Cytoscape3.7.2软件进行可视化，得到黄芪-有效成分-疾病靶基因网络图如图2。

2.4蛋白质相互作用网络的构建及拓扑分析

将疾病与中药取交集得到的基因进行蛋白互作分析，得到PPI网络图见图3。其中“边”代表作用靶基因之间的关联，粗细用combine score表示，“边”越粗，combine score值越大，代表其结合度越大；节点代表作用靶基因。共有143个节点、2353条边，平均节点度为32.9，平均局部聚类系数为0.619。

绘制degree最大的前20个靶基因的柱状图如图4，Degree值排名前10位的核心节点分别为AKT1、IL6、VEGFA、CASP3、JUN、MYC、MAPK8、EGF、EGFR、MAPK1。

拓扑分析筛选的关键节点绘制得到黄芪治疗STC的核心调控网络图见图5。degree值排在前10位的分别为JUN、VEGFA、AKT1、CASP3、IL6、MAPK1、EGF、MAPK8、MYC、EGFR。

2.5GO富集分析和KEGG通路分析

GO功能富集分析得到GO条目2173个，其中生物过程、细胞组成、分子功能各1943、77、153个。生物过程主要与对脂多糖的回应、对细菌起源的分子的反应、对抗生素的反应等信号通路有关；在细胞成分主要与薄膜筏、膜微区、膜区、突触后膜的组成部分等信号通路有关；在分子功能主要与核受体结合、类固醇激素受体活性、泛素样蛋白连接酶结合、组蛋白激酶活性、泛素蛋白连接酶结合等信号通路有关。表明黄芪可通过调节多个生物过程、细胞组成、分子功能治疗STC，BP、CC、MF富集分析结果详见表1。

KEGG富集分析共得到163条信号通路，可知黄芪可通过乙型肝炎、前列腺癌、IL-17信号通路、TNF信号通路等作用通路发挥对STC的治疗作用，KEGG通路富集分析结果前20个条形图和气泡图详见图6（条形图颜色代表富集程度，条形长度代表富集的基因数目，气泡图颜色代表富集程度，气泡大小代表富集的基因数目）。

2.6分子对接

利用SwissDock在线分子对接工具，对JUN、VEGFA、AKT1、CASP3、IL6五个靶蛋白和黄芪20种活性成分进行分子对接,分子对接结果显示叶酸、1,7-二羟基-3,9-二甲氧基紫檀素、山奈酚、槲皮素四种成分与5个靶点结合比较紧密，详细结果见表2,蛋白质与小分子结合越紧密,Full Fitness值越低，叶酸与IL6对接结果见图7。

表2分子对接结果

慢传输型便秘是由于大肠传导功能紊乱而导致的排便困难。祖国医学归属不大便、脾约等范畴。其病机主要为气虚大肠传道无力。本病属本虚标实，主要以脾气虚为主。《症因脉治》中指出气虚则大肠传道失司，肠失濡润致糟粕传送失常，形成腑气不通，而致便秘。黄芪有补气健脾、升阳举陷、利水消肿的功效，临床上用于治疗肺脾气虚证等。有研究表明，黄芪可明显升高脾虚小鼠Hb值、淋巴因子(LAK)细胞活性、白介素-2(IL-2)活性及CD4/CD8细胞比值，并对其免疫功能低下有改善作用。余克强等通过对26个研究，1955例患者进行meta分析证明以黄芪为主的中成药或方剂治疗气虚型功能性便秘安全有效。

本研究从化合物-靶点网络图可知，黄芪中山柰酚、叶酸、槲皮素、异鼠李素、芒柄花素等20个化合物作用于STC疾病相关的143个靶点，说明了黄芪是通过多成分、多靶点协同发挥治疗STC的作用的。String在线数据库所得PPI网络显示节点度排名前的核心节点有AKT1、IL6、VEGFA、CASP3、JUN。核心基因调控网络图拓扑分析结果中degree值靠前的有JUN、VEGFA、AKT1、CASP3、IL6。二者取交集得到黄芪治疗STC的关键核心靶点AKT1、JUN、IL6、VEGFA、CASP3。

AKT1是PI3K/Akt信号通路下游的主要靶点，主要参与调节炎症反应。张高松通过实验研究推导出AKT在STC大鼠表达的减少可能是STC发病机制之一。现有研究认为JUN主要与癌症的发生有关，JUN癌基因激活后其蛋白以JUN或JUN同源二聚体的形式结合到转录激活蛋白1启动子结合位点,诱导癌基因表达,参与细胞恶性转化。IL6具有广泛的促炎症作用，其过度表达会造成肠道炎性病变，并可导致胃肠运动及分泌功能紊乱，其作用机制可能与激活胃肠道黏膜下层的神经元影响肾上腺应激反应和胃肠蠕动有关。VEGFA基因的表达可使肠道细胞间的紧密连接松弛，损伤肠黏膜屏障从而导致肠道微生物入侵肠壁。李超报道CASP3是细胞凋亡途径中的关键角色，肠细胞的非正常调亡可能会引起肠黏膜屏障破坏受损，进一步造成STC的发生。GO富集分析预测出黄芪对STC的治疗机制可能是通过对脂多糖的回应、对细菌起源的分子的反应、薄膜筏、膜微区、丝氨酸/苏氨酸蛋白激酶复合物、核受体结合、类固醇激素受体活性产生作用的。

KEGG富集分析预测出黄芪对STC的治疗机制可能是通过流体剪切应力与动脉粥样硬化、乙型肝炎、IL-17信号通路产生作用的，结合黄芪主要活性成分和其作用靶点，我们推测黄芪可能通过抗炎、提高免疫、保护肠黏膜屏障、调节肠道菌群和神经递质等作用机制，从而促进肠道运动，达到益气健脾、通阳利水的功效，改善STC症状。

本研究得到的相关通路虽然与目前研究较多的胶质细胞源性神经营养因子（GDNF）/GFRα1/神经细胞黏附因子（NCAM）、SCF/c-kit、5-羟色胺（5-HT）、PKA和MAPK、PLC[1]γ1/PLC-γ2等典型的治疗STC信号通路不大相同，但正说明了黄芪治疗STC是以多靶点、多通路的形式来发挥作用的，也恰恰体现了黄芪多功效治疗STC的特点。黄芪有效成分与STC疾病靶点分子对接结果显示山柰酚、叶酸、槲皮素、3,9-二-O-甲基尼森香豌豆紫檀酚与JUN、VEGFA、AKT1、CASP3、IL6核心靶点结合较好。相关研究显示槲皮素能可抑制肠炎小鼠血清中髓过氧化物酶、丙二醛和一氧化氮浓度水平升高，减轻细胞氧化损伤，从而改善肠道微生物群失调。汤利华等研究发现，山柰酚可通过抑制DN大鼠NLRP3炎症小体的激活，从而抑制炎症因子的表达。叶酸可增强机体免疫能力，因此可用于免疫功能低的便秘患者治疗，但目前尚无公开发表的叶酸与STC相关的研究文献。关于3,9-二-O-甲基尼森香豌豆紫檀酚与慢传输型便秘之间作用机制的目前尚无相关文献报道，待后期进一步实验验证。

黄芪可能通过抗炎、提高免疫、保护肠黏膜屏障、调节肠道菌群和神经递质，从而促进肠道运动等过程改善STC症状；通过调控对脂多糖的回应、对细菌起源的分子的反应、薄膜筏、膜微区、丝氨酸/苏氨酸蛋白激酶复合物、核受体结合、类固醇激素受体活性、乙型肝炎、IL-17信号通路等信号通路发挥治疗STC作用。本研究借助网络药理学方法与分子对接技术，探讨了黄芪活性成分、靶点、通路之间相互作用的关系，发现黄芪能通过多靶点、多通路发挥治疗STC的作用，证明了黄芪是通过多途径发挥行气健脾，通阳利水的功效，达到治疗STC的目的，为黄芪用于治疗STC的进一步实验验证及临床研究提供了思路与理论基础。