β-拉帕醌衍生物靶向利什曼原虫蛋白揭示一个有前景的多靶点候选物

《Archives of Microbiology》：Targeting of Leishmania proteins by β-Lapachone derivatives reveals a promising multitarget candidate

　　

摘要：利什曼病（Leishmaniasis）是全球主要公共卫生问题之一。犬利什曼病（Canine leishmaniasis, CanL）在这一疾病的持续传播中扮演着重要角色，因为现有的治疗方案有限、成本高昂、需要长期治疗，并且与日益增长的耐药性有关。在此背景下，β-拉帕醌（β-lapachone），一种从拉帕醇（lapachol）衍生而来的生物活性萘醌，因其药理特性（包括抗寄生虫活性）而成为一种潜在的治疗替代方案，这使其成为开发新型抗利什曼病药物的有希望的候选物。在本研究中，研究人员采用整合计算方法评估了一组β-拉帕醌衍生物的抗利什曼病潜力和药代动力学特征。分子对接（Molecular docking）模拟表明，与现有药物相比，这些化合物对利什曼原虫（Leishmania spp.）靶点具有更高的选择性。在分析的衍生物中，化合物LPI3脱颖而出，成为先导化合物（lead compound），显示出有前景的多靶点（multitarget）特征。基于简正模分析（Normal Mode Analysis, NMA）的分子动力学分析显示，涉及LPI3的蛋白质-配体复合物具有高的构象稳定性，与在参考药物中观察到的稳定性相似。关于吸收、分布、代谢、排泄和毒性（ADMET）的药代动力学描述符的预测表明，该先导化合物具有良好的药代动力学特征，符合药物相似性和安全性标准。因此，研究结果确定LPI3是未来在利什曼病治疗中进行实验评估的潜在候选化合物。

**β-拉帕醌衍生物作为抗利什曼病多靶点候选物的计算研究**

**1. 研究背景与研究必要性**

利什曼病（Leishmaniasis）是由利什曼原虫（*Leishmania* spp.）引起的、在全球热带和亚热带地区广泛流行的被忽视热带病，对人类和动物健康构成重大威胁。犬作为利什曼原虫的主要贮存宿主，其犬利什曼病（Canine leishmaniasis, CanL）在维持疾病传播中起关键作用。目前临床上使用的抗利什曼病药物，如五价锑剂、两性霉素B（amphotericin B, AmphoB）、米替福新（miltefosine, MTF）和巴龙霉素等，普遍存在疗效有限、成本高昂、需要长期治疗以及耐药性日益增加等局限。因此，迫切需要开发更安全、更有效和更具选择性的新型治疗药物。

天然产物因其结构多样性和广谱生物活性，成为药物开发的重要来源。β-拉帕醌（β-lapachone）是一种来源于拉帕醇的邻萘醌，具有抗肿瘤、抗炎和抗寄生虫等多种药理活性，显示出作为抗利什曼病药物先导化合物的潜力。同时，计算方法的快速发展为药物发现提供了强有力的工具，可以深入探究配体与靶点蛋白的相互作用及复合物稳定性，并预测化合物的药代动力学和毒理学性质，从而有效降低研发成本和周期。

本研究旨在通过整合计算策略，系统评估一组β-拉帕醌衍生物与利什曼原虫关键蛋白靶点的相互作用、药代动力学特征及药理活性，筛选出具有多靶点作用潜力的先导化合物，为抗利什曼病药物研发提供理论依据和候选分子。该研究成果发表在《Archives of Microbiology》上。

**2. 主要关键技术方法**

研究人员采用了一套整合的计算方法流程。首先，通过分子对接（molecular docking）模拟技术，使用AutoDock Vina软件，将β-拉帕醌衍生物与三种利什曼原虫关键蛋白靶点——二氢乳清酸脱氢酶（dihydroorotate dehydrogenase, DHODH）、N-肉豆蔻酰基转移酶（N-myristoyltransferase, NMT）和木瓜蛋白酶样半胱氨酸蛋白酶B（cathepsin B, CatB）——进行对接，预测结合模式和亲和力。其次，基于简正模分析（Normal Mode Analysis, NMA），利用iMODS服务器对优选出的复合物进行构象动力学分析，评估其稳定性、柔性和变形性。最后，利用多种在线服务器（如SwissADME、pkCSM、ADMETlab 3.0）对先导化合物进行吸收、分布、代谢、排泄和毒性（ADMET）属性的计算机模拟（in silico）预测。

