smina-盲对接
思路: 把整个蛋白质装进对接盒子,允许在整个盒子中搜索小分子构象
对接中心确认:
pymol中切换-chains;选中整个蛋白质;命令行计算质心:
- centerofmass sele
输出案例如下图:
得到的三个值分别为 smina对接命令里的 --center_x / --center_y / --center_z 后跟的值
对接盒子尺寸确认
有现成的脚本:原链接: DrawGridBox - PyMOLWiki
下载后将脚本“drawgridbox.py” 拖入pymol
- load drawgridbox.py
- # 命令的参数可见链接DrawGridBox - PyMOLWiki;一般改一下 ${your_protein_name}就行,其他参数默认
- drawgridbox ${your_protein_name}, nx=5, ny=5, nz=5, lw=1, g=0, b=0
- drawgridbox ${your_protein_name}, nx=5, ny=5, nz=5, padding=5, lw=1, r=0, g=0
输出案例如下图所示:
脚本会自动生成这个红色的盒子,对smina对接而言,有用的是红框里的三个数值,代表盒子尺寸
即 smina对接命令里的 --size_x / --size_y / --size_z 后跟的值
对接
与常规smina对接无异;只是 --center_x / --size_x 20 (y;z) 后面的值换成前两步得到的值
- # 对接命令案例 smina -h 可查看参数
- smina --seed 1 --num_modes 10 -r 6lk0fix-md.pdb -l 0415.pdb --center_x 48.997 --center_y 44.261 --center_z 37.907 --size_x 20 --size_y 20 --size_z 20 -o dock_0415-t13.pdb