g_sdf

Langue: en

Autres versions - même langue

Version: 259238 (debian - 07/07/09)

Section: 1 (Commandes utilisateur)

NAME

g_sdf - calculate spatial distribution function

SYNOPSIS

g_sdf -s topol.tpr -f traj.xtc -n index.ndx -o gom_plt.dat -r refmol.gro -[no]h -nice int -b time -e time -dt time -mode int -triangle vector -dtri vector -bin real -grid vector

DESCRIPTION

g_sdf calculates the spatial distribution function (SDF) of a set of atoms within a coordinate system defined by three atoms. There is single body, two body and three body SDF implemented (select with option -mode). In the single body case the local coordinate system is defined by using a triple of atoms from one single molecule, for the two and three body case the configurations are dynamically searched complexes of two or three molecules (or residues) meeting certain distance consitions (see below).

OPTIONS

-s topol.tpr
Structure+mass(db) (Input, optional): tpr tpb tpa gro g96 pdb xml
-f traj.xtc
Generic trajectory (Input): xtc trr trj gro g96 pdb
-n index.ndx
Index file (Input, Output)
-o gom_plt.dat
Generic data file (Output)
-r refmol.gro
Generic structure (Output, optional): gro g96 pdb xml
-[no]h
Print help info and quit (default: no).
-nice int
Set the nicelevel (default: 19).
-b time
First frame (ps) to read from trajectory (default: 0).
-etime
Last frame (ps) to read from trajectory (default: end of file).
-mode int
SDF in [1,2,3] particle mode (default: 1).
-triangle vector
r(1,3), r(2,3), r(1,2) (default: 0 0 0).
-dtri vector
dr(1,3), dr(2,3), dr(1,2) (default: 0.03 0.03 0.03).
-binreal
Binwidth for the 3D-grid (nm) (default: 0.05).
-grid vector
Size of the 3D-grid (nm,nm,nm) (default: 1 1 1).

USAGE

The program needs a trajectory, a GROMACS run input file and an index file to work. You have to setup 4 groups in the index file before using g_sdf:

The first three groups are used to define the SDF coordinate system. The programm will dynamically generate the atom tripels according to the selected -mode: In -mode 1 the triples will be just the 1st, 2nd, 3rd, ... atoms from groups 1, 2 and 3. Hence the nth entries in groups 1, 2 and 3 must be from the same residue. In -mode 2 the triples will be 1st, 2nd, 3rd, ... atoms from groups 1 and 2 (with the nth entries in groups 1 and 2 having the same res-id). For each pair from groups 1 and 2 group 3 is searched for an atom meeting the distance conditions set with -triangle and -dtri relative to atoms 1 and 2. In -mode 3 for each atom in group 1 group 2 is searched for an atom meeting the distance condition and if a pair is found group 3 is searched for an atom meeting the further conditions. The triple will only be used if all three atoms have different res-id's.

The local coordinate system is always defined using the following scheme: Atom 1 will be used as the point of origin for the SDF. Atom 1 and 2 will define the principle axis (Z) of the coordinate system. The other two axis will be defined inplane (Y) and normal (X) to the plane through Atoms 1, 2 and 3. The fourth group contains the atoms for which the SDF will be evaluated.

For -mode 2 and 3 you have to define the distance conditions for the 2 resp. 3 molecule complexes to be searched for using -triangle and -dtri.

The SDF will be sampled in cartesian coordinates. Use '-grid x y z' to define the size of the SDF grid around the reference molecule. The Volume of the SDF grid will be V=x*y*z (nm^3). Use -bin to set the binwidth for grid.

The output will be a binary 3D-grid file (gom_plt.dat) in the .plt format that can be be read directly by gOpenMol. The option -r will generate a .gro file with the reference molecule(s) transfered to the SDF coordinate system. Load this file into gOpenMol and display the SDF as a contour plot (see http://www.csc.fi/gopenmol/index.phtml for further documentation).

For further information about SDF's have a look at: A. Vishnyakov, JPC A, 105, 2001, 1702 and the references cited within.

ALTERNATIVES

g_spatial also produces spatial distribution functions. g_spatial differs from g_sdf in that g_spatial gives the user full control of the alignment procedure. The downside is that g_spatial requires more time and more disk space.

SEE ALSO

gromacs(7), g_spatial(1)

More information about the GROMACS suite is available in /usr/share/doc/gromacs or at http://www.gromacs.org/.

Further information about SDFs is published in A. Vishnyakov, JPC A, 105, 2001, 1702 and the references cited within.