Conventional MMPBSA is not suitable as the forcefield will vary in amber. I do not have access to CHARMM software in our lab. Is there any other way I could compute the strength of interaction between my protein and ligand along the trajectory apart from running 3 explicit simulations for protein, ligand and protein-ligand?