缺血性卒中(ischemic stroke，IS)是由高血压、糖尿病、高脂血症及体内高凝状态等多种危险因素导致脑血流量降低、脑部缺血缺氧引发神经功能缺损的一种疾病^{[1, 2]}。根据中国2019年流行病学调查显示，缺血性卒中的发病率、患病率、复发率、死亡率逐年攀升，总体发病风险仍位居世界首位^[3]，这种持续上升的形势依然严峻。健脾补肾活血方(Jianpi Bushen Huoxue Prescription，JPBSHXP)由河南中医药大学第一附属医院脑病专家团队拟定，已在临床应用十余年。该方针对缺血性卒中本虚标实的基本病机，以健脾补肾，活血通络为治疗原则。方中黄芪、党参、白术健脾益气；牛膝、杜仲、桑寄生、肉苁蓉补肾固本；丹参，川芎活血通络；石菖蒲化痰开窍；炙甘草调和诸药。前期研究表明，JPBSHXP可提高缺血性卒中患者脑缺血区灌注水平，减轻神经功能缺损，改善中医证候、临床症状和体征，提高日常生活能力^[4]，但其作用机制尚不明确。网络药理学是从分子层次阐释中药方剂作用原理的有效方法，本文旨在应用网络药理学的分析方法，通过靶点的筛选、通路分析以及分子对接技术，为验证该方治疗IS的有效性提供合理依据，同时为此方进一步研究提供新的方向。

1 材料与方法 1.1 材料

本文所需数据库及相关应用程序，如Tab 1。

1.2 方法 1.2.1 健脾补肾活血方化学成分及靶点的获取

利用TCMSP中药系统药理数据库及分析平台检索健脾补肾活血方11味中药的有效活性成分及对应靶点信息，根据药代动力学参数吸收、分布、代谢和排泄(ADME)的原则，以口服生物利用度(oral bioavailability，OB)≥30%，类药性(Drug-like，DL)≥0.18为条件筛选潜在化合物，OB、DL越高，说明药物在体内发挥效能的作用越强^[5]。使用Perl语言将11味中药整理合并，得到健脾补肾活血方的中药-活性成分-靶点的总体信息表。

1.2.2 收集疾病靶点及预测药物作用靶点

通过Genecards、OMIM、DrugBank、TTD、PharmGkb五种疾病数据库，以“Ischemic stroke”为关键词，查找IS的疾病基因靶点，用R语言程序将其合并；Genecards数据库获取基因分值越高与疾病相关性越强，选择分值大于1与IS相关性高的靶点进行研究；在Uniprot平台将所有关于人的及被证实的基因下载，并通过perl语言将药物靶点转换为基因ID编码的形式；取药物和疾病的共同靶点并绘制venn图，预测药物靶点。

1.2.3 中药复方网络的构建

设计perl语言程序对健脾补肾活血方与缺血性卒中相同靶点进行整理，应用Cytoscape 3.8.0软件构建JPBSHXP-活性成分-靶点和IS相互作用的调控信息网络图，应用拓扑插件计算介数中心性(betweenness centrality，BC)，紧密中心性(closeness centrality，CC)，度值(degree centrality，DC)，特征向量中心性(eigenvector centrality，EC)，局部平均连通度中心性(local average connectivity centrality，LAC)，通过DC、BC筛选该方核心活性成分。

1.2.4 蛋白互作(PPI)网络的构建

通过String平台构建PPI蛋白互作网络，将药物与疾病交集所得到的190个基因导入List Of Names，设置highest confidence为0.9，去除游离靶点后得到170个靶点，应用Cytoscape3.8.0进行PPI网络分析，拓扑参数取大于BC、CC、DC、EC、LAC的中位值得到核心靶点。

1.2.5 GO功能分析和KEGG信号通路富集分析

应用David数据库进行GO分析及KEGG信号通路富集分析，根据P值提取前20名的生物学过程(biological process，BP)、细胞组成(cellular components，CC)、分子功能(molecular function，MF)及信号通路的数据，运用微生信在线制图进行气泡图的绘制。

1.2.6 分子对接

使用分子对接技术研究JPBSHXP的活性成分与治疗IS的相关靶点。首先，从pubchem数据库下载化合物结构，保存为SDF格式，并通过Chem 3D软件将SDF格式转化为mol2格式，从RCSB数据库下载相关蛋白，使用Pymol软件去除溶剂分子与配体，使用AutoDock Tools 1.5.6软件计算电荷、加氢等。最后运行Autodock vina 1.1.2进行分子对接, 采用Discovery Studio 2020可视化分析对接构象。

2 结果 2.1 成分及靶点的筛选

