ledock盒子的确定
本帖最后由 eming 于 2014-7-28 04:03 编辑Ledock是论坛fireflying博士期间写的一个分子对接软件，其对接准确性高于Autodock Vina. 其中lepro程序提供受体的准备，包括加氢，输出文件仅含坐标。以1BMA为例。
1. 抽取受体。
直接运行 lepro 1bma.pdb
结果可以发现，lepro可以自动去除杂原子，并加氢，残基的名字也改为charmm的命名方式。Arg226在晶体结构中构象不确定，lepro自动保存第一个；
2. 对接的盒子
lepro还产生一个dock.in的对接配置文件：Receptor
pro.pdb
RMSD
1.0
Binding pocket
xmin xmax
ymin ymax
zmin zmax
Number of binding poses
20
Ligands list
ligands
END其中Binding pocket将是确定对接位置的关键信息，其定义方式与普通软件不同的是，一般软件以立方体中心和三边长定义，这里以立方体的一个最小顶点，分别从xyz三个方向延伸，定义xyz三个方向的最大顶点值。这里借助vmd,我们可以确定盒子的大概位置.
需借助一个小的脚本文件boxer.tclproc boxer {args} {
set nargs
puts $nargs
if { $nargs == 1 } {
puts "Generate a box from selection"
set sel
puts "your selection is $sel"
# get the min and max values for each of the directions
set coords ]
set minx
set maxx 0]
set coords ]
set miny
set maxy 0]
set coords ]
set minz
set maxz 0]
} elseif{$nargs == 6} {
set minx
set maxx
set miny
set maxy
set minz
set maxz
} else {
puts "I'm sure your input has some problem, and you will have some errors like:"
}
# and draw the lines
draw materials off
# line in x direction
draw color red
draw line "$minx $maxy $minz" "$maxx $maxy $minz"
draw line "$minx $miny $maxz" "$maxx $miny $maxz"
draw line "$maxx $maxy $maxz" "$minx $maxy $maxz"
draw line "$minx $miny $minz" "$maxx $miny $minz"
# line in y direction
draw color green
draw line "$maxx $miny $minz" "$maxx $maxy $minz"
draw line "$minx $miny $maxz" "$minx $maxy $maxz"
draw line "$maxx $maxy $maxz" "$maxx $miny $maxz"
draw line "$minx $miny $minz" "$minx $maxy $minz"
# line in z direction
draw color blue
draw line "$maxx $miny $minz" "$maxx $miny $maxz"
draw line "$minx $maxy $minz" "$minx $maxy $maxz"
draw line "$maxx $maxy $maxz" "$maxx $maxy $minz"
draw line "$minx $miny $minz" "$minx $miny $maxz"
# print result
puts"xmin and xmax: "
puts"ymin and ymax: "
puts"zmin and zmax: "
}脚本提供两个方法定义顶点，一是根据选择的原子的坐标来定义，原子可以选择某些残基，或者配体周围，二是自定义顶点的位置，例如，我们以晶体中配体极其周围4埃范围内的残基，包括水，离子等，来定义盒子，打开tk consolesource boxer.tcl
setactive_site
boxer $active_site
得到box的信息xmin and xmax: 27.5648.12
ymin and ymax: 13.6331.20
zmin and zmax: 28.0448.87
另外也可以自己定义xmin xmax ymin ymax zmin zmax来输入，看盒子是否合适，并且进行调整draw delete all
boxer 28 50 13 30 28 50这六个数值根据自己定义，来看盒子是否在自己所需要的位置
3.利用pymol的autodock插件定义盒子
这里，我们对autodock插件进行了稍微处理，尽管加了一个Ledock Box的按钮，但其实，跟show box是一样的功能。只要在调整参数的时候，就会自动输出ledock的参数,由于操作比较简单，这里不做详细阐述，只要安装这个插件即可
dengzho5068 发表于 2014-10-25 10:09
我想问一下怎么样在mac osx上实现pymol的插件功能，我的系统是osx x，pymol版本是1.2r1 ...
X11 Hybrid
MacPyMOL comes with a native OSX interface which doesn't require X11 (XQuartz). However, some features, like the Plugin menu, are only available in the X11 interface of PyMOL. There are two ways to launch the X11 interface:
rename /Applications/MacPyMOL.app to /Applications/MacPyMOLX11Hybrid.app
launch the unix executable with the -m flag: /Applications/MacPyMOL.app/Contents/MacOS/MacPyMOL -m 许青良 发表于 2017-3-29 12:29
请问，假如我的对接口袋是一个长方体形状，而且是水平方向的长方体，在vmd中执行完source boxer.tcl seta ...
在vmd中生成的box形状是左图，在pymol中的box是右图中的方框中的区域，对接的口袋类似一条“沟”，对接的分子在里面伸展开来。假如能把vmd中的box顺时针旋转一定的角度，使之能和“沟”相匹配，后期的对接效果应该能不错。如果单纯把box大小扩大，使之能够将“沟”给包围起来，感觉box的改变就太大了，不太好。怎样才能解决这个问题？
请问，假如我的对接口袋是一个长方体形状，而且是水平方向的长方体，在vmd中执行完source boxer.tcl setactive_site 以后出现的box是一个倾斜的box,请问能否对box进行旋转，使其能够和受体本身的口袋方向相一致。 呵呵，飞天出手果真不一样！刚动手试了一下，挺好用的。:lol
我还是第一次知道可以用VMD来显示一个box，然后来设置结合位点。以前都是用autodock的adt软件设置结合口袋，不是特别方便，而且adt的图形化效果比VMD要差很多。估计这个脚本修改一下也可以用来设置Autodock或者vina的口袋吧。
要是输出结果小数点只保留一位就好了，那么多数字看得眼花啊。:lol
学习了，谢谢分享 我想问一下怎么样在mac osx上实现pymol的插件功能，我的系统是osx x，pymol版本是1.2r1 dengzho5068 发表于 2014-10-25 10:09
我想问一下怎么样在mac osx上实现pymol的插件功能，我的系统是osx x，pymol版本是1.2r1 ...
好像pymol在mac osx上不支持插件。有解决办法，网上找一下吧 跟师兄学习了，谢谢~ 你好，我在lephar research官网上面下载不了，能不能提供下载软件包？ 谢谢分享！ 安装了基于pymol的插件后，调整盒子大小就报错。怎么回事？请教解决办法
页:
[1]
2