Computational studies of ligand-receptor interactions in bitter taste receptors

Phenylthiocarbamide tastes intensely bitter to some individuals, but others find it completely tasteless. Recently, it was suggested that phenylthiocarbamide elicits bitter taste by interacting with a human G protein-coupled receptor (hTAS2R38) encoded by the PTC gene. The phenylthiocarbamide nontaster trait was linked to three single nucleotide polymorphisms occurring in the PTC gene. Using the crystal structure of bovine rhodopsin as template, we generated the 3D structure of hTAS2R38 bitter taste receptor. We were able to map on the receptor structure the amino acids affected by the genetic polymorphisms and to propose molecular functions for two of them that explained the emergence of the nontaster trait. We used molecular docking simulations to find that phenylthiocarbamide exhibited a higher affinity for the target receptor than the structurally similar molecule 6-n-propylthiouracil, in line with recent experimental studies. A 3D model was constructed for the hTAS2R16 bitter taste receptor as well, by applying the same protocol. We found that the recently published experimental ligand binding affinity data for this receptor correlated well with the binding scores obtained from our molecular docking calculations.     

Insights into the Binding of Phenyltiocarbamide (PTC) Agonist to Its Target Human TAS2R38 Bitter Receptor

Humans'' bitter taste perception is mediated by the hTAS2R subfamily of the G protein-coupled membrane receptors (GPCRs). Structural information on these receptors is currently limited. Here we identify residues involved in the binding of phenylthiocarbamide (PTC) and in receptor activation in one of the most widely studied hTAS2Rs (hTAS2R38) by means of structural bioinformatics and molecular docking. The predictions are validated by site-directed mutagenesis experiments that involve specific residues located in the putative binding site and trans-membrane (TM) helices 6 and 7 putatively involved in receptor activation. Based on our measurements, we suggest that (i) residue N103 participates actively in PTC binding, in line with previous computational studies. (ii) W99, M100 and S259 contribute to define the size and shape of the binding cavity. (iii) W99 and M100, along with F255 and V296, play a key role for receptor activation, providing insights on bitter taste receptor activation not emerging from the previously reported computational models.     

The molecular basis of individual differences in phenylthiocarbamide and propylthiouracil bitterness perception

Individual differences in perception are ubiquitous within the chemical senses: taste, smell, and chemical somesthesis . A hypothesis of this fact states that polymorphisms in human sensory receptor genes could alter perception by coding for functionally distinct receptor types . We have previously reported evidence that sequence variants in a presumptive bitter receptor gene (hTAS2R38) correlate with differences in bitterness recognition of phenylthiocarbamide (PTC) . Here, we map individual psychogenomic pathways for bitter taste by testing people with a variety of psychophysical tasks and linking their individual perceptions of the compounds PTC and propylthiouracil (PROP) to the in vitro responses of their TAS2R38 receptor variants. Functional expression studies demonstrate that five different haplotypes from the hTAS2R38 gene code for operatively distinct receptors. The responses of the three haplotypes we also tested in vivo correlate strongly with individuals' psychophysical bitter sensitivities to a family of compounds. These data provide a direct molecular link between heritable variability in bitter taste perception to functional variations of a single G protein coupled receptor that responds to compounds such as PTC and PROP that contain the N-C=S moiety. The molecular mechanisms of perceived bitterness variability have therapeutic implications, such as helping patients to consume beneficial bitter-tasting compounds-for example, pharmaceuticals and selected phytochemicals.     

Probenecid inhibits the human bitter taste receptor TAS2R16 and suppresses bitter perception of salicin

