Prediction of the odorant binding site of olfactory receptor proteins by human-mouse comparisons

Decoding Molecular Mechanism Underlying Human Olfactory Receptor OR8D1 Activation by Sotolone Enantiomers

Sotolone, a chiral compound, plays an important role in the food industry. Herein, (R)-/(S)-sotolone were separated to determine their odor characteristics and thresholds in air (R-form: smoky, burned, herb, and green aroma, 0.0514 μg/m3; S-form: sweet, milk, acid, and nutty aroma, 0.0048 μg/m3). OR8D1 responses to (R)-/(S)-sotolone were detected in a HEK293 cell-based luminescence assay. (S)-Sotolone was a more potent agonist than (R)-sotolone (EC50 values of 84.98 ± 1.05 and 167.20 ± 0.25 μmol/L, respectively). Molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area analyses confirmed that the combination of (S)-sotolone and OR8D1 was more stable than that of (R)-sotolone. Odorant docking, multiple sequence alignments, site-directed mutagenesis, and functional studies with recombinant odorant receptors (ORs) in a cell-based luminescence assay identified 11 amino-acid residues that influence the enantioselectivity of OR8D1 toward sotolone significantly and that N2065.46 was indispensable to the activation of OR8D1 by (S)-sotolone.

Tracking the Diversity and Chromosomal Distribution of the Olfactory Receptor Gene Repertoires of Three Anurans Species

Olfaction is a crucial capability for most vertebrates and is realized through olfactory receptors in the nasal cavity. The enormous diversity of olfactory receptors has been created by gene duplication, following a birth-and-death model of evolution. The olfactory receptor genes of the amphibians have received relatively little attention up to now, although recent studies have increased the number of species for which data are available. This study analyzed the diversity and chromosomal distribution of the OR genes of three anuran species (Engystomops pustulosus, Bufo bufo and Hymenochirus boettgeri). The OR genes were identified through searches for homologies, and sequence filtering and alignment using bioinformatic tools and scripts. A high diversity of OR genes was found in all three species, ranging from 917 in B. bufo to 1194 in H. boettgeri, and a total of 2076 OR genes in E. pustulosus. Six OR groups were recognized using an evolutionary gene tree analysis. While E. pustulosus has one of the highest numbers of genes of the gamma group (which detect airborne odorants) yet recorded in an anuran, B. bufo presented the smallest number of pseudogene sequences ever identified, with no pseudogenes in either the beta or epsilon groups. Although H. boettgeri shares many morphological adaptations for an aquatic lifestyle with Xenopus, and presented a similar number of genes related to the detection of water-soluble odorants, it had comparatively far fewer genes related to the detection of airborne odorants. This study is the first to describe the complete OR repertoire of the three study species and represents an important contribution to the understanding of the evolution and function of the sense of smell in vertebrates.

Modeling the Orthosteric Binding Site of the G Protein-Coupled Odorant Receptor OR5K1

With approximately 400 encoding genes in humans, odorant receptors (ORs) are the largest subfamily of class A G protein-coupled receptors (GPCRs). Despite its high relevance and representation, the odorant-GPCRome is structurally poorly characterized: no experimental structures are available, and the low sequence identity of ORs to experimentally solved GPCRs is a significant challenge for their modeling. Moreover, the receptive range of most ORs is unknown. The odorant receptor OR5K1 was recently and comprehensively characterized in terms of cognate agonists. Here, we report two additional agonists and functional data of the most potent compound on two mutants, L104^3.32 and L255^6.51. Experimental data was used to guide the investigation of the binding modes of OR5K1 ligands into the orthosteric binding site using structural information from AI-driven modeling, as recently released in the AlphaFold Protein Structure Database, and from homology modeling. Induced-fit docking simulations were used to sample the binding site conformational space for ensemble docking. Mutagenesis data guided side chain residue sampling and model selection. We obtained models that could better rationalize the different activity of active (agonist) versus inactive molecules with respect to starting models and also capture differences in activity related to minor structural differences. Therefore, we provide a model refinement protocol that can be applied to model the orthosteric binding site of ORs as well as that of GPCRs with low sequence identity to available templates.

Odor Discrimination by Lipid Membranes

Odor detection and discrimination in mammals is known to be initiated by membrane-bound G-protein-coupled receptors (GPCRs). The role that the lipid membrane may play in odor discrimination, however, is less well understood. Here, we used model membrane systems to test the hypothesis that phospholipid bilayer membranes may be capable of odor discrimination. The effect of S-carvone, R-carvone, and racemic lilial on the model membrane systems was investigated. The odorants were found to affect the fluidity of supported lipid bilayers as measured by fluorescence recovery after photobleaching (FRAP). The effect of odorants on surface-supported lipid multilayer microarrays of different dimensions was also investigated. The lipid multilayer micro- and nanostructure was highly sensitive to exposure to these odorants. Fluorescently-labeled lipid multilayer droplets of 5-micron diameter were more responsive to these odorants than ethanol controls. Arrays of lipid multilayer diffraction gratings distinguished S-carvone from R-carvone in an artificial nose assay. Our results suggest that lipid bilayer membranes may play a role in odorant discrimination and molecular recognition in general.

Advances in the Production of Olfactory Receptors for Industrial Use

In biological olfactory systems, olfactory receptors (ORs) can recognize and discriminate between thousands of volatile organic compounds with very high sensitivity and specificity. The superior properties of ORs have led to the development of OR-based biosensors that have shown promising potential in many applications over the past two decades. In particular, newly designed technologies in gene synthesis, protein expression, solubilization, purification, and membrane mimetics for membrane proteins have greatly opened up the previously inaccessible industrial potential of ORs. In this review, gene design, expression and solubilization strategies, and purification and reconstitution methods available for modern industrial applications are examined, with a focus on ORs. The limitations of current OR production technology are also estimated, and future directions for further progress are suggested.

Large-Scale G Protein-Coupled Olfactory Receptor-Ligand Pairing

