Structure and ligand-binding site characteristics of the human P2Y11 nucleotide receptor deduced from computational modelling and mutational analysis

Abstract
The P2Y11-R (P2Y11 receptor) is a less explored drug target. We computed an hP2Y11-R (human P2Y11) homology model with two templates, bovine-rhodopsin (2.6 Å resolution; 1 Å=0.1 nm) and a hP2Y1–ATP complex model. The hP2Y11-R model was refined using molecular dynamics calculations and validated by virtual screening methods, with an enrichment factor of 5. Furthermore, mutational analyses of Arg106, Glu186, Arg268, Arg307 and Ala313 confirmed the adequacy of our hP2Y11-R model and the computed ligand recognition mode. The E186A and R268A mutants reduced the potency of ATP by one and three orders of magnitude respectively. The R106A and R307A mutants were functionally inactive. We propose that residues Arg106, Arg268, Arg307 and Glu186 are involved in ionic interactions with the phosphate moiety of ATP. Arg307 is possibly also H-bonded to N6 of ATP via the backbone carbonyl. Activity of ATP at the F109I mutant revealed that the proposed π-stacking of Phe109 with the adenine ring is a minor interaction. The mutation A313N, which is part of a hydrophobic pocket in the vicinity of the ATP C-2 position, partially explains the high activity of 2-MeS-ATP at P2Y1-R as compared with the negligible activity at the P2Y11-R. Inactivity of ATP at the Y261A mutant implies that Tyr261 acts as a molecular switch, as in other G-protein-coupled receptors. Moreover, analysis of cAMP responses seen with the mutants showed that the efficacy of coupling of the P2Y11-R with Gs is more variable than coupling with Gq. Our model also indicates that Ser206 forms an H-bond with Pγ (the γ-phosphate of the triphosphate chain of ATP) and Met310 interacts with the adenine moiety.