gmx chi [-s [<.gro/.g96/...>]] [-f [<.xtc/.trr/...>]] [-ss [<.dat>]]
[-o [<.xvg>]] [-p [<.pdb>]] [-jc [<.xvg>]] [-corr [<.xvg>]]
[-g [<.log>]] [-ot [<.xvg>]] [-oh [<.xvg>]] [-rt [<.xvg>]]
[-cp [<.xvg>]] [-b <time>] [-e <time>] [-dt <time>] [-[no]w]
[-xvg <enum>] [-r0 <int>] [-[no]phi] [-[no]psi] [-[no]omega]
[-[no]rama] [-[no]viol] [-[no]periodic] [-[no]all] [-[no]rad]
[-[no]shift] [-binwidth <int>] [-core_rotamer <real>]
[-maxchi <enum>] [-[no]normhisto] [-[no]ramomega]
[-bfact <real>] [-[no]chi_prod] [-[no]HChi] [-bmax <real>]
[-acflen <int>] [-[no]normalize] [-P <enum>] [-fitfn <enum>]
[-beginfit <real>] [-endfit <real>]
gmx chi computes phi, psi, omega, and chi dihedrals for all your amino acid backbone and sidechains. It can compute dihedral angle as a function of time, and as histogram distributions. The distributions (histo-(dihedral)(RESIDUE).xvg) are cumulative over all residues of each type.
If option -corr is given, the program will calculate dihedral autocorrelation functions. The function used is C(t) = <cos(chi(tau)) cos(chi(tau+t))>. The use of cosines rather than angles themselves, resolves the problem of periodicity. (Van der Spoel & Berendsen (1997), Biophys. J. 72, 2032-2041). Separate files for each dihedral of each residue (corr(dihedral)(RESIDUE)(nresnr).xvg) are output, as well as a file containing the information for all residues (argument of -corr).
With option -all, the angles themselves as a function of time for each residue are printed to separate files (dihedral)(RESIDUE)(nresnr).xvg. These can be in radians or degrees.
A log file (argument -g) is also written. This contains
- information about the number of residues of each type.
- The NMR ^3J coupling constants from the Karplus equation.
- a table for each residue of the number of transitions between rotamers per nanosecond, and the order parameter S^2 of each dihedral.
- a table for each residue of the rotamer occupancy.
All rotamers are taken as 3-fold, except for omega and chi dihedrals to planar groups (i.e. chi_2 of aromatics, Asp and Asn; chi_3 of Glu and Gln; and chi_4 of Arg), which are 2-fold. “rotamer 0” means that the dihedral was not in the core region of each rotamer. The width of the core region can be set with -core_rotamer
The S^2 order parameters are also output to an .xvg file (argument -o ) and optionally as a .pdb file with the S^2 values as B-factor (argument -p). The total number of rotamer transitions per timestep (argument -ot), the number of transitions per rotamer (argument -rt), and the ^3J couplings (argument -jc), can also be written to .xvg files. Note that the analysis of rotamer transitions assumes that the supplied trajectory frames are equally spaced in time.
If -chi_prod is set (and -maxchi > 0), cumulative rotamers, e.g. 1+9(chi_1-1)+3(chi_2-1)+(chi_3-1) (if the residue has three 3-fold dihedrals and -maxchi >= 3) are calculated. As before, if any dihedral is not in the core region, the rotamer is taken to be 0. The occupancies of these cumulative rotamers (starting with rotamer 0) are written to the file that is the argument of -cp, and if the -all flag is given, the rotamers as functions of time are written to chiproduct(RESIDUE)(nresnr).xvg and their occupancies to histo-chiproduct(RESIDUE)(nresnr).xvg.
The option -r generates a contour plot of the average omega angle as a function of the phi and psi angles, that is, in a Ramachandran plot the average omega angle is plotted using color coding.
Options to specify input files:
Options to specify output files:
Other options: