基于网络药理学探讨缓衰方对慢性肾脏病心肾保护作用的机制

基于网络药理学探讨缓衰方对慢性肾脏病心肾保护作用的机制
作者：刘童童王宇阳王倩冒慧敏占永立
来源：《世界中医药》2021年第20期
摘要目的：探討缓衰方对慢性肾脏病（CKD）心肾保护作用的机制。

方法：基于中药系统药理学数据库与分析平台（TCMSP）检索缓衰方中5味中药的有效活性成分并获得相应的预测靶点;通过GEO数据库、GeneCards数据库和人类孟德尔遗传数据库（OMIM）检索与CKD和心血管疾病（CVD）相关的靶点，并对药物靶点和疾病靶点进行映射取交集作为缓衰方对CKD患者心肾保护作用的关键靶点。

基于String数据库对关键靶点进行蛋白质-蛋白质相互作用（PPI）网络分析，并通过拓扑分析得到缓衰方对CKD患者心肾保护作用的核心靶点，基于autodock软件对缓衰方的活性成分和核心靶点作分子对接，分析成分和靶点的结合能力。

并对关键靶点进行基因本体（GO）富集分析和京都基因和基因组百科全书（KEGG）通路富集分析。

结果：共得到缓衰方有效活性成分105种，预测靶点135个，得到CKD疾病相关靶点1 436个，CVD疾病相关靶点744个。

共得到缓衰方对CKD心肾保护作用的关键靶点53个，拓扑分析得到核心靶点18个，分子对接结果显示缓衰方的11个活性成分和10个核心靶点平均最低结合能为-6.58 kcal/mol，其中豆甾醇、槲皮素与PTGS2靶点蛋白结合能最高。

共得到对氧含量的反应、细胞对药物的反应、循环系统中的血管过程等1 275个GO富集结果和低氧诱导因子-1（HIF-1）信号通路、肿瘤坏死因子（TNF）信号通路、白细胞介素-17（IL-17）信号通路等76条KEGG信号通路。

结论：缓衰方通过多成分、多靶点、多信号通路途径实现对CKD的心肾保护作用。

关键词网络药理学;缓衰方;慢性肾脏病;心血管疾病;活性成分;靶点;信号通路;作用机制
Abstract Objective：To explore the mechanism of the heart and kidney protective effect of Huan Shuai Formula on the chronic kidney disease（CKD）.Methods：We searched for targets related to CKD and Cardiovascular Disease（CVD） through GEO database，GeneCards database and Human Mendelian Inheritance Database（OMIM），and map drug targets and disease targets to take the intersection as a slowing down prescription as the key target of cardio-renal protection in CKD patients.The protein-protein interaction（PPI） network analysis of key targets was based on the String database，and topological analysis was used to obtain the core targets of the heart and kidney
protection of CKD patients.Based on the autodock software，molecular docking of the active ingredients and core targets of the Huan Shuai formula was performed，and the binding ability of the ingredients and the targets was analyzed.And we carried out GO（Gene Ontology） biological process enrichment analysis and Genome Encyclopedia（KEGG） pathway enrichment analysis for key targets.Results：A total of 105 active components，135 predicted targets，1 436 CKD disease related targets and 744 CVD disease related targets were obtained.A total of 53 key targets for the cardio-renal protection of CKD by Yan shuai Formula were obtained，and 18 core targets were obtained by topological analysis.The molecular docking results showed that the average minimum binding energy of 11 active ingredients and 10 core targets of Yan Shuai Formula was-6.58
kcal/mol，of which stigmasterol and quercetin had the highest binding energy to PTGS2 target protein.A total of 1275 GO enrichment results，hypoxia inducible Factor-1（HIF-1） signaling pathway，Tumor Necrosis Factor（TNF） signaling pathway，Interleukin-17（Interleukin-17，IL-17），the 76 KEGG signal pathways were obtained.Conclusion：Huan Shuai Formula can protect the heart and kidney function of CKD through multi-component，multi-target and multi signal pathway.
Keywords Network pharmacology; Huan Shuai Formula; Chronic kidney disease; Cardiovascular disease; Active components; Targets; Signaling pathway; Mechanism
中圖分类号：R242;R692文献标识码：Adoi：10.3969/j.issn.1673-7202.2021.20.004
慢性肾脏病（Chronic Kidney Disease，CKD）已成为全球共同关注的公共卫生问题，全世界超过13.4%的成年人患有CKD[1]。