G protein-coupled receptors (GPCRs) conserve common structural folds and activation mechanisms, yet their ligand spectra and functions are highly diverse. This work investigated how the amino-acid sequences of olfactory receptors (ORs)-the largest GPCR family-encode diversified responses to various ligands. We established a proteochemometric (PCM) model based on OR sequence similarities and ligand physicochemical features to predict OR responses to odorants using supervised machine learning. The PCM model was constructed with the aid of site-directed mutagenesis, in vitro functional assays, and molecular simulations. We found that the ligand selectivity of the ORs is mostly encoded in the residues up to 8 Å around the orthosteric pocket. Subsequent predictions using Random Forest (RF) showed a hit rate of up to 58%, as assessed by in vitro functional assays of 111 ORs and 7 odorants of distinct scaffolds. Sixty-four new OR-odorant pairs were discovered, and 25 ORs were deorphanized here. The best model demonstrated a 56% deorphanization rate. The PCM-RF approach will accelerate OR-odorant mapping and OR deorphanization.

Two entry tunnels in mouse TAAR9 suggest the possibility of multi-entry tunnels in olfactory receptors

Orthosteric binding sites of olfactory receptors have been well understood for ligand-receptor interactions. However, a lack of explanation for subtle differences in ligand profile of olfactory receptors even with similar orthosteric binding sites promotes more exploration into the entry tunnels of the receptors. An important question regarding entry tunnels is the number of entry tunnels, which was previously believed to be one. Here, we used TAAR9 that recognizes important biogenic amines such as cadaverine, spermine, and spermidine as a model for entry tunnel study. We identified two entry tunnels in TAAR9 and described the residues that form the tunnels. In addition, we found two vestibular binding pockets, each located in one tunnel. We further confirmed the function of two tunnels through site-directed mutagenesis. Our study challenged the existing views regarding the number of entry tunnels in the subfamily of olfactory receptors and demonstrated the possible mechanism how the entry tunnels function in odorant recognition.

An evolutionary conserved olfactory receptor for foodborne and semiochemical alkylpyrazines

Molecular recognition is a fundamental principle in biological systems. The olfactory detection of both food and predators via ecological relevant odorant cues are abilities of eminent evolutionary significance for many species. Pyrazines are such volatile cues, some of which act as both human-centered key food odorants (KFOs) and semiochemicals. A pyrazine-selective odorant receptor has been elusive. Here we screened 2,3,5-trimethylpyrazine, a KFO and semiochemical, and 2,5-dihydro-2,4,5-trimethylthiazoline, an innate fear-associated non-KFO, against 616 human odorant receptor variants, in a cell-based luminescence assay. OR5K1 emerged as sole responding receptor. Tested against a comprehensive collection of 178 KFOs, we newly identified 18 pyrazines and (2R/2S)-4-methoxy-2,5-dimethylfuran-3(2H)-one as agonists. Notably, OR5K1 orthologs in mouse and domesticated species displayed a human-like, potency-ranked activation pattern of pyrazines, suggesting a domestication-led co-evolution of OR5K1 and its orthologs. In summary, OR5K1 is a specialized olfactory receptor across mammals for the detection of pyrazine-based key food odors and semiochemicals.

Functional Architecture of Deleterious Genetic Variants in the Genome of a Wrangel Island Mammoth

Woolly mammoths were among the most abundant cold-adapted species during the Pleistocene. Their once-large populations went extinct in two waves, an end-Pleistocene extinction of continental populations followed by the mid-Holocene extinction of relict populations on St. Paul Island ∼5,600 years ago and Wrangel Island ∼4,000 years ago. Wrangel Island mammoths experienced an episode of rapid demographic decline coincident with their isolation, leading to a small population, reduced genetic diversity, and the fixation of putatively deleterious alleles, but the functional consequences of these processes are unclear. Here, we show that a Wrangel Island mammoth genome had many putative deleterious mutations that are predicted to cause diverse behavioral and developmental defects. Resurrection and functional characterization of several genes from the Wrangel Island mammoth carrying putatively deleterious substitutions identified both loss and gain of function mutations in genes associated with developmental defects (HYLS1), oligozoospermia and reduced male fertility (NKD1), diabetes (NEUROG3), and the ability to detect floral scents (OR5A1). These data suggest that at least one Wrangel Island mammoth may have suffered adverse consequences from reduced population size and isolation.

Citronellal perception and transmission by Anopheles gambiae s.s. (Diptera: Culicidae) females

Anopheles gambiae s.s. is a key vector of Plasmodium parasites. Repellents, which may be a promising alternative to pesticides used to control malaria mosquitoes. Although citronellal is a known mosquito repellent, its repellency characteristics are largely unknown. Determining the specific odorant-binding proteins (OBPs) and odorant receptors (ORs) that detect and transfer the citronellal molecule in A. gambiae s.s. will help to define the mode of action of this compound. In this research, we assessed the repellent activity of citronellal in A. gambiae s.s. using a Y-tube olfactory meter, screened candidate citronellal-binding OBPs and ORs using reverse molecular docking, clarified the binding properties of predicted proteins for citronellal using fluorescence competition binding assay. Results showed that citronellal had a dosage effect on repelling A. gambiae s.s.. The 50% repellent rate was determined to be 4.02 nmol. Results of simulated molecular docking showed that the only proteins that bound tightly with citronellal were AgamOBP4 and AgamORC7. Fluorescence competitive binding assays confirmed the simulations. This research determined that citronellal was captured by AgamOBP4 and transmitted to AgamORC7 in A. gambiae s.s.. Our study will be beneficial in the further understanding the repellent mechanism of citronellal against A. gambiae s.s..

Odorant Receptor Inhibition Is Fundamental to Odor Encoding

