PyMOL命令

在线下载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目录。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,830评论 5 468
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 83,992评论 2 376
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,875评论 0 331
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,837评论 1 271
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,734评论 5 360
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,091评论 1 277
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,550评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,217评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,368评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,298评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,350评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,027评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,623评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,706评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,940评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,349评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,936评论 2 341