I have simulated multiple systems (dimer+ligand) for 200 ns and trying to calculate the MMPBSA binding energy using g_mmpbsa. However, I am getting positive binding energy in each case. The ligand is occupying the dimer interface throughout the simulation. I have used multiple snapshots as well, but the issue remains the same. Kindly suggest some solutions.