基于UPLC-Q-TOF-MS/MS结合网络药理学及分子对接研究益母草水提物治疗动脉粥样硬化的作用机制*

李梓宁 ，靳雅 ，李霄 ，张晗 ，张鹏

（1.天津中医药大学组分中药国家重点实验室，天津 301617；2.方剂学教育部重点实验室，天津 301617；3.天津中医药大学中医药研究院，天津 301617）

动脉粥样硬化（AS）是一种累及大中动脉的血管壁慢性炎症性疾病，其发病机制复杂，受遗传因素、环境因素、生活方式等多方面影响[1]。AS是心脑血管疾病与外周血管疾病的共同病理基础，常导致缺血性心脏病、缺血性中风等疾病的发生。随着西医的不断发展，以他汀类、抗血小板为主的药物治疗，以及以经皮冠状动脉介入治疗为主的外科介入治疗手段，在一定程度上减少了AS引发的心脑血管疾病的病死率。然而，他汀类药物具有一定的引发肌损伤及肝肾损伤的风险，且在他汀的降胆固醇治疗达标后，患者仍存在明显的心血管残留风险。此外，他汀单药不能满足极高危患者的治疗需求[2-4]。外科介入手术虽可临时缓解冠状动脉管腔狭窄及堵塞，但常易在术后几月或几年内出现“支架内再狭窄”情况，并未从根本上解决疾病状况，且患者术后需终身服药，增加医疗成本与负担[5]。因此，寻找可靠、有效、不良反应小的天然药物来防治AS疾病，一直以来被寄予厚望。

较之单一化学成分或单一作用靶点的药物而言，中药具有多成分、多靶点的优势，在AS的临床治疗方面已取得一定的效果[6]。益母草为唇形科植物益母草Leonurus japonicas Houtt.的新鲜或干燥地上部分，最早收录于《神农本草经》，被列为上品。其味辛、微苦，性微寒，入心包、肝经，具有活血化瘀、利尿消肿、清热解毒之效[7]。中医认为，AS疾病多以“痰、瘀、浊、毒、虚”等为主要病理因素，其“瘀”不仅贯穿其发生、发展始终，又常与痰、浊等相互交结加速疾病进程[8]，是临床防治AS的关键之一。益母草虽为治疗月经不调、痛经经闭、产后恶露不尽的妇科良药，但长于活血养血，《本草汇言》云：“益母草，行血养血，性而不伤心血，养血而不滞血，诚为血家之圣药也”[9]。因此，益母草目前已常用于心脑血管疾病的防治，尤其是复方的组方配伍之中，如冠心康[10]、醒脑治瘫胶囊[11]等，取其活血养血、通利脉道之功。现代药理研究表明，益母草化学成分主要有萜类、生物碱、黄酮类和苯丙素类等，其中萜类成分显示出抗血栓、血管舒张、抗血小板聚集、抗炎等生物活性作用；生物碱类成分不仅具有抗炎、抗氧化活性，还可调节血管生成过程，并对子宫平滑肌具有双向收缩调节作用[12]。药效学研究发现，益母草能通过降低低密度脂蛋白胆固醇、总胆固醇和三酰甘油延缓AS的发展，还可通过减少氧化应激和减少炎症因子的表达，对肿瘤坏死因子-α（TNF-α）诱导的人脐静脉内皮细胞炎症起到保护作用，但缺乏对其机制的研究[13-15]。

网络药理学采用现代化的研究方法，将医学与计算机科学结合在一起，基于“多基因、多靶点、多途径”的互作网络，为中药治疗各种疾病的作用靶点和潜在作用机制提供了预测价值[16]。为了更好地探索益母草活性成分治疗AS的潜在分子机制，本研究借助超高压液相色谱串联四级杆飞行时间质谱联用（UPLC-Q-TOF-MS/MS）技术、网络药理学和分子对接的方法，从化学成分、作用靶点和信号传导通路多角度对益母草水提物干预AS的作用机制进行系统研究，以期为进一步揭示益母草的药理作用机制提供理论基础。

1 仪器与材料

UPLC-Q-TOF-MS/MS仪UPLC Xevo G2-XS QTOF（Waters公司），ACQUITY UPLC HSS T3液相色谱柱（2.1 mm×100 mm，1.8 μm），高速冷冻离心机（赛默飞世尔科技有限公司，货号：75002440），电热套（北京中兴伟业有限公司，批号：2D4W），旋转蒸发仪（上海禾汽玻璃仪器有限公司，批号：RE5005BEX），电子天平（赛多利斯科学仪器有限公司，序列号：3137811120），色谱级甲酸（ASC公司），乙腈（Fisher公司），超纯水（屈臣氏），益母草药材2021年购于北京同仁堂，经天津中医药大学李天祥教授鉴定为唇形科植物益母草Leonurus japonicas Houtt.地上部分。