Most natural odors are complex mixtures of volatile components, competing to bind odorant receptors (ORs) expressed in olfactory sensory neurons (OSNs) of the nose. To date, surprisingly little is known about how OR antagonism shapes neuronal representations in the detection layer of the olfactory system. Here, we investigated its prevalence, the degree to which it disrupts OR ensemble activity, and its conservation across phylogenetically related ORs. Calcium imaging microscopy of dissociated OSNs revealed significant inhibition, often complete attenuation, of responses to indole-a commonly occurring volatile associated with both floral and fecal odors-by a set of 36 tested odorants. To confirm an OR mechanism for the observed inhibition, we performed single-cell transcriptomics on OSNs exhibiting specific response profiles to a diagnostic panel of odorants and identified three paralogous receptors-Olfr740, Olfr741, and Olfr743-which, when tested in vitro, recapitulated OSN responses. We screened ten ORs from the Olfr740 gene family with ∼800 perfumery-related odorants spanning a range of chemical scaffolds and functional groups. Over half of these compounds (430) antagonized at least one of the ten ORs. OR activity fitted a mathematical model of competitive receptor binding and suggests normalization of OSN ensemble responses to odorant mixtures is the rule rather than the exception. In summary, we observed OR antagonism occurred frequently and in a combinatorial manner. Thus, extensive receptor-mediated computation of mixture information appears to occur in the olfactory epithelium prior to transmission of odor information to the olfactory bulb.

Copper-mediated thiol potentiation and mutagenesis-guided modeling suggest a highly conserved copper-binding motif in human OR2M3

Sulfur-containing compounds within a physiological relevant, natural odor space, such as the key food odorants, typically constitute the group of volatiles with the lowest odor thresholds. The observation that certain metals, such as copper, potentiate the smell of sulfur-containing, metal-coordinating odorants led to the hypothesis that their cognate receptors are metalloproteins. However, experimental evidence is sparse-so far, only one human odorant receptor, OR2T11, and a few mouse receptors, have been reported to be activated by sulfur-containing odorants in a copper-dependent way, while the activation of other receptors by sulfur-containing odorants did not depend on the presence of metals. Here we identified an evolutionary conserved putative copper interaction motif CC/CSSH, comprising two copper-binding sites in TMH5 and TMH6, together with the binding pocket for 3-mercapto-2-methylpentan-1-ol in the narrowly tuned human receptor OR2M3. To characterize the copper-binding motif, we combined homology modeling, docking studies, site-directed mutagenesis, and functional expression of recombinant ORs in a cell-based, real-time luminescence assay. Ligand activation of OR2M3 was potentiated in the presence of copper. This effect of copper was mimicked by ionic and colloidal silver. In two broadly tuned receptors, OR1A1 and OR2W1, which did not reveal a putative copper interaction motif, activation by their most potent, sulfur-containing key food odorants did not depend on the presence of copper. Our results suggest a highly conserved putative copper-binding motif to be necessary for a copper-modulated and thiol-specific function of members from three subfamilies of family 2 ORs.

Concentration-Dependent Recruitment of Mammalian Odorant Receptors

A fundamental challenge in studying principles of organization used by the olfactory system to encode odor concentration information has been identifying comprehensive sets of activated odorant receptors (ORs) across a broad concentration range inside freely behaving animals. In mammals, this has recently become feasible with high-throughput sequencing-based methods that identify populations of activated ORs in vivo In this study, we characterized the mouse OR repertoires activated by the two odorants, acetophenone (ACT) and 2,5-dihydro-2,4,5-trimethylthiazoline (TMT), from 0.01% to 100% (v/v) as starting concentrations using phosphorylated ribosomal protein S6 capture followed by RNA-Seq. We found Olfr923 to be one of the most sensitive ORs that is enriched by ACT. Using a mouse line that genetically labels Olfr923-positive axons, we provided evidence that ACT activates the Olfr923 glomeruli in the olfactory bulb. Through molecular dynamics stimulations, we identified amino acid residues in the Olfr923 binding cavity that facilitate ACT binding. This study sheds light on the active process by which unique OR repertoires may collectively facilitate the discrimination of odorant concentrations.

Signature of adaptive evolution in olfactory receptor genes in Cory's Shearwater supports molecular basis for smell in procellariiform seabirds

Olfactory receptors (ORs), encoded by the largest vertebrate multigene family, enable the detection of thousands of unique odorants in the environment and consequently play a critical role in species survival. Here, we advance our knowledge of OR gene evolution in procellariiform seabirds, an avian group which relies on the sense of olfaction for critical ecological functions. We built a cosmid library of Cory's Shearwater (Calonectris borealis) genomic DNA, a model species for the study of olfaction-based navigation, and sequence OR gene-positive cosmid clones with a combination of sequencing technologies. We identified 220 OR open reading frames, 20 of which are full length, intact OR genes, and found a large ratio of partial and pseudogenes to intact OR genes (2:1), suggestive of a dynamic mode of evolution. Phylogenetic analyses revealed that while a few genes cluster with those of other sauropsid species in a γ (gamma) clade that predates the divergence of different avian lineages, most genes belong to an avian-specific γ-c clade, within which sequences cluster by species, suggesting frequent duplication and/or gene conversion events. We identified evidence of positive selection on full length γ-c clade genes. These patterns are consistent with a key role of adaptation in the functional diversification of olfactory receptor genes in a bird lineage that relies extensively on olfaction.

Convergent degeneration of olfactory receptor gene repertoires in marine mammals

Background

Olfactory receptors (ORs) can bind odor molecules and play a crucial role in odor sensation. Due to the frequent gains and losses of genes during evolution, the number of OR members varies greatly among different species. However, whether the extent of gene gains/losses varies between marine mammals and related terrestrial mammals has not been clarified, and the factors that might underlie these variations are unknown.

Results

To address these questions, we identified more than 10,000 members of the OR family in 23 mammals and classified them into 830 orthologous gene groups (OGGs) and 281 singletons. Significant differences occurred in the number of OR repertoires and OGGs among different species. We found that all marine mammals had fewer OR genes than their related terrestrial lineages, with the fewest OR genes found in cetaceans, which may be closely related to olfactory degradation. ORs with more gene duplications or loss events tended to be under weaker purifying selection. The average gain and loss rates of OR genes in terrestrial mammals were higher than those of mammalian gene families, while the average gain and loss rates of OR genes in marine mammals were significantly lower and much higher than those of mammalian gene families, respectively. Additionally, we failed to detect any one-to-one orthologous genes in the focal species, suggesting that OR genes are not well conserved among marine mammals.

Conclusions

