I was performing MD simulation for protein-ligand complex following GROMACS tutorial.
I compared the energy of two ligands (S1 and S2) to a protein, using GROMACS LJ-SR (Lennard-Jones Short Range) and Coul-SR (Coulombic Short Range) option.
Below is the plot.
As you can see the LJ-SR is complete opposite of Coul-SR. Are they expected to be the opposite? Which is more appropriate choice to represent real binding energy?