Gene expression meta-analysis in patients with schizophrenia reveals up-regulation of RGS2 and RGS16 in Brodmann Area 10.

Regulator of G-protein signalling (RGS) proteins inhibit signalling by G-protein-coupled receptors (GPCRs). GPCRs mediate the functions of several important neurotransmitters and serve as targets of many anti-psychotics. RGS2, RGS4, RGS5 and RGS16 are located on chromosome 1q23.3-31, a locus found to be associated with schizophrenia. Although previous gene expression analysis detected down-regulation of RGS4 expression in brain samples of patients with schizophrenia, the results were not consistent. In the present study, we performed a systematic meta-analysis of differential RGS2, RGS4, RGS5 and RGS16 expression in Brodmann Area 10 (BA10) samples of patients with schizophrenia and from healthy controls. Two microarray datasets met the inclusion criteria (overall, 41 schizophrenia samples and 38 controls were analysed). RGS2 and RGS16 were found to be up-regulated in BA10 samples of individuals with schizophrenia, whereas no differential expression of RGS4 and RGS5 was detected. Analysis of dorso-lateral prefrontal cortex samples of the CommonMind Consortium (258 schizophrenia samples vs. 279 controls) further validated the results. Given their central role in inactivating G-protein-coupled signalling pathways, our results suggest that differential gene expression might lead to enhanced inactivation of G-protein signalling in schizophrenia. This, in turn, suggests that additional studies are needed to further explore the consequences of the differential expression we detected, this time at the protein and functional levels.

Residue-level determinants of RGS R4 subfamily GAP activity and specificity towards the Gi subfamily.

The structural basis for the GTPase-accelerating activity of regulators of G protein signaling (RGS) proteins, as well as the mechanistic basis for their specificity in interacting with the heterotrimeric (αßγ) G proteins they inactivate, is not sufficiently understood at the family level. Here, we used biochemical assays to compare RGS domains across the RGS family and map those individual residues that favorably contribute to GTPase-accelerating activity, and those residues responsible for attenuating RGS domain interactions with Gα subunits. We show that conserved interactions of RGS residues with both the Gα switch I and II regions are crucial for RGS activity, while the reciprocal effects of "modulatory" and "disruptor" residues selectively modulate RGS activity. Our results quantify how specific interactions between RGS domains and Gα subunits are set by a balance between favorable RGS residue interactions with particular Gα switch regions, and unfavorable interactions with the Gα helical domain.

Residue-level determinants of angiopoietin-2 interactions with its receptor Tie2.

We combined computational and experimental methods to interrogate the binding determinants of angiopoietin-2 (Ang2) to its receptor tyrosine kinase (RTK) Tie2-a central signaling system in angiogenesis, inflammation, and tumorigenesis. We used physics-based electrostatic and surface-area calculations to identify the subset of interfacial Ang2 and Tie2 residues that can affect binding directly. Using random and site-directed mutagenesis and yeast surface display (YSD), we validated these predictions and identified additional Ang2 positions that affected receptor binding. We then used burial-based calculations to classify the larger set of Ang2 residues that are buried in the Ang2 core, whose mutations can perturb the Ang2 structure and thereby affect interactions with Tie2 indirectly. Our analysis showed that the Ang2-Tie2 interface is dominated by nonpolar contributions, with only three Ang2 and two Tie2 residues that contribute electrostatically to intermolecular interactions. Individual interfacial residues contributed only moderately to binding, suggesting that engineering of this interface will require multiple mutations to reach major effects. Conversely, substitutions in substantially buried Ang2 residues were more prevalent in our experimental screen, reduced binding substantially, and are therefore more likely to have a deleterious effect that might contribute to oncogenesis. Computational analysis of additional RTK-ligand complexes, c-Kit-SCF and M-CSF-c-FMS, and comparison to previous YSD results, further show the utility of our combined methodology.

Interplay between negative and positive design elements in Gα helical domains of G proteins determines interaction specificity toward RGS2.

