Scroll to navigation

GMX-DSSP(1) GROMACS GMX-DSSP(1)

NAME

gmx-dssp - Calculate protein secondary structure via DSSP algorithm

SYNOPSIS

gmx dssp [-f [<.xtc/.trr/...>]] [-s [<.tpr/.gro/...>]] [-n [<.ndx>]]


         [-o [<.dat>]] [-num [<.xvg>]] [-b <time>] [-e <time>]


         [-dt <time>] [-tu <enum>] [-fgroup <selection>]


         [-xvg <enum>] [-[no]rmpbc] [-[no]pbc] [-sf <file>]


         [-selrpos <enum>] [-seltype <enum>] [-sel <selection>]


         [-hmode <enum>] [-hbond <enum>] [-[no]nb] [-cutoff <real>]


         [-[no]clear] [-[no]pihelix] [-ppstretch <enum>]


DESCRIPTION

gmx dssp allows using the DSSP algorithm (namely, by detecting specific patterns of hydrogen bonds between amino acid residues) to determine the secondary structure of a protein.

One-symbol secondary structure designations that are used in the output file:

H — alpha-helix;

B — residue in isolated beta-bridge;

E — extended strand that participates in beta-ladder;

G — 3_10-helix;

I — pi-helix;

P — kappa-helix (poly-proline II helix);

S — bend;

T — hydrogen-bonded turn;

= — break;

~ — loop (no special secondary structure designation).

-num allows you to get a plot of the number of secondary structures of each type as a function of time at the output.

-hmode selects between using hydrogen atoms directly from the structure ("gromacs" option) and using hydrogen pseudo-atoms based on C and O atom coordinates of previous residue ("dssp" option). You should always use the "dssp" option for structures with absent hydrogen atoms!

-hbond selects between different definitions of hydrogen bond. "energy" means the calculation of a hydrogen bond using the electrostatic interaction energy and "geometry" means the calculation of the hydrogen bond using geometric criterion for the existence of a hydrogen bond.

-nb allows using GROMACS neighbor-search method to find residue pairs that may have a hydrogen bond instead of simply iterating over the residues among themselves.

-cutoff is a real value that defines maximum distance from residue to its neighbor residue used in -nb. Minimum (and also recommended) value is 0.9.

-clear allows you to ignore the analysis of the secondary structure residues that are missing one or more critical atoms (CA, C, N, O or H). Always use this option together with -hmode dssp for structures that lack hydrogen atoms!

-pihelix changes pattern-search algorithm towards preference of pi-helices.

-ppstretch defines stretch value of polyproline-helices. "shortened" means stretch with size 2 and "default" means stretch with size 3.

Note that gmx dssp currently is not capable of reproducing the secondary structure of proteins whose structure is determined by methods other than X-ray crystallography (structures in .pdb format with incorrect values in the CRYST1 line) due to the incorrect cell size in such structures.

Please note that the computation is always done in single precision, regardless of the precision for which GROMACS was configured.

OPTIONS

Options to specify input files:

Input trajectory or single configuration: xtc trr cpt gro g96 pdb tng
Input structure: tpr gro g96 pdb brk ent
Extra index groups

Options to specify output files:

Filename for DSSP output
Output file name for secondary structures statistics for the trajectory

Other options:

First frame (ps) to read from trajectory
Last frame (ps) to read from trajectory
Only use frame if t MOD dt == first time (ps)
Unit for time values: fs, ps, ns, us, ms, s
Atoms stored in the trajectory file (if not set, assume first N atoms)
Plot formatting: xmgrace, xmgr, none
-[no]rmpbc (yes)
Make molecules whole for each frame
-[no]pbc (yes)
Use periodic boundary conditions for distance calculation
Provide selections from files
Selection reference positions: atom, res_com, res_cog, mol_com, mol_cog, whole_res_com, whole_res_cog, whole_mol_com, whole_mol_cog, part_res_com, part_res_cog, part_mol_com, part_mol_cog, dyn_res_com, dyn_res_cog, dyn_mol_com, dyn_mol_cog
Default selection output positions: atom, res_com, res_cog, mol_com, mol_cog, whole_res_com, whole_res_cog, whole_mol_com, whole_mol_cog, part_res_com, part_res_cog, part_mol_com, part_mol_cog, dyn_res_com, dyn_res_cog, dyn_mol_com, dyn_mol_cog
Group for DSSP
Hydrogens pseudoatoms creating mode: gromacs, dssp
Selects between different definitions of hydrogen bond: energy, geometry
-[no]nb (yes)
Use GROMACS neighbor-search method
Distance from residue to its neighbor residue in neighbor search. Must be >= 0.9
-[no]clear (no)
Clear defective residues from the structure
-[no]pihelix (no)
Prefer Pi Helices
Stretch value for PP-helices: shortened, default

SEE ALSO

gmx(1)

More information about GROMACS is available at <http://www.gromacs.org/>.

COPYRIGHT

2024, GROMACS development team

January 30, 2024 2024