在我国，成年人CKD的发病率约为10.8%[2]。

心血管疾病（Cardiovascular Disease，CVD）为CKD常见的并发症，超过64.5%的CKD患者患有不同程度的CVD，约为健康人群的2～4倍[3-4]。

有研究表明，随着CKD的进展，CVD的发生风险显著增加，终末期肾脏疾病（End Stage Renal Disease，ESRD）患者CVD的发生率可达80%[5]。

CKD患者血脂异常、血压控制不佳等传统CVD危险因素导致动脉硬化和心肌重构是CVD形成的基本病变，临床主要表现为冠状动脉粥样硬化性心脏病、心律失常、慢性心力衰竭等[6-8]。

CKD肾素-血管紧张素-醛固酮系统（Renin-angiotensin-aldosterone System，RAAS）激活、蛋白尿、肾小球滤过率下降、贫血、全身营养不良等非传统CKD危险因素，可加速CVD进展;而心输出量的减少和静脉压的增加，使肾脏灌注减少，导致肾功能丢失加速。

CKD 和CVD的相互影响使伴有CVD的CKD患者死亡率高达50%，为健康人群的10～30倍，是CKD死亡的首要因素[9]。

由此可见，CKD患者防治CVD至关重要。

缓衰方以黄芪、丹参、当归、大黄、白花蛇舌草为核心组成，具有益气活血、解毒降浊之效。

临床和实验研究发现，缓衰方能够有效改善CKD患者的生命质量，延缓CKD患者进入终点事件（透析、全因死亡、血清肌酐双倍增加）的时间[10-11]。

在长期临床实践中我们发现，缓衰方可显著减轻CKD患者胸闷、胸痛、心悸等CVD症状。

然而目前关于缓衰方对CKD心肾保护作用的有效活性成分、作用靶点和作用机制尚未完全阐明。

因此，本研究拟采用网络药
理学的方法从分子水平方面研究和分析缓衰方对CKD患者心肾保护作用的有效活性成分、潜在靶点和作用机制，以期为缓衰方的应用和开发提供科学依据。

1 资料与方法
1.1 缓衰方活性成分及靶点的筛选通过检索中药系统药理学数据库与分析平台（TCMSP，http：///tcmsp.php），收集缓衰方中“黄芪”“丹参”“当归”“大黄”“白花蛇舌草”5味中药的活性成分，根据活性成分的药物毒代动力学参数进行筛选，以类药性（Drug-like，DL）≥0.18和口服生物利用度（Oral Bioavailability，OB）≥30%的化合物作为缓衰方的有效活性成分。

同时根据TCMSP数据库整理有效活性成分所涉及的靶点，并通过Uniprot数据库（https：//）进行靶点信息比对和靶点名校正。

1.2 CKD和CVD疾病相关靶点的检索 GEO数据库是全球最大的公共基因表达数据库，具有较全面的高通量分子数据。

以关键词“chronic kidney disease”在GEO数据库中检索与CKD相关的靶点，选择物种为人类（Homo Sapiens），对得到的基因靶点进行整理并找出差异基因靶点。

同时分别在GeneCards数据库（https：//）和人类孟德尔遗传数据库（OMIM，https：//）检索与CKD和CVD相关的疾病靶点，保留相关评分（Relevance Score）≥25的疾病相关靶点。

1.3 蛋白质-蛋白质相互作用网络构建利用半编程性质的脚本性软件R project（×64 3.6.1）对缓衰方活性成分所预测的靶点和CKD、CVD相关的疾病靶点进行映射，获得3个靶点数据的交集作为缓衰方对CKD心肾保护作用的关键靶点进行机制探讨。

