# The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding

Tammy Shim, Jody Pacalon, Won-Cheol Kim, Xiaojing Cong, Jérémie Topin,

Jérôme Golebiowski, Cheil Moon

► To cite this version:

Tammy Shim, Jody Pacalon, Won-Cheol Kim, Xiaojing Cong, Jérémie Topin, et al.. The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding. International Journal of Molecular Sciences, 2022, 23 (20), pp.12501. 10.3390/ijms232012501 . hal-04712540

HAL Id: hal-04712540

https://hal.science/hal-04712540v1

Submitted on 27 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Article

# The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding

Tammy Shim <fnref n="1" />,†,✉, Jody Pacalon <fnref n="2" />,†,✉, Won-Cheol Kim <fnref n="1" />, Xiaojing Cong <fnref n="3" />,✉, Jérémie Topin <fnref n="2" />,*,✉, Jérôme Golebiowski <fnref n="1" /> and Cheil Moon <fnref n="1" />,*,✉

**Abstract:** Mammals recognize chemicals in the air via G protein-coupled odorant receptors (ORs). In addition to their orthosteric binding site, other segments of these receptors modulate ligand recognition. Focusing on human hOR1A1, which is considered prototypical of class II ORs, we used a combination of molecular modeling, site-directed mutagenesis, and in vitro functional assays. We showed that the third extracellular loop of ORs (ECL3) contributes to ligand recognition and receptor activation. Indeed, site-directed mutations in ECL3 showed differential effects on the potency and efficacy of both carvones, citronellol, and 2-nonanone.

**Keywords:** odorant receptors; ECL3; ligand selectivity; molecular modeling; functional assays

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Shim, T.; Pacalon, J.; Kim, W.-C.; Cong, X.; Topin, J.; Golebiowski, J.; Moon, C. The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding. *Int. J. Mol. Sci.* **2022**, *23*, 12501. https://doi.org/10.3390/ijms232012501

Academic Editor: Peihua Jiang

Received: 10 August 2022

Accepted: 6 October 2022

Published: 18 October 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

Mammals rely on their sense of smell to assess the volatile chemical environment. Smell information is decoded by chemical interactions between odorants and G protein-coupled odorant receptors (ORs) expressed in olfactory neurons present in the nasal cavity. While these transmembrane proteins are the cornerstone of chemical recognition by our neurons, perireceptor events can modify the chemical composition of odorants before reaching the ORs [1]. For example, the enzyme carboxyl esterase has been shown to influence odor recognition by altering the chemical function [2]. Odorant-binding proteins also modulate the chemical signal by helping to solubilize odorants or playing the role of scavengers [3].

Nonetheless, the coding of a perceived odor by the olfactory system relies on a combinatorial code where ORs are differentially activated by odorants and where one odorant can activate multiple ORs [4,5]. The subtle interactions between ORs and odorants at the molecular level remain extremely difficult to rationalize. In fact, very subtle modifications in the chemical structure of a molecule can drastically alter its odor [6,7]. Similarly, mutations within an OR gene might strongly alter smell perception [8].

ORs represent more than 3% of the whole proteome and belong to the class A G protein-coupled receptor (GPCR) family of proteins, which are responsible for transmitting signals through the cell membrane. The large number of ORs (~400 in humans and ~1100 in mice, for example) coupled with this combinatorial activation endows mammals with incredible discriminatory power [9,10]. GPCRs are one of the largest and most diverse membrane protein families. They adopt a typical architecture consisting of seven transmembrane helices (TM1 to TM7) linked by intracellular and extracellular loops (ICLs and ECLs, respectively). ECLs, while peripheral within the GPCR tertiary structure, are involved in numerous receptor functions such as ligand or protein recognition and receptor activation [11–14]. In most class A GPCRs, the internal binding site for the ligand is found

<footnote n="1">¹ Department of Brain Sciences, Daegu Gyeongbuk Institute of Science and Technology, Daegu 711-873, Korea</footnote>

<footnote n="2">² Institut de Chimie de Nice UMR7272, Université Côte d'Azur, CNRS, 06108 Nice, France</footnote>

<footnote n="3">³ Institut de Génomique Fonctionnelle, Université de Montpellier, CNRS, INSERM, CEDEX 05, 34094 Montpellier, France</footnote>

<footnote>* Correspondence: jeremie.topin@univ-cotedazur.fr (J.T.); cmoon@dgist.ac.kr (C.M.)</footnote>

<footnote>† These authors contributed equally to this paper.</footnote>
at ca. 10 Å with respect to the extracellular side of the receptor. In the prototypical beta2-adrenergic receptor, multiple interactions have been described between the ligand and ECL3 during the process of ligand migration from the bulk solvent to the internal binding site [15]. In particular, ECL3 was suggested as a functional region important for ligand specificity. In a sub-class of mammalian ORs, namely, the class I ORs, we showed that the extracellular part of the receptor, notably ECL3, played the role of a vestibular binding site [16]. The sequence variability in the ECL3 is high between the two OR classes, which prevents us from making any conclusions on the role of ECL3 in class II ORs based on results from class I. This function of ECL3 in the class II OR has not been investigated until now.

