Latrophilin fragments behave as independent proteins that associate and signal on binding of LTX(N4C)

Heptahelical, or G-protein-coupled, receptors control many cellular functions and normally consist of one polypeptide chain. In contrast, heptahelical receptors that belong to the long N-terminus, group B (LNB) family are cleaved constitutively into two fragments. The N-terminal fragments (NTFs) resemble cell-adhesion proteins and the C-terminal fragments (CTFs) are typical G-protein-coupled receptors (GPCRs) with seven transmembrane regions. However, the functional roles of this cleavage and of any subsequent NTF-CTF interactions remain to be identified. Using latrophilin, a well-studied member of the LNB family, we now demonstrate that cleavage is critical for delivery of this receptor to the cell surface. On the plasma membrane, NTF and CTF behave as separate membrane proteins involved, respectively, in cell-surface reception and signalling. The two fragments can also internalise independently. However, separated NTF and CTF can re-associate on solubilisation. Agonist binding to NTF on the cell surface also induces re-association of fragments and provokes signal transduction via CTF. These findings define a novel principle of structural and functional organisation of the cleaved, two-subunit GPCRs.     

Novel human G protein-coupled receptors with long N-terminals containing GPS domains and Ser/Thr-rich regions

We report eight novel members of the superfamily of human G protein-coupled receptors (GPCRs) found by searches in the human genome databases, termed GPR97, GPR110, GPR111, GPR112, GPR113, GPR114, GPR115 and GPR116. Phylogenetic analysis shows that these are additional members of a family of GPCRs with long N-termini, previously termed EGF-7TM, LNB-7TM, B2 or LN-7TM. Five of the receptors form their own phylogenetic cluster, while three others form a cluster with the previously reported HE6 and GPR56 (TM7XN1). All the receptors have a GPS domain in their N-terminus and long Ser/Thr-rich regions forming mucin-like stalks. GPR113 has a hormone binding domain and one EGF domain. GPR112 has over 20 Ser/Thr repeats and a pentraxin domain. GPR116 has two immunoglobulin-like repeats and a SEA box. We found several human EST sequences for most of the receptors showing differential expression patterns, which may indicate that some of these receptors participate in reproductive functions while others are more likely to have a role in the immune system.     

Functional and biophysical analysis of the C-terminus of the CGRP-receptor; a family B GPCR

G-protein coupled receptors (GPCRs) typically have a functionally important C-terminus which, in the largest subfamily (family A), includes a membrane-parallel eighth helix. Mutations of this region are associated with several diseases. There are few C-terminal studies on the family B GPCRs and no data supporting the existence of a similar eighth helix in this second major subfamily, which has little or no sequence homology to family A GPCRs. Here we show that the C-terminus of a family B GPCR (CLR) has a disparate region from N400 to C436 required for CGRP-mediated internalization, and a proximal region of twelve residues (from G388 to W399), in a similar position to the family A eighth helix, required for receptor localization at the cell surface. A combination of circular and linear dichroism, fluorescence and modified waterLOGSY NMR spectroscopy (SALMON) demonstrated that a peptide mimetic of this domain readily forms a membrane-parallel helix anchored to the liposome by an interfacial tryptophan residue. The study reveals two key functions held within the C-terminus of a family B GPCR and presents support for an eighth helical region with striking topological similarity to the nonhomologous family A receptor. This helix structure appears to be found in most other family B GPCRs.     

There exist at least 30 human G-protein-coupled receptors with long Ser/Thr-rich N-termini

We report six novel members of the superfamily of human G-protein coupled receptors (GPCRs) found by searches in the human genome databases, termed GPR123, GPR124, GPR125, GPR126, GPR127, and GPR128. Phylogenetic analysis demonstrates that these are additional members of the family of GPCRs with long N-termini, previously termed EGF-7TM, LNB-7TM, B2 or LN-7TM, showing that there exist at least 30 such GPCRs in the human genome. Three of these receptors form their own phylogenetic cluster, while two other places in a cluster with the previously reported HE6 and GPR56 (TM7XN1) and one with EMR1-3. All the novel receptors have a GPS domain in their N-terminus, except GPR123, as well as long Ser/Thr rich regions forming mucin-like stalks. GPR124 and GPR125 have a leucine rich repeat (LRR), an immunoglobulin (Ig) domain, and a hormone-binding domain (HBD). The Ig domain shows similarities to motilin and titin, while the LRR domain shows similarities to LRIG1 and SLIT1-2. GPR127 has one EGF domain while GPR126 and GPR128 do not contain domains that are readily recognized in other proteins beyond the GPS domain. We found several human EST sequences for most of the receptors showing differential expression patterns, which may indicate that some of these receptors participate in central functions while others are more likely to have a role in the immune or reproductive systems.     

