| 查看: 2770 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助关于solvate的一些疑惑
|
|||
|
请问怎么判断溶解蛋白的水球是否足够大，我一直使用vmd的wat_sphere.tcl脚本将蛋白溶解于水溶液中，vmd 中看到的蛋白质肽链的末段好像有一部分在水球的外面。。。。。。。 有没有什么方法可以稍微增大水球的大小呀？求指点 |
回复此楼» 收录本帖的淘帖专辑推荐
 我学习计算的一些帖子 |
» 本帖已获得的红花（最新10朵）
Mars_kang» 猜你喜欢
心脉受损
已经有8人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求助求助~~~关于套磁滴一些疑问~！
已经有21人回复- 关于洪堡申请中的pulications list 中的一些疑惑
已经有8人回复
- 关于海外博后的一些疑惑
已经有8人回复
- 关于免疫组化SABC法的一些疑惑？
已经有3人回复
- 关于可用于保健食品的物品是不允许用于普通食品的一些疑惑
已经有11人回复
- 求助关于硅橡胶课题现在哪一个方向好一些！
已经有8人回复
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- 关于信号肽的一些疑惑
已经有10人回复
- 【求助】关于GPA的一些问题
已经有3人回复
- 【求助】关于计算NICS（核独立化学位移）的一些疑问A
已经有3人回复
- 【求助】关于自旋多重度（基态vs激发态）的一些疑问
已经有9人回复
- 【求助/交流】MTT实验，24孔板接种细胞的一些疑惑
已经有15人回复
- 【求助】关于PL的一点疑惑
已经有4人回复
3楼2012-07-01 21:57:25
2楼2012-07-01 16:47:22
4楼2012-07-01 23:45:27
|
### Script to immerse ubiquitin in a sphere of water just large enough ### to cover it set molname xia41h mol new ${molname}.psf mol addfile ${molname}.pdb ### Determine the center of mass of the molecule and store the coordinates set cen [measure center [atomselect top all] weight mass] set x1 [lindex $cen 0] set y1 [lindex $cen 1] set z1 [lindex $cen 2] set max 0 ### Determine the distance of the farthest atom from the center of mass foreach atom [[atomselect top all] get index] { set pos [lindex [[atomselect top "index $atom"] get {x y z}] 0] set x2 [lindex $pos 0] set y2 [lindex $pos 1] set z2 [lindex $pos 2] set dist [expr pow(($x2-$x1)*($x2-$x1) + ($y2-$y1)*($y2-$y1) + ($z2-$z1)*($z2-$z1),0.5)] if {$dist > $max} {set max $dist} } mol delete top ### Solvate the molecule in a water box with enough padding (20 A). ### One could alternatively align the molecule such that the vector ### from the center of mass to the farthest atom is aligned with an axis, ### and then use no padding package require solvate solvate ${molname}.psf ${molname}.pdb -t 20 -o del_water resetpsf package require psfgen mol new del_water.psf mol addfile del_water.pdb readpsf del_water.psf coordpdb del_water.pdb ### Determine which water molecules need to be deleted and use a for loop ### to delete them set wat [atomselect top "same residue as {water and ((x-$x1)*(x-$x1) + (y-$y1)*(y-$y1) + (z-$z1)*(z-$z1))<($max*$max)}"] set del [atomselect top "water and not same residue as {water and ((x-$x1)*(x-$x1) + (y-$y1)*(y-$y1) + (z-$z1)*(z-$z1))<($max*$max)}"] set seg [$del get segid] set res [$del get resid] set name [$del get name] for {set i 0} {$i < [llength $seg]} {incr i} { delatom [lindex $seg $i] [lindex $res $i] [lindex $name $i] } writepsf ${molname}_ws.psf writepdb ${molname}_ws.pdb mol delete top mol new ${molname}_ws.psf mol addfile ${molname}_ws.pdb puts "CENTER OF MASS OF SPHERE IS: [measure center [atomselect top all] weight mass]" puts "RADIUS OF SPHERE IS: $max" mol delete top 我对solvate ${molname}.psf ${molname}.pdb -t 20 -o del_water这一行命令不是很懂。。。。。还有版主一般使用什么工具来进行溶剂化呀，昨天试了下packmol，但只是得到了pdb文件，不知道psf文件怎么得到？ |
5楼2012-07-02 09:43:13