| 查看: 746 | 回复: 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 } |
回复此楼» 猜你喜欢
今天维护系统维护 祝所有人 高中
已经有7人回复
-
读博的好处
已经有5人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有98人回复
-
2027广东省杰青
已经有8人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有5人回复
-
投稿咨询
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有6人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