使用Bisogenet数据库（http：//）构建关键靶点的蛋白质-蛋白质相互作用（Protein Protein Interaction，PPI）网络，设定物种为人类（Homo Sapiens），基于Cytoscape（Version 3.7.2）软件的cytNCA插件计算PPI网路中各节点的度中心性（Degree Centrality，DC）、中介中心性（Betweenness Centrality，BC）和紧密中心性（Closeness Centrality，CC）等信息并进行拓扑分析，选取DC、BC和CC值均高于均值的靶点作为核心靶点。

基于Pubchem数据库（https：//）查找缓衰方药物活性成分的3D结构，并进行能量最小化处理，利用RCSB PDB数据库（http：//）查找核心靶点的相关蛋白晶体结构，利用PyMOL（Version 1.7.2）软件对靶点蛋白进行去水、加氢处理，并以对心肾血管具有保护作用的厄贝沙坦（Irbesartan）作为阳性对照，基于autodock（Version 1.5.6）软件对缓衰方药物活性成分和靶点相关蛋白进行分子模拟对接，并计算最低结合能（Binding Affinity），最低结合能≤-5 kcal/mol则说明药物活性成分和靶点蛋白对接较好，结合能的绝对值越大说明活性成分与靶点蛋白结合越好[12]。

1.4 网络分析基于bioconductor数据库（http：//）对缓衰方心肾同治作用的关键靶点进行基因本体（Gene Ontology，GO）生物学过程富集分析和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集分析，保留
P<0.05，Q<0.05的富集结果，获得关键靶点参与的生物过程（Biological Process，BP）、细胞
组分（Cellular Component，CC）、分子功能（Molecular Function，MF）和作用信号通路。

基于Cytoscape软件，构建“成分-靶点-通路”网络，对缓衰方活性成分、靶点和疾病之间的相互
关系进行分析。

Keywords Network pharmacology; Huan Shuai Formula; Chronic kidney disease; Cardiovascular disease; Active components; Targets; Signaling pathway; Mechanism
中图分类号：R242;R692文献标识码：Adoi：10.3969/j.issn.1673-7202.2021.20.004
慢性肾脏病（Chronic Kidney Disease，CKD）已成为全球共同关注的公共卫生问题，全世界超过13.4%的成年人患有CKD[1]。

在我国，成年人CKD的发病率约为10.8%[2]。

心血管疾病（Cardiovascular Disease，CVD）为CKD常见的并发症，超过64.5%的CKD患者患有不同程度的CVD，约为健康人群的2～4倍[3-4]。

有研究表明，随着CKD的进展，CVD的发生风险显著增加，终末期肾脏疾病（End Stage Renal Disease，ESRD）患者CVD的发生率可达80%[5]。

CKD患者血脂异常、血压控制不佳等传统CVD危险因素导致动脉硬化和心肌重构是CVD形成的基本病变，临床主要表现为冠状动脉粥样硬化性心脏病、心律失常、慢性心力衰竭等[6-8]。

CKD肾素-血管紧张素-醛固酮系统（Renin-angiotensin-aldosterone System，RAAS）激活、蛋白尿、肾小球滤过率下降、贫血、全身营养不良等非传统CKD危险因素，可加速CVD进展;而心输出量的减少和静脉压的增加，使肾脏灌注减少，导致肾功能丢失加速。

CKD 和CVD的相互影响使伴有CVD的CKD患者死亡率高达50%，为健康人群的10～30倍，是CKD死亡的首要因素[9]。

由此可见，CKD患者防治CVD至关重要。

缓衰方以黄芪、丹参、当归、大黄、白花蛇舌草为核心组成，具有益气活血、解毒降浊之效。

临床和实验研究发现，缓衰方能够有效改善CKD患者的生命质量，延缓CKD患者进入终点事件（透析、全因死亡、血清肌酐双倍增加）的时间[10-11]。

