整合计算研究Glycosmis pentaphylla来源化合物对Daboia russelii主要毒液蛋白的抑制作用

《Frontiers in Molecular Biosciences》：Integrated computational investigation of Glycosmis pentaphylla-derived compounds against major venom proteins of Daboia russelii

　　

蛇咬伤中毒仍是一种重要的被忽视热带病，尤其是在南亚地区，圆斑蝰（Daboia russelii）导致大量发病和死亡。当前抗蛇毒血清疗法因对局部组织损伤中和效果差、生产成本高、毒素特异性有限以及不良免疫反应而受到限制，凸显了对替代性或辅助性治疗策略的需求。研究人员采用整合计算与实验方法，评估了从Glycosmis pentaphylla叶提取物中鉴定的植物化合物作为主要毒液蛋白潜在抑制剂的可行性。通过气相色谱-质谱联用（GC-MS）鉴定的32种化合物，使用分子对接技术针对磷脂酶A2（PLA2）、蛇毒金属蛋白酶（SVMPs）、蛇毒丝氨酸蛋白酶（SVSPs）、L-氨基酸氧化酶（LAAO）和蛇C型凝集素样蛋白（SNACLECs）进行了筛选。对代表性高亲和力蛋白-配体复合物进行了100 ns分子动力学模拟。此外，还进行了药代动力学、毒性、蛋白质-蛋白质相互作用及通路富集分析，以评估所鉴定化合物的治疗潜力。初步实验验证采用胰蛋白酶-酪蛋白蛋白酶抑制试验、脱脂牛奶琼脂蛋白酶抑制试验和红细胞膜稳定试验。与N-乙酰半胱氨酸（NAC）相比，β-谷甾醇、3-氯-5-胆甾烯和维生素E对多种毒液蛋白表现出更强的结合亲和力和更稳定的相互作用。分子动力学分析（均方根偏差RMSD、均方根波动RMSF、回旋半径、溶剂可及表面积和氢键）进一步支持了所选蛋白-配体复合物的稳定性。药代动力学和毒性预测表明，大多数所选化合物具有可接受的类药性和相对安全的毒性特征。基于STRING的蛋白质-蛋白质相互作用和通路富集分析揭示了与毒液病理相关的炎症、凝血、氧化应激和细胞外基质相关通路的关联。G. pentaphylla提取物表现出浓度依赖性蛋白酶抑制活性（5 mg/mL时为56.86%）、降低的酪蛋白水解活性（抑制率28.19%）以及显著的红细胞膜稳定作用（39.23%），为计算预测提供了初步实验支持。总之，这些发现将G. pentaphylla来源的植物化合物确定为未来抗蛇毒血清开发的潜在先导化合物，并为后续针对蛇毒的特异性生化和体内验证提供了计算基础，以促进更安全、更经济的蛇咬伤管理辅助疗法的开发。

蛇咬伤中毒（snakebite envenomation）仍是被忽视的热带病，南亚地区圆斑蝰（Daboia russelii）咬伤导致严重发病与死亡。其毒液含有磷脂酶A2（phospholipase A2, PLA2）、蛇毒金属蛋白酶（snake venom metalloproteinases, SVMPs）、蛇毒丝氨酸蛋白酶（snake venom serine proteases, SVSPs）、L-氨基酸氧化酶（L-amino acid oxidase, LAAO）及蛇C型凝集素样蛋白（snake C-type lectin-like proteins, SNACLECs）等，引发凝血障碍、出血、组织坏死和肾衰竭。传统抗蛇毒血清虽为主流治疗，但存在中和局部损伤效果不佳、成本高、特异性有限及不良反应等局限。因此，植物来源抑制剂作为替代或辅助策略受到关注。Glycosmis pentaphylla（芸香科）为传统药用植物，具有抗炎、抗菌和抗毒活性。该研究首次系统探索其叶提取物来源化合物对D. russelii毒液蛋白的抑制作用，整合计算与初步实验验证，旨在为开发更安全、经济的蛇咬伤辅助治疗提供基础。

研究人员从前期GC-MS分析鉴定的32种化合物（14种来自甲醇提取物MME，13种来自石油醚提取物MPE，5种来自己烷提取物MHE）入手，使用PyRx（AutoDock Vina）将其对接到9个毒液蛋白靶标：PLA2（PDB ID: 1DPY, 1VIP, 2PVT, 2PYC）、SVSP（3S9B, 3SBK）、SVMP（2E3X）、LAAO（3KVE）和SNACLEC（3UBU）。基于对接亲和力与生物学意义，选择四个代表复合物（1VIP、2PYC、3S9B、2E3X）进行100 ns GROMACS分子动力学模拟，分析RMSD、RMSF、回旋半径、溶剂可及表面积和氢键。随后利用SwissADME和ProTox-II预测药代动力学与毒性，并通过STRING数据库进行蛋白质-蛋白质相互作用（PPI）和通路富集分析。实验验证采用胰蛋白酶-酪蛋白抑制试验、脱脂牛奶琼脂抑制试验和红细胞膜稳定试验，以胰蛋白酶作为丝氨酸蛋白酶模型，并以N-乙酰半胱氨酸（NAC）作为参照化合物。

