[gmx-users] Hydrogen bonds between new residue and protein.

[Dawid das]

Dear Gromacs experts,

This time I would like to use g_hbond to calculate hydrogen bonds between
one residue of the protein and all other residues in the protein. The thing
is that the former residue is a new one, defined by me from scratch. Now
when I run g_hbond and choose CH6 and Protein index groups (CH6 is my new
residue and Protein also includes this CH6 residue) I got 0.0 for every
time step although I strongly believe there are couple of hydrogen bonds.

I was thinking that the issue may be that CH6 is built of atomtypes which
are not recognised by g_hbond as ones able to form hydrogen bonds. Is it
possible explanation? If yes, what can I do? Is recompiling what I need to
do?

Thanks,

Dawid Grabarek