The mouse eugenol odorant receptor: structural and functional plasticity of a broadly tuned odorant binding pocket

Molecular interactions of odorants with their olfactory receptors (ORs) are of central importance for the ability of the mammalian olfactory system to detect and discriminate a vast variety of odors with a limited set of receptors. How a particular OR binds and distinguishes different odorant molecules remains largely unknown on a structural basis. Here we investigated this question for the mouse eugenol receptor (mOR-EG). By screening a large odorant library, we discovered a wide range of chemical structures activating the receptor in heterologous mammalian cells. Potent agonists comprise (i) benzene, (ii) cyclohexane, or (iii) polycyclic structures substituted with alcohol, aldehyde, keto, ether, or esterified carboxylic groups. To detect those amino acids within the receptor that are in contact with a particular bound odorant molecule, we investigated how distinct mOR-EG point mutants were activated by the different odorant agonists found for the wild-type receptor. We identified 11 amino acids as a part of the receptor's ligand binding pocket. Molecular modeling predicted 10 of these residues in transmembrane helices TM3-TM6 and one in the extracellular loop between TM2 and TM3. These amino acids participate in odorant binding with variable importance depending on the type of odorant, revealing functional "fingerprints" of ligand-receptor interactions.     

Constitutive activation of the angiotensin II type 1 receptor alters the spatial proximity of transmembrane 7 to the ligand-binding pocket

Activation of G protein-coupled receptors by agonists involves significant movement of transmembrane domains (TM) following binding of agonist. The underlying structural mechanism by which receptor activation takes place is largely unknown but can be inferred by detecting variability within the environment of the ligand-binding pocket, which constitutes a water-accessible crevice surrounded by the seven TM helices. Using the substituted cysteine accessibility method, we initially identified those residues within the seventh transmembrane domain (TM7) of wild type angiotensin II type 1 (AT1) receptor that contribute to forming the binding site pocket. We have substituted successively TM7 residues ranging from Ile276 to Tyr302 to cysteine. Treatment of A277C, V280C, T282C, A283C, I286C, A291C, and F301C mutant receptors with the charged sulfhydryl-specific alkylating agent MTSEA significantly inhibited ligand binding, which suggests that these residues orient themselves within the water-accessible binding pocket of the AT1 receptor. Interestingly, this pattern of acquired MTSEA sensitivity was greatly reduced for TM7 reporter cysteines engineered in a constitutively active mutant of the AT1 receptor. Our data suggest that upon activation, TM7 of the AT1 receptor goes through a pattern of helical movements that results in its distancing from the binding pocket per se. These studies support accumulating evidence whereby elements of TM7 of class A GPCRs promote activation of the receptor through structural rearrangements.     

Ligand binding to the human MT2 melatonin receptor: the role of residues in transmembrane domains 3, 6, and 7

To better understand the mechanism of interactions between G-protein-coupled melatonin receptors and their ligands, our previously reported homology model of human MT2 receptor with docked 2-iodomelatonin was further refined and used to select residues within TM3, TM6, and TM7 potentially important for receptor-ligand interactions. Selected residues were mutated and radioligand-binding assay was used to test the binding affinities of hMT2 receptors transiently expressed in HEK293 cells. Our data demonstrate that residues N268 and A275 in TM6 as well as residues V291 and L295 in TM7 are essential for 2-iodomelatonin binding to the hMT2 receptor, while TM3 residues M120, G121, V124, and I125 may participate in binding of other receptor agonists and/or antagonists. Presented data also hint at possible specific interaction between the side-chain of Y188 in second extracellular loop and N-acetyl group of 2-iodomelatonin.     

Dissecting the molecular mechanism of drosophila odorant receptors through activity modeling and comparative analysis

