24小时热门版块排行榜    

查看: 2143  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们,祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 研究生招生 +3 徐海涛11 2026-03-10 7/350 2026-03-12 14:26 by 徐海涛11
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
信息提示
请填处理意见