针对耐甲氧西林金黄色葡萄球菌(MRSA)中青霉素结合蛋白2a(PBP2a)这一高度柔性靶点的药物发现挑战,研究人员开发了一种整合构象系综分析与深度生成建模的综合框架。该研究通过分子动力学模拟绘制PBP2a的动态构象景观并识别隐秘结合口袋,随后利用靶标条件化的深度生成模型探索这些状态以产生新型候选骨架,并通过多目标优化和结构优先级排序筛选出化合物1。该化合物在1000纳秒的分子动力学模拟中表现出保守的结合模式和持续稳定性,其MMGBSA结合能优于天然配体。全细胞实验显示其对金黄色葡萄球菌ATCC 29213和MRSA ATCC 43300具有剂量依赖性的生长抑制作用,且24小时暴露后活性更强。碘化丙啶染色进一步证实了浓度依赖性的膜通透性增加,表明化合物诱导了包膜扰动。这项研究提出了一种针对动态柔性蛋白的药物发现通用工作流程,并为抗MRSA感染提供了新的先导化合物线索。
《Frontiers in Microbiology》:Integrated framework for targeting dynamic penicillin-binding protein 2a via ensemble structural biology and deep generative modeling
编辑推荐:
针对耐甲氧西林金黄色葡萄球菌(MRSA)在全球范围内导致的高发病率和死亡率,研究人员开发了一种整合构象系综分析与深度生成建模的综合发现框架,以解决青霉素结合蛋白2a(PBP2a)——一种由mecA基因编码的低亲和力转肽酶——这一高度柔性靶点的挑战。通过广泛的分
针对耐甲氧西林金黄色葡萄球菌(MRSA)在全球范围内导致的高发病率和死亡率,研究人员开发了一种整合构象系综分析与深度生成建模的综合发现框架,以解决青霉素结合蛋白2a(PBP2a)——一种由mecA基因编码的低亲和力转肽酶——这一高度柔性靶点的挑战。通过广泛的分子动力学(MD)模拟绘制PBP2a的动态景观并识别隐秘结合口袋,随后采用靶标条件化的生成工作流程探索这些状态以产生新型候选骨架,并进行多目标优化和结构优先级排序。在生成的分子中,化合物1显示出保守的结合模式,在1000纳秒的分子动力学模拟中保持持续稳定性,且具有比天然配体更有利的MMGBSA结合能。全细胞实验显示其对金黄色葡萄球菌ATCC 29213和MRSA ATCC 43300具有剂量依赖性的生长抑制作用,且在24小时暴露后活性更强。额外的碘化丙啶(PI)染色显示两种菌株中PI摄取呈浓度依赖性增加,与化合物诱导的包膜扰动一致。这些发现支持化合物1作为一种源自PBP2a引导发现策略的抗菌先导化合物,但其直接靶标结合机制仍有待确立。
针对耐甲氧西林金黄色葡萄球菌(MRSA)引发的全球公共卫生危机,由于细菌表达的青霉素结合蛋白2a(PBP2a)具有高度动态的闭合活性位点,传统基于静态晶体结构的药物设计往往失效。为此,研究人员在《Frontiers in Microbiology》发表研究,提出了一种整合构象系综分析与深度生成建模的综合框架,旨在突破针对此类柔性靶点的药物发现瓶颈。
为实现这一目标,研究人员采用了多项关键技术方法。首先,通过对3ZG5和5M18晶体结构进行广泛的分子动力学(MD)模拟及聚类分析,构建了包含20种代表性构象的靶标系综。其次,基于DeepDTA架构构建了亲和力预测模型,并结合进化策略进行分子生成,利用BRICS片段重组和突变初始化种群。随后,采用GNINA进行多阶段对接筛选,并利用1000纳秒的MD模拟及分子力学/广义波恩表面积(MMGBSA)方法验证结合自由能。最后,通过肉汤微量稀释法测定最低抑菌浓度(MIC),并结合碘化丙啶(PI)染色评估体外抗菌活性及膜完整性影响。
研究结果详细阐述了从计算设计到实验验证的全过程。模型训练稳定性与进化优化动力学评估显示,DeepDTA模型训练损失迅速收敛,验证集R2达到0.787,表明模型能有效捕捉构效关系(SAR)。进化轨迹分析揭示了严格的硬过滤器(如毒性警报、分子量限制)对种群分布的影响,最终从120个独特分子中筛选出82个符合标准的候选物。多目标优化与化学空间表征表明,生成分子在定量药物相似性(QED)、集合评分和合成可及性(SA)之间取得了平衡。主成分分析(PCA)证实生成分子远离种子化合物,具有高度的结构新颖性,且类药性分布符合利平斯基五规则。药物相似性与特征重要性分析指出,大多数分子无利平斯基违规,排除原因主要为毒性警报和分子量过低。随机森林分析确定QED和分子量(MW)是驱动集合评分的主要因素。分子对接验证与结合相互作用分析结果显示,化合物1在3ZG5和5M18两种受体模型中表现出高度保守的结合模式,共享ARG-241、ARG-151、VAL-277和TYR-373等关键残基,而化合物2和3的保守性差,因此化合物1被选为最佳先导。化合物1的动态稳定性与结合能学研究表明,在1000纳秒MD模拟中,化合物1的蛋白质骨架和配体原子的均方根偏差(RMSD)低于共晶配体,且MMGBSA结合自由能始终优于天然抑制剂,证实了其卓越的结构稳定性。化合物1的抗菌活性评估显示,其对金黄色葡萄球菌ATCC 29213和MRSA ATCC 43300表现出浓度和时间依赖性的生长抑制,24小时活性强于18小时。PI染色进一步观察到浓度依赖性的荧光增强,表明化合物诱导了包膜应激表型。
讨论部分总结认为,这项研究成功地将系综结构生物学与深度生成建模相结合,克服了静态结构药物设计的局限性。通过映射PBP2a的动态构象景观,研究人员识别了隐秘口袋,并利用生成模型探索了广阔的化学空间,实现了结合亲和力、药物相似性和合成可行性的平衡。化合物1被确定为最具潜力的先导化合物,其保守的结合模式、优异的MD稳定性及体外抗菌活性支持了PBP2a引导的发现策略的有效性。结论部分强调,该研究提出的整合工作流程通过明确纳入蛋白质动力学,改善了针对柔性靶点的命中发现,支持化合物1作为一种有前景的抗MRSA感染先导化合物,但直接生化验证仍需后续研究。