Elucidation of the role of peptide linker in calcium-sensing receptor activation process

Family 3 G-protein-coupled receptors (GPCRs), which includes metabotropic glutamate receptors (mGluRs), sweet and "umami" taste receptors (T1Rs), and the extracellular calcium-sensing receptor (CaR), represent a distinct group among the superfamily of GPCRs characterized by large amino-terminal extracellular ligand-binding domains (ECD) with homology to bacterial periplasmic amino acid-binding proteins that are responsible for signal detection and receptor activation through as yet unresolved mechanism(s) via the seven-transmembrane helical domain (7TMD) common to all GPCRs. To address the mechanism(s) by which ligand-induced conformational changes are conveyed from the ECD to the 7TMD for G-protein activation, we altered the length and composition of a 14-amino acid linker segment common to all family 3 GPCRs except GABA(B) receptor, in the CaR by insertion, deletion, and site-directed mutagenesis of specific highly conserved residues. Small alterations in the length and composition of the linker impaired cell surface expression and abrogated signaling of the chimeric receptors. The exchange of nine amino acids within the linker of CaR with the homologous sequence of mGluR1, however, preserved receptor function. Ala substitution for the four highly conserved residues within this amino acid sequence identified a Leu at position 606 of the CaR critical for cell surface expression and signaling. Substitution of Leu(606) for Ala resulted in impaired cell surface expression. However, Ile and Val substitutions displayed strong activating phenotypes. Disruption of the linker by insertion of nine amino acids of a random-coiled structure uncoupled the ECD from regulating the 7TMD. These data are consistent with a model of receptor activation in which the peptide linker, and particularly Leu(606), provides a critical interaction for the CaR signal transmission, a finding likely to be relevant for all family 3 GPCRs containing this conserved motif.     

Common structural requirements for heptahelical domain function in class A and class C G protein-coupled receptors

G protein-coupled receptors (GPCRs) are key players in cell communication. Several classes of such receptors have been identified. Although all GPCRs possess a heptahelical domain directly activating G proteins, important structural and sequence differences within receptors from different classes suggested distinct activation mechanisms. Here we show that highly conserved charged residues likely involved in an interaction network between transmembrane domains (TM) 3 and 6 at the cytoplasmic side of class C GPCRs are critical for activation of the gamma-aminobutyric acid type B receptor. Indeed, the loss of function resulting from the mutation of the conserved lysine residue into aspartate or glutamate in the TM3 of gamma-aminobutyric acid type B(2) can be partly rescued by mutating the conserved acidic residue of TM6 into either lysine or arginine. In addition, mutation of the conserved lysine into an acidic residue leads to a nonfunctional receptor that displays a high agonist affinity. This is reminiscent of a similar ionic network that constitutes a lock stabilizing the inactive state of many class A rhodopsin-like GPCRs. These data reveal that despite their original structure, class C GPCRs share with class A receptors at least some common structural feature controlling G protein activation.     

Structural basis of M3 muscarinic receptor dimer/oligomer formation

Class A G protein-coupled receptors (GPCRs) are known to form dimers and/or oligomeric arrays in vitro and in vivo. These complexes are thought to play important roles in modulating class A GPCR function. Many studies suggest that residues located on the "outer" (lipid-facing) surface of the transmembrane (TM) receptor core are critically involved in the formation of class A receptor dimers (oligomers). However, no clear consensus has emerged regarding the identity of the TM helices or TM subsegments involved in this process. To shed light on this issue, we have used the M(3) muscarinic acetylcholine receptor (M3R), a prototypic class A GPCR, as a model system. Using a comprehensive and unbiased approach, we subjected all outward-facing residues (70 amino acids total) of the TM helical bundle (TM1-7) of the M3R to systematic alanine substitution mutagenesis. We then characterized the resulting mutant receptors in radioligand binding and functional studies and determined their ability to form dimers (oligomers) in bioluminescence resonance energy transfer saturation assays. We found that M3R/M3R interactions are not dependent on the presence of one specific structural motif but involve the outer surfaces of multiple TM subsegments (TM1-5 and -7) located within the central and endofacial portions of the TM receptor core. Moreover, we demonstrated that the outward-facing surfaces of most TM helices play critical roles in proper receptor folding and/or function. Guided by the bioluminescence resonance energy transfer data, molecular modeling studies suggested the existence of multiple dimeric/oligomeric M3R arrangements, which may exist in a dynamic equilibrium. Given the high structural homology found among all class A GPCRs, our results should be of considerable general relevance.     

