求助，使用Gromacs得到的分子动力学模拟的结果如何进行蛋白-配体相互作用分数分析？

郑素娟

各位老师好，本人在使用Gromacs进行分子动力学模拟时，如何能得到如下图所示的蛋白质氨基酸残基与小分子的相互作用分数图呢？感谢分享经验！

student0618

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.

低调的板凳

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

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

rpestana94

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