Making sense of olfaction through predictions of the 3-D structure and function of olfactory receptors

We used the MembStruk first principles computational technique to predict the three-dimensional (3-D) structure of six mouse olfactory receptors (S6, S18, S19, S25, S46 and S50) for which experimental odorant recognition profiles are available for a set of 24 odorants (4-9 carbons aliphatic alcohols, acids, bromo-acids and diacids). We used the HierDock method to scan each predicted OR structure for potential odorant binding site(s) and to calculate binding energies of each odorant in these binding sites. The calculated binding affinity profiles are in good agreement with experimental activation profiles, validating the predicted 3-D structures and the predicted binding sites. For each of the six ORs, the binding site is located between trans-membrane domains (TMs) 3-6, with contributions from extracellular loops 2 and 3. In particular, we find six residue positions in TM3 and TM6 to be consistently involved in the binding modes of the odorants. Indeed, the differences in the experimental recognition profiles can be explained on the basis of these critical residues alone. These predictions are also consistent with mutation data on ligand binding for catecholamine receptors and sequence hypervariability studies for ORs. Based on this analysis, we defined amino acid patterns associated with the recognition of short aliphatic alcohols and mono-acids. Using these two sequence fingerprints to probe the alignment of 869 OR sequences from the mouse genome, we identified 34 OR sequences matching the fingerprint for aliphatic mono-acids and 36 corresponding to the recognition pattern for aliphatic alcohols. We suggest that these two sets of ORs might function as basic arrays for uniquely recognizing aliphatic alcohols and acids. We screened a library of 89 additional molecules against the six ORs and found that this set of ORs is likely to respond to aldehydes and esters with longer carbon chains than their currently known agonists. We also find that compounds associated with the flavor in foods are often among the best calculated binding affinities. This suggests that physiologic ligands for these ORs may be found among aldehydes and esters associated with flavor.     

An improved bioluminescence‐based signaling assay for odor sensing with a yeast expressing a chimeric olfactory receptor

The goal of this work was to improve the bioluminescence‐based signaling assay system to create a practical application of a biomimetic odor sensor using an engineered yeast‐expressing olfactory receptors (ORs). Using the yeast endogenous pheromone receptor (Ste2p) as a model GPCR, we determined the suitable promoters for the firefly luciferase (luc) reporter and GPCR genes. Additionally, we deleted some genes to further improve the sensitivity of the luc reporter assay. By replacing the endogenous yeast G‐protein α‐subunit (Gpa1p) with the olfactory‐specific Gα_olf, the optimized yeast strain successfully transduced signal through both OR and yeast Ste2p. Our results will assist the development of a bioluminescence‐based odor‐sensing system using OR‐expressing yeast. Biotechnol. Bioeng. 2012; 109: 3143–3151. © 2012 Wiley Periodicals, Inc.     

Structural determinants for the interactions between muscarinic toxin 7 and muscarinic acetylcholine receptors

Muscarinic acetylcholine receptors (mAChRs) have five subtypes and play crucial roles in various physiological functions and pathophysiological processes. Poor subtype specificity of mAChR modulators has been an obstacle to discover new therapeutic agents. Muscarinic toxin 7 (MT7) is a natural peptide toxin with high selectivity for the M1 receptor. With three to five residues substituted, M3, M4, and M5 receptor mutants could bind to MT7 at nanomolar concentration as the M1 receptor. However, the structural mechanisms explaining MT7–mAChRs binding are still largely unknown. In this study, we constructed 10 complex models of MT7 and each mAChR subtype or its mutant, performed molecular dynamics simulations, and calculated the binding energies to investigate the mechanisms. Our results suggested that the structural determinants for the interactions on mAChRs were composed of some critical residues located separately in the extracellular loops of mAChRs, such as Glu4.56, Leu4.60, Glu/Gln4.63, Tyr4.65, Glu/Asp6.67, and Trp7.35. The subtype specificity of MT7 was attributed to the non‐conserved residues at positions 4.56 and 6.67. These structural mechanisms could facilitate the discovery of novel mAChR modulators with high subtype specificity and enhance the understanding of the interactions between ligands and G‐protein‐coupled receptors. Copyright © 2015 John Wiley & Sons, Ltd.     

Mapping protein energy landscapes with amide hydrogen exchange and mass spectrometry: I. A generalized model for a two-state protein and comparison with experiment

