24小时热门版块排行榜    

查看: 3125  |  回复: 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呢?
回复此楼

» 猜你喜欢

走自己的路,让别人说去吧~不以物喜,不以己悲
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

向阳奔跑1996

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by 大姚 at 2018-07-06 09:41:32
老师你好,我现在也在跑Amber分子动力学,在运行MMPBSA这一步遇到了和你说的一样的问题,想请教你当时时怎么解决的呢?

请问你的问题解决了吗?我有同样的问题

发自小木虫Android客户端
4楼2019-10-21 18:51:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

zhaochuanfan

木虫 (正式写手)

你好!请问你的问题解决了吗?我也遇到了类似的问题,求赐教!
谢谢!
知识改变命运,改善自我,才能改善环境
2楼2017-03-14 16:52:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大姚

金虫 (小有名气)

老师你好,我现在也在跑Amber分子动力学,在运行MMPBSA这一步遇到了和你说的一样的问题,想请教你当时时怎么解决的呢?
我的地盘听我的!
3楼2018-07-06 09:41:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

syd89522

木虫 (小有名气)

当年我遇到的这个问题,是通过重装amber来解决的,同样的文件,重装以后就能运行的很好,当时我也给amber官网写过信,他们也没有解决。而现在我有遇到了同样的问题,请问谁能帮我正面解答一下?
走自己的路,让别人说去吧~不以物喜,不以己悲
5楼2020-05-25 14:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +15 硕星赴 2026-03-23 15/750 2026-03-26 02:16 by BruceLiu320
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 0856求调剂 +5 zhn03 2026-03-25 6/300 2026-03-25 22:58 by peike
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +7 niko- 2026-03-22 7/350 2026-03-25 20:14 by qingfeng258
[考研] 086000生物与医药292求调剂 +4 小小陈小小 2026-03-22 7/350 2026-03-25 19:07 by 星空星月
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
信息提示
请填处理意见