pymol分子对接

PyMOL 本身是可视化工具而非对接引擎，但在分子对接工作流里它承担着不可替的关键角色：结果检查、结合模式作图与相互作用标注。本项目在每轮对接后都用 PyMOL 做肉眼复核与发表级配图。下文聚焦 PyMOL 在对接中的实际用法，而不是泛泛的软件教程。

为什么对接离不开 PyMOL

对接软件（Vina、Glide、AutoDock）输出的是坐标和分数，但分数排名前三的 pose 未必都合理。本项目习惯把 top 10 pose 全部载入 PyMOL，逐个检查：配体是否卡在蛋白骨架外、是否与金属中心合理配位、关键氢键是否满足几何。这一步能拦掉大量”高分低质”的假阳性。

载入与对齐

# 载入蛋白与对接结果
load receptor.pdb
load poses.pdb
# 按配体名称对齐，避免整体漂移
align poses, receptor

对接结果文件常含多个 MODEL，PyMOL 用 `split_states` 可把每个 pose 拆成独立对象，方便挨个审视。

相互作用可视化

• 氢键：`dist hb, ligand, receptor, mode=2` 可标注距离；本项目偏好手动挑选占有率高的氢键，避免把瞬时接触画进图里。
• 疏水接触：用 `show sticks` + 残基着色（按疏水/极性分色），一眼看出结合由疏水还是极性主导。
• 口袋轮廓：`show surface` 显示蛋白表面，半透明处理让配体位点立体关系更清楚。

发表级配图要点

1. 背景设白（`bg_color white`），避免默认黑底在论文里突兀。
2. 视角用 `get_view` 存下，保证多图视角一致。
3. 配体单独着色（如品红），蛋白用渐变灰，关键残基高亮黄。
4. 导出用 `ray 2400,2400` 光线追踪后再 `png`，分辨率才够印刷。

一个被忽略的坑

PyMOL 的氢键距离判定默认 3.2 Å 且不考虑角度，直接用它标”氢键”会过宽松。本项目做法是：在 PyMOL 里定位候选氢键，再回 MD 轨迹统计占有率交叉验证，只有既满足几何又高占有的才写进图注。

对接流程的衔接

PyMOL 检视合格的 pose，下一步应进入分子动力学短时验证，确认结合不是对接软件的”快照幻觉”。对接软件选型与打分函数差异的更多讨论见本站的分子对接计算栏目。整体科研计算服务可查看科研学术网首页。