Protein amide hydrogen exchange (HDX) is a convoluted process, whose kinetics is determined by both dynamics of the protein and the intrinsic exchange rate of labile hydrogen atoms fully exposed to solvent. Both processes are influenced by a variety of intrinsic and extrinsic factors. A mathematical formalism initially developed to rationalize exchange kinetics of individual amide hydrogen atoms is now often used to interpret global exchange kinetics (e.g., as measured in HDX MS experiments). One particularly important advantage of HDX MS is direct visualization of various protein states by observing distinct protein ion populations with different levels of isotope labeling under conditions favoring correlated exchange (the so-called EX1 exchange mechanism). However, mildly denaturing conditions often lead to a situation where the overall HDX kinetics cannot be clearly classified as either EX1 or EX2. The goal of this work is to develop a framework for a generalized exchange model that takes into account multiple processes leading to amide hydrogen exchange, and does not require that the exchange proceed strictly via EX1 or EX2 kinetics. To achieve this goal, we use a probabilistic approach that assigns a transition probability and a residual protection to each equilibrium state of the protein. When applied to a small protein chymotrypsin inhibitor 2, the algorithm allows complex HDX patterns observed experimentally to be modeled with remarkably good fidelity. On the basis of the model we are now in a position to begin to extract quantitative dynamic information from convoluted exchange kinetics.     

Dynamics of the Hck-SH3 domain: comparison of experiment with multiple molecular dynamics simulations

Molecular dynamics calculations provide a method by which the dynamic properties of molecules can be explored over timescales and at a level of detail that cannot be obtained experimentally from NMR or X-ray analyses. Recent work (Philippopoulos M, Mandel AM, Palmer AG III, Lim C, 1997, Proteins 28:481-493) has indicated that the accuracy of these simulations is high, as measured by the correspondence of parameters extracted from these calculations to those determined through experimental means. Here, we investigate the dynamic behavior of the Src homology 3 (SH3) domain of hematopoietic cell kinase (Hck) via 5N backbone relaxation NMR studies and a set of four independent 4 ns solvated molecular dynamics calculations. We also find that molecular dynamics simulations accurately reproduce fast motion dynamics as estimated from generalized order parameter (S2) analysis for regions of the protein that have experimentally well-defined coordinates (i.e., stable secondary structural elements). However, for regions where the coordinates are not well defined, as indicated by high local root-mean-square deviations among NMR-determined structural family members or high B-factors/low electron density in X-ray crystallography determined structures, the parameters calculated from a short to moderate length (less than 5-10 ns) molecular dynamics trajectory are dependent on the particular coordinates chosen as a starting point for the simulation.     

Molecular organization and dynamics of the melatonin MT1 receptor/RGS20/Gi protein complex reveal asymmetry of receptor dimers for RGS and Gi coupling

Functional asymmetry of G‐protein‐coupled receptor (GPCR) dimers has been reported for an increasing number of cases, but the molecular architecture of signalling units associated to these dimers remains unclear. Here, we characterized the molecular complex of the melatonin MT₁ receptor, which directly and constitutively couples to G_i proteins and the regulator of G‐protein signalling (RGS) 20. The molecular organization of the ternary MT₁/G_i/RGS20 complex was monitored in its basal and activated state by bioluminescence resonance energy transfer between probes inserted at multiple sites of the complex. On the basis of the reported crystal structures of G_i and the RGS domain, we propose a model wherein one G_i and one RGS20 protein bind to separate protomers of MT₁ dimers in a pre‐associated complex that rearranges upon agonist activation. This model was further validated with MT₁/MT₂ heterodimers. Collectively, our data extend the concept of asymmetry within GPCR dimers, reinforce the notion of receptor specificity for RGS proteins and highlight the advantage of GPCRs organized as dimers in which each protomer fulfils its specific task by binding to different GPCR‐interacting proteins.     

The structure of the major transition state for folding of an FF domain from experiment and simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.     

Distinct second extracellular loop structures of the brain cannabinoid CB(1) receptor: implication in ligand binding and receptor function

