硕博科研从0到1|Gromacs下载安装&计算
- 2026-09-20 13:40:10
硕博科研从0到1|Gromacs下载安装&计算一、下载&安装 Conda一键下载安装. 具体命令（linux）： 二、计算 1. 准备工作 （1）力场选择和力场文件下载 ① CHARMM36 力场：由MacKerell等开发的全原子生物大分子力场，旨在高精度模拟蛋白质、核酸、脂质和碳水化合物等生物大分子。 下载链接：http://mackerell.umaryland.edu/CHARMM_ff_params.html （2）溶剂（除水外）的力场文件准备 SwissParam 为溶剂生成拓扑文件（生成与 CHARMM36 力场兼容的 .itp 拓扑文件）。 LigParGe 是OPLS-AA力场参数生成工具。由Jorgensen团队开发。 CGenFF力场在CGenFF 的参数生成网站 https://www.paramchem.org/生成。 需要用到的几种文件格式： 2. 模拟计算 （1）建立盒子→填充溶剂（水、有机溶剂）→添加离子（中和电荷） gmx editconf -f GG.gro -o box.gro -c -d 1.0 -bt cubic （将溶质GG置于盒子中心，与边缘保持 1.0 nm 距离） gmx insert-molecules -f box.gro -ci 溶剂.gro -o box_溶剂.gro -nmol N （插入指定数量的其它溶剂分子） gmx solvate -cp box_溶剂.gro -cs spc216.gro -o solv.gro -p topol.top（用水填满盒子剩余空间，自动更新拓扑） gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr （插入离子前的准备） gmx genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -neutral （插入离子替换部分溶剂） （2）模拟：能量最小化 （升温至目标温度）→NPT 平衡（调整密度和压力）→生产模拟（采集数据） gmx grompp -fem.mdp -c solv_ions.gro -p topol.top -o em.tpr -maxwarn1（建立输入文件） gmx mdrun -v -deffnm em （运行能量最小化） gmx grompp -f nvt.mdp -c em.gro -p topol.top -o nvt.tpr -maxwarn 1 （准备 NVT输入文件） gmx mdrun -v -deffnm nvt -ntmpi 1 （运行NVT平衡，升温至目标温度） gmx grompp -f npt.mdp -c nvt.gro -t nvt.cpt -p topol.top -o npt.tpr -maxwarn1 （准备 NPT输入文件） gmx mdrun -v-deffnm npt -ntmpi1 （运行NPT平衡，调整密度和压力） gmx grompp -f md.mdp -c npt.gro -t npt.cpt -p topol.top -o md.tpr -maxwarn 1 （生产模拟） gmx mdrun -v -deffnm md -ntmpi 1 （运行生产模拟（时间根据需要设定）） 输出： （3）数据分析 gmx rms -s md.tpr -f md.xtc -o rmsd.xvg （构象变化分析（RMSD）RMSD 值越低且平稳 ） gmx gyrate -s md.tpr -f md.xtc -o gyrate.xvg （计算回转半径（Rg）：Rg衡量分子尺寸和紧凑程度。Rg 值越小，分子结构越紧凑） gmx hbond -s md.tpr -f md.xtc -num hbnum.xvg （氢键分析） gmx energy -f md.edr -o energy.xvg（能量分析） gmx rdf -s md.tpr -f md.xtc -o rdf.xvg（计算径向分布函数：分析溶剂在溶质周围的分布密度） gmx sasa -s md.tpr -f md.xtc -o sasa.xvg（计算溶剂可及表面积：分析溶质与溶剂的接触面积） vmd md.gro md.xtc（vmd可视化轨迹） 输出：.xvg （作图数据） 对比不同数据时取稳定段的平均值 ± 标准差。
conda install -c conda-forge gromacs=2025.2 （安装指定版本）
conda create -n gromacs_env （创建一个名为 gromacs_env 的新环境）
conda activate gromacs_env （激活该环境）
conda install -c conda-forge gromacs （在新环境中安装 GROMACS）
conda search -c conda-forge gromacs （安装支持 GPU 的版本）
gmx --version
② OPLS-AA 力场：由Jorgensen等开发的全原子力场，其参数优化目标是精确再现纯有机液体的热力学和结构性质。
下载链接：https://traken.chem.yale.edu/oplsaam.html
③ CGenFF力场：主要由 Kenno Vanommeslaeghe 和 Alexander D. MacKerell Jr. 及其团队开发专门为类药分子和广泛有机分子提供参数，使其能与CHARMM力场完美兼容。
下载链接：http://dogmans.umaryland.edu/~kenno/cgenff/download.html
④ AMBER力场 ：其中通用小分子力场 GAFF 由 Junmei Wang（王军梅）、David A. Case、Peter A. Kollman 等开发（2004 年）；而蛋白质/核酸等生物大分子力场（如 ff14SB、ff99）则由 AMBER 开发团队历代成员持续贡献。
下载链接：https://ambermd.org/GetAmber.php，也可以conda一键下载
整个模拟体系（溶质 + 所有溶剂）应使用 同一力场家族 的参数（如都使用 GAFF/AMBER 体系，或都使用 CHARMM/CGenFF 体系），混用不同家族的力场会导致原子类型冲突。
Acpype 默认调用 AmberTools 的antechamber，其默认力场是 GAFF。.mol2：转化gromacs的输入文件。.top (或 .topol.top)主拓扑文件，引用所有 .itp 文件，定义体系组成。
.itp：最核心的gromacs计算输入文件，直接放到topol.top里引用。.pdb：可作为坐标输入，但建议转成.gro格式再使用。坐标文件，定义原子的初始位置。.rtf、.prm、.psf、.crd：过程文件，可以保留作为备份或检查，但对 GROMACS 模拟不是必需的。.trr/.xtc轨迹文件，记录原子坐标随时间的变化
.edr能量文件，记录体系的能量、温度、压力等随时间的变化
.mdp模拟参数文件，定义模拟条件（温度、压力、步长等）
.tpr最终的二进制输入文件，由 gmx grompp 生成，是 gmx mdrun 的唯一输入，把.gro,
.top,.mdp完整合并、压缩、二进制化成一个.tpr文件
.xvg计算输出文件，分析提取可读的数值数据，通过gmx energy/rmsd等命令从.edr/.xtc中提取，用于绘图
nano topol.top （建立top文件）md.xtc（轨迹）、md.edr（能量）、md.log（日志）
本文来自网友投稿或网络内容，如有侵犯您的权益请联系我们删除，联系邮箱：wyl860211@qq.com 。