2 方法

2.1 益母草水提物的制备及处理 称取适量益母草饮片加入蒸馏水回流提取3次，每次加入蒸馏水用量分别为10倍量、8倍量、8倍量，分别回流提取30 min、20 min、20 min，煮沸后计时，合并 3 次滤液，旋蒸，减压浓缩冻干。

2.2 益母草水提物的UPLC-Q-TOF-MS/MS分析 鉴定称取益母草水提物药粉适量，采用50%甲醇超声提取20 min，提取液在转速13 000 r/min，离心半径8.6 cm，温度20℃下离心15 min，取上清液进样。液相色谱条件流速0.4 mL/min；柱温35℃；进样量10 μL；流动相为0.1%甲酸水溶液（A）-乙腈（B），梯度洗脱。洗脱时间见表1。

表1 梯度洗脱时间Tab.1 Gradient elution time

2.3 质谱条件 采用电喷雾离子源（ESI），质谱扫描模式采用正、负离子2种扫描模式，毛细管电压为ESI+3.0 kV/ESI-2.5 kV，锥孔电压为40 V。离子源温度110℃，雾化器温度为450℃，雾化器流速为800 L/h，锥孔气流量为50 L/h。MSE Continuum扫描模式检测，质量扫描范围为m/z 100～1 200，碰撞气为氩气，碰撞电压为 15～30 V、20～40 V、30～50 V，用亮氨酸脑啡肽LE（m/z 556.2677[M＋H]+）作为内标进行质量实时校正，体积流量为0.5 μL/min，进样量10 μL，数据采集软件为MassLynx 4.2，数据处理软件为UNIFI 1.9.4。

2.4 药物活性成分筛选及作用 靶点预测利用有机小分子生物活性数据库（PubChem数据库，https：//pubchem.ncbi.nlm.nih.gov）收集益母草水提物各成分化学结构。借助SwissADME数据库（http：//www.swissadme.ch），根据益母草水提物各成分的化学结构或SMILES字符，预测其成分的ADME参数、药代动力学特性、类药物性质等，采用胃肠道吸收（GI-absorption=high）、 类 药 性 Druglikeness 中Lipinski、Ghose、Veber、Egan 和 Muegge 为 yes且≥3为筛选标准，确定益母草水提物的主要活性成分。利用 SwissTarget Prediction数据平台（http：//www.swisstargetprediction.ch/），预测益母草水提物主要活性成分的潜在作用靶点。

2.5 疾病潜在作用靶点筛选及药物-疾病作用 靶点获取通过 GeneCards（https：//www.genecards.org/）、OMIM（https：//www.ncbi.nlm.nih.gov/omim）和 DisGeNET（https：//www.disgenet.org/search） 等数据库挖掘与AS相关的靶基因；利用Excel对收集到的疾病靶点进行去重筛选，得到疾病相关的作用靶点。将AS相关的作用靶点映射至益母草水提物活性成分的潜在作用靶点上，借助Draw Venn Diagram在线绘图工具绘制益母草活性成分潜在作用靶点与AS疾病靶点的Venn图并得到交集靶点。

2.6 药物防治疾病的调控网络构建将AS相关的作用 靶点映射至益母草水提物活性成分的潜在作用靶点上，得到成分-疾病靶点文件。借助Cytoscape 3.7.1软件进行网络可视化，构建活性成分-交集靶点调控网络。

2.7 靶点蛋白互作网络（PPI）构建及核心靶点选筛 使用STRING数据库构建益母草抗AS的PPI网络图，生物体设置为“Homo”，得到PPI.tsv数据，借助Cytoscape3.7.1软件绘制并构建“共同作用靶点”的关系网络，PPI网络由代表目标蛋白质的节点和代表蛋白质-蛋白质相互作用的边组成。边缘的粗细代表组合分数。度是指直接连接到一个节点的其他节点的数量。度数越高，节点越重要。同时利用Cytoscape 3.7.1软件中 CytoNCA和 BisoGenet对PPI网络进行拓扑分析，以介度中心性（BC）、对接中心性（CC）、度自由性（DC）和Degree值等为筛选标准，筛选出核心靶点。

2.8 GO和京都基因和基因组百科全书（KEGG）通路富集分析使用 DAVID在线分析软件将交集靶点进行GO生物过程和KEGG通路富集分析，得到通路富集文件，以P≤0.05为标准筛选通路，使用在线工具微生信进行可视化，得到GO柱状图和KEGG气泡图。