The G-protein-coupled receptor (GPCR) second extracellular loop (E2) is known to play an important role in receptor structure and function. The brain cannabinoid (CB(1)) receptor is unique in that it lacks the interloop E2 disulfide linkage to the transmembrane (TM) helical bundle, a characteristic of many GPCRs. Recent mutation studies of the CB(1) receptor, however, suggest the presence of an alternative intraloop disulfide bond between two E2 Cys residues. Considering the oxidation state of these Cys residues, we determine the molecular structures of the 17-residue E2 in the dithiol form (E2(dithiol)) and in the disulfide form (E2(disulfide)) of the CB(1) receptor in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer, using a combination of simulated annealing and molecular dynamics simulation approaches. We characterize the CB(1) receptor models with these two E2 forms, CB(1)(E2(dithiol)) and CB(1)(E2(disulfide)), by analyzing interaction energy, contact number, core crevice, and cross correlation. The results show that the distinct E2 structures interact differently with the TM helical bundle and uniquely modify the TM helical topology, suggesting that E2 of the CB(1) receptor plays a critical role in stabilizing receptor structure, regulating ligand binding, and ultimately modulating receptor activation. Further studies on the role of E2 of the CB(1) receptor are warranted, particularly comparisons of the ligand-bound form with the present ligand-free form.     

Structure-based identification of binding sites,native ligands and potential inhibitors for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.     

Cyclic pentapeptides of chiral sequence DLDDL as scaffold for antagonism of G‐protein coupled receptors: Synthesis,activity and conformational analysis by NMR and molecular dynamics of ITF 1565 a substance P inhibitor

Under the hypotheses of a structurally related binding site for antagonists of G‐protein coupled receptors and the ability of cyclic pentapeptides of chiral sequence D ¹L ²D ³D ⁴L ⁵ to form rigid structures with which probe the pharmacophoric specificity of these receptors, inhibitors of substance P were designed based on available structure–activity relationships. ITF 1565, cyclo[D ‐Trp¹‐Pro²‐D ‐Lys³‐D ‐Trp⁴‐Phe⁵], antagonized substance P activity mediated by type 1 neurokinin receptor (NK1) whereas it acted weakly against NK2 and did not inhibit endothelin at all. The preferential conformation of the peptide was obtained from nmr spectroscopy and computer calculations, and shown to contain the same βII‐turn and γ′‐turn found in other cyclic pentapeptides with the same chiral sequence. The structure of the peptide was compared with that of the β‐D ‐glucose molecule that has been proposed as a semirigid scaffold for antagonists of G‐protein coupled receptors. The γ′‐turn of the cyclic peptide superimposed well with β‐D ‐glucose in the chair conformation. Furthermore, when the side chains were considered, the aromatic groups of the two molecules were found to generally overlap. These results support the view of G‐protein coupled receptors as possessing structurally similar binding sites for antagonists and suggest that cyclic pentapeptides of chiral sequence D ¹L ²D ³D ⁴L ⁵ may be useful as semirigid scaffolds for the design of antagonists of this family of receptors. © 1999 John Wiley & Sons, Inc. Biopoly 50: 211–219, 1999     

Progress toward heterologous expression of active G‐protein‐coupled receptors in Saccharomyces cerevisiae: Linking cellular stress response with translocation and trafficking

High-level expression of mammalian G-protein-coupled receptors (GPCRs) is a necessary step toward biophysical characterization and high-resolution structure determination. Even though many heterologous expression systems have been used to express mammalian GPCRs at high levels, many receptors are improperly trafficked or are inactive in these systems. En route to engineering a robust microbial host for GPCR expression, we have investigated the expression of 12 GPCRs in the yeast Saccharomyces cerevisiae, where all receptors are expressed at the mg/L scale. However, only the human adenosine A₂a (hA₂aR) receptor is active for ligand-binding and located primarily at the plasma membrane, whereas other tested GPCRs are mainly retained within the cell. Selective receptors associate with BiP, an ER-resident chaperone, and activated the unfolded protein response (UPR) pathway, which suggests that a pool of receptors may be folded incorrectly. Leader sequence cleavage of the expressed receptors was complete for the hA₂aR, as expected, and partially cleaved for hA₂bR, hCCR5R, and hD_2LR. Ligand-binding assays conducted on the adenosine family (hA₁R, hA₂aR, hA₂bR, and hA₃R) of receptors show that hA₂aR and hA₂bR, the only adenosine receptors that demonstrate leader sequence processing, display activity. Taken together, these studies point to translocation as a critical limiting step in the production of active mammalian GPCRs in S. cerevisiae.     

Membrane orientation and oligomerization of the melanocortin receptor accessory protein 2

