Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis

Calcium-Driven In Silico Inactivation of a Human Olfactory Receptor

Conformational changes as well as molecular determinants related to the activation and inactivation of olfactory receptors are still poorly understood due to the intrinsic difficulties in the structural determination of this GPCR family. Here, we perform, for the first time, the in silico inactivation of human olfactory receptor OR51E2, highlighting the possible role of calcium in this receptor state transition. Using molecular dynamics simulations, we show that a divalent ion in the ion binding site, coordinated by two acidic residues at positions 2.50 and 3.39 conserved across most ORs, stabilizes the receptor in its inactive state. In contrast, protonation of the same two acidic residues is not sufficient to drive inactivation within the microsecond timescale of our simulations. Our findings suggest a novel molecular mechanism for OR inactivation, potentially guiding experimental validation and offering insights into the possible broader role of divalent ions in GPCR signaling.

Machine Learning-Based Modeling of Olfactory Receptors in Their Inactive State: Human OR51E2 as a Case Study

Atomistic-level investigation of olfactory receptors (ORs) is a challenging task due to the experimental/computational difficulties in the structural determination/prediction for members of this family of G-protein coupled receptors. Here, we have developed a protocol that performs a series of molecular dynamics simulations from a set of structures predicted de novo by recent machine learning algorithms and apply it to a well-studied receptor, the human OR51E2. Our study demonstrates the need for simulations to refine and validate such models. Furthermore, we demonstrate the need for the sodium ion at a binding site near D^2.50 and E^3.39 to stabilize the inactive state of the receptor. Considering the conservation of these two acidic residues across human ORs, we surmise this requirement also applies to the other ∼400 members of this family. Given the almost concurrent publication of a CryoEM structure of the same receptor in the active state, we propose this protocol as an in silico complement to the growing field of ORs structure determination.

Inhibiting a promiscuous GPCR: iterative discovery of bitter taste receptor ligands

The human GPCR family comprises circa 800 members, activated by hundreds of thousands of compounds. Bitter taste receptors, TAS2Rs, constitute a large and distinct subfamily, expressed orally and extra-orally and involved in physiological and pathological conditions. TAS2R14 is the most promiscuous member, with over 150 agonists and 3 antagonists known prior to this study. Due to the scarcity of inhibitors and to the importance of chemical probes for exploring TAS2R14 functions, we aimed to discover new ligands for this receptor, with emphasis on antagonists. To cope with the lack of experimental structure of the receptor, we used a mixed experimental/computational methodology which iteratively improved the performance of the predicted structure. The increasing number of active compounds, obtained here through experimental screening of FDA-approved drug library, and through chemically synthesized flufenamic acid derivatives, enabled the refinement of the binding pocket, which in turn improved the structure-based virtual screening reliability. This mixed approach led to the identification of 10 new antagonists and 200 new agonists of TAS2R14, illustrating the untapped potential of rigorous medicinal chemistry for TAS2Rs. 9% of the ~ 1800 pharmaceutical drugs here tested activate TAS2R14, nine of them at sub-micromolar concentrations. The iterative framework suggested residues involved in the activation process, is suitable for expanding bitter and bitter-masking chemical space, and is applicable to other promiscuous GPCRs lacking experimental structures.

Structural basis for strychnine activation of human bitter taste receptor TAS2R46

Taste sensing is a sophisticated chemosensory process, and bitter taste perception is mediated by type 2 taste receptors (TAS2Rs), or class T G protein–coupled receptors. Understanding the detailed molecular mechanisms behind taste sensation is hindered by a lack of experimental receptor structures. Here, we report the cryo–electron microscopy structures of human TAS2R46 complexed with chimeric mini–G protein gustducin, in both strychnine-bound and apo forms. Several features of TAS2R46 are disclosed, including distinct receptor structures that compare with known GPCRs, a new “toggle switch,” activation-related motifs, and precoupling with mini–G protein gustducin. Furthermore, the dynamic extracellular and more-static intracellular parts of TAS2R46 suggest possible diverse ligand-recognition and activation processes. This study provides a basis for further exploration of other bitter taste receptors and their therapeutic applications.