To gain insight into the molecular mechanism of odorant receptors (ORs) in Drosophila species, we developed a Quantitative Structure Activity Relationship (QSAR) model that predicts experimentally measured electrophysiological activities between 24 D. melanogaster ORs and 108 odorants. Although the model is limited by the tested odorants,analyzing the model allowed dissection of specific topological and chemical properties necessary for an odorant to elicit excitatory or inhibitory receptor response. Linear odorants with five to eight nonhydrogen atoms at the main chain and hydrogen‐bond acceptor and/or hydrogen‐bond donor at its ends were found to stimulate strong excitatory response. A comparative sequence analysis of 90 ORs in 15 orthologous groups identified 15 putative specificity‐determining residues (SDRs) and 15 globally conserved residues that we postulate as functionally key residues. Mapping to a model of secondary structure resulted in 14 out of 30 key residues locating to the transmembrane (TM) domains. Twelve residues, including six SDRs and six conserved residues, are located at the extracellular halves of the TM domains. Combining the evidence from the QSAR modeling and the comparative sequence analysis, we hypothesize that the Drosophila ORs accept odorants into a binding pocket located on the extracellular halves of its TM domains. The QSAR modeling suggests that the binding pocket is around 15 Å in depth and about 6 Å in width. Twelve mainly polar or charged key residues, both SDRs and conserved, are located inthis pocket and postulated to distinguish docked odorants via primarily geometry fitting and hydrogen‐bond interaction. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Structural activation pathways from dynamic olfactory receptor-odorant interactions

We have simulated an odor ligand's dynamic behavior in the binding region of an olfactory receptor (OR). Our short timescale computational studies (up to 200 ps) have helped identify unprecedented postdocking ligand behavior of ligands. From in vacuo molecular dynamics simulations of interactions between models of rat OR I7 and 10 aldehyde ligands, we have identified a dissociative pathway along which the ligand exits and enters the OR-binding pocket--a transit event. The ligand's transit through the receptor's binding region may mark the beginning of a signal transduction cascade leading to odor recognition. We have graphically traced the rotameric changes in key OR amino acid side chains during the transit. Our results have helped substantiate or refute previously held notions of amino acid contribution to ligand stability in the binding pocket. Our observations of ligand activity when compared to those of experimental (electroolfactogram response) OR-activation studies provide a view to predicting the stability of ligands in the binding pocket as a precursor to OR activation by the ligand.     

A role for transmembrane domains V and VI in ligand binding and maturation of the angiotensin II AT1 receptor

Several studies have proposed that angiotensin II (Ang II) binds to its receptor AT1 through interactions with residues in helices V and VI, suggesting that the distance between these helices is crucial for ligand binding. Based on a 3D model of AT1 in which the C-terminus of Ang II is docked, we identified the hydrophobic residues of TM V and VI pointing towards the external face of the helices, which may play a role in the structure of the binding pocket and in the structural integrity of the receptor. We performed a systematic mutagenesis study of these residues and examined the binding, localization, maturation, and dimerization of the mutated receptors. We found that mutations of hydrophobic residues to alanine in helix V do not alter binding, whereas mutations to glutamate lead to loss of binding without a loss in cell surface expression, suggesting that the external face of helix V may not directly participate in binding, but may rather contribute to the structure of the binding pocket. In contrast, mutations of hydrophobic residues to glutamate in helix VI lead to a loss in cell surface expression, suggesting that the external surface of helix VI plays a structural role and ensures correct folding of the receptor.     

Vasopressin V2 Receptor/Bioligand Interactions

We predict some essential interactions between the V2 vasopressin renal receptor (V2R) and its agonists [Arg⁸]vasopressin (AVP) and [D-Arg⁸]vasopressin (DAVP), and the non-peptide antagonist OPC-31260. V2R controls antidiuresis and belongs to the superfamily of heptahelical transmembrane (7TM) G-protein-coupled receptors (GPCRs). The receptor was built, the ligands were docked and the structures relaxed using advanced molecular modeling techniques. Docked agonists and antagonists appear to prefer similar V2R compartments. A number of receptor amino acid residues are indicated, mainly in the TM3–TM7 helices, as potentially important in ligand binding. Many of these residues are invariant for either the GPCR superfamily or the subfamily of related (vasopressin V2R, V1aR and V1bR and oxytocin OR) receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for ligand affinity [Mouillac et al., J. Biol. Chem., 270 (1995) 25771].     

Predicted 3-D structures for mouse I7 and rat I7 olfactory receptors and comparison of predicted odor recognition profiles with experiment

