摘要
利用TCMSP和PharmMapper完成对“杜仲-牛膝”药对有效成分和有效成分靶点的筛选预测;通过相关疾病数据库获取LDH相关的作用靶点;利用Venny分析获取此药对治疗LDH的潜在靶点,将获得的潜在治疗靶点基因上传至STRING数据库进行分析,并通过可视化软件分析蛋白质-蛋白质相互作用网络(PPI),获得药物治疗疾病的关键靶点;利用DAVID数据库对关键靶点蛋白进行生物富集分析。
获得“杜仲-牛膝”的有效成分共45个,其中共同有效成分3个,有效成分作用靶点共287个;LDH疾病靶点基因506个,药物作用于疾病的可能潜在靶点62个;PPI网络分析表明IL6、TNF、IL1B、ALB、FN1、MMP9、AKT1、JUN、PTGS2、TP53等可能是“杜仲-牛膝”治疗LDH的核心靶点;主要通过RNA聚合酶II启动子转录正调控、调节炎症及免疫反应、调节细胞增殖和细胞凋亡、调节血管生成及对缺氧反应等生物学进程,参与TNF信号通路、Toll样受体信号通路、破骨细胞分化、PI3K-Akt信号通路、MAPK信号通路。
腰椎间盘突出症(Lumbar disc herniation,LDH)是骨科临床的一种常见疾病,是由椎间盘退变、突出、压迫或刺激临近组织或相应平面的神经根、马尾神经后表现出来的一种综合
网络药理学是近年来新兴的基于系统生物学,利用计算机进行生物网络分析,构建“药材-成分-靶标”网络,分析通路机制的一门学科,其为药物治疗疾病从靶点分子层面提供了导向。本文利用网络药理学的相关原理和方法,探讨“杜仲-牛膝”治疗LDH的分子机制,为进一步实验研究及临床应用提供思
通过检索中药系统药理学数据库TCMS
通过TCMSP及PharmMapper数据库(http://lilab-ecust.cn/pharmmapper/
通过OMIM(https://omim.org)、DrugBan
通过Venny分析获得药物与疾病共同作用的靶点,绘制韦恩图,将获得的中药有效成分、药物与疾病作用靶点导入Cytoscape 3.7.1软件,构建“有效成分-靶点-疾病”的可视化网络图。
将初步得到的潜在作用靶点输入STRING数据
利用DAVI
通过TCMSP数据库挖掘杜仲及牛膝的药物有效成分,获得“杜仲-牛膝”潜在有效成分。根据ADME参数筛选标准,设置OB≥30%、DL≥0.18为筛选条件,结果:牛膝有效成分共176个,符合条件的有20个;杜仲有效成分共147个,符合条件的有28个;其中两药共同有效成分有3个。见
通过TCMSP及PharmMapper数据库预测所筛选出的“杜仲-牛膝”有效成分的作用靶点,利用Uniprot数据库进行靶点基因名称(Gene name)的标准转换,得到杜仲有效成分作用靶点共260个,牛膝作用靶点共238个,其中两种药物的共同作用靶点211个。利用网络可视化软件Cytoscape 3.7.1进行“杜仲-牛膝”有效成分及作用靶点的网络构建,利用软件中的network analyzer插件进行网络特征分析,以明确“杜仲-牛膝”中有效成分与作用靶点的相互作用关系,用节点(node)表示药物、有效成分及靶蛋白,用边(edge)表示其中的作用关系,节点(node)共334个,作用关系边(edge)共1164条,见

图1 “杜仲-牛膝”药对-有效成分-作用靶点网络图
通过对OMIM、DrugBank、GeneCards等疾病数据库进行筛选,合并去重后通过Uniprot数据库完成靶点基因名称(Gene name)的标准转换,得到目前已知LDH疾病靶点506个。利用Venny分析绘制Venny图获得药物与疾病共同作用的靶点,得到药物与疾病交集62个,代表“杜仲-牛膝”治疗LDH的潜在靶点数量62个,同时通过Cytoscape 3.7.1进行药物-疾病-靶点的可视化分析,见

图2 药物-疾病-靶点韦恩图

图3 药物-有效成分-疾病-靶点网络图
将所得到的“杜仲-牛膝”治疗LDH潜在靶点通过STRING数据库进行PPI分析(DDX6与其他靶点蛋白均未关联,故将其移除),见

图4 STRING绘制的“杜仲-牛膝”药对治疗LDH的PPI网络图

图5 药物治疗疾病的PPI网络图
将药物治疗疾病的潜在靶点导入DAVID数据库,进行GO富集分析,得到生物学进程(Biological Process,BP)366条,细胞组分(Cellular Component,CC)36条,分子功能(Molecular Function,MF)71条。靶点主要参与RNA聚合酶II启动子转录正调控、参与炎症反应、调节细胞增殖和细胞凋亡、免疫应答、脂多糖应答反应、调节血管生成及对缺氧反应等生物学进程,主要在胞外空间、细胞核线粒体等区域发挥作用,主要参与蛋白质、酶和转录因子的结合、细胞因子及蛋白质二聚化活性等MF。见

图6 GO富集分析气泡图
潜在靶点KEGG富集分析得到100条相关信号通路,药物-疾病共同靶点主要参与:类风湿性关节炎、结直肠癌、癌症中的蛋白多糖、癌症中的途径、HTLV-I感染、甲型流感、乙型肝炎、结核、阿米巴病、弓形虫、恰加斯病(美国锥虫病)、利什曼病、百日咳、沙门氏菌感染、TNF信号通路、Toll样受体信号通路、破骨细胞分化、PI3K Akt信号通路、细胞因子-细胞因子-受体相互作用、MAPK信号通路。见

图7 KEGG富集分析
将筛选得到的“杜仲-牛膝”有效成分及药物与LDH的共同作用靶点和靶点参与数量前20的信号通路导入Cytoscape 3.7.1进行可视化分析,同时利用软件内数据分析插件进行网络特征分析,节点颜色及大小反映重要程度,得到节点共119个,相互作用关系623条。见

图8 成分-疾病-靶点-信号通路网络图
腰椎间盘突出症归属于中医学“腰痛”“痹病”的范畴,西医对症治疗的不良反应较多,而疾病早期及时的中医药干预能有效改善症状,减轻患者心理及身体负担,改善生活质量。其中杜仲-牛膝作为治疗LDH的常用药对,临床使用频率高且疗效显著。本文基于网络药理学的相关理论及技术,研究了“杜仲-牛膝”联合治疗LDH的分子作用机制。
通过筛选得到“杜仲-牛膝”有效成分共45种,有效成分作用靶点共287个,其中治疗LDH的潜在靶点62个。根据构建的“药物-疾病-靶点-信号通路”网络图,杜仲-牛膝治疗LDH的主要作用成分有槲皮素(quercetin)、山奈酚(aempferol)、汉黄芩素、β-胡萝卜素、去氢双丁香酚等。现代药理学研究表明,槲皮素具有抗肿瘤、抗氧化、抗炎、免疫调节等药理活性,能促进小鼠脾淋巴细胞的增殖,增强小鼠免疫功能,又能抑制其T、B细胞的增殖,抑制TNF-α、IL-6、INF-γ等促炎因子的分
药物-疾病共同靶点PPI网络图结果中,Degree值前10分别为IL6、TNF、IL1B、ALB、FN1、MMP9、AKT1、JUN、PTGS2、TP53。白细胞介素(interleukin,IL)能通过影响基质金属蛋白酶(MMPs)及抑制MMPs中的蛋白多糖合成而作用于椎间
生物富集分析结果显示,“杜仲-牛膝”治疗LDH主要涉及的信号通路包括TNF信号通路、Toll样受体信号通路、破骨细胞分化、PI3K-Akt信号通路、MAPK信号通路等,主要涉及细胞分化和凋亡、炎症反应相关及代谢(如内分泌、氧化应激等)等相关生物学进程。TNF(肿瘤坏死因子)主要引起细胞凋亡,目前主要研究的有TNF-α和TNF-β两大类,TNF相关信号通路主要参与细胞凋亡、调节机体免疫及炎症反应等进程,椎间盘中髓核细胞产生的TNF-α是引起椎间盘退变的主要炎性物质之一,其作用主要表现在可增加椎间盘细胞中基质金属蛋白酶(MMPs)的有效表达,并使基质中的蛋白多糖、弹性蛋白及胶原等大量减
本文通过网络药理学分析,发现“杜仲-牛膝”通过多成分、多靶点、多通路发挥抗炎、抗氧化、调节细胞增殖分化及细胞凋亡的药理作用治疗LDH,从靶点分子层面为药物作用于疾病提供了参考。此次研究既分析了“杜仲-牛膝”的主要成分及作用靶点,为进一步研究药物治疗不同疾病做基础;同时分析了LDH疾病的相关靶点基因,可进一步研究疾病的相关信号通路,基于信号通路挖掘相关治疗药物,同时两者结合,为杜仲-牛膝治疗LDH进一步的实验研究及临床应用提供参考思路。
参考文献
向 熙,司群超,成伟益,等.不同手术方式治疗腰椎间盘突出症的网状Meta分析[J].中国组织工程研究,2020,24(27):4398-4405. [百度学术]
洪炳杰,王道明,张沈屹婷.腰椎间盘突出症中医证素分布与组合规律文献研究[J].中医药临床杂志,2020,32(11):2109-2113. [百度学术]
秦晓宽,孙 凯,方圣杰,等.基于文献数据挖掘和网络药理学分析的腰椎间盘突出症用药规律及作用机制研究[J].中药新药与临床药理,2021,32(10):1490-1499. [百度学术]
侯宇龙,王晶石,王旭凯.牛膝治疗腰椎间盘突出症潜在作用机制的网络药理学分析[J].中国组织工程研究,2021,25(17):2734-2739. [百度学术]
王 猛,刘 刚,路聊东,等.基于网络药理学和分子对接探讨“杜仲-牛膝”药对治疗骨质疏松的分子机制[J].世界中医药,2021,16(20):2966-2974. [百度学术]
JINLONG RU,PENG LI,WANGJINAN,et al.TCMSP:a database of systems pharmacology for drug discovery from herbal medicines[J].J Cheminformatics,2014,6(1):13. [百度学术]
温明韬,梁学振,许 波,等.基于网络药理学探讨杜仲-牛膝治疗股骨头坏死的分子机制[J].西部中医药,2021,34(8):76-83. [百度学术]
LIUXIAOFENG,OUYANGSISHENG,YUBIAO,et al.PharmMapper Server:a web server for potential drug target identification via pharmacophore mapping approach[J].Nucleic Acids Res,2010,38(Web Server issue):W609-614. [百度学术]
WANG X,PAN C,GONG J,et al.Enhancing the Enrichment of Pharmacophore-Based Target Prediction for the Polypharmacological Profiles of Drugs[J].J Chem Inf Model,2016,56(6):1175-1183. [百度学术]
XIA WANG,SHENYIHANG,WANGSHIWEI,et al.PharmMapper 2017 update:a web server for potential drug target identification with a comprehensive target pharmacophore database[J].Nucleic Acids Res,2017,45(W1):W356-W360. [百度学术]
CONSORTIUMUNIPROT.UniProt:the universal protein knowledgebase in 2021[J]. Nucleic Acids Res,2021,49(D1):D480-D489. [百度学术]
WISHART DS,FEUNANG YD,GUO AC,et al.DrugBank 5.0:a major update to the DrugBank database for 2018[J].Nucleic Acids Res,2018,46(D1):D1074-D1082. [百度学术]
SZKLARCZYK D,GABLE AL,NASTOU KC,et al.The STRING database in 2021:customizable protein-protein networks,and functional characterization of user-uploaded gene/measurement sets[J].Nucleic Acids Res,2021,49(D1):D605-D612. [百度学术]
HUANG DA W,SHERMAN BT,LEMPICKI RA.Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources[J].Nat Protoc,2009,4(1):44-57. [百度学术]
司丽君,王 雪,王林林,等.槲皮素的抗炎免疫及部分机制研究[J].中国医药导报,2021,18(27):26-29,34. [百度学术]
郑泽陆,张瀚文,梁 栋,等.川芎-牛膝药对治疗膝骨关节炎的网络药理学探究[J].世界中西医结合杂志,2021,16(5):890-896. [百度学术]
蔡美云,庄文欣,吕 娥,等.山奈酚通过抑制p38 MAPK通路减轻6-羟多巴胺(6-OHDA)诱导的PC12细胞炎症[J].细胞与分子免疫学杂志,2020,36(7):583-589. [百度学术]
邹雪莲,胡 雯,雷子贤,等.山奈酚减轻人角质形成细胞氧化应激损伤[J].基础医学与临床,2021,41(8):1103-1108. [百度学术]
王甲任,顾留歌,秦 燕.基于网络药理学及分子对接研究穿心莲抗炎的分子机制[J].中国医药工业杂志,2021,52(4):519-533. [百度学术]
李若楠,洪 盼,郎伍营,等.β-胡萝卜素对脂多糖刺激的IPEC-J2细胞紧密连接蛋白表达的影响[J].中国免疫学杂志,2017,33(11):1611-1615. [百度学术]
BITTENCOURT-MERNAK MI,PINHEIRO NM,DA SILVA RC,et al.Effects of eugenol and dehydrodieugenol B from Nectandra leucantha against lipopolysaccharide(LPS)-Induced experimental acute lung inflammation[J].J Nat Prod,2021,84(8):2282-2294. [百度学术]
SANTANA FPR,DA SILVA RC,PONCI V,et al.Dehydrodieugenol improved lung inflammation in an asthma model by inhibiting the STAT3/SOCS3 and MAPK pathways[J].Biochem Pharmacol,2020,180:114175. [百度学术]
李洪涛,魏心然,孙晓伟,等.脊痛消胶囊对腰椎间盘突出症的疗效及对IL-1β、IL-6、TNF-α因子的影响[J].中医药信息,2020,37(2):95-99. [百度学术]
吕立江,毛凌宇,李景虎,等.杠杆定位手法结合脉冲电场对腰椎间盘突出症患者镇痛效应及IL-1β、TNF-α的影响[J].中国骨伤,2021,34(8):780-784. [百度学术]
张潇潇,姜 宏,张志刚,等.肿瘤坏死因子与腰椎间盘退变的相关 性研究进展[J].中国中医骨伤科杂志,2017,25(12):79-82. [百度学术]
WANG SJ,LIU C,SUN ZY,et al.IL-1p increases asporin expression via the NF-KB p65 pathway in nucleus pulposus cells dunng intervertebral disc degeneration[J].SciRep,2017,7(1):4112. [百度学术]
谢小妤,王海疆.红细胞分布宽度及血清ALT、AST、ALB水平与慢性乙型肝炎肝细胞炎症的关系分析[J].检验医学与临床,2020,17(12):1684-1687. [百度学术]
马 蓉,张新爱,杨 洋,等.腺苷酸环化酶相关蛋白2、纤维连接蛋白1表达与胃癌病理类型及预后关系探讨[J].陕西医学杂志,2021,50(2):241-244. [百度学术]
李玉华.Toll样受体信号转导通路与急性肺损伤的研究进展[J].中国免疫学杂志,2021,37(1):115-118,123. [百度学术]
WANG P F,XIONG X Y,CHEN J,et al.Function and mechanism of toll-like receptors in cerebral ischemic tolerance: from preconditioning to treatment[J].J Neuroinflammation,2015,12(1):80. [百度学术]
王链链,郭晓英.破骨细胞分化过程中的信号通路及信号因子的研究进展[J].中国骨质疏松杂志,2015(6):742-748. [百度学术]
左 斌,夏晓枫,车 彪,等.虎杖苷对小鼠腰椎间盘退变髓核细胞凋亡及SIRT1/mTOR通路的影响[J].中国组织工程研究,2021,25(35):5619-5625. [百度学术]
陈仁场,李念虎,管华鹏,等.椎间盘退变中的炎症反应与通路作用机制的研究进展[J].中华骨科杂志,2021,41(20):1509-1518. [百度学术]