24小时热门版块排行榜    

查看: 1755  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 18/900 2026-08-25 14:50 by Weed114114
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 +4 ziyangfang 2026-08-25 4/200 2026-08-25 14:41 by ZJTJZ
[基金申请] 什么时候开奖? +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +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
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +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
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见