The first step in the perception of an odor is the activation of one or more olfactory receptors (ORs) following binding of the odorant molecule to the OR. In order to initiate the process of determining how the molecular level receptor-odorant interactions are related to odor perception, we used the MembStruk computational method to predict the three-dimensional (3-D) structure of the I7 OR for both mouse and rat. We then used the HierDock ligand docking computational method to predict the binding site and binding energy for the library of 56 odorants to these receptors for which experiment response data are now available. We find that the predicted 3-D structures of the mouse and rat I7 OR lead to predictions of odorant binding that are in good agreement with the experimental results, thus validating the accuracy of both the 3-D structure and the predicted binding site. In particular we predict that heptanal and octanal both bind strongly to both mouse and rat I7 ORs, which conflicts with the older literature but agrees with recent experiments. To provide the basis of additional validations of our 3-D structures, we also report the odorant binding site for a new odorant (8-hydroxy-octanal) with a novel functionality designed to bind strongly to mouse I7. Such validated computational methods should be very useful in predicting the structure and function of many other ORs.     

Vasopressin V2 Receptor/Bioligand Interactions

Summary We predict some essential interactions between the V2 vasopressin renal receptor (V2R) and its agonists [Arg⁸]vasopressin (AVP) and [D-Arg⁸]vasopressin (DAVP), and the non-peptide antagonist OPC-31260. V2R controls antidiuresis and belongs to the superfamily of heptahelical transmembrane (7TM) G-protein-coupled receptors (GPCRs). The receptor was built, the ligands were docked and the structures relaxed using advanced molecular modeling techniques. Docked agonists and antagonists appear to prefer similar V2R compartments. A number of receptor amino acid residues are indicated, mainly in the TM3-TM7 helices, as potentially important in ligand binding. Many of these residues are invariant for either the GPCR superfamily or the subfamily of related (vasopressin V2R, V1aR and V1bR and oxytocin OR) receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for ligand affinity [Mouillac et al., J. Biol. Chem., 270 (1995) 25771].     

Binding modes of dihydroquinoxalinones in a homology model of bradykinin receptor 1

We report the first homology model of human bradykinin receptor B1 generated from the crystal structure of bovine rhodopsin as a template. Using an automated docking procedure, two B1 receptor antagonists of the dihydroquinoxalinone structural class were docked into the receptor model. Site-directed mutagenesis data of the amino acid residues in TM1, TM3, TM6, and TM7 were incorporated to place the compounds in the binding site of the homology model of the human B1 bradykinin receptor. The best pose in agreement with the mutation data was selected for detailed study of the receptor-antagonist interaction. To test the model, the calculated antagonist-receptor binding energy was correlated with the experimentally measured binding affinity (K(i)) for nine dihydroquinoxalinone analogs. The model was used to gain insight into the molecular mechanism for receptor function and to optimize the dihydroquinoxalinone analogs.     

Molecular similarities in the ligand binding pockets of an odorant receptor and the metabotropic glutamate receptors

The 5.24 odorant receptor is an amino acid sensing receptor that is expressed in the olfactory epithelium of fish. The 5.24 receptor is a G-protein-coupled receptor that shares amino acid sequence identity to mammalian pheromone receptors, the calcium-sensing receptor, the T1R taste receptors, and the metabotropic glutamate receptors (mGluRs). It is most potently activated by the basic amino acids L-lysine and L-arginine. In this study we generated a homology model of the ligand binding domain of the 5.24 receptor based on the crystal structure of mGluR1 and examined the proposed lysine binding pocket using site-directed mutagenesis. Mutants of truncated glycosylated versions of the receptor containing only the extracellular domain were analyzed in a radioligand binding assay, whereas the analogous full-length membrane-bound mutants were studied using a fluorescence-based functional assay. In silico analysis predicted that aspartate 388 interacts with the terminal amino group on the side chain of the docked lysine molecule. This prediction was supported by experimental observations demonstrating that mutation of this residue caused a 26-fold reduction in the affinity for L-lysine but virtually no change in the affinity for the polar amino acid L-glutamine. In addition, mutations in four highly conserved residues (threonine 175, tyrosine 223, and aspartates 195 and 309) predicted to establish interactions with the alpha amino group of the bound lysine ligand greatly reduced or eliminated binding and receptor activation. These results define the essential features of amino acid selectivity within the 5.24 receptor binding pocket and highlight an evolutionarily conserved motif required for ligand recognition in amino acid activated receptors in the G-protein-coupled receptor superfamily.     

