摘要
通过使用TCMSP、ETCM和BATMAN-TCM数据库筛选陈皮、山楂的活性成分和作用靶点;通过Genecards、Disgenet和OMIM数据库获取HLP潜在作用靶点;通过Venny网站获取“陈皮-山楂”药对与疾病的交集靶点;使用Cytoscape软件构建活性成分-疾病-交集靶点网络图;使用String数据库和Cytoscape软件构建交集靶点蛋白质-蛋白质相互作用(PPI)网络;利用Metascape数据库对交集靶点进行基因本体(GO)生物功能分析及京都基因和基因组百科全书(KEGG)通路富集分析;最后利用AutoDock Vina软件对主要活性成分和关键靶点进行分子对接。
得到“陈皮-山楂”药对潜在活性成分71个、潜在作用靶点88个,其中核心靶点12个,GO生物过程(BP)作用于激素反应、细胞对激素刺激的反应等602个条目,GO细胞组分(CC)富集作用于转录调节复合物、囊泡腔等23个条目,GO分子功能(MF)富集作用于核受体活性、配体激活转录因子活性等47个条目,KEGG通路分析共得到127条通路,与脂肪细胞因子信号通路、PPAR信号通路和AMPK信号通路等有关。分子对接结果证实主要活性成分槲皮素、柚皮素与关键靶点PPARA、TP53结合能力良好。
高脂血症(Hyperlipidemia,HLP),通常表现为血浆或血清中总胆固醇(TC)、甘油三酯(TG)、低密度脂蛋白胆固醇(LDL-C)水平升高和高密度脂蛋白胆固醇(HDL-C)水平下
使用TCMSP数据库(https: //tcmspw.com/tcmsp.php),以口服生物利用度(OB)≥20%,类药性(DL)≥0.1作为阈值,对陈皮、山楂活性成分进行筛选,在TCMSP数据库的“related targets”栏找出活性成分对应的靶蛋白,通过UniProt数据库(https://www.uniprot.org/)查询靶点对应的基因名称。并检索ETCM(http://www.tcmip.cn/ETCM/index.php/Home/)数据库,筛选分值大于中位数的成分以及使用BATMAN-TCM(http://bionet.ncpsb.org.cn/batman-tcm/)数据库,筛选分值大于80的成分对其活性成分进行补充,合并三种数据库的基因。
在Genecards数据库(https://www.genecards.org/)、Disgenet数据库(https://www.disgenet.org/)和OMIM数据库(https://omim.org/)中以“Hyperlipidemia”为关键词检索高脂血症的相关靶点。为提高准确度,在Genecards中根据“Relevance score”值,筛选大于中位数的靶点,合并三种数据库的靶点,删除重复值,最终得到高血脂疾病作用靶点。
将取得的交集靶标导入String数据库(https://cn.string-db.org/),物种设为“Homo Sapiens”,最低相互作用阈值选择“highest confidence(.900)”,隐藏网络中游离的节点,下载tsv文件,将文件导入Cytoscape3.7.1中构建蛋白相互作用(PPI)网络。
将药物与疾病的交集靶点导入Metascape基因富集分析工具(https://metascape.org/),设置物种为“H.sapiens”,其他参数保持默认设置,对交集靶点分别进行基因本体(Gene Ontology,GO)的生物学过程(Biological Process,BP)、细胞组分(Cellular Component,CC)、分子功能(Molecular Functions,MF)的分析和京都基因和基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)通路富集分析。
通过三种数据库检索共得到陈皮活性成分36个,山楂活性成分36个。通过化合物信息数据库获取作用靶标,经UniProt数据库标准化后,最终得到陈皮55个活性成分作用靶点,山楂292个活性成分作用靶点。“陈皮-山楂”药对共71个活性成分,88个活性成分潜在作用靶点。
通过Venny 2.1.0在线绘图工具,最终得到成分与疾病交集靶点共88个,见

图1 陈皮-山楂活性成分与疾病交集靶点维恩图
左侧圈为陈皮-山楂共同活性成分作用靶点数量,右侧圈为HLP作用靶点数量,中间交集处为交集靶点数量
使用Cytoscape 3.7.1软件构建“陈皮-山楂”的“成分-疾病-交集靶点”网络,见

图2 成分-疾病-交集靶点网络图
图中外圈紫色三角为山楂活性成分,内圈粉色三角为陈皮活性成分,中间绿色菱形为交集靶点
Medicine | Compound | Average ShortestPath Length | BetweennessCentrality | Closeness Centrality | Degree |
---|---|---|---|---|---|
陈皮 | Naringenin(柚皮素) | 2.65957447 | 0.37310318 | 0.37600000 | 23 |
山楂 | Linoleic Acid(亚油酸) | 3.15602837 | 0.23848168 | 0.31685393 | 21 |
陈皮 | Nobiletin(川皮苷) | 2.91489362 | 0.11844473 | 0.34306569 | 14 |
山楂 | Quercetin(槲皮素) | 3.02836879 | 0.08202619 | 0.33021077 | 12 |
山楂 | Methylheptenone(甲基庚烯酮) | 3.43971631 | 0.10223312 | 0.29072165 | 10 |
陈皮 | Tangeretin(橘皮素) | 3.08510638 | 0.08539722 | 0.32413793 | 8 |
山楂 | Eburical(齿孔醛) | 2.94326241 | 0.07088254 | 0.33975904 | 8 |
山楂 | (-)-Epicatechin(表儿茶素) | 3.46808511 | 0.02050152 | 0.28834356 | 7 |
陈皮 | Hepta-3(3,5,6,7,8,3',4'-七甲氧基黄酮) | 3.02836879 | 0.02138338 | 0.33021077 | 6 |
陈皮 | Lauric Acid(月桂酸) | 3.36879433 | 0.06174753 | 0.29684211 | 6 |
山楂 | Ergosterol(麦角甾醇) | 3.04255319 | 0.04165902 | 0.32867133 | 6 |
陈皮 |
5,7-dihydroxy-2-(3-hydroxy-4-methoxyphenyl)chroman-4-one (橙皮素) | 3.46808511 | 0.00876946 | 0.28834356 | 5 |
陈皮 | DIBP(邻苯二甲酸二异丁酯) | 3.78014184 | 0.03085261 | 0.26454034 | 5 |
陈皮 | Hesperidin(橙皮苷) | 3.62411348 | 0.03421568 | 0.27592955 | 5 |
陈皮 | Beta-Myrcene(月桂烯) | 4.20567376 | 0.00216075 | 0.23777403 | 4 |
陈皮 | Perillaldehyde(紫苏醛) | 3.42553191 | 0.05005381 | 0.29192547 | 4 |
陈皮 | Limonene(柠檬烯) | 4.20567376 | 0.00216075 | 0.23777403 | 4 |
陈皮 | Caryophyllene(石竹烯) | 4.20567376 | 0.00216075 | 0.23777403 | 4 |
陈皮 | Alpha-Pinene[(±)-2-蒎烯)] | 4.20567376 | 0.00216075 | 0.23777403 | 4 |
陈皮 | Gamma-Terpinene(γ-松油烯) | 4.20567376 | 0.00216075 | 0.23777403 | 4 |
在String中导入交集靶点基因,物种设置为“Homo Sapiens”,最低阈值≧0.900,排除与其他蛋白无相互作用的孤立蛋白(2个),将从String数据库得到的数据导入Cytoscape 3.7.1软件中,建立可视化PPI网络,构建一个79个节点286条边的PPI网络图,其中,节点的大小和颜色深浅变化与度值(Degree)呈正相关,见

图3 交集靶点PPI网络图
Protein name | Degree | EPC | EcCentricity | Closeness | Betweenness |
---|---|---|---|---|---|
AKT1 | 26 | 17.739 | 0.31646 | 48.5 | 1124.11746 |
RXRA | 21 | 16.858 | 0.23734 | 44.16667 | 308.6323 |
NCOA1 | 21 | 16.637 | 0.18987 | 41.43333 | 198.71931 |
TP53 | 19 | 16.716 | 0.23734 | 43.83333 | 379.55181 |
PPARA | 19 | 16.648 | 0.23734 | 43.41667 | 472.63736 |
NCOA2 | 18 | 16.297 | 0.18987 | 39.93333 | 116.54937 |
HSP90AA1 | 17 | 15.657 | 0.23734 | 42.5 | 351.07983 |
TNF | 17 | 15.712 | 0.23734 | 43.08333 | 369.72758 |
ESR1 | 16 | 15.166 | 0.23734 | 41.91667 | 257.33722 |
INS | 16 | 15.402 | 0.23734 | 42.41667 | 377.19589 |
PPARG | 16 | 16.252 | 0.23734 | 42.5 | 250.59434 |
NR3C1 | 16 | 15.699 | 0.23734 | 42.08333 | 179.54841 |
SREBF1 | 15 | 14.409 | 0.23734 | 40.83333 | 245.31159 |
MAPK1 | 14 | 15.395 | 0.23734 | 41.25 | 297.9127 |
PTGS2 | 13 | 11.966 | 0.23734 | 39 | 457.57823 |
PIK3CA | 13 | 12.537 | 0.23734 | 38.83333 | 189.15269 |
APOB | 13 | 11.524 | 0.18987 | 38.31667 | 262.6285 |
PRKACA | 12 | 12.400 | 0.23734 | 38.16667 | 547.11396 |
AR | 12 | 12.538 | 0.23734 | 38.5 | 126.421 |
RXRG | 11 | 12.817 | 0.23734 | 37.83333 | 53.64969 |
ADIPOQ | 11 | 13.674 | 0.18987 | 37.98333 | 48.73605 |
TLR4 | 10 | 11.264 | 0.23734 | 37.41667 | 96.46271 |
HMGCR | 9 | 10.958 | 0.18987 | 35.26667 | 75.90774 |
MMP9 | 9 | 9.461 | 0.23734 | 35.66667 | 115.81856 |
ESR2 | 9 | 12.016 | 0.23734 | 37.16667 | 55.99061 |
CAT | 9 | 9.544 | 0.23734 | 36.83333 | 407.94433 |
RARG | 8 | 11.258 | 0.23734 | 36.41667 | 35.78216 |
VDR | 7 | 10.931 | 0.18987 | 34.48333 | 9.68526 |
NOS2 | 7 | 10.762 | 0.23734 | 37.25 | 33.12429 |
FASN | 7 | 9.300 | 0.18987 | 34.4 | 149.56719 |

图4 陈皮-山楂核心富集模块
GO生物功能分析涉及BP、CC、MF,将各自基因数排名前10的绘制成柱状图,见

图5 交集靶点GO生物功能分析图

图6 交集靶点KEGG富集分析图
利用Cytoscape 3.7.1软件构建陈皮-山楂活性成分-靶点-信号通路网络图,见

图7 活性成分-靶点-信号通路网络图
利用Autodock Vina对排名靠前的活性成分与靶点进行分子对接。通常,结合能≤-5.0 kcal/mol说明活性成分与作用靶点之间有较好的结合活性,见
主要活性成分 | 关键靶标 | 结合能(kcal/mol) |
---|---|---|
柚皮素
| AKT1 | -5.9 |
RXRA | -8.6 | |
NCOA1 | -9.6 | |
TP53 | -9.9 | |
PPARA | -10.8 | |
川皮苷
| AKT1 | -5.5 |
RXRA | -8.3 | |
NCOA1 | -9.1 | |
TP53 | -9.7 | |
PPARA | -11.0 | |
槲皮素
| AKT1 | -6.2 |
RXRA | -9.4 | |
NCOA1 | -10.4 | |
TP53 | -9.6 | |
PPARA | -11.6 | |
甲基庚烯酮
| AKT1 | -2.9 |
RXRA | -4.1 | |
NCOA1 | -5.2 | |
TP53 | -5.7 | |
PPARA | -5.8 | |
橘皮素
| AKT1 | -5.6 |
RXRA | -8.3 | |
NCOA1 | -9.1 | |
TP53 | -9.4 | |
PPARA | -10.4 |


图8 活性成分与关键靶点的分子对接
注: 上面左为槲皮素与PPARA,右为柚皮素与TP53;下面左为槲皮素与TP53,右为柚皮素与PPARA
随着人们生活水平的提高,高脂血症的发病率也越来越高,同时高脂血症也会引起一系列疾病的发生,如糖尿
本研究基于网络药理学筛选出“陈皮-山楂”药对71个活性成分及88个潜在靶点;与HPL取得交集后共获得88个交集靶点,通过筛选最后获得20个主要活性成分,包括柚皮素、川皮苷、槲皮素、甲基庚烯酮与橘皮素等;柚皮素具有降血脂、抗炎、保肝的作
GO分析中BP分析表明“陈皮-山楂”药对可能通过激素反应、细胞对激素刺激的反应、细胞内受体信号通路、细胞对脂质的反应、对类固醇激素的反应、小分子代谢过程的调节等作用于HPL;CC分析结果表明“陈皮-山楂”药对可能通过转录调节复合物、RNA聚合酶II转录调节复合物、囊泡腔、内质网腔、细胞质小泡腔、细胞质核周区和受体复合物等功能调控;MF分析结果表明“陈皮-山楂”药对可能通过核受体活性、配体激活转录因子活性、转录因子结合、RNA聚合酶II特异性DNA结合转录因子结合、核受体结合、染色质结合、蛋白激酶结合和酶结合等分子功能调控。有研究表明,莲子抗性淀粉低剂量组能够显著调控与甘油代谢和类固醇激素生物合成相关的基
KEGG富集分析结果表明“陈皮-山楂”药对主要富集在脂肪细胞因子信号通路、催乳素信号通路、甲状腺激素信号通路、内分泌抗性、雌激素信号通路、VEGF信号通路、PPAR信号通路和AMPK信号通路等信号通路上。实验表明,青钱柳多糖通过调节脂肪细胞因子信号通路和AMPK信号通路可以降低高脂小鼠肝脏全基因组DNA甲基化水平,起到降脂作
分子对接对“陈皮-山楂”药对治疗HPL排名靠前的5个活性成分与排名前5的核心靶点进行验证。结果表明,柚皮素、川皮苷、槲皮素、甲基庚烯酮和橘皮素与AKT1、RXRA、NCOA1、TP53和PPARA靶点具有良好的结合亲和力。表明“陈皮-山楂”药对可能通过柚皮素、川皮苷、槲皮素等活性成分作用于AKT1、RXRA、NCOA1等关键靶点起到治疗HPL的作用。
综上,本研究基于网络药理学对“陈皮-山楂”药对治疗HPL的有效成分、潜在靶点以及可能的通路进行挖掘,结果表明“陈皮-山楂”药对可能通过柚皮素、川皮苷和槲皮素等影响AKT1、RXRA、NCOA1等关键蛋白表达从而调节脂肪细胞因子信号通路、VEGF信号通路、PPAR信号通路和AMPK信号通路等来调节脂质代谢进而起到降低血脂的作用,分子对接结果显示部分有效成分与关键靶点具有良好的结合亲和力。反映了中药治疗HPL多成分、多靶点、多通路的治疗特点,未来可对关键成分及核心靶点进行体内外实验验证。
参考文献补充并核对卷号、期号、页码起止页
[1]王 琛,郭姗姗,宋云龙,等.心脉通胶囊合非诺贝特对饮食控制疗法效果不佳的高脂血症患者血脂、血液流变学和血管内皮功能的影响[J].现代生物医学进展,2022,22(14):2695-2699.
[2]杨园园,吴圣贤,赵 颖,等.从湿浊论治原发性高脂血症[J].中华中医药杂志,2019,34(12):5602-5604.
[3]安冬青,吴宗贵,梁 春,等.血脂异常中西医结合诊疗专家共识[J].中国全科医学,2017,20(3):262-269.
[4]俞静静,苏 洁,颜美秋,等.陈皮降脂药效与黄酮类成分的相关性研究[J].中国中药杂志,2019,44(15):3335-3342.
[5]于舒雁,田 硕.山楂中黄酮的提取工艺及其调节糖脂代谢的活性研究[J].中国调味品,2022,47(2):182-185.
[6]ZHOU X, ZHANG W, LIU X,et al. Interrelationship between diabetes and periodontitis: role of hyperlipidemia[J]. Arch Oral Biol. 2015,60(4):667-74.
[7]GORDTS PLSM, ESKO JD. The heparan sulfate proteoglycan grip on hyperlipidemia and atherosclerosis[J]. Matrix Biol. 2018,71(72):262-282.
[8]FLORA GD, NAYAK MK. A Brief Review of Cardiovascular Diseases, Associated Risk Factors and Current Treatment Regimes[J]. Curr Pharm Des. 2019,25(38):4063-4084.
[9]OKERSON T, PATEL J, DIMARIO S, et al. Effect of 2013 ACC/AHA Blood Cholesterol Guidelines on Statin Treatment Patterns and Low-Density Lipoprotein Cholesterol in Atherosclerotic Cardiovascular Disease Patients[J]. J Am Heart Assoc. 2017,6(3):e004909.
[10]陈克永.瑞舒伐他汀联合烟酸治疗对高脂血症患者血脂、颈动脉内中膜厚度及踝肱指数的影响[J].中国实用医药,2019,14(3):1-3.
[11]CAMACHO-MUÑOZ D, KIEZEL-TSUGUNOVA M, KISS O, et al. Omega-3 carboxylic acids and fenofibrate differentially alter plasma lipid mediators in patients with non-alcoholic fatty liver disease[J]. FASEB J. 2021,35(11):e21976.
[12]OSCARSSON J, ÖNNERHAG K, RISÉRUS U, et al. Effects of free omega-3 carboxylic acids and fenofibrate on liver fat content in patients with hypertriglyceridemia and non-alcoholic fatty liver disease: A double-blind, randomized, placebo-controlled study[J]. J Clin Lipidol. 2018,12(6):1390-1403.
[13]HAO P, JIANG F, CHENG J, et al. Traditional Chinese Medicine for Cardiovascular Disease: Evidence and Potential Mechanisms[J]. J Am Coll Cardiol. 2017,69(24):2952-2966.
[14]张思琪. 基于数据挖掘分析游运舸治疗代谢相关脂肪性肝病的用药规律[D].长沙:湖南中医药大学,2022.
[15]CUI W, HE Z, ZHANG Y, et al. Naringenin Cocrystals Prepared by Solution Crystallization Method for Improving Bioavailability and Anti-hyperlipidemia Effects[J]. AAPS PharmSciTech. 2019,20(3):115.
[16]马兆臣,陈奎奎,潘琦雪,等.丹荷颗粒治疗高脂血症质量标志物发现研究[J].分析测试学报,2021,40(1):43-49.
[17]LIU H, ZHAO H, CHE J, et al. Naringenin Protects against Hypertension by Regulating Lipid Disorder and Oxidative Stress in a Rat Model[J]. Kidney Blood Press Res. 2022,47(6):423-432.
[18]PENG L, ZHANG Q, ZHANG Y, et al. Effect of tartary buckwheat, rutin, and quercetin on lipid metabolism in rats during high dietary fat intake[J]. Food Sci Nutr. 2019,8(1):199-213.
[19]CHEN PY, CHAO TY, HSU HJ, et al. The Lipid-Modulating Effect of Tangeretin on the Inhibition of Angiopoietin-like 3 (ANGPTL3) Gene Expression through Regulation of LXRα Activation in Hepatic Cells[J]. Int J Mol Sci. 2021,22(18):9853.
[20]HE S, XIONG Z, LI L, et al.Lotus seed resistant starch ameliorates high-fat diet induced hyperlipidemia by fatty acid degradation and glycerolipid metabolism pathways in mouse liver. Int J Biol Macromol[J]. 2022,215:79-91.
[21]张彩虹,刘宏峰,朱 斌.微粒体三酰甘油转运蛋白研究进展[J].国际检验医学杂志,2017,38(10):1384-1386.
[22]INOUE J, IKEDA S, KANAYAMA T, et al. The flavonoid derivative 4'-nitro-6-hydroxyflavone suppresses the activity of HNF4α and stimulates the degradation of HNF4α protein through the activation of AMPK[J]. Biosci Biotechnol Biochem. 2017,81(8):1548-1552.
[23]YANG Z, ZHAO J, LI J, et al. Genome-wide DNA methylation profiling of high-fat emulsion-induced hyperlipidemia mice intervened by a polysaccharide from Cyclocarya paliurus (Batal) Iljinskaja[J]. Food Chem Toxicol. 2021,152:112230.
[24]HU N, CHEN C, WANG J, et al. Atorvastatin Ester Regulates Lipid Metabolism in Hyperlipidemia Rats via the PPAR-signaling Pathway and HMGCR Expression in the Liver[J]. Int J Mol Sci. 2021,22(20):11107.
[25]周 玉,寿鑫甜,柴若宁,等.基于CiteSpace的PPAR与脂代谢相关可视化研究[J/OL].心脏杂志,2022(6):699-704[2022-10-14]. http://kns.cnki.net/kcms/detail/61.1268.R.20220914.1749.015.html.
[26]JI J, XUE TF, GUO XD, et al. Antagonizing peroxisome proliferator-activated receptor γ facilitates M1-to-M2 shift of microglia by enhancing autophagy via the LKB1-AMPK signaling pathway[J]. Aging Cell. 2018,17(4):e12774.
[27]史 婧,关 键.PPARγ相关自噬在急性肺损伤中的作用机制研究[J].临床肺科杂志,2020,25(12):1817-1821.
[28]黄小强,李宣宣,吴水生,等.泽泻三萜对小鼠急性肺损伤的保护作用[J].中成药,2022,44(9):2979-2984.
[29]SU XIN,CHEN XIANG,WANG BIN. Relationship between the development of hyperlipidemia in hypothyroidism patients.[J]. Molecular biology reports,2022.
SU X, CHEN X, WANG B. Relationship between the development of hyperlipidemia in hypothyroidism patients. Mol Biol Rep[J]. 2022,49(11):11025-11035.