结果部分：
配体选择与制备：共选定32种提取物来源化合物，结构取自PubChem并优化，适用于对接。蛋白选择与制备：选定9个毒液相关蛋白，覆盖PLA2、SVSP、SVMP、LAAO、SNACLEC五类毒素家族。分子对接分析：32种化合物对9种毒液蛋白的结合能范围为?3.9至?9.3 kcal/mol，显著优于NAC（?4.0至?5.0 kcal/mol）。β-谷甾醇（β-sitosterol）、维生素E（vitamin E）和3-氯-5-胆甾烯（3-chloro-5-cholestene）表现出广谱强结合，其中β-谷甾醇对PLA2（1DPY）和LAAO（3KVE）的结合能均达?9.3 kcal/mol。对接结果显示这些化合物与PLA2催化残基His48、Asp49及SVSP的Arg62、Lys73等关键残基形成氢键和疏水接触，提示可能抑制酶活性。RMSD分析：3-氯-5-胆甾烯和维生素E复合物较NAC具有更低的构象漂移（如1VIP复合物RMSD为0.15–0.18 nm），表明配体能稳定蛋白整体构象。RMSF分析：这些配体降低催化位点残基的波动，例如围绕PLA2催化二联体His48、Asp49、Tyr52的波动从>0.12 nm降至0.08–0.10 nm，提示可限制活性位点运动。回旋半径（Rg）分析：配体结合使蛋白结构更紧凑，如SVMP（2E3X）复合物Rg为1.79–1.82 nm，而NAC为1.88–1.91 nm，表明配体诱导致密化。氢键分析：3-氯-5-胆甾烯可与催化残基形成3–5个持久氢键，而NAC通常仅形成1–2个瞬时氢键，支持更强结合。溶剂可及表面积（SASA）分析：3-氯-5-胆甾烯使SASA降至112–115 nm2，显著低于NAC的约125–130 nm2，说明配体深埋于活性口袋，减少溶剂暴露和底物进入。实验验证：G. pentaphylla提取物在胰蛋白酶-酪蛋白试验中呈浓度依赖性抑制（1 mg/mL时19.61%至5 mg/mL时56.86%）；脱脂牛奶琼脂试验中使酪蛋白水解区从17.17 mm减至12.33 mm，抑制率28.19%；红细胞膜稳定试验中在5 mg/mL时达39.23%膜稳定率，均具有统计学显著性。药代动力学与毒性预测：SwissADME预测显示大多数化合物具有可接受的类药性和生物利用度分数（约0.55），但甾醇和萜类化合物亲脂性高（LogP 8.02–10.60），存在Lipinski规则违例，可能限制全身生物利用度。ProTox-II预测维生素E和角鲨烯毒性较低（LD₅₀ 5000 mg/kg），而6-硝基-2-甲基吡咯并[2,3-g]喹啉显示致突变和致癌风险。总体而言，β-谷甾醇和维生素E安全性较好。蛋白质-蛋白质相互作用网络分析：STRING网络显示TNF、IL6、MMP2、MMP9、F2、F10、SERPINE1、PLAT等靶点紧密关联，涉及炎症、凝血和细胞外基质降解。KEGG通路富集分析：显著富集补体与凝血级联、TNF信号通路、IL-17信号通路及细胞因子介导的炎症通路。Reactome通路富集分析：显著富集细胞因子信号转导、止血、血小板活化、细胞外基质组织和免疫系统调节通路。基因本体生物学过程分析：显著富集炎症反应、凝血调节、细胞外基质分解、氧化应激反应和细胞因子介导的信号通路，与蛇毒病理机制一致。

讨论部分总结了该研究的几点意义：植物来源化合物具有多靶点抑制潜力，分子动力学模拟弥补了静态对接的不足；NAC作为参照具有抗氧化和细胞保护作用，但并非选择性毒素抑制剂；实验采用胰蛋白酶模型和红细胞膜稳定模型，虽不能完全模拟真实蛇毒中毒，但为计算预测提供了初步功能学支持。未来需使用纯化毒液酶、全毒液中和试验和体内模型进一步验证。研究结论翻译如下：本研究首次提供了G. pentaphylla成分作为D. russelii主要毒液蛋白未来抗蛇毒开发潜在先导化合物的综合计算证据。分子对接和分子动力学分析显示，β-谷甾醇、3-氯-5-胆甾烯和维生素E与磷脂酶A2、蛇毒金属蛋白酶和丝氨酸蛋白酶等多个毒液相关靶标形成稳定且有利