You cannot accurately compute binding FREE energy only by re-evaluating your MD trajectories. To estimate binding free energy of two molecules in explicit solvent, you need to carry out enhanced sampling simulation methods, either umbrella sampling or free energy perturbation. For details on how to use these methods with GROMACS, please refer to the following link: