24小时热门版块排行榜    

查看: 3452  |  回复: 5

mojie1987

铜虫 (小有名气)

[求助] amber的tleap出错

有一个小配体是GLP,tleap之后,出错啦。。

Done adding ions.
> savepdb MODEL model.pdb
Writing pdb file: model.pdb
Converting N-terminal residue name to PDB format: NGLY -> GLY
Converting C-terminal residue name to PDB format: CGLU -> GLU
> saveamberparm MODEL model.prmtop model.inpcrd
Checking Unit.
WARNING: There is a bond of 3.201702 angstroms between:
-------  .R.A and .R.A
WARNING: There is a bond of 4.241588 angstroms between:
-------  .R.A and .R.A
FATAL:  Atom .R.A does not have a type.
FATAL:  Atom .R.A

does not have a type.
FATAL:  Atom .R.A does not have a type.
FATAL:  Atom .R.A does not have a type.
FATAL:  Atom .R.A does not have a type.
Failed to generate parameters
Parameter file was not saved.
> quit
        Quit
在tleap之前,蛋白质大分子已经除氢了。。。那这个该肿么解决呢?

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

yyuan8658

金虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
mojie1987: 金币+1, 有帮助 2013-04-10 19:26:25
GLP609中的几个原子没有力场参数
搜一下以前的关于做力场参数的帖子吧,
2楼2013-04-10 07:23:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan10

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
mojie1987: 金币+3, ★★★★★最佳答案, 那个帖子很好哈。。。已经解决啦。。谢谢啦。。。 2013-04-10 19:27:07
你这个用tleap处理的蛋白里面有amber力场不识别的原子类型。可以参看amber的手册进行操作。
http://ambermd.org/tutorials/basic/tutorial4b/
另外,如果用amber自带的量化计算模块不能给出正确的计算结果,可以结合使用gaussian。论坛里面有amber+gaussian构建小分子力场的方法,可以参看一下。
http://muchong.com/bbs/viewthread.php?tid=3052985
努力的生活着
3楼2013-04-10 08:13:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaozhq

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
mojie1987: 金币+1, 有帮助 2013-04-10 19:27:32
楼上说的很明白了 你没有给GLP设定力场参数 需要你用antechamebr或者gaussian来处理一下
4楼2013-04-10 10:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangdafa

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
antechamber处理未定义的residue,如果分配resp电荷就gaussian生成out文件导入antechamber,如果是用bcc,就之接用pdb的结构文件,希望可以帮到你
5楼2013-04-10 11:15:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mojie1987

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by jiangdafa at 2013-04-10 11:15:44
antechamber处理未定义的residue,如果分配resp电荷就gaussian生成out文件导入antechamber,如果是用bcc,就之接用pdb的结构文件,希望可以帮到你

金币木有啦。。。不过很感谢哈!~~~
6楼2013-04-10 19:28:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mojie1987 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 01:03 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:41 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见