24小时热门版块排行榜    

查看: 1051  |  回复: 3

林祎

木虫 (正式写手)

[求助] psf 不生成C-H键

我想将蛋白子的PDB利用VMD生成PSF文件,我想不生成类似C-H键,图1为我生成的,图2为目标,显然图2中没有C-H这类的氢原子的干扰,
实在没辙了,请求高人指点,谢谢

1.png



2.png
回复此楼

» 猜你喜欢

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

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

qianxiong

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2013-03-06 08:35:59
林祎: 金币+5, 我想到办法了,就是生成psf文件后在选择原子导出 2013-03-06 12:22:13
林祎: 金币+2, ★★★很有帮助 2013-03-06 12:24:05
psf 是根据top文件(default 全原子的)中各个氨基酸的结构定义而来,psfgen 根据 top将通常原始的PDB(x-ray) 中缺少的原子补全。所以即使你要生成 不含H,只能去修改top. 可以参考 charmm 力场官网中 top 压缩包 中notcharmm 文件夹,其中包括使用联合原子的gromos的top.
戒微博
2楼2013-03-06 08:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
林祎: 金币+5, 回复同上,就是生成psf文件后在选择原子导出 2013-03-06 12:23:46
林祎: 金币+2, ★★★很有帮助 2013-03-06 12:24:14
引用回帖:
2楼: Originally posted by qianxiong at 2013-03-06 08:23:31
psf 是根据top文件(default 全原子的)中各个氨基酸的结构定义而来,psfgen 根据 top将通常原始的PDB(x-ray) 中缺少的原子补全。所以即使你要生成 不含H,只能去修改top. 可以参考 charmm 力场官网中 top 压缩包 ...

我有个比较愚蠢的办法,VMD 的 InorganicBuilder 的插件中,Task下面有个Add bonds
利用这个功能可以实现,不添加与H 有关的键。
里面可以选择根据原子类型来添加键,如果你知道其他你需要添加键的原子类型,把他们都列出来即可。
3楼2013-03-06 08:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qianxiong

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by qianxiong at 2013-03-06 08:23:31
psf 是根据top文件(default 全原子的)中各个氨基酸的结构定义而来,psfgen 根据 top将通常原始的PDB(x-ray) 中缺少的原子补全。所以即使你要生成 不含H,只能去修改top. 可以参考 charmm 力场官网中 top 压缩包 ...

补充一下:namd 常用的 charmm27/31/36 都是全原子的。而gromos的联合原子模型就不含非极性氢,并且 “charmm 力场官网中 top 压缩包 中notcharmm 文件夹” 已是 charmm的书写格式,可以直接用
戒微博
4楼2013-03-06 08:51:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 林祎 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见