The melanocortin receptor accessory protein 2 (MRAP2) plays a pivotal role in the regulation of several G protein–coupled receptors that are essential for energy balance and food intake. MRAP2 loss-of-function results in obesity in mammals. MRAP2 and its homolog MRAP1 have an unusual membrane topology and are the only known eukaryotic proteins that thread into the membrane in both orientations. In this study, we demonstrate that the conserved polybasic motif that dictates the membrane topology and dimerization of MRAP1 does not control the membrane orientation and dimerization of MRAP2. We also show that MRAP2 dimerizes through its transmembrane domain and can form higher-order oligomers that arrange MRAP2 monomers in a parallel orientation. Investigating the molecular details of MRAP2 structure is essential for understanding the mechanism by which it regulates G protein–coupled receptors and will aid in elucidating the pathways involved in metabolic dysfunction.     

Delineation of ligand binding and receptor signaling activities of purified P2Y receptors reconstituted with heterotrimeric G proteins

P2Y receptors are G protein coupled receptors that respond to extracellular nucleotides to promote a multitude of signaling events. Our laboratory has purified several P2Y receptors with the goal of providing molecular insight into their: (1) ligand binding properties, (2) G protein signaling selectivities, and (3) regulation by RGS proteins and other signaling cohorts. The human P2Y₁ receptor and the human P2Y₁₂ receptor, both of which are intimately involved in ADP-mediated platelet aggregation, were purified to near homogeneity and studied in detail. After high-level expression from recombinant baculovirus infection of Sf9 insect cells, approximately 50% of the receptors were successfully extracted with digitonin. Purification of nearly homogeneous epitope-tagged P2Y receptor was achieved using metal-affinity chromatography followed by other traditional chromatographic steps. Yields of purified P2Y receptors range from 10 to 100 g/l of infected cells. Once purified, the receptors were reconstituted in model lipid vesicles along with their cognate G proteins to assess receptor function. Agonist-promoted increases in steady-state GTPase assays demonstrated the functional activity of the reconstituted purified receptor. We have utilized this reconstitution system to assess the action of various nucleotide agonists and antagonists, the relative G protein selectivity, and the influence of other proteins, such as phospholipase C, on P2Y receptor-promoted signaling. Furthermore, we have identified the RGS expression profile of platelets and have begun to assess the action of these RGS proteins in a reconstituted P2Y receptor/G protein platelet model.     

Molecular dynamics simulations of transitions for ECD epidermal growth factor receptors show key differences between human and drosophila forms of the receptors

Recent X‐ray structural work on the Drosophila epidermal growth factor receptor (EFGR) has suggested an asymmetric dimer that rationalizes binding affinity measurements that go back decades (Alvarado et al., Cell 2010;142:568–579; Dawson et al., Structure 2007;15:942–954; Lemmon et al., Embo J 1997;16:281–294; Mattoon et al., Proc Natl Acad Sci USA 2004;101:923–928; Mayawala et al., Febs Lett 2005;579:3043–3047; Ozcan et al., Proc Natl Acad Sci USA 2006;103:5735–5740). This type of asymmetric structure has not been seen for the human EGF receptor family and it may or may not be important for function in that realm. We hypothesize that conformational changes in the Drosophila system have been optimized for the transition, whereas the barrier for the same transition is much higher in the human forms. To address our hypothesis we perform dynamic importance sampling (DIMS) (Perilla et al., J Comput Chem 2010;32:196–209) for barrier crossing transitions in both Drosophila and human EFGRs. For each set of transitions, we work from the hypothesis, based on results from the AdK system, that salt‐bridge pairs making and breaking connections are central to the conformational change. To evaluate the effectiveness of the salt‐bridges as drivers for the conformational change, we use the effective transfer entropy based on stable state MD calculations (Kamberaj and Der Vaart, Biophys J 2009;97:1747–1755) to define a reduced subset of degrees of freedom that seem to be important for driving the transition (Perilla and Woolf, J Chem Phys 2012;136:164101). Our results suggest that salt‐bridge making and breaking is not the dominant factor in driving the symmetric to asymmetric transition, but that instead it is a result of more concerted and correlated functional motions within a subset of the dimer structures. Furthermore, the analysis suggests that the set of residues involved in the transitions from the Drosophila relative to the human forms differs and that this difference in substate distributions relates to why the asymmetric form may be more common to Drosophila than to the human forms. We close with a discussion about the residues that may be changed in the human and the Drosophila forms to potentially shift the kinetics of the symmetric to asymmetric transition. Proteins 2013; 81:1113–1126. © 2013 Wiley Periodicals, Inc.     

