基于网络药理学及分子对接技术探讨头花蓼治疗糖尿病的预防作用

Mechanism of Cephalophyllum cephalus Intervention in Diabetes Based on Network Pharmacology and Molecular Docking Technology

摘要: 目的：借助网络药理学及分子对接方法探讨头花蓼治疗糖尿病的作用机制，为扩大头花蓼用药范围提供依据。方法：通过文献研究及Swiss数据库筛选得到头花蓼活性成分并分析相关作用蛋白靶点，借助OMIM、GeneCards、DRUG BANK、DisGeNet等数据库分析糖尿病基因靶点，并取交集；借助STRING数据库及Cytoscape软件构建PPI网络关系并得到头花蓼干预糖尿病关键靶点；借助metascape数据库进行生物富集分析；使用autodock软件进行核心化合物–蛋白分子对接。结果：得到头花蓼核心化合物为49个，涉及45个蛋白靶点，其中与糖尿病相关核心靶点为ALB、EGFR、SRC、CASP3等；在60组分子对接结果中，Affinity < −3 kcal∙mol⁻¹的对接结果有45个，提示多数化合物–蛋白有较大结合能；KEGG富集分析结果显示关键通路为Lipid and atherosclerosis、Pathways in cancer、IL-17 signaling pathway等。结论：头花蓼干预糖尿病依靠多通路、多靶点协调作用，其主要活性成分为槲皮素、红景天甘、木犀草素、齐墩果酸等，通过ALB、EGFR、SRC、CASP3等蛋白作用与Lipid and atherosclerosis、Pathways in cance、IL-17 signaling pathway等通路发挥治疗糖尿病作用。

Abstract: Objective: To explore the mechanism of action of Polygonum cuspidatum in the treatment of diabetes mellitus with the help of network pharmacology and molecular docking methods, and to provide a basis for expanding the scope of drug use of Polygonum cuspidatum. Methods: We obtained the active ingredients of Polygonum cuspidatum through literature research and screening of Swiss database and analyzed the related protein targets, analyzed the diabetes gene targets with the help of OMIM, GeneCards, DRUG BANK, DisGeNet and other databases, and took the intersection; constructed PPI network relationships with the help of STRING database and Cytoscape software and obtained the key targets of Polygonum cuspidatum for diabetes intervention were obtained by using STRING database and Cytoscape software; bioenrichment analysis was performed by using metascape database; core compound-protein molecular docking was performed by using autodock software. Results: 49 core compounds were obtained, involving 45 protein targets, including ALB, EGFR, SRC and CASP3, which are related to diabetes. Among the 60 groups of molecular docking results, 45 docking results with Affinity < −3 kcal∙mol−1 suggested that most of the compounds had large binding energy to proteins; the results of KEGG enrichment analysis showed that the key pathways were Lipid and atherosclerosis, Pathways in cancer, IL-17 signaling pathway, etc. Conclusion: Polygonum cuspidatum intervenes in diabetes mellitus relying on multi-pathway and multi-target coordination, and its main active ingredients are quercetin, rhodopsin, lignan, oleanolic acid, etc., which act with Lipid and atherosclerosis, Pathways in cance, IL-17 signaling pathway through ALB, EGFR, SRC, CASP3 and other proteins signaling pathway and other pathways to play a role in the treatment of diabetes.

1. 引言

糖尿病是一种以多饮、多食、多尿、消瘦为主要临床表现、血糖代谢紊乱为病理特征的代谢障碍性疾病 [1] 。主要分为1型和2型糖尿病，其中以2型糖尿病最为普遍。糖尿病的发病率逐年升高，据国际糖尿病联合会统计，中国约有1.164亿人患有糖尿病，居全球糖尿病患者之首 [2] 。临床降糖药物治疗疗效快、作用机制明确，但存在副作用大、费用高昂的缺点 [3] 。中草药头花蓼(Polygonum capitatum Buch.-Ham. ex D. Don)作为我国民间的常用中草药，其资源丰富、价格低廉、毒副作用小、不易产生耐药性，相关基础及临床研究已证实其具有抗炎、抗菌、抗氧化损伤等作用，因此推测其可为优质降糖药物的研制开发提供广阔空间 [4] [5] 。由于目前尚缺乏头花蓼对糖尿病治疗的作用机制，物质基础的具体研究，本文利用网络药理学研究方法，筛选头花蓼干预糖尿病的化学活性物质及作用靶点、机制，为后期该药物临床应用提供科学理论依据。

2. 材料与方法

2.1. 头花蓼化学成分及靶点筛选

