Molecular modelling of odorant/olfactory receptor complexes

Modeling of Olfactory Receptors

Olfactory receptors (ORs) form the largest subfamily within class A G protein-coupled receptors (GPCRs). No experimental structural data of any OR is available to date. Homology modeling has become a popular strategy to propose plausible OR models, in order to study the structure-function relationships of the receptors and to aid the discovery and development of ligands capable of modulating receptor activity. In this chapter, we provide a general guideline for OR structure construction, including the collection of candidate templates, structure-based sequence alignment, 3D structure construction, ligand docking, and molecular dynamic simulation.

Hot Spot Mutagenesis Improves the Functional Expression of Unique Mammalian Odorant Receptors

Vertebrate animals detect odors through olfactory receptors (ORs), members of the G protein-coupled receptor (GPCR) family. Due to the difficulty in the heterologous expression of ORs, studies of their odor molecule recognition mechanisms have progressed poorly. Functional expression of most ORs in heterologous cells requires the co-expression of their chaperone proteins, receptor transporting proteins (RTPs). Yet, some ORs were found to be functionally expressed without the support of RTP (RTP-independent ORs). In this study, we investigated whether amino acid residues highly conserved among RTP-independent ORs improve the functional expression of ORs in heterologous cells. We found that a single amino acid substitution at one of two sites (N_BW3.39 and 3.43) in their conserved residues (E and L, respectively) significantly improved the functional expression of ORs in heterologous cells. E^3.39 and L^3.43 also enhanced the membrane expression of RTP-dependent ORs in the absence of RTP. These changes did not alter the odorant responsiveness of the tested ORs. Our results showed that specific sites within transmembrane domains regulate the membrane expression of some ORs.

Bitter Taste and Olfactory Receptors: Beyond Chemical Sensing in the Tongue and the Nose

Abstract

The Up-and-Coming-Scientist section of the current issue of the Journal of Membrane Biology features the invited essay by Dr. Mercedes Alfonso-Prieto, Assistant Professor at the Forschungszentrum Jülich (FZJ), Germany, and the Heinrich-Heine University Düsseldorf, Vogt Institute for Brain Research. Dr. Alfonso-Prieto completed her doctoral degree in chemistry at the Barcelona Science Park, Spain, in 2009, pursued post-doctoral research in computational molecular sciences at Temple University, USA, and then, as a Marie Curie post-doctoral fellow at the University of Barcelona, worked on computations of enzyme reactions and modeling of photoswitchable ligands targeting neuronal receptors. In 2016, she joined the Institute for Advanced Science and the Institute for Computational Biomedicine at the FZJ, where she pursues research on modeling and simulation of chemical senses. The invited essay by Dr. Alfonso-Prieto discusses state-of-the-art modeling of molecular receptors involved in chemical sensing – the senses of taste and smell. These receptors, and computational methods to study them, are the focus of Dr. Alfonso-Prieto’s research. Recently, Dr. Alfonso-Prieto and colleagues have presented a new methodology to predict ligand binding poses for GPCRs, and extensive computations that deciphered the ligand selectivity determinants of bitter taste receptors. These developments inform our current understanding of how taste occurs at the molecular level.

Graphic Abstract

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.

Simulative structure and binding sites of lyral with olfactory receptor 10J5 using computational prediction methods

Olfactory receptor (OR) genes are extensively distributed throughout the human organism. Although these receptors are predominantly located in the olfactory epithelium, binding between odorant chemicals and corresponding ORs initiates downstream events in other tissues. In particular, exposure to allergen fragrances results in the induction of contact dermatitis. At present, current methodologies are limited in their ability to predict the consequences of fragrancy chemicals on humans. The aim of this study was designed to simulate the bindingstructure between lyral and OR10J5, a known allergen which produces contact dermatitis, and its corresponding OR OR10J5 in an effort to predict dermal outcomes using computational methods. Results demonstrated that binding between lyral and OR10J5 involved amino acid residues Phe104, Val105, Cyx178, Ile180, and Tyr258, respectively, which were located on binding sites of the receptor transmembrane 3(TM3), TM3, extracellular loop 2(EL2), EL2, TM6. Evidence indicates that computer simulating binding interactions occurred between an odorant chemical and its receptors which initiated downstream alterations accounting possibly for the observed in vivo contact dermatitis.

