Mechanism of Cistanches Herba in Treatment of Parkinson’s Disease Based on Network Pharmacology and Molecular Docking
-
摘要: 基于网络药理学及分子对接对肉苁蓉治疗帕金森疾病的作用机制进行研究。利用 TCMSP筛选肉苁蓉中的活性成分作为候选化合物，再利用PubChem和sea数据库获取肉苁蓉的靶点信息。在GeneCard数据库中获取帕金森病靶点信息，通过比对中药和疾病靶点，得到共有靶点基因信息。运用STRING进行靶点蛋白互作（protein-protein interaction，PPI）网络的构建，在DAVID6.8数据库中对关键靶基因进行基因本体（gene ontology，GO）富集和基因组百科全书（Kyoto encyclopedia of genes and geomes，KEGG）通路富集得出治疗靶基因，利用cytoscape软件构建中药-化合物-靶点网络图。并将该网络中的基因靶点以degree值高低进行排序，选择前5个关键靶点，最后利用AutoDock Vina工具将5个关键靶点与肉苁蓉有效活性成分进行分子对接。结果表明：从肉苁蓉中共筛选得到6个有效活性成分，分别为β-谷甾醇、花生四烯酸、脱水内酯、丁香树脂二甲醚、槲皮素、枫马京。交集获得308个预测靶点，将其进行PPI分析，并从中筛选出38个关键靶基因，对关键靶基因进行GO与KEGG通路富集。最终发现，GO富集中生物过程（biologica process，BP）包括细胞内信号转导等；细胞成分（cellular component，CC）分析主要为蛋白质复合物、线粒体等；分子功能（molecular function，MF）结果主要为腺嘌呤核苷三磷酸结合、蛋白质复合物结合等。KEGG通路富集主要包括磷脂酰肌醇-3-激酶（phosphatidylin-ositol-3-kinase，PI3K）/蛋白激酶B（Protein Kinase B，PKB/AKT）、微型核醣核酸（microRNAs，miRNAs）、缺氧诱导因子-1（hypoxia inducible factor-1，HIF-1）、肿瘤坏死因子（tumor necrosis factor，TNF）、叉形头转录因子O亚型（forkhead box O，FOXO）等相关信号通路。筛选前5个关键靶点分别为MAPK8、ESR1、SRC、JAK2、AKT1。分子对接结果显示肉苁蓉6个有效活性成分与5个关键靶点具有较好的结合活性，结合能均小于−5 kcal/mol。本研究得出肉苁蓉通过保护线粒体功能、调节蛋白质活动、调节自噬、抑制神经炎症反应等治疗帕金森病，发挥多成分、多通路、多靶点协同作用。Abstract: The study aimed to explore the mechanism of Cistanches Herba in treatment of Parkinson’s disease based on network pharmacology and molecular docking. The active components in Cistanches Herba were selected as candidate compounds by the TCMSP. And the target information of Cistanches Herba was obtained by PubChem and sea database. The target information of Parkinson’s disease was obtained by GeneCard database and the common target gene information was obtained by comparing traditional Chinese medicine and disease targets. PPI network was constructed by means of STRING. The key target genes were enriched by GO and KEGG pathway in DAVID 6.8 database. The Cytoscape software was applied to construct the network diagram of traditional Chinese medicine-chemical compound-target. The gene targets in the network were sorted according to the degree value, the first five key targets were selected. Finally, the five key targets were in connection with the effective active components of Cistanches Herba by using AutoDock Vina tool. In the end, the results showed that 6 active components were screened from Cistanches Herba beta-sitosterol, arachidonate, suchilactone, Yangambin, quercetin, Marckine. 308 prediction targets were obtained by intersection, which were analyzed by PPI. 38 key target genes were filtered out, and the key target genes were enriched by GO and KEGG pathway. Finally, it was found that GO rich BP includes intracellular signal transduction; CC analysis mainly included protein complex, mitochondrion; MF results mainly included ATP binding, protein complex binding. KEGG pathway enrichment mainly included PI3K/Akt, miRNAs, HIF-1, TNF, FOXO signal pathways, etc. The first five key targets were MAPK8, ESR1, SRC, JAK2 and AKT1. The results of molecular docking showed that the six active components of Cistanches Herba had good binding activity with five key targets. This study found that Cistanches Herba could play a multi-component, multi-channel and multi-target synergistic role in the treatment of Parkinson's disease by protecting mitochondrial function, regulating protein activity, regulating autophagy and inhibiting neuroinflammatory response.
-
Keywords:
- Cistanches Herba /
- Parkinson's disease /
- network pharmacology /
- molecular docking
-
帕金森病（Parkinson’s disease，PD）是一种常见的中老年神经系统退行性疾病，临床上以运动迟缓、静止性震颤、肌强直、姿势平衡障碍等运动症状以及睡眠障碍、嗅觉障碍、自主神经功能障碍、认知和精神障碍等非运动症状为显著特征[1]。PD主要病理特征为黑质中多巴胺能神经元的丢失和由α-突触核蛋白（α-Synuclein，α-Syn）组成的路易小体的存在，但发病机制尚未完全明确。目前我国65岁以上人群中PD发病率为1.7%，且随着老龄化程度加剧，据推测我国PD患病人数到2030年将达到500万人，约占到全球PD患病人数的一半[2]。目前PD的药物治疗多以左旋多巴制剂为主，虽然能够在一定程度上改善症状，但长期使用存在药效减退、剂末现象、异动症等运动并发症，且无法从根本上阻止疾病进展。
肉苁蓉，俗称大芸，始载于《神农本草经》，列属上品，味甘、咸，性温，归肾、大肠经，中医认为其具有补肾阳，益精血，润肠通便之效。研究表明，肉苁蓉具有抗凋亡、抗氧化、抗衰老、抗脑缺血再灌注损伤等神经保护作用[3]，其化学成分有苯苷类、环烯醚萜及其苷类、木脂素及其苷类、多糖以及单萜苷类、生物碱等，其中又以苯乙醇苷类和多糖为主[4]，药理学研究发现肉苁蓉在PD及阿尔兹海默病等神经退行性疾病中发挥的作用显著[5-6]。据统计，肉苁蓉在抗衰老延年类古方中的出现率仅次于人参[7]。目前，肉苁蓉常被应用于PD的临床治疗，但其发挥疗效的具体机制并不明确，本文主要从网络药理学和分子对接的角度，利用相关数据库筛选肉苁蓉的主要活性化合物，并将其作用靶点与PD靶点比对，得出关键作用靶点，探析其治疗PD的物质基础与作用机制，以期为深入研究提供指导。
1. 资料与方法
1.1 有效化合物
本研究使用中药系统药理数据库（Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform TCMSP，http://lsp.nwu.edu.cn/tcmsp.php）检索肉苁蓉成分，口服生物利用度（oral bioavailability，OB）≥30%，类药性（drug-likeness，DL）≥0.18[8]为阈值，因血脑屏障功能障碍与PD病程有关[9]，故设定血脑屏障通透性（blood brain barrier，BBB）≥−0.85，通过以上限定条件筛选肉苁蓉中的活性成分作为有效化合物。
1.2 查询PD相关疾病靶点
利用PubChem数据库（www.ncbi.nlm.nih.gov/pccompound/）、sea数据库（http://sea.bkslab.org/）获得肉苁蓉所有靶点基因名称。应用GeneCard数据库（https://www.genecards.org/）以“Parkinson’s disease”为检索词获取PD靶点基因名称。过比对中药和疾病靶点，利用VENNY2.1软件作交集饼图得到PD与肉苁蓉共有靶点基因信息。
1.3 靶点蛋白互作（protein-protein interaction，PPI）网络的构建
将1.1及1.2获得的中药、疾病共有基因输入STRING（https://string-db.org/），选择“Homo sapiens”，将最小相互作用阈值设为中等“Highest confidence>0.9”，其余设置为默认值[10]，进行PPI网络构建。
1.4 基因本体（gene ontology，GO）富集和基因组百科全书（Kyoto encyclopedia of genes and geomes，KEGG）通路富集
筛选关键靶基因（大于2倍degree），将筛选出的靶基因导入DAVID6.8数据库（https://david.ncifcrf.gov/summary.jsp），限定物种为“Homo Sapiens”，以P<0.05为阈值，选择其中BP、CC、MF 3对关键靶基因进行GO功能富集分析，筛选具有显著差异的生物过程及可靠的靶点通路，按照P值从小到大排序，前20个作为重要通路用Omicshare （http://www.omicshare.com/） 制作成可视化气泡图。通过KEGG通路富集分析靶点信号通路，以P<0.05为阈值，依据count值排序，选取前20条信号通路选用Omicshare （ http://www.omicshare.com/） 制作成可视化气泡图。
1.5 构建“肉苁蓉-化合物-靶点”网络图
选取前20条KEGG信号通路上的基因作为治疗靶基因，将筛选的治疗靶基因对应的化合物和中药导入到cytoscape软件中，构建“肉苁蓉-化合物-靶点”网络图。
1.6 “活性成分-靶点”分子对接
选取网络药理学筛选出的肉苁蓉有效化合物与帕金森病的关键靶点，利用AutoDock Vina工具进行分子对接，了解靶点蛋白质与活性成分的结合方式。
2. 结果与分析
2.1 肉苁蓉有效化合物及靶点基因筛选
使用TCMSP网站筛选出肉苁蓉中药共6个有效化合物，分别为beta-sitosterol（β-谷甾醇）、arachidonate（花生四烯酸）、suchilactone（脱水内酯）、Yangambin（丁香树脂二甲醚）、quercetin（槲皮素）、Marckine（枫马京）。见表1。根据有效成分筛选出435个靶点基因名称。
表 1 肉苁蓉有效成分Table 1. Effective components of Cistanches Herba编号 化合物名称 OB（%） DL BBB MOL000358 beta-sitosterol 36.91 0.75 0.99 MOL005320 arachidonate 45.57 0.20 0.58 MOL005384 suchilactone 57.52 0.56 0.28 MOL007563 Yangambin 57.53 0.81 0.01 MOL000098 quercetin 46.43 0.28 -0.77 MOL008871 Marckine 37.05 0.69 -0.10 2.2 帕金森病靶点基因筛选
检索GeneCard数据库（https://www.genecards.org/）得到PD相关基因，共计7846个，如LRRK2、SNCA、PRKN、PARK7、GBA、PINK1、MAPT等。
2.3 肉苁蓉活性成分治疗帕金森病预测靶点的筛选
通过TCMSP收集肉苁蓉6个活性化合物对应靶点，去重得到435个蛋白靶标，与GeneCard数据库去重所得7846个蛋白靶标进行比对取交集，运用VENNY2.1软件对疾 病靶点与药物靶点整理后得到308个交集靶点，见图1。蓝色代表PD靶点，黄色代表肉苁蓉靶点，重叠部分为共同靶点。
![]()
图 1 帕金森病-肉苁蓉靶点交集韦恩图Figure 1. Wayne diagram of target intersection of Parkinson’s disease-Cistanches Herba2.4 PPI网络的构建及关键靶点筛选
将上述筛选出的靶基因导入STRING（https://string-db.org/）数据库，物种选“Homo sapiens”进行PPI分析，得到肉苁蓉作用于PD的PPI网络图，见图2。共得到283个节点，3509条边。
2.5 GO功能富集分析
筛选关键靶基因（大于2倍degree），将筛选出的靶基因导入DAVID6.8数据库，限定物种为“Homo Sapiens”，对38个靶基因进行GO富集分析。结果包括生物过程（biologica process，BP）305条、细胞成分（cellular component，CC）25条、分子功能（molecular function，MF） 69条，按照P值从小到大排序，前20个重要通路用Omicshare作气泡图进行可视化分析，如图3。图中气泡越大富集基因个数越多，气泡颜色越红，P值越小，显著性越高。图3A中BP包括RNA聚合酶II启动子的转录负调控、以DNA作为模板转录、细胞增殖、凋亡过程的正调控、细胞内信号转导等；图3B可见CC分析中包括膜筏、神经元突触、细胞膜、细胞核膜、蛋白质复合物、线粒体等；图3C可见MF富集结果主要为腺嘌呤核苷三磷酸结合、蛋白异源二聚体活性、蛋白质复合物结合、蛋白质酪氨酸激酶活性等。
![]()
图 3 GO功能富集分析注：A：GO-BP；B：GO-CC；C：GO-MF。Figure 3. Gene ontology function enrichment analysis2.6 KEGG通路分析及构建“肉苁蓉-化合物-靶点“网络图
通过KEGG通路富集分析靶点信号通路，以P<0.05为阈值，依据count值排序，选取前20条信号通路，其中主要包括PI3K/AKT信号通路、miRNAs调控、HIF-1、TNF、FOXO信号通路等，依据结果绘制气泡图，见图4。图中气泡越大富集基因个数越多，气泡颜色越红，P值越小，显著性越高。选取前20条KEGG信号通路上的基因作为治疗靶基因，将91个治疗靶基因对应的化合物和中药导入到cytoscape软件中，构建中药-化合物-靶点网络图，见图5。网络中的节点表示靶标或化合物，而边则表示活性成分与PD靶标之间的相互作用，颜色越深，degree值越高。
![]()
图 4 KEGG信号通路富集分析Figure 4. Enrichment analysis of Kyoto encyclopedia of genes and geomes signal pathway2.7 “活性成分-靶点”分子对接
从中药-化合物-靶点网络中将基因靶点以degree值高低进行排序，选择前5个关键靶点，分别为丝裂原活化蛋白激酶8（mitogen-activated protein kinase 8，MAPK8）、雌激素受体1（estrogen receptor alpha 1，ESR1）、类固醇受体辅助活化因子（steroid receptor coactivator，SRC）、蛋白酪氨酸激酶2（janus kinase signal transducers 2，JAK2）、丝氨酸/苏氨酸激酶1（serine/threonine kinases 1，AKT1）。利用AutoDock Vina[11]分别与肉苁蓉的6个活性成分（β-谷甾醇、花生四烯酸、脱水内酯、丁香树脂二甲醚、槲皮素和枫马京）进行分子对接。结合能<−4.25 kcal/mol提示配体小分子与受体蛋白之间有一定的结合活性；结合能<−5.0 kcal/mol表明二者之间有较好的结合活性；结合能<−7.0 kcal/mol表示配体与受体具有强烈的结合活性[12]。分子对接结果见表2，β-谷甾醇、花生四烯酸、脱水内酯、丁香树脂二甲醚、槲皮素和枫马京6种化合物与关键作用靶点MAPK8、ESR1、SRC、JAK2和AKT1的结合能均小于−5 kcal/mol，说明药物与靶点的结合活性较好[10]。其中β-谷甾醇、脱水内酯和丁香树脂二甲醚与JAK2、MAPK8和SRC，以及槲皮素和枫马京与ESR1、JAK2、MAPK8和SRC的结合能均小于−7 kcal/mol，说明药物与靶点之间具有较强的结合活性。随后利用PyMol（www.pymol.org）对6个肉苁蓉有效成分与其结合活性最高靶点的分子对接结果进行可视化分析，如图6所示。其中β-谷甾醇与活性位点附近的ASN-397形成氢键结合到SRC，花生四烯酸与活性位点附近的ARG-394形成氢键结合到ESR1，脱水内酯与活性位点附近的SER-936和ARG-980形成氢键结合到JAK2，丁香树脂二甲醚与活性位点附近的ARG-980形成氢键结合到JAK2，槲皮素与活性位点附近的THR-338、GLU-339、MET-341和ASP-404形成氢键结合到SRC，枫马京与活性位点附近的PRO-361、GLN-362、GLU517、PHE520和SER-522形成氢键结合到SRC。
表 2 肉苁蓉有效活性成分与对应核心靶点的分子对接Table 2. Molecular docking of effective active components of Cistanches Herba with corresponding core targets有效活性成分 MAPK8 ESR1 SRC JAK2 AKT1 beta-sitosterol（β-谷甾醇） −7.8 −6.8 −8.1 −7.4 −6.3 arachidonate（花生四烯酸） −5.4 −6.3 −6.3 −6.1 −5.3 suchilactone（脱水内酯） −7.9 −5.9 −7.7 −8.6 −5.9 Yangambin（丁香树脂二甲醚） −7.4 −6.2 −7.5 −7.9 −6.2 quercetin（槲皮素） −8.1 −7.2 −9.2 −9 −6.6 Marckine（枫马京） −9 −7.9 −10.2 −10 −7.2 注：单位：结合能（kcal/mol）。 ![]()
图 6 分子对接模式图注：A：β-谷甾醇和SRC的对接图；B：花生四烯酸和ESR1的对接图；C：脱水内酯和JAK2的对接图；D：丁香树脂二甲醚和JAK2的对接图；E：槲皮素和SRC的对接图；F：枫马京和SRC的对接图。Figure 6. Molecular docking mode3. 讨论
研究表明肉苁蓉主要成分苯乙醇苷可通过抑制Caspase-3的活性，对抗1-甲基-4-苯基吡啶离子 （1-methyl-4- phylpyridiniumion，MPP+） 导致的多巴胺能神经细胞凋亡[13]。研究也发现松果菊苷作为肉苁蓉活性成分之一，能抑制PD小鼠黑质多巴胺能神经元的减少，且通过GFRa1/AKT信号通路上调胶质细胞源性神经营养因子（glial cell-derived neurotrophic factor，GDNF）表达[14]，通过ROS/ATF3/CHOP通路抑制神经细胞凋亡[15]，增加纹状体多巴胺神经递质含量，改善PD小鼠的行为学异常，表现出良好的神经保护和修复作用[16]。但关于肉苁蓉其他有效成分在治疗PD中的作用机制仍未探究。
本研究通过筛选发现肉苁蓉的有效活性成分一共6种，分别是beta-sitosterol（β-谷甾醇）、arachidonate（花生四烯酸）、suchilactone（脱水内酯）、Yangambin（丁香树脂二甲醚）、quercetin（槲皮素）、Marckine（枫马京）。研究发现β-谷甾醇具有抗氧化、抑制细胞凋亡等作用[17]，而氧化应激损伤是PD发病的关键环节之一[18]，推测β-谷甾醇可能通过抗氧化应激治疗PD。花生四烯酸可以通过激活过氧化物酶增殖体激活受体（peroxisome proliferators-activated receptors，PPAR）及促进新蛋白合成来增强神经组织内部的抗氧化系统，发挥抗氧化应激的作用，保护神经元细胞[19]。有研究发现槲皮素能够通过抑制AGEs/RAGE/NF-κB炎症通路，从而降低脑组织炎性因子的生成[20]。而另一项研究发现中药提取物可通过抑制晚期糖基化终末产物（advanced glycation end products，AGEs），阻断晚期糖基化终产物受体（receptor of advanced glycation endproducts，RAGE）表达，从而抑制人神经母细胞瘤细胞（human neuroblastoma cells，SH-SY5Y）中甲基乙二醛（methylglyoxal，MGO）触发的核因子κB（nuclear factor kappa-B，NF-κB）激活和易位，而已知MGO是一种二羰基代谢物，与神经退行性疾病的发展有关，提示可通过调节AGEs/RAGE/NF-κB通路减轻反应性二羰基化合物的神经毒性[21]。此外，另一项研究发现具有高氧化还原活性的铁可催化大脑氧化反应，还参与产生高毒性自由基（Fenton反应），从而对细胞造成氧化应激损伤，最终导致多巴胺能神经元死亡，而槲皮素可充当强有力的铁螯合剂发挥抗氧化作用，表现出一定的神经保护作用[22]。
GO功能富集分析表明肉苁蓉具有促进ATP合成、保护线粒体呼吸链的功能完整的作用，而线粒体功能障碍是引起多巴胺能神经元变性、死亡机制之一[23]。肉苁蓉还可以调节蛋白质的合成、折叠、降解等活动，维持“蛋白质内稳态”。而研究表明蛋白质构象异常或过度沉积可抑制突触功能，干扰细胞内信号传导，导致由泛素蛋白酶体系统介导的蛋白质降解功能紊乱，最终导致多巴胺能神经细胞死亡[24]。
通过KEGG富集得到活性成分与PD靶标之间相互作用的主要信号通路，包括PI3K/AKT信号通路、miRNAs调控、HIF-1、TNF及FOXO信号通路等，表明肉苁蓉通过多条通路作用于PD发病的调控过程。AKT，又名蛋白激酶B，是一种促进细胞存活的激酶，其存在于磷脂酰肌醇-3-激酶（phosphatidylin-ositol-3-kinase，PI3K）信号途径中。PI3K激活后的产物与AKT的PH结构域相结合，使其Ser473和Thr308位点磷酸化，从而激活AKT。研究发现PD患者黑质致密部（substantia nigra compacta，SNc）中总AKT水平及活性磷酸化AKT水平均显著降低[25]。林瑶等[26]发现肉苁蓉可使PD大鼠黑质纹状体中PI3K、ATK的表达升高及其下游凋亡相关蛋白Bcl-2表达升高、Bax表达降低，提示肉苁蓉通过激活PI3K/AKT通路，从而升高Bcl-2/Bax比值，减少神经细胞凋亡来发挥神经保护作用。糖原合成酶激酶-3β（glycogen synthase kinase-3β，GSK-3β）活性增高，诱导多巴胺能神经元凋亡，而研究发现载有碱性成纤维细胞生长因子（basic fibroblast growth factor，bFGF）的脂质体通过激活PI3K/AKT通路增加GSK-3β的磷酸化，从而抑制GSK-3β活性保护神经元[27]。此外，哺乳动物雷帕霉素靶蛋白（mammalian target of rapamycin，mTOR）作为PI3K/AKT下游信号分子，是自噬过程的主要调节因子，而α-Syn的聚集是PD发病机制的核心之一，激活PI3K/ATK/mTOR通路诱发自噬途径使聚集的α-Syn降解，从而发挥治疗PD的作用[28]。由此可见肉苁蓉有效成分可能通过激活PI3K/AKT信号通路促进神经细胞存活发挥治疗PD的作用。同样涉及α-Syn聚集导致PD的发病机制，Pino等[29]发现叉形头转录因子O亚型（forkhead box O，FOXO）家族重要的自噬调控因子FOXO3有效地减少可溶性α-Syn含量，且FOXO3控制神经元细胞中的自噬流量。由此可以推测，肉苁蓉有效成分可能通过FOXO信号通路减少α-Syn的积累保护多巴胺神经元。
miRNAs作为一簇微小的内源性非编码小RNA分子，能够通过序列特异性的方式来调节靶基因表达，参与调控多巴胺神经元增殖、凋亡及自噬等生物过程[30]。研究发现α-Syn（其编码基因为SNCA）的异常表达会诱导多巴胺神经元凋亡，而miRNA-7及miRNA-153可降低PD模型中SNCA的表达水平，改善PD症状[31]；同时miR-425可以恢复功能失调的多巴胺能神经院并改善PD小鼠行为缺陷[32]。这些也验证了肉苁蓉通过调控miRNA影响PD发病过程，但具体机制尚需进一步研究。
缺氧诱导因子-1（hypoxia inducible factor-1，HIF-1）是一种对氧敏感的转录激活因子，由组成型表达的β亚基和低氧诱导型的α亚基构成，参与血管生成、铁代谢、葡萄糖能量代谢和细胞增殖及存活的基因转录诱导[33]。有研究发现鱼藤酮或MPP+处理的PD模型SH-SY5Y细胞中HIF-1α的蛋白质水平降低，铁铁螯合剂去铁胺可引HIF-1α的积累，伴随着SH-SY5Y细胞中自噬的增强；当HIF-1α基因被抑制时，去铁胺诱导的自噬也相应受到抑制，表明去铁胺在鱼藤酮诱导的PD细胞凋亡中的保护作用可能是通过HIF-1α积累以及HIF-1α介导的自噬诱导[34]。食欲素A（Orexin-A）是下丘脑神经元分泌的一种神经肽，可在神经病变中发挥保护作用。有研究发现Orexin-A作用于经MPP+处理SH-SY5Y细胞后，不仅诱导HIF-1α的产生，同时激活HIF-1下游靶标如血管内皮生长因子，减轻细胞损伤，起到神经保护作用[35]。可见，肉苁蓉有效成分可能通过增加HIF-1靶基因的表达激活诱导受损细胞自噬，减轻细胞损伤。
越来越多证据支持神经炎症反应参与PD的发病机制[36]，而肿瘤坏死因子（tumor necrosis factor，TNF）通路作为重要的炎症反应通路，介导激活炎性小体，进而激活小胶质细胞，活化的小胶质细胞通过增殖和分泌炎症因子可损伤多巴胺能神经元[37]。Liang等[38]用肉苁蓉提取物松果菊苷处理亚急性PD小鼠模型后发现，松果菊苷可减少小胶质细胞和星形胶质细胞的活化，降低炎症因子TNF-α水平，上调GDNF表达，改善PD小鼠行为学症状。
通过构建中药-化合物-靶点网络图，根据degree值高低进行排序，选择前5个关键靶点，分别为MAPK8、ESR1、SRC、JAK2、AKT1。Liu等[39]发现中药提取物可以通过激活丝裂原活化蛋白激酶8（mitogen-activated protein kinase 8，MAPK8），从而促进抗凋亡蛋白Bcl-2 磷酸化，进而清除受损线粒体，间接抑制细胞凋亡，最终通过调节自噬、凋亡再平衡而发挥神经保护作用，可能也是肉苁蓉的作用机制之一。Chakrabarti等[40]发现雌激素受体（estrogen receptor alpha，ESR）激动剂可通过抑制小胶质细胞活性减轻炎症反应发挥神经保护作用；而另一项研究发现雌激素作用于ESR抑制NADPH氧化酶（reduced nicotinamide adenine dinucleotide phosphate oxidase）活性，减少神经炎症标志物，从而减轻神经损伤[41]。由此可以得出肉苁蓉有效成分可能作用于ESR1靶点减轻神经炎症反应治疗PD。与ESR相反，类固醇受体辅助活化因子（steroid receptor coactivator，SRC）作为促炎因子会造成PD小鼠多巴胺能神经元损伤[42]。丝氨酸/苏氨酸激酶1（AKT1）属于AKT家族，是参与细胞存活和凋亡调控信号转导的中心分子，磷酸化时完全激活[43]。有研究显示在PD中，维持AKT1中半胱氨酸296和310的硫醇状态对于AKT1激酶活性和防止其被蛋白酶体降解至关重要[44]。可见，上述四个分子均可能是肉苁蓉保护效应的靶点。再利用AutoDock Vina将以上5个关键靶点与肉苁蓉的6个活性成分（β-谷甾醇、花生四烯酸、脱水内酯、丁香树脂二甲醚、槲皮素和枫马京）分别进行分子对接，均成功对接，结合能均小于−5 kcal/mol，提示药物与靶点的结合活性较好。其中槲皮素和SRC具有强烈的结合活性，结合能为−9.2 kcal/mol，可从抗神经炎症、抗凋亡角度进一步深入研究肉苁蓉的作用机制。
综上，本研究通过网络药理学及分子对接探讨了临床常用于治疗PD的有效药物肉苁蓉作用靶点，发现肉苁蓉通过多组分、多通路、多靶点的协同作用发挥其在PD中的疗效，为进一步研究肉苁蓉神经保护作用机理和深入探讨PD发病机制提供了依据。另外，本研究的筛选标准与中药各组分的口服吸收度及生物利用度均有一定关系，所以并非对所有有效成分都能进行充分分析，这对我们天然药物的开发也有一定启示作用，即应结合中药成分的吸收、代谢等特点，同时考虑不同给药途径及不同个体或组织器官的分布情况进行深入研究。
-
![]()
图 1 帕金森病-肉苁蓉靶点交集韦恩图
Figure 1. Wayne diagram of target intersection of Parkinson’s disease-Cistanches Herba
![]()
图 3 GO功能富集分析
注：A：GO-BP；B：GO-CC；C：GO-MF。
Figure 3. Gene ontology function enrichment analysis
![]()
图 4 KEGG信号通路富集分析
Figure 4. Enrichment analysis of Kyoto encyclopedia of genes and geomes signal pathway
![]()
图 6 分子对接模式图
注：A：β-谷甾醇和SRC的对接图；B：花生四烯酸和ESR1的对接图；C：脱水内酯和JAK2的对接图；D：丁香树脂二甲醚和JAK2的对接图；E：槲皮素和SRC的对接图；F：枫马京和SRC的对接图。
Figure 6. Molecular docking mode
表 1 肉苁蓉有效成分
Table 1 Effective components of Cistanches Herba
编号 化合物名称 OB（%） DL BBB MOL000358 beta-sitosterol 36.91 0.75 0.99 MOL005320 arachidonate 45.57 0.20 0.58 MOL005384 suchilactone 57.52 0.56 0.28 MOL007563 Yangambin 57.53 0.81 0.01 MOL000098 quercetin 46.43 0.28 -0.77 MOL008871 Marckine 37.05 0.69 -0.10 
下载: 导出CSV
表 2 肉苁蓉有效活性成分与对应核心靶点的分子对接
Table 2 Molecular docking of effective active components of Cistanches Herba with corresponding core targets
有效活性成分 MAPK8 ESR1 SRC JAK2 AKT1 beta-sitosterol（β-谷甾醇） −7.8 −6.8 −8.1 −7.4 −6.3 arachidonate（花生四烯酸） −5.4 −6.3 −6.3 −6.1 −5.3 suchilactone（脱水内酯） −7.9 −5.9 −7.7 −8.6 −5.9 Yangambin（丁香树脂二甲醚） −7.4 −6.2 −7.5 −7.9 −6.2 quercetin（槲皮素） −8.1 −7.2 −9.2 −9 −6.6 Marckine（枫马京） −9 −7.9 −10.2 −10 −7.2 注：单位：结合能（kcal/mol）。
下载: 导出CSV
-
[1] 中华医学会神经病学分会帕金森病及运动障碍学组, 中国医师协会神经内科医师分会帕金森病及运动障碍学组. 中国帕金森病治疗指南(第四版)[J]. 中华神经科杂志,2020,53(12):973−986. [Chinese Society of Parkinson's disease and Movement Disorders, Parkinson’s Disease and Movement Disorder Section of Neurologist Branch of Chinese Medical Doctor Association. Chinese guidelines for the treatment of Parkinson's disease (fourth edition)[J]. Chin J Neurol,2020,53(12):973−986. doi: 10.3760/cma.j.cn113694-20200331-00233 [2] LI G, MA J, CUI S, et al. Parkinson's disease in China: A forty-year growing track of bedside work[J]. Translational Neurodegeneration,2019,8(1):22. doi: 10.1186/s40035-019-0162-z
[3] DENG M, ZHAO J Y, TU P F, et al. Echinacoside rescues the SHSY5Y neuronal cells from TNFα-induced apoptosis[J]. European Journal of Pharmacology,2004,505(1-3):11−18. doi: 10.1016/j.ejphar.2004.09.059
[4] 毕萃萃, 刘银路, 魏芬芬, 等. 肉苁蓉的主要化学成分及生物活性研究进展[J]. 药物评价研究,2019,42(9):1896−1900. [BI C C, LIU Y L, WEI F F, et al. Research progress on main chemical constituents and biological activities of Cistanche deserticola[J]. Drug Evaluation Research,2019,42(9):1896−1900. [5] 陆艳, 张亚杰, 阮杰, 等. 肉苁蓉颗粒剂对帕金森病大鼠模型黑质纹状体多巴胺能神经元的保护作用研究[J]. 中华中医药学刊,2016,34(12):2927−2931. [LU Y, ZHANG Y J, RUAN J, et al. Nigrostriatal dopaminergic protection of Cistanche Granule on rat model of Parkinson's disease[J]. Chinese Archives of Traditional Chinese Medicine,2016,34(12):2927−2931. [6] 尹若熙, 李刚, 俞腾飞, 等. 肉苁蓉多糖对东莨菪碱所致学习记忆障碍模型小鼠在突触可塑性方面的保护作用[J]. 中国药理学通报,2014(6):801−806,807. [YIN R X, LI G, YU T F, et al. Effect of polysaccharide from Cistanche deserticola on learning and memory deficits induced by scopolamine under improving synaptic plasticity in mice[J]. Chinese Pharmacological Bulletin,2014(6):801−806,807. [7] 周志锦. 肉苁蓉的研究进展[J]. 中医药信息,1995,12(5):28−30. [ZHOU Z J. Research progress of Cistanche deserticola[J]. Information of Traditional Chinese Medicine,1995,12(5):28−30. [8] LIU H, WANG J, ZHOU W, et al. Systems approaches and polypharmacology for drug discovery from herbal medicines: An example using licorice[J]. Journal of Ethnopharmacology,2013,146(3):773−793. doi: 10.1016/j.jep.2013.02.004
[9] GRAY M T, WOULFE J M. Striatal blood-brain barrier permeability in Parkinson’s disease[J]. J Cereb Blood Flow Metab,2015,35(5):747−750. doi: 10.1038/jcbfm.2015.32
[10] 陈文娜, 朱爱松, 丛培玮, 等. 基于网络药理学的二陈汤治疗冠心病分子机制研究[J]. 中草药,2019(2):441−448. [CHEN W N, ZHU A S, CONG P W, et al. Molecular mechanism of Erchen decoction in treatment of coronary heart disease based on network pharmacology[J]. Chinese Traditional and Herbal Drugs,2019(2):441−448. doi: 10.7501/j.issn.0253-2670.2019.02.025 [11] TROTT O, OLSON A J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading[J]. Journal of Computational Chemistry,2010,31(2):455−461.
[12] 龚普阳, 郭瑜婕, 李晓朋, 等. 基于网络药理学与分子对接技术的金花清感颗粒防治新型冠状病毒肺炎的潜在药效物质研究[J]. 中草药,2020,51(7):1685−1693. [GONG P Y, GUO Y J, LI X P, et al. Exploring active compounds of Jinhua Qinggan Granules for prevention of COVID-19 based on network pharmacology and molecular docking[J]. Chinese Traditional and Herbal Drugs,2020,51(7):1685−1693. doi: 10.7501/j.issn.0253-2670.2020.07.002 [13] 郑晓珂, 刘媛媛, 冯卫生, 等. 天然苯乙醇苷类化合物研究进展[J]. 中国新药杂志,2011,20(3):230−234. [ZHENG X K, LIU Y Y, FENG W S, et al. Development in research of natural phenylethanoid glycosides[J]. Chinese Journal of New Drugs,2011,20(3):230−234. [14] 赵卿, 张鹏, 白宇, 等. 松果菊苷通过调控GFRα1/AKT信号通路抑制MPP+诱导的多巴胺能细胞SH-SY5Y凋亡[J]. 中成药,2016(6):1225−1231. [ZHAO Q, ZHANG P, BAI Y, et al. Echinacoside inhibits apoptosis of SH-SY5Y cell through GFRα1/AKT pathway[J]. Chinese Traditional Patent Medicine,2016(6):1225−1231. [15] ZHAO Q, YANG X Y, CAI D F, et al. Echinacoside protects against MPP(+)-induced neuronal apoptosis via ROS/ATF3/CHOP pathway regulation[J]. Neuroscience Bulletin,2016,32(4):349−362. doi: 10.1007/s12264-016-0047-4
[16] ZHAO Q, GAO J, LI W, et al. Neurotrophic and neurorescue effects of echinacoside in the subacute MPTP mouse model of Parkinson's disease[J]. Brain research,2010,1346:224−236. doi: 10.1016/j.brainres.2010.05.018
[17] CHENG Y F, CHEN Y P, LI J, et al. Dietary β-sitosterol improves growth performance, meat quality, antioxidant status, and mitochondrial biogenesis of breast muscle in broilers[J]. Animals an Open Access Journal from Mdpi,2019,9(3):71.
[18] SUBRAMANIAM S R, CHESSELET M F. Mitochondrial dysfunction and oxidative stress in Parkinson's disease[J]. Progress in Neurobiology,2013,106-107(3):17−32.
[19] WANG Z J, LIANG C L, LI G M, et al. Neuroprotective effects of arachidonic acid against oxidative stress on rat hippocampal slices[J]. Chemico-Biological Interactions,2006,163(3):207−217. doi: 10.1016/j.cbi.2006.08.005
[20] 王德杰, 郑傲. 槲皮素对D-半乳糖致衰老小鼠模型学习记忆及脑内炎症的影响[J]. 中华老年多器官疾病杂志,2018,17(9):691−695. [WANG D J, ZHENG A. Effects of quercetin on learning and memory and cerebral inflammation in the D-galactose-induced aging mice[J]. Chin J Mult Organ Dis Elderly,2018,17(9):691−695. doi: 10.11915/j.issn.1671-5403.2018.09.158 [21] WANG Y H, YU H T, PU X P, et al. Myricitrin alleviates methylglyoxal-induced mitochondrial dysfunction and AGEs/RAGE/NF-κB pathway activation in SH-SY5Y cells[J]. Journal of Molecular Neuroscience:MN,2014,53(4):562−570. doi: 10.1007/s12031-013-0222-2
[22] BOYINA H K, GEETHAKHRISHNAN S L, PANUGANTI S, et al. In silico and in vivo studies on quercetin as potential anti-Parkinson agent[J]. Advances in Experimental Medicine and Biology,2020,1195:1−11.
[23] LUO Y, HOFFER A, HOFFER B, et al. Mitochondria: A therapeutic target for Parkinson’s disease?[J]. International Journal of Molecular Sciences,2015,16(9):20704−20730. doi: 10.3390/ijms160920704
[24] GANGULY G, CHAKRABARTI S, CHATTERJEE U, et al. Proteinopathy, oxidative stress and mitochondrial dysfunction: Cross talk in Alzheimer’s disease and Parkinson's disease[J]. Drug Design, Development and Therapy,2017,2017(11):797−810.
[25] TIMMONS S, COAKLEY M F, MOLONEY A M, et al. Akt signal transduction dysfunction in Parkinson's disease[J]. Neuroscience Letters,2009,467(1):30−35. doi: 10.1016/j.neulet.2009.09.055
[26] 林瑶, 杨莎莎, 刘婷, 等. 肉苁蓉对帕金森病模型大鼠PI3K、Akt、Bcl-2、Bax蛋白表达的影响[J]. 江西中医药大学学报,2019,31(1):83−86. [LIN Y, YANG S S, LIU T, et al. Effects of HerbaCistanche on the protein expression of PI3K, Akt, Bcl-2 and Bax in rats model of Parkinson's disease[J]. Journal of Jiangxi University of TCM,2019,31(1):83−86. [27] YANG P H, ZHU J X, HUANG Y D, et al. Human basic fibroblast growth factor inhibits Tau phosphorylation via the PI3K/Akt-GSK3β signaling pathway in a 6-Hydroxydopamine-induced model of Parkinson's disease[J]. Neurodegenerative Diseases,2016,16(5-6):357−369. doi: 10.1159/000445871
[28] QIAO C M, SUN M F, JIA X B, et al. Sodium butyrate causes α-synuclein degradation by an Atg5-dependent and PI3K/Akt/mTOR-related autophagy pathway[J]. Experimental Cell Research,2019,387(1):111772.
[29] PINO E, AMAMOTO R, ZHENG L, et al. FOXO3 determines the accumulation of α-synuclein and controls the fate of dopaminergic neurons in the substantia nigra[J]. Human Molecular Genetics,2014,23(6):1435−1452. doi: 10.1093/hmg/ddt530
[30] BRIGGS C E, WANG Y, KONG B, et al. Midbrain dopamine neurons in Parkinson's disease exhibit a dysregulated miRNA and target-gene network[J]. Brain Research,2015,1618:111−121. doi: 10.1016/j.brainres.2015.05.021
[31] DOXAKIS E. Post-transcriptional regulation of alpha-Synuclein expression by mir-7 and mir-153[J]. Journal of Biological Chemistry,2010,285(17):12726−12734. doi: 10.1074/jbc.M109.086827
[32] HU Y B, ZHANG Y F, WANG H, et al. MiR-425 deficiency promotes necroptosis and dopaminergic neurodegeneration in Parkinson's disease[J]. Cell Death & Disease,2019,10(8):589.
[33] SCHUMACKER P T. Hypoxia-inducible factor-1 (HIF-1)[J]. Critical Care Medicine,2005,33(5):423−425.
[34] WU Y C, LI X Q, XIE W J, et al. Neuroprotection of deferoxamine on rotenone-induced injury via accumulation of HIF-1α and induction of autophagy in SH-SY5Y cells[J]. Neurochemistry International,2010,57(3):198−205. doi: 10.1016/j.neuint.2010.05.008
[35] FENG Y, LIU T, LI X Q, et al. Neuroprotection by orexin-A via HIF-1α induction in a cellular model of Parkinson's disease[J]. Neuroscience Letters,2014,579:35−40. doi: 10.1016/j.neulet.2014.07.014
[36] ALAM Q, ALAM M Z, MUSHTAQ G, et al. Inflammatory process in Alzheimer’s and Parkinson’s diseases: Central role of cytokines[J]. Current Pharmaceutical Design,2016,22(5):541−548. doi: 10.2174/1381612822666151125000300
[37] YANG L, MAO K, YU H, et al. Neuroinflammatory responses and Parkinson’ disease: Pathogenic mechanisms and therapeutic targets[J]. Journal of Neuroimmune Pharmacology,2020,15(4):830−837. doi: 10.1007/s11481-020-09926-7
[38] LIANG Y, CHEN C, XIA B, et al. Neuroprotective effect of echinacoside in subacute mouse model of Parkinson's disease[J]. BioMed Research International,2019,2019:1−8.
[39] LIU J, LIU W, LU Y, et al. Piperlongumine restores the balance of autophagy and apoptosis by increasing BCL2 phosphorylation in rotenone-induced Parkinson disease models[J]. Autophagy,2018,14(5):845−861. doi: 10.1080/15548627.2017.1390636
[40] CHAKRABARTI M, HAQUE A, BANIK N L, et al. Estrogen receptor agonists for attenuation of neuroinflammation and neurodegeneration[J]. Brain Research Bulletin,2014,109:22−31. doi: 10.1016/j.brainresbull.2014.09.004
[41] LABANDEIRA-GARCIA J L, RODRIGUEZ-PEREZ A I, VALENZUELA R, et al. Menopause and Parkinson's disease. Interaction between estrogens and brain renin-angiotensin system in dopaminergic degeneration[J]. Frontiers in Neuroendocrinology,2016,43:44−59. doi: 10.1016/j.yfrne.2016.09.003
[42] YANG H, WANG L, ZANG C, et al. Src inhibition attenuates neuroinflammation and protects dopaminergic neurons in Parkinson’s disease models[J]. Frontiers in Neuroscience,2020,14:45. doi: 10.3389/fnins.2020.00045
[43] RAN C, WESTERLUND M, ANVRET A, et al. Genetic studies of the protein kinase AKT1 in Parkinson's disease[J]. Neuroscience Letters,2011,501(1):41−44. doi: 10.1016/j.neulet.2011.06.038
[44] AHMAD F, NIDADAVOLU P, DURGADOSS L, et al. Critical cysteines in Akt1 regulate its activity and proteasomal degradation: Implications for neurodegenerative diseases[J]. Free Radical Biology & Medicine,2014,74:118−128.
-
期刊类型引用(1)
1. 刀梅，宫树森，韦祖粉，王青芬，杨自云，吴田. 诱导剂对海巴戟莨菪亭累积的影响. 西北农林科技大学学报(自然科学版). 2024(08): 133-142 . 
百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 374
- HTML全文浏览量: 114
- PDF下载量: 55
- 被引次数: 1
