Pungent chemical compounds originating from decaying tissue are strong drivers of animal behavior. Two of the best-characterized death smell components are putrescine (PUT) and cadaverine (CAD), foul-smelling molecules produced by decarboxylation of amino acids during decomposition. These volatile polyamines act as 'necromones', triggering avoidance or attractive responses, which are fundamental for the survival of a wide range of species. The few studies that have attempted to identify the cognate receptors for these molecules have suggested the involvement of the seven-helix trace amine-associated receptors (TAARs), localized in the olfactory epithelium. However, very little is known about the precise chemosensory receptors that sense these compounds in the majority of organisms and the molecular basis of their interactions. In this work, we have used computational strategies to characterize the binding between PUT and CAD with the TAAR6 and TAAR8 human receptors. Sequence analysis, homology modeling, docking and molecular dynamics studies suggest a tandem of negatively charged aspartates in the binding pocket of these receptors which are likely to be involved in the recognition of these small biogenic diamines.

BACKGROUND: Protein sequence alignments and database search methods use standard scoring matrices calculated from amino acid substitution frequencies in general sets of proteins. These general-purpose matrices are not optimal to align accurately sequences with marked compositional biases, such as hydrophobic transmembrane regions found in membrane proteins. In this work, an amino acid substitution matrix (GPCRtm) is calculated for the membrane spanning segments of the G protein-coupled receptor (GPCR) rhodopsin family; one of the largest transmembrane protein family in humans with great importance in health and disease. RESULTS: The GPCRtm matrix reveals the amino acid compositional bias distinctive of the GPCR rhodopsin family and differs from other standard substitution matrices. These membrane receptors, as expected, are characterized by a high content of hydrophobic residues with regard to globular proteins. On the other hand, the presence of polar and charged residues is higher than in average membrane proteins, displaying high frequencies of replacement within themselves. CONCLUSIONS: Analysis of amino acid frequencies and values obtained from the GPCRtm matrix reveals patterns of residue replacements different from other standard substitution matrices. GPCRs prioritize the reactivity properties of the amino acids over their bulkiness in the transmembrane regions. A distinctive role is that charged and polar residues seem to evolve at different rates than other amino acids. This observation is related to the role of the transmembrane bundle in the binding of ligands, that in many cases involve electrostatic and hydrogen bond interactions. This new matrix can be useful in database search and for the construction of more accurate sequence alignments of GPCRs.

MOTIVATION: Integral polytopic membrane proteins contain only two types of folds in their transmembrane domains: α-helix bundles and β-barrels. The increasing number of available crystal structures of these proteins permits an initial estimation of how sequence variability affects the structure conservation in their transmembrane domains. We, thus, aim to determine the pairwise sequence identity necessary to maintain the transmembrane molecular architectures compatible with the hydrophobic nature of the lipid bilayer. RESULTS: Root-mean-square deviation (rmsd) and sequence identity were calculated from the structural alignments of pairs of homologous polytopic membrane proteins sharing the same fold. Analysis of these data reveals that transmembrane segment pairs with sequence identity in the so-called 'twilight zone' (20-35%) display high-structural similarity (rmsd < 1.5 Å). Moreover, a large group of β-barrel pairs with low-sequence identity (<20%) still maintain a close structural similarity (rmsd < 2.5 Å). Thus, we conclude that fold preservation in transmembrane regions requires less sequence conservation than for globular proteins. These findings have direct implications in homology modeling of evolutionary-related membrane proteins. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.