Regulators of G protein signaling (RGS) proteins inactivate Gα subunits, thereby controlling G protein-coupled signaling networks. Among all RGS proteins, RGS2 is unique in interacting only with the Gαq but not with the Gαi subfamily. Previous studies suggested that this specificity is determined by the RGS domain and, in particular, by three RGS2-specific residues that lead to a unique mode of interaction with Gαq This interaction was further proposed to act through contacts with the Gα GTPase domain. Here, we combined energy calculations and GTPase activity measurements to determine which Gα residues dictate specificity toward RGS2. We identified putative specificity-determining residues in the Gα helical domain, which among G proteins is found only in Gα subunits. Replacing these helical domain residues in Gαi with their Gαq counterparts resulted in a dramatic specificity switch toward RGS2. We further show that Gα-RGS2 specificity is set by Gαi residues that perturb interactions with RGS2, and by Gαq residues that enhance these interactions. These results show, for the first time, that the Gα helical domain is central to dictating specificity toward RGS2, suggesting that this domain plays a general role in governing Gα-RGS specificity. Our insights provide new options for manipulating RGS-G protein interactions in vivo, for better understanding of their 'wiring' into signaling networks, and for devising novel drugs targeting such interactions.

Structural motifs in the RGS RZ subfamily combine to attenuate interactions with Gα subunits.

Regulators of G-protein Signaling (RGS) proteins inactivate heterotrimeric G proteins, thereby setting the duration of active signaling. In particular, the RGS RZ subfamily, which consists of RGS17, RGS19, and RGS20, mediates numerous physiological functions and human pathologies - mostly by functioning as GTPase Activating Proteins (GAPs) towards the Gαi subfamily. Yet, which RZ subfamily members mediate particular functions and how their GAP activity and specificity are governed at the amino acid level is not well understood. Here, we show that all RZ subfamily members have similar and relatively low GAP activity towards Gαo. We characterized four RZ-specific structural motifs that mediate this low activity, and suggest they perturb optimal interactions with the Gα subunit. Indeed, inserting these RZ-specific motifs into the representative high-activity RGS16 impaired GAP activity in a non-additive manner. Our results provide residue-level insights into the specificity determinants of the RZ subfamily, and enable to study their interactions in signaling cascades by using redesigned mutants such as those presented in this work.

Engineering a monomeric variant of macrophage colony-stimulating factor (M-CSF) that antagonizes the c-FMS receptor.

Enhanced activation of the signaling pathways that mediate the differentiation of mononuclear monocytes into osteoclasts is an underlying cause of several bone diseases and bone metastasis. In particular, dysregulation and overexpression of macrophage colony-stimulating factor (M-CSF) and its c-FMS tyrosine kinase receptor, proteins that are essential for osteoclast differentiation, are known to promote bone metastasis and osteoporosis, making both the ligand and its receptor attractive targets for therapeutic intervention. With this aim in mind, our starting point was the previously held concept that the potential of the M-CSFC31S mutant as a therapeutic is derived from its inability to dimerize and hence to act as an agonist. The current study showed, however, that dimerization is not abolished in M-CSFC31S and that the protein retains agonistic activity toward osteoclasts. To design an M-CSF mutant with diminished dimerization capabilities, we solved the crystal structure of the M-CSFC31S dimer complex and used structure-based energy calculations to identify the residues responsible for its dimeric form. We then used that analysis to develop M-CSFC31S,M27R, a ligand-based, high-affinity antagonist for c-FMS that retained its binding ability but prevented the ligand dimerization that leads to receptor dimerization and activation. The monomeric properties of M-CSFC31S,M27R were validated using dynamic light scattering and small-angle X-ray scattering analyses. It was shown that this mutant is a functional inhibitor of M-CSF-dependent c-FMS activation and osteoclast differentiation in vitro Our study, therefore, provided insights into the sequence-structure-function relationships of the M-CSF/c-FMS interaction and of ligand/receptor tyrosine kinase interactions in general.

Compartmentalization and Ca2+ buffering are essential for prevention of light-induced retinal degeneration.

