I have used MM/PBSA with a similar system and obtained a similarly huge
binding free energy. However, I was interested in how mutations affected the
binding free energy, the "delta-delta-G". This ddG is what's
"quantitatively" correct, not, usually, the dG as I believe you have
reported below.