Global fold of human cannabinoid type 2 receptor probed by solid‐state 13C‐, 15N‐MAS NMR and molecular dynamics simulations

The global fold of human cannabinoid type 2 (CB₂) receptor in the agonist‐bound active state in lipid bilayers was investigated by solid‐state ¹³C‐ and ¹⁵N magic‐angle spinning (MAS) NMR, in combination with chemical‐shift prediction from a structural model of the receptor obtained by microsecond‐long molecular dynamics (MD) simulations. Uniformly ¹³C‐ and ¹⁵N‐labeled CB₂ receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. ¹³C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of C_α, C_β, and C?O bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the ¹³C chemical shift distribution of C_α resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB₂. Thus the shape of the C_α band can be used as an indicator of CB₂ global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥1.5 ppm for carbons and ≥5.0 ppm for nitrogens). Simulated two‐dimensional ¹³C_α(i)? ¹³C?O(i) and ¹³C?O(i)? ¹⁵NH(i + 1) dipolar‐interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB₂ by solid‐state MAS NMR. Proteins 2014; 82:452–465. © 2013 Wiley Periodicals, Inc.     

Predicted structure of the extracellular region of ligand-gated ion-channel receptors shows SH2-like and SH3-like domains forming the ligand-binding site.

Fast synaptic neurotransmission is mediated by ligand-gated ion-channel (LGIC) receptors, which include receptors for acetylcholine, serotonin, GABA, glycine, and glutamate. LGICs are pentamers with extracellular ligand-binding domains and form integral membrane ion channels that are selective for cations (acetylcholine and serotonin 5HT3 receptors) or anions (GABAA and glycine receptors and the invertebrate glutamate-binding chloride channel). They form a protein superfamily with no sequence similarity to any protein of known structure. Using a 1D-3D structure mapping approach, we have modeled the extracellular ligand-binding domain based on a significant match with the SH2 and SH3 domains of the biotin repressor structure. Refinement of the model based on knowledge of the large family of SH2 and SH3 structures, sequence alignments, and use of structure templates for loop building, allows the prediction of both monomer and pentamer models. These are consistent with medium-resolution electron microscopy structures and with experimental structure/function data from ligand-binding, antibody-binding, mutagenesis, protein-labeling and subunit-linking studies, and glycosylation sites. Also, the predicted polarity of the channel pore calculated from electrostatic potential maps of pentamer models of superfamily members is consistent with known ion selectivities. Using the glycine receptor alpha 1 subunit, which forms homopentamers, the monomeric and pentameric models define the agonist and antagonist (strychnine) binding sites to a deep crevice formed by an extended loop, which includes the invariant disulfide bridge, between the SH2 and SH3 domains. A detailed binding site for strychnine is reported that is in strong agreement with known structure/function data. A site for interaction of the extracellular ligand-binding domain with the activation of the M2 transmembrane helix is also suggested.     

Design of novel quinazolinone derivatives as inhibitors for 5HT7 receptor

To study the pharmacophore properties of quinazolinone derivatives as 5HT₇ inhibitors, 3D QSAR methodologies, namely Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) were applied, partial least square (PLS) analysis was performed and QSAR models were generated. The derived model showed good statistical reliability in terms of predicting the 5HT₇ inhibitory activity of the quinazolione derivative, based on molecular property fields like steric, electrostatic, hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. This is evident from statistical parameters like q² (cross validated correlation coefficient) of 0.642, 0.602 and r² (conventional correlation coefficient) of 0.937, 0.908 for CoMFA and CoMSIA respectively. The predictive ability of the models to determine 5HT₇ antagonistic activity is validated using a test set of 26 molecules that were not included in the training set and the predictive r² obtained for the test set was 0.512 & 0.541. Further, the results of the derived model are illustrated by means of contour maps, which give an insight into the interaction of the drug with the receptor. The molecular fields so obtained served as the basis for the design of twenty new ligands. In addition, ADME (Adsorption, Distribution, Metabolism and Elimination) have been calculated in order to predict the relevant pharmaceutical properties, and the results are in conformity with required drug like properties.     

Structure of a GRK5-Calmodulin Complex Reveals Molecular Mechanism of GRK Activation and Substrate Targeting

1. Download : Download high-res image (234KB)
2. Download : Download full-size image