In this article, we report on the role of ECL3 in class II ORs and show that ECL3 modulates ligand binding independent of the sequence variability between the two classes. We consider hOR1A1 as the prototypical class II OR. We thus conclude that ECL3 is involved in odorant selectivity in all mammalian ORs.

## 2. Results

### 2.1. *The Highly Variable ECL3 Sequence Acts as a Vestibular Binding Pocket for Ligands*

In class A GPCRs, ECL3 connects the extracellular parts of TM6 and TM7. Up to now, no structure of hOR1A1 has been experimentally solved. A 3D modeled structure was built using a previously established protocol [17,18]. An alternative model was obtained using AlphaFold [19]. Figure 1a depicts the two modeled 3D structures of hOR1A1. The overall structure of the two models is conserved with an RMSD of 3.5 Å on the backbone. In general, the main difference occurs in the unstructured segment of the receptor. In both models, ECL3 was located ca. 4 Å toward the extracellular side of the orthosteric binding cavity. In hOR1A1, ECL3 comprised five/six residues (P261/L262 to S266) and showed an unfolded 3D structure. The vicinity of ECL3 and the binding cavity was consistent with a potential role of ECL3 as a vestibular binding site as already observed by us on class I ORs, and by others on the beta2-adrenergic receptor and muscarinic M2 receptors [15,16,20].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** (a) Structure of hOR1A1 from homology modeling (Modeller) compared to that obtained from AlphaFold2 (AlphaFold). In both structures, the third extracellular loop (ECL3) (shown in red) was predicted to be close to the orthosteric binding cavity, shown as a cyan surface. (b) Conservation analysis of ECL3 sequences of both classes of human odorant receptors and the highlight of hOR1A1 specific ECL3 sequence.</figcaption></figure>

From a sequence point of view, the conservation analysis and its decomposition into class I and class II ORs revealed that the positions in ECL3 were conserved in class II ORs and formed a P261xSxxS motif (Figure 1b). In human ORs, the proline (here P261) was conserved at 72.3%, S263 at 46.5%, and S266 at 44.1%. In hOR1A1, the sequence reads PLTNYS. They were not only different from those found in class I but also different from any other class A GPCRs.

The involvement of ECL3 in the binding process was assessed by a molecular dynamics (MD) simulation where (–)-carvone, a strong ligand for hOR1A1 [21,22], was initially
located within the bulk solvent. When all replicas were aggregated, the binding mechanism
could be decomposed into a three-step process. Figure 2 depicts the main events identified
during the trajectories and highlights the density of (−)-carvone in the protein with respect
to both residue Y251<sup>6.48</sup> (bottom of the orthosteric cavity) and S266<sup>ECL3</sup>. Starting from the
bulk solvent, the ligand rapidly approached the extracellular segment of the receptor and
initiated contact with ECL3 before reaching the orthosteric binding site. The density map
confirmed the three regions in which the ligand spent more time. Upon the binding process,
the ligand interacted with various residues in ECL3, highlighting this loop's involvement
in the ligand's entry within the receptor.

<figure><img src="image_6.png" /></figure>

**Figure 2.** Entry of (−)-carvone inside receptor hOR1A1. The ligand is initially located outside the receptor (1). It then migrates to the cradle of the orthosteric binding cavity (2,3), as indicated by Y251<sup>6.48</sup>. During this process, the ligand interacts with several residues from ECL3 (indicated in red). Contour map of (−)-carvone migration as the minimum distance from S266 (taken as the distance from ECL3) and minimum distance from Y251<sup>6.48</sup> (taken as the distance from the cradle of the cavity). All replicas were considered. The three highlighted basins show the ligand's largest density.

In our homology model, P261<sup>ECL3</sup> is the first residue of ECL3 and acts as a (−)-carvone contact point with the receptor after ca. 100 nanoseconds. The ligand then interacted with the hydrophobic residues Leu262<sup>ECL3</sup> and diffused towards the orthosteric binding site, engaging H-bonds with Ser266<sup>ECL3</sup> (Figure 2, position 2). Finally, the ligand interacted with Asp269<sup>7.34</sup> before entering the binding cavity (Figure 2, on the left, position 3), i.e., at less than 2.5 Å of Y251<sup>6.48</sup> (6.48 refers to the Ballesteros–Weinstein notation) [23,24]. To assess the role of amino acids from ECL3 in the receptor recognition process, in vitro experiments were also performed.

## 2.2. *In Vitro Functional Assays Highlight the Differential Interactions with Diverse Ligands*

Based on the results of both the conservation analysis and the MD simulations, several amino acids from ECL3 seem to be involved in the binding process. To investigate the specific chemical functions of amino acids in ECL3 upon ligand binding, we designed mutant ORs with various (small, charged, lipophilic, or aromatic) properties, and assessed their response to four different ligands (Figure 3).
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Chemical structure of four agonists of hOR1A1.</figcaption></figure>

