You can try to use ProLIF, a python library to get different types of interaction and the frames where the interaction was present, with that you can make something similar to the picture

看图应该是某商业软件的MD分析结果，对应不同颜色数据时氢键，疏水作用，水桥，离子相互作用。

不过对于gmx产生的轨迹，目前没有什么特别方便的方法直接拿到这种数据。氢键的话你可以通过gmx hbond或者VMD对于不同残基分组后进行统计。其他的相互作用需要比较复杂的脚本才能分析，比较麻烦。 如果对于残基作用想进行仔细分析的话，可以考虑楼上的方案，算MMPBSA，之后直接分析残基的能量。 而不纠结具体相互作用数量的统计。

Probably trajectory analysis to count contacts and H-bonds etc. Then get the fraction with some data analysis tools. Not sure what type of interaction represented by each color.

Another way to do something similar is to do energy analysis with e.g. gmx_MMPBSA (https://valdes-tresanco-ms.github.io/gmx_MMPBSA) and use the per-residue energy decomposition option. This would tell the energy contribution of each residue.