应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】*重金悬赏*静电能何以如何离谱？
51/1返回列表
查看: 856  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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]
1楼 2009-06-21 15:59:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lccbat

金虫 (初入文坛)
★ ★ ★ ★ ★ ★
gwdavid(金币+3,VIP+0):谢谢！ 7-14 10:41
lei0736(金币+3,VIP+0):谢谢 7-14 11:37
这个方法有意识地忽略了protein-solvent相互作用，可能是认为变化不大。我对这一点是有疑问的，这可能造成很大的误差。你的结果看起来合理，ligand在复合物中是-200多，单独是-800多，我说合理是指相对而言。
另，实际能结合的抑制剂，计算得到正值，也经常发生吧。呵呵
赞一下(13人) 回复此楼
8楼2009-07-14 01:09:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

recoli

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

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

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

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★ ★
lei0736(金币+1,VIP+0):谢谢 6-21 22:04
gwdavid(金币+5,VIP+0):谢谢，请再帮忙看一下 6-22 10:00
我不是很懂。这样相减的物理意义是什么？
赞一下(15人) 回复此楼
3楼2009-06-21 20:50:24
已阅   回复此楼   关注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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售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 6/300 2026-11-11 04:32 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 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.55.1.O.54,科目全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:59 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一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:32 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