Our previously derived models of the active state of the β2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The β2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures. The general features of the microswitches were well reproduced, but with minor differences, partly because of the unexpected X-ray results for the rotamer toggle switch. In addition, most of the interacting residues between the receptor and the G-protein were identified. This analysis of the modelling has also given important additional insight into GPCR dimerization: re-analysis of results on photoaffinity analogues of rhodopsin provided additional evidence that TM4 (transmembrane helix 4) resides at the dimer interface and that ligands such as bivalent ligands may pass between the mobile helices. A comparison, and discussion, is also carried out between the use of implicit and explicit solvent for active-state modelling.

In recent years there has been an increasing number of reports describing G protein-coupled receptor (GPCR) dimerization and heterodimerization. However, the evidence on the nature of the dimers and their role in GPCR activation is inconclusive. Consequently, we present here a review of our computational studies on G protein-coupled receptor dimerization and domain swapping. The studies described include molecular dynamics simulations on receptor monomers and dimers in the absence of ligand, in the presence of an agonist, and in the presence of an antagonist (or more precisely an inverse agonist). Two distinct sequence-based approaches to studying protein interfaces are also described, namely correlated mutation analysis and evolutionary trace analysis. All three approaches concur in supporting the proposal that the dimerization interface includes transmembrane helices 5 and 6. These studies cannot distinguish between domain swapped dimers and contact dimers as the models used were restricted to the helical part of the receptor. However, it is proposed that for the purpose of signalling, the domain swapped dimer and the corresponding contact dimer are equivalent. The evolutionary trace analysis suggests that every GPCR family and subfamily (for which sufficient sequence data is available) has the potential to dimerize through this common functional site on helices 5 and 6. The evolutionary trace results on the G protein are briefly described and these are consistent with GPCR dimerization. In addition to the functional site on helices 5 and 6, the evolutionary trace analysis identified a second functional site on helices 2 and 3. Possible roles for this site are suggested, including oligomerization.