在长期临床实践中我们发现，缓衰方可显著减轻CKD患者胸闷、胸痛、心悸等CVD症状。

然而目前关于缓衰方对CKD心肾保护作用的有效活性成分、作用靶点和作用机制尚未完全阐明。

因此，本研究拟采用网络药理学的方法从分子水平方面研究和分析缓衰方对CKD患者心肾保护作用的有效活性成分、潜在靶点和作用机制，以期为缓衰方的应用和开发提供科学依据。

1 资料与方法
1.1 缓衰方活性成分及靶点的筛选通过检索中药系统药理学数据库与分析平台（TCMSP，http：///tcmsp.php），收集缓衰方中“黄芪”“丹参”“当归”“大黄”“白花蛇舌草”5味中药的活性成分，根据活性成分的药物毒代动力学参数进行筛选，以类药性（Drug-like，DL）≥0.18和口服生物利用度（Oral Bioavailability，OB）≥30%的化合物作为缓衰方的有效活性成分。

同时根据TCMSP数据库整理有效活性成分所涉及的靶点，并通过Uniprot数据库（https：//）进行靶点信息比对和靶点名校正。

1.2 CKD和CVD疾病相关靶点的检索 GEO数据库是全球最大的公共基因表达数据库，具有较全面的高通量分子数据。

以关键词“chronic kidney disease”在GEO数据库中检索与CKD相关的靶点，选择物种为人类（Homo Sapiens），对得到的基因靶点进行整理并找出差异基因靶点。

同时分别在GeneCards数据库（https：//）和人类孟德尔遗传数据库（OMIM，https：//）检索与CKD和CVD相关的疾病靶点，保留相关评分（Relevance Score）≥25的疾病相关靶点。

1.3 蛋白质-蛋白质相互作用网络构建利用半编程性质的脚本性软件R project（×64 3.6.1）对缓衰方活性成分所预测的靶点和CKD、CVD相关的疾病靶点进行映射，获得3个靶点数据的交集作为缓衰方对CKD心肾保护作用的关键靶点进行机制探讨。

使用Bisogenet数据库（http：//）构建关键靶点的蛋白质-蛋白质相互作用（Protein Protein Interaction，PPI）网络，设定物种为人类（Homo Sapiens），基于Cytoscape（Version 3.7.2）软件的cytNCA插件计算PPI网路中各节点的度中心性（Degree Centrality，DC）、中介中心性（Betweenness Centrality，BC）和紧密中心性（Closeness Centrality，CC）等信息并进行拓扑分析，选取DC、BC和CC值均高于均值的靶点作为核心靶点。

基于Pubchem数据库（https：//）查找缓衰方药物活性成分的3D结构，并进行能量最小化处理，利用RCSB PDB數据库（http：//）查找核心靶点的相关蛋白晶体结构，利用PyMOL（Version 1.7.2）软件对靶点蛋白进行去水、加氢处理，并以对心肾血管具有保护作用的厄贝沙坦（Irbesartan）作为阳性对照，基于autodock（Version 1.5.6）软件对缓衰方药物活性成分和靶点相关蛋白进行分子模拟对接，并计算最低结合能（Binding Affinity），最低结合能≤-5 kcal/mol则说明药物活性成分和靶点蛋白对接较好，结合能的绝对值越大说明活性成分与靶点蛋白结合越好[12]。

1.4 网络分析基于bioconductor数据库（http：//）对缓衰方心肾同治作用的关键靶点进行基因本体（Gene Ontology，GO）生物学过程富集分析和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集分析，保留
P<0.05，Q<0.05的富集结果，获得关键靶点参与的生物过程（Biological Process，BP）、细胞组分（Cellular Component，CC）、分子功能（Molecular Function，MF）和作用信号通路。

基于Cytoscape软件，构建“成分-靶点-通路”网络，对缓衰方活性成分、靶点和疾病之间的相互
关系进行分析。