All the mutant ORs investigated in this study were confirmed to be expressed at the membrane surface and showed some basal activities, confirming that they remain functional (Figure 4b,c and Supplementary Materials). However, it appeared that the expression level of the mutant ORs was differentially affected compared to the *wt* OR (Figure 4b,c). The presence of a phenylalanine residue systematically resulted in a decrease in surface expression for all four mutant ORs. Mutation to an alanine residue at all three hydrophilic positions (T263ECL3, S266ECL3 and D269) also decreased surface expression of the receptor, whereas this mutation had no effect on expression for the P261A mutant OR. Surprisingly, the presence of a positive charge at positions P261ECL3, T263ECL3, and S266ECL3 did not affect surface expression. The same modification at position D269 induced a decrease in surface expression. Overall, we observed that bulky amino acids in ECL3 resulted an increased basal activity (Figure S3) independent of surface expression.

The *wt* OR exhibited a dose-dependent increase in cAMP levels in response to four ligands ((-)-carvone, (+)-carvone, citronellol, and 2-nonanone) known as bona fide agonists in vitro (Figures 3 and 4a) [21,23,25]. Mutations at positions P261, T263, S266, and D269 differentially affected the receptor response to agonists (Figure 4a).

MD simulations revealed interactions between the (–)-carvone and P261<sup>ECL3</sup> early in the binding process. P261X mutant ORs (X = A, R, L, or F) did not show a dose-dependent response for all four ligands: their efficacy and potency were significantly reduced (Figure 4a). Mutations affected the agonist-induced response differently: P261A and P261R remained sensitive to the most potent ligand, (–)-carvone, whereas the receptors’ responses were negligible for the other three ligands. Mutations to a more hydrophobic and bulkier residue, such as leucine or phenylalanine (Figure 4a pink and green, respectively), eliminated the response to all four ligands. These data show that mutations at position P261 differentially affect the response to ligands and reinforce the hypothesis of the initial interaction between P261 and ligands during the binding process.

T263 and S266 are two small and polar conserved residues in the ECL3 P*S*S*S motif, which strongly interacted with (–)-carvone during the MD simulations. Notably, mutations at position 263 have been reported to affect OR2AG1 response to agonist [23]. Interestingly, mutations at these positions had a different impact on the receptor response to the 4 agonists. The mutations into apolar (A), charged (R), or aromatic (F) residues were—as for P261—also associated with differential modification of the receptor response to ligands in vitro.

The T263F mutation abolished the receptor response to all four ligands (Figure 4a, second row, green curves). The presence of positively charged residues systematically decreased the efficacy of all four ligands (Figure 4a, second row, blue curve). Finally, all four ligands were differentially affected by introducing an alanine at position 263. The response to 2-nonanone and citronellol remained unchanged compared to *wt* OR, while the efficacy of the two enantiomers of carvone decreased. This observation suggested a specific interaction of residue T263 with both carvone enantiomers.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** In vitro data of hOR1A1 and mutant ORs. (a) In vitro dose-response curves of four ligands (–)-carvone, (+)-carvone, citronellol, and 2-nonanone towards *wt* hOR1A1 and mutant ORs at positions P261, T263, S266, and D269. (*) indicates the response value is significantly different compared to *wt* hOR1A1 (one-way ANOVA, followed by a Dunnett test; \* *p* < 0.05). (RLU = relative luminescence unit) (b) The fluorescence intensity of the R-phycoerythrin (R-PE) signal of *wt* hOR1A1 and mutant ORs at positions P261, T263, S266, and D269. (c) Normalized graph of cell-surface expression of ECL3 mutant ORs against *wt* hOR1A1 (one-way ANOVA, followed by a Dunnett test; \* *p* < 0.05, \*\* *p* < 0.01, and \*\*\* *p* < 0.001).</figcaption></figure>
Mutations at position S266 also induced changes in efficacy. Surprisingly, (-)-carvone efficacy for the S266R mutant OR remained unchanged (Figure 4a, 3rd row, black curve vs. blue curve), whereas it was systematically reduced for the other ligands. The S266A and S266F mutant ORs decreased the efficacy of all four ligands. As for the other three residues investigated, introducing a phenylalanine residue almost abrogated the receptor response.

D269<sup>7.34</sup> was not predicted to belong to ECL3 in both the Modeller and AlphaFold models. It was located at the upper extremity of TM7. One can observe in Figure 4c that surface expression of the mutant ORs was systematically decreased except for D269L, in line with the role of the extracellular part of TM7 in the activation mechanism [26]. Compared to the mutant ORs at ECL3, the impact of the mutation on receptor response to agonist stimulation was more systematic. The mutations to A, R, or F induced a similar decrease in the receptor response for the four studied ligands. The D269L mutant OR showed a much more pronounced decrease in response to agonist stimulation except for citronellol (Figure 4a, last line, pink color). This systematic trend suggests that D269 seems to be involved in the dynamics of receptor activation rather than in ligand-binding recognition.