Bitter taste stimuli are detected by a diverse family of G protein-coupled receptors (GPCRs) expressed in gustatory cells. Each bitter taste receptor (TAS2R) responds to an array of compounds, many of which are toxic and can be found in nature. For example, human TAS2R16 (hTAS2R16) responds to β-glucosides such as salicin, and hTAS2R38 responds to thiourea-containing molecules such as glucosinolates and phenylthiocarbamide (PTC). While many substances are known to activate TAS2Rs, only one inhibitor that specifically blocks bitter receptor activation has been described. Here, we describe a new inhibitor of bitter taste receptors, p-(dipropylsulfamoyl)benzoic acid (probenecid), that acts on a subset of TAS2Rs and inhibits through a novel, allosteric mechanism of action. Probenecid is an FDA-approved inhibitor of the Multidrug Resistance Protein 1 (MRP1) transporter and is clinically used to treat gout in humans. Probenecid is also commonly used to enhance cellular signals in GPCR calcium mobilization assays. We show that probenecid specifically inhibits the cellular response mediated by the bitter taste receptor hTAS2R16 and provide molecular and pharmacological evidence for direct interaction with this GPCR using a non-competitive (allosteric) mechanism. Through a comprehensive analysis of hTAS2R16 point mutants, we define amino acid residues involved in the probenecid interaction that result in decreased sensitivity to probenecid while maintaining normal responses to salicin. Probenecid inhibits hTAS2R16, hTAS2R38, and hTAS2R43, but does not inhibit the bitter receptor hTAS2R31 or non-TAS2R GPCRs. Additionally, structurally unrelated MRP1 inhibitors, such as indomethacin, fail to inhibit hTAS2R16 function. Finally, we demonstrate that the inhibitory activity of probenecid in cellular experiments translates to inhibition of bitter taste perception of salicin in humans. This work identifies probenecid as a pharmacological tool for understanding the cell biology of bitter taste and as a lead for the development of broad specificity bitter blockers to improve nutrition and medical compliance.     

Characterization of the ��-d-Glucopyranoside Binding Site of the Human Bitter Taste Receptor hTAS2R16

G-protein-coupled receptors mediate the senses of taste, smell, and vision in mammals. Humans recognize thousands of compounds as bitter, and this response is mediated by the hTAS2R family, which is one of the G-protein-coupled receptors composed of only 25 receptors. However, structural information on these receptors is limited. To address the molecular basis of bitter tastant discrimination by the hTAS2Rs, we performed ligand docking simulation and functional analysis using a series of point mutants of hTAS2R16 to identify its binding sites. The docking simulation predicted two candidate binding structures for a salicin-hTAS2R16 complex, and at least seven amino acid residues in transmembrane 3 (TM3), TM5, and TM6 were shown to be involved in ligand recognition. We also identified the probable salicin-hTAS2R16 binding mode using a mutated receptor experiment. This study characterizes the molecular interaction between hTAS2R16 and β-d-glucopyranoside and will also facilitate rational design of bitter blockers.     

Modeling the human PTC bitter-taste receptor interactions with bitter tastants

We employed the first principles computational method MembStruk and homology modeling techniques to predict the 3D structures of the human phenylthiocarbamide (PTC) taste receptor. This protein is a seven-transmembrane-domain G protein-coupled receptor that exists in two main forms worldwide, designated taster and nontaster, which differ from each other at three amino-acid positions. 3D models were generated with and without structural similarity comparison to bovine rhodopsin. We used computational tools (HierDock and ScanBindSite) to generate models of the receptor bound to PTC ligand to estimate binding sites and binding energies. In these models, PTC binds at a site distant from the variant amino acids, and PTC binding energy was equivalent for both the taster and nontaster forms of the protein. These models suggest that the inability of humans to taste PTC is due to a failure of G protein activation rather than decreased binding affinity of the receptor for PTC. Amino-acid substitutions in the sixth and seventh transmembrane domains of the nontaster form of the protein may produce increased steric hindrance between these two α-helices and reduce the motion of the sixth helix required for G protein activation.     

Genetic variation in the TAS2R38 taste receptor contributes to the oral microbiota in North and South European locations: a pilot study

Background

Microbial communities are influenced by environmental factors including host genetics. We investigated the relationship between host bitter taste receptor genotype hTAS2R38 and oral microbiota, together with the influence of geographical location.

Methods

hTAS2R38 polymorphisms and 16S bacterial gene sequencing from oral samples were analyzed from a total of 45 healthy volunteers from different geographical locations.

Results