Keywords Network pharmacology; Huan Shuai Formula; Chronic kidney disease; Cardiovascular disease; Active components; Targets; Signaling pathway; Mechanism
中图分类号：R242;R692文献标识码：Adoi：10.3969/j.issn.1673-7202.2021.20.004
慢性肾脏病（Chronic Kidney Disease，CKD）已成为全球共同关注的公共卫生问题，全世界超过13.4%的成年人患有CKD[1]。

在我国，成年人CKD的发病率约为10.8%[2]。

心血管疾病（Cardiovascular Disease，CVD）为CKD常见的并发症，超过64.5%的CKD患者患有不同
程度的CVD，约为健康人群的2～4倍[3-4]。

有研究表明，随着CKD的进展，CVD的发生风险显著增加，终末期肾脏疾病（End Stage Renal Disease，ESRD）患者CVD的发生率可达80%[5]。

CKD患者血脂异常、血压控制不佳等传统CVD危险因素导致动脉硬化和心肌重构是CVD形成的基本病变，临床主要表现为冠状动脉粥样硬化性心脏病、心律失常、慢性心力衰竭等[6-8]。

CKD肾素-血管紧张素-醛固酮系统（Renin-angiotensin-aldosterone System，RAAS）激活、蛋白尿、肾小球滤过率下降、贫血、全身营养不良等非传统CKD危险因素，可加速CVD进展;而心输出量的减少和静脉压的增加，使肾脏灌注减少，导致肾功能丢失加速。

CKD 和CVD的相互影响使伴有CVD的CKD患者死亡率高达50%，为健康人群的10～30倍，是CKD死亡的首要因素[9]。

由此可见，CKD患者防治CVD至关重要。

缓衰方以黄芪、丹参、当归、大黄、白花蛇舌草为核心组成，具有益气活血、解毒降浊之效。

临床和实验研究发现，缓衰方能够有效改善CKD患者的生命质量，延缓CKD患者进入终点事件（透析、全因死亡、血清肌酐双倍增加）的时间[10-11]。

在长期临床实践中我们发现，缓衰方可显著减轻CKD患者胸闷、胸痛、心悸等CVD症状。

然而目前关于缓衰方对CKD心肾保护作用的有效活性成分、作用靶点和作用机制尚未完全阐明。

因此，本研究拟采用网络药理学的方法从分子水平方面研究和分析缓衰方对CKD患者心肾保护作用的有效活性成分、潜在靶点和作用机制，以期为缓衰方的应用和开发提供科学依据。

1 资料与方法
1.1 缓衰方活性成分及靶点的筛选通过检索中药系统药理学数据库与分析平台（TCMSP，http：///tcmsp.php），收集缓衰方中“黄芪”“丹参”“当归”“大黄”“白花蛇舌草”5味中药的活性成分，根据活性成分的药物毒代动力学参数进行筛选，以类药性（Drug-like，DL）≥0.18和口服生物利用度（Oral Bioavailability，OB）≥30%的化合物作为缓衰方的有效活性成分。

同时根据TCMSP数据库整理有效活性成分所涉及的靶点，并通过Uniprot数据库（https：//）进行靶点信息比对和靶点名校正。

1.2 CKD和CVD疾病相关靶点的检索 GEO数据库是全球最大的公共基因表达数据库，具有较全面的高通量分子数据。

以关键词“chronic kidney disease”在GEO数据库中检索与CKD相关的靶点，选择物种为人类（Homo Sapiens），对得到的基因靶点进行整理并找出差异基因靶点。

同时分别在GeneCards数据库（https：//）和人类孟德尔遗传数据库（OMIM，https：//）检索与CKD和CVD相关的疾病靶点，保留相关评分（Relevance Score）≥25的疾病相关靶点。

1.3 蛋白质-蛋白质相互作用网络构建利用半编程性质的脚本性软件R project（×64 3.6.1）对缓衰方活性成分所预测的靶点和CKD、CVD相关的疾病靶点进行映射，获得3个靶点数据的交集作为缓衰方对CKD心肾保护作用的关键靶点进行机制探讨。

