Molecular modeling of agonist binding to the human A(2A) adenosine receptor (AR) was assessed and extended in light of crystallographic structures. Heterocyclic adenine nitrogens of cocrystallized agonist overlaid corresponding positions of the heterocyclic base of a bound triazolotriazine antagonist, and ribose moiety was coordinated in a hydrophilic region, as previously predicted based on modeling using the inactive receptor. Automatic agonist docking of 20 known potent nucleoside agonists to agonist-bound A(2A)AR crystallographic structures predicted new stabilizing protein interactions to provide a structural basis for previous empirical structure activity relationships consistent with previous mutagenesis results. We predicted binding of novel C2 terminal amino acid conjugates of A(2A)AR agonist CGS21680 and used these models to interpret effects on binding affinity of newly synthesized agonists. d-Amino acid conjugates were generally more potent than l-stereoisomers and free terminal carboxylates more potent than corresponding methyl esters. Amino acid moieties were coordinated close to extracellular loops 2 and 3. Thus, molecular modeling is useful in probing ligand recognition and rational design of GPCR-targeting compounds with specific pharmacological profiles.