根据TCMSP得到JPBSHXP中药-活性成分-靶点(单位：个)，其中黄芪-17个-462个、党参-17个-214个、白术-4个-23个、杜仲-26个-519个、桑寄生-2个-157个、怀牛膝-17个-444个、肉苁蓉-6个-223个、丹参-59个-932个、川芎-6个-42个、石菖蒲-4个-103个、甘草87个-1 761个。共有活性成分245个，靶点4 880个，剔重后取得该方有效成分201个，靶点263个。

2.2 疾病靶点及预测药物作用靶点

通过R语言设计程序将疾病基因合并、剔除重复项并取唯一值，得到2 808个疾病靶点。药物靶点与疾病靶点相交集得到共同基因190个。(详见Fig 1)。

2.3 健脾补肾活血方网络构建及拓扑分析

通过Sytoscape 3.8.0构建中药-活性成分-靶点-疾病网络图。Network Analyzer计算分析得到398个节点，2 118条边。拓扑研究显示健脾类药物主要化学成分槲皮素(quercetin)、山柰酚(kaempferol)、木犀草素(luteolin)等。其中槲皮素、山柰酚、豆甾醇为健脾类、补肾类共有活性成分，其余补肾类药物显著化学成分有汉黄芩素(wogonin)、β-谷甾醇(beta-sitosterol)等。活血类药物为丹参、川芎，研究显示发挥活血作用的药物为丹参，主要有丹参酮IIA(tanshinone iia)、隐丹参酮(cryptotanshinone)等化合物(Tab 2)。

2.4 PPI网络构建及分析

将190个交集基因导入String数据库，置信度设置为0.9，并剔除游离靶点，共得到170个基因，将170个基因通过sytoscape 3.8.0进行网络可视化分析和拓扑分析。根据拓扑分析结果，取BC>58.4、DC>7、CC>0.11、EC>0.03、LAC>3(中位值)的靶点作为该方的核心基因。STAT3、Tp53、MAPK1、MAPK3、JUN、APP、MAPK14、MAPK8、AKT1、EGFR、VEGFA等基因可能是该方防治缺血性卒中的潜在靶点(Fig 2)。

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

GO分析共得到880个分析条目，其中生物学过程有662条，最多对应交集基因有40个，最少2个；细胞组成有82条，最多对应交集基因有87个，最少2个；分子功能有136条；最多对应交集基因有150个，最少2个。根据排名前20的GO条目分析显示，健脾补肾活血方治疗缺血性卒中可能与对药物的反应、对缺氧的反应、对乙醇的反应、对脂多糖的反应、对雌二醇的反应、老化等生物学过程密切相关，分子生物学功能集中在酶结合、蛋白质均二聚活性、血红素结合、类固醇结合、转录因子结合、多巴胺结合、蛋白激酶结合、多巴胺神经递质受体活性等方面，细胞组分主要有细胞外空间、膜筏、质膜等。

KEGG富集分析共得到125条通路，主要与癌症通路、乙型肝炎、弓形虫病、大肠癌、TNF信号通路、膀胱癌、胰腺癌、蛋白聚糖在癌症中、HIF-1信号通路、前列腺癌、局灶性粘连、丙型肝炎、PI3K-Akt信号通路、细胞凋亡、Toll样受体信号通路等密切相关(Fig 3)。

2.6 分子对接

Quercetin、Tanshinone IIa、Luteolin、Kaempferol与STAT3、Tp53、MAPK3、MAPK1、JUN、APP、AKT1关键靶蛋白进行28次对接，结果显示靶点与化合物均有较强的结合能(Tab 3)。Quercetin主要是健脾补肾药物共有活性成分，β淀粉样前体蛋白(APP)与Quercetin结合较强，Tanshinone IIa为活血药物主要成分，主要与转录因子AP-1(JUN)、丝裂原活化蛋白激酶3(MAPK3)、丝裂原活化蛋白激酶1(MAPK1)等化合物结合具有较好的稳定性(Fig 4)。另外，氢键使结构更具稳定性。例如，APP与Quercetin的a环Arg381、Pro385形成氢键作用；b环与Arg500、Tyr369形成氢键作用；c环与Ala334、Glu362、Tyr501形成氢键作用。

3 讨论

课题组成员前期多角度、多层次阐释健脾补肾活血方治疗缺血性卒中的药理机制。实验研究表明, 该方通过抑制CD3⁺、CD8⁺和调节CD4⁺淋巴细胞在缺血/再灌注损伤后脑组织内的浸润，减少免疫应答所介导的神经细胞凋亡，增加大鼠外周血中内皮组细胞(EPCs)的含量，增强脑组织CD31抗原的表达，促进大鼠脑组织新生血管的形成、有效建立侧支循环，从而减轻脑损伤。临床试验也显示，该方明显增加患者外周血EPCs含量，提高缺血区血液灌注水平、改善细胞免疫、拮抗炎症反应，减少卒中后氧化应激反应，增强机体抗氧化应激能力，同时修复受损的血管内皮细胞^{[4, 6-8]}。