Genetic variation in the bitter taste receptor TAS2R38 reflected in the microbial composition of oral mucosa in Finnish and Spanish subjects. Multivariate analysis showed significant differences in the microbial composition between country and also dependent on taste genotype. Oral microbiota was shown to be more stable to the geographical location impact among AVI-homozygotes than PAV-homozygotes or heterozygotes (PAV/AVI).

Conclusion

Geographical location and genetic variation in the hTAS2R38 taste receptor impact oral mucosa microbial composition. These findings provide an advance in the knowledge regarding the interactions between taste receptor genes and oral microbiota. This study suggests the role of host-microbiota interactions on the food taste perception in food choices, nutrition, and eating behavior.

     

Structure-function relationships of the human bitter taste receptor hTAS2R1: insights from molecular modeling studies

Bitter taste receptors (T2Rs) belong to G-protein-coupled receptors (GPCRs). Despite extensive studies, the precise mechanisms of GPCR activation are still poorly understood. In this study, the models of the human bitter taste receptor hTAS2R1 alone and in complex with various ligands were constructed on the basis of template-based modeling and molecular docking. Then these models were subjected to all-atom molecular dynamics (MD) simulations in explicit lipid bilayers. The binding pocket of hTAS2R1 is mainly formed by transmembrane helix (TM) III, TM V, TM VI, and TM VII. Most of the residues contributing to ligand binding are positionally conserved comparing with other hTAS2Rs. By comparing the final conformations obtained by extensive MD simulations, we identified the changes in the transmembrane helices and the intra- and extracellular loops, which were expected to initiate the activation of the receptor. The intracellular loop II (ICL2) and TM III were found to play prominent roles in the process of activation. We proposed that a set of interactions between the aromatic Phe115 in the middle of ICL2 and three residues (Tyr103, Lys106, and Val107) at the cytoplasmic end of TM III may serve as a conformational switch of hTAS2R1 activation. All of the residues involved in the switch are highly conserved among T2Rs. This indicates that the control switch we proposed may be universal in T2Rs. Besides, our results also suggest that the formation of a short helical segment in ICL2 may be necessary for the activation of hTAS2R1.     

The human bitter taste receptor, hTAS2R16, discriminates slight differences in the configuration of disaccharides

Sweetness and bitterness are key determinants of food acceptance and rejection, respectively. Sugars, such as sucrose and fructose, are generally recognized as sweet. However, not all sugars are sweet, and even anomers may have quite different tastes. For example, gentiobiose is bitter, whereas its anomer, isomaltose, is sweet. Despite this unique sensory character, the molecular basis of the bitterness of gentiobiose remains to be clarified. In this study, we used calcium imaging analysis of human embryonic kidney 293T cells that heterologously expressed human taste receptors to demonstrate that gentiobiose activated hTAS2R16, a bitter taste receptor, but not hT1R2/hT1R3, a sweet taste receptor. In contrast, isomaltose activated hT1R2/hT1R3. As a result, these anomers elicit different taste sensations. Mutational analysis of hTAS2R16 also indicated that gentiobiose and β-d-glucopyranosides, such as salicin share a common binding site of hTAS2R16.     

Structure–function relationships of the human bitter taste receptor hTAS2R1: insights from molecular modeling studies

Bitter taste receptors (T2Rs) belong to G-protein-coupled receptors (GPCRs). Despite extensive studies, the precise mechanisms of GPCR activation are still poorly understood. In this study, the models of the human bitter taste receptor hTAS2R1 alone and in complex with various ligands were constructed on the basis of template-based modeling and molecular docking. Then these models were subjected to all-atom molecular dynamics (MD) simulations in explicit lipid bilayers. The binding pocket of hTAS2R1 is mainly formed by transmembrane helix (TM) III, TM V, TM VI, and TM VII. Most of the residues contributing to ligand binding are positionally conserved comparing with other hTAS2Rs. By comparing the final conformations obtained by extensive MD simulations, we identified the changes in the transmembrane helices and the intra- and extracellular loops, which were expected to initiate the activation of the receptor. The intracellular loop II (ICL2) and TM III were found to play prominent roles in the process of activation. We proposed that a set of interactions between the aromatic Phe115 in the middle of ICL2 and three residues (Tyr103, Lys106, and Val107) at the cytoplasmic end of TM III may serve as a conformational switch of hTAS2R1 activation. All of the residues involved in the switch are highly conserved among T2Rs. This indicates that the control switch we proposed may be universal in T2Rs. Besides, our results also suggest that the formation of a short helical segment in ICL2 may be necessary for the activation of hTAS2R1.     

