I have run a 30 ns MD of receptor -ligand complex. Now I want to calculate the distance a particular non-bond interaction (hydrophoibc) over the simulation period, so how can I don that? Finally how can I calculate the MMGBSA or MMPBSA binding energy?