All these experiments underline the importance of ECL3 residues during binding. It can be observed that the presence of a phenylalanine at the position studied systematically altered the expression level. Furthermore, the response of the phenylalanine mutant ORs to the four ligands was abrogated, except for D269F, which remained sensitive to agonist stimulation.

In general, the mutations did not alter the EC50. When a dose-response curve was observed, the EC50 never shifted more than one order of magnitude in concentration. The mutations clearly affect the efficacy of the ligands. Finally, the major takeaway of our study is that from a general point of view, the efficacy and potency of ligands are differentially affected by the mutations in ECL3.

### 3. Discussion

Odorant receptors, as all class A GPCRs, are structured with a seven-transmembrane domain where helices are connected by three extra- and three intra-cellular loops. Among these loops, the short ECL3 connects TM6 to TM7 and shows a large sequence variability. Several works have already evaluated the role of this very poorly conserved ECL3. In particular, they assessed its essential contribution to the ligand-binding process of the beta2-adrenergic receptor or the muscarinic M2 receptor [15,20]. The active role of ECL3 in ligand binding is consistent with the close vicinity of ECL3 and the orthosteric binding site. This has been, for example, illustrated in some peptide-binding class A GPCRs, where the peptide directly interacts with both ECL3 and the orthosteric binding pocket [27,28]. Note that ECL3 structures in various peptide-bound class A GPCR are much more variable than ECL1 or ECL2, suggesting that ECL3 plays a modulation role rather than a ligand recognition role [29]. All in all, the role of ECL3 in class A GPCRs is not only structural but it seems to modulate the ligand's recognition in combination with the binding site.

In class I ORs, a sub-family of mammalian ORs, we previously proposed that ECL3 functions as a vestibule site and contributes to ligand binding [16]. In this work, we hypothesized a similar involvement of ECL3, although the conserved residues are significantly different between the two sub-families of ORs. Through a joint approach combining molecular modeling, heterologous functional expression, site-directed mutagenesis, and functional assays on hOR1A1—a prototypical class II OR—we have evaluated the role of ECL3 in ligand binding. In class II ORs, ECL3 is based on the conserved P<sup>261</sup> xS<sup>263</sup> xxS<sup>266</sup> motif. Various mutant ORs at positions 261, 263, and 266 showed differential responses to four agonists belonging to various chemical families compared to *wt*.

Meanwhile, additional mutant ORs at position D269<sup>7.34</sup> (at the junction between ECL3 and TM7) showed a conserved modulation regardless of the type of substituted amino acid. Molecular dynamics simulations of a homology model of hOR1A1 in interaction with (-)-carvone sampled structures where the ligand was in regular contact with ECL3 prior to entering the orthosteric binding cavity.
hOR1A1 has been extensively studied for its property of discriminating the two enantiomers of carvone [23]. In a former study combining homology modeling, site-directed mutagenesis, and functional expression studies, the authors identified eleven positions on the transmembrane segments involved in the chiral recognition of carvone. However, our results reveal that the two enantiomers of carvone are similarly affected by mutations in ECL3. The conserved residues in ECL3, although involved in agonist recognition, are apparently not involved in chiral discrimination in hOR1A1.

Finally, our results suggest that ECL3 in hOR1A1 plays a similar role in ligand binding as observed in the beta-2-adrenergic receptor, the M2 muscarinic receptor, and class I ORs. This study therefore reconciles all previous studies and demonstrates a similar function of ECL3 on ligand binding in all class A GPCRs, independent of the ECL3 sequence. The large variability in the ECL3 sequence among class A GPCR subfamilies is consistent with the diverse chemical space associated with these receptors, which can bind highly polar or highly lipophilic ligands. It could become an interesting and alternative target for allosteric ligand design.

# 4. Materials and Methods

## 4.1. *In Silico* Experiments

### 4.1.1. Molecular Modeling

hOR1A1 models were created using Modeller 9.25 [30]. Our previous protocol was described by De March et al. with four experimental structures as templates (PDB codes: 1U19, 3ODU, 2YDV, and 2LNL) [18]. Cysteine residues 74 with 156, and 146 with 166 were assigned to form two disulfide bonds [31]. The best-generated models were selected based on the DOPE score, and one model was finally chosen based on visual inspection. The hOR1A1 AlphaFold model was retrieved from the AlphaFold structure database [32,33].

Modeller and AlphaFold models have a 3.5 Å RMSD (calculated on the backbone of TMs). Both models nicely superimpose concerning TM6 and TM7. The distance between Y251⁶.⁴⁸ and ECL3's amino acids is comparable in the hOR1A1 Modeller and AlphaFold model (15.7/17.4 Å, 19.0/20.2 Å, 19.0/19.4 Å and 14.7/14.3 Å for residues P261, T263, S266 and D269, respectively, for Modeller/AlphaFold). While ECL3 is made up of almost identical residues, its structure is predicted with low confidence by AlphaFold (pLDDT between 50 and 70). The model built by homology was kept for the molecular dynamics simulations.

### 4.1.2. Molecular Dynamics Preparation

