Nature of H bonding in gromacs

User discussions
,
1

GROMACS version:2023
GROMACS modification: Yes/No

What is the nature of H-bonding in GROMACS ? Is it primarily electrostatic , or does it calculate the number of H-bonds based on the angle and distance criteria?
Can i relate the electrostatic interaction to H bonding?

(I used the CHARMM36 force field for my MD simulation)

2

It is based on angle and distance (https://manual.gromacs.org/current/reference-manual/analysis/hydrogen-bonds.html).

There isn’t any command that relates the h-bonds to electrostatic interactions directly (at least not one I could find), but what you could do is use the hbond -hbn flag to retrieve an index file of all the atoms involved with hydrogen bonding, as is done in this tutorial (How to study Hydrogen bonds using GROMACS - Compchems) than use that index file to specify the atoms to be used as the energy groups when extracting information on electrostatic interactions using gmx energy, similar to the analysis done in this tutorial (Protein-Ligand Complex)

3

Thanks for your valuable feedback

1 Like
4

Fig 59. shows a water bridge hydrogen bond, but without any explicit statements.
Can gmx hbond calculate water-bridged hydrogen bonds?

5

I haven’t been able to find anything in the documentation for gmx hbond for identifying water bridges specifically, but the python library mdanalysis does have an analysis module for water bridge analysis: 4.3.3. Water Bridge analysis — MDAnalysis.analysis.hydrogenbonds.WaterBridgeAnalysis — MDAnalysis 2.7.0 documentation