24小时热门版块排行榜    

查看: 2041  |  回复: 11

YAYA_Z

新虫 (初入文坛)

[求助] 将糖链的PDB生成PSF文件 已有1人参与

我是刚开始接触用NAMD动力学分析这块,现在需要把一个糖链生成psf文件,可是生成后放到VMD里看有问题,各个残基断开了,就是不是一条链了,不是只加氢吗,为什么也多加了氧把链拆开了?就是糖生成psf需要注意什么呢?谢谢高手解答我的疑问
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gn02530640

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
YAYA_Z: 金币+5, ★★★很有帮助, 感谢帮助 2014-10-27 14:41:10
完成了,请楼主先看一下我做出来的文件是否可行。

压缩文件中的 tang_mod.pdb 和 top_all36_carb_mod.rtf 是我将原本的 tang.pdb 还有 top_all36_carb.rtf 部分内容修改过后的档案,而 tang_new.pdb 和 tang_new.psf 则是用 psfgen 所生成的档案。

关于残基断开的问题:
在利用 psfgen 生成 psf 文件时,必须注意到,pdb 的 atom name 和 residue name 必须和 top 文件里的定义完全相同。
举例来说,必须将 pdb 文件中的 BMA 残基名称改成 top 文件所定义的 BMAN,
此外,pdb 文件中 NAG 残基的 C7 C8 N2 O7 这几个原子在 top 文件中 BGLCNA 的命名分别是 C CT N O,
因此也有必要将 atom name 作修改。

另外,由于 pdb 的 residue name 最长似乎只能用4个字母代表,因此我将top文件 BGLCNA 残基的名称改成 BGL,基本上只要注意别跟其他残基名称重复就行了。

大致上就这样

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : tang.rar
  • 2014-10-27 13:39:50, 72.95 K
4楼2014-10-27 14:13:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gn02530640

银虫 (小有名气)

引用回帖:
7楼: Originally posted by 916756502 at 2014-10-27 14:50:32
我和发帖的是同学,我俩一个研究组的,我刚才看了你生成的pdb和psf,与我之前生成的相比NAG那部分正常了,可是你放大了一看,在两个环相连的氧原子部分还是分开的,也就是说糖链还是断开的...

我这看到的是连起来的没错,不过会有一部分奇怪的键结。

补氢的话可以用MS或者其他软件补上,之后再将pdb中氢原子的atom name改成top文件中定义的形式(EX:HO61 HO62 etc….),也许就能解决(糖链断开)这个问题。
8楼2014-10-28 21:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gn02530640

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
Psf?n?????Y???Y?????]???k??????@???????????????psfgen ????a????pdb?n?

?????韩??pdb?????top??????P???朵???x?Y???N????????????NAMD?????????????????
2楼2014-10-25 10:05:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

YAYA_Z

新虫 (初入文坛)

???????:
2?: Originally posted by gn02530640 at 2014-10-25 10:05:59
Psf?n?????Y???Y?????]???k??????@???????????????psfgen ????a????pdb?n?

?????韩??pdb?????top??????P???朵???x?Y???N????????????NAMD?????????????????

?????????????л????????top????????л??????г???
BMA---BMAN
MAN---AMAN
NAG?????????BGLCNA?????top??????кü?????????塣??????????NAG??????6?????????????????????????5????????????????????????????????????????????
???????????????????????psf??????????pdb???

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : top_all36_carb.rtf
  • 2014-10-27 08:52:06, 274.31 K
  • 附件 2 : tang.pdb
  • 2014-10-27 08:52:23, 8.36 K
3楼2014-10-27 08:56:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gn02530640

银虫 (小有名气)

但这里有个问题(用VMD打开tang_new.pdb应该就有发现),利用 psfgen 自动补氢有时候会产生不合理的键结,因此我个人并不建议用psfgen自动补氢,最好的方法还是一开始就把所有的原子(包括氢原子)都写进pdb文件后,才用psfgen生成psf文件。
5楼2014-10-27 14:37:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

916756502

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by gn02530640 at 2014-10-27 14:37:36
但这里有个问题(用VMD打开tang_new.pdb应该就有发现),利用 psfgen 自动补氢有时候会产生不合理的键结,因此我个人并不建议用psfgen自动补氢,最好的方法还是一开始就把所有的原子(包括氢原子)都写进pdb文件后,才用 ...

是在pdb里面自己把氢用键盘敲上去吗?
6楼2014-10-27 14:41:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

916756502

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by gn02530640 at 2014-10-27 14:37:36
但这里有个问题(用VMD打开tang_new.pdb应该就有发现),利用 psfgen 自动补氢有时候会产生不合理的键结,因此我个人并不建议用psfgen自动补氢,最好的方法还是一开始就把所有的原子(包括氢原子)都写进pdb文件后,才用 ...

我和发帖的是同学,我俩一个研究组的,我刚才看了你生成的pdb和psf,与我之前生成的相比NAG那部分正常了,可是你放大了一看,在两个环相连的氧原子部分还是分开的,也就是说糖链还是断开的
7楼2014-10-27 14:50:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

916756502

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by gn02530640 at 2014-10-28 21:29:22
我这看到的是连起来的没错,不过会有一部分奇怪的键结。

补氢的话可以用MS或者其他软件补上,之后再将pdb中氢原子的atom name改成top文件中定义的形式(EX:HO61 HO62 etc….),也许就能解决(糖链断开)这个问题。...

你说为什么会有这些奇怪的键结呢?你平时用AUTODOCK这个软件不
9楼2014-10-29 08:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

916756502

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by gn02530640 at 2014-10-28 21:29:22
我这看到的是连起来的没错,不过会有一部分奇怪的键结。

补氢的话可以用MS或者其他软件补上,之后再将pdb中氢原子的atom name改成top文件中定义的形式(EX:HO61 HO62 etc….),也许就能解决(糖链断开)这个问题。...

生成psf过程中产生了一个pdb和一个psf文件,只在VMD里打开pdb会出现你说的奇怪键结的现象,不过你先打开psf再导入pdb,就会看到分子之间是断开的了。我还在纠结这个断开问题,这个问题不解决我没法继续往下做
10楼2014-10-29 15:28:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 YAYA_Z 的主题更新
信息提示
请填处理意见