24小时热门版块排行榜    

查看: 1774  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 16:25 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:57 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见