24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2219  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +15 熊二想上岸 2026-04-04 15/750 2026-04-08 22:52 by 猪会飞
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 求调剂 +14 张zic 2026-04-05 15/750 2026-04-08 16:14 by 一只好果子?
[考研] 266求调剂 +18 阳阳哇塞 2026-04-07 18/900 2026-04-08 12:59 by grayjzr
[考研] 求助071001调剂!!! +6 黄守松 2026-04-05 7/350 2026-04-08 11:54 by 猪会飞
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
[考研] 化工调剂303分,过四级 +34 栖梧待风 2026-04-02 34/1700 2026-04-07 12:30 by 1018329917
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 专硕304找调剂,一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 338求调剂,一志愿能源动力,外语是日语203 +5 zzz,,r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
信息提示
请填处理意见