从TCMSP（https://tcmspw.com/tcmsp.php）^[11]、TCMID（http://www.megabionet.org/tcmid/）、ETCM（http://www.tcmip.cn/ETCM/）^[12]和BATMAN（http://bionet.ncpsb.org.cn/batman-tcm/）^[13]等数据库中获取冰片、蟾酥、人工牛黄、人参提取物、肉桂、人工麝香、苏合香等七味药材所含化合物及其对应靶点。吸收、分布、代谢、排泄（ADME）指标是筛选活性化合物的重要指标，将全方所含化合物在TCMSP数据库中进行ADME筛选，将OB≥20%且DL≥0.1的化合物视为有效活性成分^[14]，对不符合ADME筛选的化合物进行检视，将已被证明为生物活性成分的化合物进行回收，去除重复项，即得麝香保心丸全方有效成分及其潜在靶点。为了提高获取数据的可信度，将所获得的靶点信息进行进一步筛选，剔除BATMAN数据库中score低于48的靶点和TCMID中score低于400的靶点，并去除仅有一个成分作用的靶点，得到全方靶点信息。

从GeneCards（https://www.genecards.org/）^[15]、OMIM（https://omim.org/）^[16]、TCMSP（https://tcmspw.com/tcmsp.php）、DrugBank（https://go.drugbank.com/）和DisGeNet（https://www.disgenet.org/search）^[17]等数据库中获取动脉粥样硬化疾病靶点。剔除DisGeNet中靶点分值低于0.02的靶点和GeneCards中分值低于1的靶点，合并各数据库靶点，得到AS相关的靶点蛋白。

将获取的麝香保心丸潜在靶点和AS相关靶点导入STRING（https://string-db.org/）数据库，统一转换为基因ID。取两者交集，得到麝香保心丸治疗AS相关靶点及其对应的化合物，将获取的数据导入Cytoscape3.8.2软件中进行进一步分析，构建化合物-靶点网络，以便于分析网络中的关键节点。

通过STRING数据库对麝香保心丸治疗AS的靶点进行分析，构建蛋白-蛋白相互作用（PPI）网络，并将PPI网络导入CytoScape3.8.2软件进一步进行分析，探究麝香保心丸治疗AS靶点之间的相互作用关系。

使用Metascape（http://metascape.org）^[18]数据库对获取的麝香保心丸抗AS靶点蛋白进行GO富集分析和KEGG通路分析，设定P<0.01，取生物学过程（GOBP）、分子功能（GOMF）和细胞成分（GOCC）条目的前10项，取KEGG富集通路的前20项，分析麝香保心丸治疗AS的主要相关通路及其功能。

各个数据库中共检索到麝香保心丸中474个化合物，对应的潜在靶点共有6158个。经筛选后共得到麝香保心丸中化合物124个，对应的潜在靶点为452个。

结合GeneCards、OMIM、TCMSP、DrugBank和DisGeNet等数据库，经过筛选，得到AS相关靶点1786个。

取麝香保心丸潜在靶点与AS相关靶点绘制维恩图（图1），取二者交集，即得麝香保心丸治疗AS的175个靶点，相关成分114个。采用STRING数据库对获取的175个靶点进行分析，得到一个具有175个点和2303个边的PPI网络（图2），通过对该网络图进行分析，根据各靶点在网络中的度值判断该靶点在网络中的重要程度。在麝香保心丸治疗AS的过程中，ALB、INS、ACTB、AKT1、IL-6、TNF、IL-1B、PPARG、MAPK3、CREB1等靶点的度值排在前10位，表明这些靶点作为关键节点起到最主要的作用。

图  1  麝香保心丸与动脉粥样硬化靶点维恩图

图  2  麝香保心丸治疗动脉粥样硬化蛋白之间相互作用网络

用CytoScape3.8.2软件构建麝香保心丸治疗AS的化合物及其对应的靶点网络（图3），得到具有289个节点和1166条边的复合网络，图中化合物节点为红色，靶点节点为绿色。依据PPI网络中各靶点度值对各化合物进行评分，每个化合物分值为其潜在作用靶点的度值之和，以分值高低为依据评判麝香保心丸治疗AS过程中的贡献度，筛选出关键的20个化合物（表1）。鹅去氧胆酸、人参皂苷R_b1、肉桂醛、胆酸、熊去氧胆酸等化合物作为关键化合物，分值为1486、925、827、809和805，分别作用63、17、42、49和31个靶点，体现了中药的多成分、多靶点协同作用特点。