Fly photoreceptors are polarized cells, each of which has an extended interface between its cell body and the light-signaling compartment, the rhabdomere. Upon intense illumination, rhabdomeric calcium concentration reaches millimolar levels that would be toxic if Ca(2+) diffusion between the rhabdomere and cell body was not robustly attenuated. Yet, it is not clear how such effective attenuation is obtained. Here we show that Ca(2+) homeostasis in the photoreceptor cell relies on the protein calphotin. This unique protein functions as an immobile Ca(2+) buffer localized along the base of the rhabdomere, separating the signaling compartment from the cell body. Generation and analyses of transgenic Drosophila strains, in which calphotin-expression levels were reduced in a graded manner, showed that moderately reduced calphotin expression impaired Ca(2+) homeostasis while calphotin elimination resulted in severe light-dependent photoreceptor degeneration. Electron microscopy, electrophysiology, and optical methods revealed that the degeneration was rescued by prevention of Ca(2+) overload via overexpression of CalX, the Na(+)-Ca(2+) exchanger. In addition, Ca(2+)-imaging experiments showed that reduced calphotin levels resulted in abnormally fast kinetics of Ca(2+) elevation in photoreceptor cells. Together, the data suggest that calphotin functions as a Ca(2+) buffer; a possibility that we directly demonstrate by expressing calphotin in a heterologous expression system. We propose that calphotin-mediated compartmentalization and Ca(2+) buffering constitute an effective strategy to protect cells from Ca(2+) overload and light-induced degeneration.

Regulator of G-protein signaling-21 (RGS21) is an inhibitor of bitter gustatory signaling found in lingual and airway epithelia.

The gustatory system detects tastants and transmits signals to the brain regarding ingested substances and nutrients. Although tastant receptors and taste signaling pathways have been identified, little is known about their regulation. Because bitter, sweet, and umami taste receptors are G protein-coupled receptors (GPCRs), we hypothesized that regulators of G protein signaling (RGS) proteins may be involved. The recent cloning of RGS21 from taste bud cells has implicated this protein in the regulation of taste signaling; however, the exact role of RGS21 has not been precisely defined. Here, we sought to determine the role of RGS21 in tastant responsiveness. Biochemical analyses confirmed in silico predictions that RGS21 acts as a GTPase-accelerating protein (GAP) for multiple G protein α subunits, including adenylyl cyclase-inhibitory (Gα(i)) subunits and those thought to be involved in tastant signal transduction. Using a combination of in situ hybridization, RT-PCR, immunohistochemistry, and immunofluorescence, we demonstrate that RGS21 is not only endogenously expressed in mouse taste buds but also in lung airway epithelial cells, which have previously been shown to express components of the taste signaling cascade. Furthermore, as shown by reverse transcription-PCR, the immortalized human airway cell line 16HBE was found to express transcripts for tastant receptors, RGS21, and downstream taste signaling components. Over- and underexpression of RGS21 in 16HBE cells confirmed that RGS21 acts to oppose bitter tastant signaling to cAMP and calcium second messenger changes. Our data collectively suggests that RGS21 modulates bitter taste signal transduction.

Cooperative phosphoinositide and peptide binding by PSD-95/discs large/ZO-1 (PDZ) domain of polychaetoid, Drosophila zonulin.

PDZ domains are well known protein-protein interaction modules that, as part of multidomain proteins, assemble molecular complexes. Some PDZ domains have been reported to interact with membrane lipids, in particular phosphatidylinositol phosphates, but few studies have been aimed at elucidating the prevalence or the molecular details of such interactions. We screened 46 Drosophila PDZ domains for phosphoinositide-dependent cellular localization and discovered that the second PDZ domain of polychaetoid (Pyd PDZ2) interacts with phosphatidylinositol 4,5-bisphosphate (PtdIns(4,5)P(2)) at the plasma membrane. Surface plasmon resonance binding experiments with recombinant protein established that Pyd PDZ2 interacts with phosphatidylinositol phosphates with apparent affinities in the micromolar range. Electrostatic interactions involving an extended positively charged surface of Pyd PDZ2 are crucial for the PtdIns(4,5)P(2)-dependent membrane interactions as shown by a combination of three-dimensional modeling, mutagenesis, binding, and localization studies. In vivo localization studies further suggested that both lipid and peptide binding contribute to membrane localization. We identified the transmembrane protein Crumbs as a Pyd PDZ2 ligand and probed the relation between peptide and PtdIns(4,5)P(2) binding. Contrary to the prevalent view on PDZ/peptide/lipid binding, we did not find competition between peptide and lipid ligands. Instead, preloading the protein with the 10-mer Crb3 peptide increased the apparent affinity of Pyd PDZ2 for PtdIns(4,5)P(2) 6-fold. Our results suggest that membrane localization of Pyd PDZ2 may be driven by a combination of peptide and PtdIns(4,5)P(2) binding, which raises the intriguing possibility that the domain may coordinate protein- and phospholipid-mediated signals.

