药学学报  2020, Vol. 55 Issue (2): 256-264     DOI: 10.16438/j.0513-4870.2019-0521   PDF    
基于系统的化合物-靶点相互作用预测模型的消栓通络方抗脑卒中网络药理学研究
郑一夫1,2, 孔令雷1, 贾皓1, 张宝月1, 王喆1, 许律捷1, 刘艾林1, 杜冠华1     
1. 中国医学科学院、北京协和医学院药物研究所, 北京100050;
2. 武汉大学药学院, 湖北 武汉 430072
摘要: 消栓通络方临床上主治脑血栓引起的精神呆滞、言语迟涩等症状,疗效显著,但其作用机制尚不明确。本文通过搜集消栓通络方中的化学成分和治疗脑卒中相关靶点,获得1251个化学成分和10个脑卒中相关靶点,采用朴素贝叶斯和递归分割等机器学习算法,基于分子指纹和分子描述符,结合分子对接方法,构建了脑卒中10个相关靶点的18个化合物-靶点相互作用预测模型。应用这些模型预测了消栓通络方的活性化学成分及其作用靶点,发现了153个潜在活性化学成分,其中22个可以与2个以上脑卒中药物靶点相互作用。在此基础上,利用网络构建专业软件,构建了化学成分-靶点网络,并通过Gene-Ontology(GO)富集分析,确证了靶点重要的生物过程,如凝血(blood coagulation)、血管生成正调控(positive regulation ofangiogenesis)和离子转运正调控(positive regulation of ion transport)等。在此基础上,构建了靶点-通路网络。本研究利用机器学习、分子对接、虚拟筛选、数据挖掘及网络构建等技术方法,探索并部分揭示了消栓通络方抗脑卒中的活性物质基础及其化学成分-靶点-通路的网络作用机制,为消栓通络方的深入研究提供重要信息依据。
关键词: 消栓通络方    脑卒中    药物靶点    机器学习    分子对接    虚拟筛选    网络药理学    
Network pharmacology study on anti-stroke of Xiaoshuan Tongluo formula based on systematic compound-target interaction prediction models
ZHENG Yi-fu1,2, KONG Ling-lei1, JIA Hao1, ZHANG Bao-yue1, WANG Zhe1, XU Lü-jie1, LIU Ai-liu1, DU Guan-hua1     
1. Institute of Materia Medica, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing 100050, China;
2. College of Pharmacy, Wuhan University, Wuhan 430072, China
Abstract: Xiaoshuan Tongluo formula is effective in treating mental retardation and speech astringency caused by cerebral thrombosis, but its mechanism remains unclear. In this investigation, by collecting the chemical constituents from Xiaoshuan Tongluo formula and the targets related to stroke, we obtained 1 251 constituents from the formula and 10 drug targets related with stroke. We established 18 prediction models of compound-target interaction for 10 stroke-related targets, using molecular docking method and machine learning methods includes Naive Bayesian and recursive partitioning based on the input of molecular fingerprints and molecular descriptors. Using these models, we predicted the active chemical constituents from Xiaoshuan Tongluo formula and their drug targets, 153 potential active constituents were discovered, 22 among them could interact with at least two drug targets related with stroke. On this basis, the chemical constituent-target network was constructed using network construction software, and then the important metabolic pathways of the targets were identified by using Gene-Ontology (GO) enrichment analysis, such as blood coagulation, positive regulation of angiogenesis, positive regulation of ion transport and so on. On this basis, a target-pathway network was constructed. In conclusion, using machine learning, molecular docking, virtual screening, data mining and network construction, this study explored and partially revealed the active chemical constituents and chemical constituent-target-pathway network action mechanism of Xiaoshuan Tongluo formula against stroke, which will provide important information for its further study.
Key words: Xiaoshuan Tongluo formula    stroke    drug target    machine learning    molecular docking    virtual screening    network pharmacology    