As the N-terminal of GPCR is not resolved in the template we used, the N-terminal part of the hOR1A1 Modeller model was discarded until residue E24. Propka3 was used to predict protonation states of the protein at a target pH 6.5 [34]. The extremities of the model were capped accordingly. hOR1A1 orientation in its membrane was determined using the OPM server [35]. Three (–)-carvone molecules were added in different orientations on the extracellular side at 5–10 Å of the top of the receptor. The system was embedded into a POPC-only model membrane using PACKMOL-Memgen [36]. The simulation box was completed using TIP3P water molecules and neutralized using $K^+$ and $Cl^-$ ions with a final concentration of 0.15 M. The total system comprised 49,536 atoms in a $7 \times 10^5$ Å³ periodic box.

### 4.1.3. Molecular Dynamics Protocol

Molecular dynamics simulations were performed with the sander and pmemd.cuda modules of AMBER18 [37], with the ff14SB force field for the proteins and the lipid14 forcefield for the membrane. (–)-carvone parameters were generated by calculating partial atomic charges with the HF/6–31G* basis set using Gaussian 09 [38]. The obtained electrostatic potential was fitted by the RESP program [39]. The other parameters were taken from the general Amber force field 2 (gaff2). Bonds involving hydrogen atoms were constrained using the SHAKE algorithm and long-range electrostatic interactions were
handled using particle mesh Ewald. The cut-off for non-bonded interaction was set to 10 Å. Each system was first minimized with the AMBER sander module, with 5000 steps of steepest descent and then 5000 steps of the conjugate gradient with a 50 kcal·mol⁻¹·Å² harmonic potential restraint on the protein. A second minimization of the same length without restraint was applied. The systems were then thermalized from 100 to 310 K for 10,000 steps (restraining the protein and ligands with a 200 kcal·mol⁻¹·Å² harmonic potential). Each system underwent 50,000 steps of equilibration in the NPT ensemble and 1 bar (restraining the protein and ligands with a 15 kcal·mol⁻¹·Å² harmonic potential) before the production phase. During the equilibration and production phase, the temperature was kept constant in the system at 310 K using a Langevin thermostat with a collision frequency of 5 ps⁻¹. To increase sampling, all 3 (−)-carvone molecules were constrained in a sphere of 30–40 Å radius, centered on the center of mass of the Thr89 hOR1A1 (with a potential of 10 kcal·mol⁻¹·Å²).

The system stability was evaluated from the root mean square deviation (RMSD) evolution computed on the receptor backbone. For the 8 replicas, the receptors underwent small fluctuations (RMSD < 4 Å), showing that they remained correctly folded during microsecond simulations (Figure S1). Binding events occurred in Rep1 (stable) and Rep3 (partial entry) (Figure S2). Cpptraj and Pytraj v2.02.dev0 were used for distance and RMSD analysis [40].

## 4.2. *In Vitro* Experiments

### 4.2.1. Cell Culture