BitterMatch: recommendation systems for matching molecules with bitter taste receptors

Bitterness is an aversive cue elicited by thousands of chemically diverse compounds. Bitter taste may prevent consumption of foods and jeopardize drug compliance. The G protein-coupled receptors for bitter taste, TAS2Rs, have species-dependent number of subtypes and varying expression levels in extraoral tissues. Molecular recognition by TAS2R subtypes is physiologically important, and presents a challenging case study for ligand-receptor matchmaking. Inspired by hybrid recommendation systems, we developed a new set of similarity features, and created the BitterMatch algorithm that predicts associations of ligands to receptors with ~ 80% precision at ~ 50% recall. Associations for several compounds were tested in-vitro, resulting in 80% precision and 42% recall. The encouraging performance was achieved by including receptor properties and integrating experimentally determined ligand-receptor associations with chemical ligand-to-ligand similarities.

BitterMatch can predict off-targets for bitter drugs, identify novel ligands and guide flavor design. The novel features capture information regarding the molecules and their receptors, which could inform various chemoinformatic tasks. Inclusion of neighbor-informed similarities improves as experimental data mounts, and provides a generalizable framework for molecule-biotarget matching.

Functional molecular switches of mammalian G protein-coupled bitter-taste receptors

Bitter taste receptors (TAS2Rs) are a poorly understood subgroup of G protein-coupled receptors (GPCRs). The experimental structure of these receptors has yet to be determined, and key-residues controlling their function remain mostly unknown. We designed an integrative approach to improve comparative modeling of TAS2Rs. Using current knowledge on class A GPCRs and existing experimental data in the literature as constraints, we pinpointed conserved motifs to entirely re-align the amino-acid sequences of TAS2Rs. We constructed accurate homology models of human TAS2Rs. As a test case, we examined the accuracy of the TAS2R16 model with site-directed mutagenesis and in vitro functional assays. This combination of in silico and in vitro results clarifies sequence-function relationships and proposes functional molecular switches that encode agonist sensing and downstream signaling mechanisms within mammalian TAS2Rs sequences.

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

Ligand Affinities within the Open-Boundary Molecular Mechanics/Coarse-Grained Framework (I): Alchemical Transformations within the Hamiltonian Adaptive Resolution Scheme

Our recently developed Open-Boundary Molecular Mechanics/Coarse Grained (OB-MM/CG) framework predicts ligand poses in important pharmaceutical targets, such as G-protein Coupled Receptors, even when experimental structural information is lacking. The approach, which is based on GROMOS and AMBER force fields, allows for grand-canonical simulations of protein-ligand complexes by using the Hamiltonian Adaptive Resolution Scheme (H-AdResS) for the solvent. Here, we present a key step toward the estimation of ligand binding affinities for their targets within this approach. This is the implementation of the H-AdResS in the GROMACS code. The accuracy of our implementation is established by calculating hydration free energies of several molecules in water by means of alchemical transformations. The deviations of the GROMOS- and AMBER-based H-AdResS results from the reference fully atomistic simulations are smaller than the accuracy of the force field and/or they are in the range of the published results. Importantly, our predictions are in good agreement with experimental data. The current implementation paves the way to the use of the OB-MM/CG framework for the study of large biological systems.

In Silico Investigation of Bitter Hop-Derived Compounds and Their Cognate Bitter Taste Receptors

The typical bitter taste of beer is caused by adding hops (Humulus lupulus L.) during the wort boiling process. The bitter taste of hop-derived compounds was found to be mediated by three bitter taste receptors: TAS2R1, TAS2R14, and TAS2R40. In this work, structural bioinformatics analyses were used to characterize the binding modes of trans-isocohumulone, trans-isohumulone, trans-isoadhumulone, cis-isocohumulone, cis-isohumulone, cis-isoadhumulone, cohumulone, humulone, adhumulone, and 8-prenylnaringenin into the orthosteric binding site of their cognate receptors. A conserved asparagine in transmembrane 3 was found to be essential for the recognition of hop-derived compounds, whereas the surrounding residues in the binding site of the three receptors encode the ligand specificity. Hop-derived compounds are renowned bioactive molecules and are considered as potential hit molecules for drug discovery to treat metabolic diseases. A chemoinformatics analysis revealed that hop-derived compounds cluster in a different region of the chemical space compared to known bitter food-derived compounds, pinpointing hop-derived compounds as a very peculiar class of bitter compounds.

