基于网络药理学探究滇黄精抗疲劳作用及其机制*

张 楠，杨 帆，马仲帅，马舒韵，李曦明，温燕龙，李凌飞

(云南农业大学 食品科学技术学院，云南 昆明 650201)

随着社会的飞速发展，疲劳已成为常见的健康问题，体力疲劳、心理疲劳、脑力疲劳和病理疲劳[1]会严重影响人们的工作效率并危害人体健康，主要表现为注意力不集中、嗜睡、情绪低落、从事正常活动能力减弱、体力不支等症状[2]，更有甚者会增加多发性硬化症、脑出血、神经衰弱等疾病的患病率[3]，长期疲劳会降低人体免疫力，促使病原微生物入侵机体，引发疾病[4]，过度疲劳甚至会危及生命。滇黄精 (Polygonatum kingianum) 是百合科(Liliaceae)黄精属(Polygonatum)多年生草本植物，主要分布于云南、贵州、四川等地；其主要成分滇黄精多糖具有增强免疫能力、抗肿瘤、减缓衰老、抗炎等生物活性[5-6]，滇黄精皂苷具有降血糖、促凝血等作用[7]，滇黄精黄酮具有抗氧化、抗疲劳等作用[8-9]。尽管已有报道表明滇黄精具有抗疲劳功效，但目前对滇黄精的抗疲劳作用机制尚不清楚。

网络药理学于2007 年由药理学专家ANDREW 提出，随后应用于众多研究领域。网络药理学旨在针对药物中的某种或多种成分对某种或多种疾病进行分析，它以拓扑网络的形式，强调药物对信号通路进行多途径调控，旨在改善药物的有效性，减少其不良反应，降低新药的研发成本和缩短药物的研发时间[10]。近年来，网络药理学广泛应用于药物活性成分的作用机制预测研究。李非洲等[11]通过网络药理学分析方法探讨人参养荣汤治疗慢性疲劳的作用机制，发现其主要通过作用于MAPK14、NOS2、PTGS2 等多个靶点调控PI3K/Akt 信号通路发挥抗疲劳作用；高寒等[12]也利用网络药理学方法探究了刺五加总苷抗疲劳的潜在作用机制，发现刺五加总苷中的7 种活性成分通过作用于HSP90AA1、STAT3、VEGFA 等靶点而发挥抗疲劳作用。本研究基于药物—靶点—途径的研究思路，通过多平台数据库获取滇黄精活性成分及其抗疲劳作用靶点，使用Cytoscape 和Venny 2.1.0 软件构建滇黄精—活性成分—抗疲劳拓扑网络，采用网络药理学分析技术及分子对接技术探索滇黄精抗疲劳的作用靶点及其机制，以期为滇黄精的抗疲劳功效提供理论依据。

1 材料与方法

1.1 检索滇黄精的活性成分

通过中药系统药理学数据库与分析平台 (Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform，TCMSP) 和CNKI 数据库查询滇黄精的活性成分。首先，使用TCMSP 数据库，设置“Herb Name”为“滇黄精”进行搜索，导出其活性成分；其次，使用CNKI 数据库搜寻与滇黄精活性成分相关的文献并导出。去除上述步骤所获得的重复活性成分，即获得滇黄精的活性成分。在活性成分筛选时，通常需要满足口服利用度大于30%和生物类药性大于0.18 的条件[13]。在本研究中，虽然滇黄精中的一些化学成分(如滇黄精皂苷、薯蓣皂苷、重楼皂苷等)不满足口服利用度和生物类药性的条件，但相关文献[14]报道了这些活性成分具有显著的药理作用，因此，也将这些成分保留作为滇黄精的活性成分。

1.2 滇黄精活性成分的靶点预测

SwissTargetPrediction 数据库是基于已知化合物的二维和三维结构的相似度，对化合物治疗靶点展开预测分析的工具。本研究通过小分子有机生物活性数据库 (PubChem) 查询滇黄精活性成分的分子线性输入规范 (SMILES) 结构式，将其输入到SwissTargetPrediction 数据库中，并将种属设置为“homo sapiens”，选取Probability＞0 的靶点，得到滇黄精活性成分潜在靶点的预测结果。

1.3 滇黄精抗疲劳靶点的筛选

