I used QM\MM method to optimize the protein-ligand complex. protein was in the MM region and ligand is in the QM region. After optimization, I want to how to calculate the interaction energy between ligand and protein. (include the electrostatic and van der Waals interaction energies) Is anyone can help me?