24小时热门版块排行榜    

查看: 2227  |  回复: 6

顽灵

新虫 (初入文坛)

[求助] AMBER MMPBSA 计算时建立prmtop文件问题

急求帮忙
问题一:
我用tleap建立prmtop文件时,检查电荷如下:
ligand:
Total unperturbed charge:   0.002000
Total perturbed charge:     0.002000
protein:
Total unperturbed charge: -20.000000
Total perturbed charge:   -20.000000
因此,将ligand和protein combine之后的complex电荷为:
Total unperturbed charge: -19.998000
Total perturbed charge:   -19.998000
这样,我用addions给complex加Na+中和,tleap只加19个Na+,为什么?

问题二:
其他体系觉得prmtop文件建立没有问题,但是MMPBSA计算还是出错,过程如下:
Reading command-line arguments and input files...
Loading and checking parameter files for compatibility...
mmpbsa_py_energy found! Using /home/user/amber12/bin/mmpbsa_py_energy
cpptraj found! Using /home/user/amber12/bin/cpptraj
Preparing trajectories for simulation...
51 frames were processed by cpptraj for use in calculation.

Beginning PB calculations with /home/user/amber12/bin/mmpbsa_py_energy
  calculating complex contribution...
CalcError: /home/user/amber12/bin/mmpbsa_py_energy failed with prmtop complex.prmtop!
Exiting. All files have been retained.

这是什么原因呢?还是prmtop文件建立的不对?
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

顽灵

新虫 (初入文坛)

为什么没有人理我呢。。金币太少了么?我可以加送金币的哦。。
2楼2012-09-05 10:13:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluewhale

铁杆木虫 (正式写手)

【答案】应助回帖


jiaoyixiong: 金币+1, 鼓励交流 2012-09-23 22:28:16
可能大家没有amber12,所以没法帮助你。
3楼2012-09-23 10:01:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lujunyan1118

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
jiaoyixiong: 金币+3, 鼓励交流 2012-09-23 22:28:34
第一个问题:
tleap就是这样,只舍不入。 你在家一个 addions Na+ 1 就好了

第二个问题:
我觉得是你的prmtop文件建立的有问题

你用complex的prmtop(没有水的),用ptraj处理一段去掉水的轨迹,比如算个rmsd什么的,看看你这个prmtop是不是好的,或者用vmd载入这个拓扑,打开一段去掉水的轨迹
4楼2012-09-23 22:15:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaochuanfan

木虫 (正式写手)

送红花一朵
你好!你的问题最后怎么解决的?我也遇到类似的问题,求赐教!!!谢谢!
知识改变命运,改善自我,才能改善环境
5楼2017-03-14 16:49:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsbaa1937

禁虫 (初入文坛)

本帖内容被屏蔽

6楼2017-05-01 07:11:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大姚

金虫 (小有名气)

送红花一朵
老师你好!我现在也在跑Amber分子动力学,在计算MMPBSA的时候出现了跟你说的一样的错误,请教你时怎么解决的呢?
我的地盘听我的!
7楼2018-07-06 09:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 顽灵 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年什么时候会评啊 +9 lancet0903 2024-06-24 9/450 2024-06-27 13:22 by yrzhu1984
[有机交流] 大佬们,打的核磁氢谱与chemdraw预测的有偏差 +6 来了个 2024-06-26 7/350 2024-06-27 13:09 by 张小生子
[找工作] offer选择 +10 羡渔樵 2024-06-26 10/500 2024-06-27 12:19 by chenzhikai
[教师之家] 高校辞职,要求赔偿,这到底合不合理 +24 传动_海神 2024-06-23 36/1800 2024-06-27 11:26 by f1yfish
[硕博家园] 考博英语 +6 198新用户 2024-06-25 11/550 2024-06-27 10:13 by Bletilla
[硕博家园] 联培博士文章第一单位署名问题交流~ +10 橙成成c 2024-06-23 23/1150 2024-06-27 10:11 by 橙成成c
[基金申请] 基金评审时,评审专家还回去看申请人代表作之外的文章吗? +10 孤单12站 2024-06-26 15/750 2024-06-27 08:41 by 孤单12站
[论文投稿] 第一篇论文投稿进程记录 +7 慎独的小花卷 2024-06-20 12/600 2024-06-27 00:49 by 玩耍绝缘体
[基金申请] 博后网上办公系统维护 +14 取名字烦人 2024-06-25 22/1100 2024-06-26 22:15 by tuanbao2023
[基金申请] 青年基金E02口青基去年几个函评专家? +6 他山攻玉之石 2024-06-25 9/450 2024-06-26 15:09 by 他山攻玉之石
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +6 babu2015 2024-06-25 6/300 2024-06-26 14:52 by felicity6056
[基金申请] 厅级项目出校却没中 +13 Iwould 2024-06-23 20/1000 2024-06-26 06:14 by foolishmani
[基金申请] 这样的说辞是上会了吗 +12 学员d3zYCz 2024-06-24 15/750 2024-06-25 22:05 by 3115321
[考博] 没读上博,好焦虑! +6 wangzhe_bs 2024-06-24 8/400 2024-06-25 21:15 by wangzhe_bs
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[基金申请] 说博后基金7月出的真打电话了吗? +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 国自然资助比率是不是要下降了?? +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[基金申请] 国自然青年基金,1A4B能上会吗?青年和面上的上会标准是一样的吗? +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[博后之家] 在国内某高校做全职博士后2年,现在找到新的单位,出站或退站对新工作有什么影响? +10 nxplfcc 2024-06-20 10/500 2024-06-22 07:52 by 徐长安
信息提示
请填处理意见