计算化学公社

 找回密码 Forget password
 注册 Register
Views: 67|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 模拟蛋白和多聚体糖类形成复合物的过程,有许多LINCS WARNING报错问题

[复制链接 Copy URL]

23

帖子

0

威望

175

eV
积分
198

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好:
我需要模拟蛋白和多聚体糖类形成复合物的过程。多糖通过charmm-gui获得pdb文件,然后在B3LYP-D3(BJ)/def2-SVP级别下对结构进行优化,然后在B3LYP-D3(BJ)/def2-TZVP级别下计算单点能同时产生分子表面静电势数据,最后调用Multiwfn产生RESP电荷。接着使用sobtop建立gaff力场下的拓扑结构。蛋白力场选用了amber14sb,但是在nvt时就出现了许多LINCS WARNING:
relative constraint deviation after LINCS:
bonds that rotated more than 30 degrees:
于是把蛋白和配体分别在水中进行模拟,发现单独模拟配体时还是出现该问题,于是把配体的模拟过程进行了多种尝试:

尝试了:
缩小了模拟步长1 fs乃至0.5 fs不加约束,采用不同方法能量最小化,采用双精度版本进行能量最小化,gen-temp = 310改为50K,设置退火,都没有解决问题。

下面是配体相关文件,还请老师们帮忙解决困惑,谢谢各位老师帮助。

topol.top

543 Bytes, 下载次数 Times of downloads: 1

lig_RESP.itp

190.67 KB, 下载次数 Times of downloads: 3

lig_RESP.mol2

21.25 KB, 下载次数 Times of downloads: 1

em_steep.mdp

1.09 KB, 下载次数 Times of downloads: 0

nvt.mdp

2.41 KB, 下载次数 Times of downloads: 0

npt.mdp

2.64 KB, 下载次数 Times of downloads: 0

5万

帖子

99

威望

5万

eV
积分
112476

管理员

公社社长

2#
发表于 Post on 2024-10-28 12:23:13 | 只看该作者 Only view this author
短的配体跑成功了再跑长的
另外,建议每个重复单元单独产生RESP电荷而别直接对整体产生,http://sobereva.com/soft/Sobtop里"关于给聚合物用RESP原子电荷"部分专门提到了
体系中有的是NH3有的是NH2,注意确保质子化状态没错
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 16:15 , Processed in 0.161444 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list