Marine mammals have experienced large numbers of OR gene losses compared with their related terrestrial lineages, which may result from the frequent birth-and-death evolution under varied functional constrains. Due to their independent degeneration, OR genes present in each lineage are not well conserved among marine mammals. Our study provides a basis for future research on the olfactory receptor function in mammals from the perspective of evolutionary trajectories.

A Highly Selective Biosensor Based on Peptide Directly Derived from the HarmOBP7 Aldehyde Binding Site

This paper presents the results of research on determining the optimal length of a peptide chain to effectively bind octanal molecules. Peptides that map the aldehyde binding site in HarmOBP7 were immobilized on piezoelectric transducers. Based on computational studies, four Odorant Binding Protein-derived Peptides (OBPPs) with different sequences were selected. Molecular modelling results of ligand docking with selected peptides were correlated with experimental results. The use of low-molecular synthetic peptides, instead of the whole protein, enabled the construction OBPPs-based biosensors. This work aims at developing a biomimetic piezoelectric OBPPs sensor for selective detection of octanal. Moreover, the research is concerned with the ligand binding affinity depending on different peptides’ chain lengths. The authors believe that the chain length can have a substantial influence on the type and effectiveness of peptide–ligand interaction. A confirmation of in silico investigation results is the correlation with the experimental results, which shows that the highest affinity to octanal is exhibited by the longest peptide (OBPP4 – KLLFDSLTDLKKKMSEC-NH₂). We hypothesized that the binding of long chain aldehydes to the peptide, mimicking the binding site of HarmOBP7, induced a conformational change in the peptide deposited on a selected transducer. The constructed OBPP4-based biosensors were able to selectively bind octanal in the gas phase. It was also shown that the sensors were characterized by high selectivity with respect to octanal, as well as to acetaldehyde and benzaldehyde. The results indicate that the OBPP4 peptide, mimicking the binding domain in the Odorant Binding Protein, can provide new opportunities for the development of biomimicking materials in the field of odor biosensors.

Is It Possible to Predict the Odor of a Molecule on the Basis of its Structure?

The olfactory sense is the dominant sensory perception for many animals. When Richard Axel and Linda B. Buck received the Nobel Prize in 2004 for discovering the G protein-coupled receptors’ role in olfactory cells, they highlighted the importance of olfaction to the scientific community. Several theories have tried to explain how cells are able to distinguish such a wide variety of odorant molecules in a complex context in which enantiomers can result in completely different perceptions and structurally different molecules. Moreover, sex, age, cultural origin, and individual differences contribute to odor perception variations that complicate the picture. In this article, recent advances in olfaction theory are presented, and future trends in human olfaction such as structure-based odor prediction and artificial sniffing are discussed at the frontiers of chemistry, physiology, neurobiology, and machine learning.

Carbon chain shape selectivity by the mouse olfactory receptor OR-I7

The rodent OR-I7 is an olfactory receptor exemplar activated by aliphatic aldehydes such as octanal. Normal alkanals shorter than heptanal bind OR-I7 without activating it and hence function as antagonists in vitro. We report a series of aldehydes designed to probe the structural requirements for aliphatic ligand chains too short to meet the minimum approximate 6.9 Å length requirement for receptor activation. Experiments using recombinant mouse OR-I7 expressed in heterologous cells show that in the context of short aldehyde antagonists, OR-I7 prefers binding aliphatic chains without branches, though a single methyl on carbon-3 is permitted. The receptor can accommodate a surprisingly large number of carbons (e.g. ten in adamantyl) as long as the carbons are part of a conformationally constrained ring system. A rhodopsin-based homology model of mouse OR-I7 docked with the new antagonists suggests that small alkyl branches on the alkyl chain sterically interfere with the hydrophobic residues lining the binding site, but branch carbons can be accommodated when tied back into a compact ring system like the adamantyl and bicyclo[2.2.2]octyl systems.

Olfactory receptor 43 reduces hepatic lipid accumulation and adiposity in mice

Olfactory receptors are primarily expressed in nasal olfactory epithelium, but these receptors are also ectopically expressed in diverse tissues. In this study, we investigated the biological functions of Olfr43, a mouse homolog of human OR1A1, in cultured hepatocytes and mice to assess its functionality in lipid metabolism. Olfr43 was expressed in mouse hepatocytes, and Olfr43 activation by a known ligand, (-)-carvone, stimulated cAMP response element-binding protein (CREB) activity. In ligand-receptor binding studies using site-directed mutagenesis, (-)-carvone binding required two residues, M257 and Y258, in Olfr43. In the mouse study, oral administration of (-)-carvone for 5 weeks in high-fat diet-fed mice improved energy metabolism, including reductions in hepatic steatosis and adiposity, and improved glucose and insulin tolerance. In mouse livers and cultured mouse hepatocytes, Olfr43 activation simulated the CREB-hairy and enhancer of split 1 (HES1)-peroxisome proliferator-activated receptor (PPAR)-γ signaling axis, leading to a reduction in hepatic triglyceride accumulation in the mouse liver. Thus, long-term administration of (-)-carvone reduces hepatic steatosis. The knockdown of Olfr43 gene expression in cultured hepatocytes negated these effects of (-)-carvone. In conclusion, an ectopic olfactory receptor, hepatic Olfr43, regulates energy metabolism via the CREB-HES1-PPARγ signaling axis.

Agonists of G-Protein-Coupled Odorant Receptors Are Predicted from Chemical Features

Predicting the activity of chemicals for a given odorant receptor is a longstanding challenge. Here the activity of 258 chemicals on the human G-protein-coupled odorant receptor (OR)51E1, also known as prostate-specific G-protein-coupled receptor 2 (PSGR2), was virtually screened by machine learning using 4884 chemical descriptors as input. A systematic control by functional in vitro assays revealed that a support vector machine algorithm accurately predicted the activity of a screened library. It allowed us to identify two novel agonists in vitro for OR51E1. The transferability of the protocol was assessed on OR1A1, OR2W1, and MOR256-3 odorant receptors, and, in each case, novel agonists were identified with a hit rate of 39-50%. We further show how ligands' efficacy is encoded into residues within OR51E1 cavity using a molecular modeling protocol. Our approach allows widening the chemical spaces associated with odorant receptors. This machine-learning protocol based on chemical features thus represents an efficient tool for screening ligands for G-protein-coupled odorant receptors that modulate non-olfactory functions or, upon combinatorial activation, give rise to our sense of smell.

