摘要

目的 采用孟德尔随机化（MR）及分子动力学模拟（MD）分析附子通过调控免疫细胞治疗慢性心力衰竭（CHF）的机制。方法 通过中药系统药理学数据库分析平台（TCMSP）及Swiss ADME在线工具获得附子有效成分，采用SwissTargetPrediction数据库获取附子靶点。采用双样本MR分析附子靶点、免疫细胞、CHF之间的因果关联，将MR-Egger回归、加权中位数法作为逆方差加权法（IVW）的补充方法。采用Cochran's Q检验、MR-Egger回归的截距项及留一法进行敏感性分析；采用MD评估附子活性成分与附子核心靶点蛋白的结合模式及稳定性。结果 IVW分析结果显示：（1）27个附子靶点中仅内质网α-甘露糖苷酶1（MAN1B1）与CHF发生风险呈负相关（P=0.005）。MR-Egger回归分析结果显示，MAN1B1与CHF发生风险无因果关系（P>0.05），但其β值与IVW的β值方向一致。加权中位数法分析结果显示，MAN1B1与CHF发生风险呈负相关（P=0.004）。（2）6种免疫细胞亚群中仅CD25hi%CD4+与CHF发生风险呈正相关（P=0.016），仅CD64onCD14-CD16+monocytes与CHF发生风险呈负相关（P=0.024）。MR-Egger回归、加权中位数法分析结果显示，CD25hi%CD4+、CD64onCD14-CD16+monocytes与CHF发生风险无因果关系（P>0.05），且MR-Egger回归的β值与IVW的β值方向不一致，而加权中位数法的β值与IVW的β值方向一致。（3）MAN1B1与CD25hi%CD4+呈负相关（P=0.040），与CD64onCD14-CD16+monocytes呈正相关（P=0.044）。MR-Egger回归、加权中位数法分析结果显示，MAN1B1与CD25hi%CD4+、CD64onCD14-CD16+monocytes无因果关系（P>0.05），但其β值与IVW的β值方向一致。Cochran's Q检验分析结果显示，与MAN1B1、CD25hi%CD4+、CD64onCD14-CD16+monocytes强相关的单核苷酸多态性（SNP）间不存在统计学异质性（P>0.05）；MR-Egger回归的截距项分析结果显示，与MAN1B1、CD25hi%CD4+和CD64onCD14-CD16+monocytes强相关的SNP间不存在水平多效性（P>0.05）。留一法分析结果显示，剔除单个SNP后，MR分析结果无明显改变。MD分析结果显示，附子活性成分中的惰碱与MAN1B1的结合能<-7 kcal/mol，且惰碱-MAN1B1复合物在MD过程中保持稳定。结论 附子活性成分惰碱可与其靶点MAN1B1稳定结合，并通过抑制CD25hi%CD4+表达、促进CD64onCD14-CD16+monocytes表达而延缓CHF病情进展。

关键词