Prediction of characteristic chromatogram for Abri Herba based on network pharmacology and molecular docking

赛默飞Vanquish Core型液相色谱仪（赛默飞世尔科技有限公司）；岛津LC-20A型液相色谱仪（日本岛津株式会社）；SK7200H超声波清洗器（上海科导超声仪器有限公司）；RH-600A高速多功能粉碎机（浙江荣洁工资有限公司）；XS105DU电子天平（梅特勒-托利多公司）；鼓风式干燥箱（上海般诺生物科技有限公司）；色谱柱：Thermo AcclaimTM 120 C₁₈（4.6 mm×250 mm，5 μm，120Å，赛默飞世尔科技有限公司）、Agilent ZORBAX SB-C₁₈（4.6 mm×250 mm，5 μm，安捷伦科技公司）、Diamonsil® C₁₈（4.6 mm×250 mm，5 μm，北京迪马科技有限公司）。

鸡骨草（经海军军医大学药学系黄宝康教授鉴定为豆科植物广州相思子Abrus cantoniensis Hance的干燥全草），药材市场采购或产地采集，根据采集地不同分别设置批号为JGC-xx（采集序号)-xxx（样品序号），共13个批次；相思子碱对照品（中国食品药品检定研究院，批号：111808-202003）、刺桐碱对照品（中国食品药品检定研究院，批号：112058-202001）、夏佛塔苷对照品（中国食品药品检定研究院，批号：111912-202204，纯度：94.9%）、异夏佛塔苷对照品（四川维克奇生物科技有限公司，纯度：98%）；乙腈、甲醇（色谱纯，Merck公司）；甲酸（分析级，国药集团）；蒸馏水（广州屈臣氏食品饮料有限公司）。

检索中药网络药理学常用数据库与分析平台TCMSP（https://old.tcmsp-e.com/tcmsp.php）、TCMID（http://www.megabionet.org/tcmid/）、ETCM（http://www.tcmip.cn/ETCM/index.php/Home/Index/）等，结合鸡骨草主要活性化合物文献检索，根据OB值、OD值以及类药5原则等，获取鸡骨草中活性明确且具有成药性的化合物，再通过 Pubchem查询相关化合物Canonical SMILES 编号，输入Swiss Target Prediction（http://www.swisstargetprediction.ch/）获取化合物所对应的靶标蛋白的Uniprot ID。共选取Soyasaponin Bb、Kaikasaponin Ⅲ、Abrisapogenol A、Abrisapogenol D、Sophoradiol、Kudzusapogenol A、 Abrisaponin I、Abrine、Hypaphorine、Schaftoside、Isoschaftoside、butin等12种不同结构的化合物，筛选其药理活性作用靶点，筛去重复靶标，最终获得237个靶标，通过Cytoscape3.8.2软件构建“鸡骨草-成分-潜在靶点”可视化网络图（图1）。网络中包括250个节点和440条边，平均节点度值3.520。

将“2.1.1”项下筛选所得靶标导入STRING 11.0 在线数据库（https://string-db.org/cgi/input.pl）进行PPI 网络分析。选择物种为“Homo sapiens”，蛋白交互参数评分值为“high confidence>0.9”，隐藏网络中无联系的节点，其余参数设置不变，获得化合物靶点PPI网络图（图2）。结果共获得237个节点，360条边，平均节点度为3.04，预期边数135，PPI富集P值<1.0e-16。再将PPI结果以TSV文本格式导入Cytoscape 3.8.2软件中，进行拓扑属性分析，选取亲密度、间隔度、自由度3个重要参数，共计筛选21个网络参数大于均值（0.322、9.662、5.180）的核心靶标（图3）。包括信号传导与转录激活因子3（STAT3）、基质金属蛋白酶2、9（MMP2、9）、蛋白激酶（AKT1）、表面活性蛋白C（SRC）、淋巴细胞特异蛋白酪氨酸激酶（LCK）、雌激素受体alpha（ESRα）、过氧化物酶体增生激活受体γ（PPARγ）、一氧化氮合酶2（NOS2）、缺氧诱导因子-1alpha（HIF1A）、生长因子受体结合蛋白2（GRB2）、双微体2蛋白（MDM2）、糖原合酶激酶3Β（GSK3B）、酪氨酸蛋白激酶受体B2（EPHB2）、α-突触核蛋白（SNCA）等靶标。

