Modeling of mammalian olfactory receptors and docking of odorants

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.

Binding site identification of G protein-coupled receptors through a 3D Zernike polynomials-based method: application to C. elegans olfactory receptors

Studying the binding processes of G protein-coupled receptors (GPCRs) proteins is of particular interest both to better understand the molecular mechanisms that regulate the signaling between the extracellular and intracellular environment and for drug design purposes. In this study, we propose a new computational approach for the identification of the binding site for a specific ligand on a GPCR. The method is based on the Zernike polynomials and performs the ligand-GPCR association through a shape complementarity analysis of the local molecular surfaces. The method is parameter-free and it can distinguish, working on hundreds of experimentally GPCR-ligand complexes, binding pockets from randomly sampled regions on the receptor surface, obtaining an Area Under ROC curve of 0.77. Given its importance both as a model organism and in terms of applications, we thus investigated the olfactory receptors of the C. elegans, building a list of associations between 21 GPCRs belonging to its olfactory neurons and a set of possible ligands. Thus, we can not only carry out rapid and efficient screenings of drugs proposed for GPCRs, key targets in many pathologies, but also we laid the groundwork for computational mutagenesis processes, aimed at increasing or decreasing the binding affinity between ligands and receptors.

A two-stage computational approach to predict novel ligands for a chemosensory receptor

Olfactory receptor (OR) 1A2 is the member of largest superfamily of G protein-coupled receptors (GPCRs). OR1A2 is an ectopically expressed receptor with only 13 known ligands, implicated in reducing hepatocellular carcinoma progression, with enormous therapeutic potential. We have developed a two-stage screening approach to identify novel putative ligands of OR1A2. We first used a pharmacophore model based on atomic property field (APF) to virtually screen a library of 5942 human metabolites. We then carried out structure-based virtual screening (SBVS) for predicting the potential agonists, based on a 3D homology model of OR1A2. This model was developed using a biophysical approach for template selection, based on multiple parameters including hydrophobicity correspondence, applied to the complete set of available GPCR structures to pick the most appropriate template. Finally, the membrane-embedded 3D model was refined by molecular dynamics (MD) simulations in both the apo and holo forms. The refined model in the apo form was selected for SBVS. Four novel small molecules were identified as strong binders to this olfactory receptor on the basis of computed binding energies.

Therapeutic potential of ectopic olfactory and taste receptors

Olfactory and taste receptors are expressed primarily in the nasal olfactory epithelium and gustatory taste bud cells, where they transmit real-time sensory signals to the brain. However, they are also expressed in multiple extra-nasal and extra-oral tissues, being implicated in diverse biological processes including sperm chemotaxis, muscle regeneration, bronchoconstriction and bronchodilatation, inflammation, appetite regulation and energy metabolism. Elucidation of the physiological roles of these ectopic receptors is revealing potential therapeutic and diagnostic applications in conditions including wounds, hair loss, asthma, obesity and cancers. This Review outlines current understanding of the diverse functions of ectopic olfactory and taste receptors and assesses their potential to be therapeutically exploited.

Structure-function relationship between a natural cosmetic active ingredient and the olfactory receptor OR2AT4

OBJECTIVE

Although the olfactory receptor OR2AT4 was described as involved in epidermal renewal, there is no data about a cosmetic active ingredient activating this receptor. The aim of this research work was thus to identify a natural molecule binding to this receptor in order to stimulate keratinocyte migration.

METHODS

For this purpose, natural molecules were extracted from Cocos nucifera flour. Then, efficacy of this natural extract was evaluated on keratinocyte migration in vitro. Molecules of the Cocos nucifera flour extract were then identified by UPLC-MS/MS. Molecular docking was finally conducted to investigate the potential interaction between identified molecules and the olfactory receptor OR2AT4.

RESULTS

The Cocos nucifera flour extract significantly increased keratinocyte migration and results demonstrated that this effect was mediated by the olfactory receptor OR2AT4. Metabolomic analysis revealed two molecules, nonioside D and butyl 4-O-alpha-D-glucopyranosyl-beta-D-glucopyranoside, as significantly present in the Cocos nucifera flour extract compared to both Cocos nucifera oil and water. Finally, molecular docking revealed that butyl 4-O-alpha-D-glucopyranosyl-beta-D-glucopyranoside could interact with the extracellular domain 2 of the OR2AT4.

CONCLUSION

This study highlighted for the first time a natural molecule, extracted from Cocos nucifera flour, able to interact with the olfactory receptor OR2AT4 and promote the keratinocyte migration and thus the epithelialization.

The role of secondary interactions on the preferred conformers of the fenchone-ethanol complex

New atomic-level experimental data on the intermolecular non-covalent interactions between a common odorant and a relevant residue at odorant binding sites are reported. The preferred arrangements and binding interactions of fenchone, a common odorant and ethanol, a mimic of serine's side chain, have been unambiguously identified using a combination of high resolution rotational spectroscopy and computational methods. The observed conformers include homochiral (RR) and heterochiral (RS) conformers, with a slight preference for a heterochiral form, and exhibit primary OH-O hydrogen bonds between fenchone and ethanol. Secondary interactions play a key role in determining the relative configurations of fenchone and ethanol, and in shaping quite a flat potential energy surface, with many conformers close in energy and small barriers for interconversion.

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.

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.

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.