热淋清颗粒是由传统苗药头花蓼为原料的单方制剂，本文通过对头花蓼相关化学成分实验归纳总结，借助chemdraw软件绘制各化合物结构并导入SwissADME数据库，以胃肠吸率“高”，五项类药性评价中至少获得2个“是”为原则，若有研究已证实与糖尿病相关性较高的化合物，则仍保留入选头花蓼干预糖尿病潜在化合物。将所得化合物逐一导入PharmMapper数据库得到各化合物核心作用靶点，以NF ≥ 0.9为标准筛选关键靶点，并结合相关文献研究，最终获得头花蓼核心基因靶点。

2.2. 糖尿病相关靶点筛选

以“diabetes mellitus”、“glycuresis”、“diabetes”等为检索关键词，检索GeneCards数据库、DRUGBANK数据库、OMIM数据库、TTD数据库获取糖尿病核心基因靶点信息。在这些数据库中，DRUGBANK数据库以临床实验结果构建，GeneCards、OMIM、Therapeutic Target Database数据库基于文献研究构建，因此选用多个数据库进行糖尿病基因靶点获取有益于得到更加系统的靶点数据。

2.3. 头花蓼成分–糖尿病靶点

将“2.2.”项所得靶点信息导入STRING分析平台，以物种为“智人”，置信度为0.4为标准，分析建立PPI网络拓扑关系，借助CytoScape 3.7.2软件进一步优化结果，获得头花蓼干预糖尿病关键基因靶点信息。

2.4. GO功能富集及KEGG通路富集

将“2.3.”项结果导入Metascape平台进行核心基因靶点的生物信息数据富集分析，其中包括GOBP (生物过程) GOMF (分子功能) GOCC (细胞成分)及KEGG通路。将筛选结果构建富集条形分析图及通路气泡图。

2.5. 成分–靶点–通路网络图的构建

整理“2.4.”项分析结果及头花蓼–糖尿病交集基因数据并导入Cytoscape3.7.2软件建立：头花蓼成分–靶点–通路网络拓扑图，根据介度、连接度等分析得到该网络中的核心靶点及药物干预疾病的核心成分。

2.6. 分子对接计算验证

通过对2.3.项结果分析得到该网络中度值排名最高的6个靶点蛋白，并结合RCSB平台下载该蛋白pdb格式，借助Openbabel软件将头花蓼活性前10化合物转化为pdb格式。借助pymol及Autodock软件对蛋白及化合物进行预处理并进行分子对接，以结合能(Affinity)作为判断二者结合活性的参考。

3. 结果

3.1. 头花蓼主要化学成分基因靶点预测

通过文献研究共收集到100个头花蓼所含化合物，其中主要以黄酮、木质素及挥发油类为主。利用SwissADME平台共筛选出49个化合物，利用PharmMapper平台进行预测并去除所有重复基因后共计得到114个预测靶点。

3.2. 头花蓼与糖尿病相关基因靶点获取

利用Genecards平台共收集到糖尿病相关靶点17,703个，选取Score ≥ 3.88的结果结合TTD、DRUGBANK、OMIM数据库进行补充，共获得1386个核心疾病靶点。将头花蓼核心靶点与糖尿病核心靶点取交集，VENN分析得到45个核心作用靶点，见图1。

Figure 1. Drug-disease target VENN diagram

图1. 药物–疾病靶点VENN图

Figure 2. Polygonum cephalicum-diabetes target PPI network

图2. 头花蓼–糖尿病靶点PPI网络

3.3. PPI网络构建

将“3.2.”项所得45个交集靶点导入String11.5平台做PPI网络分析，见图2。该网络共44个节点数，225条边，根据度值推测ALB、EGFR、SRC、CASP3等为核心靶点。

3.4. GO功能富集和KEGG通路富集

利用Metascape平台对交集基因进行生物信息富集分析，GO分析结果包括GO BP 441条，GO CC 35条，GO MF 63条；KEGG富集分析得到103条信号通路；取GO及KEGG分析结果排名前20建立生物信息富集分析图，见图3。

(A) (B) (C) (D)

Figure 3. Enrichment analysis of potential targets of major components of Polygonum cephalicum (Note: A: GO-BP analysis; B: GO-CC analysis; C: GO-MF analysis; D: KEGG analysis)

图3. 头花蓼主要成分潜在靶点的富集分析(注：A：GO-BP分析；B：GO-CC分析；C：GO-MF分析；D：KEGG分析)

3.5. 头花蓼成分–靶点–通路网络的构建

使用CytoScape3.7.2软件建立“头花蓼成分–靶点–通路”网络拓扑关系图，并分析得到头花蓼干预糖尿病靶点网络信息数据、核心化学成分及作用靶点，见图4。此网络包含117个节点553条边，其中红色节点表示头花蓼潜在作用靶点，橙色节点表示潜在信号通路，绿色节点表示活性化学成分，连线代表不同节点间的作用关系，图中节点面积越大、颜色越深表明其对糖尿病的影响越大。