应用metascape 数据库（https://metascape.org/gp/index.html）对21个主要潜在的靶点进行基因本体（GO）功能和基因组百科全书（KEGG）通路富集分析。其中，靶点经生物过程（BP）相关条目463 条，主要涉及细胞凋亡正调控（positive regulation of cell death）、蛋白质水解调控（regulation of proteolysis）、细胞对有机氮化合物的反应（cellular response to organonitrogen compound）、TGF-β响应（response to transforming growth factor-beta）等；细胞组分（CC）有34个条目，主要涉及复合物转录调节（transcription regulator complex）、谷氨酸突触（glutamatergic synapse）、细胞质核周区（perinuclear region of cytoplasm）等；分子功能（MF）有42个条目，主要涉及DNA-结合转录因子结合（DNA-binding transcription factor binding）、转录因子结合（transcription factor binding）、蛋白质结构域特异性结合（protein domain specific binding）、蛋白激酶活性（protein kinase activity）、分子适配器活性（molecular adaptor activity）等，各选取显著性前10 条目展示（图4）。KEGG通路富集分析得到79个通路条目，选取显著性前20的条目通过桑基图呈现（图5）。Y轴为信号通路，X轴为该通路靶向基因占总基因的比率，气泡颜色表示该通路基因富集的显著性，气泡大小表示该通路的基因数量。其中AKT1、STAT3、HIF1A、GRB2、MMP9等主要靶点蛋白与松弛素信号通路（relaxin signaling pathway）、甲状腺激素信号通路（thyroid hormone signaling pathway）、B细胞受体信号通路（B cell receptor signaling pathway）、HIF-1信号通路（HIF-1 signaling pathway）、ErbB信号途径（ErbB signaling pathway）、PI3K-Akt信号通路（PI3K-Akt signaling pathway）、JAK-STAT 信号通路（JAK-STAT signaling pathway）、乙型肝炎（hepatitis B）等具有较强的关联性。

通过PubChem获取相思子碱、刺桐碱、夏佛塔苷等化合物3D结构，根据自由度选取“成分-靶点-通路”网络中靠前的靶标，检索PDB数据库（http://www.rcsb.org/）获得靶标蛋白结构，再将蛋白与化合物结构文件导入 pymol 软件进行去水、加氢、删除重复链等预处理，模拟化合物与靶标作用模式，计算结合能，结合能小于0表明可以自由结合，结果如表1、图6所示。

成分分析采用C₁₈色谱柱（Thermo、Agilent、Diamonsil）；流动相采用梯度洗脱，A相为乙腈，B相为0.2%甲酸溶液，梯度洗脱程序：0～10 min（5%～10%B）；10～20 min（10%～12%B）；20～30 min（12%～13%B）；30～35 min（13%～14%B）；35～40 min（14%～16%B）；40～55 min（16%～18%B）；55～60 min（18%～20%B）；60～65 min（20%～22%B）；65～70 min（22%～30%B）；70～75 min（30%～40%B）；75～80 min（40%～60%B）；80～85 min（60%～75%B）；85～95 min（75%～75%B）。柱温设置为25℃、30℃、35℃。流速：0.8 ml/min、1.0 ml/min、1.2 ml/min。检测波长278 nm。进样量：10 μl。

①对照品储备液的制备：精密称取相思子碱、刺桐碱、夏佛塔苷、异夏佛塔苷对照品5~10 mg，分别置10 ml容量瓶中，加甲醇定容，摇匀，制备对照品储备液。②供试品溶液的制备：精密称取样品粉末1 g置锥形瓶中，加甲醇50 ml，称定重量，加热回流提取2 h，取出，称定，补足损失的量，滤过，回收溶剂，残渣加甲醇5 ml使溶解，微孔滤膜过滤，制得供试品溶液。

精密量取混合对照品溶液10 μl注入液相色谱仪，测定对照品保留时间及峰面积，平行测定5次，计算相思子碱、刺桐碱、夏佛塔苷、异夏佛塔苷保留时间及峰面积RSD，保留时间RSD分别为0.6%、0.7%、0.4%、0.3%，峰面积RSD均为0.3%（图7、表2）。

取供试品JGC01-010按照2.3.2项下供试品溶液制备法制备供试品溶液，分别在0、2、4、12、24 h测定相思子碱、刺桐碱、夏佛塔苷、异夏佛塔苷对照品保留时间，计算RSD，分别为0.09%、0.10%、0.05%、0.05%（表3）。

测定13批样品，记录样品中相思子碱、刺桐碱、夏佛塔苷、异夏佛塔苷等色谱峰的保留时间及峰面积，并计算其相对保留时间及偏差（表4 ）。根据出峰时间及样品化合物峰型特征分析，部分供试品中异夏佛塔苷未检测到或含量较低超出仪器检出范围，不具有共有峰特征。因此选择相思子碱、刺桐碱、夏佛塔苷等3种化合物的色谱峰作为特征峰，相对保留时间为1.00、1.21、2.39，可控制偏差范围±5%。参考《中华人民共和国药典》（2020年版）通则指导原则的有关规定，拟定供试品色谱中应呈现相思子碱、刺桐碱、夏佛塔苷3种化合物的特征峰，并应与对照品色谱峰保留时间相对应（图8）。

