- 摘要Abstract
- 关键词Keywords
-
1 资料与方法
-
1.1 雷公藤重要活性成分筛选及靶点基因的预测
-
1.2 构建与分析活性成分-预测靶点网络
-
1.3 IgAN相关靶点基因的收集
-
1.4 治疗IgAN靶点的预测及构建与分析药物-活性成分-靶点-疾病网络
-
1.5 相关靶点基因功能分析及通路分析
-
2 结果
-
2.1 雷公藤的重要活性成分及对应的靶点基因
-
2.2 IgA肾病相关靶点收集
-
2.3 雷公藤重要活性成分-靶点基因网络构建与分析
-
2.4 药物活性成分-靶点-疾病（IgAN）的PPI网络构建与分析
-
2.5 药物-活性成分-靶点-疾病网络的构建与分析
-
2.6 相关作用靶点基因功能分析及通路分析
-
3 讨论
- 参考文献
目的
基于网络药理学的理念，运用网络分析的技术，拟探讨雷公藤干预IgAN（immunoglobulin Anephropathy）的潜在靶点预测及其作用机制。
方法
利用TCMSP平台收集雷公藤的重要活性成分及潜在相关靶点信息，分别从 GeneCards数据库、OMIM数据库中获取与IgAN相关联的靶点；采用STRING平台构建PPI作用网络，使用Cytoscape 3.7.0对药物-活性成分-潜在靶点-疾病网络进行构建；通过DAVID数据库对相关疾病-药物靶点基因进行GO生物功能富集分析和KEGG代谢通路分析。
结果
检索出雷公藤成分共144个，筛选获得重要活性成分共29个；收集到潜在作用靶点基因122个。检索出IgAN相关靶点333个，基因功能GO富集分析及KEGG代谢通路富集结果显示，雷公藤中的9个有效成分干预26个疾病靶点，通过对20个生物过程、12个分子功能、13个细胞成分产生作用。得到通路共50条，通过可视化分析，取靠前的30条通路，其中6条通路可能与IgAN发生及发展密切相关。
结论
该研究基于网络药理学方法，表明了雷公藤中药多成分、多靶点、多功能、多通路干预IgAN的作用，为后续临床及研究开展相应工作提供了较为可观的信息和理论依据。
Objective
Based on the concept of network pharmacology, using network analysis techniques, it is proposed to explore the potential target prediction of Tripterygium wilfordii IgAN and its mechanism of action.
METHODS
The important active components and potential related target information of Tripterygium wilfordii were collected by TCMSP platform. The targets associated with IgAN were obtained from GeneCards database and OMIM database respectively. The PPI network was constructed by STRING platform, using Cytoscape 3.7.0. The drug-active ingredient-potential target-disease network was constructed, and the GO-function enrichment analysis and KEGG metabolic pathway analysis were performed on the relevant disease-drug target genes by the DAVID database.
RESULTS
A total of 144 components of Tripterygium wilfordii were retrieved, and 29 important active ingredients were screened. 122 potential target genes were collected. 333 IgAN-related targets were retrieved, and gene function GO enrichment analysis and KEGG metabolic pathway enrichment showed that 9 active components in Tripterygium wilfordii interfered with 26 disease targets through 20 biological processes and 12 molecular functions. 13 cellular components work. A total of 50 pathways were obtained. Through visual analysis, the top 30 pathways were taken, and 6 of them may be closely related to the occurrence and development of IgAN.
Conclusion
This study based on the network pharmacology method, clarified the role of multi-component, multi-target, multi-functional and multi-channel intervention of IgAN in Tripterygium wilfordii Hook. It provided a considerable information and theoretical basis for the subsequent clinical and research work.
1 资料与方法
1.1 雷公藤重要活性成分筛选及靶点基因的预测
使用TCMSP平台（http://lsp.nwu.edu.cn/tcmsp.php）收集雷公藤的主要活性成分，对成分药物动力学(absorption and distribution and metabolism and excretion，ADME)参数进行设置以筛选出在人体内生物活性较高的成分，筛选条件为：药物相似性（drug likeness，DL） > 0.18，口服生物利用度（Oral bioavailability，OB） > 30，氢键给体（hydrogen bond donors，Hdon） ≤ 5，氢键受体（hydrogen bond acceptors，Hacc） ≤ 10。根据筛选的活性成分，从TCMSP平台的“Targets Infomation”中获取该成分作用的靶蛋白名（target name = protein names），剔除其中没有对应靶点蛋白的活性成分，利用UniProt（https://www.uniprot.org/）数据库对TCMSP平台中采集到的靶点信息进行规范化处理，然后将每种蛋白名称输入其中，得到对应的靶点基因名称。
1.2 构建与分析活性成分-预测靶点网络
利用网络可视化软件Cytoscape 3.7.0构建雷公藤的活性成分-作用靶点网络，并通过network analyzer插件对该网络的拓扑学特征进行分析。在网络中，活性成分与靶点基因由节点（node）来表示，节点之间相连的为边（edge）。基于网络中节点拓扑结构特征值“自由度(degree)”与最大自由度筛选核心网络，筛选出雷公藤中发挥治疗作用的重要活性成分，其中自由度越大，则说明介入的生物功能越多，生物学重要性权重越大。
1.3 IgAN相关靶点基因的收集
以“IgA nephropathy”为关键词，分别检索GeneCards数据库（https://www.genecards.org/）、OMIM数据库（http://omim.org/）获得IgAN疾病有关的靶点基因信息，通过相关性得分（Relevance score） ≥ 4筛选与IgAN关联性较高的靶点基因，将两个数据库获取的靶点基因进行合并。
1.4 治疗IgAN靶点的预测及构建与分析药物-活性成分-靶点-疾病网络
选用Rx64 3.5.3软件中的venn Diagram插件将雷公藤预测的相关靶点映射到IgAN相关靶点中，重叠的基因部分即可以预测为雷公藤成分治疗IgAN的作用靶点，通过STRING对雷公藤治疗IgAN相关靶点进行PPI网络分析。为了进一步可视化呈现雷公藤活性化合物使用Cytoscape3.7.0软件构建雷公藤-成分-靶点-疾病网络。借用network analyzer插件计算该网络拓扑参数，包括自由度、介数。例如靶点的自由度越大，表明与其相关的生物功能越多，所以生物学重要性占比更大；介数值与节点在网络中的重要性呈正相关。
1.5 相关靶点基因功能分析及通路分析
通过DAVID数据库（https://david.ncifcrf.gov/）对“1.4”结果中的预测靶点予以基因本体（GO）生物过程分析和京都基因与基因组百科全书（KEGG）代谢通路分析。使用GraphPad Prism 8.0软件对生物过程及利用OmicShare数据库（http://www.omicshare.com）对代谢通路均进行可视化。
2 结果
2.1 雷公藤的重要活性成分及对应的靶点基因
通过在TCMSP平台检索中输入雷公藤，共检索到已报导的雷公藤有效成分共144个。依据设置的筛选条件，共获得预测的活性成分42个，见表1。其中剔除没有对应靶点蛋白14个，剩余预测活性成分28个，并对各预测的活性成分进行编号（MOL ID）。通过TCMSP数据库中的靶点预测模型获得的上述活性成分的作用靶点蛋白，利用UniProt数据库得到对应的靶点基因名称，并获得靶点-成分关系对。
表1 雷公藤的预测活性成分
Mol ID Molecule Name MW Hdon Hacc OB (%) DL MOL000211 Mairin （丁子香萜） 456.78 2 3 55.38 0.78 MOL000296 Hederagenin （常春藤皂苷） 414.79 1 1 36.91 0.75 MOL000358 Beta-sitosterol （β-谷甾醇） 414.79 1 1 36.91 0.75 MOL000422 Kaempferol （山奈酚） 286.25 4 6 41.88 0.24 MOL000449 Stigmasterol（豆甾醇） 412.77 1 1 43.83 0.76 MOL002058 40957-99-1 388.45 2 7 57.2 0.62 MOL003182 (+)-Medioresinoldi-O-beta-D- glucopyranoside_qt（（+）-树脂醇二-O-β-D-吡喃葡萄糖苷） 388.45 2 7 60.69 0.62 MOL003184 81827-74-9 342.47 1 4 45.42 0.53 MOL003185 (1R,4aR,10aS)-5-hydroxy-1-(hydroxymethyl)-7-isopropyl-8-methoxy-1,4a-dimethyl-4,9,10,10a-tetrahydro-3H-phenanthren-2-one（（1R，4AR，10AS）-5-羟基-1-（羟甲基）-7-异丙基-8-甲氧基-1,4a-二甲基-4,9,10,10a-四氢-3H-菲仑-2-酮） 346.51 2 4 48.84 0.38 MOL003187 Triptolide（雷公藤内酯醇） 360.44 1 6 51.29 0.68 MOL003196 Tryptophenolide（胰托酚酸盐） 312.44 1 3 48.5 0.44 MOL003199 5,8-Dihydroxy-7-(4-hydroxy-5-methyl-coumarin-3)-coumarin（5,8-二羟基-7-（4-羟基-5-甲基-香豆素-3）-香豆素） 352.31 3 7 61.85 0.54 MOL003208 Celafurine（青蒿素） 369.51 2 6 72.94 0.44 MOL003209 Celallocinnine（青蒿素） 405.59 2 5 83.47 0.59 MOL003217 Isoxanthohumol（异蛇麻素） 354.43 2 5 56.81 0.39 MOL003224 Tripdiotolnide（雷公藤甲素） 360.44 2 6 56.4 0.67 MOL003225 Hypodiolide A（山海棠二萜内酯 A） 318.5 1 3 76.13 0.49 MOL003229 Triptinin B（雷公藤素B） 314.46 2 3 34.73 0.32 MOL003231 Triptoditerpenic acid B（雷藤二萜酸 B） 328.49 1 3 40.02 0.36 MOL003245 Triptonoditerpenic acid（雷藤二萜酸） 344.49 2 4 42.56 0.39 MOL003248 Triptonoterpene（雷公藤萜烯） 300.48 1 2 48.57 0.28 MOL003266 21-Hydroxy-30-norhopan-22-one 428.77 1 2 34.11 0.77 MOL003280 TRIPTONOLIDE （雷公藤内酯） 326.42 1 4 49.51 0.49 MOL003283 (2R,3R,4S)-4-(4-hydroxy-3-methoxy-phenyl)-7-methoxy-2,3-dimethylol-tetralin-6-ol
（（2R，3R，4S）-4-（4-羟基-3-甲氧基-苯基）-7-甲氧基-2,3-二甲基-四氢萘-6-醇）
360.44 4 6 66.51 0.39 MOL004443 Zhebeiresinol（浙贝素） 280.3 1 6 58.72 0.19 MOL005828 Nobiletin（川陈皮素） 402.43 0 8 61.67 0.52 MOL007415 [(2S)-2-[[(2S)-2-(benzoylamino)-3-phenylpropanoyl]amino]-3-phenylpropyl] acetate
（[（2s）-2-[（2s）-2-（苯甲酰氨基）-3-苯丙酰基]氨基]-3-苯丙基]乙酸盐）
444.57 2 6 58.02 0.52 MOL007535 (5S,8S,9S,10R,13R,14S,17R)-17-[(1R,4R)-4-ethyl-1,5-dimethylhexyl]-10,13-dimethyl-2,4,5,7,8,9,11,12,14,15,16,17-dodecahydro-1H-cyclopenta[a]phenanthrene-3,6-dione（（5S，8S，9S，10R，13R，14S，17R）-17-[（1R，4R）-4-乙基-1,5-二甲基己基]-10,13-二甲基-2,4,5,7,8,9,11,12,14,15,16,17-十二氢-1H-环戊烷[A]菲-3,6-二酮） 428.77 0 2 33.12 0.79 MOL009386 3,3’-bis-(3,4-dihydro-4-hydroxy-6-methoxy)-2H-1-benzopyran（3,3’-双-（3,4-二氢-4-羟基-6-甲氧基）-2H-1-苯并吡喃） 358.42 2 6 52.11 0.54 2.2 IgA肾病相关靶点收集
从GeneCards数据库中获取330条与IgAN相关的靶点信息、OMIM数据库中获取10条IgAN相关的靶点信息，合并2个数据库数据并且运用 EXCEL剔除重复靶点基因，结果共检索出333个相关基因信息，见表2。
表2 IgAN疾病相关靶点
数据库 基因 GeneCards数据库 CD79A、ACE、CD40LG、REN、AGTR1、ALB、TGFB1、B2M、FCAR、SPRY2、IGAN1、IGAN2、IL6、NPHS1、CCL2、SERPINB7、NPHS2、UMOD、CD40、HPRT1、IL1RN、TNF、COL4A1、NAGLU、MBL2、AGT、C1GALT1、HLA-DQA1、MAFB、VEGFA、SCGB1A1、CFP、MASP1、IL10、HFE、CRP、TNFRSF13B、EPO、IL4、NLRP3、SAA1、MME、TP53、CD19、WAS、MVK、IFNG、ADIPOQ、AICDA、HP、INS、IKBKG、CR2、TF、NEU1、MIR21、IL1B、SELE、CCND1、BTK、MPO、IL2、EDN1、TNFRSF1A、TNFSF13B、CXCL8、ICAM1、PPARG、FASLG、CASP10、PTH、AMBP、CCR5、BLM、NOS3、HLA-DQB1、WT1、IL2RA、LRP2、MUC1、HRAS、TNFRSF11B、NRAS、LIG4、HELLS、IGHE、MMP9、IL5、KRAS、CAV1、PRKCD、PMM2、VCAM1、PIGR、JCHAIN、COL4A4、HLA-DRB1、TFRC、APOA1、SERPINE1、SPP1、IGF1、TTR、RECQL4、TLR9、COL4A3、PDGFD、FN1、SOD2、COL4A5、CCN2、PVT1、ICOS、AKR1B1、CCL5、MTHFR、AGPAT2、BSCL2、WIPF1、AIRE、APOE、COL7A1、C3、COL4A2、IGHA1、LTA、FCAMR、FAS、IL1R1、FGFR3、PLA2R1、UBD、FGF2、CD2AP、H19、MEFV、PON1、SELP、GAS6、MTR、CCR1、STAT3、LAP3、TPO、IL18、PGM3、TNFSF13、FGA、NPPA、TGM2、PIK3R1、HAVCR1、AGER、ITGA4、CFH、CST3、MTOR、CLU、LCN2、CD36、CFHR5、PTPN11、MIR25、VDR、VWF、MIR93、SLC2A1、BRAF、HGF、CDK4、IDH1、IDH2、DMD、MIR106B、MIR181A2、MIR32、MIR19A、MIR15A、IFNA1、C4A、MYH9、AGTR2、CAVIN1、CASP8、SYNPO、VTN、CCL3、CASP1、CSF2、ADA、FABP1、AXL、TNFSF11、APOL1、SELL、CXCL12、IL6R、LYZ、MALAT1、CD44、SRC、FCGRT、CTLA4、CASP3、CASR、SERPINC1、TLR4、HLA-A、IL15、MIRLET7B、SOD1、HLA-B、HMOX1、ST6GALNAC2、FOS、LEPQTL1、CYP11B2、RBP4、XBP1、HNP1、HNF4A、APOH、ITGA5、ADM、TG、CXCL10、MYC、ACTN4、CDKN2A、BMP6、PSMB8、MBP、LEP、NOD2、PDGFA、THBD、NAGA、PARP1、IL1A、TNFRSF17、CXCR4、BMP7、CD27、IRF5、HLA-DPB1、TLR2、PTPRC、MYOM2、C1GALT1C1、PSTPIP1、MMD、CDC42、IL12A、SMAD2、LTF、SMAD3、SNX9、CYBA、NOS2、WRN、BCL2L11、RECQL、RECQL5、ENPP1、FLT1、TNFRSF1B、F2、MIR146A、MIF、SGK1、KCNJ11、SPG7、PCNA、LPA、IL13、GSN、PTX3、ACTR3、TIMP1、ACE2、CAT、GRB2、AKT1、CFHR1、CXCL9、HMGCR、MAPK1、HNF1A、INPP5E、EGF、CD38、HSPG2、RNASE3、ACTR2、HLA-DPA1、NCKIPSD、PLG、KL、JAK2、CASP9、RASGRP1、DNMT1、FYN、NCK1、CTTN、WASF1、TFE3、WASF2、WASL、SPN、WASF3、TRIP10、
PES1、SCHIP1、WASH3P、WASH6P、WASF4P、WASF5P、TGFB2、AQP2、PRDM1、ABCA1、CYP3A4、LPL、NSD2、ITGAM、IGHA2、PDGFB、ABCB1
OMIM数据库 FCAMR、HNFJ3、IGAN2、MBNP、IGAN1、HNP1、NEDE、CD79A、FCAR、IGA 2.3 雷公藤重要活性成分-靶点基因网络构建与分析
雷公藤活性成分-靶点基因网络由150个节点组成，由28个雷公藤成分节点和122个靶点基因节点，图中28个绿色三角形节点表示雷公藤活性成分节点，图中红色椭圆形122节点表示靶点基因节点，节点间的边表示活性成分与靶点基因之间的相互作用关系，最大自由度值51，最小自由度值1，平均自由度4.973，结果见图1。根据最大自由度与平均自由度之间的值筛选核心成分-靶点基因网络，结果见图2，筛选出的核心成分按最大自由度由大到小筛选分别为：山奈酚（51）、雷公藤内酯醇（31）、川陈皮素（27）、豆甾醇（26）、β-谷甾醇（26）、雷藤二萜酸 B（21）、雷公藤素B（18）、雷公藤萜烯（17）、胰托酚酸盐（17）、常春藤皂苷（16）、81827-74-9（16）、（2R，3R，4S）-4-（4-羟基-3-甲氧基-苯基）-7-甲氧基-2,3-二甲基-四氢萘-6-醇（14）、（1R，4AR，10AS）-5-羟基-1-（羟甲基）-7-异丙基-8-甲氧基-1,4a-二甲基-4,9,10,10a-四氢-3H-菲仑-2-酮（14）、雷公藤内酯（13）、异黄花蒿醇（12）、5,8-二羟基-7-（4-羟基-5-甲基-香豆素-3）-香豆素（12）、雷藤二萜酸（10）、浙贝素（6）、40957-99-1（6）、3,3 -双-（3,4-二氢-4-羟基-6-甲氧基）-2H-1-苯并吡喃（5）。
2.4 药物活性成分-靶点-疾病（IgAN）的PPI网络构建与分析
通过Rx64 3.5.3软件映射出药物与疾病相互作用的靶点交集，交集的靶点可能为雷公藤治疗IgAN的预测靶点，用韦恩图表示，如图3，得到共同靶点数26个。将作用的预测靶点导入String数据库，经过数据库运行获得靶点相互作用网络图（PPI），如图4，图中有26个节点，200条边，平均节点度为15.4，平均局部聚类系数为0.773。图示说明雷公藤的活性成分作用于多个靶点基因，且各个靶点之间相互作用、协同起效。
2.5 药物-活性成分-靶点-疾病网络的构建与分析
借用Cytoscape 3.7.0软件对该网络进行分析和可视化处理，利用软件中的Network Analyzer插件分析其节点的网络拓扑学特征，根据其自由度（Degree）、介数（Betweenness Centrality）等网络拓扑特征值。筛选网络中的核心靶点基因。“药物-活性成分-靶点-疾病”网络由37个节点组成，由9个雷公藤活性成分节点、26个靶点节点、1个疾病节点和1个药物节点组成，黄色椭圆形代表活性成分节点，紫色菱形代表蛋白靶点节点。这里的26个靶点节点既是雷公藤成分作用靶点，亦是IgAN疾病发生相互作用的靶点。见图5。网络分析表明，雷公藤内酯醇自由度为12，介数为0.110211，预测为雷公藤治疗IgAN中发挥的重要成分，其次为山奈酚（自由度为11、介数为0.08184897）、川陈皮素（自由度为7、介数为0.03074641）等。PPARG为活性成分中7个活性成分的共有目标靶点，因此可以认为PPARG是雷公藤治疗IgAN的最重要靶点，其余NOS2、CASP9、CASP3等亦为相对重要靶点，见表3。
表3 雷公藤作用IgAN重要活性成分及相关作用靶点自由度和介数
MOL ID 自由度 介数 靶点 自由度 介数 MOL003187 12 0.11021095 PPARG 7 0.05139536 MOL000422 11 0.08184897 NOS2 5 0.05139536 MOL005828 7 0.03074641 CASP3 4 0.02087731 MOL000358 5 0.02333455 CASP9 3 0.0137202 MOL003199 3 0.00386054 AKR1B1 2 0.01986017 MOL000296 2 0.00541572 AKT1 2 0.00264871 MOL003217 2 0.00153345 C3 2 0.00200734 MOL000449 2 0.00541572 CASP8 2 0.00615646 MOL003283 2 0.00153345 CD40 2 0.00200734 靶点 自由度 介数 靶点 自由度 介数 CXCL8 2 0.00200734 MMP9 2 0.00532571 CXCR4 2 0.00200734 PON1 2 0.00615646 CYP3A4 2 0.00264871 SELE 2 0.00264871 FOS 2 0.00200734 STAT3 2 0.00200734 HMOX1 2 0.00264871 TIMP1 2 0.00532571 ICAM1 2 0.00264871 VCAM1 2 0.00264871 IFNG 2 0.00200734 VEGFA 2 0.00200734 IL2 2 0.00200734 LYZ 2 0.01986017 IL4 2 0.00200734 2.6 相关作用靶点基因功能分析及通路分析
通过David数据库对活性成分-疾病相互作用的26个关键靶点分别从生物过程（Bioprocess）、分子功能（Molecular Function）、细胞成分（Cell Component）3个水平进行分析，设置P < 0.05，筛选出具有统计学意义的生物学过程，进行生物学过程分析，使用GraphPad Prism 8.0.1软件将这些条目做成可视化条形图，如图6。对GO富集的生物学过程进行分析，预测靶点在生物过程（BP）中一共获得20个条目，在分子功能（MF）中一共获得12个条目，细胞成分（CC）中一共获得13个条目。
对雷公藤作用于IgAN的26个关键靶点进行KEGG代谢通路分析，设置P < 0.05，筛选出具有统计学意义的代谢通路，得到通路共50条，借用OmicShare软件进行可视化分析（通路条目过多，取靠前的30条通路，其中6条通路与IgAN发生及发展相关，如图7。
3 讨论
目前IgAN的发病机制尚不明确，有研
究[8] 表明IgAN发病过程主要与自身免疫失调、机体感染（特别是上呼吸道感染）或遗传易感性等方面密切相关，并且IgA1分子异常糖基化[9] 和免疫复合物系膜沉积导致的补体激活，激发免疫炎症[10] 反应加剧了IgAN的临床表现。近年来，多项研究证实了雷公藤的多种成分具有抗氧化、抗炎、免疫抑制及抗增殖等作用[11] ，但雷公藤成分改善IgAN理化指标的潜在机制研究较少或不足以系统阐述生物过程及代谢通路。本研究通过“药物-活性成分-作用靶点-疾病”网络的构建、拓扑学参数分析以及基因功能和通路的信息开发对雷公藤治疗IgAN的活性成分、潜在靶点、生物功能以及分子机制进行了剖析。本研究结果表明雷公藤中环氧二萜内酯（雷公藤内酯醇）、黄酮类（山奈酚、川陈皮素、异蛇麻素）、甾醇类（β-谷甾醇、谷甾醇）、香豆素衍生物（（5,8-二羟基-7-（4-羟基-5-甲基-香豆素-3）-香豆素）等化合物显示了治疗IgAN的潜在活性。同时研究中剖析出的122个预测靶标中，26个是多个活性化合物的共同靶标，其中自由度最高的靶标为PPARG、NOS2、CASP3、CASP9、AKR1B，同时表明它们可能是雷公藤干预IgAN的关键靶点。通过String软件可视化显示IgAN靶点蛋白与化合物有着直接或间接作用的关系，并且大部分成分均能直接作用于PPARG这个靶点上而发挥关键治疗作用，同时也与其他主要靶点相互作用共同干预IgAN进展。过氧化物酶体增殖物激活受体-γ（PPARG）是一种已被证明是配体激活且同时可调节自噬信号的受体，而自噬本身就是一种抗炎机制，而晚期糖基化终末产物AGEs可以通过增加炎症介质水平从而下调PPAR
G[12] ，研究[13] 显示AGEs受体（RAGE）的表达引发IL-1，IL-6，TNF-α和MMP-13的表达,炎症细胞因子似乎突出下调PPARG。本研究通过KEGG通路分析及靶点基因GO功能分析以进一步探讨雷公藤治疗IgAN的潜在通路及作用机制，结果显示这些靶点涉及的生物过程（BP）、分子功能（MF）、细胞成分（CC）相关信息。如图6所示，P值（PValue） < 0.05的20个生物过程中，排名靠前如凋亡过程的负调节、对脂多糖的反应、炎症反应、对缺氧的反应、蛋白磷酸化的正调控、免疫反应的调节、免疫应答等过程，可通过减少氧化应激反应、炎症反应及调控免疫反应的产生与进展等方式，最终影响IgAN疾病的发生与转归；在P值（PValue） < 0.05的12个分子功能中，其中IgAN发病相关的功能有细胞因子活性、半胱氨酸型内肽酶活性参与细胞凋亡过程、死亡受体结合、转录调控区DNA结合等，干预方式主要涉及IgAN病程中的炎症反应、细胞免疫调节等因素；在P值（PValue） < 0.05的13个细胞成分中，与IgAN病程密切相关的细胞成分有细胞外空间、质膜外侧、膜筏、细胞表面、细胞外外泌体、细胞质的核周区域、死亡诱导信号复合体、细胞质、细胞内膜结合细胞器、RNA聚合酶II转录因子复合物，这些细胞成分可能参与了炎症反应，氧化应激及体液免疫等过程，同时表明雷公藤可能是通过作用在这些生物过程从而干预IgAN以达到治疗的目的；同时对雷公藤干预IgAN相关的26个靶点进行KEGG代谢通路分析，结果如图7所示。在可视化的30条通路中，与IgAN的病程密切相关的可能是TNF信号通路、HIT-1信号通路、T细胞受体信号通路、Toll样受体信号通路、Jak-STAT信号通路及NF-κB信号通路这6条，其中雷公藤作用于TNF信号通路干预IgAN病程的显著性最大，可能是因为该通路是雷公藤干预IgAN疾病过程中的重要通路，该通路可介入各种信号的传导，主要是通过调节细胞凋亡或控制炎症反应以释放IL-4、IL-2等促炎因子和NO等炎症介质，造成炎症反应与激活更多的炎症基因，造成IgAN的补体系统的激活，最后通过旁路途径和凝集素途径导致对IgAN的肾损
伤[14,15,16,17] 。同时本研究显示HIF-1信号通路的表达与IgAN发生及发展密切相关，其很有可能就是雷公藤通过干预HIF-1信号通路的表达实现抗氧化作用。Tanaka T等人的研
究[18] 表明肾小管间质缺血或缺氧被认为是肾纤维化发生的关键微环境因子。介导应答的分子机制ES对缺氧的作用是复杂的，其中缺氧诱导因子(HIF)起着关键作用，HIF表达与细胞外基质增多呈正相关，从而加快了肾间质纤维化的进程[19] 。JAK-STAT通路是细胞间信号传导的主要通道之一，也是代表细胞因子信息内传导的重要通路之一[20] 。蔡婷等人[21] 通过构建IgAN模型肾组织，发现JAK2、STAT3的表达显著上调，尤以肾小管间质中为甚，在使用雷公藤后，二者的表达明显下降，因此猜测JAK-STAT信号通路极大概率参与IgAN肾组织损伤，以肾小管间质损伤为主，雷公藤通过干预JAK-STAT信号通路对IgAN的肾小管间质进行保护。陈明喆[22] 通过降低JAK-STAT信号通路中IgAN大鼠模型肾脏组织的蛋白STAT3表达，发现IgAN大鼠肾功能具有显著的改善作用，可能是由于降低蛋白STAT3表达，使机体抗氧化能力得以提高。Hye-Ju Han等[23] 构建TCRδ-/- 小鼠实验发现T细胞可以控制体液免疫反应，T细胞具有重要的调节功能，尤其是对维持正常黏膜和全身IgA反应有重要作用。Toll样受体是直接参与与炎症反应相关的慢性肾病的单个跨膜非催化性蛋白质，在人类中，Toll样受体在多种细胞类型上表达，包括肾小球系膜细胞和肾小管上皮细胞，TLR4已被证实与系膜细胞损伤和肾纤维化有关[24,25,26] 。哺乳动物核因子NF-κB信号通路是细胞内的基本转录因子系统，是针对各种诱导的细胞外刺激来源[27] 。NF-κB的指标激活是二聚体Rel蛋白8的核转位，它调节许多涉及NF-κB依赖性基因，包括炎症，免疫，细胞凋亡，细胞增殖和分化。有研究[28] 显示miＲ-21的低表达可能会使p65磷酸化达到激活NF-κB通路以此增强IgAN病理的炎性反应。综上所述，雷公藤可能通过干预诱导炎症发生的关键靶点基因的转录和表达、阻断促肾间质损伤信号的传导、减少氧化应激反应的产生及调控免疫反应等方面实现改善IgAN发生发展过程的作用。本研究通过网络拓扑学分析和大数据挖掘的联合使用，建立了“雷公藤-活性成分-靶点-IgAN”关系网络，紧接着使用DAVID数据库进行基因功能（GO）富集分析及KEGG代谢通路富集分析，最后预测了雷公藤干预IgAN的关键靶点及其相关作用机制，这将为雷公藤治疗IgAN的药效活性成分、探讨相关谱效、作用机制原理及生物评价相应指标的筛选等探索提供理论依据，为雷公藤治疗IgAN的机制研究提供了科学依据与实用方法，为接下来进一步的机制探讨提供了理论基础。
参考文献
2
- 1
Ouyang Y, Zhu L, Shi M, et al.A Rare Genetic Defect of Increased the Risk for Progression of IgA Nephropathy. Frontiers in immunology, 2019, 10:537. DOI:10.3389/fimmu.2019.00537.
- 2
Knoppova B,Reily C,Maillard N,et al.The Origin and Activities of IgA1-Containing Immune Complexes in IgAN ephropathy. Frontier in immunology, 2016, 7:117. DOI:10.3389/fimmu.2016.00117.
- 3
朱娇丽. IgA肾病的发病机制及治疗进展. 临床医药文献电子杂志, 2017, 4(94): 18430.
- 4
卢晓梅, 唐雪莲, 秦登优, 等. 雷公藤多苷联合RAS阻断剂治疗CKD2-3期IgA肾病. 实用医学杂志, 2016, 32(01): 137-139.
- 5
朱彩凤, 朱斌, 包自阳, 等. 雷公藤内酯醇对IgA肾病大鼠蛋白尿和nephrin、podocin蛋白及mRNA表达的影响. 中国中西医结合肾病杂志, 2015, 16(02): 106-109+1.
- 6
李文文, 刘晓丽, 吴颢, 等. 雷公藤治疗IgA肾病疗效和安全性的Meta分析. 中国循证医学杂志, 2015, 15(02): 206-214.
- 7
汪力, 樊均明, 王少清, 等. 雷公藤治疗IgA肾病的系统评价. 实用药物与临床, 2015, 18(03): 293-297.
- 8
Suzuki Y, Suzuki H, Sato D, et al. Reevaluation of the mucosa-bone marrow axis in IgA nephropathy with animal models . Advances in oto-rhino-laryngology, 2011, 72: 64-67.
- 9
Di Lernia V. IgA nephropathy during treatment with TNF-alpha blockers: Could it be predicted?. Medical hypotheses, 2017, 107: 12-13.
- 10
汪力, 王少清. IgA肾病的发病机制及治疗进展. 医学综述, 2015, 21(15): 2712-2714.
- 11
朱林波, 王娟, 刘林林, 等. 雷公藤多甙联合血管紧张素受体拮抗剂治疗中等量蛋白尿IgA肾病疗效研究. 中国实用内科杂志, 2017, 37(06): 549-551.
- 12
ChenC, MaC, Zhang Y, et al.Pioglitazone inhibits advanced glycation end product-induced TNF-α and MMP-13 expression via the antagonism of NF-κB activation in chondrocytes. Pharmacology, 2014, 94(5-6): 265-272.
- 13
Wang ZJ, Zhang HB, Chen C, et al. Effect of PPARG on AGEs-induced AKT/MTOR signaling-associated human chondrocytes autophagy. Cell biology International, 2018, 42(7): 841-848.
- 14
Lee CH, Kim HJ, Lee YS, et al. Hypothalamic Macrophage Inducible Nitric Oxide Synthase Mediates Obesity-Associated Hypothalamic Inflammation. Cell reports, 2018, 25(4): 934-946.e5.
- 15
Chen J, Liu GZ, Sun Q, et al.Protective effects of ginsenoside Rg3 on TNF-α-induced human nucleus pulposus cells through inhibiting NF-κB signaling pathway. Life sciences, 2019, 216: 1-9.
- 16
Maillard N, Wyatt RJ, Julian BA, et al. Current Understanding of the Role of Complement in IgA Nephropathy Journal of the American Society of Nephrology:JASN, 2015, 26(7): 1503-1512.
- 17
Onda K, Ohi H, Tamano M, et al. Hypercomplementemia in adult patients with IgA nephropathy. Journal of clinical laboratory analysis, 2007, 21(2): 77-84.
- 18
Tanaka T.The role of hypoxia, increased oxygen consumption, and hypoxia-inducible factor-1 alpha in progression of chronic kidney disease. Current opinion in nephrology and hypertension, 2010, 19(1): 43-50.
- 19
Fang Y,Yu X,Liu Y,et al.miR-29c is downregulated in renal interstitial fibrosis in humans and rats and restored by HIF-α activation. American journal of physiology. Renal physiology, 2013, 304(10): F1274-82.
- 20
戎赞华. PDGF-D和JAK/STAT信号通路在IgA肾病发病中的作用以及雷帕霉素对IgA肾病的影响. 河北: 河北医科大学, 2008.
- 21
蔡婷. JAK2/STAT3信号通路在大鼠IgA肾病模型中的表达及雷公藤对其影响. 辽宁: 辽宁医学院, 2013.
- 22
陈明喆, 李易, 何伟. 橙皮苷对IgA肾病大鼠肾脏组织中JAK/STAT信号通路相关蛋白STAT3和p-STAT3表达的影响. 解放军医药杂志, 2017, 29(03): 10-14.
- 23
Han HJ, Jang YS, Seo GY, et al. MurineγδT Cells Render B Cells Refractory to Commitment of IgA Isotype Switching. Immune network,2018, 18(4): e25.
- 24
Zou JN, Xiao J, Hu SS, et al. Toll-like Receptor 4 Signaling Pathway in the Protective Effect of Pioglitazone on Experimental Immunoglobulin A Nephropathy. Chinese medical journal, 2017, 130(8): 906-913.
- 25
Anders HJ. Toll-like receptors: emerging concepts in kidney disease. Current opinion in nephrology and hypertension, 2007, 16(3): 177-183.
- 26
Lim BJ,Lee D,Hong SW.Toll-like receptor 4 signaling is involved in IgA-stimulated mesangial cell activation. Yonsei medical journal, 2011, 52(4): 610-615.
- 27
Zhang H, Mao X, Sun Y, et al. NF-κB upregulates ubiquitin C-terminal hydrolase 1 in diseased podocytes in glomerulonephritis. Molecular medicine reports, 2015, 12(2): 2893-2901.
- 28
张莹, 龚剑, 李晓蕾, 等. 下调的miR-21通过激活NF-κB通路促进IgA肾病的研究. 吉林医学, 2016, 37(06): 1302-1304.
2
摘要
Abstract
＊ 国家自然科学基金委员会地区科学基金项目（81760807）：基于SIRT3调控能量敏感通路AMPK\PCG-1\CPT1A促进肾纤维化PTCs生成机制研究“益肾活血法”抗肾纤维化效应；负责人：钟建；广西自然科学基金委员会项目（2016JJA140633）：从周细胞-内皮细胞cross-talk调控肾纤维化机制研究CKD“久病入络”实质及益肾活血法干预效应；负责人：钟建;广西中医药大学2019年研究生教育创新计划项目立项课题（YCSY20190066）:基于数据挖掘的IgA肾病中医证候及用药特点研究；负责人：陈虹彤。
免疫球蛋白A肾病（IgAN）是一种自身免疫疾病，其特征在于半乳糖缺乏的IgA1免疫复合物在肾小球系膜区沉
雷公藤是一种卫矛科植物，其主要功效为杀虫、消炎、解毒。雷公藤相对于免疫抑制药物价格较为低廉，且由于中药多靶点、多途径的优势，该药具有的多种有效成分制剂在治疗IgAN中效果明显。雷公藤中的有效成分雷公藤内酯醇具有抗炎、免疫抑制等作用，在减少蛋白尿、抑制肾功能转归方面具有可靠的疗效，因此在免疫性肾损伤疾病方面有广泛应
为了进一步探寻雷公藤干预IgAN的作用机制，给予雷公藤制剂治疗IgAN临床应用提供证据，运用网络药理学拓扑分析技术，对中药材的活性成分及其作用的靶点基因网络、药物-疾病靶点相互作用（PPI）网络等多层次多网络的构建，探索雷公藤有效成分作用于IgAN靶点之间的关联性网络，契合中医药作用的整体性、复杂性和系统性的特点。本文借助中药系统药理学技术平台（TCMSP）结合GeneCards、OMIM疾病靶点数据库，通过构建雷公藤-活性成分-作用靶点-IgAN疾病网络相互映射来挖掘雷公藤治疗IgAN疾病的作用靶点，利用DAVID数据库进行GO生物功能富集分析和KEGG代谢通路分析，探索雷公藤干预IgAN的作用机制及代谢通路研究，为雷公藤药效物质治疗和干预IgAN发病机制的研究提供理论依据。
