24小时热门版块排行榜    

查看: 1828  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

nufang19a

金虫 (正式写手)


[交流] 【讨论】水分子的两个HT原子之间不应该成键 已有5人参与

我知道产生问题的原因是和坛子里的一个同学一样:水分子的两个HT原子之间不应该成键。在产生psf是出现的问题,可我不会解决,望大家帮忙
在运行后出现这个问题:UNABLE TO FIND ANGLE PARAMETERS FOR HT HT OT(如图)

这是我的命令:
去掉硫酸根:
mol load pdb 1KIM.pdb
set sel [atomselect top "not resname SO4"]
$sel writepdb 1kims.pdb
$sel writepsf 1kims.psf
exit
去掉THM:
mol load pdb 1kims.pdb
set sel [atomselect top "not resname THM"]
$sel writepdb 1kimst.pdb
exit
分链:
mol load pdb 1kimst.pdb
set a {atomselect top "protein and chain A"]
set b {atomselect top "protein and chain B"]
$a writepdb 1kimstA.pdb
$b writepdb 1kimstB.pdb
exit
mol load pdb 1kimst.pdb
set a {atomselect top "water and chain A"]
set b {atomselect top "water and chain B"]
$a writepdb crystwatA.pdb
$b writepdb crystwatB.pdb
exit

产生psf:
package require psfgen
topology top_all27_prot_lipid.inp
pdbalias residue HIS HSE
pdbalias atom ILE CD1 CD
pdbalias atom HOH O OH2
pdbalias residue HOH TIP3
foreach S { A B } {
segment $S {
pdb 1kimst$S.pdb
}
coordpdb 1kimst$S.pdb $S
regenerate angles dihedrals
segment WC$S {
auto none
pdb crystwat$S.pdb
}
coordpdb crystwat$S.pdb WC$S
}
guesscoord
writepdb 1kim.pdb
writepsf 1kim.psf

加水:
set all [atomselect top all]
measure minmax $all
package require solvate
solvate 1kim.psf 1kim.pdb -o 1kim_wb -b 1.5 -minmax {{-10 48 8} {70 120 80}}
平衡动力学:
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


自己顶上去
2楼2010-12-17 15:10:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


我在立场和psf文件都发现了HT OT啊
4楼2010-12-17 15:40:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


好的,谢谢你,但我不会加。能否指点一下
6楼2010-12-17 17:15:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nufang19a

金虫 (正式写手)


★ ★
zh1987hs(金币+2):鼓励自己解决问题!呵呵 2010-12-17 22:19:06
解决了,谢谢coolrainbow 的提醒,十分感谢。
把top..27.inp中的H1和H2删掉就可以了
7楼2010-12-17 17:56:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nufang19a 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见