The theoretical calculation of protein-protein binding free energy is a grand challenge in computational biology. Accurate prediction of critical residues along with their specific and quantitative contributions to protein-protein binding free energy is extremely helpful to reveal binding mechanisms and identify drug-like molecules that alter protein-protein interactions. In this paper, we propose an interaction entropy approach combined with the molecular mechanics/generalized Born surface area (MM/GBSA) method for solvation to compute residue-specific protein-protein binding free energy. In the current approach, the entropic loss in binding free energy of individual residues is explicitly computed from moledular dynamics (MD) simulation by using the interaction entropy method. In this approach the entropic contribution to binding free energy is determined from fluctuation of the interaction in MD simulation. Studies for an extensive set of realistic protein-protein interaction systems showed that by including the entropic contribution, the computed residue-specific binding free energies are in better agreement with the corresponding experimental data.