Sex Differences in the Effects of Inherited Bitter Thiourea Sensitivity on Body Weight in 4–6‐Year‐Old Children

Previous studies have shown that inherited taste blindness to bitter compounds like 6‐n‐propylthiouracil (PROP) may be a risk factor for obesity, but this literature has been highly controversial. The objectives of this study were (i) to confirm findings that show an interaction between PROP status and sex on BMI z‐score, and (ii) to determine if sex also interacts with variations in TAS2R38 (phenylthiocarbamide (PTC) genotype) to influence weight status in 4–6 year olds. Also, we tested whether nontaster children consumed more fat and total energy at laboratory‐based meals. Seventy‐two ethnically diverse children who ranged in weight status were classified as tasters (N = 52) or nontasters (N = 20) using a standard PROP screening solution. Anthropometric measures were taken, and at the end of each visit, children ate ad libitum from test meals intended for exploratory purposes. Genomic DNA was extracted from saliva and alleles at TAS2R38 were genotyped for A49P polymorphisms. In 75.8% of children, PTC genotype predicted PROP phenotype, whereas in 24.4%, genotype did not predict phenotype. PROP nontaster males had higher BMI z‐scores than taster‐males and females in both groups (P < 0.05), but due to a three‐way interaction between PROP phenotype, TAS2R38 genotype, and sex, this relationship was only true for children who were homozygous for the bitter‐insensitive allele (P < 0.0005). There were no differences in test‐meal intake as a function of PROP phenotype or TAS2R38 genotype. These results suggest that the TAS2R38 variation, PROP phenotype, and sex interact to impact obesity risk in children. Future studies should be done to determine how this trait influences energy balance.     

Supertasting and PROP bitterness depends on more than the TAS2R38 gene

Polymorphisms in the TAS2R38 gene provide insight to phenotypes long associated 6-n-propylthiouracil (PROP) and phenylthiocarbamide bitterness. We tested relationships between TAS2R38 genotype, taste phenotype, and fungiform papillae (FP) number in 139 females and 59 males (age range 21-60 years), primarily of European ancestry. DNA was analyzed for 3 polymorphic sites, identifying common (alanine-valine-isoleucine [AVI/AVI], heterozygotes, proline-alanine-valine [PAV/PAV]) and rare (proline-valine-isoleucine, alanine-alanine-valine, AAI) forms. Individuals with PROP threshold >0.15 mM were almost exclusively AVI/AVI; those with threshold <0.1 mM could have any genotype. PAV/PAVs were more difficult to identify with PROP taste measures, although perceived bitterness of moderate PROP concentrations (0.32, 1 mM) had better correspondence with genotype than did threshold. For AVI/AVIs, increases in bitterness from 1 to 3.2 mM PROP nearly paralleled those of TAS2R38 heterozygotes and PAV/PAVs. Some bitterness gains were related to FP number sampled from a standard area on the tongue tip, yet the PROP bitterness-FP relationship differed across genotype. Among homozygotes, FP was a significant determinant of PROP bitterness; heterozygotes showed a flat relationship. Those tasting concentrated PROP as more bitter also tasted concentrated sucrose, citric acid, sodium chloride, and quinine as more intense, even after statistically controlling for TAS2R38 genotype, FP, and intensity of tones (nonoral standard). To summarize, although PROP threshold generally exhibited single-gene complete dominance, PROP bitterness may involve additional bitter receptors as evidenced by misclassification of some nontaster homozygotes and the bitterness functions for concentrated PROP. Variability in receptor expression may explain attenuated bitterness-FP relationships. PROP bitterness does associate with heightened taste sensations (i.e., supertasting), but this is not due to TAS2R38 polymorphisms.     