分别采用Thermo、Agilent、Diamonsil色谱柱考察色谱柱的影响；柱温在20℃、25℃、30℃考察柱温的影响；流速为0.8、1.0、1.2 ml/min考察流速的影响。根据测定数据分析，色谱柱、柱温、流速均对特征峰相对保留时间具有一定的影响，其中流速及色谱柱型号对峰3、峰4影响较大，各色谱峰的影响RSD小于5.0%，属于可接受的范围（表5）。

运用网络药理学的方法对中药成分进行研究，然后预测其药理活性作用测靶点、分析通路，再通过实验验证，探索中药的作用机制^[21]，进而预测中药作用的物质基础，寻找中药质量标志物，能够较好地结合中药的药效建立质量控制标准^[22]。通过网络药理学预测发现，鸡骨草主要成分相思子碱、刺桐碱、夏佛塔苷等药理活性与STAT3、MMP2、AR、AKT1、SPC、LCK、ESR1、CTNNB1、PPARγ、NOS2、HIF1A等靶标有关，可能作用于松弛素信号通路、HIF-1信号通路、ErbB信号途径、PI3K-Akt信号通路、JAK-STAT 信号通路等多种信号通路。实验研究发现，相思子碱可以通过抑制NOS、VEGF的表达，促进IL-1β刺激的人软骨细胞C28/I2进行细胞增殖并抑制细胞凋亡，进而起到治疗骨性关节炎的作用^[23]。刺桐碱可以通过PI3K-Akt信号通路调控PPARγ起到抗炎的作用^[24]。夏佛塔苷及异夏佛塔苷等能够通过调控PPARα及其下游蛋白，减少促炎性细胞因子生成，产生抗炎、抗氧化等作用，具有潜在的非酒精性脂肪肝治疗作用^[25]。异夏佛塔苷可有效抑制脂多糖诱导的iNOS生成和促炎性细胞因子（TNF-α、IL-1β和COX2）的表达，能够显著降低脂多糖诱导的HIF1A、HK2和PFKFB3蛋白的表达^[26]。文献报道的实验研究结果进一步验证了网络药理学的预测结果。因此，网络药理学在基于药理活性的鸡骨草质量标志物预测上具有一定价值，能够为鸡骨草质量标准的建立提供药效预测及研究基础。

鸡骨草药用为全草，但市场上鸡骨草药材叶损失严重，这是否会造成药材成分及其含量的变化进而影响药材质量？该实验分别取鸡骨草（JGC-003-012）供试品根、茎、叶等不同部位粉末1 g，参照2.3.2项下制备供试品溶液，测定药材不同部位相思子碱、刺桐碱、夏佛塔苷等化合物的含量，记录保留时间及峰面积。计算含量发现根、茎、叶等不同部位所含3种成分含量分别为根：0.12%、0.08%、未检出；茎：0.02%、0.01%、0.06%；叶：未检出、0.01%、0.20%，相思子碱及刺桐碱主要存在于根、茎中，夏佛塔苷主要存在于叶和茎中，与文献报道^[27]一致。因此保证根、茎、叶的相对存量是保障鸡骨草药材质量的重要因素。另外该实验同时测定了13批次鸡骨草药材的水分、总灰分、酸不溶性灰分及浸出物的量，对鸡骨草质量标准的研究具有一定意义。

鸡骨草除药用外，岭南多地常用于食疗保健中。由于相思子属鸡骨草、毛鸡骨草等基源相近、形态相似、民间及中医临床上常同功效应用，市场上常有毛鸡骨草等混淆或掺伪为鸡骨草药材，而现行版《中华人民共和国药典》（一部）2020年版及部分收载鸡骨草的中药标准等对鸡骨草质量的控制主要集中在性状和相思子碱的薄层鉴别上，而相思子碱在相思子属相思子、毛鸡骨草、美丽相思子等植物中均有发现，仅以此为指标化合物专属性不强。目前基于化学成分与药理活性的关联为依据制定中药质量仍然是中药质量标准研究的主要方向。本实验通过网络药理学与分子对接预测鸡骨草主要活性成分，筛选出药理活性较强的成分，并以此为基础，建立鸡骨草质量控制的HPLC特征图谱，最终选择相思子碱、刺桐碱、夏佛塔苷等作为鸡骨草特征图谱的质量标志物。