摘要:
目的:基于中药系统药理学方法探讨六味地黄丸治疗骨质疏松症(osteoporosis,OP)的药效成分、作用靶点及作用特点.方法:在中药综合数据库中检索六味地黄丸中的已知化学成分,然后使用中医药分子机制生物信息学分析工具预测六味地黄丸已知化学成分的作用靶点.通过检索DisGeNET数据库和人工阅读文献,查找OP相关靶基因和蛋白,通过与预测得到的作用靶点进行交叠,筛选六味地黄丸治疗OP的作用靶点.利用Cytoscape 3.7软件自带的BiNGO和ClueGO插件,分别对筛选出的靶点进行基因注释和信号通路富集分析.最后选取前期构建的绝经后OP疾病模型中的87个差异表达蛋白(基于血清蛋白质组学筛选的绝经后OP女性与绝经后骨量正常女性的血清差异表达蛋白),与筛选出的六味地黄丸治疗OP的作用靶点进行交叠,筛选六味地黄丸治疗OP的关键靶点,通过分析原始数据,确定相应的中药及其药效成分.结果:通过检索发现六味地黄丸化学成分130个,其中熟地黄化学成分8个、山药化学成分20个、山茱萸化学成分46个、泽泻化学成分21个、牡丹皮化学成分18个、茯苓化学成分21个;茯苓与山茱萸、山茱萸与牡丹皮各含有1个相同化学成分,泽泻与山茱萸含有2个相同化学成分.根据六味地黄丸化学成分共预测到1169个作用靶点,其中君药熟地黄作用靶点42个,臣药山茱萸和山药作用靶点997个,佐药泽泻、牡丹皮和茯苓作用靶点共622个;君药与臣药共有12个共同作用靶点,君药与佐药有14个共同作用靶点,臣药与佐药有472个共同作用靶点,君药、臣药、佐药共有6个共同作用靶点.经检索DisGeNET数据库和人工阅读文献,共检索到OP相关靶点687个,通过与预测到的六味地黄丸作用靶点交叠,共筛选出六味地黄丸治疗OP的作用靶点142个.靶点基因注释发现,六味地黄丸治疗OP的142个作用靶点主要参与了20种生物学过程,分布于8种细胞成分中,具有10种分子功能,其作用涉及骨重建、骨化、血管重塑、软骨发育、肌肉发育、骨骼发育、维生素D代谢、维生素应答、雌激素刺激应答、破骨细胞分化、成骨细胞分化、骨矿化、Wnt信号通路和Samd信号通路等生物学行为的调控.信号通路富集分析发现,119个信号通路被显著富集,破骨细胞分化、NF-κB信号通路和卵巢类固醇生成等24个功能组被显著富集,其中Wnt信号通路和与绝经后OP密切相关的雌激素调控信号通路等在24个功能组中广泛分布.通过将87个差异表达蛋白与筛选出的六味地黄丸治疗OP的作用靶点交叠,筛选出血管紧张素原、载脂蛋白E、Dickkopf相关蛋白3、RAS相关蛋白7a和蛋白质二硫键异构酶共5个关键靶点,对应山茱萸、山药、茯苓的6个药效成分,包括马兜铃酮、α-考绕咖烯、豆甾醇、尿囊素、鞘氨醇A及齿孔醇.结论:六味地黄丸治疗OP的药效成分包括山药、山茱萸和茯苓所含化合物马兜铃酮、α-考绕咖烯、豆甾醇、尿囊素、鞘氨醇A及齿孔醇,主要作用靶点包括血管紧张素原、载脂蛋白E、Dickkopf相关蛋白3、RAS相关蛋白7 a和蛋白质二硫键异构酶.六味地黄丸治疗OP具有多靶点、多系统协同作用的特点.