使用Bisogenet数据库（http：//）构建关键靶点的蛋白质-蛋白质相互作用（Protein Protein Interaction，PPI）网络，设定物种为人类（Homo Sapiens），基于Cytoscape（Version 3.7.2）软件的
cytNCA插件计算PPI网路中各节点的度中心性（Degree Centrality，DC）、中介中心性（Betweenness Centrality，BC）和紧密中心性（Closeness Centrality，CC）等信息并进行拓扑分析，选取DC、BC和CC值均高于均值的靶点作为核心靶点。

基于Pubchem数据库（https：//）查找缓衰方药物活性成分的3D结构，并进行能量最小化处理，利用RCSB PDB数据库（http：//）查找核心靶点的相关蛋白晶体结构，利用PyMOL（Version 1.7.2）软件对靶点蛋白进行去水、加氢处理，并以对心肾血管具有保护作用的厄贝沙坦（Irbesartan）作为阳性对照，基于autodock（Version 1.5.6）软件对缓衰方药物活性成分和靶点相关蛋白进行分子模拟对接，并计算最低结合能（Binding Affinity），最低结合能≤-5 kcal/mol则说明药物活性成分和靶点蛋白对接较好，结合能的绝对值越大说明活性成分与靶点蛋白结合越好[12]。

1.4 网络分析基于bioconductor数据库（http：//）對缓衰方心肾同治作用的关键靶点进行基因本体（Gene Ontology，GO）生物学过程富集分析和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集分析，保留
P<0.05，Q<0.05的富集结果，获得关键靶点参与的生物过程（Biological Process，BP）、细胞组分（Cellular Component，CC）、分子功能（Molecular Function，MF）和作用信号通路。

基于Cytoscape软件，构建“成分-靶点-通路”网络，对缓衰方活性成分、靶点和疾病之间的相互
关系进行分析。

Keywords Network pharmacology; Huan Shuai Formula; Chronic kidney disease; Cardiovascular disease; Active components; Targets; Signaling pathway; Mechanism
中图分类号：R242;R692文献标识码：Adoi：10.3969/j.issn.1673-7202.2021.20.004
慢性肾脏病（Chronic Kidney Disease，CKD）已成为全球共同关注的公共卫生问题，全世界超过13.4%的成年人患有CKD[1]。

在我国，成年人CKD的发病率约为10.8%[2]。

心血管疾病（Cardiovascular Disease，CVD）为CKD常见的并发症，超过64.5%的CKD患者患有不同程度的CVD，约为健康人群的2～4倍[3-4]。

有研究表明，随着CKD的进展，CVD的发生风险显著增加，终末期肾脏疾病（End Stage Renal Disease，ESRD）患者CVD的发生率可达80%[5]。

CKD患者血脂异常、血压控制不佳等传统CVD危险因素导致动脉硬化和心肌重构是CVD形成的基本病变，临床主要表现为冠状动脉粥样硬化性心脏病、心律失常、慢性心力衰竭等[6-8]。

CKD肾素-血管紧张素-醛固酮系统（Renin-angiotensin-aldosterone System，RAAS）激活、蛋白尿、肾小球滤过率下降、贫血、全身营养不良等非传统CKD危险因素，可加速CVD进展;而心输出量的减少和静脉压的增加，使肾脏灌注减少，导致肾功能丢失加速。

CKD 和CVD的相互影响使伴有CVD的CKD患者死亡率高达50%，为健康人群的10～30倍，是CKD死亡的首要因素[9]。

由此可见，CKD患者防治CVD至关重要。

缓衰方以黄芪、丹参、当归、大黄、白花蛇舌草为核心组成，具有益气活血、解毒降浊之效。

临床和实验研究发现，缓衰方能够有效改善CKD患者的生命质量，延缓CKD患者进入终
点事件（透析、全因死亡、血清肌酐双倍增加）的时间[10-11]。

