| 查看: 727 | 回复: 6 | |||
[交流]
【求助】选取the water molecules inside the channels to maintain the structural s 已有2人参与
|
|
这是pdb文件（在附件中）： 在产生psf文件前，我不太清楚怎样选取通道中的水分子，我自己选了一下，大家看一下我选的对不对，期待大家的指点，先谢谢了： 之前我是利用替换矩阵的方法构建了一个四聚体为1YMG-all.pdb mol load pdb 1YMG-all.pdb foreach S { A B C D } { set seg [atomselect top "segname $S and chain A"] $seg writepdb seg$S.pdb $seg delete } foreach S { A B C D } { set wat [atomselect top "segname $S and resname HOH and resid 524 579 417 433 447 448 453 462 469 470 477 480 490 491 492 493 501 507 513 515 518 519 520 522 523 529 531 532 539 542 544 546 551 553 558 559 562 563 564 566 567 569 572 573 576 578 581"] $wat writepdb crystwat$S.pdb $wat delete } 或者： mol load pdb 1YMG-all.pdb foreach S { A B C D } { set seg [atomselect top "segname $S and chain A"] $seg writepdb seg$S.pdb $seg delete } foreach S { A B C D } { set wat [atomselect top "segname $S and resname HOH and within 12.14 of chain A"] $wat writepdb crystwat$S.pdb $wat delete } |
回复此楼» 猜你喜欢
微信指数没变化，科研之友没阅读
已经有4人回复
-
这种情况还有戏吗
已经有9人回复
-
同事接到电话了，我却没有
已经有4人回复
-
准备明年的基金了
已经有5人回复
-
今年的WR进展到哪一步了？
已经有6人回复
-
时间戳他又来了
已经有11人回复
-
小木虫看见有人已经知道结果了
已经有15人回复
-
这个自发加氧反应的机理是什么？
已经有6人回复
-
J-J-W不为人知的一面
已经有18人回复
-
没消息就是被刷了呗
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
2楼2011-02-16 09:23:00
3楼2011-02-16 09:32:57
4楼2011-02-16 11:48:34
5楼2011-02-16 12:35:13
6楼2011-02-16 15:01:25
★ ★ ★ ★ ★
zh1987hs(金币+5): 鼓励以上自己解决问题 2011-02-18 22:51:01
zh1987hs(金币+5): 鼓励以上自己解决问题 2011-02-18 22:51:01
|
已解决，单独生成所需水的pdb文件和蛋白的pdb文件（因为没有金属离子等），然后合并这两个pdb文件即可。 生成所需水的pdb文件小脚本： mol load pdb 1YMG.pdb set model [atomselect top "water and (resid 411 or resid 415 or resid 455 or resid 426 or resid 407 or resid 426)"] $model writepdb wat.pdb exit |
7楼2011-02-16 15:03:48