Despite a considerable diversity in size, sequence, secondary, and tertiary structure of the nine peptide ligands, we observed a significant overlap in the receptor region they bind to, particularly in a binding pocket between the outer leaflet portions of transmembrane helices (
Figure 2). Using Rosetta [
70,
71], we calculated the per-residue ΔΔG of the interacting residues on the transmembrane helices, two conserved ECL1 residues (23.49 and 23.50), and three conserved ECL2 residues (45.50, 45.51, and 45.52). The details of the structure optimization and ΔΔG analysis protocols are listed in the
Supplementary Material, and the ΔΔG values of all residues are listed in the
Supplement Table S1. To make the optimization and ΔΔG analysis possible for the apelin/ApelinR (PDB ID: 5VBL) and the PMX53/C5aR (PDB ID: 6C1Q) structures, which contain a non-nature peptide backbone, we generated 5VBL* and 6C1Q* as natural-backbone peptide analogs of those structures. More specifically, the 5VBL* peptide ligand has the native apelin sequence, and the covalent bond between ornithine (ORN) at position 2 and the N-terminal acetyl group is omitted in the 6C1Q* peptide (
Figure S3). Then, the GPCR residues were ranked based on their calculated ΔΔG. We selected 14 common residues with ΔΔG of less than −1 and contact peptide ligands in at least seven out of nine GPCR-peptide complexes. The details of the list and their locations on a GPCR structure are mapped in the structure of the ET-1/ETB receptor complex, as shown in
Figure 6. This list of the top 14 residues implies a potential common peptide-binding mechanism among class A GPCRs. This common binding pocket encompasses two residues of TM2 (2.60 and 2.63), one from TM3 (3.32), three from TM6 (6.51, 6.55, and 6.58), five from TM7 (7.28, 7.32, 7.35, 7.36, and 7.39), and all three conserved ECL2 residues. More specifically, the common peptide engagement mechanism starts from the end of the β-hairpin of ECL2, extends to the tip of TM2, touches the extracellular half of TM7 and TM6, then ends at the core of TM3. The
Supplementary Table S2 summarizes the non-Van Der Waal interactions between these 14 residues and the corresponding peptides. Although additional structures of peptide-GPCR complexes are still needed to validate our hypothesis of the common peptide binding pocket, this finding could help guide future structural studies of this family of GPCRs.
Herein, we examine whether the common binding mechanism agrees with the models of three class A GPCRs—Y
1 [
56], Y
2 [
14], and Ghrelin receptor [
72]—and their endogenous peptide ligands—NPY and Ghrelin. In those studies, the peptide docking experiments were conducted using FlexPepDock [
73] with constraints from mutagenesis, cross-linking, and NMR data. For each complex, the ΔΔG analysis was performed on an ensemble of docking models. The per-residue ΔΔG values were assigned to the interacting residues of the GPCR targets. The peptides’ binding pockets contain all of the 14 common residues, except for ghrelin, which does not contact the residue 7.36. Furthermore, most of the interactions between the common residues and NPY or ghrelin are favorable or at least neutral, except for the high ΔΔG value of residue 7.32 from Y
2 (
Figure 8). These results imply that the observation of the common peptide engagement pocket can also be applied to the docking study of peptide class A GPCRs, especially with limited experimental data.
Figure 8. Models of peptide/class A GPCR complexes show that the peptides interact with the top 14 common residues. (From
left to
right): A table lists the ΔΔGs values of the 14 common residues of Y
1 [
56], Y
2 [
14], and ghrelin receptors [
72], as well as their sum and average values. The absence of the ΔΔG values indicates that the corresponding residues do not interact with the peptide ligands. The residue ΔΔG cells are colored based on the ΔΔG values (negative: Blue, neutral: White, and positive: Red). The blank cells indicate that the residues do not interact with the peptide ligands. Models of NPY (cyan) bind with the Y
1 receptor (grey) and the Y
2 receptor (orange), and ghrelin (magenta) binds with the ghrelin receptor (green).
A GPCR pharmacogenomics study has extracted polymorphism data for the coding-region of the 108 GPCR drug targets [
74]. From the data provided by the authors, we found around 30 relevant GPCR mutants that were predicted to be deleterious by the sorting intolerant from tolerant (SIFT) [
75] or Polyphen [
76]. Those 30 genetic invariants have population allele frequencies of around 1 to 28 over 120,000 individuals and are related to the shared peptide interacting residues or are close to those residues. The table containing the information regarding the relevant mutants of peptide and protein binding class A GPCRs is summarized in the
Supplementary Table Peptide_binding_pocket_genetic_variants.xlsx. The data suggest the great potential of the proposed common peptide-binding pocket as drug targets for class A GPCRs.