《Journal of Molecular Graphics and Modelling》：Quinolone-Inspired MCM2-7 Helicase Inhibitors: Computational Design, Dynamic Stability, and Preclinical Promise for Targeted Anticancer Therapy
编辑推荐：
本研究通过结构虚拟筛选、分子动力学模拟和MM-PBSA计算，设计并筛选出12个新型喹诺酮类化合物抑制MCM2-7复合体的AAA+ ATP酶域。结果显示化合物1、7、8结合亲和力最高（-8.66 kcal/mol），其中化合物8在稳定OB域和ATP酶域方面效果最佳，化合物1形成稳定氢键网络（均7.2±1.8键）。MM-PBSA分析表明化合物1结合自由能最优（-51.27±8.46 kJ/mol），ADMET预测显示低突变风险和半衰期延长（≥3h），但口服生物利用度有限。该研究构建了高效计算药物设计流程，为靶向复制压力的抗癌治疗提供候选分子。
Seifeldin Elabed | Wael M. Elshemey | Ammar Y. Mohamed | Medhat W. Shafaa
埃及赫尔万大学理学院物理系，医学生物物理学部门
摘要
迷你染色体维护（MCM）2–7解旋酶复合体是真核生物DNA复制机制中的关键组成部分，由于其在复制许可和基因组稳定性中的作用，也成为肿瘤学领域的一个新兴治疗靶点。尽管其生物学重要性显著，但针对该复合体的小分子抑制剂却很少被合理设计出来，尤其是针对人类MCM2的抑制剂。在这项研究中，我们采用了基于结构的虚拟筛选策略，结合分子动力学（MD）模拟、MM-PBSA结合自由能计算以及全面的药代动力学建模，以鉴定出新型喹诺酮类抑制剂，这些抑制剂能够靶向MCM2的AAA+ ATP酶结构域。
十二种经过合理设计的化合物被对接到经过优化的AlphaFold3衍生的MCM2结构上。配体1、7和8表现出最强的结合亲和力（高达–8.66 kcal/mol），优于参考化合物环丙沙星。100纳秒内的MD模拟显示，配体8最有效地稳定了OB折叠结构域和ATP酶结构域，而配体1形成了高度持久的氢键网络（平均7.2 ± 1.8个氢键）。MM-PBSA能量分解表明，配体1具有最有利的结合自由能（–51.27 ± 8.46 kJ/mol），这归因于其强烈的极性溶剂化作用。计算机模拟的ADMET预测显示这些化合物具有较低的突变风险、较长的半衰期（≥3小时）以及良好的药物类似性，尽管其口服生物利用度仍然有限。
总体而言，这项工作建立了一个用于解旋酶靶向药物设计的强大计算流程，并鉴定出了具有优异结合力和药理特性的有前景的MCM2抑制剂。这些发现为配体诱导的复制机制稳定提供了机制上的见解，并支持进一步在体外和体内探索这些化合物作为利用复制压力脆弱性的抗癌治疗手段的可能性。
引言
基因组的准确和及时复制对于细胞增殖和生物体的生存至关重要。这一过程的核心是一系列蛋白质复合体的协同作用，它们确保了细胞周期中DNA复制的准确性。其中，迷你染色体维护（MCM）复合体作为复制机制的关键组成部分，起着主要解旋酶的作用，负责解开DNA双螺旋结构——这是复制叉进展的前提条件。
MCM2-7复合体是由六个不同但相关的亚基组成的异六聚体，它通过促进复制起始位点上预复制复合体（pre-RC）的形成来协调DNA复制的启动。在G1阶段，复制起始位点通过加载MCM2-7复合体而被“许可”，这一过程受到多个细胞周期检查点和激酶信号级联的严格调控。在S阶段被激活后，MCM2-7的解旋酶功能驱动双链DNA的分离，为DNA聚合酶和其他复制因子的组装提供单链模板。这一基本作用不仅强调了它在正常细胞功能中的重要性，也在复制准确性受损的情况下（如肿瘤发生过程中）显得尤为关键。
新的证据进一步揭示了调控MCM2-7解旋酶的多方面机制。这些机制包括其亚基之间的复杂相互作用、翻译后修饰以及与各种调节蛋白的相互作用，这些因素共同微调其活性。任何这些调节层中的紊乱或改变都可能导致复制压力——这是基因组不稳定性的已知因素，也是癌症发展的驱动力。因此，MCM2-7复合体作为细胞增殖潜力的关键标志物以及治疗干预的有希望的目标，特别是在那些以细胞分裂失控为特征的癌症中，受到了越来越多的关注。
然而，开发针对真核生物解旋酶（包括MCM2-7）的有效且特异的抑制剂面临着重大挑战。例如，现有的用于解旋酶抑制的化合物（如抗菌化合物环丙沙星）存在局限性，包括对MCM2-7靶点的效力较低，以及由于其主要药理活性而产生的不良副作用，这使得它们不适合用于专门的抗癌治疗。这凸显了发现和开发新型化学类型抑制剂的重要需求，这些抑制剂能够以改进的效力和安全性特异性调节MCM2-7解旋酶活性。
MCM2-7解旋酶的复杂性也体现在其结构动力学上。该复合体包含负责ATP结合、核酸相互作用和亚基寡聚化的不同结构域，构成了一个高度适应性强的分子机器。这种适应性使其能够在复制周期中精确调节，并为通过靶向破坏关键功能域来调节其活性提供了潜在途径。值得注意的是，研究已经开始揭示解旋酶构象的微妙变化如何对复制完整性产生深远影响，为旨在操控这一关键复制因子的创新药物发现工作开辟了新的方向。
针对真核生物DNA解旋酶的小分子抑制剂的研发是一个活跃且不断发展的研究领域，其动力在于调节DNA复制在过度增殖性疾病中的治疗潜力。该领域的努力越来越集中在识别能够选择性干扰这些关键酶的ATP酶活性或DNA解旋功能的化合物上，从而为癌症治疗提供新的策略。
总之，MCM2-7解旋酶复合体在控制DNA复制中占据着关键地位，整合了多层的调节机制以确保基因组稳定性。它在DNA复制的启动和进展中的核心作用不仅使其对细胞存活不可或缺，还使其成为开发新型治疗策略的理想靶点。随着研究的深入，以及对靶向解旋酶抑制剂开发的持续努力，对其机制基础的更深入理解有望为那些以异常细胞增殖为特征的疾病提供新的干预措施。
在这项工作中，为了克服现有抑制剂的局限性，我们描述了十二种新型喹诺酮类衍生物的设计和计算评估，这些衍生物靶向人类MCM2的AAA+ ATP酶结构域。这一重点的提出是基于现有药物（如环丙沙星）所面临的挑战：尽管它们能优先阻断解旋酶活性，但其IC50值高达约632 μM，并且具有不适合抗癌治疗的抗菌活性。通过基于结构的对接、长时间尺度的分子动力学模拟以及计算机模拟的ADMET预测，我们鉴定出了具有增强结合亲和力、选择性和药物样特性的先导化合物，以便进一步进行实验验证。
结构预测、优化与评估
本研究的初始阶段涉及准确预测和优化人类MCM2蛋白的三维结构，以便进行基于结构的药物设计。全长MCM2序列（UniProtKB: P49736）使用AlphaFold 3进行建模，选择AlphaFold 3的原因包括：（i）其在多结构域打包方面的出色性能；（ii）在使用Pairformer Transformer架构预测生物分子结构时，没有实验确定的完整长度的人类MCM2单体结构。
MCM2蛋白的同源性建模
AlphaFold 3（AF3）对人类MCM2单体（UniProt P49736）的预测结果为pTM得分0.63和ipTM得分0.56，这些指标表明整体折叠准确性和结构域间打包信心处于中等至高水平（pTM > 0.5通常表示正确的整体拓扑结构，ipTM > 0.5支持合理的界面排列）。QMEANDisCo全球得分为0.68 ± 0.05，进一步证实了模型的可靠性，因为得分高于0.6表明与距离分布的一致性良好。
局限性、结论与未来工作
局限性：本研究完全基于计算，解释结果时应考虑几个局限性。（i）所使用的受体是预测的完整长度的人类MCM2结构；尽管通过多种信心指标和外部评分进行了验证，但在灵活/无序区域可能存在局部错误，这可能会影响对接姿势和动力学表现。（ii）对接和分子动力学模拟是在MCM2单体及其孤立的ATP口袋模型上进行的；而完整的MCM2–7六聚体和复制体模型并未包含在内。
CRediT作者贡献声明
Seifeldin Elabed：撰写——原始草稿、可视化、方法论、数据整理、概念化。
Wael M Elshemey：撰写——审稿与编辑、可视化。
Ammar Mohamed：撰写——审稿与编辑。
Medhat Shafaa：撰写——审稿与编辑、监督。
未引用的参考文献
25.; 26.
利益冲突声明
作者声明没有已知的竞争性财务利益或个人关系可能影响本文所述的工作。
