摘要
为了进行有效的分子动力学模拟,研究人员需要精确的计算工具和算法。
GROMACS 是其中一个广泛使用的分子动力学模拟软件包,它特别优化了性能以处理大规模的生物分子系统,如蛋白质和脂质体。本文将介绍 GROMACS 中的 gmx solvate 工具,这是一个用于在溶剂中溶解溶质分子或生成溶剂盒子的实用程序,展示了如何使用它在进行复杂模拟前准备分子系统。
gmx 溶剂化物 核心指令
gmx solvate [-cp [<.gro/.g96/…>]] [-cs [<.gro/.g96/…>]]
[-p [<.top>]] [-o [<.gro/.g96/...>]] [-box <vector>]
[-radius <real>] [-scale <real>] [-shell <real>]
[-maxsol <int>] [-[no]vel]
描述
gmx solvate可以做以下两件事之一:
- 生成溶剂的盒子。指定
-cs
和-box
,或者指定-cs
和带有盒子但没有原子的结构文件的-cp
。 - 在溶剂分子的浴中溶解一个溶质配置,例如蛋白质。指定 -cp(溶质)和 -cs(溶剂)。除非设置了
-box
,否则使用溶质坐标文件(-cp)
中指定的盒子。如果你希望溶质在盒子中居中,程序gmx editconf
有复杂的选项可以改变盒子尺寸并将溶质居中。从盒子中移除任何溶质分子的原子与任何溶剂分子的原子之间距离小于两者原子的缩放范德华半径之和的溶剂分子。程序会读取一个范德华半径数据库(vdwradii.dat
),并按-scale
缩放得到的半径。如果数据库中找不到半径,则这些原子会被赋予(预缩放的)距离-radius
。注意,这些半径的有效性依赖于原子名称,因此随力场的不同而有很大差异。
默认的溶剂是简单点电荷水(SPC
),坐标来自 $GMXLIB/spc216.gro
。这些坐标也可用于其他三位点水模型,因为短暂的平衡会消除模型之间的细微差异。也支持其他溶剂以及混合溶剂。溶剂类型的唯一限制是一个溶剂分子只能包含一个残基。坐标文件中的残基信息被使用,因此应该或多或少保持一致。实际上,这意味着在溶剂坐标文件中两个连续的溶剂分子应该有不同的残基编号。通过从坐标文件中读取的坐标堆叠来构建溶质的盒子。这意味着这些坐标应在周期边界条件下进行平衡,以确保在堆叠接口上分子的良好对齐。-maxsol 选项
只添加最初的 -maxsol
溶剂分子,剩下的可以放入盒子中的将被留下。这可能会创建一个可能在后面造成问题的空洞。明智地选择你的体积。
如果设置 -shell
大于零,将在溶质周围放置一层指定厚度(纳米)的水层。提示:首先使用 gmx editconf
将蛋白质放在盒子中心是一个好主意。
最后,gmx solvate
将可选地从你的拓扑文件中删除已经添加了若干溶剂分子的行,并添加一行,其中包含你的坐标文件中的溶剂分子总数。
选项
-cp[<.gro/.g96/…>] (protein.gro) (可选)
结构文件:gro g96 pdb brk ent esp tpr
-cs[<.gro/.g96/…>] (spc216.gro) (图书馆)
结构文件:gro g96 pdb brk ent esp tpr
指定输入/输出文件的选项:
-p[<.top>] (topol.top) (可选)
拓扑文件
指定输出文件的选项:
-o[<.gro/.g96/…>] (出.gro)
结构文件:gro g96 pdb brk ent esp
其他选择:
-box<向量> (0 0 0)
盒子尺寸(单位:纳米)
-radius<真实> (0.105)
默认范德华距离
-scale<实际> (0.57)
比例因子,用于将 share/gromacs/top/vdwradii.dat 数据库中的范德华半径相乘。默认值 0.57 表示水中蛋白质的密度接近 1000 g/l。
-shell<真实> (0)
溶质周围可选水层的厚度
-maxsol<整数> (0)
如果溶剂分子适合盒子,则最多添加溶剂分子数。如果为零(默认),则忽略此值
-[no]vel (不)
保持输入溶质和溶剂的速度
注意
分子在初始配置中必须是完整的哦。
返回主要文章:《【分子动力学】 分子动力学新手入门:一文读懂GROMACS使用全流程,轻松开启模拟之旅! GROMACS初学者了解资料 GROMACS安装之前的准备 windows如何安装Gromacs》