Prediction of heterodimerization interfaces of G-protein coupled receptors with a new subtractive correlated mutation method

Recent studies employing differential epitope tagging, selective immunoprecipitation of receptor complexes and fluorescence or bioluminescence resonance energy transfer techniques provide direct evidence for heterodimerization between both closely and distantly related members of the G-protein coupled receptor (GPCR) family. Since, heterodimerization appears to play a role in modulating agonist affinity, efficacy and/or trafficking properties, the molecular models of GPCRs required to understand receptor function must consider these oligomerization hypotheses. To advance knowledge in this field, we present here a computational approach based on correlated mutation analysis and the structural information contained in three-dimensional molecular models of the transmembrane regions of GPCRs built using the rhodopsin crystal structure as a template. The new subtractive correlated mutation method reveals likely heterodimerization interfaces amongst the different alternatives for the positioning of two tightly packed bundles of seven transmembrane domains next to each other in contact heterodimers of GPCRs. Predictions are applied to GPCRs in the class of opioid receptors. However, in the absence of a known structure of any GPCR dimer, the features of the method and predictions are also illustrated and analyzed for a dimeric complex of known structure.