24小时热门版块排行榜    

查看: 4270  |  回复: 9

张筱宇

金虫 (小有名气)

[求助] gromacs模拟分子在二氧化硅表面的吸附

我的体系是用gromacs模拟计算荧光分子在二氧化硅表面的吸附,用ms导出pdb的文件(附件是pdb文件),用pdb2gmx处理时,出现Residue 'O' not found in residue topology database
的错误提示,我该如何修改残基文件呢?谢谢
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : si.pdb
  • 2012-06-14 21:54:10, 65.39 K

» 收录本帖的淘帖专辑推荐

高分子学习--qlh

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
2楼2012-06-14 22:32:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
chaizhm: 金币+3, 谢谢~ 2012-06-15 10:37:03
张筱宇: 金币+20, 有帮助 2012-06-19 20:47:13
pdb2gmx处理时你选择的是哪个力场?
你选择的力场文件中有没有对二氧化硅的定义?

若没有
你需要添加,然后修改一致

若有,你直接把这里面的残基名与力场里修改一致即可。

具体仔细看看gromacs手册吧,手册里讲解的比任何人都详细。
3楼2012-06-15 07:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zh1987hs

金虫 (著名写手)

分子模拟新手


chaizhm: 金币+1, 谢谢~ 2012-06-15 14:12:18
引用回帖:
3楼: Originally posted by jiaoyixiong at 2012-06-15 07:30:53
pdb2gmx处理时你选择的是哪个力场?
你选择的力场文件中有没有对二氧化硅的定义?

若没有
你需要添加,然后修改一致

若有,你直接把这里面的残基名与力场里修改一致即可。

具体仔细看看gromacs手册吧,手 ...

我估计gromacs应该没有对SiO2的定义~不过最新的g53a6有si的参数,应该还是能做出itp文件的
4楼2012-06-15 11:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
chaizhm: 金币+2, 谢谢~ 2012-06-15 19:28:24
张筱宇: 金币+30, 有帮助 2012-06-19 20:41:44
引用回帖:
4楼: Originally posted by zh1987hs at 2012-06-15 11:09:25
我估计gromacs应该没有对SiO2的定义~不过最新的g53a6有si的参数,应该还是能做出itp文件的...

我估计也没有,
但是gromacs 和namd这两个软件方便的地方就是可以根据自己,依据要模拟的物质的力场信息,把没有的力场部分,自己编辑好了之后添加进去。从建模以及参数化来说,比lammps方便一些。
5楼2012-06-15 14:12:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张筱宇

金虫 (小有名气)

引用回帖:
3楼: Originally posted by jiaoyixiong at 2012-06-15 07:30:53
pdb2gmx处理时你选择的是哪个力场?
你选择的力场文件中有没有对二氧化硅的定义?

若没有
你需要添加,然后修改一致

若有,你直接把这里面的残基名与力场里修改一致即可。

具体仔细看看gromacs手册吧,手 ...

对不起,现在才回复您。我用的力场是G53a6,没有对二氧化硅残基的定义,现在就难在这里了,我该如何在力场文件中添加二氧化硅的参数呢,不知道如何修改itp和top,希望您能给建议,谢谢
6楼2012-06-19 20:41:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张筱宇

金虫 (小有名气)

引用回帖:
4楼: Originally posted by zh1987hs at 2012-06-15 11:09:25
我估计gromacs应该没有对SiO2的定义~不过最新的g53a6有si的参数,应该还是能做出itp文件的...

请问我该如何写这个itp文件呢,G53a6力场中,Bond Streching Parameters中Si-OA Force consant Kb 4.72, Ideal bond length 0.163,就是说我该怎样利用这些参数写自己的itp呢,初学gromacs,不是很懂,希望您能给些细节指导,谢谢啦
7楼2012-06-19 20:46:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张筱宇

金虫 (小有名气)

内容已删除
8楼2012-06-19 20:50:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

草莓米粑

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
chaizhm: 金币+2, 谢谢~ 2012-06-20 14:10:30
张筱宇: 金币+30, ★★★很有帮助 2012-06-20 16:25:31
这个没试过,不过对于蛋白质在表面的吸附,表面一般都是固定不动的,所以,这种情况下你要是考虑各种表面内部分子键的话其实没有太大的意义,所以假设表面有3层,你需要做的就是把下面2层固定,只要有二氧化硅的坐标信息即可,然后给每个原子附上LJ参数以及带电量,我看到的文献一般处理第一层的二氧化硅原子会做不固定处理,这时才需要你说的一些键的相互作用信息。然后在.itp文件中加入键的作用信息。至于如何加,我是通过一个GMX的例子学习的,推荐你试试。http://www.bevanlab.biochem.vt.e ... _protein/index.html
9楼2012-06-20 10:51:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjr

铁杆木虫 (正式写手)

你好,我想问一下,看了您发的帖子,用的是MS建二氧化硅晶体,我想问您用的空间群以及空间参数,我在网上查了,在建二氧化硅不知道怎么建不成?
10楼2012-07-05 10:08:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 张筱宇 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见