Dynamical Binding Modes Determine Agonistic and Antagonistic Ligand Effects in the Prostate-Specific G-Protein Coupled Receptor (PSGR)

We analysed the ligand-based activation mechanism of the prostate-specific G-protein coupled receptor (PSGR), which is an olfactory receptor that mediates cellular growth in prostate cancer cells. Furthermore, it is an olfactory receptor with a known chemically near identic antagonist/agonist pair, α- and β-ionone. Using a combined theoretical and experimental approach, we propose that this receptor is activated by a ligand-induced rearrangement of a protein-internal hydrogen bond network. Surprisingly, this rearrangement is not induced by interaction of the ligand with the network, but by dynamic van der Waals contacts of the ligand with the involved amino acid side chains, altering their conformations and intraprotein connectivity. Ligand recognition in this GPCR is therefore highly stereo selective, but seemingly lacks any ligand recognition via polar contacts. A putative olfactory receptor-based drug design scheme will have to take this unique mode of protein/ligand action into account.

Structural determinants of a conserved enantiomer-selective carvone binding pocket in the human odorant receptor OR1A1

Chirality is a common phenomenon within odorants. Most pairs of enantiomers show only moderate differences in odor quality. One example for enantiomers that are easily discriminated by their odor quality is the carvones: humans significantly distinguish between the spearmint-like (R)-(-)-carvone and caraway-like (S)-(+)-carvone enantiomers. Moreover, for the (R)-(-)-carvone, an anosmia is observed in about 8% of the population, suggesting enantioselective odorant receptors (ORs). With only about 15% de-orphaned human ORs, the lack of OR crystal structures, and few comprehensive studies combining in silico and experimental approaches to elucidate structure-function relations of ORs, knowledge on cognate odorant/OR interactions is still sparse. An adjusted homology modeling approach considering OR-specific proline-caused conformations, odorant docking studies, single-nucleotide polymorphism (SNP) analysis, site-directed mutagenesis, and subsequent functional studies with recombinant ORs in a cell-based, real-time luminescence assay revealed 11 amino acid positions to constitute an enantioselective binding pocket necessary for a carvone function in human OR1A1 and murine Olfr43, respectively. Here, we identified enantioselective molecular determinants in both ORs that discriminate between minty and caraway odor. Comparison with orthologs from 36 mammalian species demonstrated a hominid-specific carvone binding pocket with about 100% conservation. Moreover, we identified loss-of-function SNPs associated with the carvone binding pocket of OR1A1. Given carvone enantiomer-specific receptor activation patterns including OR1A1, our data suggest OR1A1 as a candidate receptor for constituting a carvone enantioselective phenotype, which may help to explain mechanisms underlying a (R)-(-)-carvone-specific anosmia in humans.

Twenty years of ModelDB and beyond: building essential modeling tools for the future of neuroscience

Neuron modeling may be said to have originated with the Hodgkin and Huxley action potential model in 1952 and Rall's models of integrative activity of dendrites in 1964. Over the ensuing decades, these approaches have led to a massive development of increasingly accurate and complex data-based models of neurons and neuronal circuits. ModelDB was founded in 1996 to support this new field and enhance the scientific credibility and utility of computational neuroscience models by providing a convenient venue for sharing them. It has grown to include over 1100 published models covering more than 130 research topics. It is actively curated and developed to help researchers discover and understand models of interest. ModelDB also provides mechanisms to assist running models both locally and remotely, and has a graphical tool that enables users to explore the anatomical and biophysical properties that are represented in a model. Each of its capabilities is undergoing continued refinement and improvement in response to user experience. Large research groups (Allen Brain Institute, EU Human Brain Project, etc.) are emerging that collect data across multiple scales and integrate that data into many complex models, presenting new challenges of scale. We end by predicting a future for neuroscience increasingly fueled by new technology and high performance computation, and increasingly in need of comprehensive user-friendly databases such as ModelDB to provide the means to integrate the data for deeper insights into brain function in health and disease.

Genomic organization and evolution of olfactory receptors and trace amine-associated receptors in channel catfish, Ictalurus punctatus

BACKGROUND

Channel catfish (Ictalurus punctatus) live in turbid waters with limited visibility to chase prey within a certain distance. This can be compensated through detecting specific water-soluble substances by the olfactory receptors (ORs) and trace amine associated receptors (TAARs) expressed on the olfactory epithelium.

METHODS

We identified the OR and TAAR repertoires in channel catfish, and characterized the genomic organizations of these two gene families by data mining available genomic resources.

RESULTS

A total of 47 putative OR genes and 36 putative TAAR genes were identified in the channel catfish genome, including 27 functional OR genes and 28 functional TAAR genes. Phylogenetic and orthogroup analyses were conducted to illustrate the evolutionary dynamics of the vertebrate ORs and TAARs. Collinear analysis revealed the presence of two conserved orthologous blocks that contain OR genes between the catfish genome and zebrafish genome. The complete loss of a conserved motif in fish OR family H may contribute to the divergence of family H from other families. The dN/dS analysis indicated that the highest degree of selection pressure was imposed on TAAR subfamily 14 among all fish ORs and TAARs.

CONCLUSIONS

The present study provides understanding of the evolutionary dynamics of the two gene families (OR and TAAR) associated with olfaction in channel catfish.

GENERAL SIGNIFICANCE

This is the first systematic study of ORs and TAARs in catfish, which could provide valuable genomic resources for further investigation of olfactory mechanisms in teleost fish.

A Family of non-GPCR Chemosensors Defines an Alternative Logic for Mammalian Olfaction