**3. 研究结果**

**蛋白质结构验证与口袋预测**：通过拉马钱德兰图（Ramachandran plot）验证了DHODH、NMT和CatB三个靶点蛋白的三维结构，其立体化学质量良好，无残基位于不允许区域。DoGSite3预测的口袋分析显示，DHODH和NMT具有较大且化学组成多样的结合口袋，而CatB的最佳口袋体积较小且功能组织更为复杂。

**分子对接模拟**：所有β-拉帕醌衍生物与靶点蛋白的结合都表现出低均方根偏差（RMSD<2.0 ?）和低亲和能（E_A），优于或媲美参考药物（MTF和AmphoB）及共结晶抑制剂。在七个衍生物中，化合物LPI3在所有三个靶点中均表现出最低的RMSD和E_A值。RMSD与E_A的相关性四象限图显示，LPI3主要分布在第三象限（Q3），表明其兼具高构象一致性和高理论亲和力，是选择性最高的衍生物。

**相互作用分析**：LPI3与DHODH主要形成疏水相互作用和π-阳离子相互作用；与NMT形成疏水相互作用、氢键和π-π堆积作用；与CatB形成疏水相互作用和氢键。与其他衍生物及参考药物相比，LPI3的结合模式更优，且能与关键活性位点残基发生特异性相互作用。

**基于NMA的构象动力学分析**：简正模分析显示，与游离蛋白和参考药物复合物相比，LPI3与DHODH、NMT和CatB形成的复合物在整体上具有更低的柔性（RMSF）和变形性。所有蛋白-配体复合物的变形性谱图高度相似，且LPI3复合物的变形性值略低于对照组。构象趋近分析表明，LPI3复合物能与实验结构的靶蛋白良好趋近，说明其构象稳定。

**多参数优化分析与ADMET预测**：LPI3具有分子量为378.48 g/mol、拓扑极性表面积（TPSA）为37.91 ?²、高脂溶性（logP=6.32）的理化性质。其口服生物利用度评分为0.55，符合利平斯基五规则（Lipinski’s Rule），但未满足辉瑞规则的两项标准。多参数优化（MPO）得分为3.6，落入可接受范围。ADMET预测表明，LPI3的水溶性较低，但具有良好的肠道吸收率（79.255%），可能与P-糖蛋白（Pgp）发生相互作用。代谢预测显示，该化合物是CYP1A2、CYP2C9和CYP2C19的抑制剂，同时是CYP2D6和CYP3A4的底物，其二期代谢可能通过尿苷二磷酸葡萄糖醛酸转移酶（UGT）介导发生葡萄糖醛酸化。毒性预测结果表明，LPI3无肝毒性和皮肤致敏性，但具有中等程度的hERG（human ether-a-go-go-related gene）通道抑制风险，并表现出一定的口服急性毒性和致突变性。

**4. 讨论与结论总结**

讨论部分指出，利什曼原虫特异的靶点DHODH和NMT因其较大的口袋和多样的化学组成，为衍生物提供了多种相互作用模式。而CatB因其最佳口袋与催化位点分离，结合位点的差异并未导致复合物稳定性的显著变化。LPI3在分子对接和构象动力学分析中均表现出优于参考药物的性能，体现了其作为多靶点抑制剂的潜力。尽管其药代动力学和毒性特征较为有利，但中等程度的hERG风险和口服急性毒性提示了进一步优化的必要性。

**研究结论翻译**：“本研究结果表明，所采用的整合计算策略在评估β-拉帕醌衍生物的抗利什曼病潜力方面是有效的，能够对与利什曼原虫（*Leishmania* spp.）必需靶点的结构、药代动力学和分子相互作用特性进行全面表征。在评估的化合物中，LPI3衍生物脱颖而出，成为最有希望的候选化合物，其表现出与药物相似性和安全性标准相一致的有利的理化和药代动力学特征。此外，对接和基于简正模分析的构象动力学分析表明，该化合物与分子靶点建立了稳定且能量有利的相互作用，其性能优于参考药物，尤其是与利什曼原虫属特异性蛋白的结合，暗示其具有更强的选择性和生物活性潜力。” 这些发现证实了LPI3作为抗利什曼病候选药物的潜力，强调了计算机模拟（in silico）策略在生物活性化合物筛选和优化中的重要性，但仍需进行体内外实验验证。