Multidimensional scaling reveals the main evolutionary pathways of class A G-protein-coupled receptors

Class A G-protein-coupled receptors (GPCRs) constitute the largest family of transmembrane receptors in the human genome. Understanding the mechanisms which drove the evolution of such a large family would help understand the specificity of each GPCR sub-family with applications to drug design. To gain evolutionary information on class A GPCRs, we explored their sequence space by metric multidimensional scaling analysis (MDS). Three-dimensional mapping of human sequences shows a non-uniform distribution of GPCRs, organized in clusters that lay along four privileged directions. To interpret these directions, we projected supplementary sequences from different species onto the human space used as a reference. With this technique, we can easily monitor the evolutionary drift of several GPCR sub-families from cnidarians to humans. Results support a model of radiative evolution of class A GPCRs from a central node formed by peptide receptors. The privileged directions obtained from the MDS analysis are interpretable in terms of three main evolutionary pathways related to specific sequence determinants. The first pathway was initiated by a deletion in transmembrane helix 2 (TM2) and led to three sub-families by divergent evolution. The second pathway corresponds to the differentiation of the amine receptors. The third pathway corresponds to parallel evolution of several sub-families in relation with a covarion process involving proline residues in TM2 and TM5. As exemplified with GPCRs, the MDS projection technique is an important tool to compare orthologous sequence sets and to help decipher the mutational events that drove the evolution of protein families.     

Molecular dynamics simulations reveal insights into key structural elements of adenosine receptors

The crystal structure of the human A(2A) adenosine receptor, a member of the G protein-coupled receptor (GPCR) family, is used as a starting point for the structural characterization of the conformational equilibrium around the inactive conformation of the human A(2) (A(2A) and A(2B)) adenosine receptors (ARs). A homology model of the closely related A(2B)AR is reported, and the two receptors were simulated in their apo form through all-atom molecular dynamics (MD) simulations. Different conditions were additionally explored in the A(2A)AR, including the protonation state of crucial histidines or the presence of the cocrystallized ligand. Our simulations reveal the role of several conserved residues in the ARs in the conformational equilibrium of the receptors. The "ionic lock" absent in the crystal structure of the inactive A(2A)AR is rapidly formed in the two simulated receptors, and a complex network of interacting residues is presented that further stabilizes this structural element. Notably, the observed rotameric transition of Trp6.48 ("toggle switch"), which is thought to initiate the activation process in GPCRs, is accompanied by a concerted rotation of the conserved residue of the A(2)ARs, His6.52. This new conformation is further stabilized in the two receptors under study by a novel interaction network involving residues in transmembrane (TM) helices TM5 (Asn5.42) and TM3 (Gln3.37), which resemble the conformational changes recently observed in the agonist-bound structure of β-adrenoreceptors. Finally, the interaction between Glu1.39 and His7.43, a pair of conserved residues in the family of ARs, is found to be weaker than previously thought, and the role of this interaction in the structure and dynamics of the receptor is thoroughly examined. All these findings suggest that, despite the commonalities with other GPCRs, the conformational equilibrium of ARs is also modulated by specific residues of the family.     

GPCRs heterodimerization: a new way towards the discovery of function for the orphan receptors?

G protein-coupled receptors (GPCRs), also called seven transmembrane domain (7TM) proteins, represent the largest family of cell surface receptors. GPCRs control a variety of physiological processes, are involved in multiple diseases and are major drug targets. Despite a vast effort of academic and industrial research, more than one hundred receptors remain orphans. These orphan GPCRs offer a great potential for drug discovery, as almost 60% of currently prescribed drugs target GPCRs. Deorphenization strategies have concentrated mainly on the identification of the natural ligands of these proteins. Recent advances have shown that orphan GPCRs, similar to orphan nuclear receptors, can regulate the function of non-orphan receptors by heterodimerization. These findings not only help to better understand the extraordinary diversity of GPCRs, but also open new perspectives for the identification of the function of these orphan receptors that hold great therapeutic potential.     

The Secretin GPCRs descended from the family of Adhesion GPCRs