Hybrid MM/CG Webserver: Automatic Set Up of Molecular Mechanics/Coarse-Grained Simulations for Human G Protein-Coupled Receptor/Ligand Complexes

Hybrid Molecular Mechanics/Coarse-Grained (MM/CG) simulations help predict ligand poses in human G protein-coupled receptors (hGPCRs), the most important protein superfamily for pharmacological applications. This approach allows the description of the ligand, the binding cavity, and the surrounding water molecules at atomistic resolution, while coarse-graining the rest of the receptor. Here, we present the Hybrid MM/CG Webserver (mmcg.grs.kfa-juelich.de) that automatizes and speeds up the MM/CG simulation setup of hGPCR/ligand complexes. Initial structures for such complexes can be easily and efficiently generated with other webservers. The Hybrid MM/CG server also allows for equilibration of the systems, either fully automatically or interactively. The results are visualized online (using both interactive 3D visualizations and analysis plots), helping the user identify possible issues and modify the setup parameters accordingly. Furthermore, the prepared system can be downloaded and the simulation continued locally.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

Understanding and Controlling Food Protein Structure and Function in Foods: Perspectives from Experiments and Computer Simulations

The structure and interactions of proteins play a critical role in determining the quality attributes of many foods, beverages, and pharmaceutical products. Incorporating a multiscale understanding of the structure-function relationships of proteins can provide greater insight into, and control of, the relevant processes at play. Combining data from experimental measurements, human sensory panels, and computer simulations through machine learning allows the construction of statistical models relating nanoscale properties of proteins to the physicochemical properties, physiological outcomes, and tastes of foods. This review highlights several examples of advanced computer simulations at molecular, mesoscale, and multiscale levels that shed light on the mechanisms at play in foods, thereby facilitating their control. It includes a practical simulation toolbox for those new to in silico modeling.

Highly conserved intracellular H208 residue influences agonist selectivity in bitter taste receptor T2R14

Bitter taste receptors (T2Rs) are a specialized class of cell membrane receptors of the G protein-coupled receptor family and perform a crucial role in chemosensation. The 25 T2Rs in humans are activated by structurally diverse ligands of plant, animal and microbial origin. The mechanisms of activation of these receptors are poorly understood. Therefore, identification of structural determinants of T2Rs that regulate its efficacy could be beneficial in understanding the molecular mechanisms of T2R activation. In this work, we characterized a highly conserved histidine (H208), present at TM5-ICL3 region of T2R14 and its role in agonist-induced T2R14 signaling. Surprisingly, mutation of the conserved H208 (H208A) did not result in increased basal activity of T2R14, in contrast to similar H206A mutation in T2R4 that showed constitutive or basal activity. However, H208A mutation in T2R14 resulted in an increase in agonist-induced efficacy for Flufenamic acid (FFA). Interestingly, H208A did not affect the potency of another T2R14 agonist Diphenhydramine (DPH). The H208R compensatory mutation showed FFA response similar to wild-type T2R14. Molecular modeling suggests a FFA-induced shift in TM3 and TM5 helices of H208A, which changes the network of interactions connecting TM5-ICL3-TM6. This report identifies a crucial residue on the intracellular surface of T2Rs that is involved in bitter ligand selectivity. It also highlights the varied roles carried out by some conserved residues in different T2Rs.

Toward the Identification of Extra-Oral TAS2R Agonists as Drug Agents for Muscle Relaxation Therapies via Bioinformatics-Aided Screening of Bitter Compounds in Traditional Chinese Medicine