Evolution of homo-oligomerization of methionine S-adenosyltransferases is replete with structure-function constrains.

Homomers are prevalent in bacterial proteomes, particularly among core metabolic enzymes. Homomerization is often key to function and regulation, and interfaces that facilitate the formation of homomeric enzymes are subject to intense evolutionary change. However, our understanding of the molecular mechanisms that drive evolutionary variation in homomeric complexes is still lacking. How is the diversification of protein interfaces linked to variation in functional regulation and structural integrity of homomeric complexes? To address this question, we studied quaternary structure evolution of bacterial methionine S-adenosyltransferases (MATs)-dihedral homotetramers formed along a large and conserved dimeric interface harboring two active sites, and a small, recently evolved, interdimeric interface. Here, we show that diversity in the physicochemical properties of small interfaces is directly linked to variability in the kinetic stability of MAT quaternary complexes and in modes of their functional regulation. Specifically, hydrophobic interactions within the small interface of Escherichia coli MAT render the functional homotetramer kinetically stable yet impose severe aggregation constraints on complex assembly. These constraints are alleviated by electrostatic interactions that accelerate dimer-dimer assembly. In contrast, Neisseria gonorrhoeae MAT adopts a nonfunctional dimeric state due to the low hydrophobicity of its small interface and the high flexibility of its active site loops, which perturbs small interface integrity. Remarkably, in the presence of methionine and ATP, N. gonorrhoeae MAT undergoes substrate-induced assembly into a functional tetrameric state. We suggest that evolution acts on the interdimeric interfaces of MATs to tailor the regulation of their activity and stability to unique organismal needs.

A guide to adhesion GPCR research.

Adhesion G protein-coupled receptors (aGPCRs) are a class of structurally and functionally highly intriguing cell surface receptors with essential functions in health and disease. Thus, they display a vastly unexploited pharmacological potential. Our current understanding of the physiological functions and signaling mechanisms of aGPCRs form the basis for elucidating further molecular aspects. Combining these with novel tools and methodologies from different fields tailored for studying these unusual receptors yields a powerful potential for pushing aGPCR research from singular approaches toward building up an in-depth knowledge that will facilitate its translation to applied science. In this review, we summarize the state-of-the-art knowledge on aGPCRs in respect to structure-function relations, physiology, and clinical aspects, as well as the latest advances in the field. We highlight the upcoming most pressing topics in aGPCR research and identify strategies to tackle them. Furthermore, we discuss approaches how to promote, stimulate, and translate research on aGPCRs 'from bench to bedside' in the future.

Structural design principles for specific ultra-high affinity interactions between colicins/pyocins and immunity proteins.

The interactions of the antibiotic proteins colicins/pyocins with immunity proteins is a seminal model system for studying protein-protein interactions and specificity. Yet, a precise and quantitative determination of which structural elements and residues determine their binding affinity and specificity is still lacking. Here, we used comparative structure-based energy calculations to map residues that substantially contribute to interactions across native and engineered complexes of colicins/pyocins and immunity proteins. We show that the immunity protein α1-α2 motif is a unique structurally-dissimilar element that restricts interaction specificity towards all colicins/pyocins, in both engineered and native complexes. This motif combines with a diverse and extensive array of electrostatic/polar interactions that enable the exquisite specificity that characterizes these interactions while achieving ultra-high affinity. Surprisingly, the divergence of these contributing colicin residues is reciprocal to residue conservation in immunity proteins. The structurally-dissimilar immunity protein α1-α2 motif is recognized by divergent colicins similarly, while the conserved immunity protein α3 helix interacts with diverse colicin residues. Electrostatics thus plays a key role in setting interaction specificity across all colicins and immunity proteins. Our analysis and resulting residue-level maps illuminate the molecular basis for these protein-protein interactions, with implications for drug development and rational engineering of these interfaces.

Integrative Structural Biology in the Era of Accurate Structure Prediction.

Characterizing the three-dimensional structure of macromolecules is central to understanding their function. Traditionally, structures of proteins and their complexes have been determined using experimental techniques such as X-ray crystallography, NMR, or cryo-electron microscopy-applied individually or in an integrative manner. Meanwhile, however, computational methods for protein structure prediction have been improving their accuracy, gradually, then suddenly, with the breakthrough advance by AlphaFold2, whose models of monomeric proteins are often as accurate as experimental structures. This breakthrough foreshadows a new era of computational methods that can build accurate models for most monomeric proteins. Here, we envision how such accurate modeling methods can combine with experimental structural biology techniques, enhancing integrative structural biology. We highlight the challenges that arise when considering multiple structural conformations, protein complexes, and polymorphic assemblies. These challenges will motivate further developments, both in modeling programs and in methods to solve experimental structures, towards better and quicker investigation of structure-function relationships.

