| 查看: 1640 | 回复: 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或是其它? |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有210人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
gwdavid
木虫 (著名写手)
小木虫散兵坑坑长
- 应助: 0 (幼儿园)
- 贵宾: 0.522
- 金币: 4097.7
- 散金: 1548
- 红花: 57
- 帖子: 1768
- 在线: 151.2小时
- 虫号: 514002
- 注册: 2008-02-28
- 性别: GG
- 专业: 生化分析及生物传感
★ ★ ★ ★ ★ ★
csfn(金币+6,VIP+0):3q
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项，首先是字面翻译，其次是我认为分子间的结合自由能计算不应包括分子内能，愿与楼主讨论此判断结果！ |
2楼2008-07-01 20:49:59
chemozhang
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3431.1
- 散金: 462
- 红花: 12
- 帖子: 615
- 在线: 356小时
- 虫号: 282954
- 注册: 2006-10-08
- 性别: GG
- 专业: 药物学
★ ★ ★
csfn(金币+3,VIP+0):3q
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
weishenme
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 981.3
- 帖子: 128
- 在线: 52.6小时
- 虫号: 45892
- 注册: 2004-05-15
- 性别: GG
- 专业: 生物大分子结构与功能
4楼2008-07-02 04:37:21
luai071
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 3000.5
- 散金: 100
- 红花: 1
- 帖子: 1348
- 在线: 126.8小时
- 虫号: 351005
- 注册: 2007-04-21
- 专业: 合成药物化学
5楼2008-07-02 07:15:13
chemozhang
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 3431.1
- 散金: 462
- 红花: 12
- 帖子: 615
- 在线: 356小时
- 虫号: 282954
- 注册: 2006-10-08
- 性别: GG
- 专业: 药物学
★
csfn(金币+1,VIP+0):3q
csfn(金币+1,VIP+0):3q
|
感谢楼上的两位 但我还是觉得2楼说的很有道理的,AutoDocK中评价函数中不应该包含内能(差值),设想若对接的起始构象不同,理论上如果对接好的化,它们会得到同一构象和取向,那么这里internal energy(差值)是不同的,显然不能包含在AutoDocK中评价函数中 我查阅了徐莜杰的书及相关资料以及tutarial中的解释,得出一下结论 △G(bind)(评价函数)= △G(vdw)+△G(H-bond)+△G(ele)+△G(del-sol)+ △G(tor) 经过计算证实前四项包含在intermolecular energy中, 另外,free energy=intermolecular energy+torsion entropy penalty 正如2楼兄弟所说 freeEnergy: 自由能＝分子间能量项+扭转熵因子罚分项（这一项的值是一个常数乘以配体内可旋转键数目，这一项的值是那个unbound extended energy,可能性比较大，因为这两项加起来很接近freeEnergy的值，但又不是绝对相等，小数点后总差点。） 对接结果中 Torsional Energy : 3.84 Unbound Extended Energy: -0.83 这两个值是固定的,即不论你排名如何,构象如何,他们的值都是这个.这两个能量表示的是同一个东西. 所以推测Unbound Extended Energy=△G(tor) 这样 △G(bind)(评价函数)= △G(vdw)+△G(H-bond)+△G(ele)+△G(del-sol)+ △G(tor) =intemolecar energy+ △G(tor)=free energy约=binding energy 而docked energy在这里对我们来说没有多少价值 |
6楼2008-07-02 10:42:49
gwdavid
木虫 (著名写手)
小木虫散兵坑坑长
- 应助: 0 (幼儿园)
- 贵宾: 0.522
- 金币: 4097.7
- 散金: 1548
- 红花: 57
- 帖子: 1768
- 在线: 151.2小时
- 虫号: 514002
- 注册: 2008-02-28
- 性别: GG
- 专业: 生化分析及生物传感
7楼2008-07-08 17:16:41
