24小时热门版块排行榜    

查看: 4022  |  回复: 1

gn02530640

银虫 (小有名气)

[求助] 关于Material Studio 的pdb输出档 已有1人参与

先说一下,我想问的问题其实跟分子模拟其实没什么关系…

情况是这样的。我个人目前使用NAMD作模拟动力学计算,我的做法是先在Material Studio 或者是Gaussian 下创建我所需要的单体分子/聚合物以及作前处理(能量优化,计算电荷,etc……),甚至是用Amorphous Tools把分子放进盒子后,再输出成NAMD所需要的pdb文件。总之流程是这样:

用 Gaussian 创建分子并且做结构优化和电荷计算后,输出成pdb文件
-> 用 MS 读取该 pdb文件,用 Amorphous Tools 把分子放进盒子里
-> 将得到的 xtd 文件输出成pdb文件,用于NAMD
(这个过程可能不完全正确甚至有点麻烦,但至少我目前的确是这样做的)

问题来了,我从MS转换输出最后的pdb文件,分子的残基序号都是相同的,举例来说:

HETATM    1  OH2 HOH W   1  
HETATM    2  H1  HOH W   1  
HETATM    3  H2  HOH W   1  
TER       4
HETATM    5  OH2 HOH W   1  
HETATM    6  H1  HOH W   1  
HETATM    7  H2  HOH W   1  
TER       8
HETATM    9  OH2 HOH W   1  
HETATM   10  H1  HOH W   1  
HETATM   11  H2  HOH W   1  
TER      12

可以发现3个水分子的残基序号都是1,这样的pdb文件在产生psf文件的过程中肯定会出错。相反的,如果是在MS下直接画好水分子并且输出成pdb文件的话,就不会有这种情况发生:

ATOM      1  O1  MOL     1
ATOM      2  H2  MOL     1
ATOM      3  H3  MOL     1
ATOM      4  O4  MOL     2
ATOM      5  H5  MOL     2
ATOM      6  H6  MOL     2
ATOM      7  O7  MOL     3
ATOM      8  H8  MOL     3
ATOM      9  H9  MOL     3
TER      10

回到主题,我目前的解决方式是用文本编辑器开启pdb文件后手动修改残基序号,但效率太低…只有几个分子的话还好处理,如果是几千个分子……

因此我的问题很简单:有没有什么办法让MS在输出pdb文件时自动修改残基序号? 或者是有没有什么办法可以快速修改残基序号(比方说写脚本?)

有点太啰嗦了,总之想请各位虫友指点一下迷津,感激不尽!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

qchem

铁杆木虫 (著名写手)


【答案】应助回帖

感谢参与,应助指数 +1
这个可能是不同的软件输出的格式差别
用第三方软件转转?

或者自己写个小程序搞搞

[ 发自小木虫客户端 ]
2楼2014-12-13 18:08:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gn02530640 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 02:25 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 02:23 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 02:20 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 02:12 by 9fOTqj074tIg
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 01:52 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-12 01:50 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-12 01:41 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-12 01:28 by 9fOTqj074tIg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 01:28 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 01:25 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 YHiw1gou1Eda 2026-11-10 4/200 2026-11-12 01:17 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 01:08 by 9fOTqj074tIg
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 01:06 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 01:03 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 00:33 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 4/200 2026-11-11 18:20 by TgvkbZZl9tqB
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-11 16:41 by TgvkbZZl9tqB
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
信息提示
请填处理意见