2.9 分子对接分析 将益母草水提物有效成分与核心蛋白进行分子对接，从Pubchem数据库（https：//pub-chem.ncbi.nlm.nih.gov/）和 Chemspider数据库（http：//www.8chemspider.com/）下载益母草水提物活性成分的SDF格式文件，使用Chem3D软件转化为Mol2格式文件。从Protein DataBank数据库（PDB，https：//www.rcsb.org/）获取 PPI网络中核心靶点的3D结构的PBD格式文件[白细胞介素（IL）-6（PDB：1ALU）；肿瘤坏死因子（TNF，PDB：2E7A）；IL-1β（PDB：1LZH）；血管内皮生长因子 A（VEGFA）（PDB：1MKK）；CASP3（PDB：1NMS）；myc（PDB：1NKP）；表皮生长因子受体（EGFR，PDB：1XKK）；SRC（PDB：1O41）；丝裂原活化蛋白激酶 3（MAPK3）（PDB：4QTB）；HSP90AA1（PDB：1OSF）]，导入 PyMol软件去除水分子和异质分子。用Auto Dock Tools分子对接模拟软件对蛋白质进行添加氢原子等操作，再将化学成分和靶点文件均转化为pdbqt格式的文件，运用Grid模块设置蛋白原配体为对接盒子中心，利用Auto Dock Vina软件将受体（靶点）与配体（化学成分）进行半柔性对接，结合能小于0 k J/mol表明受体和配体具有自发结合潜力，结合能小于等于-5.0 kJ/mol说明化学成分与靶点对接效果较好[17]。

3 结果与分析

3.1 药物活性成分筛选及作用 靶点预测益母草水提物UPLC-Q-TOF-MS/MS鉴定结果得到46个成分，具体信息见表2，正离子及负离子模式图见图1。借助SwissADME数据库，通过ADME筛选标准，得到7个活性成分，分别是益母草碱、益母草叟苷、Leoheteronin C、汉黄芩素、葫芦巴碱、N，N，N-三甲基-色氨酸和芹菜醇。利用Swiss Target Prediction数据平台，预测各成分的潜在作用靶点。其中得到13个益母草碱预测靶点，64个益母草叟苷潜在预测靶点，74个Leoheteronin C潜在靶点，103个汉黄芩素潜在靶点，34个葫芦巴碱潜在靶点，62个N，N，N-三甲基-色氨酸潜在靶点，92个芹菜醇潜在靶点。整合去重后共得到益母草水提物潜在作用靶点357个[1]。

图1 益母草水提物活性成分正离子（A）和负离子（B）模式质谱离子流图Fig.1 Positive(A)and negative(B)total ion chromatograms of LHAE active ingredients

表2 益母草水提物活性成分UPLC-Q-TOF-MS/MS鉴定结果Tab.2 Compounds identified in LHAE active ingredients by UPLC-Q-TOF-MS/MS

3.2 疾病潜在作用靶点筛选及药物成分-疾病作用靶点获取 通过GeneCards、OMIM及DisGeNET数据库分别得到AS相关的疾病基因4 717个、188个、2 044个，整合去重后共得到5 237个疾病基因。随后与2.3得到的357个益母草成分作用靶点取交集，得到259个交集靶点并绘制Venn图，见图2。进一步将交集靶点一一映射至成分潜在作用靶点上，得到药物成分-疾病靶点文件。

图2 药物靶点和疾病靶点VENN图Fig.2 VENN diagram of drug targets and disease targets

3.3 药物防治疾病的调控 利用Cytoscape软件构建“药物-成分-靶点”网络，网络中有7个活性成分，259个靶点，相互作用的边共349条。见图3。

图3 益母草水提物治疗AS的“相关活性成分-潜在作用靶点”网络图Fig.3 Network of“related active components-potential targets”of LHAE in the treatment of AS

3.4 靶点蛋白互作网络构建及核心靶点筛选 将得到的259个作用靶点导入STRING数据库进行PPI分析，利用Cytoscape3.7.1软件进行可视化处理，同时用软件中CytoNCA和BisoGenet插件对PPI网络进行拓扑分析，以 BC、CC、DC、degree值等2倍中位数为筛选标准结果见图4。并进一步筛选出核心靶点同时对核心靶点构建蛋白互作网络图，见图5。

图4 药物-疾病靶点筛选策略Fig.4 Screening strategies of drug-disease targets