在长期临床实践中我们发现，缓衰方可显著减轻CKD患者胸闷、胸痛、心悸等CVD症状。

然而目前关于缓衰方对CKD心肾保护作用的有效活性成分、作用靶点和作用机制尚未完全阐明。

因此，本研究拟采用网络药理学的方法从分子水平方面研究和分析缓衰方对CKD患者心肾保护作用的有效活性成分、潜在靶点和作用机制，以期为缓衰方的应用和开发提供科学依据。

1 资料与方法
1.1 缓衰方活性成分及靶点的筛选通过检索中药系统药理学数据库与分析平台（TCMSP，http：///tcmsp.php），收集缓衰方中“黄芪”“丹参”“当归”“大黄”“白花蛇舌草”5味中药的活性成分，根据活性成分的药物毒代动力学参数进行筛选，以类药性（Drug-like，DL）≥0.18和口服生物利用度（Oral Bioavailability，OB）≥30%的化合物作为缓衰方的有效活性成分。

同时根据TCMSP数据库整理有效活性成分所涉及的靶点，并通过Uniprot数据库（https：//）进行靶点信息比对和靶点名校正。

1.2 CKD和CVD疾病相关靶点的检索 GEO数据库是全球最大的公共基因表达数据库，具有较全面的高通量分子数据。

以关键词“chronic kidney disease”在GEO数据库中检索与CKD相关的靶点，选择物种为人类（Homo Sapiens），对得到的基因靶点进行整理并找出差异基因靶点。

同时分别在GeneCards数据库（https：//）和人类孟德尔遗传数据库（OMIM，https：//）检索与CKD和CVD相关的疾病靶点，保留相关评分（Relevance Score）≥25的疾病相关靶点。

1.3 蛋白质-蛋白质相互作用网络构建利用半编程性质的脚本性软件R project（×64 3.6.1）对缓衰方活性成分所预测的靶点和CKD、CVD相关的疾病靶点進行映射，获得3个靶点数据的交集作为缓衰方对CKD心肾保护作用的关键靶点进行机制探讨。

使用Bisogenet数据库（http：//）构建关键靶点的蛋白质-蛋白质相互作用（Protein Protein Interaction，PPI）网络，设定物种为人类（Homo Sapiens），基于Cytoscape（Version 3.7.2）软件的cytNCA插件计算PPI网路中各节点的度中心性（Degree Centrality，DC）、中介中心性（Betweenness Centrality，BC）和紧密中心性（Closeness Centrality，CC）等信息并进行拓扑分析，选取DC、BC和CC值均高于均值的靶点作为核心靶点。

基于Pubchem数据库（https：//）查找缓衰方药物活性成分的3D结构，并进行能量最小化处理，利用RCSB PDB数据库（http：//）查找核心靶点的相关蛋白晶体结构，利用PyMOL（Version 1.7.2）软件对靶点蛋白进行去水、加氢处理，并以对心肾血管具有保护作用的厄贝沙坦（Irbesartan）作为阳性对照，基于autodock（Version 1.5.6）软件对缓衰方药物活性成分和靶点相关蛋白进行分子模拟对接，并计算最低结合能（Binding Affinity），最低结合能≤-5 kcal/mol则说明药物活性成分和靶点蛋白对接较好，结合能的绝对值越大说明活性成分与靶点蛋白结合越好[12]。

1.4 网络分析基于bioconductor数据库（http：//）对缓衰方心肾同治作用的关键靶点进行基因本体（Gene Ontology，GO）生物学过程富集分析和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）通路富集分析，保留
P<0.05，Q<0.05的富集结果，获得关键靶点参与的生物过程（Biological Process，BP）、细胞组分（Cellular Component，CC）、分子功能（Molecular Function，MF）和作用信号通路。

基于Cytoscape软件，构建“成分-靶点-通路”网络，对缓衰方活性成分、靶点和疾病之间的相互关系进行分析。