24小时热门版块排行榜    

查看: 608  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268求调剂 +8 简单点0 2026-03-17 8/400 2026-03-21 13:43 by 南粤飞花
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +6 墨墨漠 2026-03-20 6/300 2026-03-21 09:23 by 天青等烟雨
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 297求调剂 +9 戏精丹丹丹 2026-03-17 9/450 2026-03-21 01:49 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 0703化学调剂 ,六级已过,有科研经历 +13 曦熙兮 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见