24小时热门版块排行榜    

查看: 1867  |  回复: 6
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chemozhang

木虫 (正式写手)

[交流] [交流]关于AutoDock输出结果中的能量值

下面是autodock运行 打分值最好的输出结果
Rank: 1_1
Binding Energy:             -15.1
kI :                        8.58pM
Intermolecular Energy :     -14.39
Internal Energy :           -5.38
Torsional Energy :          3.84
Unbound Extended Energy:    -0.83
Cluster RMS:                 0.0
Ref RMS:                     6.79

tutorial中说Dock Energy=intermolecular energy+internal energy
那么在此例中Dock Energy=-14.39+(-5.38)=-19.77,而Binding Energy=-15.1,两者并不相等,
但两者是平行的,按这两个能量值进行排序都得到相同的排序结果

我很困惑的是 如果Binding Energy就是我们关心的Binding Free Energy(结合自由能△G),那么Dock Energy又是什么含义?
Dock Energy and Binding Energy有什么不同,他们都是怎么定义的?
还有这里Binding Energy的单位是kcal/mol还是kJ/mol或是其它?
回复此楼
故地如重游月圆更寂寞
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemozhang

木虫 (正式写手)

★ ★ ★
csfn(金币+3,VIP+0):3q
分子间的结合自由能计算不应包括分子内能,这点我非常的认同.
我的理解是文献中引用的能量值△G(结合自由能,binding free energy)就是输出中的binding energy,即 △G=binding free energy=binding energy,AutoDocK结果中的△G应该直接计算得出,不可能再经过计算得到.
freeEnergy与Binding energy 虽然很接近,但不绝对相等.很可能是它们的定义略有不同,而不能将他们等同起来,因为我觉得AutoDocK那么广泛使用的程序,不可能在这点上没有做到精确.
故地如重游月圆更寂寞
3楼2008-07-01 22:22:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

★ ★ ★ ★ ★ ★
csfn(金币+6,VIP+0):3q
楼主,我愿意与您讨论这个问题,autodock里边的这些能量项确实有些令人迷惑,但是请您先看一下这个,我自己提取并按自己的理解翻译总结的:
The upper panel displays information about the current
conformation. This includes its overall rank, for example the
best result is always 1_1: lowest energy cluster_best individual
in cluster.

Docked Energy is the sum of the intermolecular
and internal energy components.  

Cluster RMS is the root
mean square difference rms between this individual and the
seed for the cluster. 1_1 is the seed for the first cluster so its
Cluster RMS is 0.0. Ref RMS is the rms between the specified
reference structure.  If no reference structure is specified in the
DPF, the input structure is used as the reference. freeEnergy
is the sum of the intermolecular energy plus the torsion entropy
penalty which is a constant times the number of rotatable bonds
in the ligand, kI calculated from the Docked Energy
.


译文:(此为自己翻译如有不确切,敬请谅解,仅供参考理解)
上部面板显示当前构象的信息。
包括此构象的综合评价,例如最好的结果总是(1_1:最低能量类群_此群中的最好构象)

关于docked Energy: 对接能量=分子间能量+分子内能量
关于freeEnergy:    自由能=分子间能量项+扭转熵因子罚分项(这一项的值是一个常数乘以配体内可旋转键数目,此处我怀疑这一项的值要么是那个unbound extended energy,是这一项值的可能性比较大,因为这两项加起来很接近freeEnergy的值,但又不是绝对相等,小数点后总差点。)
按我的理解总结是:结合自由能应为freeEnergy项,首先是字面翻译,其次是我认为分子间的结合自由能计算不应包括分子内能,愿与楼主讨论此判断结果!
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
2楼2008-07-01 20:49:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weishenme

金虫 (小有名气)

★ ★
csfn(金币+2,VIP+0):3q
Binding Energy=Intermolecular Energy +Internal Energy +Torsional Energy- Unbound Extended Energy,单位是kcal/mol。
internal energy可能是指配体结合前后内能的差值,不是内能的绝对值。

[ Last edited by weishenme on 2008-7-2 at 04:53 ]
4楼2008-07-02 04:37:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luai071

木虫 (著名写手)


csfn(金币+1,VIP+0):3q
引用回帖:
Originally posted by weishenme at 2008-7-2 04:37:
Binding Energy=Intermolecular Energy +Internal Energy +Torsional Energy- Unbound Extended Energy,单位是kcal/mol。
internal energy可能是指配体结合前后内能的差值,不是内能的绝对值。

[ Last edite ...

不错,是这回事,楼主手里有没有徐莜杰编的书,翻一下就是了。
5楼2008-07-02 07:15:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
信息提示
请填处理意见