Human P2Y receptors encompass at least eight subtypes of Class A G protein-coupled receptors (GPCRs), responding to adenine and/or uracil nucleotides. Using a BLAST search against the Homo sapiens subset of the SWISS-PROT and TrEMBL databases, we identified 68 proteins showing high similarity to P2Y receptors. To address the problem of low sequence identity between rhodopsin and the P2Y receptors, we performed a multiple-sequence alignment of the retrieved proteins and the template bovine rhodopsin, combining manual identification of the transmembrane domains (TMs) with automatic techniques. The resulting phylogenetic tree delineated two distinct subgroups of P2Y receptors: Gq-coupled subtypes (e.g., P2Y1) and those coupled to Gi (e.g., P2Y12). On the basis of sequence comparison we mutated three Tyr residues of the putative P2Y1 binding pocket to Ala and Phe and characterized pharmacologically the mutant receptors expressed in COS-7 cells. The mutation of Y306 (7.35, site of a cationic residue in P2Y12) or Y203 in the second extracellular loop selectively decreased the affinity of the agonist 2-MeSADP, and the Y306F mutation also reduced antagonist (MRS2179) affinity by 5-fold. The Y273A (6.48) mutation precluded the receptor activation without a major effect on the ligand-binding affinities, but the Y273F mutant receptor still activated G proteins with full agonist affinity. Thus, we have identified new recognition elements to further define the P2Y1 binding site and related these to other P2Y receptor subtypes. Following sequence-based secondary-structure prediction, we constructed complete models of all the human P2Y receptors by homology to rhodopsin. Ligand docking on P2Y1 and P2Y12 receptor models was guided by mutagenesis results, to identify the residues implicated in the binding process. Different sets of cationic residues in the two subgroups appeared to coordinate phosphate-bearing ligands. Within the P2Y1 subgroup these residues are R3.29, K/R6.55, and R7.39. Within the P2Y12 subgroup, the only residue in common with P2Y1 is R6.55, and the role of R3.29 in TM3 seems to be fulfilled by a Lys residue in EL2, whereas the R7.39 in TM7 seems to be substituted by K7.35. Thus, we have identified common and distinguishing features of P2Y receptor structure and have proposed modes of ligand binding for the two representative subtypes that already have well-developed ligands.