24小时热门版块排行榜    

查看: 574  |  回复: 6

nufang19a

金虫 (正式写手)


[交流] 【求助】选取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
}
回复此楼

» 猜你喜欢

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

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

qphll

金虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+2): 谢谢 2011-02-16 13:03:21
引用回帖:
Originally posted by nufang19a at 2011-02-16 08:26:23:
这是pdb文件(在附件中):
在产生psf文件前,我不太清楚怎样选取通道中的水分子,我自己选了一下,大家看一下我选的对不对,期待大家的指点,先谢谢了:
之前我是利用替换矩阵的方法构建了一个四聚体为1YMG-al ...

没有明白你要做什么.

你是要将结构(你的PDB)中的水分子做一些筛选, 只留下在通道中的水分子嘛?
Life, Love, Laugh.
2楼2011-02-16 09:23:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


刚才一位高人说,通道中的水分子应该在生成四聚体前保留的,而不是在生成四聚体后。看他操作了一下,自己还在没看懂时就结束了。他时间很紧,没好意思再接着问。郁闷
3楼2011-02-16 09:32:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


是的,现在我用vmd已选择好通道中的水分子,采用的是最古老的方法。就是在vmd中直接看的,以供选择了6个水分子
4楼2011-02-16 11:48:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


采用最古老的方法已经选择出在通道中的水分子,他们的编号是411 415 455 426 407 436 ,可生成还有这些水分子的pdb时,却无论如何不能生成,大家看看我的小脚本写的哪里对不对好吗??
#mol load pdb 1YMG-all.pdb
#set model [atomselect top "protein and water and index 411 415 455 426 407 436"]
#$model writepdb abc.pdb
#exit
5楼2011-02-16 12:35:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


已解决,单独生成水的pdb文件和蛋白的pdb文件(因为没有金属离子等),然后合并这两个pdb文件即可
6楼2011-02-16 15:01:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nufang19a 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见