The Adhesion G-protein-coupled receptors (GPCRs) are the most complex gene family among GPCRs with large genomic size, multiple introns, and a fascinating flora of functional domains, though the evolutionary origin of this family has been obscure. Here we studied the evolution of all class B (7tm2)-related genes, including the Adhesion, Secretin, and Methuselah families of GPCRs with a focus on nine genomes. We found that the cnidarian genome of Nematostella vectensis has a remarkably rich set of Adhesion GPCRs with a broad repertoire of N-terminal domains although this genome did not have any Secretin GPCRs. Moreover, the single-celled and colony-forming eukaryotes Monosiga brevicollis and Dictyostelium discoideum contain Adhesion-like GPCRs although these genomes do not have any Secretin GPCRs suggesting that the Adhesion types of GPCRs are the most ancient among class B GPCRs. Phylogenetic analysis found Adhesion group V (that contains GPR133 and GPR144) to be the closest relative to the Secretin family in the Adhesion family. Moreover, Adhesion group V sequences in N. vectensis share the same splice site setup as the Secretin GPCRs. Additionally, one of the most conserved motifs in the entire Secretin family is only found in group V of the Adhesion family. We suggest therefore that the Secretin family of GPCRs could have descended from group V Adhesion GPCRs. We found a set of unique Adhesion-like GPCRs in N. vectensis that have long N-termini containing one Somatomedin B domain each, which is a domain configuration similar to that of a set of Adhesion-like GPCRs found in Branchiostoma floridae. These sequences show slight similarities to Methuselah sequences found in insects. The extended class B GPCRs have a very complex evolutionary history with several species-specific expansions, and we identified at least 31 unique N-terminal domains originating from other protein classes. The overall N-terminal domain structure, however, concurs with the phylogenetic analysis of the transmembrane domains, thus enabling us to track the origin of most of the subgroups.     

Structure-based identification of binding sites,native ligands and potential inhibitors for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.     

A novel gene family in Arabidopsis encoding putative heptahelical transmembrane proteins homologous to human adiponectin receptors and progestin receptors

A novel seven-transmembrane receptor family, that is comprised of human adiponectin receptors (AdipoRs) and membrane progestin receptors (mPRs) that share little sequence homology with all known G protein-coupled receptors (GPCRs), has been identified recently. Although a fish mPR has been suggested to be a GPCR, human AdipoRs seem to be structurally and functionally distinct from all known GPCRs. The identification of a novel gene family, the heptahelical protein (HHP) gene family, encoding proteins in Arabidopsis predicted to have a heptahelical transmembrane topology is reported here. There are at least five HHP genes in Arabidopsis whose encoded amino acid sequences have significant similarities to human AdipoRs and mPRs.The expression and regulation of the Arabidopsis HHP gene family has been studied here. The expression of the HHP gene family is differentially regulated by plant hormones. Steady-state levels of HHP1 mRNA are increased by treatments with abscisic acid and gibberellic acid, whereas levels of HHP2 mRNA are increased by abscisic acid and benzyladenine treatments. In addition, the expression of the HHP gene family is up-regulated by the presence of sucrose in the medium. Temperature and salt stress treatments also differentially affect the expression of the HHP genes. These novel seven-transmembrane proteins previously described in yeast and animals, and now identified in plants, may represent a new class of receptors that are highly conserved across kingdoms.     

Membrane topology of a metabotropic glutamate receptor

The metabotropic glutamate receptors (mGluRs) have been predicted to have a classical seven transmembrane domain structure similar to that seen for members of the G-protein-coupled receptor (GPCR) superfamily. However, the mGluRs (and other members of the family C GPCRs) show no sequence homology to the rhodopsin-like GPCRs, for which this seven transmembrane domain structure has been experimentally confirmed. Furthermore, several transmembrane domain prediction algorithms suggest that the mGluRs have a topology that is distinct from these receptors. In the present study, we set out to test whether mGluR5 has seven true transmembrane domains. Using a variety of approaches in both prokaryotic and eukaryotic systems, our data provide strong support for the proposed seven transmembrane domain model of mGluR5. We propose that this membrane topology can be extended to all members of the family C GPCRs.     

Membrane Environment Imposes Unique Selection Pressures on Transmembrane Domains of G Protein-Coupled Receptors

We have investigated the influence of the plasma membrane environment on the molecular evolution of G protein-coupled receptors (GPCRs), the largest receptor family in Metazoa. In particular, we have analyzed the site-specific rate variation across the two primary structural partitions, transmembrane (TM) and extramembrane (EM), of these membrane proteins. We find that TM domains evolve more slowly than do EM domains, though TM domains display increased rate heterogeneity relative to their EM counterparts. Although the majority of residues across GPCRs experience strong to weak purifying selection, many GPCRs experience positive selection at both TM and EM residues, albeit with a slight bias towards the EM. Further, a subset of GPCRs, chemosensory receptors (including olfactory and taste receptors), exhibit increased rates of evolution relative to other GPCRs, an effect which is more pronounced in their TM spans. Although it has been previously suggested that the TM’s low evolutionary rate is caused by their high percentage of buried residues, we show that their attenuated rate seems to stem from the strong biophysical constraints of the membrane itself, or by functional requirements. In spite of the strong evolutionary constraints acting on the TM spans of GPCRs, positive selection and high levels of evolutionary rate variability are common. Thus, biophysical constraints should not be presumed to preclude a protein’s ability to evolve.     

Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β₂-adrenergic receptor (β₂AR) and the cholecystokinin (CCK) receptor subfamily. The β₂AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β₂AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK₁R shows functional sensitivity to membrane cholesterol while CCK₂R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β₂AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-T_m) shows that mutation of a conserved CRAC sequence on TM7 of the β₂AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.     

