24小时热门版块排行榜    

查看: 2444  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 5/250 2026-08-04 04:44 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 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
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见