Piezoelectric olfactory biosensor: ligand specificity and dose-dependence of an olfactory receptor expressed in a heterologous cell system

An olfactory receptor protein of rats, I7, was expressed on the surface of human embryonic kidney (HEK)-293 cells. For targeting and detecting the protein, rho-tag import sequence was fused with the I7 protein. The olfactory receptor was expressed on the plasma membrane of HEK-293 cells, and stable cell lines regulated by an inducer were obtained. The expression on the cell surface was confirmed by immunocytochemical and Western blotting methods, and the binding of specific odorant molecules to the olfactory receptor was measured using quartz crystal microbalance (QCM). The results for QCM coated with cells containing the olfactory receptor showed that the expressed protein I7 strongly interacted with octyl aldehyde (octanal), which is an odorant specific to the I7 protein. Several other odorants were tested, and the results showed that I7 interacted differently with them. The QCM response to the serial concentrations of octyl aldehyde showed that the response is dose dependent. All these results indicate that the I7 receptor protein expressed on the surface of the heterologous cell system is sensitive to the specific odorant and can be used for the quantitative measurement of the odorant.     

An olfactory receptor pseudogene whose function emerged in humans: a case study in the evolution of structure–function in GPCRs

Human olfactory receptor, hOR17-210, is identified as a pseudogene in the human genome. Experimental data has shown however, that the gene product of frame-shifted, cloned hOR17-210 cDNA was able to bind an odorant-binding protein and is narrowly tuned for excitation by cyclic ketones. Supported by experimental results, we used the bioinformatics methods of sequence analysis (genome-wide and pair-wise), computational protein modeling and docking, to show that functionality in this receptor is retained due to sequence-structure features not previously observed in mammalian ORs. This receptor does not possess the first two transmembrane helical domains (of seven typically seen in GPCRs). It however, possesses an additional TM that has not been observed in other human olfactory receptors. By incorporating these novel structural features, we created two putative models for this receptor. We also docked odor ligands that were experimentally shown to bind hOR17-210. We show how and why structural modifications of OR17-210 do not hinder this receptor's functionality. Our studies reveal that novel gene rearrangements that result in sequence and structural diversity may have a bearing on OR and GPCR function and evolution.     

An evaluation of automated in silico ligand docking of amino acid ligands to Family C G-protein coupled receptors

Family C G-protein coupled receptors (GPCRs) consist of the metabotropic glutamate receptors (mGluRs), the calcium-sensing receptor (CaSR), the T1R taste receptors, the GABA(B) receptor, the V2R pheromone receptors, and several chemosensory receptors. A common feature of Family C receptors is the presence of an amino acid binding pocket. The objective of this study was to evaluate the ability of the automatic docking program FlexX to predict the favored amino acid ligand at several Family C GPCRs. The docking process was optimized using the crystal structure of mGluR1 and the 20 amino acids were docked into homology models of the CaSR, the 5.24 chemosensory receptor, and the GPRC6A amino acid receptor. Under optimized docking conditions, glutamate was docked in the binding pocket of mGluR1 with a root mean square deviation of 1.56 angstroms from the co-crystallized glutamate structure and was ranked as the best ligand with a significantly better FlexX score compared to all other amino acids. Ligand docking to a homology model of the 5.24 receptor gave generally correct predictions of the favored amino acids, while the results obtained with models of GPRC6A and the CaSR showed that some of the favored amino acids at these receptors were correctly predicted, while a few other top scoring amino acids appeared to be false positives. We conclude that with certain caveats, FlexX can be successfully used to predict preferred ligands at Family C GPCRs.     

Molecular Basis for Benzodiazepine Agonist Action at the Type 1 Cholecystokinin Receptor

