《Frontiers in Bioinformatics》：Targeting KasA: isosakuranetin derivatives as promising scaffolds for novel anti-tuberculosis agents against drug-resistant Mycobacterium tuberculosis
编辑推荐：
为应对多重耐药和广泛耐药结核病带来的治疗困境，研究人员聚焦于结核分枝杆菌细胞壁合成关键酶KasA，采用计算机辅助药物设计手段，对源自药用植物Chromolaena odorata的香橙素(Isosakuranetin)衍生物进行筛选。研究发现，衍生物Isn_96在分子对接、动力学模拟及ADMET预测中均表现出最优的结合亲和力、稳定性与药代动力学特性，是极具前景的KasA抑制剂候选药物。此项研究为开发新型抗结核药物提供了强有力的先导化合物和理论依据。
在人类与传染病的漫长斗争中，结核病（Tuberculosis, TB）始终是挥之不去的阴影。由结核分枝杆菌（Mycobacterium tuberculosis, Mtb）引起的这种古老疾病，至今仍在全球范围内肆虐，每年造成数百万人感染和上百万人死亡。其“顽固”不仅在于高感染性，更在于其病原体拥有一件强大的“盔甲”——一层由分枝菌酸（mycolic acid）构成的、异常致密且富含脂质的细胞壁。这层“盔甲”不仅能帮助细菌逃避免疫系统的攻击，还构成了抵御多种抗生素的物理屏障，是结核分枝杆菌产生固有耐药性和持久性的关键。
更棘手的是，随着抗生素的长期和不当使用，能够抵抗多种一线药物的多重耐药结核病（Multi-Drug Resistant TB, MDR-TB）以及对更多药物耐药的广泛耐药结核病（Extensively Drug Resistant TB, XDR-TB）不断涌现。目前世界卫生组织（WHO）推荐的标准治疗方案（包含利福平、异烟肼、吡嗪酰胺和乙胺丁醇）疗程长达6-9个月，不仅易导致患者依从性差，还存在肝毒性等副作用，对耐药菌株的疗效更是有限。尽管出现了贝达喹啉、德拉马尼等新药，但耐药性问题依然存在。因此，寻找作用于全新靶点、能有效穿透或破坏细菌“盔甲”的新型药物，已成为抗击结核病，特别是耐药结核病的迫切需求。
科学家们将目光投向了细菌合成“盔甲”核心材料——分枝菌酸的“工厂”。在结核分枝杆菌中，这个“工厂”由两个系统协同运作：FAS-I系统（哺乳动物型I型脂肪酸合成）负责生产短链脂肪酸前体，而FAS-II系统（细菌型II型脂肪酸合成）则负责将其“加长”为构成分枝菌酸骨架的长链分支菌酸。其中，FAS-II系统中的一个关键“装配工”——β-酮脂酰载体蛋白合酶A（β-ketoacyl acyl carrier protein synthase A, KasA）成为了理想的药物靶点。KasA催化分枝菌酸链延伸的核心步骤，且是人类所不具备的细菌特有酶，这意味着针对它设计的药物可能具有高选择性，在杀死细菌的同时对宿主细胞影响较小。抑制KasA，就等于切断了“盔甲”原料的供应链，有望导致细菌因细胞壁缺陷而死亡。
与此同时，大自然一直是新药发现的重要宝库。药用植物Chromolaena odorata因其抗菌、抗结核活性而被传统医学所认知。从该植物中分离出的黄酮类化合物香橙素（Isosakuranetin）已被报道具有中度的抗结核活性。这为药物研发提供了宝贵的“种子”化合物。然而，天然产物的活性往往有待提高，需要通过化学修饰来优化其与靶点蛋白的结合能力、稳定性和成药性。那么，能否以香橙素为“蓝图”，通过计算机设计并筛选出更强效、更安全的KasA抑制剂呢？
为了回答这一问题，一项发表于《Frontiers in Bioinformatics》的研究应运而生。研究人员运用了一套综合的计算生物学“组合拳”，旨在从一系列香橙素衍生物中，筛选出靶向结核分枝杆菌KasA酶的最佳候选药物。
研究人员主要运用了以下几个关键技术方法：首先，他们采用结构生物学方法，从蛋白质数据库（RCSB-PDB）中获取了结核分枝杆菌KasA同源二聚体与已知抑制剂TLM4的共晶结构（PDB: 4C70）作为受体模型。其次，利用计算机辅助药物设计（CADD）工具，以香橙素为母核，通过片段库合成和结构修饰，生成了一系列衍生物库。接着，通过分子对接（Molecular Docking）模拟，评估了这些化合物与KasA活性位点的结合模式和亲和力，并利用MM-GBSA（Molecular Mechanics with Generalized Born and Surface Area）方法计算结合自由能。然后，对筛选出的高亲和力化合物进行了ADMET（吸收、分布、代谢、排泄和毒性）性质预测，以评估其成药潜力。最后，对最有前景的化合物与KasA的复合物进行了长达250纳秒的分子动力学模拟（Molecular Dynamics Simulation, MDS），以评估复合物在接近生理条件下的动态稳定性，并通过轨迹分析计算了结合自由能、残基波动、溶剂可及表面积等关键参数。
分子对接筛选高亲和力化合物
研究人员首先验证了分子对接流程的可靠性，将共晶配体TLM4重新对接到KasA活性位点，取得了0.59 ?的低均方根偏差（RMSD）。随后，他们对香橙素及其六个衍生物进行了对接筛选。结果显示，所有化合物的对接分数在-5.712至-7.921 kcal/mol之间。其中，衍生物Isn_96表现最为突出，其对接分数最低，为-7.921 kcal/mol，预示着最强的结合亲和力。分子相互作用分析表明，Isn_96能与KasA活性位点的关键残基，如HIS311、HIS345、ASP273、PHE404和ASN408等，形成稳定的氢键、π-π堆积和疏水相互作用，这表明它能有效占据并稳定在KasA的底物结合隧道中。
ADMET性质预测评估成药潜力
在进入更耗时的动力学模拟前，研究人员对筛选出的化合物进行了ADMET性质预测，以提前排除可能存在严重药代动力学或毒性问题的分子。分析显示，Isn_96展现出较为平衡的ADMET特征：具有相对较高的水溶性、较低的血浆蛋白结合率（34.4%）、适中的半衰期以及较低的预测毒性风险。相比之下，原型化合物香橙素预测有较高的肝代谢率和血液毒性风险，而某些其他衍生物则在渗透性、P-糖蛋白（P-gp）底物特性或特定毒性指标上存在短板。综合来看，Isn_96、Isn_184_A67、香橙素原型以及标准抑制剂被选入后续的分子动力学模拟进行深入评估。
分子动力学模拟验证结合稳定性
为了在动态、接近生理的环境中验证复合物的稳定性，研究人员对上述四个化合物与KasA的复合物进行了250纳秒的分子动力学模拟。
- •
均方根偏差（RMSD）分析：RMSD反映了蛋白质骨架结构随时间变化的偏差。Isn_96-KasA复合物的Cα原子RMSD平均值最低（2.65 ± 0.45 ?），且波动最小，表明该复合物在整个模拟过程中结构最为稳定。而标准抑制剂复合物的RMSD最高（4.75 ± 0.85 ?），稳定性相对较差。
- •
均方根涨落（RMSF）分析：RMSF揭示了蛋白质各个残基的柔性。分析发现，Isn_96复合物在催化关键残基（如HIS311, HIS345）所在区域表现出较低的涨落，说明配体的结合使活性位点区域更加刚性稳定。而在一些非活性区的环状结构处，所有体系都表现出固有的柔性，这与预期一致。
- •
蛋白质-配体相互作用分析：轨迹分析显示，Isn_96能与KasA的催化三联体残基HIS311和HIS345保持长时间（分别约86%和92%的模拟时间）的强相互作用。同时，它与GLN171、ASP273、ASN408等残基也形成了稳定的氢键网络。这些持久的相互作用确保了Isn_96能牢固地“锁”在KasA的活性口袋中。
- •
结合自由能计算：基于分子动力学轨迹的MM-GBSA计算给出了更精确的结合自由能估计。Isn_96的结合自由能（ΔGbind）为-56.20 ± 6.90 kcal/mol，是所有测试化合物中最负的，表明其与KasA的结合在热力学上最为有利。能量分解显示，其强大的结合力主要得益于显著的库仑静电相互作用贡献（-96.98 ± 16.85 kcal/mol）。
研究结论与讨论
综上所述，这项系统的计算生物学研究成功地筛选并验证了香橙素衍生物Isn_96作为靶向结核分枝杆菌KasA酶的高潜力抑制剂。Isn_96在分子对接中显示出最强的初始结合亲和力，在长达250纳秒的分子动力学模拟中表现出优异的复合物稳定性，其药代动力学（ADMET）预测特征也较为理想，并且计算得出的结合自由能远优于原型化合物和参考抑制剂。
这项研究的意义在于：首先，它明确指出了源自植物的香橙素衍生物，特别是Isn_96，是开发新型抗结核药物的一个极具前景的“先导结构”（scaffold）。这为基于天然产物的抗结核药物研发提供了新的思路和具体的候选分子。其次，研究证实了KasA作为一个抗结核药物靶点的“可成药性”，并详细阐释了高效抑制剂与靶点之间的相互作用模式，为后续的理性药物设计提供了关键信息。最后，研究展示了一套整合分子对接、ADMET预测和分子动力学模拟的计算工作流程，在药物发现早期阶段进行高效筛选和优化方面的有效性与价值。
当然，论文作者也明确指出，目前的所有发现均基于计算机模拟。Isn_96的真正效力与安全性，需要通过后续的化学合成、体外抗菌活性测试、细胞毒性实验以及更复杂的体内动物模型研究来最终验证。然而，这项计算研究无疑为抗击耐药结核病的新药探索之旅点亮了一盏明灯，指明了前进的方向，并为后续的实验验证奠定了坚实的理论基础。
