Nucleotides are emerging as an ubiquitous family of extracellular signaling molecules. These effects are mediated through a specific class of plasma membrane receptors called P2 receptors that, according to the molecular structure, are further subdivided into two subfamilies: P2Y and P2X. Specifically, P2X-receptors are ligand-gated ion channels, whereas P2Y-receptors belong to the superfamily of G-protein-coupled receptors. In this review, we focus our attention to GPCRs molecular architecture, with the special emphasis on our work on the human P2Y1 receptor. In fact, despite an enormous amount of research on the structure and function of these receptors, fundamental understanding of the molecular details of ligand / GPCR interactions remains very rudimentary. How agonist binding transforms a resting GPCR into its active form and the microscopic basis of binding site blockade by an antagonist are generally still unclear. In the absence of high-resolution structural knowledge of GPCRs, such questions only can be addressed by building models, which are tested through pharmacological and biochemical studies. In this review, we underline how different molecular modeling approaches can help the investigation of both receptor architecture and ligand / receptor molecular recognition.