拓扑分析得出槲皮素、木犀草素、山柰酚、丹参酮IIA可能是JPBSHXP治疗IS的关键化学成分。Oh等^[9]将槲皮素应用于缺血/再灌注诱导小鼠模型中发现，抑制在细胞活化期间由糖蛋白VI介导的血小板信号转导，明显破坏胶原相关肽诱导的血小板聚集，减少颗粒分泌及活性氧生成，降低细胞内钙动员水平，抑制动脉血栓形成。木犀草素通过抑制基质金属蛋白酶(MMP)9和激活PI3K/Akt信号通路增强细胞活力，并下调细胞凋亡，减少了脑梗塞体积^[10]。山柰酚促进血管内皮生长因子(VEGF)诱导血管内皮生长因子受体2(VEGFR2)磷酸化、内皮一氧化氮合酶(eNOS)和细胞外信号调节激酶(Erk)表达，从而促使内皮细胞增殖、迁移和血管新生^[11]。徐璇实验得出丹参酮IIA可使脑缺血大鼠的PTEN表达增强，降低自噬通路PI3K/AKT/mTOR的表达水平，减轻脑缺血后炎症反应，保护大鼠的脑组织^[12]。该方治疗IS的核心靶点有STAT3、TP53、MAPK1、MAPK3、APP、JUN、AKT1等。转录激活因子3(STAT3)广泛存在于中枢神经系统中，参与细胞生长、代谢、凋亡，血管生成和氧化应激等多种生物学途径，同时介导机体免疫失调、炎症反应等过程^[13]。余录等^[14]研究发现山柰酚糖苷可能通过STAT3信号通路减少卒中后TNF-α、IL-1β炎性因子，降低白质区胼胝体和纹状体的β淀粉样前体蛋白(APP)阳性表达，改善大鼠神经行为障碍，缩小脑梗死范围，减轻神经元损伤。

健脾补肾活血方作用于缺血性卒中的炎症通路有肿瘤坏死因子(tumor necrosis factor，TNF)信号通路、Toll样受体(Toll-like receptors，TLRs)信号通路；促血管新生通路主要是低氧诱导因子-1(hypoxia inducible factor-1，HIF-1α)信号通路、磷脂酰肌醇-3-激酶/丝苏氨酸蛋白激酶(phosphatei-dylinositol-3-kinase/serine-threonine kinase，PI3K /Akt)通路(Tab 4，Fig 5)。脑梗死急性期常出现大量的炎症因子，TNF信号通路富集基因显示可能调节IL-1β及CCL2、CXCL10、NFKBIA等致炎因子来对抗脑缺血。Toll样受体信号通路参与神经炎症的反应过程^[15]，可能是通过MAPK信号通路，调控转录因子AP-1，调节下游基因IL-8、IL-1β、NFKBIA、IKBKB、CXCL10、CXCL11等炎性因子参与中风后炎性反应。这也为未来临床及实验研究炎症信号通路方向提供新的研究指标。PI3K/Akt通路，参与多种通路的调节，与炎症免疫反应、细胞凋亡、HIF-1α等通路均有相关性^[16-17]。PI3K-Akt通路得出该方可能是调控GF、RTK、ECM等上游蛋白，使PI3K激活，AKT磷酸化启动下游PKCs、eNOS、ERK、GSK3、p21、Bcl-2、Bax、caspase-9等靶蛋白来作用于缺血性卒中。以上得知健脾补肾活血方可能通过免疫、炎症反应、氧化应激、血管新生、细胞凋亡等途径作用于缺血性卒中。

网络药理学通过生物信息学、计算机技术、药理学等对方剂进行现代化阐释，为新药创制研发提供思路，也给予关键的技术支撑，对后续实验研究指明有希望的验证方向。本研究通过该方法得出缺血性卒中的潜在化学成分、小分子蛋白以及信号通路，在分子生物学层面为该方进一步实验研究及临床应用提供了理论基础。中医药在传承中发展，在发展中创新，中药新药以证候为干预对象，以辨证论治为理论基础，具有异病同治的特征。临床研究中本方对脾肾两虚型脑梗死、卒中后认知障碍和脑白质病变具有明显疗效，是否通过调控这些通路、靶蛋白发挥作用，有待于课题组下一步实验和临床验证。