Odor perception in mammals is mediated by parallel sensory pathways that convey distinct information about the olfactory world. Multiple olfactory subsystems express characteristic seven-transmembrane G-protein-coupled receptors (GPCRs) in a one-receptor-per-neuron pattern that facilitates odor discrimination. Sensory neurons of the "necklace" subsystem are nestled within the recesses of the olfactory epithelium and detect diverse odorants; however, they do not express known GPCR odor receptors. Here, we report that members of the four-pass transmembrane MS4A protein family are chemosensors expressed within necklace sensory neurons. These receptors localize to sensory endings and confer responses to ethologically relevant ligands, including pheromones and fatty acids, in vitro and in vivo. Individual necklace neurons co-express many MS4A proteins and are activated by multiple MS4A ligands; this pooling of information suggests that the necklace is organized more like subsystems for taste than for smell. The MS4As therefore define a distinct mechanism and functional logic for mammalian olfaction.

Odorant receptors of Drosophila are sensitive to the molecular volume of odorants

Which properties of a molecule define its odor? This is a basic yet unanswered question regarding the olfactory system. The olfactory system of Drosophila has a repertoire of approximately 60 odorant receptors. Molecules bind to odorant receptors with different affinities and activate them with different efficacies, thus providing a combinatorial code that identifies odorants. We hypothesized that the binding affinity of an odorant-receptor pair is affected by their relative sizes. The maximum affinity can be attained when the molecular volume of an odorant matches the volume of the binding pocket. The affinity drops to zero when the sizes are too different, thus obscuring the effects of other molecular properties. We developed a mathematical formulation of this hypothesis and verified it using Drosophila data. We also predicted the volume and structural flexibility of the binding site of each odorant receptor; these features significantly differ between odorant receptors. The differences in the volumes and structural flexibilities of different odorant receptor binding sites may explain the difference in the scents of similar molecules with different sizes.

In Silico Modeling of Novel Drug Ligands for Treatment of Concussion Associated Tauopathy

The objective of this study was to develop an in silico screening model for characterization of potential novel ligands from commercial drug libraries able to functionally activate certain olfactory receptors (ORs), which are members of the class A rhodopsin-like family of G protein couple receptors (GPCRs), in the brain of murine models of concussion. We previously found that concussions may significantly influence expression of certain ORs, for example, OR4M1 in subjects with a history of concussion/traumatic brain injury (TBI). In this study, we built a 3-D OR4M1 model and used it in in silico screening of potential novel ligands from commercial drug libraries. We report that in vitro activation of OR4M1 with the commercially available ZINC library compound 10915775 led to a significant attenuation of abnormal tau phosphorylation in embryonic cortico-hippocampal neuronal cultures derived from NSE-OR4M1 transgenic mice, possibly through modulation of the JNK signaling pathway. The attenuation of abnormal tau phosphorylation was rather selective since ZINC10915775 significantly decreased tau phosphorylation on tau Ser202/T205 (AT8 epitope) and tau Thr212/Ser214 (AT100 epitope), but not on tau Ser396/404 (PHF-1 epitope). Moreover, no response of ZINC10915775 was found in control hippocampal neuronal cultures derived from wild type littermates. Our in silico model provides novel means to pharmacologically modulate select ubiquitously expressed ORs in the brain through high affinity ligand activation to prevent and eventually to treat concussion induced down regulation of ORs and subsequent cascade of tau pathology. J. Cell. Biochem. 117: 2241-2248, 2016. © 2016 Wiley Periodicals, Inc.

Exchanging ligand-binding specificity between a pair of mouse olfactory receptor paralogs reveals odorant recognition principles

A multi-gene family of ~1000 G protein-coupled olfactory receptors (ORs) constitutes the molecular basis of mammalian olfaction. Due to the lack of structural data its remarkable capacity to detect and discriminate thousands of odorants remains poorly understood on the structural level of the receptor. Using site-directed mutagenesis we transferred ligand specificity between two functionally related ORs and thereby revealed amino acid residues of central importance for odorant recognition and discrimination of the two receptors. By exchanging two of three residues, differing at equivalent positions of the putative odorant binding site between the mouse OR paralogs Olfr73 (mOR-EG) and Olfr74 (mOR-EV), we selectively changed ligand preference but remarkably also signaling activation strength in both ORs. Computer modeling proposed structural details at atomic resolution how the very same odorant molecule might interact with different contact residues to induce different functional responses in two related receptors. Our findings provide a mechanistic explanation of how the olfactory system distinguishes different molecular aspects of a given odorant molecule, and unravel important molecular details of the combinatorial encoding of odorant identity at the OR level.

G protein-coupled odorant receptors: From sequence to structure

Odorant receptors (ORs) are the largest subfamily within class A G protein-coupled receptors (GPCRs). No experimental structural data of any OR is available to date and atomic-level insights are likely to be obtained by means of molecular modeling. In this article, we critically align sequences of ORs with those GPCRs for which a structure is available. Here, an alignment consistent with available site-directed mutagenesis data on various ORs is proposed. Using this alignment, the choice of the template is deemed rather minor for identifying residues that constitute the wall of the binding cavity or those involved in G protein recognition.

Computational Approaches for Decoding Select Odorant-Olfactory Receptor Interactions Using Mini-Virtual Screening

Olfactory receptors (ORs) belong to the class A G-Protein Coupled Receptor superfamily of proteins. Unlike G-Protein Coupled Receptors, ORs exhibit a combinatorial response to odors/ligands. ORs display an affinity towards a range of odor molecules rather than binding to a specific set of ligands and conversely a single odorant molecule may bind to a number of olfactory receptors with varying affinities. The diversity in odor recognition is linked to the highly variable transmembrane domains of these receptors. The purpose of this study is to decode the odor-olfactory receptor interactions using in silico docking studies. In this study, a ligand (odor molecules) dataset of 125 molecules was used to carry out in silico docking using the GLIDE docking tool (SCHRODINGER Inc Pvt LTD). Previous studies, with smaller datasets of ligands, have shown that orthologous olfactory receptors respond to similarly-tuned ligands, but are dramatically different in their efficacy and potency. Ligand docking results were applied on homologous pairs (with varying sequence identity) of ORs from human and mouse genomes and ligand binding residues and the ligand profile differed among such related olfactory receptor sequences. This study revealed that homologous sequences with high sequence identity need not bind to the same/ similar ligand with a given affinity. A ligand profile has been obtained for each of the 20 receptors in this analysis which will be useful for expression and mutation studies on these receptors.