通过对网络分析发现，该网络平均度值为6.698，平均介度为0.014，平均紧密的为0.317，且头花蓼每种活性成分对应多个靶点，每个靶点连接多种成分，体现了头花蓼内含的多种成分可能通过多靶点干预糖尿病。其中THL27 (槲皮素)度值为39，介度为0.151，紧密度为0.423，推测THL27可能是头花蓼干预糖尿病的核心化学成分(表1)。MAPK10、CA2、MAPK14、MAPK1、MAPK8靶点基因度值排名靠前(表2)。故预测MAPK通路为头花蓼干预糖尿病的核心通路。

3.6. 头花蓼分子对接筛选

将头花蓼中排名前十名的化合物与ALB、SRC、EGFR、ESR1、CASP3、MAPK1 6个核心蛋白进行分子对接，从而得到60组化合物–蛋白对接结果。其中Affinity < −3 kcal∙mol⁻¹的对接结果有45个，Affinity < −10 kcal∙mol⁻¹的有5个，其中对接得分最高的为SRC-THL37，分值为−16.03 kcal∙mol⁻¹。本文对接结果为进一步进行头花蓼相关实验提供理论支撑，对接结果见图5，核心对接模式见图6。

Table 1. Characterization parameters of the network nodes of the main active ingredients of Polygonum cephalicum

表1. 头花蓼主要活性成分网络节点特征参数

Table 2. Characterization parameters of the nodes of the target network of the main active ingredients of Polygonum cephalicum

表2. 头花蓼主要活性成分靶点网络节点特征参数

Figure 4. Component-target-pathway diagram of Polygonum cephalicum intervention for diabeticization

图4. 头花蓼干预糖尿病化成分–靶点–通路图

Figure 5. Molecular docking results

图5. 分子对接结果

Figure 6. Molecular docking pattern of some core compounds of Polygonum cephalicum

图6. 头花蓼部分核心化合物分子对接模式

4. 讨论

2型糖尿病是由机体胰岛素抵抗造成的体内葡萄糖代谢障碍 [6] 。长期持续性的血糖升高会引发心、脑血管、糖尿病肾病、糖尿病足等病变 [7] 并增加多种癌症的风险，包括肝癌、胰腺癌、子宫内膜癌、结直肠癌和绝经后乳腺癌 [8] ，严重残害人们身体健康、影响人们的生活质量。

中医中糖尿病属于“消渴”的范畴，首载于《素问·奇病论篇》。唐代《外台秘要》云“渴而饮水不能多，小便数，阴痿弱，但腿肿 [9] ”指出渴不多饮、小便数、水肿是消渴病的重要表现；《皇帝内经》云“此人必数食甘美而多肥也，肥者令人内热，甘者令人中满，故其气上溢，转为消渴 [10] ”、《灵枢·五变》云“五脏皆柔弱者，善病消……血气逆留，腹皮充胀，血脉不行，转而为热，热则消肌，故为消瘫 [11] ”均认为消渴发病与“热”密切相关。《金匮》认为以“大便难、脉浮数、消谷善饥、小便频数”为主症的消渴属胃热炽盛，应荡涤胃肠，泻热存津 [12] ；《素问吴注》云：“膏粱之人，内多滞热，故其变病，能生大疔 [13] ”，《中藏经》云：“恣饮醇酒、多嗜甘肥、毒鱼醋酱，畜其毒邪，浸渍脏腑，久不摅散，始变为丁 [14] ”认为过食肥甘、膏粱醇酒可致脾胃气机失调、湿热毒邪内蕴以致毒聚肌肤、血肉腐败引发外科疮疡，即对应现代医学中的糖尿病足。消渴病病机复杂，变症繁多，以阴虚为本，燥热为标的阴津亏损、燥热偏盛为其基本病机，根据“多饮”、“多食”、“多尿”并结合“三焦学说”，中医将消渴分为上、中、下三消，以清热解毒、养阴清热、活血补气为其主要治则，

头花蓼为蓼科植物头花蓼的干燥全草或地上部分，是贵州省的特色苗药之一 [15] ，具有清热解毒、利尿通淋的功效，主要用于治疗肾孟肾炎、膀胱炎、尿道感染、疮疡湿疹等疾病 [16] ，具有降糖 [17] 、抗菌 [18] 、抗炎镇痛 [19] 、抗氧化 [20] 等药理活性。以头花蓼为主的单方制剂“热淋清颗粒”已收载于2020版《中华人民共和国药典》一部，具有清热泻火，利尿通淋的功效 [21] 。

