求助计算分子的亲脂性logP的值是实验数据的10倍

weixue · 发表于 Post on 2022-5-12 10:07:36

大家好！想计算一下实验分子的亲脂性logp，采用了https://doi.org/10.26434/chemrxiv.11320130 这篇文章的计算方法，如下图所示。但是计算出来的值是实验数据的10倍，请问这合理吗？用这种计算方法可靠吗？

管理员注：本帖作者严重不当引用北京科音基础量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）的ppt，已删除处理。ppt见12L

chands · 发表于 Post on 2022-5-12 10:24:59

是logP差了10倍，还是P差了10倍？

weixue · 发表于 Post on 2022-5-12 10:26:12

chands 发表于 2022-5-12 10:24
是logP差了10倍，还是P差了10倍？

您好，算出来logP比实验值大十倍

chands · 发表于 Post on 2022-5-12 10:30:43

weixue 发表于 2022-5-12 10:26
您好，算出来logP比实验值大十倍

你做优化时是各自带着溶剂模型吗，还有这个分子在溶液中不要有太复杂得过程，比如解离。

weixue · 发表于 Post on 2022-5-12 10:32:13

chands 发表于 2022-5-12 10:30
你做优化时是各自带着溶剂模型吗，还有这个分子在溶液中不要有太复杂得过程，比如解离。

是的，我的分子是带正电荷的金属离子，直接用的m052x/lanl2dz/6-31g(d) 加上smd溶剂模型计算的

chands · 发表于 Post on 2022-5-12 10:45:38

weixue 发表于 2022-5-12 10:32
是的，我的分子是带正电荷的金属离子，直接用的m052x/lanl2dz/6-31g(d) 加上smd溶剂模型计算的

你干脆把你的体系贴出来，说一下过程步骤，优化用的级别，算溶解自由能用的级别，让大家看看。

weixue · 发表于 Post on 2022-5-12 11:07:24

大家好！想计算一下实验分子的亲脂性logp，采用了https://doi.org/10.26434/chemrxiv.11320130 这篇文章的计算方法，如下图所示。我的分子是带两个正电荷的金属配合物，但是计算出来的值是实验数据的10倍，实验值是0.5，我计算的是5.4，请问这合理吗？用这种计算方法可靠吗？具体的输入文件如下图和附件所示。

weixue · 发表于 Post on 2022-5-12 11:08:33

chands 发表于 2022-5-12 10:45
你干脆把你的体系贴出来，说一下过程步骤，优化用的级别，算溶解自由能用的级别，让大家看看。

您好！我又发了一个新的帖子，谢谢您的帮助

weixue · 发表于 Post on 2022-5-12 11:26:53

这是计算的在两种溶剂中的结果文件

chands · 发表于 Post on 2022-5-12 14:12:15

weixue 发表于 2022-5-12 11:26
这是计算的在两种溶剂中的结果文件

你的文件一个是G16 A.03(Linux)，一个是G16 C.01(Linux), 我的GV是6.0.16Win，不知道为什么两个文件都打不开，奇怪的是我能打开我的G16 B.01(linux)，所以更详细的东西我看不到。这是过渡金属体系，你用这个级别优化结构未必最优，可以试用其他泛函加溶剂模型优化结构，再用M05-2X算溶解自由能。
（1）我不知道你的logP在哪个量级，咋眼一看这个分子在水和辛醇中的溶剂度都不会太小，相差10倍具体是多少，溶剂模型是有一定误差的，这个公式两个溶解自由能相减，不知误差是抵消了还是叠加更大了。
（2）第二个公式实际上假设分子在水中和辛醇中的结构是相同的（或者在水中的结构和在辛醇中的结构有相同的单点能），这个分子有一定柔性，在水中和辛醇中的结构是未必相同的，搞不好还要做构象平均啥的。
我暂时想到这两点，看看其他人还有什么想法。

wzkchem5 · 发表于 Post on 2022-5-12 14:53:47

因为实验上这个阳离子在辛醇/水之间分配的时候，必然带着阴离子一起，不可能只有阳离子从一相跑到另一相而阴离子不动。所以实验测得的不是这个阳离子自己的logP，而是阳离子加上抗衡阴离子这个总体系的logP。所以计算必须考虑抗衡阴离子，其中如果已知抗衡阴离子基本是完全电离开的，那么应该单独算抗衡阴离子再把自由能加进来；如果已知抗衡阴离子基本是完全没有电离开的，那么应该把抗衡阴离子放在阳离子附近一块优化、一块算（这个要做得严格的话还涉及到构型搜索的问题）

sobereva · 发表于 Post on 2022-5-12 20:39:40

完全相同的问题，发两个主题帖子干嘛？这次扣50 eV予以警告处理。这次我把帖子合并了，下次再发现积分扣光

引用我的培训班ppt的时候不得把顶端的版权信息截掉

已经把你帖子里的图片删除处理，下次再发现此行为将删帖处理

sobereva · 发表于 Post on 2022-5-12 20:43:32

weixue 发表于 2022-5-12 11:08
您好！我又发了一个新的帖子，谢谢您的帮助

这种情况不得重新发帖！
别人提示后，有要补充的内容直接在回帖里发，如果没别人回复过直接编辑原帖！仔细看置顶的新社员必读贴
相同问题发两个帖子会严重扰乱论坛秩序！

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“采用公式计算分子的亲脂性logP”改了，以后务必注意，下次再违规扣分处理

weixue · 发表于 Post on 2022-5-14 09:00:45

wzkchem5 发表于 2022-5-12 14:53
因为实验上这个阳离子在辛醇/水之间分配的时候，必然带着阴离子一起，不可能只有阳离子从一相跑到另一相而 ...

您好！假如我不用公式计算，有什么理论计算软件能够比较好的计算logP吗，主要是想和实验进行对应

weixue · 发表于 Post on 2022-5-14 09:00:57

chands 发表于 2022-5-12 14:12
你的文件一个是G16 A.03(Linux)，一个是G16 C.01(Linux), 我的GV是6.0.16Win，不知道为什么两个文件都打 ...

您好！假如我不用公式计算，有什么理论计算软件能够比较好的计算logP吗，主要是想和实验进行对应

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 求助计算分子的亲脂性logP的值是实验数据的10倍

采用公式计算分子的亲脂性logP ，计算值比实验值大10倍

评分 Rate

浏览过的版块