Conserved Residues Control Activation of Mammalian G Protein-Coupled Odorant Receptors

Odorant receptor (OR) genes and proteins represent more than 2% of our genome and 4% of our proteome and constitute the largest subgroup of G protein-coupled receptors (GPCRs). The mechanism underlying OR activation remains poorly understood, as they do not share some of the highly conserved motifs critical for activation of non-olfactory GPCRs. By combining site-directed mutagenesis, heterologous expression, and molecular dynamics simulations that capture the conformational change of constitutively active mutants, we tentatively identified crucial residues for the function of these receptors using the mouse MOR256-3 (Olfr124) as a model. The toggle switch for sensing agonists involves a highly conserved tyrosine residue in helix VI. The ionic lock is located between the "DRY" motif in helix III and a positively charged "R/K" residue in helix VI. This study provides an unprecedented model that captures the main mechanisms of odorant receptor activation.

The molecular biology of vertebrate olfaction

The importance of chemosensation for vertebrates is reflected in the vast and variable nature of their chemosensory tissues, neurons, and genes, which we explore in this review. Immense progress has been made in elucidating the molecular biology of olfaction since the discovery of the olfactory receptor genes by Buck and Axel, which eventually won the authors the Nobel Prize. In particular, research linking odor ligands to olfactory receptors (ORs) is truly revolutionizing our understanding of how a large but limited number of chemosensory receptors can allow us to perceive the massive diversity of odors in our habitat. This research is providing insight into the evolution of genomes and providing the raw data needed to explore links between genotype and phenotype, still a grand challenge in biology. Research into olfaction is still developing and will no doubt continue until we have a clear understanding of how all odors are detected and the evolutionary forces that have molded the chemosensory subgenome in vertebrates. This knowledge will not only be a huge step in elucidating olfactory function, advancing scientific knowledge and techniques, but there are also commercial applications for this research. This review focuses on the molecular basis of chemosensation, particularly olfaction, its evolution across vertebrates and the recent molecular advances linking odors to their cognate receptors.

Molecular recognition of ketamine by a subset of olfactory G protein-coupled receptors

Ketamine elicits various neuropharmacological effects, including sedation, analgesia, general anesthesia, and antidepressant activity. Through an in vitro screen, we identified four mouse olfactory receptors (ORs) that responded to ketamine. In addition to their presence in the olfactory epithelium, these G protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptors (GPCRs) are distributed throughout the central nervous system. To better understand the molecular basis of the interactions between ketamine and ORs, we used sequence comparison and molecular modeling to design mutations that (i) increased, reduced, or abolished ketamine responsiveness in responding receptors, and (ii) rendered nonresponding receptors responsive to ketamine. We showed that olfactory sensory neurons (OSNs) that expressed distinct ORs responded to ketamine in vivo, suggesting that ORs may serve as functional targets for ketamine. The ability to both abolish and introduce responsiveness to ketamine in GPCRs enabled us to identify and confirm distinct interaction loci in the binding site, which suggested a signature ketamine-binding pocket that may guide exploration of additional receptors for this general anesthetic drug.

Receptive range analysis of a mouse odorant receptor subfamily

Mammals deploy a large array of odorant receptors (ORs) to detect and distinguish a vast number of odorant molecules. ORs vary widely in the type of odorant structures recognized and in the breadth of molecular receptive range (MRR), with some ORs recognizing a small group of closely related molecules and other ORs recognizing a wide range of structures. While closely related ORs have been shown to have similar MRRs, the functional relationships among less closely related ORs are unclear. We screened a small group of ORs with a diverse odorant panel to identify a new odorant-OR pairing (unsaturated aldehydes and MOR263-3). We then extensively screened MOR263-3 and a series of additional MORs related to MOR263-3 in various ways. MORs related by phylogenetic analysis (several other members of the MOR263 subfamily) had MRRs that overlapped with the MRR of MOR263-3, even with amino acid identity as low as 48% (MOR263-2). MOR171-17, predicted to be functionally related to MOR263-3 by an alternative bioinformatic analysis, but with only 39% amino acid identity, had a distinct odorant specificity. Our results support the use of phylogenetic analysis to predict functional relationships among ORs with relatively low amino acid identity. We screened a small group of mouse odorant receptors (MORs) with a diverse odorant panel to identify a new odorant-OR pairing (unsaturated aldehydes and MOR263-3), then extensively screened a series of additional MORs related to MOR263-3 in various ways. MORs related by phylogenetic analysis had odorant specificities that overlapped with that of MOR263-3, but MOR171-17, predicted to be functionally related to MOR263-3 by an alternative bioinformatic analysis, had a distinct odorant specificity.

Mammalian olfactory receptors: molecular mechanisms of odorant detection, 3D-modeling, and structure-activity relationships

This chapter describes the main characteristics of olfactory receptor (OR) genes of vertebrates, including generation of this large multigenic family and pseudogenization. OR genes are compared in relation to evolution and among species. OR gene structure and selection of a given gene for expression in an olfactory sensory neuron (OSN) are tackled. The specificities of OR proteins, their expression, and their function are presented. The expression of OR proteins in locations other than the nasal cavity is regulated by different mechanisms, and ORs display various additional functions. A conventional olfactory signal transduction cascade is observed in OSNs, but individual ORs can also mediate different signaling pathways, through the involvement of other molecular partners and depending on the odorant ligand encountered. ORs are engaged in constitutive dimers. Ligand binding induces conformational changes in the ORs that regulate their level of activity depending on odorant dose. When present, odorant binding proteins induce an allosteric modulation of OR activity. Since no 3D structure of an OR has been yet resolved, modeling has to be performed using the closest G-protein-coupled receptor 3D structures available, to facilitate virtual ligand screening using the models. The study of odorant binding modes and affinities may infer best-bet OR ligands, to be subsequently checked experimentally. The relationship between spatial and steric features of odorants and their activity in terms of perceived odor quality are also fields of research that development of computing tools may enhance.

Discrimination between olfactory receptor agonists and non-agonists

