| 查看: 582 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
wuhanhgf2002金虫 (正式写手)
|
[交流]
【求助】关于氢键
|
||
|
请问有没有查询mol2结构中氢键供体和受体数目的软件？最好是linux下的。谢谢 [ Last edited by zdhlover on 2009-11-25 at 13:09 ] |
回复此楼» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有14人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
5楼2009-03-11 09:36:38
wuhanhgf2002
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 323.8
- 帖子: 435
- 在线: 43.1小时
- 虫号: 121302
- 注册: 2005-12-03
- 性别: GG
- 专业: 有机合成
2楼2009-03-10 09:23:40
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★ ★
xuefei06(金币+6,VIP+0):thanks！ 3-10 23:41
xuefei06(金币+6,VIP+0):thanks！ 3-10 23:41
|
vmd 提供了计算HB的接口，但具体用法需要自己写脚本。 下面是一个计算水中每个分子同周围形成氢键的数目。最关键的是 measure hbonds $dist $ang 这一行，但全文我还没有看明白。 set all [atomselect top all] # determine the dimensions of the water box set mes [measure minmax $all] set xmin [lindex [lindex $mes 0] 0] set ymin [lindex [lindex $mes 0] 1] set zmin [lindex [lindex $mes 0] 2] set xmax [lindex [lindex $mes 1] 0] set ymax [lindex [lindex $mes 1] 1] set zmax [lindex [lindex $mes 1] 2] set O [atomselect top "name OH2"] # set the number of atoms for which you want to count H-bonds set atomcount 300.0 set totalbondnum 0 for {set i 1} {$i <= $atomcount} {incr i} { # pick a random water molecule (oxygen atom) and record its position set rand [expr int(rand()*[expr [$O num] - 1])] set oatom [atomselect top "index [lindex [$O get index] $rand]"] set xpos [lindex [lindex [$oatom get {x y z}] 0] 0] set ypos [lindex [lindex [$oatom get {x y z}] 0] 1] set zpos [lindex [lindex [$oatom get {x y z}] 0] 2] set buffer 3 # the selected atom should be away from the surface (by if {$xpos < [expr $xmin+$buffer] || $xpos > [expr $xmax-$buffer] || $ypos < [expr $ymin+$buffer] || $ypos > [expr $ymax-$buffer] || $zpos < [expr $zmin+$buffer] || $zpos > [expr $zmax-$buffer]} then { set i [expr $i-1] } else { # set the distance of a H-bond to 3.0 Angstroms and the angle to 30 degrees (default values) set dist 3.0 set ang 30.0 set temp1 [atomselect top "same residue as index [$oatom get index]"] set temp2 [atomselect top "same residue as exwithin 10 of index [$temp1 get index]"] # count the number of H-bonds for th molecule set count1 [llength [lindex [measure hbonds $dist $ang $temp1 $temp2] 0]] set count2 [llength [lindex [measure hbonds $dist $ang $temp2 $temp1] 0]] # and add it to the total so far set totalbondnum [expr $totalbondnum+$count1 + $count2] } } puts "Each water molecule makes [expr double($totalbondnum/$atomcount)] hydrogen bonds on average with its neighbors." |
4楼2009-03-10 17:58:13