Maybe you can try including indexs of those atoms in the index file with different names and pass it with gmx hbond command.

Dear Sara Ahmadi

GROMACS prints out the trajectory of the system. As long as you set an appropriate output frequency and you save the positions of the groups you are interested in, then it is just a matter of analysis.

For hbonds, take a look at the tool gmx hbond.

https://manual.gromacs.org/documentation/5.1/onlinehelp/gmx-hbond.html

The selection is, as always, coming from an index file. If you write the indexes for each OH group of your sugars then you can analyse the hbonding properties of each one of them separately.

Hope this helps,

Nicola