A general computational scheme to evaluate the effects of single point mutations on ligand binding is reported. This scheme is applied to characterize agonist binding to the A_2A adenosine receptor, and is found to accurately explain how point mutations of different nature affect the binding affinity of a potent agonist.