本研究利用网络药理学方法，构建成分–靶点–信号通路网络，通过分析相关网络要素(节点、度值等)，探讨了头花蓼治疗糖尿病主要物质基础。

4.1. 潜在有效成分分析

通过网络药理学方法，发现头花蓼中槲皮素、红景天甘，木犀草素，齐墩果酸为潜在治疗糖尿病的有效成分。铁死亡可导致胰腺β细胞(PBC)丢失和功能障碍，研究结果表明，槲皮素除有抗炎、抗氧化作用外 [22] ，还可能通过抑制PBC铁死亡和胰腺铁沉积对糖尿病发挥作用 [23] 。红景天苷是一种对羟基苯乙基-β-葡萄糖苷类化合物，是红景天的主要活性成分，有降血糖、降血脂作用以及神经保护、缓解糖尿病并发症等作用，研究发现，红景天甘可通过激活AMPK及AMPK相关信号通路，调节糖脂代谢、氧化应激反应和炎症反应相关的代谢应激，从而缓解糖尿病 [24] 。木犀草素可通过调节Nrf2介导的氧化应激和NF-κB介导的炎症反应保护糖尿病小鼠的心脏组织 [25] 。齐墩果酸具有广泛的生物活性，研究表明其可通过增强胰岛素反应、保持β细胞的功能和存活以及预防糖尿病并发症来有效对抗糖尿病、代谢综合征与血脂异常 [26] 。

4.2. 靶点分析

网络拓扑分析结果显示，头花蓼活性成分49个，对应114个蛋白靶点，其中与糖尿病相关的共同靶点45个。通过构建PPI网络发现，ALB、EGFR、SRC、CASP3等的度值较大，可能是头花蓼最关键的治疗糖尿病靶点。Caspase-3是一种凋亡蛋白酶，研究发现，糖尿病患者血浆中Caspase-3增加 [27] ，且Caspase-3依赖性细胞凋亡介导的胰岛β细胞功能障碍是糖尿病的潜在病理机制 [28] 。丝裂原活化蛋白激酶(MAPK)信号通路在糖尿病、胰岛素抵抗和炎症的过程中发挥关键作用 [29] 。此外，糖化反应导致的晚期糖基化终末产物(AGEs)，是糖尿病和糖尿病并发症进展的关键因素。AGEs与其主要细胞受体RAGE的结合激活MAPK信号通路，导致氧化应激和炎症增强 [30] 。白蛋白(ALB)是糖尿病的诊断标志物之一，研究证明其与糖尿病肾脏损伤程度呈正相关 [31] 。

4.3. 生物富集分析

KEGG信号通路分析发现头花蓼主要是通过Lipid and atherosclerosis (脂质和动脉粥样硬化)，Pathways in cancer，IL-17 signaling pathway (IL-17信号通路)等通路抵抗糖尿病。高胰岛素血症及血脂异常是糖尿病患癌的重要因素 [32] 。胰岛素与胰岛素受体(IR)结合激活PI3K/Akt通路，IR/PI3K/Akt信号通路激活mTOR及其下游信号通路，参与调节癌细胞的存活、增殖、侵袭、迁移、分化、血管生成和转移 [33] 。此外，高脂血症导致非酯化脂肪酸(NEFA)水平升高并激活Akt/mTOR和Akt/3β (GSK3β)/β-连环蛋白致癌信号通路 [34] 。研究发现，脂质和动脉粥样硬化通路可通过高血糖、氧化应激和炎症增加促进糖尿病的发生发展 [35] 。此外，糖尿病可改变口腔微生物组的组成并增加牙周炎症和骨质流失 [36] ，而IL-17是一种与免疫保护和免疫病理学相关的细胞因子，在慢性牙周炎中发现IL-17升高并刺激促炎介质(如IL-6和RANKL)的产生，且促进破骨细胞生成 [37] 。

5. 总结

本研究借助网络药理学及分子对接技术，利用数据库进行靶点预测、蛋白互作关系、生物富集分析，初步探明头花蓼作用与糖尿病的多通路、多靶点协调作用。头花蓼可能通过槲皮素、红景天甘，木犀草素，齐墩果酸等主要活性成分，作用于ALB、EGFR、SRC、CASP3等蛋白靶点，参与Akt/mTOR、PI3K/Akt、IL-17信号通路等经典通路发挥治疗糖尿病作用。本研究采用网络药理学及分子对接方法初步探讨了头花蓼干预糖尿病的作用机制，为相关新药开发提供理论依据，但本研究仍存在一定的局限性，后续研究团队将在此基础上开展动物体内实验，进一步探明头花蓼干预糖尿病的作用机制，为中医药防治糖尿病提供新思路。

基金项目

贵阳市科技计划项目筑科合同[2020]-10-12号。

NOTES

^*通讯作者。