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.     

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.     

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.     

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.     

Natural selection and molecular evolution in PTC, a bitter-taste receptor gene

The ability to taste phenylthiocarbamide (PTC) is a classic phenotype that has long been known to vary in human populations. This phenotype is of genetic, epidemiologic, and evolutionary interest because the ability to taste PTC is correlated with the ability to taste other bitter substances, many of which are toxic. Thus, variation in PTC perception may reflect variation in dietary preferences throughout human history and could correlate with susceptibility to diet-related diseases in modern populations. To test R. A. Fisher's long-standing hypothesis that variability in PTC perception has been maintained by balancing natural selection, we examined patterns of DNA sequence variation in the recently identified PTC gene, which accounts for up to 85% of phenotypic variance in the trait. We analyzed the entire coding region of PTC (1,002 bp) in a sample of 330 chromosomes collected from African (n=62), Asian (n=138), European (n=110), and North American (n=20) populations by use of new statistical tests for natural selection that take into account the potentially confounding effects of human population growth. Two intermediate-frequency haplotypes corresponding to "taster" and "nontaster" phenotypes were found. These haplotypes had similar frequencies across Africa, Asia, and Europe. Genetic differentiation between the continental population samples was low (FST=0.056) in comparison with estimates based on other genes. In addition, Tajima's D and Fu and Li's D and F statistics demonstrated a significant deviation from neutrality because of an excess of intermediate-frequency variants when human population growth was taken into account (P<.01). These results combine to suggest that balancing natural selection has acted to maintain "taster" and "nontaster" alleles at the PTC locus in humans.     

Specific alleles of bitter receptor genes influence human sensitivity to the bitterness of aloin and saccharin

Variation in human taste is a well-known phenomenon. However, little is known about the molecular basis for it. Bitter taste in humans is believed to be mediated by a family of 25 G protein-coupled receptors (hT2Rs, or TAS2Rs). Despite recent progress in the functional expression of hT2Rs in vitro, up until now, hT2R38, a receptor for phenylthiocarbamide (PTC), was the only gene directly linked to variations in human bitter taste. Here we report that polymorphism in two hT2R genes results in different receptor activities and different taste sensitivities to three bitter molecules. The hT2R43 gene allele, which encodes a protein with tryptophan in position 35, makes people very sensitive to the bitterness of the natural plant compounds aloin and aristolochic acid. People who do not possess this allele do not taste these compounds at low concentrations. The same hT2R43 gene allele makes people more sensitive to the bitterness of an artificial sweetener, saccharin. In addition, a closely related gene's (hT2R44's) allele also makes people more sensitive to the bitterness of saccharin. We also demonstrated that some people do not possess certain hT2R genes, contributing to taste variation between individuals. Our findings thus reveal new examples of variations in human taste and provide a molecular basis for them.     

Gut bitter taste receptor signalling induces ABCB1 through a mechanism involving CCK

T2Rs (bitter taste-sensing type 2 receptors) are expressed in the oral cavity to prevent ingestion of dietary toxins through taste avoidance. They are also expressed in other cell types, including gut enteroendocrine cells, where their physiological role is enigmatic. Previously, we proposed that T2R-dependent CCK (cholecystokinin) secretion from enteroendocrine cells limits absorption of dietary toxins, but an active mechanism was lacking. In the present study we show that T2R signalling activates ABCB1 (ATP-binding cassette B1) in intestinal cells through a CCK signalling mechanism. PTC (phenylthiocarbamide), an agonist for the T2R38 bitter receptor, increased ABCB1 expression in both intestinal cells and mouse intestine. PTC induction of ABCB1 was decreased by either T2R38 siRNA (small interfering RNA) or treatment with YM022, a gastrin receptor antagonist. Thus gut ABCB1 is regulated through signalling by CCK/gastrin released in response to PTC stimulation of T2R38 on enteroendocrine cells. We also show that PTC increases the efflux activity of ABCB1, suggesting that T2R signalling limits the absorption of bitter tasting/toxic substances through modulation of gut efflux membrane transporters.     

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.     

Multiplex minisequencing screening for PTC genotype associated with bitter taste perception

Sensitivity to phenylthiocarbamide (PTC) has a bimodal distribution pattern and the genotype of the TAS2R38 gene, which is composed of combinations of three coding single nucleotide polymorphisms (SNPs), p.A49P (c.145G>C), p.V262A (c.785T>C) and p.I296 V (c.886A>G), determines the ability or inability to taste PTC. In this study, we developed a tool for genotyping of these SNPs in the TAS2R38 gene using SNaPshot minisequencing and investigated the accuracy of the tool in 100 subjects who were genotyped by Sanger sequencing. The minor allele frequencies of the three SNPs were 0.39, and these genotypes corresponded to those determined by direct sequencing. In conclusion, we successfully developed a precise and rapid genetic tool for analysis of PTC genotype associated with bitter taste perception.     

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.     

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.

     

Bitter taste study in a sardinian genetic isolate supports the association of phenylthiocarbamide sensitivity to the TAS2R38 bitter receptor gene

Recently, a major locus on chromosome 7q was found in association with the taste sensitivity to phenylthiocarbamide (PTC) in humans. This region contains the TAS2R38 gene that encodes a member of the TAS2R bitter taste receptor family. Three SNPs within this gene demonstrated a strong association with taster status in Utah families and in an additional sample of 85 unrelated individuals. We studied a small isolated village in eastern Sardinia and carried out a genome-wide scan to map the genetic basis of PTC perception in this population. We performed both qualitative and quantitative PTC-taste linkage analysis. Qualitative analysis was carried out by defining a cut-off from the bimodal distribution of the trait and classifying subjects as tasters and non-tasters (75 and 25%, respectively). Linkage analysis on 131 subjects belonging to a unique large multi-generation pedigree comprising 239 subjects confirmed significant evidence for linkage at 7q35 also in our population. Haplotype analyses of the three SNPs inside the PTC gene allowed us to identify only two haplotypes that were associated with the non-taster phenotype (80% AVI homozygous) and to taster phenotype (40% PAV homozygous and 56% PAV/AVI heterozygous). Sex, age and haplotype effect explained 77.2 % of the total variance in PTC sensitivity.     

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.     

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.     

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.