Molecular evolution of the transmembrane domains of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are a superfamily of integral membrane proteins vital for signaling and are important targets for pharmaceutical intervention in humans. Previously, we identified a group of ten amino acid positions (called key positions), within the seven transmembrane domain (7TM) interhelical region, which had high mutual information with each other and many other positions in the 7TM. Here, we estimated the evolutionary selection pressure at those key positions. We found that the key positions of receptors for small molecule natural ligands were under strong negative selection. Receptors naturally activated by lipids had weaker negative selection in general when compared to small molecule-activated receptors. Selection pressure varied widely in peptide-activated receptors. We used this observation to predict that a subgroup of orphan GPCRs not under strong selection may not possess a natural small-molecule ligand. In the subgroup of MRGX1-type GPCRs, we identified a key position, along with two non-key positions, under statistically significant positive selection.     

Comparative Sequence and Structural Analyses of G-Protein-Coupled Receptor Crystal Structures and Implications for Molecular Models

Background

Up until recently the only available experimental (high resolution) structure of a G-protein-coupled receptor (GPCR) was that of bovine rhodopsin. In the past few years the determination of GPCR structures has accelerated with three new receptors, as well as squid rhodopsin, being successfully crystallized. All share a common molecular architecture of seven transmembrane helices and can therefore serve as templates for building molecular models of homologous GPCRs. However, despite the common general architecture of these structures key differences do exist between them. The choice of which experimental GPCR structure(s) to use for building a comparative model of a particular GPCR is unclear and without detailed structural and sequence analyses, could be arbitrary. The aim of this study is therefore to perform a systematic and detailed analysis of sequence-structure relationships of known GPCR structures.

Methodology

We analyzed in detail conserved and unique sequence motifs and structural features in experimentally-determined GPCR structures. Deeper insight into specific and important structural features of GPCRs as well as valuable information for template selection has been gained. Using key features a workflow has been formulated for identifying the most appropriate template(s) for building homology models of GPCRs of unknown structure. This workflow was applied to a set of 14 human family A GPCRs suggesting for each the most appropriate template(s) for building a comparative molecular model.

Conclusions

The available crystal structures represent only a subset of all possible structural variation in family A GPCRs. Some GPCRs have structural features that are distributed over different crystal structures or which are not present in the templates suggesting that homology models should be built using multiple templates. This study provides a systematic analysis of GPCR crystal structures and a consistent method for identifying suitable templates for GPCR homology modelling that will help to produce more reliable three-dimensional models.     

Homology modeling and molecular dynamics simulations of the active state of the nociceptin receptor reveal new insights into agonist binding and activation

The opioid receptor-like receptor, also known as the nociceptin receptor (NOP), is a class A G protein-coupled receptor (GPCR) in the opioid receptor family. Although NOP shares a significant homology with the other opioid receptors, it does not bind known opioid ligands and has been shown to have a distinct mechanism of activation compared to the closely related opioid receptors mu, delta, and kappa. Previously reported homology models of the NOP receptor, based on the inactive-state GPCR crystal structures, give limited information on the activation and selectivity features of this fourth member of the opioid receptor family. We report here the first active-state homology model of the NOP receptor based on the opsin GPCR crystal structure. An inactive-state homology model of NOP was also built using a multiple template approach. Molecular dynamics simulation of the active-state NOP model and comparison to the inactive-state model suggest that NOP activation involves movements of transmembrane (TM)3 and TM6 and several activation microswitches, consistent with GPCR activation. Docking of the selective nonpeptidic NOP agonist ligand Ro 64-6198 into the active-state model reveals active-site residues in NOP that play a role in the high selectivity of this ligand for NOP over the other opioid receptors. Docking the shortest active fragment of endogenous agonist nociceptin/orphaninFQ (residues 1-13) shows that the NOP extracellular loop 2 (EL2) loop interacts with the positively charged residues (8-13) of N/OFQ. Both agonists show extensive polar interactions with residues at the extracellular end of the TM domain and EL2 loop, suggesting agonist-induced reorganization of polar networks, during receptor activation.