Ectopic activation of GABAB receptors inhibits neurogenesis and metamorphosis in the cnidarian Nematostella vectensis.

The metabotropic gamma-aminobutyric acid B receptor (GABABR) is a G protein-coupled receptor that mediates neuronal inhibition by the neurotransmitter GABA. While GABABR-mediated signalling has been suggested to play central roles in neuronal differentiation and proliferation across evolution, it has mostly been studied in the mammalian brain. Here, we demonstrate that ectopic activation of GABABR signalling affects neurogenic functions in the sea anemone Nematostella vectensis. We identified four putative Nematostella GABABR homologues presenting conserved three-dimensional extracellular domains and residues needed for binding GABA and the GABABR agonist baclofen. Moreover, sustained activation of GABABR signalling reversibly arrests the critical metamorphosis transition from planktonic larva to sessile polyp life stage. To understand the processes that underlie the developmental arrest, we combined transcriptomic and spatial analyses of control and baclofen-treated larvae. Our findings reveal that the cnidarian neurogenic programme is arrested following the addition of baclofen to developing larvae. Specifically, neuron development and neurite extension were inhibited, resulting in an underdeveloped and less organized nervous system and downregulation of proneural factors including NvSoxB(2), NvNeuroD1 and NvElav1. Our results thus point to an evolutionarily conserved function of GABABR in neurogenesis regulation and shed light on early cnidarian development.

The European Research Network on Signal Transduction (ERNEST): Toward a Multidimensional Holistic Understanding of G Protein-Coupled Receptor Signaling.

G protein-coupled receptors (GPCRs) are intensively studied due to their therapeutic potential as drug targets. Members of this large family of transmembrane receptor proteins mediate signal transduction in diverse cell types and play key roles in human physiology and health. In 2013 the research consortium GLISTEN (COST Action CM1207) was founded with the goal of harnessing the substantial growth in knowledge of GPCR structure and dynamics to push forward the development of molecular modulators of GPCR function. The success of GLISTEN, coupled with new findings and paradigm shifts in the field, led in 2019 to the creation of a related consortium called ERNEST (COST Action CA18133). ERNEST broadens focus to entire signaling cascades, based on emerging ideas of how complexity and specificity in signal transduction are not determined by receptor-ligand interactions alone. A holistic approach that unites the diverse data and perspectives of the research community into a single multidimensional map holds great promise for improved drug design and therapeutic targeting.

Structural design principles that underlie the multi-specific interactions of Gαq with dissimilar partners.

Gαq is a ubiquitous molecular switch that activates the effectors phospholipase-C-ß3 (PLC-ß3) and Rho guanine-nucleotide exchange factors. Gαq is inactivated by regulators of G protein signaling proteins, as well as by PLC-ß3. Gαq further interacts with G protein-coupled receptor kinase 2 (GRK2), although the functional role of this interaction is debated. While X-ray structures of Gαq bound to representatives of these partners have revealed details of their interactions, the mechanistic basis for differential Gαq interactions with multiple partners (i.e., Gαq multi-specificity) has not been elucidated at the individual residue resolution. Here, we map the structural determinants of Gαq multi-specificity using structure-based energy calculations. We delineate regions that specifically interact with GTPase Activating Proteins (GAPs) and residues that exclusively contribute to effector interactions, showing that only the Gαq "Switch II" region interacts with all partners. Our analysis further suggests that Gαq-GRK2 interactions are consistent with GRK2 functioning as an effector, rather than a GAP. Our multi-specificity analysis pinpoints Gαq residues that uniquely contribute to interactions with particular partners, enabling precise manipulation of these cascades. As such, we dissect the molecular basis of Gαq function as a central signaling hub, which can be used to target Gαq-mediated signaling in therapeutic interventions.

RGS6 and RGS7 Discriminate between the Highly Similar Gαi and Gαo Proteins Using a Two-Tiered Specificity Strategy.