Bitter peptides activate hTAS2Rs, the human bitter receptors

Fermented food contains numerous peptides derived from material proteins. Bitter peptides formed during the fermentation process are responsible for the bitter taste of fermented food. We investigated whether human bitter receptors (hTAS2Rs) recognize bitterness of peptides with a heterologous expression system. HEK293 cells expressing hTAS2R1, hTAS2R4, hTAS2R14, and hTAS2R16 responded to bitter casein digests. Among those cells, the hTAS2R1-expressing cell was most strongly activated by the synthesized bitter peptides Gly-Phe and Gly-Leu, and none of the cells was activated by the non-bitter dipeptide Gly-Gly. The results showed that these bitter peptides, as well as many other bitter compounds, activate hTAS2Rs, suggesting that humans utilize these hTAS2Rs to recognize and perceive the structure and bitterness of peptides.     

Activation of human bitter taste receptors by polymethoxylated flavonoids

Tangeretin and nobiletin are polymethoxylated flavonoids in citrus peel. Both tangeretin and nobiletin are bitter; however, their bitterness has not been evaluated using human bitter taste receptors (hTAS2Rs). We screened 25 kinds of hTAS2Rs and found that hTAS2R14 and hTAS2R46 received both compounds.     

Evolution of functionally diverse alleles associated with PTC bitter taste sensitivity in Africa

Although human bitter taste perception is hypothesized to be a dietary adaptation, little is known about genetic signatures of selection and patterns of bitter taste perception variability in ethnically diverse populations with different diets, particularly from Africa. To better understand the genetic basis and evolutionary history of bitter taste sensitivity, we sequenced a 2,975 bp region encompassing TAS2R38, a bitter taste receptor gene, in 611 Africans from 57 populations in West Central and East Africa with diverse subsistence patterns, as well as in a comparative sample of 132 non-Africans. We also examined the association between genetic variability at this locus and threshold levels of phenylthiocarbamide (PTC) bitterness in 463 Africans from the above populations to determine how variation influences bitter taste perception. Here, we report striking patterns of variation at TAS2R38, including a significant excess of novel rare nonsynonymous polymorphisms that recently arose only in Africa, high frequencies of haplotypes in Africa associated with intermediate bitter taste sensitivity, a remarkably similar frequency of common haplotypes across genetically and culturally distinct Africans, and an ancient coalescence time of common variation in global populations. Additionally, several of the rare nonsynonymous substitutions significantly modified levels of PTC bitter taste sensitivity in diverse Africans. While ancient balancing selection likely maintained common haplotype variation across global populations, we suggest that recent selection pressures may have also resulted in the unusually high level of rare nonsynonymous variants in Africa, implying a complex model of selection at the TAS2R38 locus in African populations. Furthermore, the distribution of common haplotypes in Africa is not correlated with diet, raising the possibility that common variation may be under selection due to their role in nondietary biological processes. In addition, our data indicate that novel rare mutations contribute to the phenotypic variance of PTC sensitivity, illustrating the influence of rare variation on a common trait, as well as the relatively recent evolution of functionally diverse alleles at this locus.     

Snooker Structure-Based Pharmacophore Model Explains Differences in Agonist and Blocker Binding to Bitter Receptor hTAS2R39

The human bitter taste receptor hTAS2R39 can be activated by many dietary (iso)flavonoids. Furthermore, hTAS2R39 activity can be blocked by 6-methoxyflavanones, 4’-fluoro-6-methoxyflavanone in particular. A structure-based pharmacophore model of the hTAS2R39 binding pocket was built using Snooker software, which has been used successfully before for drug design of GPCRs of the rhodopsin subfamily. For the validation of the model, two sets of compounds, both of which contained actives and inactives, were used: (i) an (iso)flavonoid-dedicated set, and (ii) a more generic, structurally diverse set. Agonists were characterized by their linear binding geometry and the fact that they bound deeply in the hTAS2R39 pocket, mapping the hydrogen donor feature based on T5.45 and N3.36, analogues of which have been proposed to play a key role in activation of GPCRs. Blockers lack hydrogen-bond donors enabling contact to the receptor. Furthermore, they had a crooked geometry, which could sterically hinder movement of the TM domains upon receptor activation. Our results reveal characteristics of hTAS2R39 agonist and bitter blocker binding, which might facilitate the development of blockers suitable to counter the bitterness of dietary hTAS2R39 agonists in food applications.     