HEK293T cell line was obtained from ATCC (#CRL-3216, ATCC, Manassas, VA, USA) and cultured in Dulbecco's modified Eagle's medium (DMEM; #10-017-CV, Corning, NY, USA) plus 10% fetal bovine serum (#16000-044, Thermo Fisher Scientific, Waltham, MA, USA) and penicillin/streptomycin (#15140122, Thermo Fisher Scientific, Waltham, MA, USA) at 37 °C with 5% CO₂.

### 4.2.2. DNA Constructs, Site-Directed Mutagenesis, and Gene Transfection

The full genomic DNA sequence of *hOR1A1* was obtained from HEK293T cells using cloning primers (5'-GCA CGC GTA TGA GGG AAA ATA ACC AGT C-3' and 5'-GCG CGG CCG CTT ACG AGG AGA TTC TCT T-3'). PCR product was subcloned to LUCY-FLAG-Rho tagged pCI mammalian expression vector using MluI and NotI restriction enzyme. The pCI vector, which is the Rho sequence-tagged at N-terminal, RTP1s-pCI, $G_{olf}$-pCI, and Ric8b-pCI was a kind gift from Dr. Matsunami (Duke U., Durham, NC, USA). A cleavable leucine-rich 17 amino acid signal peptide (LUCY; MRPQILLLLALLTLGLA) and FLAG sequences were tagged in front of the Rho sequence to promote or detect the functional expression on the plasma membrane, respectively [41]. pHIV-EGFP was a gift from Bryan Welm and Zena Werb (Addgene plasmid #21373). pGloSensorTM-22F cAMP Plasmid was purchased (#E2301, Promega, Madison, WI, USA). Site-directed mutant ORs were generated through the QuikChange PCR protocol [42] using the mutagenic oligonucleotide primers, designed individually according to the desired mutation. The sequence of plasmids was confirmed through BigDye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific, Waltham, MA, USA). Furthermore, JetPRIME (#114-75, PolyPlus-transfection, Illkirch-Graffenstaden, France) was used for gene transfection. For FACS experiments, the OR or mutant OR, RTP1s, Ric8b, $G_{olf}$, and pHIV-EGFP plasmid were transfected, and for cAMP luminescence assay, the pGloSensorTM-22F cAMP plasmid was transfected instead of the pHIV-EGFP plasmid.

### 4.2.3. cAMP Luminescence Assay

The cAMP luminescence assay was performed using the GloSensorTM cAMP assay (Promega, Madison, WI, USA) by following the manufacturer's protocol. In brief, the media of the transfected HEK293T cells on the 96-well white plate were changed to CO₂ independent media (#18045088, Thermo Fisher Scientific, Waltham, MA, USA) with GloSensorTM
cAMP Reagent (#E1291, Promega, Madison, WI, USA) before exposure to solvent or odors. Each odorant was diluted into CO₂ independent media and used to treat each well. The endpoint of the luminescence level was measured with a SpectraMax L Microplate Reader (Molecular Devices, San Jose, CA, USA). Data were analyzed through Microsoft Excel and GraphPad Prism. To compare OR responses plate to plate, empty pCI vector and wild-type OR were always included as a control. To read the basal activity of each OR, the values of at least six wells were averaged in the absence of odorants. The experiments were repeated twice to measure the odorant-induced OR activity, and each condition was triplicated. The measured luminescence value was further corrected by subtracting the value for the lowest response to each odorant of that receptor. The basal activity and odorant-induced responses were normalized to that of wild-type OR.

#### 4.2.4. Functional Expression

The functional expression of OR on the plasma membrane was measured through fluorescent-activated cell sorting (FACS) as previously described [43]. To summarize, transfected HEK293T cells on the six-well plates were gently detached from the plate using Cellstripper (#25-056-CI, Thermo Fisher Scientific, Waltham, MA, USA). The cells were incubated with Rho4D2 antibody (#ab98887, Abcam, Cambridge, MA, USA) for an hour with gentle rotation. After the cells were washed twice with washing and staining solution (2% fetal bovine serum (FBS) and 15M NaN₃ in 500 mL PBS), they were incubated for 30 min with gentle rotation to attach R-PE-conjugated secondary antibody (#715-116-151, Jackson ImmunoResearch Inc., West Grove, PA, USA). Lastly, the cells were washed twice, and the dead cells were labeled with 7-AAD (#SML1633, Merck, Darmstadt, Germany). The fluorescence of immunolabeled cells was detected and analyzed using BD Accuri™ C6 Plus Flow Cytometer (BD Biosciences, San Jose, CA, USA). After removing 7-AAD signal-positive dead cells, the intensity of the R-PE signal among the GFP-positive cells was measured and plotted.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232012501/s1.

**Author Contributions:** Conceptualization, J.T., J.G. and C.M.; resources and funding acquisition: J.G. and C.M.; formal analysis: T.S., J.P., W.-C.K., X.C., J.T., J.G. and C.M.; writing: J.T., J.G. and C.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** In this work, J.P was supported by the Fondation Roudnitska under the aegis of Fondation de France and C.M was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1A6A1A03040516 and 2020R1A2C2005174). Calculations were done at the Azzurra HPC center of the Observatoire Pluridisciplinaire des Alpes-maritimes (OPAL) research facility.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

### References

1. Debat, H.; Eloit, C.; Blon, F.; Sarazin, B.; Henry, C.; Huet, J.C.; Trotier, D.; Pernollet, J.C. Identification of human olfactory cleft mucus proteins using proteomic analysis. *J. Proteome Res.* **2007**, *6*, 1985–1996. [CrossRef] [PubMed]

2. Nagashima, A.; Touhara, K. Enzymatic conversion of odorants in nasal mucus affects olfactory glomerular activation patterns and odor perception. *J. Neurosci.* **2010**, *30*, 16391–16398. [CrossRef] [PubMed]

3. Pelosi, P.; Knoll, W. Odorant-binding proteins of mammals. *Biol. Rev. Camb. Philos. Soc.* **2022**, *97*, 20–44. [CrossRef] [PubMed]

4. Buck, L.; Axel, R. A novel multigene family may encode odorant receptors: A molecular basis for odor recognition. *Cell* **1991**, *65*, 175–187. [CrossRef]

5. Malnic, B.; Godfrey, P.A.; Buck, L.B. The human olfactory receptor gene family. *Proc. Natl. Acad. Sci. USA* **2004**, *101*, 2584–2589. [CrossRef]
6. Muratore, A.; Clinet, J.C.; Dunach, E. Synthesis of new exo- and endo-3,8-dihydro-beta-santalols and other norbornyl-derived alcohols. *Chem. Biodivers.* **2010**, *7*, 623–638. [CrossRef]

7. Kraft, P.; Frater, G. Enantioselectivity of the musk odor sensation. *Chirality* **2001**, *13*, 388–394. [CrossRef]

8. Trimmer, C.; Keller, A.; Murphy, N.R.; Snyder, L.L.; Willer, J.R.; Nagai, M.H.; Katsanis, N.; Vosshall, L.B.; Matsunami, H.; Mainland, J.D. Genetic variation across the human olfactory receptor repertoire alters odor perception. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 9475–9480. [CrossRef]

9. Gerkin, R.C.; Castro, J.B. The number of olfactory stimuli that humans can discriminate is still unknown. *Elife* **2015**, *4*, e08127. [CrossRef]

10. Bushdid, C.; Magnasco, M.O.; Vosshall, L.B.; Keller, A. Humans can discriminate more than 1 trillion olfactory stimuli. *Science* **2014**, *343*, 1370–1372. [CrossRef]

11. Lawson, Z.; Wheatley, M. The third extracellular loop of G-protein-coupled receptors: More than just a linker between two important transmembrane helices. *Biochem. Soc. Trans.* **2004**, *32*, 1048–1050. [CrossRef]

12. Wheatley, M.; Wootten, D.; Conner, M.T.; Simms, J.; Kendrick, R.; Logan, R.T.; Poyner, D.R.; Barwell, J. Lifting the lid on GPCRs: The role of extracellular loops. *Br. J. Pharmacol.* **2012**, *165*, 1688–1703. [CrossRef] [PubMed]

13. Peeters, M.C.; van Westen, G.J.; Li, Q.; AP, I.J. Importance of the extracellular loops in G protein-coupled receptors for ligand recognition and receptor activation. *Trends Pharmacol. Sci.* **2011**, *32*, 35–42. [CrossRef]

14. Wheatley, M.; Simms, J.; Hawtin, S.R.; Wesley, V.J.; Wootten, D.; Conner, M.; Lawson, Z.; Conner, A.C.; Baker, A.; Cashmore, Y.; et al. Extracellular loops and ligand binding to a subfamily of Family A G-protein-coupled receptors. *Biochem. Soc. Trans.* **2007**, *35*, 717–720. [CrossRef] [PubMed]

15. Dror, R.O.; Arlow, D.H.; Maragakis, P.; Mildorf, T.J.; Pan, A.C.; Xu, H.; Borhani, D.W.; Shaw, D.E. Activation mechanism of the beta2-adrenergic receptor. *Proc. Natl. Acad. Sci. U.S.A.* **2011**, *108*, 18684–18689. [CrossRef] [PubMed]

16. Bushdid, C.; de March, C.A.; Topin, J.; Do, M.; Matsunami, H.; Golebiowski, J. Mammalian class I odorant receptors exhibit a conserved vestibular-binding pocket. *Cell Mol. Life Sci.* **2019**, *76*, 995–1004. [CrossRef]

17. Charlier, L.; Topin, J.; Ronin, C.; Kim, S.K.; Goddard, W.A., 3rd; Efremov, R.; Golebiowski, J. How broadly tuned olfactory receptors equally recognize their agonists. Human OR1G1 as a test case. *Cell Mol. Life Sci.* **2012**, *69*, 4205–4213. [CrossRef]

18. de March, C.A.; Kim, S.K.; Antonczak, S.; Goddard, W.A., 3rd; Golebiowski, J. G protein-coupled odorant receptors: From sequence to structure. *Protein Sci.* **2015**, *24*, 1543–1548. [CrossRef]

19. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Zidek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. *Nature* **2021**, *596*, 583–589. [CrossRef]

20. Jakubik, J.; Randakova, A.; Zimcik, P.; El-Fakahany, E.E.; Dolezal, V. Binding of N-methylscopolamine to the extracellular domain of muscarinic acetylcholine receptors. *Sci. Rep.* **2017**, *7*, 40381. [CrossRef]

21. Geithe, C.; Noe, F.; Kreissl, J.; Krautwurst, D. The Broadly Tuned Odorant Receptor OR1A1 is Highly Selective for 3-Methyl-2,4-nonanedione, a Key Food Odorant in Aged Wines, Tea, and Other Foods. *Chem. Senses* **2017**, *42*, 181–193. [CrossRef] [PubMed]

22. Geithe, C.; Krautwurst, D. Chirality Matters—Enantioselective Orthologous Odorant Receptors for Related Terpenoid Structures. *ACS Sym. Ser.* **2015**, *1212*, 161–181.

23. Geithe, C.; Protze, J.; Kreuchwig, F.; Krause, G.; Krautwurst, D. Structural determinants of a conserved enantiomer-selective carvone binding pocket in the human odorant receptor OR1A1. *Cell Mol. Life Sci.* **2017**, *74*, 4209–4229. [CrossRef] [PubMed]

24. de March, C.A.; Yu, Y.; Ni, M.J.; Adipietro, K.A.; Matsunami, H.; Ma, M.; Golebiowski, J. Conserved Residues Control Activation of Mammalian G Protein-Coupled Odorant Receptors. *J. Am. Chem. Soc.* **2015**, *137*, 8611–8616. [CrossRef]

25. Schmiedeberg, K.; Shirokova, E.; Weber, H.P.; Schilling, B.; Meyerhof, W.; Krautwurst, D. Structural determinants of odorant recognition by the human olfactory receptors OR1A1 and OR1A2. *J. Struct. Biol.* **2007**, *159*, 400–412. [CrossRef] [PubMed]

26. de March, C.A.; Topin, J.; Bruguera, E.; Novikov, G.; Ikegami, K.; Matsunami, H.; Golebiowski, J. Odorant Receptor 7D4 Activation Dynamics. *Angew. Chem. Int. Ed. Engl.* **2018**, *57*, 4554–4558. [CrossRef] [PubMed]

27. Robertson, M.J.; Meyerowitz, J.G.; Panova, O.; Borrelli, K.; Skiniotis, G. Plasticity in ligand recognition at somatostatin receptors. *Nat. Struct. Mol. Biol.* **2022**, *29*, 210–217. [CrossRef]

28. Asada, H.; Horita, S.; Hirata, K.; Shiroishi, M.; Shiimura, Y.; Iwanari, H.; Hamakubo, T.; Shimamura, T.; Nomura, N.; Kusano-Arai, O.; et al. Crystal structure of the human angiotensin II type 2 receptor bound to an angiotensin II analog. *Nat. Struct. Mol. Biol.* **2018**, *25*, 570–576. [CrossRef]

29. Vu, O.; Bender, B.J.; Pankewitz, L.; Huster, D.; Beck-Sickinger, A.G.; Meiler, J. The Structural Basis of Peptide Binding at Class A G Protein-Coupled Receptors. *Molecules* **2021**, *27*, 210. [CrossRef]

30. Webb, B.; Sali, A. Protein structure modeling with MODELLER. *Methods Mol. Biol.* **2014**, *1137*, 1–15. [CrossRef]

31. Cook, B.L.; Steuerwald, D.; Kaiser, L.; Graveland-Bikker, J.; Vanberghem, M.; Berke, A.P.; Herlihy, K.; Pick, H.; Vogel, H.; Zhang, S. Large-scale production and study of a synthetic G protein-coupled receptor: Human olfactory receptor 17-4. *Proc. Natl. Acad. Sci. USA* **2009**, *106*, 11925–11930. [CrossRef] [PubMed]

32. Shen, M.Y.; Sali, A. Statistical potential for assessment and prediction of protein structures. *Protein Sci.* **2006**, *15*, 2507–2524. [CrossRef]
33. Varadi, M.; Anyango, S.; Deshpande, M.; Nair, S.; Natassia, C.; Yordanova, G.; Yuan, D.; Stroe, O.; Wood, G.; Laydon, A.; et al. AlphaFold Protein Structure Database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models. *Nucleic Acids Res.* **2022**, *50*, D439–D444. [CrossRef] [PubMed]

34. Olsson, M.H.; Sondergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. *J. Chem. Theory Comput.* **2011**, *7*, 525–537. [CrossRef] [PubMed]

35. Lomize, M.A.; Pogozheva, I.D.; Joo, H.; Mosberg, H.I.; Lomize, A.L. OPM database and PPM web server: Resources for positioning of proteins in membranes. *Nucleic Acids Res.* **2012**, *40*, D370–D376. [CrossRef] [PubMed]

36. Schott-Verdugo, S.; Gohlke, H. PACKMOL-Memgen: A Simple-To-Use, Generalized Workflow for Membrane-Protein-Lipid-Bilayer System Building. *J. Chem. Inf. Model* **2019**, *59*, 2522–2528. [CrossRef] [PubMed]

37. Lee, T.S.; Cerutti, D.S.; Mermelstein, D.; Lin, C.; LeGrand, S.; Giese, T.J.; Roitberg, A.; Case, D.A.; Walker, R.C.; York, D.M. GPU-Accelerated Molecular Dynamics and Free Energy Methods in Amber18: Performance Enhancements and New Features. *J. Chem. Inf Model* **2018**, *58*, 2043–2050. [CrossRef]

38. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision A.02; Gaussian, Inc.: Wallingford, CT, USA, 2016.

39. Wang, J.M.; Cieplak, P.; Kollman, P.A. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? *J. Comput. Chem.* **2000**, *21*, 1049–1074. [CrossRef]

40. Roe, D.R.; Cheatham, T.E. 3rd. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. *J. Chem. Theory Comput.* **2013**, *9*, 3084–3095. [CrossRef]

41. Shepard, B.D.; Natarajan, N.; Protzko, R.J.; Acres, O.W.; Pluznick, J.L. A cleavable N-terminal signal peptide promotes widespread olfactory receptor surface expression in HEK293T cells. *PLoS ONE* **2013**, *8*, e68758. [CrossRef]

42. Hogrefe, H.H.; Cline, J.; Youngblood, G.L.; Allen, R.M. Creating randomized amino acid libraries with the QuikChange Multi Site-Directed Mutagenesis Kit. *Biotechniques* **2002**, *33*, 1158–1165. [CrossRef] [PubMed]

43. Zhuang, H.; Matsunami, H. Evaluating cell-surface expression and measuring activation of mammalian odorant receptors in heterologous cells. *Nat. Protoc.* **2008**, *3*, 1402–1413. [CrossRef] [PubMed]