RGS6 and RGS7 are regulators of G protein signaling (RGS) proteins that inactivate heterotrimeric (αßγ) G proteins and mediate diverse biological functions, such as cardiac and neuronal signaling. Uniquely, both RGS6 and RGS7 can discriminate between Gαo and Gαi1-two similar Gα subunits that belong to the same Gi sub-family. Here, we show that the isolated RGS domains of RGS6 and RGS7 are sufficient to achieve this specificity. We identified three specific RGS6/7 "disruptor residues" that can attenuate RGS interactions toward Gα subunits and demonstrated that their insertion into a representative high-activity RGS causes a significant, yet non-specific, reduction in activity. We further identified a unique "modulatory" residue that bypasses this negative effect, specifically toward Gαo. Hence, the exquisite specificity of RGS6 and RGS7 toward closely related Gα subunits is achieved via a two-tier specificity system, whereby a Gα-specific modulatory motif overrides the inhibitory effect of non-specific disruptor residues. Our findings expand the understanding of the molecular toolkit used by the RGS family to achieve specific interactions with selected Gα subunits-emphasizing the functional importance of the RGS domain in determining the activity and selectivity of RGS R7 sub-family members toward particular Gα subunits.

Author Correction: RGS7 is recurrently mutated in melanoma and promotes migration and invasion of human cancer cells.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Sequence-similar, structure-dissimilar protein pairs in the PDB.

It is often assumed that in the Protein Data Bank (PDB), two proteins with similar sequences will also have similar structures. Accordingly, it has proved useful to develop subsets of the PDB from which "redundant" structures have been removed, based on a sequence-based criterion for similarity. Similarly, when predicting protein structure using homology modeling, if a template structure for modeling a target sequence is selected by sequence alone, this implicitly assumes that all sequence-similar templates are equivalent. Here, we show that this assumption is often not correct and that standard approaches to create subsets of the PDB can lead to the loss of structurally and functionally important information. We have carried out sequence-based structural superpositions and geometry-based structural alignments of a large number of protein pairs to determine the extent to which sequence similarity ensures structural similarity. We find many examples where two proteins that are similar in sequence have structures that differ significantly from one another. The source of the structural differences usually has a functional basis. The number of such proteins pairs that are identified and the magnitude of the dissimilarity depend on the approach that is used to calculate the differences; in particular sequence-based structure superpositioning will identify a larger number of structurally dissimilar pairs than geometry-based structural alignments. When two sequences can be aligned in a statistically meaningful way, sequence-based structural superpositioning provides a meaningful measure of structural differences. This approach and geometry-based structure alignments reveal somewhat different information and one or the other might be preferable in a given application. Our results suggest that in some cases, notably homology modeling, the common use of nonredundant datasets, culled from the PDB based on sequence, may mask important structural and functional information. We have established a data base of sequence-similar, structurally dissimilar protein pairs that will help address this problem (http://luna.bioc.columbia.edu/rachel/seqsimstrdiff.htm).

Neofunctionalization in Ligand Binding Sites of Ant Olfactory Receptors.

Chemical communication is fundamental for the operation of insect societies. Their diverse vocabulary of chemical signals requires a correspondingly diverse set of chemosensory receptors. Insect olfactory receptors (ORs) are the largest family of chemosensory receptors. The OR family is characterized by frequent expansions of subfamilies, in which duplicated ORs may adapt to detect new signals through positive selection on their amino acid sequence. Ants are an extreme example with ∼400 ORs per genome-the highest number in insects. Presumably, this reflects an increased complexity of chemical communication. Here, we examined gene duplications and positive selection on ant ORs. We reconstructed the hymenopteran OR gene tree, including five ant species, and inferred positive selection along every branch using the branch-site test, a total of 3326 tests. We find more positive selection in branches following species-specific duplications. We identified amino acid sites targeted by positive selection, and mapped them onto a structural model of insect ORs. Seventeen sites were under positive selection in six or more branches, forming two clusters on the extracellular side of the receptor, on either side of a cleft in the structure. This region was previously implicated in ligand activation, suggesting that the concentration of positively selected sites in this region is related to adaptive evolution of ligand binding sites or allosteric transmission of ligand activation. These results provide insights into the specific OR subfamilies and individual residues that facilitated adaptive evolution of olfactory functions, potentially explaining the elaboration of chemical signaling in ant societies.