A computational method for the prediction of binding free energy between a protein receptor and its ligands