神经药理学报 ›› 2026, Vol. 16 ›› Issue (3): 1-.DOI: 10.3969/j.issn.2095-1396.2026.03.001

• 研究论文 •    下一篇

基于网络药理学 - 分子对接 - 分子动力学模拟探索赤芍治疗脓毒症相关性脑病的多靶点机制

贾丽芳,李璐璐,任钰凡,张思博,尤斯涵,王静   

  1. 1. 河北北方学院生命科学研究中心,张家口,075000,中国 

    2. 河北北方学院药学院,张家口,075000,中国 

    3. 河北北方学院研究生学院,张家口,075000,中国 

    4. 张家口学院,张家口,075000,中国 

    5. 河北北方学院微循环研究所,张家口,075000,中国

  • 出版日期:2026-06-26 发布日期:2026-08-04
  • 通讯作者: 王静,副研究员,博士;研究方向:中医药治疗脓毒症
  • 作者简介:贾丽芳,研究实习员,硕士;研究方向:中医药治疗脓毒症
  • 基金资助:
    河北北方学院省属高校基本科研业务费自然科学项目(No.JYT2025008)

Exploring the Multi-Target Mechanisms of Paeoniae Radix Rubra in Treating Sepsis-Associated Encephalopathy Based on Network Pharmacology, Molecular Docking, and Molecular Dynamics Simulation

JIA Li-fang, LI Lu-lu, REN Yu-fan, ZHANG Si-bo, YOU Si-han, WANG Jing   

  1. 1. Life Science Research Center, Hebei North University, Zhangjiakou, 075000, China 

    2. College of Pharmacy, Hebei North University, Zhangjiakou, 075000, China 

    3. Graduate School, Hebei North University, Zhangjiakou, 075000, China 

    4. Zhangjiakou University, Zhangjiakou, 075000, China 

    5. Microcirculation Research Institute, Hebei North University, Zhangjiakou, 075000, China

  • Online:2026-06-26 Published:2026-08-04

摘要:

目的:通过网络药理学、分子对接和分子动力学模拟,探究赤芍(paeoniae radix rubra, PRR)治疗脓毒症 相关性脑病(sepsis-associated encephalopathy, SAE)的潜在作用机制。方法:使用 TCMSP 数据库和 HERB 数 据库收集 PRR 活性成分,并利用 PubChem 和 Swiss Target Prediction 获取成分对应的靶点。使用 OMIM、TTD、 GeneCards 和 DrugBank 数据库获取 SAE 相关靶点,通过交集运算,得到 PRR 成分和 SAE 共有的靶点。使用 STRING 数据库构建蛋白质相互作用网络,利用 Cytoscape 软件结合 MCC、Degree、MNC、EPC 四种算法识别核 心靶点。以交集靶点为研究对象,利用 DAVID 在线数据库进行基因本体 (gene ontology,GO) 和京都基因与基 因组百科全书 (kyoto encyclopedia of genes and genomes,KEGG) 富集分析。选取 Degree 值排名前 10 位活性成 分与筛选得到的核心靶点进行分子对接验证,并进一步纳入代表性单萜苷类成分芍药苷进行补充验证。进一步 利用 GROMACS v2020.6 对结合能最低的复合物进行 100 ns 的分子动力学(molecular dynamics,MD)模拟,评 估结合稳定性。结果:共筛选出 PRR 活性成分 76 个,对应靶点 650 个,SAE 相关靶点 752 个,PRR-SAE 共同靶 点 80 个。网络拓扑分析提示 Baicalein、Kaempferol、Tryptanthrin、n-(4-Chlorobenzoyl)-melatonin、Albiflorin 等 可能为 PRR 干预 SAE 的重要活性成分。核心靶点涉及 ALB、TNF、STAT3、HIF1A、MMP9 等。GO 富集分析获 得 542 个条目,KEGG 分析显示 110 条信号通路显著富集。Degree 值前 6 位活性成分与核心靶点的分子对接结 合能均小于 -5.0 kcal·mol-1,表明具有良好的结合亲和力。此外,代表性单萜苷类成分芍药苷虽未进入 Degree 值 前 10,但仍被纳入分子对接验证,并与多个核心靶点表现出良好结合活性。分子动力学模拟进一步证实复合物 结构稳定,结合亲和力强。结论:该研究初步揭示 PRR“多成分 - 多靶点 - 多通路”干预 SAE 的潜在分子机制, 为其药理作用的深入研究提供理论支撑。

关键词: 赤芍, 脓毒症相关性脑病, 网络药理学, 分子对接, 分子动力学模拟

Abstract:

Objective: To explore the potential mechanism of paeoniae radix rubra (PRR) in the treatment of sepsis-associated encephalopathy (SAE) through network pharmacology, molecular docking and molecular dynamics simulation. Methods: The active components of PRR were collected from the TCMSP and HERB databases, and the corresponding targets were obtained from the PubChem and Swiss Target Prediction databases. SAE-related targets were collected from the GeneCards, OMIM, TTD and DrugBank databases, respectively. Common targets between PRR and SAE were identified by intersection analysis. A protein-protein interaction (PPI) network was constructed using the STRING database and core targets were identified using Cytoscape software with four algorithms: MCC, Degree, MNC and EPC. GO and KEGG enrichment analyses were performed using the DAVID database. The top 10 active components ranked by Degree value and the screened core targets were subjected to molecular docking with the core targets, and paeoniflorin, a representative monoterpene glycoside of PRR, was additionally included for supplementary validation. Further, the complex with the lowest binding energy was selected for a 100 ns molecular dynamics (MD) simulation using GROMACS v2020.6 to evaluate the binding stability. Results: A total of 76 active components and 650 corresponding targets of PRR were identified, along with 752 SAE-related targets. 80 overlapping targets between PRR and SAE were identified. Network topology analysis suggested that Baicalein, Kaempferol, Tryptanthrin, n-(4-Chlorobenzoyl)-melatonin, and Albiflorin may represent important active components involved in the therapeutic effects of PRR against SAE. Core targets included ALB, TNF, STAT3, HIF1A, MMP9, etc. GO enrichment analysis obtained 542 terms, and KEGG analysis revealed 110 significantly enriched signaling pathways. The molecular docking results showed that the binding energies between the top 10 active components and core targets were all lower than -5.0 kcal·mol-1, indicating favorable binding affinity. In addition, paeoniflorin, a representative monoterpene glycoside of PRR, was further included in molecular docking validation and exhibited favorable binding activity with multiple core targets. Molecular dynamics simulation further confirmed the structural stability of the complex. Conclusion: This study preliminarily revealed the potential molecular mechanism underlying the therapeutic effects of PRR on SAE through a multi-component, multi-target and multi-pathway mode, providing a theoretical basis for further pharmacological research.

Key words: paeoniae radix rubra, sepsis-associated encephalopathy, network pharmacology, molecular docking, molecular dynamics simulation

中图分类号: