24小时热门版块排行榜    

查看: 1758  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 3/150 2026-08-25 23:20 by cNXvBfCpiZOM
[基金申请] 今天务委会开完了,明天出结果吗 +17 angus9576 2026-08-25 21/1050 2026-08-25 22:39 by 悠悠美女
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见