cnDenis 发表于 2012-11-23 10:46:10

VMD配合APBS显示分子表面静电势

VMD配合APBS显示分子表面静电势

By cnDenis 2012年11月23日

VMD可以调用APBS来计算水溶液中的分子静电势,并进行图形显示。APBS是用来解泊松-波尔兹曼方程(PBE)的软件包。PBE是常用隐式溶剂模型之一,用于计算含盐水溶液中与溶质分子的静电作用。

http://farm9.staticflickr.com/8070/8210332046_0f181bcfac_n.jpg


APBS需要使用PQR或XML格式文件作为输入。PQR格式是在PDB格式的基础上增加了Q(电荷)和R(半径)。PQR格式文件可以通过PDB2PQR网站 http://nbcr-222.ucsd.edu/pdb2pqr_1.8/ 生成。
这个网站可以上传自己的PDB文件,也可以输入PDBID。然后选择一个力场,网站作者推荐用PARSE,据说这是为隐式溶剂优化过的,最适合用于蛋白静电势显示。如果是核酸或是有需要的话也可以选其它力场。其他选项可以参考 http://www.poissonboltzmann.org/pdb2pqr/user-guide 。然后点submit,等待几分钟后就可以下载PQR文件了,也可以点下一步直接让服务器帮你做APBS计算。不过APBS计算速度相当快的,而产生的数据却非常多,把计算结果下载下来的时间可能会比自己机子上完成计算花的时间更长,还不如自己算。

要在本地进行APBS计算,首先是到 http://sourceforge.net/projects/apbs/ 下载,然后安装。注意安装路径不要有空格。

要在VMD中计算APBS,
1. 在VMD中打开刚才生成的PQR文件。
2. Extensions->Analysis->APBS Electrostatics 调出APBS计算界面
3. 在APBS计算界面左上角点Edit->setting,Working Directory 写上计算结果保存的路径,APBS location 写上电脑上APBS程序的位置。确定。[吐槽一下:这每次用之前都要选,是多不人性化的设计啊]
4. 在APBS计算界面点一下Default,程序会自动生成APBS计算的参数。选中0,点Edit查看一下参数对不对,一般默认就可以了。
5. 然后就可以按Run APBS开始计算了。

大约过几十秒后计算完成,弹出一个对话框问把结果放在哪里。选load files into top molecule 就可以把APBS结果装到刚才打开的分子里。这时,VMD主面板上分子的vol项应该为1,表示里面有一组格点数据。

要示显分子表面的静电势分布,
1. 首先在VMD主面板打开Graphics->Representations。
2. 在Representations对话框点Create Rep新建一个图层,在Draw style选项卡中,Drawing Method选 Surf,Coloring Method选Volume,然后在Trajectory选项卡中,Color Scale Data Range 填上-10到10。完成。
如果想改变着色模式,可以在主面板 Graphic->Colors 在Color Scale选项卡中改变Method。

还有其他的静电势显示方法,可以参考官网教程原文。
(本文摘译自http://www.poissonboltzmann.org/apbs/examples/visualization/apbs-electrostatics-in-vmd ,更多信息请参考原文。)

相关链接:
APBS官方网站 http://www.poissonboltzmann.org/apbs
VMD的APBS插件页面 http://www.ks.uiuc.edu/Research/vmd/plugins/apbsrun/


http://farm9.staticflickr.com/8344/8210331506_49493edf0a.jpg
调用APBS插件

http://farm9.staticflickr.com/8067/8209242429_e9eb7c2155.jpg
APBS图形化界面

http://farm9.staticflickr.com/8206/8209242507_ed5b7d6ebd.jpg
APBS计算参数设定

http://farm9.staticflickr.com/8347/8209242533_724099e4ef_m.jpg
APBS计算完成后显示的对话框

http://farm9.staticflickr.com/8485/8210331684_123b657604.jpg
成功载入分子静电势数据后,Vol应该大于0

http://farm9.staticflickr.com/8066/8210331810_edf0732a8c.jpg
显示分子表面静电势

http://farm9.staticflickr.com/8489/8209242835_4f8eb1f380.jpg
改变静电势着色模式

萝卜 发表于 2013-4-20 11:52:59

我在运行apbs后无法产生.dx文件,报错说output files missing or unreadable,我检查过各种路径都没有空格……这是怎么回事?

君木颜 发表于 2012-11-23 11:14:29

学习了,谢啦!!~( ̄▽ ̄)~*

大工-阿里巴巴 发表于 2012-11-23 13:20:47

感谢denis兄的精彩文章,我以前做过一篇文章,使用pymol的APBS插件计算不用pH下蛋白表面电势的变化的~有机会和denis兄切磋一下

tanwen 发表于 2012-11-23 13:22:52

{:soso__10659464829293696901_2:}

墨竹晓风 发表于 2012-11-23 14:01:41

这个非常好,谢谢

川大-灰太狼 发表于 2012-11-23 15:27:27

呵呵 非常经典的教学!支持 Denis兄,学习了。

YAYA_Z 发表于 2014-11-2 20:02:10

正在学习APBS,正好看到你的帖子,太好了,很有用!赞一个!

YAYA_Z 发表于 2014-11-2 20:17:15

请问,用通过PDB2PQR网站生成的PQR文件该怎么下载下来呢

BGIlx 发表于 2019-8-19 17:09:18

萝卜 发表于 2013-4-20 11:52
我在运行apbs后无法产生.dx文件,报错说output files missing or unreadable,我检查过各种路径都没有空格 ...

大哥,问题解决没?我也这样了
页: [1] 2
查看完整版本: VMD配合APBS显示分子表面静电势