24小时热门版块排行榜    

查看: 2155  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 4/200 2026-03-16 21:52 by Losir
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
信息提示
请填处理意见