24小时热门版块排行榜    

查看: 1723  |  回复: 5

三木713

新虫 (小有名气)

[求助] 模拟第一步就出了问题,说是残基有问题,求大神帮忙看看是哪出了问题、、急急急!!!

本人是做分子模拟淀粉老化机理的,我先用MS的build建了50个葡萄糖单元的直链淀粉、然后用discover优化、保存为pdb格式文件,内容部分如下:
ATOM      1  C   14A P   1       8.407   4.105  -7.097  1.00  0.00           C   
ATOM      2  HC  14A P   1       9.232   3.388  -7.276  1.00  0.00           H   
ATOM      3  O   14A P   1       7.466   3.570  -6.145  1.00  0.00           O   
ATOM      4  C1  14A P   1       6.966   2.265  -6.409  1.00  0.00           C   
ATOM      5  H12 14A P   1       7.806   1.542  -6.475  1.00  0.00           H   
ATOM      6  C2  14A P   1       6.229   2.282  -7.764  1.00  0.00           C   
ATOM      7  H2  14A P   1       5.337   2.930  -7.649  1.00  0.00           H   
ATOM      8  O3  14A P   1       5.788   0.959  -8.108  1.00  0.00           O   
ATOM      9  HO3 14A P   1       5.534   1.054  -9.044  1.00  0.00           H   
ATOM     10  C3  14A P   1       7.106   2.907  -8.868  1.00  0.00           C   
ATOM     11  H3  14A P   1       7.972   2.233  -9.038  1.00  0.00           H   
ATOM     12  O4  14A P   1       6.309   2.910 -10.059  1.00  0.00           O   
ATOM     13  HO4 14A P   1       6.676   3.579 -10.660  1.00  0.00           H   
ATOM     14  C4  14A P   1       7.685   4.299  -8.475  1.00  0.00           C   
ATOM     15  H4  14A P   1       6.832   4.987  -8.320  1.00  0.00           H   
ATOM     16  O5  14A P   1       8.459   4.714  -9.619  1.00  0.00           O   
ATOM     17  C5  14A P   1       9.009   5.281  -6.278  1.00  0.00           C   
ATOM     18  H51 14A P   1       9.450   4.837  -5.369  1.00  0.00           H   
ATOM     19  O1  14A P   1      10.017   6.092  -6.859  1.00  0.00           O   
ATOM     20  H1  14A P   1       9.520   6.812  -7.277  1.00  0.00           H   
ATOM     21  H52 14A P   1       8.182   5.930  -5.935  1.00  0.00           H   
ATOM     22  O   14A P   1       6.151   1.929  -5.293  1.00  0.00           O   
ATOM     23  C   14A P   1       8.604   6.224 -12.245  1.00  0.00           C   
ATOM     24  HC  14A P   1       8.015   5.293 -12.339  1.00  0.00           H   
ATOM     25  O   14A P   1       9.550   6.092 -11.157  1.00  0.00           O  

然后用gromacs模拟第一步产生top文件,选择的是OPLS-AA全原子力场,在模拟的第一步,就出现了问题、、

residue ‘14A ’ not found in residue topology database.

是我优化出了问题?还是力场选择除了问题?应该怎么解决呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan10

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
月只蓝: 金币+2, 感谢指导! 2013-11-04 13:30:31
三木713: 金币+4, ★★★很有帮助 2013-11-04 14:15:59
对于这种自己构建的模型,要进行动力学模拟的话,是要进行力场参数的构建的。参看这个教程:
http://www.bevanlab.biochem.vt.e ... /complex/index.html
这个就是利用GROMACS软件,对ligand+protein的复合物体系进行动力学模拟的实例。其中对于小分子力场参数构建的方法,可以参考。
努力的生活着
2楼2013-11-04 12:51:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三木713

新虫 (小有名气)

引用回帖:
2楼: Originally posted by wangyan10 at 2013-11-04 12:51:23
对于这种自己构建的模型,要进行动力学模拟的话,是要进行力场参数的构建的。参看这个教程:
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/index.html
这个就是利用GROMACS软 ...

太感谢你了、我先研究研究、不懂的再问你啊
3楼2013-11-04 14:16:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三木713

新虫 (小有名气)

引用回帖:
2楼: Originally posted by wangyan10 at 2013-11-04 12:51:23
对于这种自己构建的模型,要进行动力学模拟的话,是要进行力场参数的构建的。参看这个教程:
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/complex/index.html
这个就是利用GROMACS软 ...

我看了看你发给我的教程,可我做的不是复合物,就只是一条单独的直链淀粉,,是力场不识别的原因么?要怎么解决呢?
4楼2013-11-06 15:02:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fmtzhangli

金虫 (小有名气)

【答案】应助回帖


月只蓝: 金币+1, 感谢指导! 2013-11-10 14:33:43
引用回帖:
4楼: Originally posted by 三木713 at 2013-11-06 15:02:20
我看了看你发给我的教程,可我做的不是复合物,就只是一条单独的直链淀粉,,是力场不识别的原因么?要怎么解决呢?...

去找葡萄糖的力场参数
5楼2013-11-10 13:11:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三木713

新虫 (小有名气)

引用回帖:
5楼: Originally posted by fmtzhangli at 2013-11-10 13:11:53
去找葡萄糖的力场参数...

在哪找呢?应该选什么力场呢?怎么添加进去呢?不好意思啊,是新手所以什么也不懂啊、、、
6楼2013-11-11 08:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 三木713 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐,不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:36 by MWRYsWbJw4DE
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
信息提示
请填处理意见