Beyond modeling: all-atom olfactory receptor model simulations

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.

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.

Sense of Smell: Structural, Functional, Mechanistic Advancements and Challenges in Human Olfactory Research

Olfaction, the sense of smell detects and discriminate odors as well as social cues which influence our innate responses. The olfactory system in human beings is found to be weak as compared to other animals; however, it seems to be very precise. It can detect and discriminate millions of chemical moieties (odorants) even in minuscule quantities. The process initiates with the binding of odorants to specialized olfactory receptors, encoded by a large family of Olfactory Receptor (OR) genes belonging to the G-protein-coupled receptor superfamily. Stimulation of ORs converts the chemical information encoded in the odorants, into respective neuronal action-potentials which causes depolarization of olfactory sensory neurons. The olfactory bulb relays this signal to different parts of the brain for processing. Odors are encrypted using a combinatorial approach to detect a variety of chemicals and encode their unique identity. The discovery of functional OR genes and proteins provided an important information to decipher the genomic, structural and functional basis of olfaction. ORs constitute 17 gene families, out of which 4 families were reported to contain more than hundred members each. The olfactory machinery is not limited to GPCRs; a number of non- GPCRs is also employed to detect chemosensory stimuli. The article provides detailed information about such olfaction machinery, structures, transduction mechanism, theories of odor perception, and challenges in the olfaction research. It covers the structural, functional and computational studies carried out in the olfaction research in the recent past.

Human Olfactory Receptors: Novel Cellular Functions Outside of the Nose

Olfactory receptors (ORs) are not exclusively expressed in the olfactory sensory neurons; they are also observed outside of the olfactory system in all other human tissues tested to date, including the testis, lung, intestine, skin, heart, and blood. Within these tissues, certain ORs have been determined to be exclusively expressed in only one tissue, whereas other ORs are more widely distributed in many different tissues throughout the human body. For most of the ectopically expressed ORs, limited data are available for their functional roles. They have been shown to be involved in the modulation of cell-cell recognition, migration, proliferation, the apoptotic cycle, exocytosis, and pathfinding processes. Additionally, there is a growing body of evidence that they have the potential to serve as diagnostic and therapeutic tools, as ORs are highly expressed in different cancer tissues. Interestingly, in addition to the canonical signaling pathways activated by ORs in olfactory sensory neurons, alternative pathways have been demonstrated in nonolfactory tissues. In this review, the existing data concerning the expression, as well as the physiological and pathophysiological functions, of ORs outside of the nose are highlighted to provide insights into future lines of research.

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.

Molecular Dynamics Simulations of Membrane-Bound STIM1 to Investigate Conformational Changes during STIM1 Activation upon Calcium Release

Calcium is involved in important intracellular processes, such as intracellular signaling from cell membrane receptors to the nucleus. Typically, calcium levels are kept at less than 100 nM in the nucleus and cytosol, but some calcium is stored in the endoplasmic reticulum (ER) lumen for rapid release to activate intracellular calcium-dependent functions. Stromal interacting molecule 1 (STIM1) plays a critical role in early sensing of changes in the ER's calcium level, especially when there is a sudden release of stored calcium from the ER. Inactive STIM1, which has a bound calcium ion, is activated upon ion release. Following activation of STIM1, there is STIM1-assisted initiation of extracellular calcium entry through channels in the cell membrane. This extracellular calcium entering the cell then amplifies intracellular calcium-dependent actions. At the end of the process, ER levels of stored calcium are reestablished. The main focus of this work was to study the conformational changes accompanying homo- or heterodimerization of STIM1. For this purpose, the ER luminal portion of STIM1 (residues 58-236), which includes the sterile alpha motif (SAM) domain plus the calcium-binding EF-hand domains 1 and 2 attached to the STIM1 transmembrane region (TM), was modeled and embedded in a virtual membrane. Next, molecular dynamics simulations were performed to study the conformational changes that take place during STIM1 activation and subsequent protein-protein interactions. Indeed, the simulations revealed exposure of residues in the EF-hand domains, which may be important for dimerization steps. Altogether, understanding conformational changes in STIM1 can help in drug discovery when targeting this key protein in intracellular calcium functions.

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/

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.

