计算化学公社

标题: 关于使用gromacs进行批量模拟的可行性与步骤问题 [打印本页]

作者
Author: Soulmate 时间: 3 hour ago
标题: 关于使用gromacs进行批量模拟的可行性与步骤问题

各位大佬好，鄙人现在手上有100对蛋白质-配体的复合物需要模拟，请问是否可以通过编写脚本的方式进行复合物的批量模拟？


      对于小分子配体，鄙人决定还是按照ORCA/Multiwfn联用计算RESP电荷，再使用Sobtop生成小分子的拓扑文件，然后加入RESP电荷。
      最后，与大分子蛋白组合生成新的复合物.gro文件和.top文件，接着进行能量极小化、预平衡和产生相。

      如果可行，请问怎样操作才能最大程度地提高工作效率，以及大致的步骤和注意事项是什么？
      谢谢大家！

作者
Author: student0618 时间: 2 hour ago
本帖最后由 student0618 于 2026-4-28 18:44 编辑

基本的注意事项是

Bash脚本把命令写进去，用各种sed awk cat head tail grep 等常用指令修改文件
非必要但让Project更整齐：自动化脚本为方便很多人习惯每一个stage 用相同名字的输入/输出文件，也可避免文件名太长。分开目录以免数据被覆盖，也方便后续处理。
确保输入文件Template正确、蛋白完整可用(或加入e.g. pdb2pqr/pdbfixer处理蛋白的步骤/脚本)
确保配体净电荷不是零的自动算对用来处理算Partial charge的量化输入文件
确保流程及mdp正确
确保不同软件的环境正确安装
理清是本地还是要在集群跑，集群的话最好有交作业的模板
根据目的，跑完MD要作什么后处理什么分析，分析用指令脚本预备好。

实作的话：

手写就是找个体系跑一遍把bash history直接复制贴上，文件名改variable。
理解自己在作什么又懂得找LLM agent的，知道自己要什么、文件放哪、权限设定好，自动就写好脚本、可能还自动跑完了。
我自己是先手写，后来直接让agent把工作流及脚本写成skill来用。

如果用agent的前提是要知道自己在做什么、可以受的风险等，不然是浪费时间及token。

作者
Author: Soulmate 时间: 1 hour ago

student0618 发表于 2026-4-28 18:18
基本的注意事项是

Bash脚本把命令写进去，用各种sed awk cat head tail grep 等常用指令修改文件

没怎么写过脚本，听着有点复杂，但是还是有所收获，谢谢大佬

作者
Author: jackshoulder 时间: halfhour ago
什么研究需要跑100对蛋白-小分子的MD，你组里有这么好的资源在短时间内跑100组？先不谈体系大小，规整的（近似球形）蛋白-小分子复合物，100ns至少也得一天吧，一台机子。多台同时跑另提。那怕是提交给服务器，这机时，成本也不容小觑呢。你如果是做的虚拟筛选，完全没必要跑100组的MD。100组MD，so crazy

欢迎光临计算化学公社 (http://bbs.keinsci.com/)