PyMOL命令

作者: 佳名 | 来源:发表于2019-06-17 12:53 被阅读0次

    在线下载pdb文件:

    先在NCBI子数据库structure检索所需要的蛋白结构,https://www.ncbi.nlm.nih.gov/structure/?term=
    记录下该蛋白的PDB ID,然后打开PyMOL,命令行输入:

    fetch 5ocn#foxn1
    fetch 3uf0#
    fetch 1si4#血红蛋白
    

    5ocn和3uf0为PDB ID


    1.png

    去除水分子

    remove solvent
    

    分离得到蛋白

    remove organic
    
    2.png

    参考自:
    https://www.jianshu.com/p/ed3446ebd2e8
    https://pymol.org/pymol-command-ref.html

    切换目录

    cd D:\demo
    

    导入文件

    load xxx.pdb
    

    隐藏所有对象

    hide
    

    鼠标的几个操作用法:

    A) 拖动鼠标左键:按住左键,任意方向拖动鼠标,使视图转动;

    B) 上下拖动鼠标右键:按住右键,上下拖动鼠标,使视图缩放;

    C) 拖动鼠标中间滚轮:按住滚轮,任意方向拖动鼠标,使视图移动;

    D)滚动鼠标中间滚轮:拉进或拉远镜头,使近处原子逐渐显示或隐藏。

    下载

    纤维素外切酶CBHI和纤维素糖链的复合物晶体结构

    fetch 7cel
    

    螺旋飘带模型(α螺旋和β折叠),在PyMOL中叫cartoon模式

    as cartoon
    

    命令显示多种模式

    show lines 
    show sticks 
    show ribbon
    

    那么怎么用命令实现着色呢?比如整个蛋白显示绿色,命令是

    color green
    

    对α螺旋、β折叠和无规则卷曲着不同的颜色怎么实现呢? 比如α螺旋(helix)显示红色,β折叠(sheet)显示绿色,无规则卷曲(loop)
    显示蓝色,命令是

    color red,ss h 
    color green,ss s 
    color blue,ss l+’’
    

    ss h,ss s的意思了,其实就是选择二级结构α-helix或β-sheet,那么ss l+”” (这里单引号和双引号都行)为什么还要写 +”” 呢? 因为二级结构的分类不仅仅只有这三种,PyMOL为了简化,就把不是α-helix和β-sheet的结构全部归为loop
    和其他无规则结构了,所以这里得用 l+”” 表示。
    文献上的结构图一般都是白色背景,怎么设置背景颜色呢? 命令是

    bg_color white
    

    比如选择7cel中的一个催化残基212位的谷氨酸GLU,命令如下

    select geng,resi 212
    

    点击geng框后S按键中的sticks,212GLU残基就会显示出sticks模式,再点
    击C按键中的oranges改为桔色,如下图。这一切也可用命令实现

    show sticks,geng 
    color orange,geng
    

    点击geng小框L按键中的residues,这样就对212GLU加上了标签,L就是Label的意思。
    突出显示

    zoom geng
    

    点击上图右下角的S按键,S指sequence,点击后会显示出蛋白序列。刚才选择残基是用命令实现的,也可以用鼠标左键在蛋白结构或蛋白序列上点击,点击中的残基就会被选中并起名为sele,然后使用A按键中的rename selection修改名称即可,你会发现被点中的残基都会显示出粉红色的小点
    击上图右下角的S按键,S指sequence,点击后会显示出蛋白序列。刚才选择残基是用命令实现的,也可以用鼠标左键在蛋白结构或蛋白序列上点击,点击中的残基就会被选中并起名为sele,然后使用A按键中的rename selection修改名称即可,你会发现被点中的残基都会显示出粉红色的小点。 可是到目前为止,
    所选择的最小单位都是残基,如果只想选择一个原子或者一条肽链,
    怎么做?鼠标操作的话,首先要修改选择的单位,仔细看一下右下角是不是有 Selecting Residues,点击Resdiues,看看都有什么。如果想选择单个原子,那么点击到显示Atoms的时候停下来,如图 然后再在结构上点击某个原子,就选中它了。此时不能在序列上点击了,因为序列只有残基名没有原子名。 举一反三,你也可以选择一个分子、一条链、一个Cα原子等等
    .3 布尔算符和标签命令
    中学时学过布尔算符and/or/not,在PyMOL的选择命令中布尔算符非常有用。比如选择7cel中212GLU残基的Cα原子,怎么选? select geng,resi 212 和select geng,name ca两个命令显然不行,因为它们分别对212
    残基的所有原子以及蛋白的所有Cα原子进行了选择,而不只是212GLU的Cα原子。 可以这样做

    select geng,resi 212 and name C
    

    那么如何选择212和214号残基的所有原子呢?

    select geng,resi 212 or resi 214
    

    选择除200-400号残基外的所有残基

    select geng,not resi 200-400  color white,geng 
    

    显示一下被选中的残基
    选择212和214号残基的非Cα原子

    select geng,resi 212+214 and not name ca hide all 
    show sticks,geng label geng,name
    

    结构打光并保存图片

    ray  
    png 001  
    

    打光命令ray就像在照相馆照相时需要专门的灯光设备,这里ray打光后,图像显示出景深和立体效果,更加逼真和精美。 保存图片命令非常简单,png后面加上你起的任意文件名即可,回车后图片自动保存在当前工作路径下,图片文件格式是png,看看后缀就清楚了。注意当前工作路径在哪里,还记得pwd目录。

    相关文章

      网友评论

        本文标题:PyMOL命令

        本文链接:https://www.haomeiwen.com/subject/pyhafctx.html