Understanding Ligand Binding to G-Protein Coupled Receptors Using Multiscale Simulations

Human G-protein coupled receptors (GPCRs) convey a wide variety of extracellular signals inside the cell and they are one of the main targets for pharmaceutical intervention. Rational drug design requires structural information on these receptors; however, the number of experimental structures is scarce. This gap can be filled by computational models, based on homology modeling and docking techniques. Nonetheless, the low sequence identity across GPCRs and the chemical diversity of their ligands may limit the quality of these models and hence refinement using molecular dynamics simulations is recommended. This is the case for olfactory and bitter taste receptors, which constitute the first and third largest GPCR groups and show sequence identities with the available GPCR templates below 20%. We have developed a molecular dynamics approach, based on the combination of molecular mechanics and coarse grained (MM/CG), tailored to study ligand binding in GPCRs. This approach has been applied so far to bitter taste receptor complexes, showing significant predictive power. The protein/ligand interactions observed in the simulations were consistent with extensive mutagenesis and functional data. Moreover, the simulations predicted several binding residues not previously tested, which were subsequently verified by carrying out additional experiments. Comparison of the simulations of two bitter taste receptors with different ligand selectivity also provided some insights into the binding determinants of bitter taste receptors. Although the MM/CG approach has been applied so far to a limited number of GPCR/ligand complexes, the excellent agreement of the computational models with the mutagenesis and functional data supports the applicability of this method to other GPCRs for which experimental structures are missing. This is particularly important for the challenging case of GPCRs with low sequence identity with available templates, for which molecular docking shows limited predictive power.

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.

Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis

Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.

ORDB, HORDE, ODORactor and other on-line knowledge resources of olfactory receptor-odorant interactions

We present here an exploration of the evolution of three well-established, web-based resources dedicated to the dissemination of information related to olfactory receptors (ORs) and their functional ligands, odorants. These resources are: the Olfactory Receptor Database (ORDB), the Human Olfactory Data Explorer (HORDE) and ODORactor. ORDB is a repository of genomic and proteomic information related to ORs and other chemosensory receptors, such as taste and pheromone receptors. Three companion databases closely integrated with ORDB are OdorDB, ORModelDB and OdorMapDB; these resources are part of the SenseLab suite of databases (http://senselab.med.yale.edu). HORDE (http://genome.weizmann.ac.il/horde/) is a semi-automatically populated database of the OR repertoires of human and several mammals. ODORactor (http://mdl.shsmu.edu.cn/ODORactor/) provides information related to OR-odorant interactions from the perspective of the odorant. All three resources are connected to each other via web-links.

Database URL: http://senselab.med.yale.edu; http://genome.weizmann.ac.il/horde/; http://mdl.shsmu.edu.cn/ODORactor/

Responsiveness of G protein-coupled odorant receptors is partially attributed to the activation mechanism

Mammals detect and discriminate numerous odors via a large family of G protein-coupled odorant receptors (ORs). However, little is known about the molecular and structural basis underlying OR response properties. Using site-directed mutagenesis and computational modeling, we studied ORs sharing high sequence homology but with different response properties. When tested in heterologous cells by diverse odorants, MOR256-3 responded broadly to many odorants, whereas MOR256-8 responded weakly to a few odorants. Out of 36 mutant MOR256-3 ORs, the majority altered the responses to different odorants in a similar manner and the overall response of an OR was positively correlated with its basal activity, an indication of ligand-independent receptor activation. Strikingly, a single mutation in MOR256-8 was sufficient to confer both high basal activity and broad responsiveness to this receptor. These results suggest that broad responsiveness of an OR is at least partially attributed to its activation likelihood.

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.

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.