ICM-Pro(MolSoft)是结构建模与药物设计平台,核心算法是内坐标力学(Internal Coordinate Mechanics)——用二面角而非笛卡尔坐标描述分子构象。这个选择带来一个实际优势:自由度大幅减少,构象搜索效率更高。
内坐标方法的意义
- 笛卡尔坐标下,一个 50 原子的分子有 150 个自由度;
- 内坐标下,固定键长键角(它们本来就变化很小),只优化可旋转的二面角,自由度可能只有十几个;
- 搜索空间小几个数量级,全局优化更容易收敛;
- 代价是键长键角被固定,某些需要键长变化的场景(如共价键形成)处理受限。
主要能力
| 功能 | 说明 |
|---|---|
| 柔性对接 | 配体柔性 + 受体侧链柔性 |
| 口袋检测 | ICM Pocket Finder,评估可成药性 |
| 虚拟筛选 | 大规模筛选流程 |
| 同源建模 | 模板搜索与建模 |
| GPCR 建模 | 针对 GPCR 的专门工具 |
| 3D 配体编辑 | 在口袋中交互式修改分子 |
| ICM-Browser | 免费的结构查看器 |
柔性对接为什么重要
绝大多数常规对接把受体当作刚性的,这与现实不符:
- 诱导契合:配体结合时,口袋的侧链会调整朝向以适应配体。用刚性受体对接一个与共晶配体化学型差别较大的分子,很可能因为某个侧链挡路而失败。
- 刚性对接的后果:假阴性——真正能结合的分子因为受体构象不合适而被打低分。这在虚拟筛选中是隐性的损失,因为你永远不知道漏掉了什么。
- 柔性对接的代价:搜索空间急剧增大,计算慢得多,且过度柔性可能引入假阳性(受体扭曲成任何形状去迁就配体)。
处理受体柔性的实际策略
| 策略 | 做法 | 成本 |
|---|---|---|
| 刚性对接 | 用单一结构 | 最低,但有假阴性 |
| 集合对接 | 用多个实验或模拟得到的构象分别对接,取最好结果 | 性价比最高 |
| 软化势 | 降低范德华排斥,容忍轻微重叠 | 低,但可能引入假阳性 |
| 侧链柔性 | 允许口袋侧链旋转 | 中等 |
| 诱导契合对接 | 对接与受体优化迭代 | 高 |
| MD 采样 + 对接 | 用 MD 生成构象集合再对接 | 最高 |
集合对接(ensemble docking)是最实用的折中:用同一靶点的多个晶体结构(或 MD/ENM 生成的构象,见 209《ProDy》)分别做刚性对接,取每个配体的最佳分数。用开源工具就能实现:
# 集合对接的开源实现
for receptor in receptors/*.pdbqt; do
for lig in ligands/*.pdbqt; do
vina --receptor $receptor --ligand $lig \
--center_x 11.2 --center_y 23.4 --center_z -5.7 \
--size_x 22 --size_y 22 --size_z 22 \
--exhaustiveness 16 --seed 42 \
--out out/$(basename $receptor .pdbqt)_$(basename $lig .pdbqt).pdbqt
done
done
# 对每个配体取所有受体构象中的最佳分数
# 同时记录「哪个受体构象给出最佳分数」——这本身是有用信息
选型判断
- 做 GPCR:ICM 在 GPCR 建模上有专门积累,值得评估。
- 需要交互式设计:在口袋里实时编辑分子并看能量变化,对药化很直观。
- 常规小分子对接:Vina + 集合对接的开源方案已能覆盖大部分需求。
- 免费查看器:ICM-Browser 免费,可作为结构查看的备选。
提示
商业软件功能与定价变动较快,本文为方向性介绍,采购前请核对官方最新信息。受体柔性是对接的真实问题,但用开源工具做集合对接已能获得大部分收益。
延伸资源
- 对照平台:275《Schrödinger Maestro》、301《MOE》;开源对接:182《AutoDock Vina》、183《GNINA》;
- 构象采样:209《ProDy》、200《OpenMM》;口袋概念见「结构与模拟」模块。