PyMol实战：从PDB下载1lEP到绘制靶点-药物相互作用图的保姆级教程
在药物研发和结构生物学领域，可视化分析靶点-药物相互作用是理解分子识别机制的关键环节。PyMol作为一款专业的分子可视化工具，能够帮助研究人员从原子层面解析蛋白质-配体复合物的相互作用细节。本文将以PDB ID为1lEP的晶体结构为例，手把手带你完成从数据获取到高质量相互作用图生成的全流程操作，特别适合刚接触计算机辅助药物设计(CADD)或结构生物学的初学者。
1. 环境准备与数据获取
1.1 PyMol安装与基础配置
PyMol提供开源版和专业版两种版本。对于学术用户，推荐从官网下载开源版本：
# Ubuntu系统安装命令示例 sudo apt-get install pymol首次启动PyMol后，建议进行以下基础设置：
- 界面语言切换为英文（避免命令兼容性问题）
- 将背景色预设为白色（
Display > Background > White） - 调整默认显示质量（
Setting > Quality > 最高）
1.2 PDB数据获取与预处理
1lEP是伊马替尼(Imatinib)与ABL激酶复合物的晶体结构。在PyMol命令行直接输入：
fetch 1lEP该命令会自动从RCSB PDB数据库下载结构文件并加载到工作区。对于网络受限环境，也可手动从PDB官网下载后通过菜单File > Open导入。
提示：使用
fetch命令时若出现连接错误，可尝试添加async=0参数强制同步下载：fetch 1lEP, async=0
2. 结构清理与关键组分保留
2.1 去除冗余结构单元
原始晶体结构通常包含溶剂分子、离子和对称单元等非必要组分。对于1lEP结构：
- 删除冗余蛋白链（本例中B链为对称单元）：
remove chain B - 清理溶剂分子时需特别注意保留关键水桥：
- 通过序列视图(
Show > Sequence)定位14.75和129号水分子 - 选择其余水分子后执行：
remove resn HOH and not (resi 14.75 or resi 129)
- 通过序列视图(
2.2 结构显示优化
执行以下操作提升可视化效果：
| 操作目标 | 命令/操作路径 | 效果说明 |
|---|---|---|
| 蛋白主链 | Action > Preset > Cartoon | 显示α螺旋/β折叠二级结构 |
| 配体分子 | Show > Sticks | 显示配体化学键结构 |
| 关键水分子 | Show > Spheres | 突出水桥的空间位置 |
# 一次性设置命令示例 set cartoon_smooth_loops, 1 # 平滑循环区域 set stick_radius, 0.2 # 调整棍状模型粗细3. 相互作用分析与可视化
3.1 氢键网络识别
- 首先显示潜在相互作用残基：
select interacting_res, byres(ligand around 4) show sticks, interacting_res - 明确氢键相互作用：
distance hbonds, ligand, interacting_res, mode=2 color red, hbonds
注意：
mode=2参数会同时显示直接氢键和水介导的间接相互作用
3.2 关键相互作用位点标注
通过PyMol内置测量工具可量化相互作用距离：
- 启用测量模式：
Wizard > Measurement - 依次点击配体原子和残基原子建立距离标注
- 调整标注样式：
set dash_radius, 0.1 set dash_color, black
对于1lEP结构，应重点关注以下相互作用：
- 吡啶环N与Met318主链NH
- 哌嗪环N与Ile360/His361侧链
- 亚氨基N与Glu286/Thr315侧链
4. 出版级图像渲染
4.1 视角优化与构图
- 使用鼠标操作调整最佳视角：
- 左键拖动：旋转结构
- 右键拖动：缩放结构
- 中键拖动：平移结构
- 通过命令精调视角：
orient ligand # 使配体居中 zoom ligand, 10 # 以配体为中心放大
4.2 高质量渲染设置
执行最终渲染前需配置以下参数：
set ray_shadows, 0 # 关闭阴影避免噪点 set ray_trace_mode, 1 # 启用光线追踪 set antialias, 2 # 抗锯齿级别 set surface_quality, 1 # 表面渲染质量渲染命令及参数：
ray 2400, 2400 # 设置分辨率 png image.png, dpi=300 # 导出高DPI图片4.3 图像后处理建议
虽然PyMol可直接生成高质量图像，但推荐进行以下后期优化：
- 使用Inkscape或Adobe Illustrator添加残基标签
- 用Photoshop微调对比度和色彩平衡
- 组合多视角图像创建拼图
5. 进阶技巧与问题排查
5.1 常见问题解决方案
| 问题现象 | 可能原因 | 解决方法 |
|---|---|---|
| 氢键不显示 | 距离阈值过小 | set h_bond_cutoff_center, 3.5 |
| 表面显示异常 | 计算精度不足 | set surface_algorithm, 1 |
| 渲染速度慢 | 内存不足 | set max_threads, 8 |
5.2 自动化脚本应用
将常用操作保存为脚本可提高工作效率：
# 示例自动化脚本 load 1lEP.pdb remove chain B remove resn HOH and not (resi 14.75 or resi 129) show cartoon show sticks, organic select contact, byres(organic around 3.5) show sticks, contact distance hbond, contact, organic ray 1600,1200 png output.png将此脚本保存为.pml文件后，可通过File > Run Script批量执行。
