ZDOCK-进行蛋白-蛋白分子对接

1. 软件下载

需要教育邮箱,下载连接https://zdock.umassmed.edu/software/

tar -zxvf zdock3.0.2_linux_x64.tar.gz
# 解压即用,已下载在~/programm/zdock/目录中
# ZDOCK 3.0.2版本是最新版本
# 将 export PATH=/home/shims/install/zdock3.0.2_linux_x64:$PATH添加到环境变量中。

2. 准备文件

从pdb数据库下载目标蛋白的.pdb格式的文件,并进行预处理

#  预处理
mark_sur receptor.pdb receptor_m.pdb
mark_sur ligand.pdb ligand_m.pdb
#利用mark_sur命令将下载的pdb文件处理成zdock可使用的文件
# uniCHARMM文件需在当前目录下,可从~/programm/zdock/zdock3.0.2_linux_x64/拷贝

3.分子对接

[命令行运行]https://zlab.umassmed.edu/zdock/how_to_run_zdock.shtml

 zdock -R receptor_m.pdb -L ligand_m.pdb -o zdock.out # 输出结果在当前目录的zdock.out文件中
 head -15 zdock.out > top10.zdock.out  #取前10个进行后续分析
 create.pl top10.zdock.out num_preds 
# 预测前十个complex, create_lig文件需在当前目录下。num_preds参数指定预测前多少个
 # -o是输出文件夹

4. 导入pymol进行可视化

上述输出的complex可利用pymol进行可视化
[pymol教学]https://www.bilibili.com/video/BV1fa4y1s7e6/?spm_id_from=333.788.recommend_more_video.1&vd_source=2530fee38a2e3828d07aaa316ccd75b9
利用pymol打开zdock对接的complex,进行可视化,同时打开 InterfaceResidues的python文件可以找到蛋白-蛋白对接的界面,视频教程https://www.bilibili.com/video/BV1tZ4y137VB/?spm_id_from=333.337.search-card.all.click&vd_source=2530fee38a2e3828d07aaa316ccd75b9 #命令 interfaceResidue 1hwg, chain A, chain B

结果解读

ZDOCK分子对接之后的输出文件解读
ZDOCK分子对接之后会将所有的信息打印到用户指定的文件中,如通过命令行:zdock -R receptor_m.pdb -L ligand_m.pdb -o zdock.out就会将所有的打印输出信息到文件zdock.out。输出文件通过zdock命令的-o参数设定。

128 1.2 0

-3.141593 1.450119 -2.739050

0.000000 0.000000 0.000000

receptor_m.pdb -38.481 1.151 -26.315

ligand_m.pdb -27.379 38.575 4.226

第一行是ZDOCK在进行分子对接过程中使用的最大NxNxN格点数目(128)。空间格点之间的空间间距(1.2Å),一般情况下都是1.2Å的值。接下来就是在zdock对接过程中将保持刚性的分子,0表示受体分子刚性,1表示配体分子刚性。

第二行是受体分子的初始Euler角度旋转参数。

第三行是配体分子的初始Euler角度旋转参数。该参数可以通过zdock命令的-S参数来设定。

第四行是受体分子文件名称,可以根据zdock命令的-R参数来设定。接着的是受体分子相对于中心来的位置参数。

第五行是配体分子的文件名称,可以根据zdock命令的-L参数来设定。接着是配体分子相对于中心的位置参数。

在前面的五行头文件信息之后的每一行都是配体分子相对于初始位置的旋转和平移位置。

0.261799 1.506530 1.988858 120 16 6 1101.029

每行的前三个数字是配体分子相对于初始位置的Euler角度的旋转。

接下来的三个数字是配体分子相对于初始位置的平移情况,使用空间格点空间来描述。

最后一列数字是当前结构的ZDOCK打分。这个打分一般情况下是从最高分到最低分排序的。

在ZDOCK的结果文件中没有具体的复合物的结构,当使用creat.pl脚本提取复合物结构的时候,会从上一次提取结构,如第六行对应第一个复合物结构,第七行对于第二个复合物结构等。每个复合物的结构都根据Euler角度旋转和格点空间的平移来构建复合物结构。
参考连接https://blog.sina.com.cn/s/blog_15d0344a30102xeh8.html

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

推荐阅读更多精彩内容