阿尔茨海默病(Alzheimer's disease, AD)是一种由多因素引起的神经退行性疾病, 其发病机制较为复杂, 包括: β-淀粉样蛋白(amyloid β-protein, Aβ)的异常沉积、tau蛋白的过度磷酸化和线粒体异常等^{[1, 2]}。据国际阿尔茨海默病协会统计, 全球AD患者约4 600万, 预计到2030年, 将达到7 470万^[3]。然而, 目前临床使用的药物仅有乙酰胆碱抑制剂和N-甲基-D-天冬氨酸受体拮抗剂, 且治疗效果都不尽如人意^[4]。因此, 对于AD这种多因素疾病, 单一靶标的干预很难得到理想的效果。治疗多因素疾病时, 研究者逐渐改变单靶点药物治疗的思路, 并向多靶点、多通路的思路转变, 中药作为天然药物的重要组成部分, 有多成分-多靶点的治疗机制, 且越来越多的天然药物通过FDA认证^[5]。因此, 阐明中药治疗多基因类疾病的物质基础及分子机制, 受到了研究者广泛的关注和重视。

远志, 始载于《神农本草经》, 性辛、温, 味苦, 归心、肾、脾经, 具有安神益智、交通心肾等功效^[6];现代药理学研究证明^[7], 远志具有改善学习记忆、抗氧化、保护神经元和保护心脑血管等作用; 临床应用与实验研究表明^[8], 远志具有明确改善阿尔茨海默病病理表现的作用。然而, 由于阿尔茨海默病发病的病因和病机比较复杂, 且远志中化学成分多样, 通过实验进行研究存在一定的困难。虽然有实验研究表明, 远志中化学成分具有保护脑活性和抗AD的作用, 但是其治疗AD的主要活性成分和分子机制尚未完全阐明。因此, 本研究运用系统药理学的思路和方法, 从分子角度及系统水平分析远志防治阿尔茨海默病的药效活性成分和可能的分子网络机制, 挖掘远志的潜在治疗作用, 以期为临床使用提供科学依据。

小分子数据库的构建  收集TCMID数据库^[9] (http://www.megabionet.org/tcmid/)中远志化学成分结合已报道的远志化学成分相关文献^{[8, 10-13]}, 共得到远志中的化学成分111个。由ChemSpider^[14] (http://www.chemspider.com/Default.aspx)、PubChem (https://pubchem.ncbi.nlm.nih.gov/)数据库中下载上述化学成分的2D分子结构文件, 截止数据分析时, 暂未被数据库收录, 但已有文献报道的成分通过Chemdraw软件进行2D分子结构的绘制。将上述所有2D分子结构文件转化为3D结构, 构建小分子数据库。

ADME评估  ADME/T (absorption, distribution, metabolism and excretion of drugs/toxicity)是由王永华团队^[15]所开发的一种借助机器学习、计算化学及现代统计学等方法来预测药物先导分子与机体相互作用的中药系统药理学方法模型。将上述搜集整理后的化学成分通过ADME/T模型, 主要参考口服生物利用度(bioavailability, OB)和类药性(drug-likeness, DL)两个关键参数进行活性化合物的筛选。继而使用Discovery Studio (DS)软件中ADME模块^[16]对筛选后的化学成分进行血脑屏障通透性(blood brain barrier, BBB)考察, 并以打分所得血脑屏障通透性水平(BBB level)为评判标准。

靶点预测及验证  Pharmmapper服务器(http://59.78.96.61/pharmmapper)是一个免费访问的Web服务器, 它运用反向药效团匹配的方法通过对TargetBank、Binding DB、DrugBank与PDTD数据库的检索, 自动为给定小分子(药品、天然产物或其他新发现的化合物与未识别的连接靶点)找到最佳的映射靶点^{[17, 18]}。将通过OB、DL和BBB筛选出的小分子, 以MOL2的格式上传到Pharmmapper在线分析平台, 其中Select Targets Set选择Human Protein Targets Only (v2010, 2241), 其余参数均为默认设置。本研究将上述预测到的靶点信息与TTD (Therapeutic Target Database, http://bidd.nus.edu.sg/group/cjttd/)、CTD (Comparative Toxicogenomics Database, http://ctdbase.org/)数据库中阿尔茨海默病的靶点进行验证, 筛选出有明确抗阿尔茨海默病的靶点。

构建成分-靶点-疾病网络模型  根据上述远志化学成分-潜在靶点预测结果, 应用Cytoscape软件分别构建靶点-疾病(target-disease, T-D)与化学成分-靶点网络(component-target, C-T)模型, 并运用软件中Merge功能, 将两个网络模型合并成化学成分-靶点-疾病(component-target-disease, C-T-D)网络模型。运用软件中Network analyer功能, 进一步对网络模型中的相互关系进行分析。

通路分析及靶点功能注释  将筛选的靶点, 进一步映射到KEGG数据库(http://www.kegg.jp/)中, 并且将对应的通路进行整理与分析, 运用Visio软件进行通路图的绘制。同时通过MAS3.0 (http://bioinfo.capitalbio.com/mas3/)生物分子功能注释平台对上述靶点进行GO富集分析, 将Species选择Homo sapiens (Human); Molecule Type选择protein; Database Symbols选择UniProtKB; Functions选择gene、mRNA、protein和pathway, 其余参数为默认设置, 并且将富集的疾病部分进行MeSH (https://www.nlm.nih.gov/mesh/)分类分析。

本研究将收集的111个远志化学成分, 分别用C1~C111表示(表 1)。通过ADME评估, 以OB≥ 30%、DL≥0.18^[19]和BBB level≥3^[20]为标准, 对111个化学成分进行筛选, 得到符合条件的成分共10个, 分别为:皂苷类化合物1个、远志口山酮类化合物3个、寡糖脂类化合物3个和其他类3个。将上述10个化学成分, 通过Pharmmapper服务器与TTD、CTD数据库验证后, 预测出关于AD的作用靶点共13个: β-分泌酶(beta-secretase 1, BACE1)、钙调蛋白(calmodulin, CALM1)、凝血素(prothrombin, F2)、糖原合酶激酶-3β (glycogen synthase kinase-3 beta, GSK3β)、胰岛素样生长因子(insulin-like growth factor Ⅰ, IGF1)、胰岛素样生长因子1受体(insulin-like growth factor 1 receptor, IGF1R)、尿激酶型纤溶酶原激活物(urokinase-type plasminogen activator, PLAU)、过氧化物酶体增殖物激活受体γ (peroxisome proliferator-activated receptor gamma, PPARγ)、超氧化物歧化酶(superoxide dismutase, SOD2)、胆碱酯酶(cholinesterase, BCHE)、一氧化氮合酶(nitric-oxide synthase, NOS3)、雌激素受体(estrogen receptor, ESR)、胱氨酸天冬氨酸蛋白酶(caspase-3, CASP3) (表 2)。将化学成分与靶点对应关系进行分析后发现, 化合物1, 3, 6, 7-tetramethoxyxanthone的口服生物利用度最高, 为102.4%, 显示此化合物口服时, 通过胃肠道中的被动吸收, 进入血液循环后的疗效可能得到较好地发挥。其次为化合物3-hydroxy-2, 8-dimeth oxyxanthone (99.07%)和tenuifoliose D (92.44%); 就类药性而言, 化合物onjisaponin E (0.7186) 则位于首位, 表明此化合物开发为药物的可能性最高; 血脑屏障通透性方面, 化合物tenuifoliose B (OB = 76.99%, DL = 0.3079, BBB level = 1) 具有较好的生物利用度与类药性的同时, 其血脑屏障通透性最高, 且实验研究表明, 其可与大鼠齿状回突触的相互作用, 能改善大鼠学习记忆的活性^[7]。分析结果显示:在10个活性成分中, 10.00%的成分对应4个靶点, 40.00%的成分对应5个靶点, 30.00%的成分对应6个靶点, 20.00%的成分对应7个靶点。这些成分与靶点的对应数量关系验证了中药通过多成分-多靶点发挥作用的机制, 同时验证了中药中化学成分间的协同作用, 显示中药治疗多因素多基因类疾病的优势。

表1(Table 1)

Table 1 Ingredients of Polygala Radix ID Ingredient C1 Onjisaponin A C2 Onjisaponin B C3 Onjisaponin E C4 Onjisaponin F C5 Onjisaponin G C6 Onjisaponin J C7 Onjisaponin L C8 Onjisaponin O C9 Onjisaponin R C10 Onjisaponin S C11 Onjisaponin T C12 Onjisaponin V C13 Onjisaponin W C14 Onjisaponin X C15 Onjisaponin Y C16 Onjisaponin Z C17 Onjisaponin Vg C18 Onjisaponin Pg C19 Onjisaponin Gg C20 Onjisaponin Fg C21 Onjisaponin Qg C22 Onjisaponin Ng C23 Onjisaponin Sg C24 Onjisaponin Ug C25 Onjisaponin Tg C26 Onjisaponin Wg C27 E-Onjisaponin H C28 Z-Ponjisaponin H C29 Sibiricasaponin A C30 Sibiricasaponin B C32 Sibiricasaponin E C33 Tenuifolisaponin A C34 Tenuifolisaponin B C35 Polygalasaponin XXVIII C36 Polygalasaponin XXXII C37 Lancerin C38 Onjixanthone Ⅰ C39 Onjixanthone Ⅱ C40 Polygalaxanthone Ⅲ C41 Polygalaxanthone Ⅳ C42 Polygalaxanthone Ⅴ C43 Polygalaxanthone Ⅵ C44 Polygalaxanthone Ⅶ C45 Sibiriphenone A C46 Sibiricaxanthone A C47 Sibiricaxanthone B C48 1, 2, 3-Trimethoxy-7-hydroxyxanthone C49 1, 3, 6-Trihydroxy-2, 7-dimethoxyxanthone C50 1, 2, 7-Trimethoxy-3-hydroxyxanthone C51 1, 2, 3, 7-Tetramethoxyxanthone C52 1, 7-Dihydroxy-3-methoxyxanthone C53 1, 7-Dihydroxy-2, 3-dimethoxyxanthone C54 1, 7-Dihydroxy-2, 3-methlendioxyxanthone C55 1, 7-Dihydroxy-2, 3-di-methylene-dioxyxanthone C56 6-Hydroxy-2, 3, 6, 7-tetramethoxyxanthone C57 6-Hydroxy-1, 2, 3, 7-tetramethoxyxanthone C58 1, 7-Dihydroxyxanthone C59 1, 7-Dimethoxyxanthone C60 1-Hydroxy-3, 7-dimethoxyxanthone C61 1-Hydroxy-3, 6, 7-trimethoxyxanthone C62 6, 8-Dihydroxy-1, 2, 4-trimethyoxyxanthone C63 6, 8-Dihydroxy-1, 2, 3-trimethyoxyxanthone C64 1, 2, 3, 6, 7-Pentamethoxyxanthone C65 1, 3, 7-Trihydroxyxanthone C66 1, 6, 7-Trihydroxy-2, 3-dimethoxyxanthone C67 3-Hydroxy-2, 8-dimethoxyxanthone C68 7-Hydroxy-1, 2, 3-trimethoxyxanthone C69 3-Hydroxy-1, 2, 7-trimethoxyxanthone C70 2, 3, 8-Trimethoxyxanthone C71 1, 3, 6, 7-Tetramethoxyxanthone C72 1, 3, 7-Trimethoxyxanthone C73 7-O-Methylmangiferin C74 Tenuifoliose A C75 Tenuifoliose B C76 Tenuifoliose C C77 Tenuifoliose D C78 Tenuifoliose E C79 Tenuifoliose F C80 Tenuifoliose G C81 Tenuifoliose H C82 Tenuifoliose I C83 Tenuifoliose J C84 Tenuifoliose K C85 Tenuifoliose L C86 Tenuifoliose M C87 Tenuifoliose N C88 Tenuifoliose O C89 Tenuifoliose P C90 Tenuifoliose Q C91 Tenuifoliside A C92 Tenuifoliside B C93 Tenuifoliside C C94 Tenuifoliside D C95 Tenuifoliside E C96 Sibiricose A1 C97 Sibiricose A2 C98 Sibiricose A3 C99 Sibiricose A4 C100 Sibiricose AS C101 Sibiricose A6 C102 Sibiricose A7 C103 3, 6'-Disinapoylsucrose C104 Polygalatenosides B C105 3, 4-Dimethoxy cinnamylic acid C106 Polygalitol C107 Polygalacic acid C108 3, 4, 5-Trimethoxy cinnamylic acid C109 Daucosterol C110 Sibiricaphenone C111 2-Methoxyhydroquinone

Table 1 Ingredients of Polygala Radix

表2(Table 2)

Table 2 Information of potential targets from active ingredients of Polygala Radix. OB: Bioavailability; DL: Drug-likeness; LogBB: Blood brain partitioning; BBB level: Blood brain barrier level; BACE1: Beta-secretase 1; CALM1: Calmodulin; F2: Prothrombin; GSK3β: Glycogen synthase kinase-3 beta; IGF1: Insulin-like growth factor 1; IGF1R: Insulin-like growth factor 1 receptor; PLAU: Urokinase-type plasminogen activator; PPARγ: Peroxisome proliferator-activated receptor gamma; SOD2: Superoxide dismutase; BCHE: Cholinesterase; NOS3: Nitric-oxide synthase; ESR: Endothelial; CASP3: Caspase-3 ID Ingredient OB/% DL LogBB BBB level Target C3 Onjisaponin E 40.41 0.718 6 -0.454 2 GSK3β, BACE1, CALM1, F2, IGF1 C59 1, 7-Dimethoxyxanthone 56.87 0.333 6 -0.513 2 BACE1, BCHE, GASP3, ESR, F2, GSK3β C67 3-Hydroxy-2, 8-dimethoxyxanthone 99.07 0.209 -0.496 2 BACE1, ESR, F2, SOD2, GSK3β, NOS3 C71 1, 3, 6, 7-Tetramethoxyxanthone 102.44 0.336 2 -0.616 3 BCHE, ESR, BACE1, GSK3β, F2, PLAU, PPARγ C75 Tenuifoliose B 76.99 0.307 9 0.054 1 BACE1, GASP3, F2, GSK3β, PLAU C77 Tenuifoliose D 92.44 0.299 9 -0.496 2 GSK3β, IGF1, BACE1, F2, IGF1R, SOD2 C101 Sibiricose A6 31.38 0.516 -0.700 3 PLAU, BACE1, GSK3β, GASP3, F2 C105 3, 4-Dimethoxycinnamic acid 39.58 0.399 4 -0.700 3 GSK3β, NOS3, BACE1, F2 C106 Polygalitol 35.47 0.249 2 -0.35 2 BACE1, F2, GSK3β, IGF1R, PLAU C107 Polygalacic acid 31.43 0.698 6 -0.238 2 GSK3β, GASP3, IGF1, BACE1, ESR, F2, PPARγ

Table 2 Information of potential targets from active ingredients of Polygala Radix. OB: Bioavailability; DL: Drug-likeness; LogBB: Blood brain partitioning; BBB level: Blood brain barrier level; BACE1: Beta-secretase 1; CALM1: Calmodulin; F2: Prothrombin; GSK3β: Glycogen synthase kinase-3 beta; IGF1: Insulin-like growth factor 1; IGF1R: Insulin-like growth factor 1 receptor; PLAU: Urokinase-type plasminogen activator; PPARγ: Peroxisome proliferator-activated receptor gamma; SOD2: Superoxide dismutase; BCHE: Cholinesterase; NOS3: Nitric-oxide synthase; ESR: Endothelial; CASP3: Caspase-3

运用Cytoscape软件构建的C-T-D网络模型中(图 1), 共包含10个化学成分, 13个靶点及1种疾病, 其中中心透明的大脑图形代表AD疾病, 粉红色三角形代表靶蛋白, 蓝色圆形代表远志中的化学成分。利用Network analyer功能对网络模型进行分析发现, 在C-T-D网络模型中包含24个节点, 74个边, 其网络半径为2, 网络直径为3, 最短路径为552, 网络的特征路径长度为1.982, 表明经过1.982步即可互相转化。从C-T-D图可以发现, 远志中10个化学成分与多个靶点间产生相互作用, 此现象符合中药多成分-多靶点的治疗机制。从化合物的网络特征分析来看, BACE1与GSK3β为C-T-D网络模型中网络度最好(核心)的靶点, 10个活性成分均与这两个靶点有潜在的结合作用, 其次分别为F2、ESR、CASP3、PPARγ、PLAU、BCHE、NOS3、IGF1、SOD2、CALM1和IGF1R。

Figure 1

Figure 1

Figure 1 omponents-targets-disease network. () Main active ingredients; () Targets; () Disease. The size of the graphic area indicates the importance of the nodes in the network

经过文献挖掘并对上述数据研究发现, 远志可能通过多种途径治疗AD:其一, 研究表明^{[21, 22]}, BACE1在Aβ生成中起着关键的作用, 而PPARγ^[23]的活化能够降低Aβ的沉积、增强海马区认知功能和逆转记忆力下降, 因此, 远志中的化学成分可能通过对这两个蛋白的调控而发挥治疗AD效果; 其二, 靶点GSK3β^[24]负责Tau蛋白的过度磷酸化, 在Tau蛋白过度磷酸化的假说中, 被认为是其中的关键分子, 并且GSK3β与细胞凋亡相关, 远志中19个活性成分均与此靶点存在相互作用, 因此推测远志可能通过调控Tau蛋白的磷酸化, 达到治疗AD的效果; 第三, 靶点NOS3与靶点SOD2均与氧化作用相关^{[25, 26]}, 而过度的氧化标志物在AD患者脑中的改变非常明显, 该结果预示远志可能通过抗氧化作用, 达到治疗AD的作用; 最后, 靶点ESR、IGF和IGF1R均与细胞的增殖分化、生长存活相关^{[27, 28]}, 这3个靶点在治疗AD的过程中可能起到协同或辅助的治疗作用。综上所述, 远志治疗AD可能通过对Aβ异常沉积与Tau蛋白过度磷酸化进行调控、抗氧化, 参与细胞增殖分化等多个途径发挥其治疗效果, 这可能与远志改善学习记忆、抗氧化和神经元保护等药理作用相关。

将上述13个靶点, 全部映射到KEGG数据库中, 并对结果进行分析整理(图 2), 图中橘色长方形为靶点所映射的位点, 灰色长方形为未映射但与通路有关的位点。通过对通路图的分析发现:首先, Tau蛋白磷酸化水平受到细胞内CASP3和GSK3β的共同调节, CASP3和GSK3β表达平衡对Tau蛋白是否能够磷酸化具有重要的意义。结合表 2发现, 远志中这10个活性成分均与GSK3β存在潜在的结合关系, 5个活性成分(onjisaponin E、1, 7-dimethoxyxanthone、tenuifoliose B、sibiricose A6和polygalacic acid)与CASP3存在相关关系, 初步推测, 上述成分可能通过3条途径治疗AD:其一, 10个活性成分与靶点GSK3β作用, 控制Tau蛋白的过度磷酸化; 其二, 与CASP3相互作用, 调节Tau蛋白的正常磷酸化。这两种途径共同发挥作用, 使Tau蛋白的磷酸化水平处在正常的生理状态, 从而控制神经纤维缠结(NFTs), 达到直接治疗AD的作用; 第三, 化合物onjisaponin E与靶点CALM1相互作用, 间接参与线粒体介导的细胞凋亡。通过对通路图分析发现, 化合物onjisaponin E可能通过控制线粒体过多分裂, 从而调节细胞损伤和神经元的死亡, 达到干预AD发病的作用。上述两条直接与AD相关的作用途径均与细胞凋亡通路相关。除细胞凋亡通路外, 还发现远志中的活性成分可能参与其他通路的调节, 进而对AD起到辅助治疗的作用:其一, 化合物onjisaponin E、1, 3, 6, 7-tetramethoxy xanthone、polygalacic acid与IGF1结合后, 通过作用IGF1R可以相继激活胰岛素信号通路(insulin signaling pathway)、PI3K-Akt-mTORC signaling pathway, 最终激活PPARγ。动物实验证明^[23], PPARγ的活化能够降低Aβ的沉积、增强海马区认知功能、逆转记忆力下降; 其二, 远志中的活性成分3-hydroxy-2, 8-dimeth oxyxanthone可直接与NOS3相互作用来调控心肌细胞的存活或使血管平滑肌松弛。综上所述, 远志中的化学成分通过对多靶点、多通路的调节, 从而达到干预或治疗AD的作用。

Figure 2

Figure 2

Figure 2 Enriched KEGG pathways of potential targets from main active ingredients of Polygala Radix. The orange rectangle () is the site where the target is mapped, the gray rectangle () is the unmapped but channel-related locus

运用MAS 3.0生物分子注释系统对远志化学成分预测出的13个作用靶点进行GO富集分析, 将其中富集的疾病部分进行MeSH分类。结果如图 3所示, 发现Alzheimer's disease所对应的靶点数目最多, 为10个; 除AD外, 有6个靶点对应了如帕金森病(Parkinson disease)、亨廷顿病(Huntington disease)、肌萎缩性侧索硬化(amyotrophic lateral sclerosis)和肌强直性营养不良(myotonic dystrophy)等其他神经退行性疾病; 所占比例最大(70.37%)的疾病分类为神经系统疾病(nervous system diseases, MeSH分类为C10)。结果表明, 远志除具有治疗AD的作用, 可能还有治疗其他神经退行性疾病乃至神经系统类疾病的潜在作用, 这与远志安神益智、解郁的功效, 促进神经干细胞增殖、抗氧化、清除自由基和保护神经胆碱能系统等药理作用有关。

Figure 3

Figure 3

Figure 3 GO biological pathway analysis and MeSH classification results

本研究通过ADMT评估、反向药效团匹配、构建C-T-D网络模型、KRGG分析、GO富集分析及MeSH分类, 明确了远志中10个活性分子, 其中包含口山酮类化合物3个、皂苷类化合物1个、寡糖酯类化合物3个和其他类3个。通过对网络模型的分析发现, BACE1与GSK3β为C-T-D网络模型中网络度最好(核心)的靶点, 而远志可能作用于BACE1、PPARγ靶点调控Aβ, 作用于靶点GSK3β调控Tau蛋白的磷酸化, 作用于靶点NOS3与SOD2调节氧化作用, 作用于靶点ESR、IGF和IGF1R调节细胞的增殖分化和生长存活等多途径防治AD。通过多途径治疗AD, 充分发挥了中药治疗疾病的多成分-多靶点机制, 但成分间存在的是协同作用还是拮抗作用暂不明确, 需要进一步的实验验证。通过KEGG分析, 发现远志治疗阿尔茨海默病主要通过调控细胞凋亡通路直接发挥疗效, 而调节细胞凋亡通路的途径可能包含两种: ① 活性成分与靶点GSK3β作用, 控制Tau蛋白的过度磷酸化; 与CASP3相互作用, 调节Tau蛋白的正常磷酸化。这两种途径共同发挥作用, 使Tau蛋白的磷酸化水平处在正常的生理状态, 从而控制神经纤维缠结, 达到直接治疗AD的作用; ② 化合物onjisaponin E可能通过控制线粒体过多分裂, 从而调节细胞损伤和神经元的死亡, 达到治疗AD的作用。KEGG结果显示, 远志中化学成分可能通过激活PPARγ, 调控心肌细胞的存活或使血管平滑肌松弛达到辅助的治疗效果, 但远志中的成分对通路是上调还是下调的具体机制还不明确, 需要进一步阐述。GO分析及MeSH分类的结果显示, 远志对神经系统类疾病存在着潜在的治疗作用。

本研究采用系统药理学的方法对远志中化学成分进行综合性分析, 阐释了远志中潜在活性成分群的分子作用机制, 进一步揭示了中药多成分-多靶点-多通路的治疗机制, 为中药防治由多基因引起的复杂疾病的治疗机制研究提供了思路。