图5 核心治疗靶点蛋白互作关系网络图Fig.5 Protein-protein interaction(PPI)network of hit related targets

3.5 “信号通路-靶点”网络图 通过DAVID在线分析平台对作用靶点进行GO功能富集分析（P<0.05），得到 258个 GOTERM-BP-DIRECT条目，259个GOTERM-CC-DIRECT条目，259个GOTERM-MFDIRECT条目，分别取其前10条作图，通过微生信得到图6。GO富集分析显示，益母草防治AS可能涉及：对药物的反应，对细胞外信号调节激酶ERK1、ERK2级联正调控、炎症反应等。

通过DAVID在线分析平台对作用靶点进行KEGG通路富集分析（P<0.05），共得到240条信号通路，取相关性较高的20条作图，利用微生信得到图6。KEGG富集分析显示，益母草可能通过影响磷酯酰肌醇 3激酶-蛋白激酶 B（PI3K-Akt）、TNF、低氧诱导因子-1（HIF-1）、肾素-血管紧张素系统（Ras）等信号通路发挥防治AS的作用。综合3.3与3.4的结果，通过Cytoscape软件对益母草7个核心成分，度值排名靠前的6个核心靶点，以及与其密切相关的15条核心信号通路进行“成分-靶点-通路”网络图的构建，见图7。

图6 生物学过程（BP）、细胞成分（CC）、分子功能（MF）及KEGG通路富集结果Fig.6 Biological process(BP)，cell component(CC)，molecular function(MF)analysis of targets and KEGG pathway enrichment results

图7 成分-靶点-通路网络图Fig.7 Component-target-pathway interaction network

3.6 分子对接分析 利用Auto Dock软件对核心药效成分和核心靶点进行分子对接模拟计算，见表3。对部分结果用Pymol进行可视化分析，见图8。对接结果显示，益母草提取物活性成分与AS核心靶点的最低结合能均小于0 kJ/mol，表示预测的成分与靶点可以自发进行结合，其中益母草叟苷、芹菜醇和 Leoheteronin C 与 IL-6、TNF、IL-1β、VEGFA、MYC、SRC点的结合能均小于-5 kJ/mol，说明其与靶点对接效果好[17]。

图8 分子对接模式图Fig.8 Molecular docking diagram

表3 益母草水提物与核心靶点分子对接结合自由能Tab.3 Binding energies of LHAE active ingredients to core targets KJ/mol

4 讨论

4.1 益母草治疗AS研究基础 本研究基于UPLC-Q-TOF-MS/MS与网络药理学活性成分筛选，发现益母草干预AS的主要活性成分可能是益母草碱、益母草叟苷、汉黄芩素、Leoheteronin C、芹菜醇、葫芦巴碱、N，N，N-三甲基-色氨酸。研究显示，益母草碱以剂量依赖性的方式显著减轻高胆固醇饮食诱导的新西兰雄兔AS的发展，高剂量益母草碱可升高高密度脂蛋白胆固醇、降低三酰甘油，减少VSMC由收缩表型向合成表型的转换和迁移，并改善巨噬细胞浸润，同时抑制血管细胞黏附分子-1（VCAM-1）、TNF-α、IL-6等炎症因子的水平。此外，益母草碱还可呈剂量依赖性地降低血清和肝脏的脂质过氧化水平[37]。Zhang等[38]从益母草中分离出的益母草叟苷A/B/C/D，在游离脂肪酸诱导的肝癌细胞中均有三酰甘油积聚的抑制作用。Ku等[39]发现汉黄芩素可以有效减少人脐静脉内皮细胞中由高糖诱导的血管炎症过程，减少血管通透性、单核细胞黏附、细胞黏附分子的表达，抑制活性氧的生成和核因子-κB（NF-κB）的激活。这些研究提示，益母草碱、益母草叟苷、汉黄芩素等可能是益母草水提物发挥抗AS药效的主要活性成分。

