Non-classical amine recognition evolved in a large clade of olfactory receptors
Abstract
Biogenic amines are important signaling molecules, and the structural basis for their recognition by G Protein-Coupled Receptors (GPCRs) is well understood. Amines are also potent odors, with some activating olfactory trace amine-associated receptors (TAARs). Here, we report that teleost TAARs evolved a new way to recognize amines in a non-classical orientation. Chemical screens de-orphaned eleven zebrafish TAARs, with agonists including serotonin, histamine, tryptamine, 2-phenylethylamine, putrescine, and agmatine. Receptors from different clades contact ligands through aspartates on transmembrane α-helices III (canonical Asp3.32) or V (non-canonical Asp5.42), and diamine receptors contain both aspartates. Non-classical monoamine recognition evolved in two steps: an ancestral TAAR acquired Asp5.42, gaining diamine sensitivity, and subsequently lost Asp3.32. Through this transformation, the fish olfactory system dramatically expanded its capacity to detect amines, ecologically significant aquatic odors. The evolution of a second, alternative solution for amine detection by olfactory receptors highlights the tremendous structural versatility intrinsic to GPCRs.
https://doi.org/10.7554/eLife.10441.001eLife digest
Many organisms make molecules called biogenic amines. These molecules, which include the human hormones adrenaline and histamine, have important roles in regulating the biology and behaviour of many animals. Some biogenic amines bind to receptor proteins called GPCRs on the surface of cells. Many drugs can affect the activity of GPCRs, so understanding how different GPCRs work is an important goal of the pharmaceutical industry. Like all proteins, GPCRs are made of chains of molecules called amino acids. The GPCRs that can detect biogenic amines use a particular amino acid named Asp3.32, and when this amino acid is mutated, these GPCRs become unable to bind to their target amine.
Trace amine-associated receptors (TAARs) are a type of GPCR that are found in many animals to detect odors. Most TAARs in mammals contain the Asp3.32 residue, and recognize amine odors. However, many fish TAARs do not contain Asp3.32, and it was not clear what molecules these fish receptors detect.
Here Li et al. find that these fish TAARs also recognize amines, and use a different amino acid called Asp5.42. Also, some TAARs contain both Asp3.32 and Asp5.42, and recognize chemicals with two amines named diamines. Some diamines that bind to TAARs are foul smelling odors; for example, cadaverine and putrescine are repulsive smells emitted by decomposing flesh. In total, the experiments identified amines that can bind to eleven zebrafish TAARs that previously had no odor partner.
Li et al. propose that some fish TAARs lost the Asp3.32 during the course of evolution to leave the Asp5.42 as the main interaction site for amines. This change dramatically altered how these TAARs interact with amines, which probably expanded the number of different amines that fish can detect. These findings open up new ways to study how the fish brain processes information about its surroundings.
https://doi.org/10.7554/eLife.10441.002Introduction
Biogenic amines are key intercellular signaling molecules that regulate physiology and behavior. In vertebrates, biogenic amines include adrenaline, serotonin, acetylcholine, histamine, and dopamine, and these biogenic amines can be detected by G Protein-Coupled Receptors (GPCRs) expressed on the plasma membrane of target cells. The biochemical and structural basis for biogenic amine recognition by GPCRs is well understood and involves an essential salt bridge between the ligand amino group and a highly conserved aspartic acid on the third transmembrane α-helix of the receptor (Asp3.32; Ballesteros-Weinstein indexing) (Katritch et al., 2013; Shi and Javitch, 2002). Mutation of Asp3.32 in adrenaline, serotonin, acetylcholine, histamine, and dopamine receptors eliminates amine recognition (Shi and Javitch, 2002; Surgand et al., 2006), and GPCRs that recognize amines through another motif have not been identified.
Trace amine-associated receptors (TAARs) are a family of vertebrate GPCRs distantly related to biogenic amine receptors (Liberles, 2015). There are 15 TAARs in mouse, 6 in human, and 112 in zebrafish, and in these or related species, all TAARs except TAAR1 function as olfactory receptors (Horowitz et al., 2014; Hussain et al., 2009; Liberles and Buck, 2006). All human TAARs (6/6), most mouse TAARs (13/15), and some zebrafish TAARs (27/112) retain Asp3.32, suggesting that many TAARs would recognize amines. High throughput chemical screens yielded ligands for TAAR1 (Borowsky et al., 2001; Bunzow et al., 2001; Revel et al., 2011) and many olfactory TAARs that retain Asp3.32, including TAAR3, TAAR4, TAAR5, TAAR7s, TAAR8s, TAAR9, and TAAR13c (Liberles and Buck, 2006; Ferrero et al., 2012; Hussain et al., 2013). Each of these receptors detects different amines, some of which evoke innate olfactory behaviors (Hussain et al., 2013; Dewan et al., 2013; Ferrero et al., 2011; Li et al., 2013; Li and Liberles, 2015). In mouse, TAAR4 detects the aversive carnivore odor 2-phenylethylamine (Ferrero et al., 2011), while TAAR5 detects the attractive and sexually dimorphic mouse odor trimethylamine (Liberles and Buck, 2006; Li et al., 2013). Knockout mice lacking TAAR4 and TAAR5 lose behavioral responses to TAAR ligands identified in cell culture studies (Dewan et al., 2013; Li et al., 2013).
The selectivity of olfactory TAARs for particular amine odors suggests structural variability within the ligand-binding pocket across the TAAR family. One clade of TAARs (including TAAR1, TAAR3, and TAAR4) detects primary amines derived by amino acid decarboxylation, while a second clade (including TAAR5, TAAR7s, TAAR8s, and TAAR9) prefers tertiary amines (Ferrero et al., 2012). Structural models and mutagenesis studies of TAAR1, TAAR7e, and TAAR7f predicted that ligand binding occurs in the membrane plane, Asp3.32 forms a salt bridge to the ligand amino group, and other transmembrane residues provide a scaffold that supports ligand binding and enables selectivity (Ferrero et al., 2012; Reese et al., 2014).
A third TAAR clade (clade III TAARs) arose in fish, and includes the majority (87/112) of zebrafish TAARs (Hussain et al., 2009). Most clade III TAARs lack Asp3.32 (85/87), and it has been proposed that these receptors may not detect amines (Hussain et al., 2009; Liberles, 2014). However, ligands have not been identified for any of these receptors. More generally, ligands are known for only one zebrafish TAAR, TAAR13c, which retains Asp3.32 and detects the carrion odor cadaverine (Hussain et al., 2013). Structure-activity analysis revealed that TAAR13c preferentially recognizes odd-chained, medium-sized diamines, and likely contains two distinct cation recognition sites within the ligand-binding pocket (Hussain et al., 2013). However, the structural basis for diamine recognition by TAAR13c was unknown.
Here, we find that TAAR13c recognizes cadaverine through two remote transmembrane aspartic acids, Asp3.32 and Asp5.42. A TAAR13c variant with a charge-neutralizing mutation at Asp5.42 has an altered ligand preference for monoamines rather than diamines. We propose that Asp3.32 and Asp5.42 together form a general di-cation recognition motif, and note this motif to be present in several zebrafish, mouse, and human TAARs. Surprisingly, Asp5.42 is retained in 84 out of 85 zebrafish TAARs that lack Asp3.32, raising the possibility that these receptors recognize amines positioned in a non-canonical inverted orientation with the amino group pointing towards transmembrane α-helix V. Using a high throughput chemical screening approach, we identified ligands for eleven additional zebrafish TAARs, and found that different receptors detect amines through salt bridge contacts at Asp3.32 and/or Asp5.42. TAAR10a, TAAR10b, TAAR12h, and TAAR12i are Asp3.32-containing clade I TAARs that detect serotonin, tryptamine, 2-phenylethylamine, and 3-methoxytyramine, while TAAR16c, TAAR16e, and TAAR16f are Asp5.42-containing clade III TAARs that detect N-methylpiperidine, N,N-dimethylcyclohexylamine, and isoamylamine. TAAR13a, TAAR13d, TAAR13e, and TAAR14d have both Asp3.32 and Asp5.42, and recognize the di-cationic molecules putrescine, histamine, and agmatine.
Based on these findings, we propose that a large clade of olfactory GPCRs evolved a novel way to detect amines that involves a non-classical salt bridge to transmembrane α-helix V. Furthermore, new TAAR ligands include several biogenic amines for which other receptors were identified, including serotonin, histamine, and 2-phenylethylamine. These findings illustrate that different members of the GPCR superfamily can evolve divergent ligand binding pockets for recognition of the same agonist.
Results
Two remote amine recognition sites in the cadaverine receptor
To gain structural insights into diamine recognition by the cadaverine receptor TAAR13c, we generated a TAAR13c homology model (Figure 1A) based on the X-ray crystal structure of agonist-bound β2 adrenergic receptor (Protein Data Bank Entry 4LDE) (Ring et al., 2013). The β2 adrenergic receptor is the closest homolog (33% amino acid identity) of TAAR13c for which an active-state structure is currently available. The homology model of TAAR13c indicated a canonical fold of seven transmembrane α-helices with ligand bound at the core. Using the homology model, we performed computational docking of cadaverine to the orthosteric binding site. One cadaverine amino group of the docked ligand was located near Asp1123.32, which corresponds to the highly conserved amine-contact site in biogenic amine receptors. The other cadaverine amino group was positioned at the opposing end of the ligand-binding pocket, 3 Å away from another aspartate, Asp2025.42. The second amino group was also within 5 Å of side chains from Thr2035.43 and Ser2766.55, Leu192 and Phe194 on extracellular loop 2, and the backbone of Ser1995.39 (additional nearby residues are depicted in Figure 1—figure supplement 1). Based on its proximity and charge, Asp2025.42 was deemed a prime candidate to function as a counterion that forms a salt bridge to the second cadaverine amino group. Position 5.42 in several other GPCRs directly contacts ligands through hydrogen bonds, salt bridges, or van der Waals interactions (Surgand et al., 2006), and our homology model suggests a role for this amino acid in ligand recognition by TAAR13c.
We used an established reporter gene assay (Liberles and Buck, 2006; Ferrero et al., 2012) to query the response properties of TAAR13c variants with charge-neutralizing mutations at Asp3.32 (TAAR13cD112A) and Asp5.42 (TAAR13cD202A). HEK-293 cells were co-transfected with plasmids encoding TAARs and a cAMP-dependent reporter gene encoding secreted alkaline phosphatase (Cre-Seap). Transfected cells were exposed to test ligands, and assayed for reporter activity using a fluorescent phosphatase substrate. As previously published with this paradigm (Hussain et al., 2013), cadaverine robustly activated TAAR13c, with a half maximal response occurring at ~10 μM. In contrast, cadaverine failed to activate either the TAAR13cD112A or TAAR13cD202A mutant at any concentration tested, up to 1 mM (Figure 1, Figure 1—figure supplement 1).
To ensure that mutation of Asp2025.42 did not impair protein folding or G protein coupling, we asked whether the mutant receptor could be activated by other ligands. We tested 50 chemicals for the ability to activate TAAR13cD202A, and identified the monoamine 3-methoxytyramine as an agonist (Figure 1D). 3-methoxytyramine failed to activate wild type TAAR13c at any concentration tested. Thus, charge-neutralizing mutation of Asp5.42 transformed TAAR13c from a diamine receptor into a monoamine receptor. The altered ligand-binding preference of the TAAR13cD202A mutant suggests that Asp5.42 contributes directly to the ligand-binding pocket, and homology modeling of TAAR13cD202A bound to 3-methoxytyramine (Figure 1—figure supplement 1) supports a direct interaction between this site and ligand. Based on modeling and mutagenesis studies, we conclude that Asp5.42 in TAAR13c forms a salt bridge to a cadaverine amino group.
Asp5.42 is widely conserved in clade III TAARs
We next asked whether Asp5.42 was present in other GPCRs, and might function more generally in amine detection. Biogenic amine receptors use Asp3.32 to contact primary amines (Shi and Javitch, 2002; Surgand et al., 2006), and Asp5.42 is not present in any of 12 serotonin, 5 dopamine, 5 acetylcholine, or 9 adrenergic receptors in mouse (Figure 2). Asp3.32 is also observed in all 4 mouse histamine receptors, while Asp5.42 is observed in one (histamine receptor H2). Histamine is di-cationic at low pH, and the polar imidazole group of the histamine ligand is predicted to contact a transmembrane α-helix 5 asparagine in histamine receptor H1 (Asn5.46) and transmembrane α-helix 5 anions in histamine receptors H2 (Asp5.42), H3 (Glu5.46), and H4 (Glu5.46) (Seifert et al., 2013).
Most mammalian TAARs retain Asp3.32, including 6/6 human TAARs, 13/15 mouse TAARs, and 15/17 rat TAARs. In mouse and rat, two TAAR7s have conservative changes of Asp3.32 to Glu3.32; the effect of this substitution is unknown as ligands have not been found for Glu3.32-containing TAARs. A small number of mammalian TAARs have an aspartate at position 5.43 including two in humans (TAAR6, TAAR8), mice (TAAR6, TAAR8b), and rats (TAAR6, TAAR8a). Sequence alignments indicate that Asp5.43 of these receptors corresponds to Asp5.42 of TAAR13c, and Asp5.43 may be similarly positioned in the agonist-binding pocket. These mammalian TAARs are candidates to function as diamine receptors; however, ligands have not been identified for these receptors, perhaps due to difficulties in achieving their functional expression in heterologous cells. Mice lack a TAAR13c ortholog, yet a role for other TAARs in diamine recognition is supported by the finding that mice lacking all olfactory TAARs fail to avoid cadaverine odor (Dewan et al., 2013).
Analysis of the zebrafish TAAR repertoire surprisingly revealed that in contrast to mammalian TAARs and biogenic amine receptors, the vast majority of zebrafish TAARs (91/112) have an aspartate at position 5.42. Some of these receptors (7/91) also retain Asp3.32, including TAAR13c, while most (84/91) have lost Asp3.32. In several other fish species such as medaka, stickleback, fugu, and salmon, most TAARs also contain Asp5.42 but lack Asp3.32 (Figure 2—figure supplement 1). Phylogenetic analysis indicates that TAARs containing Asp5.42 but not Asp3.32 are clade III TAARs found in teleosts, and constitute a large branch of the TAAR family tree (Figure 2, Figure 2—figure supplement 2). Zebrafish lack clade II TAARs, and the teleost TAAR family evolved differently from mammalian TAARs by forming the large cohort of clade III TAARs. Several TAARs with both Asp5.42 and Asp3.32 reside between clade I TAARs and clade III TAARs phylogenetically (including TAAR13s in zebrafish), suggesting that clade III TAARs derived from clade I TAARs by first gaining Asp5.42 and then losing Asp3.32.
Ligands for clade III TAARs were unknown, and since they lack the canonical amine recognition site, were proposed to detect chemicals other than amines (Hussain et al., 2009; Liberles, 2014). However, the widespread retention of Asp5.42 in clade III TAARs lacking Asp3.32 raised the possibility that the corresponding 84 olfactory receptors in zebrafish indeed detect amines, but do so with ligands positioned in a non-classical orientation.
Identifying ligands for eleven 'orphan' zebrafish TAARs
We used the high throughput reporter gene assay to identify agonists for additional zebrafish TAARs (Figure 3, Figure 3—figure supplement 1). Responses of 63 zebrafish TAARs were examined, including at least one representative from each TAAR subfamily in the zebrafish olfactory system (TAAR10 to TAAR20). Zebrafish TAARs were expressed as fusion proteins containing the N-terminal 20 amino acids of bovine rhodopsin, which enhances cell surface expression of some chemosensory receptors in heterologous cells (Krautwurst et al., 1998). As above, HEK-293 cells or Hana3A cells (Saito et al., 2004) (Hana3A cells express olfactory chaperones for promotion of receptor expression and were used for studies involving TAAR10b, TAAR12i, and TAAR16e) were transfected with plasmids encoding zebrafish TAARs and Cre-Seap, incubated with test chemicals, and assayed for reporter phosphatase activity. TAAR ligand identification in heterologous cells has been extensively validated; TAAR ligands identified by this technique also activate TAAR-containing neurons in vivo (Hussain et al., 2013; Dewan et al., 2013; Li et al., 2013;Zhang, et al., 2013 ), and behavioral responses to identified TAAR ligands are lost in TAAR knockout mice (Dewan et al., 2013; Li et al., 2013). Here, we used the heterologous expression approach to identify the first ligands for zebrafish TAAR10a, TAAR10b, TAAR12h, TAAR12i, TAAR13a, TAAR13d, TAAR13e, TAAR14d, TAAR16c, TAAR16e, and TAAR16f.
Ligands for clade I TAARs. TAAR10a, TAAR10b, TAAR12h, and TAAR12i are clade I TAARs that retain Asp3.32 but lack Asp5.42. TAAR10a sensitively detects serotonin (EC50 = ~0.5 μM), has significantly reduced affinity for the related indole amines 5-methoxytryptamine (EC50 = ~20 μM) and tryptamine (EC50 = ~70 μM), and does not detect other chemicals examined (Figure 3—figure supplement 2). The highest affinity ligand identified for TAAR12h was 2-phenylethylamine (EC50 = ~0.3 μM), and TAAR12h also detects other alkyl amines with reduced sensitivity (Figure 3—figure supplement 2). TAAR10b and TAAR12i detect the primary amines tryptamine and 3-methoxytyramine respectively (Figure 3—figure supplement 1).
Serotonin and 2-phenylethylamine are both biogenic amines that can be recognized by other GPCRs. Zebrafish lack an ortholog of mouse TAAR4, which detects 2-phenylethylamine (Ferrero et al., 2011), and mouse TAAR4 shares 41% identity with zebrafish TAAR12h. Likewise, zebrafish TAAR10a is distantly related to mouse serotonin receptors, sharing 27–35% identity. It is thought that the TAAR family derived from an ancestral duplicate of serotonin receptor subtype HTR4 (Lindemann and Hoener, 2005), and it is possible that either TAAR10a retains the ligand recognition properties of the ancestral TAAR or that serotonin responsiveness was lost and then re-gained in this lineage.
Ligands for putative di-cation receptors. Seven zebrafish TAARs retain both Asp3.32 and Asp5.42, including all five TAAR13s and two TAAR14s (TAAR14c and TAAR14d). We previously found that zebrafish TAAR13c detects the diamine cadaverine (Hussain et al., 2013), and now report that zebrafish TAAR13a, TAAR13d, TAAR13e, and TAAR14d also detect chemicals with two cations.
TAAR13d prefers medium-chained alkyl diamines with highest affinity for putrescine (EC50 = ~1 μM), and progressively reduced affinity for diamines with additional methylene groups (Figure 4A). TAAR13d has ~3,000-fold reduced sensitivity for diaminopropane, which is only one methylene group shorter than putrescine, indicating a strict minimal agonist length. The differing preference of TAAR13c and TAAR13d for cadaverine and putrescine respectively is consistent with cross-adaptation studies indicating that separate olfactory receptors detect these two diamines (Rolen, et al., 2003).
Histamine is an agonist for both TAAR13a (EC50 = ~20 μM) and TAAR13d (EC50 = ~7 μM). TAAR13d recognizes other alkyl diamines of medium length (see above), while TAAR13a is selective for histamine among tested ligands. TAAR13a and TAAR13d share only 18–28% identity with mouse histamine receptors, and the recognition of histamine by olfactory receptors is likely due to convergent evolution within the GPCR family. Agmatine also activates multiple TAARs, including TAAR13e with higher affinity (EC50 = ~1 μM) and TAAR14d with reduced affinity (EC50 = ~100 μM). Agmatine is derived from arginine by decarboxylation, reminiscent of several other TAAR ligands that are biogenic amines similarly derived from natural amino acids (Liberles, 2015). TAARs that detect agmatine, histamine, serotonin, 2-phenylethylamine, or isoamylamine (see below) are not activated by the natural amino acids from which their ligands were derived (Figure 4—figure supplement 1).
Ligands for clade III TAARs. Next, we found ligands for three clade III TAARs that lack Asp3.32. Despite the absence of the classical amine recognition motif present in all known amine-detecting GPCRs, zebrafish TAAR16c, TAAR16e, and TAAR16f detect amines. Zebrafish TAAR16c detects N-methylpiperidine (EC50 = ~10 μM), zebrafish TAAR16e detects N,N-dimethylcyclohexylamine (EC50 = ~30 μM), and zebrafish TAAR16f detects isoamylamine (EC50 = ~200 μM).
Removing the amino group of TAAR16c and TAAR16f ligands abolished responses (Figure 4, Figure 4—figure supplement 2). TAAR16c also detected N-methylpyrrolidine (EC50 = ~70 μM), and the corresponding secondary amines piperidine and pyrrolidine were partial agonists. In contrast, oxygen-containing analogs, tetrahydropyran and tetrahydrofuran did not activate TAAR16c at any concentration tested, indicating that the amine is essential for receptor agonism. Likewise, TAAR16f failed to detect isoamyl alcohol, and thus TAAR16f agonism also required a ligand amine. These studies indicate that these clade III TAARs are in fact amine receptors, and that the ligand amino group is an essential moiety recognized by the receptor.
Transmembrane aspartates are required for TAAR agonism
We showed that Asp3.32 and Asp5.42 are essential for cadaverine detection by TAAR13c (Figure 1), and here asked if the corresponding residues are required for activation of other fish TAARs (Figure 5). Charge-neutralizing mutation of Asp3.32 in clade I TAARs, TAAR10a and TAAR12h, completely eliminated responses to serotonin and 2-phenylethylamine respectively, while charge-neutralizing mutation of Asp5.42 in clade III TAARs, TAAR16c and TAAR16f, completely eliminated responses to N-methylpiperidine and isoamylamine respectively. Similarly, mutation of either Asp3.32 or Asp5.42 in TAAR13a, TAAR13d, TAAR13e, and TAAR14d reduced or abolished responses to all di-cationic ligands. The asymmetric diamines histamine and agmatine could be positioned in either orientation, with the primary amine contacting either Asp3.32 or Asp5.42. These mutagenesis studies reveal that clade I TAARs and clade III TAARs require different transmembrane aspartates for amine recognition.
Homology models of serotonin-bound TAAR10a and N-methylpiperidine-bound TAAR16c were generated (Figure 6), as described previously for cadaverine-bound TAAR13c (Figure 1A). TAAR10a is a clade I TAAR and homology modeling indicated proximity of the serotonin amino group to Asp1173.32. In contrast, TAAR16c is a clade III TAAR, and homology modeling indicated ligand inversion with the N-methylpiperidine amino group near Asp1925.42. Inverting the ligand orientation in clade III TAARs would be expected to impact other ligand-receptor contacts; to test this model, we attempted to engineer such a ligand inversion in a canonical biogenic amine receptor, serotonin receptor 6 (HTR6). We analyzed the ligand preference of an HTR6 double mutant (D3.32A; A5.42D) in which the canonical amine recognition site at Asp3.32 was neutralized and a non-canonical amine recognition site at Asp5.42 was introduced. Ligand-identification studies revealed that the HTR6 double mutant no longer recognized serotonin, but instead recognized a different amine, ethanolamine (Figure 6—figure supplement 1). Re-orienting the ligand amine towards Asp5.42 presumably alters other receptor interactions dramatically, allowing Asp5.42-containing receptors to sample structurally distinct ligands. Together, these studies indicate that Asp5.42 forms a non-classical amine interaction site in GPCRs, and this site is present throughout a large clade of olfactory receptors.
Discussion
Sensory receptors define the capacity of an organism to perceive its external environment. Olfactory and taste receptor families can rapidly evolve, with dramatic species-specific variations in repertoire size and function (Nei et al., 2008; Shi and Zhang, 2007). In some lineages, receptor families can be reduced or die out, such as vomeronasal receptors in humans (Liman, 2006; Zhang and Webb, 2003) or sweet receptors in carnivores (Jiang et al., 2012), while in other lineages new receptor families can be born, such as formyl peptide receptors in rodents (Liberles et al., 2009; Rivière et al., 2009). Functional evolution of chemoreceptor families often involves a pattern of gene duplication and subsequent mutation (Ferrero et al., 2012), although isolated functional transformations can also occur, such as taste receptor re-purposing in hummingbirds that led to a novel receptor for sugars found in nectar (Baldwin et al., 2014).
Fish are aquatic animals, and efficient detection of water-soluble odors such as amines is important for their ecological niche. Zebrafish generally use the same families of olfactory receptors as other vertebrates, including odorant receptors, vomeronasal receptors, and TAARs (Kermen et al., 2013). The TAAR repertoire in mice and humans is relatively small, representing 1–2% of all main olfactory receptors. However, zebrafish have comparable numbers of TAARs (112) and odorant receptors (143), and expansion of the teleost TAAR family suggests a relatively expanded role for TAARs in fish chemosensation. However, prior to this study, ligands were known for only one fish TAAR (Hussain, et al., 2013). Moreover, the vast majority of zebrafish TAARs lack the canonical amine recognition motif, and the classes of ligands they might detect were unknown. Here, we show that almost all zebrafish TAARs (111/112) contain either a canonical amine-detection site on transmembrane α-helix III (Asp3.32), a non-canonical amine-detection site on transmembrane α-helix V (Asp5.42), or both. Thus, teleost TAARs likely represent a very large family of diverse amine receptors.
Several zebrafish TAAR ligands are biogenic amines produced during tissue decomposition (Shalaby, 1996). Cadaverine, histamine, agmatine, putrescine, tryptamine, 2-phenylethylamine, and isoamylamine are derived in one step from the natural amino acids lysine, histidine, arginine, ornithine, tryptophan, phenylalanine, and leucine by microbe-mediated decarboxylation. Levels of each of these amines are measured in commercially sold fish as indicators of food quality and freshness (Shalaby, 1996); for example, excessive histamine in consumed fish causes human scombroid poisoning (Morrow et al., 1991). While some of these biogenic amines are aversive to mammals (Dewan et al., 2013; Ferrero et al., 2011; Wisman and Shrira, 2015), they are reported to evoke variable behaviors in different fish species. For example, agmatine, putrescine, and cadaverine increase feeding-associated pecking behavior in goldfish (Rolen, et al., 2003; Hara, 2006) while cadaverine and putrescine evoke avoidance behavior in zebrafish (Hussain et al., 2013). Future studies are needed to examine the roles of each TAAR and its ligands in fish behavior.
The evolution of divergent solutions to the problem of amine detection by GPCRs reveals the extensive landscape of change possible within the GPCR superfamily. Here, we propose a two-step model for the functional evolution of clade III TAARs that is based on phylogenetic analysis, ligand identification studies, structural modeling, and mutagenesis (Figure 6). First, an ancestral clade I TAAR gained Asp5.42, and thus the ability to recognize di-cationic chemicals like diamines. Ligands were identified for five zebrafish TAARs that retain Asp3.32 and Asp5.42, and both aspartates are important for ligand recognition in each receptor. Second, a diamine receptor lost Asp3.32, leading to clade III TAARs that display non-classical monoamine recognition through a transmembrane α-helix V salt bridge to Asp5.42.
Flipping the amine orientation dramatically changes GPCR-ligand contacts, and presumably allows clade III TAARs to sample a distinct repertoire of structurally divergent amine odors. Once the new amine recognition motif was established, clade III TAARs dramatically expanded through a pattern of gene duplication and subsequent mutation. It seems likely that the massive expansion of this olfactory receptor clade increased the capacity of the fish olfactory system to detect amines, which are water-soluble and ecologically important stimuli for aquatic animals.
Materials and methods
Cloning of receptor genes
Request a detailed protocolZebrafish Taar genes were cloned from zebrafish genomic DNA (AB strain), and inserted into a modified pcDNA3.1- (Invitrogen, Waltham, MA, USA) vector containing a 5’ DNA extension encoding the first 20 amino acids of bovine rhodopsin followed by a cloning linker (GCGGCCGCC). Sequencing Taar genes revealed several AB strain polymorphisms as compared with Tübingen strain-derived genomic sequences deposited at NCBI. Amino acid-altering polymorphisms identified were as follows: Taar10a (NM_001082898): A467G (Ile to Val) and G545A (Val to Ile); Taar10b (NM_001082903): A338G (His to Arg), G496A (Ala to Thr), G544A (Val to Ile), and G775A (Val to Met); Taar12h (NM_001082909): A23G (Ile to Val), T26C (Tyr to His), C62T (Pro to Ser), C285A (Ser to Tyr), T743A (Phe to Ile), C874T and G875A (Gly to Ile); Taar12i (NM_001083085): T23C (Val to Ala), G313T, T314A (Val to Tyr), and T583G (Phe to Val); Taar13a (NM_001083102): C427G (Thr to Arg), T680A (Asn to Lys), T803G (Ile to Met), T882C (Phe to Leu), and T948C (Phe to Leu); Taar13d (NM_001083041): A112G (Ile to Val), A196C (Thr to Pro), G382A (Val to Ile), and G559A (Val to Ile); Taar13e (NM_001083043): A643G (Ile to Val) and G796A (Val to Ile); Taar16c (XM_009307036): A797C (Asn to Thr); Taar16e (XM_009307037): T103C (Val to Ala), G711A, C712T (Ala to Ile), and C844T (Ala to Val); Taar16f (XM_009305637): A23T (Asn to Ile), A59G (Asn to Ser), C190A (His to Lys), A202T (Thr to Ser), and A951C (Leu to Phe). Htr6 was cloned from mouse brain cDNA and inserted into pcDNA3.1-. Taar mutants were generated by overlap extension PCR and the Htr6 mutant was generated using TagMaster site-directed mutagenesis kit (GM Biosciences, Rockville, MD, USA).
TAAR functional assays
Request a detailed protocolReporter gene assays in HEK-293 cells or Hana3A cells (Saito et al., 2004) (for studies involving TAAR10b, TAAR12i, and TAAR16e) were performed as previously described (Liberles and Buck, 2006; Ferrero et al., 2012) with minor modifications. Modifications included use of tissue culture-treated 96 well plates (BD Biosciences, San Jose, CA, USA) pre-incubated with polylysine (250 ng per well); transfection using Lipofectamine 2000 (Invitrogen) reagent according to manufacturer's protocols, and introduction of test stimuli (see below) four to six hours after transfection.
Test stimuli for high throughput chemical screens
Request a detailed protocolTo identify TAAR agonists, chemicals were simultaneously tested in mixtures (83.3 μM per chemical, DMEM). Mix 1 contained GABA, β-alanine, 2-aminopentane, benzylamine, butylamine, and cystamine; mix 2 contained dibutylamine, dimethylamine, N,N-dimethylcyclohexylamine, ethylamine, ethyleneamine, and ethanolamine; mix 3 contained hexylamine, isoamylamine, isobutylamine, isopropylamine, indole, and cadaverine; mix 4 contained methylamine, N-methylindole, N-methylpyrrolidine, N-methylpiperidine, pyrrolidine, and putrescine; mix 5 contained hexanal, ethyl butyrate, tryptamine, histamine, 2-phenylethylamine, and trimethylamine; mix 6 contained 2-methylbutylamine, cysteamine, 1-amino propan-2-ol, 3-methylthiopropylamine, and agmatine; mix 7 contained tyramine, N,N-dimethylethanolamine, octopamine, N,N-dimethylglycine, 5-hydroxyindole-3-acetic acid, and 3-methoxytyramine; mix 8 contained 5-methoxytryptamine, 4-methoxy phenethylamine, N,N-dimethylphenethylamine, 5-methoxy-N,N-dimethyltryptamine, N,N-dimethylaniline, N,N-dimethylisopropylamine; mix 9 contained 1-dimethylamino propanol, 2-dimethylaminoethanethiol, and 1-(2-aminoethyl)-pyrrolidine; mix 10 contained 1,7-diaminoheptane, 1,8-diaminooctane, 1,6-diaminohexane, 1,10-diaminodecane, 1,3-diaminopropane; and mix 11 contained spermine, spermidine, dopamine, serotonin, and adenine. TAARs were tested with all 11 chemical mixtures, except TAAR14, TAAR15, and TAAR16 family members which were tested with mixes 1–10.
TAAR phylogenetic analysis
Request a detailed protocolAmino acid sequences of 112 zebrafish TAARs, 15 mouse TAARs, 17 rat TAARs, 6 human TAARs, 12 biogenic amine receptors (zebrafish HTR2B, HRH2, and DRD2A; mouse HRH3, DRD3, HTR5, DRD1A, and ADRB1; and rat HTR2A, HRH2, DRD3, and ADRB2) and 5 odorant receptors (zebrafish OR103 and OR131, mouse OR121 and OR446, and rat ORI15) were aligned with MAFFT v7.017 (Katoh, 2002). For Figure 2—figure supplements 1, 2, published TAAR amino acid sequences from opossum, cow, chicken, frog, coelacanth, medaka, stickleback, tetraodon, fugu, Atlantic salmon, and elephant shark, as well as sea lamprey aminergic receptors, were analyzed (Hussain et al., 2009; Tessarolo et al., 2014). A phylogeny was inferred using a maximum likelihood approach implemented in RAxML (RAxML v7.7.1) using the JTT + Γ model of codon substitution, empirical base frequencies, and the rapid bootstrapping technique (100 replicates) (Stamatakis et al., 2008). Nodes with support values < 50 were collapsed in the program TreeGraph 2 and trees were visualized using FigTree v1.3.1 (Stöver and Müller, 2010).
TAAR structural modeling
Request a detailed protocolThe TAAR13c structural model was generated using SWISS-MODEL workspace (Biasini et al., 2014) (http://swissmodel.expasy.org/interactive) and based on the X-ray crystal structure of agonist-bound human β2 adrenergic receptor (Protein Data Bank entry 4LDE) (Ring et al., 2013). Modeling involved an active-state template, in which agonist binding is favored. The sequence alignment for homology modeling was verified by inspection to confirm proper alignment of sequence landmarks including conserved sequence motifs (DRY/DRH and NPxxY) and transmembrane proline residues. Homology models for TAAR10a (37% amino acid identity) and TAAR16c (30% amino acid identity) were prepared similarly using the same template. Models of receptor point mutants were prepared by side chain substitution in previously prepared homology models and then subjected to energy minimization prior to ligand docking. Docking of ligands into homology models was conducted with the Schrödinger suite (Schrodinger, 2015) using LigPrep version 3.4 and Glide version 6.7. First, the receptor model was prepared for docking by adding protons and assigning charges to ionizable side chains. Asp3.32 and Asp5.42 were modeled in the deprotonated state expected at physiological pH. The resulting model was then subjected to energy minimization using the Schrödinger impref utility prior to docking. Ligands were docked to the orthosteric binding pocket using Glide Extra Precision (Friesner et al., 2006), and figures were prepared in PyMol (Schrodinger, 2010).
References
-
SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary informationNucleic Acids Research 42:W252–W258.https://doi.org/10.1093/nar/gku340
-
Trace amines: identification of a family of mammalian G protein-coupled receptorsProceedings of the National Academy of Sciences of the United States of America 98:8966–8971.https://doi.org/10.1073/pnas.151105198
-
Amphetamine, 3,4-methylenedioxymethamphetamine, lysergic acid diethylamide, and metabolites of the catecholamine neurotransmitters are agonists of a rat trace amine receptorMolecular Pharmacology 60:1181–1188.
-
Detection and avoidance of a carnivore odor by preyProceedings of the National Academy of Sciences of the United States of America 108:11235–11240.https://doi.org/10.1073/pnas.1103317108
-
Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexesJournal of Medicinal Chemistry 49:6177–6196.https://doi.org/10.1021/jm051256o
-
Olfactory receptor patterning in a higher primateThe Journal of Neuroscience 34:12241–12252.https://doi.org/10.1523/JNEUROSCI.1779-14.2014
-
High-affinity olfactory receptor for the death-associated odor cadaverineProceedings of the National Academy of Sciences of the United States of America 110:19579–19584.https://doi.org/10.1073/pnas.1318596110
-
Positive darwinian selection and the birth of an olfactory receptor clade in teleostsProceedings of the National Academy of Sciences of the United States of America 106:4313–4318.https://doi.org/10.1073/pnas.0803229106
-
Major taste loss in carnivorous mammalsProceedings of the National Academy of Sciences of the United States of America 109:4956–4961.https://doi.org/10.1073/pnas.1118360109
-
MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transformNucleic Acids Research 30:3059–3066.https://doi.org/10.1093/nar/gkf436
-
Structure-function of the g protein-coupled receptor superfamilyAnnual Review of Pharmacology and Toxicology 53:531–556.https://doi.org/10.1146/annurev-pharmtox-032112-135923
-
Neural circuits mediating olfactory-driven behavior in fishFrontiers in Neural Circuits 7:62.https://doi.org/10.3389/fncir.2013.00062
-
Aversion and attraction through olfactionCurrent Biology 25:R120–R129.https://doi.org/10.1016/j.cub.2014.11.044
-
Formyl peptide receptors are candidate chemosensory receptors in the vomeronasal organProceedings of the National Academy of Sciences of the United States of America 106:9842–9847.https://doi.org/10.1073/pnas.0904464106
-
Mammalian pheromonesAnnual Review of Physiology 76:151–175.https://doi.org/10.1146/annurev-physiol-021113-170334
-
Trace amine-associated receptors: ligands, neural circuits, and behaviorsCurrent Opinion in Neurobiology 34:1–7.https://doi.org/10.1016/j.conb.2015.01.001
-
Use it or lose it: molecular evolution of sensory signaling in primatesPflügers Archiv : European Journal of Physiology 453:125–131.https://doi.org/10.1007/s00424-006-0120-3
-
A renaissance in trace amines inspired by a novel GPCR familyTrends in Pharmacological Sciences 26:274–281.https://doi.org/10.1016/j.tips.2005.03.007
-
Evidence that histamine is the causative toxin of scombroid-fish poisoningThe New England Journal of Medicine 324:716–720.https://doi.org/10.1056/NEJM199103143241102
-
The evolution of animal chemosensory receptor gene repertoires: roles of chance and necessityNature Reviews. Genetics 9:951–963.https://doi.org/10.1038/nrg2480
-
TAAR1 activation modulates monoaminergic neurotransmission, preventing hyperdopaminergic and hypoglutamatergic activityProceedings of the National Academy of Sciences of the United States of America 108:8485–8490.https://doi.org/10.1073/pnas.1103029108
-
Polyamines as olfactory stimuli in the goldfish carassius auratusJournal of Experimental Biology 206:1683–1696.https://doi.org/10.1242/jeb.00338
-
The PyMOL molecular graphics system, version 1.3r1
-
Schrodinger release 2015-2
-
Molecular and cellular analysis of human histamine receptor subtypesTrends in Pharmacological Sciences 34:33–58.https://doi.org/10.1016/j.tips.2012.11.001
-
Significance of biogenic amines to food safety and human healthFood Research International 29:675–690.https://doi.org/10.1016/S0963-9969(96)00066-X
-
The binding site of aminergic g protein-coupled receptors: the transmembrane segments and second extracellular loopAnnual Review of Pharmacology and Toxicology 42:437–467.https://doi.org/10.1146/annurev.pharmtox.42.091101.144224
-
A rapid bootstrap algorithm for the RAxML web serversSystematic Biology 57:758–771.https://doi.org/10.1080/10635150802429642
-
The smell of death: evidence that putrescine elicits threat management mechanismsFrontiers in Psychology 6:1274.https://doi.org/10.3389/fpsyg.2015.01274
-
Ultrasensitive detection of amines by a trace amine-associated receptorJournal of Neuroscience 33:3228–3239.https://doi.org/10.1523/JNEUROSCI.4299-12.2013
-
Evolutionary deterioration of the vomeronasal pheromone transduction pathway in catarrhine primatesProceedings of the National Academy of Sciences of the United States of America 100:8337–8341.https://doi.org/10.1073/pnas.1331721100
Article and author information
Author details
Funding
National Institute for Health Research (R01 DC013289)
- Stephen D Liberles
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Xianchi Dong for assistance with GPCR modeling, Hiroaki Matsunami for Hana3A cells, and Sean Megason for zebrafish genomic DNA. This work was supported by a grant from the NIH (SDL, Award Number R01 DC013289).
Copyright
© 2015, Li et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,971
- views
-
- 672
- downloads
-
- 41
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Cell Biology
Activation of the Wnt/β-catenin pathway crucially depends on the polymerization of dishevelled 2 (DVL2) into biomolecular condensates. However, given the low affinity of known DVL2 self-interaction sites and its low cellular concentration, it is unclear how polymers can form. Here, we detect oligomeric DVL2 complexes at endogenous protein levels in human cell lines, using a biochemical ultracentrifugation assay. We identify a low-complexity region (LCR4) in the C-terminus whose deletion and fusion decreased and increased the complexes, respectively. Notably, LCR4-induced complexes correlated with the formation of microscopically visible multimeric condensates. Adjacent to LCR4, we mapped a conserved domain (CD2) promoting condensates only. Molecularly, LCR4 and CD2 mediated DVL2 self-interaction via aggregating residues and phenylalanine stickers, respectively. Point mutations inactivating these interaction sites impaired Wnt pathway activation by DVL2. Our study discovers DVL2 complexes with functional importance for Wnt/β-catenin signaling. Moreover, we provide evidence that DVL2 condensates form in two steps by pre-oligomerization via high-affinity interaction sites, such as LCR4, and subsequent condensation via low-affinity interaction sites, such as CD2.
-
- Biochemistry and Chemical Biology
The development of proteolysis targeting chimeras (PROTACs), which induce the degradation of target proteins by bringing them into proximity with cellular E3 ubiquitin ligases, has revolutionized drug development. While the human genome encodes more than 600 different E3 ligases, current PROTACs use only a handful of them, drastically limiting their full potential. Furthermore, many PROTAC development campaigns fail because the selected E3 ligase candidates are unable to induce degradation of the particular target of interest. As more and more ligands for novel E3 ligases are discovered, the chemical effort to identify the best E3 ligase for a given target is exploding. Therefore, a genetic system to identify degradation-causing E3 ligases and suitable target/E3 ligase pairs is urgently needed. Here, we used the well-established dimerization of the FKBP12 protein and FRB domain by rapamycin to bring the target protein WDR5 into proximity with candidate E3 ligases. Strikingly, this rapamycin-induced proximity assay (RiPA) revealed that VHL, but not Cereblon, is able to induce WDR5 degradation - a finding previously made by PROTACs, demonstrating its predictive power. By optimizing the steric arrangement of all components and fusing the target protein with a minimal luciferase, RiPA can identify the ideal E3 for any target protein of interest in living cells, significantly reducing and focusing the chemical effort in the early stages of PROTAC development.