| 查看: 3007 | 回复: 9 | ||
syd89522木虫 (小有名气)
|
[求助]
运行mmpbsa时显示复合物prmtop文件出错
|
|
我计算mmpbsa的时候，总是出现如下错误： [shangyd@bioserver hc]$ $AMBERHOME/bin/MMPBSA.py -O -i mmpbsa.in -o final_results_mmpbsa.dat -sp mod_wat.prmtop -cp rec-lig.prmtop -rp rec.prmtop -lp lig.prmtop -y mod_wat_md2.mdcrd Loading and checking parameter files for compatibility... mmpbsa_py_energy found! Using /home/shangyd/amber14/bin/mmpbsa_py_energy cpptraj found! Using /home/shangyd/amber14/bin/cpptraj Preparing trajectories for simulation... 100 frames were processed by cpptraj for use in calculation. Running calculations on normal system... Beginning GB calculations with /home/shangyd/amber14/bin/mmpbsa_py_energy calculating complex contribution... /home/shangyd/amber14/bin/mmpbsa_py_energy: symbol lookup error: /home/shangyd/amber14/bin/mmpbsa_py_energy: undefined symbol: __libm_error_support File "/home/shangyd/amber14/bin/MMPBSA.py", line 104, in <module> app.run_mmpbsa() File "/home/shangyd/amber14/lib/python2.7/site-packages/MMPBSA_mods/main.py", line 218, in run_mmpbsa self.calc_list.run(rank, self.stdout) File "/home/shangyd/amber14/lib/python2.7/site-packages/MMPBSA_mods/calculation.py", line 82, in run calc.run(rank, stdout=stdout, stderr=stderr) File "/home/shangyd/amber14/lib/python2.7/site-packages/MMPBSA_mods/calculation.py", line 157, in run self.prmtop)) CalcError: /home/shangyd/amber14/bin/mmpbsa_py_energy failed with prmtop rec-lig.prmtop! Exiting. All files have been retained. 已经试过了各种方法，却始终不能正常进行运算，即使用官网教程中的pdb结构来构建prmtop，仍然会出现这种错误，但是用教程现成的prmtop文件来跑mmpbsa，却能够顺利计算出结果，我就是用这个命令：rec-lig= loadpdb "hc.pdb" saveamberparm rec-lig rec-lig.prmtop rec-lig.inpcrd 来生成top文件的 在动力学之前应该如何修改pdb文件呢？我已经按照教程的步骤进行过pdb的前处理了：pdb4amber -i hc.pdb -o hc1.pdb 另外，除了mmpbsa.py、mmpbsa.pl以外，其他涉及到prmtop的模块在计算过程中都没有出现异常。 我用vmd跑出了一段含水的复合物轨迹，一共1000帧只显示出267帧，请问究竟该如何解决这种top发生错误的问题并计算mmpbsa呢？ |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
zhaochuanfan
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 4295.8
- 散金: 1110
- 红花: 17
- 帖子: 765
- 在线: 267.2小时
- 虫号: 1235025
- 注册: 2011-03-16
- 专业: 应用有机化学
2楼2017-03-14 16:52:07
3楼2018-07-06 09:41:32
4楼2019-10-21 18:51:11
syd89522
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 10107.3
- 散金: 72
- 红花: 3
- 帖子: 154
- 在线: 111小时
- 虫号: 778712
- 注册: 2009-05-24
- 性别: GG
- 专业: 理论和计算化学
5楼2020-05-25 14:44:28
syd89522
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 10107.3
- 散金: 72
- 红花: 3
- 帖子: 154
- 在线: 111小时
- 虫号: 778712
- 注册: 2009-05-24
- 性别: GG
- 专业: 理论和计算化学
|
终于知道该怎么解决上述问题了，以前我每次运行前都是把amber重装一遍，运行就正常了，汗~现在是这样： 如下命令利用加水复合物体系生成真空的配体、受体、复合物文件，并在命令后加上后缀 --: $AMBERHOME/bin/ante-MMPBSA.py -p mod_wat.prmtop -s ":WAT,Na+" -c rec-lig.prmtop -r rec.prmtop -l lig.prmtop -m ":1-288" --radii=mbondi2 radii的半径取值取决于in文件中的igb值，一般为2或5时，上述半径值即可，其他值可在教程中查询，有详细说明。 在in文件中，只计算GBSA，不计算PBSA，显然错误就是计算PBSA时带来的，目前结合自由能和熵都已运行完毕，亲测有效！ |
6楼2020-05-28 10:08:25
冷血孤峰
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 299.1
- 帖子: 241
- 在线: 27.6小时
- 虫号: 7226251
- 注册: 2017-09-24
- 性别: GG
- 专业: 生物大分子结构与功能
7楼2020-06-20 09:35:29
8楼2020-07-29 21:17:07
syd89522
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 10107.3
- 散金: 72
- 红花: 3
- 帖子: 154
- 在线: 111小时
- 虫号: 778712
- 注册: 2009-05-24
- 性别: GG
- 专业: 理论和计算化学
9楼2020-10-03 19:55:12
syd89522
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 10107.3
- 散金: 72
- 红花: 3
- 帖子: 154
- 在线: 111小时
- 虫号: 778712
- 注册: 2009-05-24
- 性别: GG
- 专业: 理论和计算化学
10楼2020-10-03 19:57:08