Significant advances have been made in the past decade in mapping the distributions and the physiological functions of extra-oral bitter taste receptors (TAS2Rs) in non-gustatory tissues. In particular, it has been found that TAS2Rs are expressed in various muscle tissues and activation of TAS2Rs can lead to muscle cell relaxation, which suggests that TAS2Rs may be important new targets in muscle relaxation therapy for various muscle-related diseases. So far, however, there is a lack of potent extra-oral TAS2R agonists that can be used as novel drug agents in muscle relaxation therapies. Interestingly, traditional Chinese medicine (TCM) often characterizes a drug’s property in terms of five distinct flavors (bitter, sweet, sour, salty, and pungent) according to its taste and function, and commonly regards “bitterness” as an intrinsic property of “good medicine.” In addition, many bitter flavored TCM are known in practice to cause muscle relaxation after long term use, and in lab experiments the compounds identified from some bitter flavored TCM do activate TAS2Rs and thus relax muscle cells. Therefore, it is highly possible to discover very useful extra-oral TAS2R agonists for muscle relaxation therapies among the abundant bitter compounds used in bitter flavored TCM. With this perspective, we reviewed in literature the distribution of TAS2Rs in different muscle systems with a focus on the map of bitter flavored TCM which can regulate muscle contractility and related functional chemical components. We also reviewed the recently established databases of TCM chemical components and the bioinformatics software which can be used for high-throughput screening and data mining of the chemical components associated with bitter flavored TCM. All together, we aim to present a knowledge-based approach and technological platform for identification or discovery of extra-oral TAS2R agonists that can be used as novel drug agents for muscle relaxation therapies through screening and evaluation of chemical compounds used in bitter flavored TCM.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

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.

Multiscale simulations on human Frizzled and Taste2 GPCRs

Recently, molecular dynamics simulations, from all atom and coarse grained to hybrid methods bridging the two scales, have provided exciting functional insights into class F (Frizzled and Taste2) GPCRs (about 40 members in humans). Findings include: (i) The activation of one member of the Frizzled receptors (FZD4) involves a bending of transmembrane helix TM7 far larger than that in class A GPCRs. (ii) The affinity of an anticancer drug targeting another member (Smoothened receptor) decreases in a specific drug-resistant variant, because the mutation ultimately disrupts the binding cavity and affects TM6. (iii) A novel two-state recognition mechanism explains the very large agonist diversity for at least one member of the Taste2 GPCRs, hTAS2R46.

Beyond the Flavour: The Potential Druggability of Chemosensory G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) belong to the largest class of drug targets. Approximately half of the members of the human GPCR superfamily are chemosensory receptors, including odorant receptors (ORs), trace amine-associated receptors (TAARs), bitter taste receptors (TAS2Rs), sweet and umami taste receptors (TAS1Rs). Interestingly, these chemosensory GPCRs (csGPCRs) are expressed in several tissues of the body where they are supposed to play a role in biological functions other than chemosensation. Despite their abundance and physiological/pathological relevance, the druggability of csGPCRs has been suggested but not fully characterized. Here, we aim to explore the potential of targeting csGPCRs to treat diseases by reviewing the current knowledge of csGPCRs expressed throughout the body and by analysing the chemical space and the drug-likeness of flavour molecules.

Open-Boundary Molecular Mechanics/Coarse-Grained Framework for Simulations of Low-Resolution G-Protein-Coupled Receptor-Ligand Complexes

G-protein-coupled receptors (GPCRs) constitute as much as 30% of the overall proteins targeted by FDA-approved drugs. However, paucity of structural experimental information and low sequence identity between members of the family impair the reliability of traditional docking approaches and atomistic molecular dynamics simulations for in silico pharmacological applications. We present here a dual-resolution approach tailored for such low-resolution models. It couples a hybrid molecular mechanics/coarse-grained (MM/CG) scheme, previously developed by us for GPCR-ligand complexes, with a Hamiltonian-based adaptive resolution scheme (H-AdResS) for the solvent. This dual-resolution approach removes potentially inaccurate atomistic details from the model while building a rigorous statistical ensemble-the grand canonical one-in the high-resolution region. We validate the method on a well-studied GPCR-ligand complex, for which the 3D structure is known, against atomistic simulations. This implementation paves the way for future accurate in silico studies of low-resolution ligand/GPCRs models.

