gmx editconf命令

文章目录

命令行功能
参数解释
- -f （指定输入的结构文件）
- -o(指定输出结果文件)
- -bt (盒子类型)
- -d （溶液与盒子的距离）
- -c（放在盒子中心）
- -box
- -center（指定盒子中心的位置）
- -angles

命令行功能

该命令的主要功能是将pdb或者gro文件中获取得到的分子放在一个盒子中（该盒子作为模拟的单位）

参数解释

-f （指定输入的结构文件）

-o(指定输出结果文件)

-bt (盒子类型)

相当于box type的简写，用来控制要添加的盒子的类型，可选参数：

triclinic:三斜方体
cubic：正方体
dodeahedron:十二面体
octahedron:八面体

-d （溶液与盒子的距离）

由于盒子是等边长的立方体，因此，我们只用-d参数设定一个距离即可。单位为纳米。

-c（放在盒子中心）

表明将会把分子放在盒子的中心位置。假如没有指定改参数，那么就是单纯的把分子放在盒子中。

-box

指定box的边长，直接-box a b c，a，b,c表示三个数字。指定该参数之后，再生成的文件里的最后一行能看到指定的a,b,c。这里注意，当a,b,c不同时，系统会直接取最小值，作为盒子的边长（因为已经用-bt指定了盒子是一个立方体，边长必然相同）

-center（指定盒子中心的位置）

默认情况下盒子的中心为0,0,0,此时我们可以设置其他的坐标来作为中心，这样分子也会处在我们指定的这个中心上。
比如
我用以下命令加的盒子

gmx editconf -f RIB.gro -o newbox_RIB.gro -d 1.0 -c -bt cubic

那么它生成的gro文件最后的位置信息如下

(Abs_Met_Oxydation) [fanxuezhe@fanxuezhe test_gro_pdb]$ tail newbox_RIB.gro 227ALA     CA 3373   4.453   2.372   1.658227ALA     HA 3374   4.354   2.387   1.665227ALA     CB 3375   4.499   2.399   1.514227ALA    HB1 3376   4.457   2.333   1.453227ALA    HB2 3377   4.473   2.492   1.487227ALA    HB3 3378   4.599   2.390   1.509227ALA      C 3379   4.501   2.233   1.697227ALA    OT1 3380   4.454   2.135   1.632227ALA    OT2 3381   4.584   2.222   1.7918.02702   8.02702   8.02702

此外我利用-center命令指定中心之后

gmx editconf -f RIB.gro -o newbox_RIB.gro -d 1.0 -bt cubic -center 2 2 2

文件最后的信息如下

(Abs_Met_Oxydation) [fanxuezhe@fanxuezhe test_gro_pdb]$ tail newbox_RIB.gro 227ALA     CA 3373   2.439   0.358  -0.355227ALA     HA 3374   2.340   0.373  -0.348227ALA     CB 3375   2.485   0.385  -0.499227ALA    HB1 3376   2.443   0.319  -0.560227ALA    HB2 3377   2.459   0.478  -0.526227ALA    HB3 3378   2.585   0.376  -0.504227ALA      C 3379   2.487   0.219  -0.316227ALA    OT1 3380   2.440   0.121  -0.381227ALA    OT2 3381   2.570   0.208  -0.2228.02702   8.02702   8.02702

这里可以看出位置信息发生了巨大的变化，但是盒子的大小并未发生改变。

-angles

坐标空间中的坐标之间的角度，默认为90,90,90。

gmx editconf命令相关推荐

Gromacs的文件、程序、命令介绍
Gromacs的文件.程序.命令介绍 2017收集,来源于网络,整理如下: 1.文件格式 (1).pdb PDB数据库使用的文件格式: (2).top 拓扑文件,包含所有力场参数: (3).gro 分 ...
windows下gromacs中文教程（simulate chain A of insulin (PDB ID: 1ZNI).
1. 蛋白质准备文件从蛋白质文件中删除其他链和链,包括水分子(HOH).在Ubuntu和Windows中,您可以使用vim.notepad++等编辑器.也可以使用以下命令删除hetatom: #li ...
基于Gromacs的蛋白分子动力学模拟（RMSD、RMSF及蛋白的回旋半径）
一.实验要求实验对象:目标体系为modeller或其他方法建模的结果中评价最好的模型. 软件: Gromacs-5.1.2 二.实验步骤加立场 gmx pdb2gmx –h 打开帮助菜单. 选力场 ...
CADD课程学习（13）-- 研究蛋白小分子动态相互作用-I(GROMACS)
CADD课程学习(13)-- 研究蛋白小分子动态相互作用-I(GROMACS) 分子动力学基本原理分子动力学(Molecular Dynamics-MD)一门结合物理,数学和化学的综合技术. 分子动 ...
GROMACS Tutorial 1-Lysozyme in Water-translated with notes
Lysozyme Tutorial 一些GROMACS基础知识 pdb2gmx The Topology Solvation Adding Ions Energy Minimization Equil ...
Gromacs分子动力学模拟流程概述
Gromacs分子动力学模拟主要可以分为以下几个步骤,不同的体系步骤可能略有不同. 在开始之前,先简单了解一下预平衡: 分子动力学模拟的最终目的是对体系进行抽样,然后计算体系的能量,各种化学键,成分分 ...
生成小分子力场TOP
生成小分子力场TOP 自动化工具是非常受欢迎的.对于每一个力场,都有一些方法或软件程序可以提供与各种力场兼容的参数.并不是所有的结果都一样准确.有几个例子,请参考下表: 力场自动化工具备注 AMB ...
分子模拟软件amber_[gromacs使用教程] 基于amber力场模拟蛋白小分子复合物
祥请参考官网教程,使用其中的mdp参数文件(均100ps),案例只考虑模拟顺利,暂不考虑合理性. 平台:windows 软件:gaussina16, ambertools, gromacs2019.6 ...
基于Gromacs的蛋白水溶液分子动力学模拟
1. 检查结构文件有些结构文件存在少几个氢原子或者侧链的情况,所以先用spdbv软件打开结构文件,该软件可自动补加缺失的分子,用这个软件打开结构文件,再另存一下结构即可. Spdbv软件是windo ...
GROMACS 教程--水中的溶菌酶
GROMACS 教程–水中的溶菌酶中文教程:http://jerkwin.github.io/ 英文教程:http://www.mdtutorials.com 此示例将指导新用户完成模拟系统的设置过 ...