24小时热门版块排行榜    

查看: 735  |  回复: 0

codyliu

铁杆木虫 (正式写手)

[求助] 想计算配体和受体之间的氢键 能运行

群里的高手们,我写了一个vmd的tcl,想计算配体和受体之间的氢键
能运行,但是一个氢键都没有计算出来,请大家帮我看看怎么回事呗
set outfile [open hbond.dat w]
set nf [molinfo top get numframes]
set selA [atomselect top "index 577"]
# frames loop
for { set i 0 } { $i <= $nf } { incr i } {
for {set t 1} {$t < 577 } {incr t } {
set selB [atomselect top "index $t"]
$selA frame $i
$selB frame $i
puts $outfile "[measure hbonds 3.5 30 $selB  $selA ]"
}
}
close $outfile
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 codyliu 的主题更新
信息提示
请填处理意见