4.2 靶点分析 通过网络药理学成分筛选、靶点预测和网络分析，可发现占据网络中心位置的主要为IL-6、TNF、IL-1β、VEGFA、c-myc等靶点，可能是益母草干预冠状动脉粥样硬化的关键作用靶点。既往研究表明，在AS的发展过程中，IL-6具有影响炎症细胞的聚集的作用，内皮细胞通过分泌趋化因子和增加细胞间黏附分子-1（ICAM-1）的表达来影响IL-6和可溶性IL-6受体复合物的结合，共同促进白细胞的募集和向内膜迁移。Ha等[40]发现胡枝子提取物可显著抑制TNF-α诱导的小鼠动脉内膜单核细胞的浸润、ICAM的表达，以及炎性细胞因子IL-1β、IL-6、TNF-α 和 VCAM1 的表达，从而有效预防早期AS。Zhang等[41]发现TNF-α可能通过增加低密度脂蛋白跨内皮细胞的转运并促进NF-κB和过氧化物酶体增殖剂激活受体-γ（PPAR-γ）的激活，从而促进低密度脂蛋白在血管壁的滞留，加速AS。斑块内血管生成是增加斑块糜烂、破裂风险的重要因素。血管生成由各种生长因子诱导，例如VEGF家族。Mao等[42]发现，高胆固醇血症增加了炎性细胞因子TNF-α的释放，然后通过增强VEGFA/VEGFR-2和成纤维细胞生长因子-2（FGF-2）/纤维细胞生长因子受体-1（FGFR-1）表达促进斑块内血管生成，共同促进高胆固醇新西兰兔AS病变的不稳定性。血管平滑肌细胞过度迁移与病理性内膜增厚相关，Luo等[43]发现川芎内酯对病理性内膜增厚有保护作用，表现为抑制VSMC在体内和体外的迁移，并与下调c-myc蛋白和基质金属蛋白酶-2（MMP-2）蛋白的表达有关。可见，网络药理学预测的益母草水提物抗AS的关键靶点主要与炎症反应、血管新生及平滑肌细胞的增殖迁移有关，可为后续进一步实验提供参考。

4.3 生物富集分析 通过GO功能富集分析可知，生物学进程包括对药物的反应，对ERK1、ERK2级联正调控、炎症反应等，主要涉及蛋白酪氨酸激酶活性、类固醇结合药物结合等方面的功能。通过KEGG富集分析预测，益母草水提物干预AS的信号通路主要是PI3K/Akt信号通路、HIF-1信号通路、肿瘤坏死因子信号通路、Ras信号通路。

PI3K/Akt信号已被证明参与许多生理和病理过程，包括细胞增殖、细胞凋亡、炎症、缺血性损伤和肿瘤进展[44]。巨噬细胞衍生的泡沫细胞是早期AS病变的主要细胞类型，PI3K/Akt信号通路介导巨噬细胞M1和M2极性改变和凋亡影响AS的各个阶段，包括早期病变、斑块进展和斑块稳定性[45]。

HIF-1信号通路可通过上调NF-κB、血管内皮生长因子及其受体、活性氧、自由基和一氧化氮合酶（iNOS）的表达，导致内皮细胞功能障碍、增殖、迁移和血管生成，从而促进AS的发展。此外，HIF-1通过促进巨噬细胞M1表型分化和骨桥蛋白、TNF-α、iNOS、IL-6和VEGF的表达，促进巨噬细胞的炎症和迁移。HIF-1α还通过减少胆固醇外流机制，促进泡沫细胞的形成[46]。Shan等[47]发现HIF-1α参与了缺氧诱导的KLF4表达上调，从而促进了血管平滑肌细胞的迁移，影响血管重塑。

Ras是多种细胞外信号通路的上游。Ras调节的信号通路控制包括增殖、生长和细胞周期进程在内的生物过程。Yu等[48]发现Ras在VSMCs中的作用包括 ERK、FAK、PI3K/Akt的激活、MMPs的分泌以及其他小G蛋白的调节，阿卡波糖对A7r5细胞中Ras的表达有抑制作用，从而减少A7r5细胞的增殖、迁移。Lin等[49]发现哌嗪二酮衍生化合物TW-01能抑制Ras的激活，抑制MAPK及其下游c-fos、c-jun和c-myc mRNA表达，显著抑制血管平滑肌细胞的增殖以及大鼠颈总动脉球囊损伤后内膜增生。

5 结语

本研究基于UPLC-Q-TOF-MS/MS技术定性表征益母草水提物的化学成分，并结合网络药理学、分子对接筛选并验证其主要活性成分的作用靶点，初步探讨了益母草水提物治疗AS的潜在作用机制。结果表明，益母草水提物可能通过调控IL-6、TNF、IL-1β、VEGFA 和 c-myc等核心靶点及 PI3K/Akt信号通路、TNF信号通路、HIF-1信号通路发挥治疗AS作用。本研究为进一步阐明益母草水提物治疗AS机制提供了一定的研究方向和数据支持，但存在以下不足：1）仅对益母草水提物化学成分进行了定性分析，并未进行定量分析。2）通过网络药理学及分子对接研究得到的仅是预测性结果，益母草提取物及其活性成分防治AS的确切作用与机制仍需进一步实验研究验证。