图  3  化合物-靶点网络

对GO富集条目进行分析（图4），麝香保心丸调控的靶点在质膜外侧、细胞质囊泡腔、内质网内腔、细胞顶端部分等区域富集，影响核受体活性、氧化还原酶活性、蛋白质结构域特异性结合、一氧化氮合酶调节活性等分子功能，调节血液循环、细胞对脂质的反应、细胞对激素刺激的反应、对脂多糖的反应等生物过程。KEGG通路富集显示（图5），麝香保心丸主要调节流体剪切应力与动脉粥样硬化、神经活性配体-受体相互作用、钙离子信号通路、胆汁分泌等通路，实现对AS的治疗效果。

图  4  GO富集分析

图  5  KEGG通路富集

本研究基于网络药理学的方法探讨麝香保心丸治疗AS相关化合物及其潜在靶点和通路。本研究发现鹅去氧胆酸、人参皂苷Rb1、肉桂醛、胆酸、熊去氧胆酸可能在治疗AS时起到关键的作用。鹅去氧胆酸能够降低机体炎症反应、提高机体的抗氧化能力和缓解细胞自噬^[19]。熊去氧胆酸能够促进AS模型ApoE(-/-)小鼠的胆固醇流出，降低炎症反应水平，减少AS斑块面积^[20]。人参皂苷Rb1能够减轻炎症反应和氧化应激，并通过抑制细胞凋亡和促进自噬达到抑制ApoE(-/-)小鼠的早期AS的发展^[21]。肉桂醛^[22]能够调节IκB/NF-κB通路从而抑制炎症和氧化达到对AS的治疗效果。同时，这些关键成分也是麝香保心丸的入血成分，从侧面验证了网络药理学分析麝香保心丸治疗AS机制的可靠性^[23]。

对麝香保心丸治疗AS的261个靶点进行分析，得到最主要的10个靶点。有研究证实，ALB与AS的严重程度具有显著的相关性^[24]。AKT1参与代谢、细胞生长与存活、血管新生等多个过程进而防治AS^[25]。INS是调节机体能量代谢的关键蛋白，胰岛素抵抗是导致AS的关键危险因素之一^[26]。TNF能够调控细胞生长与凋亡^[27]，且TNF能与IL-6和IL-1B协同作用诱导VEGF的生成，促进血管新生^[28]。PPARG是降血脂药的关键受体之一，能够调节脂质代谢和抑制炎症^[29]，对AS的防治起到重要的作用。这些关键靶点从各方面对AS进行治疗，且相互之间存在协调作用，达到了更好的治疗效果。

麝香保心丸治疗AS的靶点主要通过调节内分泌系统和信号转导系统，对炎症、免疫和氧化应激等过程进行调控，最终达到对AS的治疗作用。血管内皮剪切应力是心血管系统中的关键调控因子，低血管内皮剪切应力会诱发炎症并且使内皮细胞更容易被低密度脂蛋白和高血糖等危险因素影响，促进AS的发展^[30]。钙离子作为第二信使在体内发挥重要作用，能够介导信号转导并调节生理功能，且钙离子能够促进AS斑块的形成，研究证实AS斑块中的钙离子浓度显著增加^[31]。胆汁酸代谢是调节机体胆固醇水平的关键途径，且胆固醇还能通过与核受体协同作用进一步调节免疫和炎症反应，在AS的发生发展过程中发挥重要作用^[32]。对这些通路进行分析发现麝香保心丸可能通过调节内分泌，调控炎症、免疫、脂质代谢，保护血管内皮细胞等作用达到对AS的治疗效果。

本研究通过网络药理学的方法全方面、多角度对麝香保心丸治疗动脉粥样硬化整体协同作用机制进行研究，发现了治疗过程中可能发挥主要作用的化合物、靶点和通路。本研究为麝香保心丸抗AS的机制研究提供了思路，为临床长期应用提供理论依据。