Purinergic receptors are membrane proteins that regulate numerous cellular functions by catalyzing reactions involving purine nucleotides or nucleosides. Among the three receptor families, i.e., P1, P2X, and P2Y, the P1 and P2Y receptors share common structural features of class A GPCR. Comprehensive sequence and structural analysis revealed that the P1 and P2Y receptors belong to two distinct groups. They exhibit different ligand-binding site features that can distinguish between specific activators. These specific amino acid residues in the binding cavity may be involved in the selectivity and unique pharmacological behavior of each subtype. In this study, we conducted a structure-based analysis of purinergic P1 and P2Y receptors to identify their evolutionary signature and obtain structural insights into ligand recognition and selectivity. The structural features of the P1 and P2Y receptor classes were compared based on sequence conservation and ligand interaction patterns. Orthologous protein sequences were collected for the P1 and P2Y receptors, and sequence conservation was calculated based on Shannon entropy to identify highly conserved residues. To analyze the ligand interaction patterns, we performed docking studies on the P1 and P2Y receptors using known ligand information extracted from the ChEMBL database. We analyzed how the conserved residues are related to ligand-binding sites and how the key interacting residues differ between P1 and P2Y receptors, or between agonists and antagonists. We extracted new similarities and differences between the receptor subtypes, and the results can be used for designing new ligands by predicting hotspot residues that are important for functional selectivity.