A joint approach combining free-energy calculations and calcium-imaging assays on the broadly tuned human 1G1 olfactory receptor is reported. The free energy of binding of ten odorants was computed by means of molecular-dynamics simulations. This state function allows separating the experimentally determined eight agonists from the two non-agonists. This study constitutes a proof-of-principle for the computational deorphanization of olfactory receptors.

A computational microscope focused on the sense of smell

In this article, we review studies of the protagonists of the perception of smell focusing on Odorant-Binding Proteins and Olfactory Receptors. We notably put forward studies performed by means of molecular modeling, generally combined with experimental data. Those works clearly emphasize that computational approaches are now a force to reckon with. In the future, they will certainly be more and more used, notably in the framework of a computational microscope meant to observe how the laws of physics govern the biomolecular systems originating our sense of smell.

The state of the art of odorant receptor deorphanization: a report from the orphanage

The odorant receptors (ORs) provide our main gateway to sensing the world of volatile chemicals. This involves a complex encoding process in which multiple ORs, each of which detects its own set of odorants, work as an ensemble to produce a distributed activation code that is presumably unique to each odorant. One marked challenge to decoding the olfactory code is OR deorphanization, the identification of a set of activating odorants for a particular receptor. Here, we survey various methods used to try to express defined ORs of interest. We also suggest strategies for selecting odorants for test panels to evaluate the functional expression of an OR. Integrating these tools, while retaining awareness of their idiosyncratic limitations, can provide a multi-tiered approach to OR deorphanization, spanning the initial discovery of a ligand to vetting that ligand in a physiologically relevant setting.

Genetic signatures for enhanced olfaction in the African mole-rats

The Olfactory Receptor (OR) superfamily, the largest in the vertebrate genome, is responsible for vertebrate olfaction and is traditionally subdivided into 17 OR families. Recent studies characterising whole-OR subgenomes revealed a 'birth and death' model of evolution for a range of species, however little is known about fine-scale evolutionary dynamics within single-OR families. This study reports the first assessment of fine-scale OR evolution and variation in African mole-rats (Bathyergidae), a family of subterranean rodents endemic to sub-Saharan Africa. Because of the selective pressures of life underground, enhanced olfaction is proposed to be fundamental to the evolutionary success of the Bathyergidae, resulting in a highly diversified OR gene-repertoire. Using a PCR-sequencing approach, we analysed variation in the OR7 family across 14 extant bathyergid species, which revealed enhanced levels of functional polymorphisms concentrated across the receptors' ligand-binding region. We propose that mole-rats are able to recognise a broad range of odorants and that this diversity is reflected throughout their OR7 gene repertoire. Using both classic tests and tree-based methods to test for signals of selection, we investigate evolutionary forces across the mole-rat OR7 gene tree. Four well-supported clades emerged in the OR phylogeny, with varying signals of selection; from neutrality to positive and purifying selection. Bathyergid life-history traits and environmental niche-specialisation are explored as possible drivers of adaptive OR evolution, emerging as non-exclusive contributors to the positive selection observed at OR7 genes. Our results reveal unexpected complexity of evolutionary mechanisms acting within a single OR family, providing insightful perspectives into OR evolutionary dynamics.

Preferential binding of an odor within olfactory receptors: a precursor to receptor activation

Using computational methods, which allow mechanistic insights at a molecular level, we explored the olfactory receptor (OR)-odor interactions for 2 mouse ORs, S79 and S86. Both ORs have been previously experimentally, functionally characterized. The odors used were mostly carboxylic acids, which differed in chain length, substituents on the primary carbon atom-chain and degree of unsaturation. These odors elicited varied activation responses from both ORs. Our studies revealed that both receptors have 2 distinct binding sites. Preferential binding in 1 of the 2 sites is correlated with OR activation. The activating odorants: nonanedioic acid, heptanoic acid, and octanoic acid for OR S79 and nonanoic acid for OR S86 preferentially bind in the region bound by transmembranes (TMs [helical domains]) III, IV, V, and VI. The non excitatory odorants heptanol for S79 and heptanoic acid for S86 showed a greater likelihood of binding in the region bound by TMs I, II, III, and VII. Nanosecond-scale molecular dynamics simulations of the physiologically relevant conditions of docked OR-odorant complexes enabled us to quantitatively assess the roles of individual OR amino acids in odor binding. Amino acid-odorant contact maps and distance determinations over the course of the simulations lend support to our conclusions.

Molecular modelling of odorant/olfactory receptor complexes

Providing a rationale that associates a chemical structure of an odorant to its induced perception has been sought for a long time. To achieve this, a detailed atomic structure of both the odorant and the olfactory receptor must be known. State-of-the-art techniques to model the 3D structure of an olfactory receptor in complex with various odorants are presented here. These range from sequence alignment with known structures to molecular dynamics simulations in a realistic environment.

HORDE: comprehensive resource for olfactory receptor genomics

Olfactory receptors (ORs) constitute the largest gene family in the mammalian genome. The existence of these proteins underlies the nature of, and variability in, odorant perception. The Human Olfactory Receptor Data Explorer (HORDE, http://genome.weizmann.ac.il/horde/ ) is a free online resource, which presents a complete compendium of all OR genes and pseudogenes in the genome of human and four other vertebrates. HORDE includes three parts: (1) an automated pipeline, which mines OR gene and pseudogene sequences out of complete genomes, and generates gene symbols based on sequence similarity; (2) a card generator that obtains and displays annotative information on individual ORs retrieved from external databases and relevant studies; and (3) a search engine that allows user retrieval of OR information. For human ORs, HORDE specifically addresses the universe of interindividual variation, as obtained from several sources, including whole genome sequences made possible by next-generation sequencing. This encompasses single nucleotide polymorphisms (SNP) and copy number variation (CNV), including deleterious mutational events. HORDE also hosts a number of tools designed specifically to assist in the study of OR evolution and function. In this chapter, we describe the status of HORDE (build #43). We also discuss plans for future enhancements and a road map for HORDE to become a better community-based bioinformatics tool. We highlight HORDE's role as a major research tool in the study of an expanding cohort of OR repertoires.