24小时热门版块排行榜    

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

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

[交流] 【求助】*重金悬赏*静电能何以如何离谱?

Energy                      Average       RMSD     Fluct.      Drift  Tot-Drift
-------------------------------------------------------------------------------

[蛋白-配体复合物在水中的模拟]
Coul-SR: Protein-LIG        -358.097    25.0043    24.3992 -0.0633488   -18.9414
LJ-SR: Protein-LIG          -145.162    15.8001    15.6884 -0.0217339   -6.49847
Coul-SR: LIG-SOL             -273.128   31.8265    31.5378 -0.0495503   -14.8157
LJ-SR: LIG-SOL              -36.9597      14.52    14.2306  0.0334187    9.99226


[配体单独在水中的模拟]
Coul-SR: LIG-SOL           -824.378    47.1395    47.1384 -0.00357615   -1.07214
LJ-SR: LIG-SOL              -94.3993    23.6678    22.2444 -0.0934096   -28.0044



各位朋友,请帮忙!我要算蛋白质-配体之间的结合自由能,上面绿色的是复合物在水中的模拟,下面蓝色的是单独配体在水中的模拟,最后的算式里需要用[复合物]-[配体],结果,你看我红色标注的部分,是静电作用项,(-358.097)+(-273.128) - (-824.378) 配体单独在水中的静电能怎么这么大呀!我跑的都是阳性抑制剂,最后的能量由于这一离谱的项变成正的了!问题就出在这个-824.378身上,我的配体确实有两个羧基,并且都是脱去H的质子化状态,请问该如何解决?

[ Last edited by lei0736 on 2009-11-13 at 21:49 ]
回复此楼
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★ ★ ★
lei0736(金币+1,VIP+0):谢谢 6-21 22:04
gwdavid(金币+5,VIP+0):谢谢,请再帮忙看一下 6-22 10:00
我不是很懂。这样相减的物理意义是什么?
3楼2009-06-21 20:50:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

recoli

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+3,VIP+0):谢谢 6-21 22:03
gwdavid(金币+5,VIP+0):谢谢,请再帮忙看一下 6-22 10:00
这样的能量能直接相减吗。。。我没做过这一类计算,只是有点疑问。。。

[复合物]-[配体],听起来有点像算自由能?但是算自由能也不是直接把势能相减的。

羧基负离子和水分子肯定有很强的氢键,所以我觉得模拟结果应该是合理的,是否需要重新考虑一下计算公式呢,当然我只是推测。。。
2楼2009-06-21 19:45:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

谢谢二位的解答,我本为了图方便,没有写计算公式,这个是LIE方法计算结合自由能的算法,就是将静电项相减乘以一个系数,再加上范德华项相减乘上一个系数。

可是我的配体分子在水中的静电作用力要远远高于结合在蛋白质上的,从这里看来这么强的水合作用力貌似不可能与蛋白质结合,但事实上这是一个成药辛伐他汀的。

怎么这么高呢。。。。。。。
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
4楼2009-06-22 08:48:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★
erylingjet(金币+2,VIP+0):感谢交流 6-22 17:38
gwdavid(金币+5,VIP+0):谢谢,请再帮忙看一下 6-23 09:25
我觉得LIE的自由能也不是这样算的。

比如要研究 [蛋白+配体] 这个整体,A状态是  [蛋白+配体] 在水溶液中,B状态是 [配体] 在水溶液中,能量差值应该是:A状态中 [蛋白+配体] 整体与水的作用  减去  B状态中 [配体] 与水的作用。我觉得你的A状态的能量取得不正确。

另外,因为蛋白比配体大很多,所以更合适的方法是将配体消去,所以B状态最好是 [蛋白] 在水溶液中。这样结果可能会更好。
5楼2009-06-22 15:06:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:40 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:29 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 04:01 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
请填处理意见