Evaluation of the bitterness of green tea catechins by a cell-based assay with the human bitter taste receptor hTAS2R39

Catechins have a broad range of physiological functions and act as the main taste ingredient of green tea. Although catechins show a strong bitterness, the bitter taste receptor for catechins has not been fully understood. The objective of this study was to identify the receptor for the major green tea catechins such as (−)-epicatechin (EC), (−)-epicatechin gallate (ECg), (−)-epigallocatechin (EGC), and (−)-epigallocatechin gallate (EGCg). By the cell-based assay using cultured cells expressing human bitter taste receptor, a clear response of hTAS2R39-expressing cells was observed to 300 μM of either ECg or EGCg, which elicit a strong bitterness in humans. The response of hTAS2R39-expressing cells to ECg was the strongest among the tested catechins, followed by EGCg. Because the cellular response to EC and EGC is much weaker than those of ECg and EGCg, galloyl groups was strongly supposed to be involved in the bitter intensity. This finding is similar to the observations of taste intensity obtained from a human sensory study. Our results suggest the participation of hTAS2R39 in the detection of catechins in humans, indicating the possibility that bitterness of tea catechins can be evaluated by using cells expressing hTAS2R39.     

Genetic variation in the hTAS2R38 taste receptor and food consumption among Finnish adults

Genetic variation in bitter taste receptors, such as hTAS2R38, may affect food preferences and intake. The aim of the present study was to investigate the association between bitter taste receptor haplotypes and the consumption of vegetables, fruits, berries and sweet foods among an adult Finnish population. A cross-sectional design utilizing data from the Cardiovascular Risk in Young Finns cohort from 2007, which consisted of 1,903 men and women who were 30–45 years of age from five different regions in Finland, was employed. DNA was extracted from blood samples, and hTAS2R38 polymorphisms were determined based on three SNPs (rs713598, rs1726866 and rs10246939). Food consumption was assessed with a validated food frequency questionnaire. The prevalence of the bitter taste-sensitive (PAV/PAV) haplotype was 11.3 % and that of the insensitive (AVI/AVI) haplotype was 39.5 % among this Finnish population. PAV homozygotic women consumed fewer vegetables than did the AVI homozygotic women, 269 g/day (SD 131) versus 301 g/day (SD 187), respectively, p = 0.03 (multivariate ANOVA). Furthermore, the intake of sweet foods was higher among the PAV homozygotes of both genders. Fruit and berry consumption did not differ significantly between the haplotypes in either gender. Individuals perceive foods differently, and this may influence their patterns of food consumption. This study showed that the hTAS2R38 taste receptor gene variation was associated with vegetable and sweet food consumption among adults in a Finnish population.     

Human Bitter Taste Receptors hTAS2R8 and hTAS2R39 with Differential Functions to Recognize Bitter Peptides

The strong bitter peptide, Phe-Phe-Pro-Arg, activated cultured cells expressing either of the known human bitter taste receptors, hTAS2R8 and hTAS2R39. The partial structure of Pro-Arg activated hTAS2R39, but did not activate hTAS2R8. These receptors may not indiscriminately recognize bitter peptides, but have a differential function in their recognition.     

Human bitter taste receptors hTAS2R8 and hTAS2R39 with differential functions to recognize bitter peptides

The strong bitter peptide, Phe-Phe-Pro-Arg, activated cultured cells expressing either of the known human bitter taste receptors, hTAS2R8 and hTAS2R39. The partial structure of Pro-Arg activated hTAS2R39, but did not activate hTAS2R8. These receptors may not indiscriminately recognize bitter peptides, but have a differential function in their recognition.