24小时热门版块排行榜    

查看: 607  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 0703化学调剂 ,六级已过,有科研经历 +13 曦熙兮 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
信息提示
请填处理意见