脑卒中是由于脑部血管突然破裂或因血管阻塞导致血液不能流入大脑而引起脑组织损伤的一组疾病。据国家卫生部第三次全国死亡原因调查显示, 中国脑卒中年发病率超过200万, 每年150万人因脑卒中死亡, 死亡率已经是西方发达国家的4~5倍, 是日本的3.5倍[1], 跃升为我国第一大致死原因。与此相应的是, 获准脑卒中治疗的药物寥寥无几, 目前临床上仅存在唯一单靶点作用药物组织型纤溶酶原激活剂(tissue plasminogen activator, tPA)[2]治疗脑卒中。但由于其治疗窗狭窄, 对患者发病情况有严苛条件要求, 仅3%~5%患者可以接受治疗[3], 且有导致脑出血的风险[4]。因此, 有关脑卒中的治疗亟待提高。脑卒中分为缺血性脑卒中和出血性脑卒中, 其中缺血性脑卒中占70%左右[5, 6]。因此, 本文将针对缺血性脑卒中开展相关研究。

中药消栓通络方为临床常用的经验方, 自1990年起为历版中国药典所收录, 主要用于脑卒中恢复期的治疗。消栓通络方包含川芎、丹参、黄芪、泽泻、三七、槐花、桂枝、郁金、木香、冰片、山楂共11味中药配伍而成[6], 具有活血化瘀, 温经通络之功效, 主治血脂增高, 脑血栓引起的精神呆滞、舌质发硬、言语迟涩、发音不清、手足发凉和活动疼痛。但此常用药方鲜有研究, 且其发挥功效的活性成分及作用机制尚未见深入的研究报道。

中药复方使用多味草药配伍而成, 作用成分复杂, 可以通过与脑卒中相关的多个靶点产生药效作用。通过生物学实验手段来确定中药中有效成分的抗脑卒中机制不仅消耗大量的人力、物力和财力, 且任务艰巨[7]。因此, 急需发展新的、有效的方法来研究中药抗脑卒中的网络作用机制。随着中药现代化发展, 国内外建立了多个中药成分数据库, 如TCMID、TCMSP、HIT、TCMDB@Taiwan等[8]。随着生命科学的发展和蛋白质结构解析技术的进步, 越来越多的蛋白结构得到解析, 在此基础上构建了许多重要的药物靶点数据库, 如Thomson Reuters Integrity[9]、Therapeutic Target Database (TTD)等[10]。基于现有数据库和机器学习方法, 进行数据挖掘及相关生物信息学分析, 结合计算模拟方法和网络构建技术, 使中药活性成分发现和网络机制研究成为可能。数据库的快速发展推动了国内外大量对于中药复方有效成分机制的研究[11, 12], 本文采用类似的方法进行消栓通络方治疗脑卒中网络药理学机制研究。

本文通过中药成分数据库及上市药物靶点数据库, 结合文献分析, 搜集消栓通络方中的化学成分及治疗脑卒中相关靶点信息, 基于分子指纹和分子描述符, 采用朴素贝叶斯和递归分割等机器学习算法, 结合分子对接计算, 系统地构建抗脑卒中的化合物-靶点相互作用预测模型。应用预测模型, 预测消栓通络方潜在活性成分及其作用靶点。在此基础上, 利用生物网络构建软件, 构建化学成分-靶点网络和靶点-通路网络, 以揭示消栓通络方抗脑卒中的活性物质基础及其网络作用机制, 为消栓通络方的深入研究提供信息依据。

材料与方法

消栓通络方中药成分的搜集  利用中药药理学数据库TCMSP (lsp.nwu.edu.cn/tcmsp.php)、TCMID (183.129.215.33/tcmid/)和Batman-TCM (bionet.ncpsb.org/batman-tcm/), 分别以消栓通络方中11味药材名称为关键词, 检索并保存化学成分结构等信息。

脑卒中靶点及靶点预测模型的选择  在Thomson Reuters Integrity和Therapeutic Target Database (TTD)数据库中, 针对stroke搜索临床试验3期及上市药物, 进一步检索药物作用靶点, 经过整理, 结合文献分析, 确定脑卒中相关靶点。然后进入BindingDB数据库搜索靶点的已知活性化合物, 由于机器学习模型都设置了9个描述符, 当活性化合物的个数是描述符个数的3倍以上, 建立的模型才具有统计学意义, 因此, 以至少50个活性化合物作为基于活性化合物结构建立预测模型的标准, 以保证机器学习模型可靠性。若靶点的活性化合物个数大于50, 则采用机器学习方法建模; 若活性化合物个数不足50, 则采用分子对接方法建模。

