【软件使用】PyMOL绘图

本篇基于PyMol 2.5.3版本,实测与当前最新版3.0.3基本通用

推荐阅读:PyMOL中文教程PyMOL常用命令

省流:安装完点点GUI,直接看绘图实例

安装

之前的文章【CentOS 7.5】利用conda本地安装蛋白质结构预测工具 里写过简略版,附有AutoDock安装教程,以后详细写AutoDock+PyMOL联用。

教育版可以申请License,见https://pymol.org/edu/

Windows

本方法本质为下载一个conda环境。

  1. 下载 exe installer。官网地址:PyMOL | pymol.org

  2. 安装位置可修改,其他一路默认即可

  3. (实测可不做,注释后依然提示no license file - for evaluation only)找到安装位置\Lib\site-packages\pymol\licensing.py文件,推荐使用everything搜索licensing.py,右键记事本打开,最后三行前面添加“#”注释符号。

  4. 在安装目录下打开PyMOLWin.exe,可以看到此时显示剩余使用时间30天。
    找到C:\Users\(用户名)\.pymoltimestamp文件,推荐使用everything搜索pymoltimestamp,右键记事本打开,数字修改为1945494007(不固定,任一大数字即可)。

  5. 重新打开PyMOLWin.exe即可正常使用。实测本方法修改一次即可,安装新版后正常使用。

Conda安装

推荐新建独立环境后安装,完成后在环境文件夹下找.exe文件。

conda install -c conda-forge -c schrodinger pymol-bundle

常用命令

select

用法:select 变量名(任意的字母或数字组合), 选择的原子

属性选择符有segi(segment),chain(链),resn(残基名),resi(残基号),name(原子名),symbol(元素符号)。逻辑运算符有 and 或者 & 、|(和)、or (或)、not(非)。

select ligand,segi lig # 选择segment为lig的所有原子,赋给变量ligand;
select alian,chain a # 选择A链的所有原子,赋给变量alian;
select aspglu,resn asp+glu # 选择残基名为ASP和GLU的所有原子,赋给aspglu;单个的残基名或残基号间需用+号,不能加空格;
select Y24,resi 24 # 选择残基号24,赋给Y24;
select loop, resi 21-29 # 选择残基号21到29的所有残基,赋给loop;
select Y24-OH,resi 24 and name OH # 选择24号残基中的原子名为OH的原子,赋给变量Y24-OH;
select Y24E27,resi 24 or resi 27 # 选择24号残基和27残基;
select site, not resi 100-200 # 除100到200号残基外的所有残基;

另外还有运算符 around,expand,byres,neighbor:

select coor,name fe around 4 # 选择铁原子周围4埃内的所有原子(不包括铁原子);
select coor,subs expand 4 # 选择变量subs周围4埃内的所有原子(包含subs);
select Y24,byres resi 24 and name OH # 扩充s1到残基; 
select key,neighbor s1 # 选择直接以化学键连接s1的原子,赋给变量key;

show/hide

用法:show 显示形式, 变量或原子

常用显示形式有lines,sticks,cartoon,ribbon,sphere,surface,mesh

show sticks,Y24  # 显示变量为Y24的所有原子为棍棒形;
show sphere, name fe # 显示铁原子为球形;
show cartoon     # 不加操作对象默认操作所有原子;

hide命令操作同show

color/label/align/png

  1. color 颜色, 变量
    color green  # 整个蛋白着色为绿色;
    color red, ss h      # α螺旋着色,ss表示二级结构,h表示helix
    color red, ss s       # β片着色,s表示sheet
    color blue,ss l+"   # loop着色;
    color grey, elem c  # 所有碳原子着色为灰色;
    bg_color white # 背景色
    
  2. label 变量, 标签名
    label resi 24 and name OH,"T24" # 在24号残基的OH原子附近标记为T24;
    label resi 24 and name OH,"" # 可清除标签T24;
    label  # 清除所有标签;
    
  3. align source, target
    Align mol1 & resi N1, mol2 & resi N2 # 按两个残基叠合;
    Align prot1 and chain A, prot2 and chain B # 按两条链叠合;
    Align mol1 & resi n1-n2 & name n+ca+c+o, mol2 & resi n3-n4 & name n+ca+c+o # 按两段残基的主链进行叠合;
    
    无指定的话可直接选择Plugin-Alignment。
  4. png 保存名称
    set ray_shadows,0 # 渲染时不显示投影;
    ray 3000,3000 # 图片像素设置;
    png name  # 输出name.png图片;
    

命令集成

输入@name.pml,就可以执行name.pml文本中的命令。

绘图实例

蛋白-底物相互作用

参考高质量PyMOL作图教程 - 知乎 (zhihu.com)

准备:蛋白与底物结合的pdb文件,确认残基和作用。注意本方法为对已知作用的可视化,不涉及计算预测关键残基。

注意命令中的分隔符。

  1. 导入文件(或右键选择目标文件,直接拖入)

    cd [工作目标]
    load [准备好的pdb文件]
    
  2. 设置配体:选中底物,在<sele>点击A-rename selection,改名lig

    设置关键残基:

    select res, resi 193+191+289+290
    
  3. 显示底物与残基

    hide
    show sticks, lig | res
    hide sticks, h.
    zoom lig | res
    
    # 将lig的碳原子(以C开头的原子名)设为绿色
    color green, lig & name C*
    
    # 设置字体30磅
    set label_size, 30
    # 设置字体为Serif
    # 编号与字体的对应见https://pymolwiki.org/index.php/Label_font_id
    set label_font id, 9
    # 设置label位置,(x,y,z)为离默认位置的三维偏移值; 
    set label_position,(x,y,z) 
    # 建议最后调整label位置,mouse mode改为editing之后右键直接拖动
    

    在res点击L-residues即可添加所有残基标签,氨基酸类型默认三字母缩写

  4. 绘制相互作用:点击PyMOL菜单栏的Wizard -> Measurement,选择作用两端的原子,自动生成虚线和距离对象。

    或者命令行:

    bond atom1,atom2  # 可提前定义变量atom1和atom2,也可以在这里直接输入选择;
    unbond atom1,atom2 # 删除键
    set stick_radius=0.15  # 改变sticks显示半径;
    distance atom1,atom2 # 添加虚键(测量原子间距离)
    
  5. 常用修饰命令

    命令注释
    bg_color white设置背景颜色
    set dash_length, 0.3设置虚线长度
    set dash_radius, 0.08设置虚线粗细
    set dash_gap, 0.3设置虚线间隔
    color gray70, measure09设置作用颜色
    show cartoon color green, !lig & name C* set cartoon transparency, 0.5显示蛋白半透明结构
  6. (可选)使用光线追踪美化图象,生成图片

    save focus.pse
    png label.png
    hide labels
    ray
    png focus.png,16.93cm,16.93cm,dpi=300
    zoom
    set cartoon_transparency, 0
    ray
    png overview.png,16.93cm,16.93cm,dpi=300
    save overview.pse
    
  • 19
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值