Binding free energy

GROMACS version:
GROMACS modification: Yes/No
Dear Gromacs user
I need your kind suggestion in setting up a relative free energy calculation for mutation of an amino acid from Arginine to Histidine in the presence of cofactor.
Could you please suggest how should I work with topology and thermodynamic cycle.
Thanks in advance.