Fix g_chi -omega
authorMark Abraham <mark.j.abraham@gmail.com>
Mon, 29 Apr 2013 22:48:24 +0000 (00:48 +0200)
committerGerrit Code Review <gerrit@gerrit.gromacs.org>
Wed, 22 May 2013 13:11:45 +0000 (15:11 +0200)
commitb0dd1667859afa216a08319d23c9e560c3d149ee
tree9ba250466f5aad58fa8455372f5903b20f497f67
parent9dc3fe4ab0329965bf093661572ff3a4306b4aa3
Fix g_chi -omega

Per IUPAC, we should calculate peptide omega angles from
Calpha-C-N-Calpha, and there are only n-1 such angles defined for a
polypeptide of n amino acid residues.

There are some other IUPAC violations in make_chi_ind(), but solving
those is deferred to 5.0, not least because there will be a
behavioural change when we conform to IUPAC.

The problem that led to the report had several causes. The hard-coding
of the string names for hydrogen bound to peptide nitrogen (which in
CHARMM27 is "HN") meant that an H atom was not found. has_dihedral()
then checked whether a Calpha was found, which was true. Then
make_chi_ind() computed the angle based on H without doing any further
check. So some garbage coordinate was looked up when computing the
angle.

A proper fix for this kind of stuff for 5.0 requires finding all the
places we try to guess what atoms are from their names and using
properly implemented selections.

Also removed orphaned header file.

Fixes #953, details and links there

Change-Id: I63257a1d6d8ed46d251b8bc0cb308dfbd2864a97
include/gstat.h
src/tools/dlist.c
src/tools/gmx_chi.c
src/tools/pp2shift.h [deleted file]