GMX-DSSP(1)

πŸ’– GMX-DSSP πŸ’–

okay bestie, welcome to the GROMACS manual page for gmx-dssp β€” yes it's a real page in the GROMACS manual, section 1, and yes we're gonna get through it together 🫢


✨ 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>]
         [-[no]polypro]

πŸ¦‹ DESCRIPTION πŸ¦‹

so gmx dssp lets you use the DSSP algorithm β€” and like, specifically the way it detects particular patterns of hydrogen bonds between amino acid residues β€” to figure out the secondary structure of a protein. it's giving structural biology main character energy fr fr πŸ’œ

here are the one-symbol secondary structure designations that show up in the output file, and you're gonna want to know all of them:

  • 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 is your bestie for getting 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 (that's the "gromacs" option) and using hydrogen pseudo-atoms based on C and O atom coordinates of previous residue (that's the "dssp" option). and listen β€” you should always use the "dssp" option for structures with absent hydrogen atoms! no cap πŸ™ˆ


-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 the GROMACS neighbor-search method to find residue pairs that may have a hydrogen bond, instead of simply iterating over the residues among themselves. lowkey a whole different vibe πŸ‘‘


-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. periodt πŸŽ€


-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.


-polypro enables the search for polyproline helices (default behavior, equivalent to DSSP v4). disabling this option will result in disabling the search for polyproline helices, reproducing the behavior of DSSP v2 🌈


⚠️ okay important one: 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. it's not her fault, it's the cell size's fault πŸ’€
⚠️ also please note that the computation is always done in single precision, regardless of the precision for which GROMACS was configured. single precision only, no exceptions, iykyk 😭

✨ OPTIONS ✨

πŸŽ€ Options to specify input files:

  • -f [<.xtc/.trr/...>] (traj.xtc) (Optional) β€” input trajectory or single configuration: xtc trr cpt gro g96 pdb tng h5md
  • -s [<.tpr/.gro/...>] (topol.tpr) (Optional) β€” input structure: tpr gro g96 pdb brk ent
  • -n [<.ndx>] (index.ndx) (Optional) β€” extra index groups

🧚 Options to specify output files:

  • -o [<.dat>] (dssp.dat) β€” filename for DSSP output
  • -num [<.xvg>] (num.xvg) (Optional) β€” output file name for secondary structures statistics for the trajectory

πŸ‘‘ Other options:

  • -b <time> (0) β€” first frame (ps) to read from trajectory
  • -e <time> (0) β€” last frame (ps) to read from trajectory
  • -dt <time> (0) β€” only use frame if t MOD dt == first time (ps)
  • -tu <enum> (ps) β€” unit for time values: fs, ps, ns, us, ms, s
  • -fgroup <selection> β€” atoms stored in the trajectory file (if not set, assume first N atoms)
  • -xvg <enum> (xmgrace) β€” plot formatting: xmgrace, xmgr, none
  • -[no]rmpbc (yes) β€” make molecules whole for each frame
  • -[no]pbc (yes) β€” use periodic boundary conditions for distance calculation
  • -sf <file> β€” provide selections from files
  • -selrpos <enum> (atom) β€” 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
  • -seltype <enum> (atom) β€” 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
  • -sel <selection> β€” group for DSSP
  • -hmode <enum> (gromacs) β€” hydrogens pseudoatoms creating mode: gromacs, dssp
  • -hbond <enum> (energy) β€” selects between different definitions of hydrogen bond: energy, geometry
  • -[no]nb (yes) β€” use GROMACS neighbor-search method
  • -cutoff <real> (0.9) β€” 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
  • -ppstretch <enum> (default) β€” stretch value for PP-helices: shortened, default
  • -[no]polypro (yes) β€” perform a search for polyproline helices

🌈 SEE ALSO 🌈

gmx(1)

more information about GROMACS is available at <http://www.gromacs.org/> πŸ’–


2026, GROMACS development team


and that's the whole page bestie, you ate that πŸ’…βœ¨ this was GMX-DSSP(1) from the GROMACS manual, version 2026.4, dated Oct 02, 2026 β€” go calculate some secondary structure, you absolute icon πŸ‘‘πŸŽ€

πŸ“¦ where this page comes from