Predicted 3D structures of olfactory receptors with details of odorant binding to OR1G1

Olfactory receptors (ORs) are responsible for mediating the sense of smell; they allow humans to recognize an enormous number of odors but the connection between binding and perception is not known. We predict the ensemble of low energy structures for the human OR1G1 (hOR1G1) and also for six other diverse ORs, using the G protein-coupled receptor Ensemble of Structures in Membrane BiLayer Environment complete sampling method that samples 13 trillion different rotations and tilts using four different templates to predict the 24 structures likely to be important in binding and activation. Our predicted most stable structures of hOR1G1 have a salt-bridge between the conserved D3.49 and K6.30 in the D(E)RY region, that we expect to be associated with an inactive form. The hOR1G1 structure also has specific interaction in transmembrane domains (TMD) 3-6 (E3.39 and H6.40), which is likely an important conformational feature for all hORs because of the ~94 to 98 % conservation among all hOR sequences. Of the five ligands studied (nonanal, 9-decen-1-ol, 1-nonanol, camphor, and n-butanal), we find that the 4 expected to bind lead to similar binding energies with nonanol the strongest.

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.

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.

Databases in SenseLab for the genomics, proteomics, and function of olfactory receptors

We present here, the salient aspects of three databases: Olfactory Receptor Database (ORDB) is a repository of genomics and proteomics information of ORs; OdorDB stores information related to odorous compounds, specifically identifying those that have been shown to interact with olfactory rectors; and OdorModelDB disseminates information related to computational models of olfactory receptors (ORs). The data stored among these databases is integrated. Presented in this chapter are descriptions of these resources, which are part of the SenseLab suite of databases, a discussion of the computational infrastructure that enhances the efficacy of information storage, retrieval, dissemination, and automated data population from external sources.

Peeping into human renal calcium oxalate stone matrix: characterization of novel proteins involved in the intricate mechanism of urolithiasis

BACKGROUND

The increasing number of patients suffering from urolithiasis represents one of the major challenges which nephrologists face worldwide today. For enhancing therapeutic outcomes of this disease, the pathogenic basis for the formation of renal stones is the need of hour. Proteins are found as major component in human renal stone matrix and are considered to have a potential role in crystal-membrane interaction, crystal growth and stone formation but their role in urolithiasis still remains obscure.

METHODS

Proteins were isolated from the matrix of human CaOx containing kidney stones. Proteins having MW>3 kDa were subjected to anion exchange chromatography followed by molecular-sieve chromatography. The effect of these purified proteins was tested against CaOx nucleation and growth and on oxalate injured Madin-Darby Canine Kidney (MDCK) renal epithelial cells for their activity. Proteins were identified by Matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF MS) followed by database search with MASCOT server. In silico molecular interaction studies with CaOx crystals were also investigated.

RESULTS

Five proteins were identified from the matrix of calcium oxalate kidney stones by MALDI-TOF MS followed by database search with MASCOT server with the competence to control the stone formation process. Out of which two proteins were promoters, two were inhibitors and one protein had a dual activity of both inhibition and promotion towards CaOx nucleation and growth. Further molecular modelling calculations revealed the mode of interaction of these proteins with CaOx at the molecular level.

CONCLUSIONS

We identified and characterized Ethanolamine-phosphate cytidylyltransferase, Ras GTPase-activating-like protein, UDP-glucose:glycoprotein glucosyltransferase 2, RIMS-binding protein 3A, Macrophage-capping protein as novel proteins from the matrix of human calcium oxalate stone which play a critical role in kidney stone formation. Thus, these proteins having potential to modulate calcium oxalate crystallization will throw light on understanding and controlling urolithiasis in humans.