Pymol学习

PyMOL命令

2019-06-17  本文已影响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目录。

上一篇下一篇

猜你喜欢

热点阅读