“…S1P1 was used as a template due to its 33% sequence homology (similarity) with GPR6 and because both share essential features such as (1) the absence of helix kinking proline residues in TMH2 (2.58 or 2.59) and TMH5 (5.50); (2) both possess an acidic residue E1.49 before the highly conserved N1.50 in TMH1 (most Class A GPCRs have a G1.49 here); and (3) the presence of an internal disulfide bridge in the EC2 loop. Neither the Lysophosphatidic acid receptor 1 (LPA1) [PDB identifier: 4Z35] [36] or CB1 [PDB ID: 5TGZ and 5U09] [58,59], the other closely related crystallized GPCRs, were used because the LPA1 shares less essential structural features than S1P1 and the CB1 structures suffer from crystal packing problems [60]. Before mutating, the sequence of GPR6 was aligned first with the sequences of other class A GPCRs such as Rhodopsin (Rho), CB1, CB2, δ-opioid receptor (DOR), β2-AR, LPA1 and S1P1 receptors using the highly conserved residues/motifs as alignment guides (N1.50, D2.50, R3.50, W4.50, Y5.58, CWXP in TMH6, and NPXXY in TMH7) (blue residues in Figure 1).…”