Allosteric sodium binding cavity in GPR3: a novel player in modulation of Aβ production

The orphan G-protein coupled receptor 3 (GPR3) belongs to class A G-protein coupled receptors (GPCRs) and is highly expressed in central nervous system neurons. Among other functions, it is likely associated with neuron differentiation and maturation. Recently, GPR3 has also been linked to the production of Aβ peptides in neurons. Unfortunately, the lack of experimental structural information for this receptor hampers a deep characterization of its function. Here, using an in-silico and in-vitro combined approach, we describe, for the first time, structural characteristics of GPR3 receptor underlying its function: the agonist binding site and the allosteric sodium binding cavity. We identified and validated by alanine-scanning mutagenesis the role of three functionally relevant residues: Cys267^6.55, Phe120^3.36 and Asp^2.50. The latter, when mutated into alanine, completely abolished the constitutive and agonist-stimulated adenylate cyclase activity of GPR3 receptor by disrupting its sodium binding cavity. Interestingly, this is correlated with a decrease in Aβ production in a model cell line. Taken together, these results suggest an important role of the allosteric sodium binding site for GPR3 activity and open a possible avenue for the modulation of Aβ production in the Alzheimer’s Disease.

Reengineering the ligand sensitivity of the broadly tuned human bitter taste receptor TAS2R14

BACKGROUND

In humans, bitterness perception is mediated by ~25 bitter taste receptors present in the oral cavity. Among these receptors three, TAS2R10, TAS2R14 and TAS2R46, exhibit extraordinary wide agonist profiles and hence contribute disproportionally high to the perception of bitterness. Perhaps the most broadly tuned receptor is the TAS2R14, which may represent, because of its prominent expression in extraoral tissues, a receptor of particular importance for the physiological actions of bitter compounds beyond taste.

METHODS

To investigate how the architecture and composition of the TAS2R14 binding pocket enables specific interactions with a complex array of chemically diverse bitter agonists, we carried out homology modeling and ligand docking experiments, subjected the receptor to point-mutagenesis of binding site residues and performed functional calcium mobilization assays.

RESULTS

In total, 40 point-mutated receptor constructs were generated to investigate the contribution of 19 positions presumably located in the receptor's binding pocket to activation by 7 different TAS2R14 agonists. All investigated positions exhibited moderate to pronounced agonist selectivity.

CONCLUSIONS

Since numerous modifications of the TAS2R14 binding pocket resulted in improved responses to individual agonists, we conclude that this bitter taste receptor might represent a suitable template for the engineering of the agonist profile of a chemoreceptive receptor.

GENERAL SIGNIFICANCE

The detailed structure-function analysis of the highly promiscuous and widely expressed TAS2R14 suggests that this receptor must be considered as potentially frequent target for known and novel drugs including undesired off-effects.

Orthosteric and benzodiazepine cavities of the α1β2γ2 GABAA receptor: insights from experimentally validated in silico methods

γ-aminobutyric acid-type A (GABA_A) receptors mediate fast synaptic inhibition in the central nervous system of mammals. They are modulated via several sites by numerous compounds, which include GABA, benzodiazepines, ethanol, neurosteroids and anaesthetics among others. Due to their potential as targets of novel drugs, a detailed knowledge of their structure-function relationships is needed. Here, we present the model of the α₁β₂γ₂ subtype GABA_A receptor in the APO state and in complex with selected ligands, including agonists, antagonists and allosteric modulators. The model is based on the crystallographic structure of the human β₃ homopentamer GABA_A receptor. The complexes were refined using atomistic molecular dynamics simulations. This allowed a broad description of the binding modes and the detection of important interactions in agreement with experimental information. From the best of our knowledge, this is the only model of the α₁β₂γ₂ GABA_A receptor that represents altogether the desensitized state of the channel and comprehensively describes the interactions of ligands of the orthosteric and benzodiazepines binding sites in agreement with the available experimental data. Furthermore, it is able to explain small differences regarding the binding of a variety of chemically divergent ligands. Finally, this new model may pave the way for the design of focused experimental studies that will allow a deeper description of the receptor.

