CML 双抑制剂分子动力学模拟——项目原理说明

研究问题

慢性粒细胞白血病（CML）的驱动因子是 BCR-ABL1 融合激酶。临床上已有两类抑制剂：

ATP 竞争性抑制剂（如奥雷巴替尼/Olverembatinib）：占据激酶催化口袋，阻断磷酸化功能
别构抑制剂（如阿思尼布/Asciminib）：结合肉豆蔻酰口袋，锁定激酶非活性构象

两药靶向不同位点，理论上可同时结合形成”双重锁定”。本研究旨在从分子层面验证：双药联合是否产生协同增效？

技术路线

晶体结构 (PDB: 5MO4, 2.5Å)
    │
    ├─ 蛋白准备 (pdbfixer修复缺失残基 → AMBER99SB-ILDN力场)
    │
    ├─ 配体参数化
    │   ├─ 奥雷巴替尼: AM1-BCC电荷 + GAFF2力场
    │   └─ 阿思尼布: Gasteiger电荷 + GAFF2力场
    │
    ├─ 分子对接 (AutoDock Vina)
    │   └─ 奥雷巴替尼 → ATP位点 (最佳打分: -8.6 kcal/mol)
    │
    ├─ 体系构建 (3组 × 3重复 = 9个独立模拟)
    │   ├─ olv_only:      蛋白 + 奥雷巴替尼 + 水 + 离子
    │   ├─ asc_only:      蛋白 + 阿思尼布 + 水 + 离子
    │   └─ dual_complex:  蛋白 + 双药 + 水 + 离子
    │
    └─ 分子动力学模拟 (GROMACS 2023, GPU加速)
        ├─ 能量最小化 → NVT平衡 → NPT平衡
        └─ 生产模拟: 200 ns × 9 runs (总计 1.8 μs)

核心原理——为什么用分子动力学？

分子动力学（MD）的本质是：给每个原子赋予初始速度，按牛顿力学求解运动方程，观察蛋白-药物体系在生理条件（310K、1 atm、显式水环境）下的动态行为。

时间步长：2 飞秒（10⁻¹⁵ s）
模拟时长：200 纳秒（10⁸ 步）
体系规模：约 79,500 个原子

这相当于给蛋白-药物复合物拍了一部”原子级电影”，能观察到： 1. 药物是否稳定待在结合口袋？还是逐渐脱离？ 2. 蛋白结构是否因药物结合而变得更刚性？ 3. 两药共存时，蛋白构象是否额外稳定？

分析方法与预期产出

分析指标	物理含义	预期结论
RMSD	蛋白骨架偏离初始结构的程度	双药体系 < 单药 → 蛋白被锁得更紧
RMSF	每个残基的运动幅度	活性位点残基在双药时波动更小
氢键分析	药物-蛋白间氢键数量与占有率	双药产生新的氢键网络
接触分析	药物与哪些残基形成稳定接触	识别协同作用的关键残基
MM-PBSA	结合自由能（ΔG）	ΔG(dual) < ΔG(olv) + ΔG(asc) → 热力学协同

实验设计的严谨性

3个独立重复（不同随机种子生成初始速度）：排除偶然性
200 ns 时长：确保体系充分平衡（最后 50 ns 用于自由能计算）
显式溶剂 + 反离子：模拟真实生理环境
AMBER99SB-ILDN + GAFF2：蛋白和小分子领域广泛验证的力场组合

当前进度

体系构建与平衡：✅ 已完成
9个 200 ns 生产模拟：🔄 运行中（V100 GPU，预计 7月16日前完成）
后分析与出稿：⏳ 待 MD 完成后执行

研究定位

本项目为临床-计算联合研究的计算部分：

临床端：教授团队在临床实践中已观察到奥雷巴替尼+阿思尼布联合用药的协同现象
计算端（本项目）：通过分子动力学模拟阐释协同现象的分子机制——为什么两药同时结合能产生”1+1>2”的效果
体外验证：补充细胞水平实验作为桥梁证据

论文叙事逻辑：临床发现 (What) → 计算解释机制 (Why) → 体外验证 (Confirm)

预期学术价值

首次从分子动力学层面阐明奥雷巴替尼+阿思尼布双靶点协同机制
为临床已观察到的协同现象提供原子分辨率的结构解释
方法可推广至其他激酶双靶点抑制剂组合筛选

技术环境：GROMACS 2023.2 (CUDA) / AutoDock Vina 1.2.5 / AmberTools / RDKit
计算资源：NVIDIA Tesla V100-SXM2-16GB