24小时热门版块排行榜    

查看: 2811  |  回复: 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文件建立的不对?
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 7 个回答

顽灵

新虫 (初入文坛)

为什么没有人理我呢。。金币太少了么?我可以加送金币的哦。。
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的回帖

zhaochuanfan

木虫 (正式写手)

送红花一朵
你好!你的问题最后怎么解决的?我也遇到类似的问题,求赐教!!!谢谢!
知识改变命运,改善自我,才能改善环境
5楼2017-03-14 16:49:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉林大学材料学硕321求调剂 +9 Ymlll 2026-03-18 12/600 2026-03-20 00:01 by 23Postgrad
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 材料专硕英一数二306 +5 z1z2z3879 2026-03-18 5/250 2026-03-19 07:43 by BruceLiu320
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
信息提示
请填处理意见