化合物-靶点相互作用机器学习模型的建立  对于活性化合物数大于50个的靶点, 并以IC50 < 10 000 nmol·L-1为条件筛选活性配体。从中随机挑选25%的配体, 上传到DUD-E [a database of useful (docking) decoys — enhanced, http://dude.docking.org/]数据库, 以产生非活性配体。

在此基础上, 按照活性化合物与非活性化合物1:3的比例, 测试集与训练集1:3的比例随机分配所有化合物。然后采用朴素贝叶斯[13]和递归分隔[14]算法, 基于ECFP_6[15]和MACCS[16]分子指纹, 以及分子描述符, 在Discovery Studio 2018平台上构建模型, 建模流程如图 1所示, 其中MACCS分子指纹利用PaDEL Descriptor软件进行计算。针对每个脑卒中相关靶点, 采用不同算法和分子指纹, 分别构建了4个相关模型, 即每个靶点有4个预测模型NB(ECFP_6)、NB(MACCS)、RP(ECFP_6)和RP(MACCS)。经过前期分析发现, 当采用朴素贝叶斯算法时, 基于ECFP_6分子指纹的模型表现更好; 当使用MACCS分子指纹时, 递归分隔算法表现更好[17]。因此, 最后每个靶点采用NB(ECFP_6)和RP(MACCS)的组合预测模型。

图 1 Workflow of the construction of machine learning models

构建的模型需经过训练集的五倍交叉验证[18]与测试集验证, 以评价预测模型预测能力。验证指标基于真阳性(true positive, TP)、真阴性(true negative, TN)、假阳性(false positive, FP)和假阴性(false negative, FN)的基础数据进行计算, 包括敏感度(sensitivity, SE)、精确度(specificity, SP)、总体预测精度(overall prediction accuracy, Q)、马修斯相关系数(Matthews correlation coefficient, MCC)和受试集工作特征曲线的曲线下面积(area under the receiver operating characteristic curve, AUC)。各指标计算公式如下:

$ {\rm{SE = }}\frac{{{\rm{TP}}}}{{{\rm{TP + FN}}}} $ (1)
$ {\rm{SP = }}\frac{{{\rm{TN}}}}{{{\rm{TN + FP}}}} $ (2)
$ Q{\rm{ = }}\frac{{{\rm{TP + TN}}}}{{{\rm{TP + TN + FP + FN}}}} $ (3)
$ {\rm{MCC = }}\frac{{{\rm{TP \times TN - FP \times FN}}}}{{\sqrt {\left( {{\rm{TP + FN}}} \right)\left( {{\rm{TP + FP}}} \right)\left( {{\rm{TN + FN}}} \right)\left( {{\rm{TN + FP}}} \right)} }} $ (4)

另外, AUC为ROC曲线下与坐标轴围成的面积, 最大不会超过1, 数值越大表示预测效果越好。

分子对接模型的建立  对于活性化合物小于50个的靶点, 在PDB (Protein Data Bank)数据库下载蛋白晶体结构, 选择分辨率 < 2.5 Å并具有配体的复合物, 以保证对接的精度。在对接前对蛋白分子进行预处理, 去除结构周围的水分子, 加上氢原子, 提取所有配体, 进行逐一重新对接, 其中对接的活性口袋由原始配体位置决定。使用Discovery Studio 2018中的CDocker分子对接方法, 对蛋白与蛋白配体进行对接[19]。对接结果通过RMSD[18] (root-mean-square deviation)进行验证, RMSD值小于2.5则认为分子对接结果有效。从所有配体中选择RMSD值最小的配体对应的活性口袋进行中药化学成分的对接。

中药化学成分-靶点相互作用预测  基于上述建立的所有机器学习模型和分子对接模型对消栓通络方中的化学成分进行靶点预测。对于机器学习模型, 预测结果为True或False, True表示化合物与靶点可以相互作用, False则相反。对于分子对接模型, 根据-CDocker energy计算预测结果, -CDocker energy可以说明配体与靶点之间结合强度, 对于-CDocker energy值比原配体大的中药化学成分, 认为其对接蛋白为化学成分的潜在靶点。然后, 对预测结果进行整合分析, 统计消栓通络方中活性化学成分与各个脑卒中相关靶点的作用关系。

作用网络分析  统计潜在活性化学成分及其作用靶点, 利用网络构建软件Cytoscape[20]构建化学成分-靶点作用网络。之后, 在DAVID数据库中, 对脑卒中靶点进行Gene-Ontology (GO)[21]富集分析, 选择P < 0.05的靶点代谢通路, 并构建靶点-通路网络[22, 23]

结果 1 消栓通络方中药成分的搜集

通过中药药理学数据库TCMSP、TCMID和Batman-TCM, 分别以药材名称为关键字, 检索并保存化合物。搜集整理得到各单味药的化学成分数目如下:冰片1个、川芎189个、丹参202个、桂枝220个、槐花27个、黄芪87个、木香106个、三七119个、山楂32个、郁金222个和泽泻46个, 总计1 251个化学成分。

2 脑卒中靶点的搜集与分类

利用Thomson Reuters Integrity和TTD数据库, 并结合文献分析, 共确定10个脑卒中相关靶点, 分别是F3、F10、F11、MMP9、P2RY1、TBXA2R、VEGFA、SERPINE1、THBD和P2RY12。其中包含凝血系统相关靶点4个, 血小板相关靶点3个、保护血管内皮靶点2个及纤溶系统靶点1个, 具体名称、UniProt序列号及其功能如表 1所示。其中THBD和P2RY12, 由于BindingDB数据库中缺少其配体相关信息, 采用分子对接方法建模, 其余8个靶点采用机器学习方法建模。

表 1 The name and classification of stroke related targets
3 机器学习模型

从BindingDB数据库和DUD-E数据库分别得到靶点活性配体及非活性配体后, 按照活性化合物与非活性化合物1:3比例, 测试集与训练集1:3的比例随机分配所有化合物, 各个模型中化合物分配情况如表 2所示。其中Tanimoto index为评价数据集中化合物空间多样性的指数, 化合物空间多样性越大, 其数值越小, 构建模型预测能力则越好。由表中数据可知, 所有靶点的数据集的Tanimoto index[24]均小于0.150, 说明活性化合物与非活性化合物相似度较低, 数据集中化合物空间多样性较高。

表 2 Diversity analysis of compounds datasets

利用Discovery Studio 2018平台, 基于ECFP_6和MACCS分子指纹, 以及分子描述符, 采用朴素贝叶斯和递归分隔算法构建模型。经过训练集的五倍交叉验证及测试集验证, 评价模型的有效性和精确度, 其结果如表 3~6所示, 其中马修斯相关系数MCC和受试集工作特征曲线的曲线下面积AUC作为评价模型预测能力的主要指标。由表可知, NB(ECFP_6)模型训练集的五倍交叉验证中AUC均大于0.996, MCC均大于0.977;测试集验证中AUC均大于0.992, MCC均大于0.778。RP(MACCS)模型训练集的五倍交叉验证中AUC均大于0.976, MCC均大于0.869;测试集验证中除F11靶点AUC为0.746外, AUC均大于0.928;除F11靶点MCC为0.583外, MCC均大于0.710。因此, 剔除F11靶点的递归分隔模型, 所有模型均具有较好的有效性和精确度, 可以用于消栓通络方中化学成分作用靶点的预测。

表 3 Performance of NB (ECFP_6) model's training set by 5-fold cross-validation. TP: True positive; FN: False negative; FP: False positive; TN: True negative; SE: Sensitivity; SP: Specificity; MCC: Matthews correlation coefficient; AUC: The area under the receiver operating characteristic curve

表 4 Performance of NB (ECFP_6) model's test set validation

表 5 Performance of RP (MACCS) model's training set by 5 fold cross-validation

表 6 Performance of RP (MACCS) model's test set validation
4 分子对接模型

应用Discovery Studio 2018中的CDocker分子对接方法, 对蛋白与蛋白配体进行对接。对接结果通过RMSD进行验证, 两个靶点THBD和P2RY12所构建模型RMSD值分别为0.758 0和0.656 6, RMSD值均小于2.5, 认为分子对接模型有效。蛋白晶体结构与蛋白原始配体相互作用如图 2所示。

图 2 The crystal structure and the interaction between proteins and their original ligands. A: Interaction between THBD and its ligand in their three dimensional (3D) structures, Graph a shows their interaction in two dimensional (2D) structures; B: Interaction between P2RY12 and its ligand in their 3D structures, Graph b shows their interaction in 2D structures
5 靶点预测结果

基于上述建立的所有机器学习模型和分子对接模型, 分别对消栓通络方中的1 251个化学成分进行虚拟筛选, 以期预测发现与靶点相互作用的活性成分。对预测结果进行统计和整理, 共得到153个活性化学成分, 其中22个化学成分可以与2个以上靶点相互作用。

6 化学成分-靶点和靶点-通路作用网络

选择活性化学成分最多的一味中药丹参, 利用Cytoscape软件构建其化学成分-靶点作用网络, 如图 3所示。在DAVID数据库中, 对脑卒中靶点进行GO富集分析, 得到靶点的代谢通路, 选择P值小于0.05的通路, 在此基础上, 构建靶点-通路网络, 如图 4所示。

图 3 The chemical constituent-target network. The rectangle means active constituents from Radix Salviae, and the oval means stroke related targets. Danshen01: (2R)-3-(3, 4-Dihydroxyphenyl)-2-[(Z)-3-(3, 4-dihydroxyphenyl)acryloyl]oxy-propionic acid; Danshen02: (2S)-2-Amino-3-[(2R)-2-amino-3-hydroxy-3-oxopropyl]disulfanylpropanoic acid; Danshen03: (2S)-3-(3, 4-Dihydroxyphenyl)-2-hydroxypropanoic acid; Danshen04: (2S, 3S)-2-(3, 4-Dihydroxyphenyl)-7-hydroxy-4-[(E)-3-hydroxy-3-oxoprop-1-enyl]-2, 3-dihydrobenzofuran-3-carboxylic acid; Danshen05: (E)-3-[2-(3, 4-Dihydroxyphenyl)-7-hydroxy-benzofuran-4-yl]acrylic acid; Danshen06: (Z)-3-[2-[(E)-2-(3, 4-Dihydroxyphenyl)vinyl]-3, 4-dihydroxy-phenyl]acrylic acid; Danshen07: 1-(3, 4-Dihydroxyphenyl)-2-hydroxyethanone; Danshen08: 6-O-Syringyl-8-O-acetyl shanzhiside methyl ester; Danshen09: (2S)-2-Amino-5-(diaminomethylideneamino) pentanoic acid

图 4 The target-pathway network based on the enrichment analysis of Gene-Ontology (GO) biological pathway
讨论

缺血性脑卒中的理想治疗模式包括血管再通和神经保护, 去除脑血栓的同时, 保护神经减轻缺血再灌注造成的损害[25], 这也意味着脑卒中的治疗需要多靶点介入。因此本研究开展网络药理学机制分析, 力图发现治疗脑卒中的多靶点有效成分组, 多个有效成分通过多靶点多通路作用于脑卒中, 以期得到更好的疗效。

本文的创新之处在于首次针对脑卒中相关药物靶点, 系统地建立了化合物-药物靶点相互作用预测模型, 并通过验证具有很好的预测能力。此外, 应用构建的化合物-靶点相互作用预测模型, 得到了消栓通络方中的活性成分, 并利用网络药理学方法, 建立了活性成分-药物靶点-通路网络, 从而揭示了消栓通络方抗脑卒中的药效物质基础和作用机制。

本研究通过机器学习算法朴素贝叶斯和递归分隔以及Cdocker分子对接方法, 对脑卒中临床应用经典处方消栓通络方中的化学成分进行作用靶点预测, 得到不错的结果。其中有3个小分子值得特别留意: ①姜黄素, 是一种植物多酚, 其具有广泛的生物学效应而引起广泛关注。本研究构建的预测模型产生的活性化合物中, 共有4个姜黄素类化合物入选, 分别是姜黄素、1, 2-二氢姜黄素、去甲氧基姜黄素[26]和双去甲氧基姜黄素, 对MMP9具有潜在活性, 提示姜黄素类化合物对脑卒中可能的治疗作用。许多研究表明, 姜黄素对脑卒中有治疗作用, 其作用机制包括抗氧化、清除氧自由基的作用、抗炎作用、降低胆固醇作用和抗凝作用[3]。此外, 还有研究表明, 姜黄素可能通过对表观遗传的影响达到抗脑卒中的作用[27]。另外, Li等[28]研究显示姜黄素通过下调炎症(NF-κB、ICAM-1、MMP-9)和凋亡(caspase-3)相关标志物, 保护脑缺血再灌注造成的损伤, 达到治疗作用。上述文献[26, 27]报道进一步确证了本研究所建模型预测结果的可靠性, 也说明姜黄素类化学成分可以作为脑卒中治疗的重要组成部分。②肉桂醛大量存在于肉桂等植物体内, 研究表明其具有抗炎、抗氧化、抗血小板聚集等多种生物学作用。模型预测其对F11、MMP9、TBXA2R和VEGFA四个靶点具有潜在活性, 提示肉桂醛对脑卒中潜在的治疗作用。Zhao等[29]研究显示, 肉桂醛可以通过抑制TLR4/TRAF6/NF-κB炎症通路, 减轻缺血脑组织的炎症反应, 从而在缺血再灌注损伤后对脑产生保护作用。Chen等[30]类似研究显示, 肉桂醛通过减弱iNOS, COX-2表达和NF-κB信号通路, 抑制神经炎症, 从而对脑产生保护作用。此外, Ballabeni等[31]以Ocotea quixos精油为实验药物, 其主要成分为反式肉桂醛, 结果显示肉桂醛具有抗血小板作用和血管舒张作用, 抑制小鼠血栓的形成。由文献[29-31]可以发现, 肉桂醛可以通过血管再通和神经保护两种方式治疗脑卒中, 是潜在优异的抗脑卒中的活性成分。③丹酚酸为丹参水溶性主要有效成分[32], 为一类结构相似的酚酸类化合物。本研究构建模型预测结果中, 多种丹酚酸入选, 包括丹酚酸A、B、C、D、E、G、J、N和异丹酚酸C, 共同作用于MMP9、SERPINE1、F11和P2RY12四个靶点, 提示丹酚酸对脑卒中可能的治疗作用。目前已有许多研究报道丹酚酸对脑损伤的保护作用, 且临床上已应用丹酚酸注射液用于脑卒中的治疗。其中, Mahmood等[33]研究显示, 丹酚酸A通过抑制eNOS解偶联和ONOO-形成神经血管保护作用。另外, 丹酚酸A可减弱血影蛋白和钙调神经磷酸酶的分解, 从而保护大脑免受缺血/再灌注损伤, 进而达到治疗脑卒中的作用。Fan等[34]研究显示, 丹酚酸B可以通过降低GFAP、Iba1、IL-1β、IL-6、TNF-α和cleaved-caspase 3的表达, 抑制星形胶质细胞和小胶质细胞的过度活化, 达到脑卒中治疗中对神经细胞的保护作用。由此可知, 多种丹酚酸化合物都存在治疗脑卒中的潜在能力, 预测得到的多种丹酚酸化合物值得同样关注。

除受到关注广为研究的化学成分外, 预测所得活性化学成分中还包含其他一些典型化学成分结构, 经整合分析, 如表 7示。

表 7 The typical examples of active constituents

综上所述, 本研究利用机器学习等算法, 借助网络药理学研究方法和技术, 构建了化学成分-靶点-通路网络, 对消栓通络方的药效物质基础和网络作用机制进行了较为深入的研究。本研究发现了许多潜在的多靶点活性成分通过血管再通和神经保护多种方式治疗脑卒中, 为消栓通络方的临床应用和后续研究提供重要信息, 同时还发现了许多潜在的多靶点活性成分用于缓解脑卒中的症状并治疗脑卒中, 也为后续开展的中药复方抗脑卒中研究提供参考。

参考文献
[1]
Jin BX. Clinical study on the prognostic factors of thrombolysis with alteplase in acute ischemic stroke[J]. Chin Heal Ind (中国卫生产业), 2013, 10: 101-102.
[2]
National institute of neurological disorders and stroke rt-PA stroke study group. Tissue plasminogen activator for acute ischemic stroke[J]. New Engl J Med, 1995, 333: 1581-1588.
[3]
Liao CY. Study on therapeutic targets and drugs of stroke[J]. Chin J Pharm Toxicol (中国药理学与毒理学杂志), 2016, 30: 1264-1272.
[4]
Tan CB, Gao M, Wang HQ, et al. Study on the selection of the effective components of Xiaoshuan Tongluo Recipe[J]. Pharmacol Clin Chinese Mater Med (中药药理与临床), 2011, 27: 75-78.
[5]
Sorensen SV, Dewilde S, Singer DE, et al. Cost-effectiveness of warfarin:trial versus "real-world" stroke prevention in atrial fibrillation[J]. Am Heart J, 2009, 157: 1064-1073. DOI:10.1016/j.ahj.2009.03.022
[6]
Grysiewicz RA, Thomas K, Pandey DK. Epidemiology of ischemic and hemorrhagic stroke:incidence, prevalence, mortality, and risk factors[J]. Neurol Clin, 2008, 26: 871-895. DOI:10.1016/j.ncl.2008.07.003
[7]
Zhang BY, Zheng YF, Pang XC, et al. Study on the material basis and network mechanism of anti-tumor effect of traditional Chinese Medicine[J]. World J Tradit Chin Med (世界中医药), 2018, 13: 1997-2009.
[8]
Xie T, Song S, Li S, et al. Review of natural product databases[J]. Cell Prolif, 2015, 48: 398-404. DOI:10.1111/cpr.12190
[9]
Bachurin SO, Bovina EV, Ustyugov AA. Drugs in clinical trials for Alzheimer's disease:the major trends[J]. Med Res Rev, 2017, 37: 1186-1225. DOI:10.1002/med.21434
[10]
Chen X, Ji ZL, Chen YZ. TTD:therapeutic target database[J]. Nucleic Acids Res, 2002, 30: 412-415. DOI:10.1093/nar/30.1.412
[11]
Liu L, Du B, Zhang H, et al. A network pharmacology approach to explore the mechanisms of Erxian decoction in polycystic ovary syndrome[J]. Chin Med, 2018, 13: 46. DOI:10.1186/s13020-018-0201-1
[12]
Lee AY, Park W, Kang TW, et al. Network pharmacology-based prediction of active compounds and molecular targets in Yijin-Tang acting on hyperlipidaemia and atherosclerosis[J]. J Ethnopharmacol, 2018, 221: 151-159. DOI:10.1016/j.jep.2018.04.027
[13]
Friedman N, Geiger D, Goldszmidt M. Bayesian network classifiers[J]. Mach Learn, 1997, 29: 131-163. DOI:10.1023/A:1007465528199
[14]
Stone M. Cross-validation and multinomial prediction[J]. Biometrika, 1974, 61: 509-515. DOI:10.1093/biomet/61.3.509
[15]
Rogers D, Brown RD, Hahn M. Using extended-connectivity fingerprints with Laplacian-modified Bayesian analysis in high-throughput screening follow-up[J]. J Biomol Screen, 2005, 10: 682-686. DOI:10.1177/1087057105281365
[16]
Cereto-Massagué A, Ojeda MJ, Valls C, et al. Molecular fingerprint similarity search in virtual screening[J]. Methods, 2015, 71: 58-63. DOI:10.1016/j.ymeth.2014.08.005
[17]
Fang J, Li Y, Liu R, et al. Discovery of multitarget-directed ligands against Alzheimer's disease through systematic prediction of chemical-protein interactions[J]. J Chem Inform Model, 2015, 55: 149-164. DOI:10.1021/ci500574n
[18]
Kirchmair J, Markt P, Distinto S, et al. Evaluation of the performance of 3D virtual screening protocols:RMSD comparisons, enrichment assessments, and decoy selection-what can we learn from earlier mistakes?[J]. J Comput Aid Mol Design, 2008, 22: 213-228. DOI:10.1007/s10822-007-9163-6
[19]
Hothorn T, Hornik K, Zeileis A. Unbiased recursive partitioning:a conditional inference framework[J]. J Comput Graph Statist, 2006, 15: 651-674. DOI:10.1198/106186006X133933
[20]
Shannon P, Markiel A, Ozier O, et al. Cytoscape:a software environment for integrated models of biomolecular interaction networks[J]. Genome Res, 2003, 13: 2498-2504. DOI:10.1101/gr.1239303
[21]
Ashburner M, Ball CA, Blake JA, et al. Gene Ontology:tool for the unification of biology[J]. Nat Gene, 2000, 25: 25-29. DOI:10.1038/75556
[22]
Pang XC, Wang Z, Fang JS, et al. Network pharmacology study of effective constituents of traditional Chinese medicine for Alzheimer's disease treatment[J]. Acta Pharm Sin (药学学报), 2016, 51: 725-731.
[23]
Xu L, Jiang W, Pang X, et al. Network pharmacology study of the effective constituents in the Compound Yizhihao against influenza disease​[J]. Acta Pharm Sin (药学学报), 2017, 52: 745-752.
[24]
Bajusz D, Rácz A, Héberger K. Why is Tanimoto index an appropriate choice for fingerprint-based similarity calculations[J]. J Cheminform, 2015, 7: 20. DOI:10.1186/s13321-015-0069-3
[25]
Nour M, Scalzo F, Liebeskind DS. Ischemia-reperfusion injury in stroke[J]. Int Neurol, 2012, 1: 185-199.
[26]
Ahmad N, Umar S, Ashafaq M, et al. A comparative study of PNIPAM nanoparticles of curcumin, demethoxycurcumin, and bisdemethoxycurcumin and their effects on oxidative stress markers in experimental stroke[J]. Protoplasma, 2013, 250: 1327-1338. DOI:10.1007/s00709-013-0516-9
[27]
Chen YH, Tang YH. The epigenetic effect of curcumin on stroke prevention[J]. Mod Med J China (中国现代医药杂志), 2016, 18: 98-101.
[28]
Li W, Suwanwela NC, Patumraj S. Curcumin prevents reperfusion injury following ischemic stroke in rats via inhibition of NF κB, ICAM-1, MMP-9 and caspase-3 expression[J]. Mol Med Rep, 2017, 16: 4710-4720. DOI:10.3892/mmr.2017.7205
[29]
Zhao J, Zhang X, Dong L, et al. Cinnamaldehyde inhibits inflammation and brain damage in a mouse model of permanent cerebral ischaemia[J]. Br J Pharmacol, 2015, 172: 5009-5023. DOI:10.1111/bph.13270
[30]
Chen YF, Wang YW, Huang WS, et al. Trans-cinnamaldehyde, an essential oil in cinnamon powder, ameliorates cerebral ischemia-induced brain injury via inhibition of neuroinflammation through attenuation of iNOS, COX-2 expression and NF-κB signaling pathway[J]. Neuromol Med, 2016, 18: 322-333. DOI:10.1007/s12017-016-8395-9
[31]
Ballabeni V, Tognolini M, Bertoni S, et al. Antiplatelet and antithrombotic activities of essential oil from wild Ocotea quixos (Lam.) Kosterm.(Lauraceae) calices from Amazonian Ecuador[J]. Pharmacol Res, 2007, 55: 23-30. DOI:10.1016/j.phrs.2006.09.009
[32]
Du GH, Zhang JT. Research progress of salvianolic acid, an effective water-soluble component of Salvia miltiorrhiza[J]. Basic Clin Med (基础医学与临床), 2000, 20: 10-14.
[33]
Mahmood Q, Wang GF, Wu G, et al. Salvianolic acid A inhibits calpain activation and eNOS uncoupling during focal cerebral ischemia in mice[J]. Phytomedicine, 2017, 25: 8-14. DOI:10.1016/j.phymed.2016.12.004
[34]
Fan Y, Luo Q, Wei J, et al. Mechanism of salvianolic acid B neuroprotection against ischemia/reperfusion induced cerebral injury[J]. Brain Res, 2018, 1679: 125-133. DOI:10.1016/j.brainres.2017.11.027