We present a comprehensive computational analysis of available 3D GPCR-G-protein complexes to inspect the structural determinants of G-protein-coupling selectivity.
Analysis of the residue contacts at interaction interfaces has revealed a network of secondary structure elements recapitulating known structural hallmarks determining G-protein-coupling specificity, including TM5, TM6 and ICLs. We coded interface contacts into generic-number fingerprints to reveal specific coupling-determinant positions. Clustering of Gs vs Gi complexes is best achieved when considering both GPCR and G-protein contacting residues rather than separated representations of the interaction partners, suggesting that coupling specificity emerges as contextual residue interactions at the interface. Interestingly, Gs-GPCR complexes contain a higher number of contacts than Gi/o-GPCR complexes, likely caused by overall higher conservation and structural constraint on the Gs interface. In contrast, Gi/o proteins adopt a wider number of alternative docking poses on cognate receptors, as assessed via structural alignments of representative 3D complexes. Furthermore, binding energy calculations demonstrate that distinct structural properties of the complexes contribute to higher stability of Gs than Gi/o complexes. AlphaFold2 predictions of experimental binary complexes confirmed several of these structural features and allowed us to augment the structural coverage of poorly characterized complexes (e.g. G12/13).
We propose that the structural properties of different G-protein complexes, such as structural restraining of Gs compared to Gi/o ones, could be instrumental in fine-tuning their activation and downstream signaling mechanisms.