The theoretical calculation of protein-ligand binding free energy is a grand challenge in computational biology.Accurate prediction of critical residues along with their specific and quantitative contributions to protein-ligand binding free energy is extremely helpful to reveal binding mechanisms and identify drug-like molecules that alter protein-ligand interactions.