在Disgenet (https://www.disgenet.org/search)数据库的搜索框中输入“Fatigue” (Name: Fatigue;CUI: C0015672::)，选择“Summary of Gene-Disease Associations”，获取与疲劳相关的靶点。由于1.2节中通过SwissTargetPrediction 预测得到的滇黄精活性成分作用靶点是所有疾病的靶点，因此，需借助Venny 2.1.0 软件进一步筛选出疲劳有关靶点与滇黄精活性成分作用靶点的交集，从而获得滇黄精抗疲劳的靶点。

1.4 建立滇黄精—活性成分—抗疲劳靶点网络

将1.1 节和1.3 节获得的数据导入Cytoscape v3.9.0 软件中，构建滇黄精—活性成分—抗疲劳靶点的网络图，其中节点度值越大，则其节点颜色越深。

1.5 蛋白质互作网络构建

将1.3 节筛选出的靶点上传到String 数据库，种群选择为“homo sapiens”，将相互作用评分调整为最高置信度 (highest confidence)=0.9，隐藏没有相互作用的蛋白靶点，获得蛋白互作 (proteinprotein interaction，PPI) 网络图，并将所获得的数据导出至Cytoscape 软件。通过节点度值设置节点大小和颜色，得到PPI 网络拓扑图。

1.6 GO 富集分析与KEGG 通路分析

为更好地探究滇黄精抗疲劳靶点涉及的生物学过程及通路，将滇黄精抗疲劳靶点导入至DAVID 数据库，在第1、2、3 列中分别选择“Official Gene Symbol”“Gene List”和“Homo Sapiens”，根据P＜0.05 从获得的数据中进一步筛选GO 富集结果和KEGG 通路。

1.7 分子对接分析

分子对接技术能直观地显示滇黄精活性成分与作用靶点的结合程度。从PubChem 中获取滇黄精主要活性成分的3D 结构，导入至Open Babel 2.3.2 软件，将文件格式转化为PDB 格式；将关键靶点导入到PDB 数据库中，得到其3D 构造的PDB 格式文件。将所获取的活性成分和靶蛋白导入至Pymol 软件，活性成分进行加氢处理，靶蛋白进行去水分子处理。使用AutoDock Vina 对主要活性成分与关键靶点进行分子对接，结合能与氢键个数反映了活性成分与关键靶蛋白结合的稳定性，结合能越低，对接的稳定性能越好；对接氢键数量越多，对接之后的化合物构象越稳定。

2 结果与分析

2.1 滇黄精的活性成分

通过TCMSP 和Drugbank 查询滇黄精的活性成分，将查询到的活性成分进行收集、整理及去除重复数据后，共获得38 种活性成分(表1)。

表1 滇黄精活性成分Tab.1 Active ingredients of Polygonatum kingianum

2.2 滇黄精潜在抗疲劳靶点预测

通过Disgenet 数据库查询到730 个疲劳相关靶点，通过SwissTargetPrediction 数据库预测得到滇黄精的所有作用靶点共129 个，使用Venny 2.1.0筛选到滇黄精的抗疲劳靶点30 个(图1)。

图1 滇黄精潜在的抗疲劳靶点数Fig.1 The number of potential antifatigue target of Polygonatum kingianum

2.3 滇黄精—活性成分—抗疲劳靶点网络

由滇黄精活性成分—抗疲劳靶点的网络图(图2)可知：包含35 个活性成分、66 个节点、187条边。除滇黄精外，3′-甲氧基大豆苷元、黄芩素和β-谷甾醇3 种成分的节点颜色最深，表明它们是滇黄精抗疲劳的主要活性成分。此外，靶点的每1 个节点均有多条边与之连接，也表明了滇黄精的抗疲劳功效是不同化合物协同作用的结果。

图2 滇黄精—活性成分—抗疲劳靶点网络Fig.2 The target network of P.kingianum-active ingredients-antifatigue targets

2.4 蛋白质互作分析

拓扑图(图3)包含22 个节点，47 条边。滇黄精抗疲劳的关键靶点为表皮生长因子受体 (EGFR)、信号转导及转录激活蛋白3 (STAT3)、蛋白酪氨酸磷酸酶非受体 (PTPN11 和JAK2)、酪氨酸激酶受体2 (ERBB2) 以及雌激素受体 (ESR1)。

图3 蛋白互作网络Fig.3 Protein-protein interaction network

2.5 滇黄精抗疲劳基因GO 分析

滇黄精抗疲劳关键靶点共涉及68 个条目。前30 个条目中(图4)，生物学过程包括RNA 聚合酶Ⅱ启动子转录正向调节、转录正向调节等13个生物学过程；细胞组成包括细胞溶质、细胞核、核质等7 个组分；分子功能包含ATP 结合、蛋白激酶结合、酶结合等10 种功能。

图4 滇黄精抗疲劳GO 富集分析Fig.4 Antiatigue GO enrichment analysis of P.kingianum

2.6 滇黄精抗疲劳KEGG 通路分析

对滇黄精抗疲劳的关键靶点进行KEGG 通路分析，得到的结果数据以P＜0.05 进行进一步筛选，共计得到39 个通路，其中基因富集数大于5 的通路主要包含癌症中的蛋白聚糖、耐药性EGFR 酪氨酸激酶抑制剂、PI3K-Akt 信号通路等10个通路(图5)。

图5 滇黄精抗疲劳KEGG 通路分析Fig.5 Antifatigue KEGG pathway analysis of P.kingianum

2.7 分子对接

将蛋白质互作网路筛选出的6 个关键靶蛋白(EGFR、STAT、PTPN11、ERBB2、ESR1 和JAK2)分别与滇黄精的3 个抗疲劳活性成分(黄芩苷元、3′-甲氧基大豆苷元和β-谷甾醇)进行分子对接，结果(表2 和图6)显示：滇黄精主要活性成分与关键靶蛋白对接的结合能均小于-5 kcal/mol，表明预测的活性成分与关键靶点结合性能良好。其中，黄芩苷元对关键靶蛋白的结合能低且结合氢键数量最多，表明黄芩苷元在滇黄精抗疲劳过程中发挥着重要作用。

图6 滇黄精活性成分(黄芩苷元)与关键靶点分子对接的可视化图Fig.6 Visualization of molecular docking between active ingredient (baicalein) and core target of P.kingianum

表2 分子对接结果Tab.2 Results of molecular docking

3 讨论

本研究通过网络药理学方法对滇黄精潜在的抗疲劳作用机制进行分析，结果表明滇黄精的抗疲劳作用具备多成分、多靶点的特性，且由不同的活性成分发挥抗疲劳功效。拓扑网络分析显示：3′-甲氧基大豆苷元、黄芩苷元和β-谷甾醇的节点最大，表明这3 种成分是滇黄精抗疲劳的主要活性成分。PPI 网络表明：滇黄精活性成分与靶蛋白之间具有较强的互作效应，EGFR、STAT3、PTPN11、ERBB2、ESR1 和JAK2 为滇黄精抗疲劳的关键靶点；分子对接结果表明：滇黄精3 种抗疲劳主要活性成分均能与抗疲劳的6 个关键靶蛋白自发结合，且黄芩苷元的抗疲劳活性最为重要。黄芩苷元属于黄酮类化合物，主要来源于黄芩根，研究表明其具有良好的抗炎、抗氧化、抗过敏以及抗癌的生物功效[15]，并能显著抑制炎症因子TNF-α 和IL-6 的表达[16]，从而缓解由炎症因子引起的疲劳。陈景等[17]研究发现：黄芩苷元能上调大鼠海马CA1 区超氧化物歧化酶和谷胱甘肽过氧化酶的表达，降低丙二醛含量，上调磷酸化核因子E2 相关因子2 及血红素氧合酶1 的蛋白表达。3′-甲氧基大豆苷元为异黄酮类物质，具有镇痛[18]、清除自由基[19]、抗炎[20]等作用。β-谷甾醇属于植物甾醇类化合物之一，具有抗炎[21]、抗菌[22]、降脂[23]以及提高人体免疫力的作用[24]。

GO 分析表明：滇黄精的抗疲劳作用涵盖细胞炎症、细胞生长、能量代谢以及自由基代谢等多个生物学过程。对细胞组分的分析发现：滇黄精抗疲劳的主要作用靶点涉及胞浆、核质、线粒体、染色质等。线粒体是骨骼肌中重要的细胞器，能为骨骼肌的舒张和收缩提供能量。WANG 等[25]研究发现：Nrf2 介导的氧化应激与疲劳有关，且血红素氧合酶1 能够促进线粒体合成[26]，加速机体产生能量，缓解疲劳。分子功能包含了ATP 结合、蛋白激酶结合、锌离子结合等作用。其中，锌是人体所必需的微量元素，也作为许多酶的金属催化剂，可通过合成RNA 聚合酶和DNA 聚合酶加速DNA、蛋白质以及胶原的合成，同时促进维生素A 的吸收、组织修复及再生，与肌肉质量关系密切[27]。此外，锌是重要的抗氧化剂，能有效减少运动过程中机体代谢产生的自由基，尤其是在高强度的激烈运动后，适量补充锌有利于运动后自由基的消除[28]。蛋白激酶存在于细胞质中，能够调节人体的神经系统，对神经系统的发育、行为的形成、长期记忆等具有正向调节作用。田振军等[29]研究证明蛋白激酶类基因可缓解运动性中枢疲劳的产生。

KEGG 通路分析表明：滇黄精抗疲劳作用主要涉及MAPK、PI3K-AKT、JAK-STAT 等信号通路。其中，MAPK 信号通路主要参与细胞炎症、细胞凋亡、免疫调节、HSCs 活化增殖和生成胶原的重要调节途径[30]。李丹等[31]通过小鼠疲劳跑台实验和悬尾实验证明艾灸能够介导JAKSTAT 的信号通路缓解疲劳。郭静等[32]利用大黄素衍生物作用于RAW264.7 炎症细胞，通过抑制PI3K-Akt 途径能够减少细胞炎症的产生。研究表明：产生疲劳的两大主要因素是细胞炎症因子水平的上升以及机体内抗氧化酶的消耗和代谢物的积累[33]，因此，本研究中滇黄精抗疲劳的主要机制可能是滇黄精通过调控MAPK 通路、PI3KAkt 途径和JAK-STAT 信号通路，降低细胞炎症水平、加快细胞周期进程及减少代谢产物的积累，从而发挥抗疲劳作用。这些结果表明滇黄精的抗疲劳作用由多个生物学过程共同参与，同时也表明了滇黄精的抗疲劳作用是多成分、多靶点、多途径共同参与。

4 结论

本研究通过网络药理学方法和分子对接探究了滇黄精的抗疲劳机制，表明滇黄精抗疲劳功效具有多成分、多靶点、多通路共同参与的特点，可通过抗炎、加速能量及自由基代谢、加速细胞周期进程等细胞生物学过程共同发挥抗疲劳作用。研究结果为滇黄精抗疲劳功效提供了理论依据。