Independent Evolution of Strychnine Recognition by Bitter Taste Receptor Subtypes

The 25 human bitter taste receptors (hT2Rs) recognize thousands of structurally and chemically diverse bitter substances. The binding modes of human bitter taste receptors hT2R10 and hT2R46, which are responsible for strychnine recognition, were previously established using site-directed mutagenesis, functional assays, and molecular modeling. Here we construct a phylogenetic tree and reconstruct ancestral sequences of the T2R10 and T2R46 clades. We next analyze the binding sites in view of experimental data to predict their ability to recognize strychnine. This analysis suggests that the common ancestor of hT2R10 and hT2R46 is unlikely to bind strychnine in the same mode as either of its two descendants. Estimation of relative divergence times shows that hT2R10 evolved earlier than hT2R46. Strychnine recognition was likely acquired first by the earliest common ancestor of the T2R10 clade before the separation of primates from other mammals, and was highly conserved within the clade. It was probably independently acquired by the common ancestor of T2R43-47 before the homo-ape speciation, lost in most T2Rs within this clade, but enhanced in the hT2R46 after humans diverged from the rest of primates. Our findings suggest hypothetical strychnine T2R receptors in several species, and serve as an experimental guide for further study. Improved understanding of how bitter taste receptors acquire the ability to be activated by particular ligands is valuable for the development of sensors for bitterness and for potential toxicity.

Plasticity of the ligand binding pocket in the bitter taste receptor T2R7

Bitter taste receptors (T2Rs) are a group of 25 G protein-coupled receptors (GPCRs) in humans. The cognate agonists and the mechanism of ligand binding to the majority of the T2Rs remain unknown. Here we report the first structure-function analysis of T2R7 and study the ability of this receptor to bind to different agonists by site-directed mutagenesis. Screening of ligands for T2R7 in calcium based assays lead to the identification of novel compounds that activate this receptor. Quinine, diphenidol, dextromethorphan and diphenhydramine showed substantial activation of T2R7. Interestingly, these bitter compounds showed different pharmacological characteristics. To investigate the structural features in T2R7 that might contribute to the observed differences in agonist specificities, molecular model guided ligand docking and site-directed mutagenesis was pursued. Amino acids D65, D86, W89, N167, T169, W170, S181, T255 and E271 in the ligand-binding pocket were replaced and the mutants characterized pharmacologically. Our results suggest D86, S181 and W170 present on the extracellular side of transmembrane 3 (TM3), TM5 and in extracellular loop 2 (ECL2) are essential for agonist binding in T2R7. Mutations of these amino acids lead to loss-of-function. We also identified gain-of-function residues that are agonist specific. These results suggest that agonists bind at an extracellular site rather than deep within the TM core involving residues present in both ECL2 and TM helices in T2R7. Similar to majority of the Class A GPCRs, ECL2 in T2R7 plays a significant role in agonist binding and activation.

Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme

The recently proposed Hamiltonian adaptive resolution scheme (H-AdResS) allows the performance of molecular simulations in an open boundary framework. It allows changing, on the fly, the resolution of specific subsets of molecules (usually the solvent), which are free to diffuse between the atomistic region and the coarse-grained reservoir. So far, the method has been successfully applied to pure liquids. Coupling the H-AdResS methodology to hybrid models of proteins, such as the molecular mechanics/coarse-grained (MM/CG) scheme, is a promising approach for rigorous calculations of ligand binding free energies in low-resolution protein models. Toward this goal, here we apply for the first time H-AdResS to two atomistic proteins in dual-resolution solvent, proving its ability to reproduce structural and dynamic properties of both the proteins and the solvent, as obtained from atomistic simulations.