Understanding the molecular basis of drug action can facilitate development of more potent and selective drugs. Here, we explore the molecular basis for action of a unique small molecule ligand that is a type 1 cholecystokinin (CCK) receptor agonist and type 2 CCK receptor antagonist, GI181771X. We characterize its binding utilizing structurally related radioiodinated ligands selective for CCK receptor subtypes that utilize the same allosteric ligand-binding pocket, using wild-type receptors and chimeric constructs exchanging the distinct residues lining this pocket. Intracellular calcium assays were performed to determine biological activity. Molecular models for docking small molecule agonists to the type 1 CCK receptor were developed using a ligand-guided refinement approach. The optimal model was distinct from the previous antagonist model for the same receptor and was mechanistically consistent with the current mutagenesis data. This study revealed a key role for Leu^7.39 that was predicted to interact with the isopropyl group in the N1 position of the benzodiazepine that acts as a “trigger” for biological activity. The molecular model was predictive of binding of other small molecule agonists, effectively distinguishing these from 1065 approved drug decoys with an area under curve value of 99%. The model also selectively enriched for agonist compounds, with 130 agonists identified by ROC analysis when seeded in 2175 non-agonist ligands of the type 1 CCK receptor (area under curve 78%). Benzodiazepine agonists in this series docked in consistent pose within this pocket, with a key role played by Leu^7.39, whereas the role of this residue was less clear for chemically distinct agonists.     

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.     

Identification of the critical residues of bradykinin receptor B1 for interaction with the kinins guided by site-directed mutagenesis and molecular modeling

We report the critical residues for the interaction of the kinins with human bradykinin receptor 1 (B1) using site-directed mutagenesis in conjunction with molecular modeling of the binding modes of the kinins in the homology model of the B1 receptor. Mutation of Lys118 in transmembrane (TM) helix 3, Ala270 in TM6, and Leu294 in TM7 causes a significant decrease in the affinity for the peptide agonists des-Arg10kallidin (KD) and des-Arg9BK but not the peptide antagonist des-Arg10Leu9KD. In contrast, mutations in TM2, TM3, TM6, and TM7 cause a significant decrease in the affinity for both the peptide agonists and the antagonist. These data indicate that the B1 bradykinin binding pocket for agonists and antagonists is similar, but the manners in which they interact with the receptor do not completely overlap. Therefore, there is a potential to influence the receptor's ligand selectivity.     

Molecular basis for binding and subtype selectivity of 1,4-benzodiazepine antagonist ligands of the cholecystokinin receptor

Allosteric binding pockets in peptide-binding G protein-coupled receptors create opportunities for the development of small molecule drugs with substantial benefits over orthosteric ligands. To gain insights into molecular determinants for this pocket within type 1 and 2 cholecystokinin receptors (CCK1R and CCK2R), we prepared a series of receptor constructs in which six distinct residues in TM2, -3, -6, and -7 were reversed. Two novel iodinated CCK1R- and CCK2R-selective 1,4-benzodiazepine antagonists, differing only in stereochemistry at C3, were used. When all six residues within CCK1R were mutated to corresponding CCK2R residues, benzodiazepine selectivity was reversed, yet peptide binding selectivity was unaffected. Detailed analysis, including observations of gain of function, demonstrated that residues 6.51, 6.52, and 7.39 were most important for binding the CCK1R-selective ligand, whereas residues 2.61 and 7.39 were most important for binding CCK2R-selective ligand, although the effect of substitution of residue 2.61 was likely indirect. Ligand-guided homology modeling was applied to wild type receptors and those reversing benzodiazepine binding selectivity. The models had high predictive power in enriching known receptor-selective ligands from related decoys, indicating a high degree of precision in pocket definition. The benzodiazepines docked in similar poses in both receptors, with C3 urea substituents pointing upward, whereas different stereochemistry at C3 directed the C5 phenyl rings and N1 methyl groups into opposite orientations. The geometry of the binding pockets and specific interactions predicted for ligand docking in these models provide a molecular framework for understanding ligand selectivity at these receptor subtypes. Furthermore, the strong predictive power of these models suggests their usefulness in the discovery of lead compounds and in drug development programs.