Improving model construction of profile HMMs for remote homology detection through structural alignment

<p>Abstract</p> <p>Background</p> <p>Remote homology detection is a challenging problem in Bioinformatics. Arguably, profile Hidden Markov Models (pHMMs) are one of the most successful approaches in addressing this important problem. pHMM packages present a relatively small computational cost, and perform particularly well at recognizing remote homologies. This raises the question of whether structural alignments could impact the performance of pHMMs trained from proteins in the <it>Twilight Zone</it>, as structural alignments are often more accurate than sequence alignments at identifying motifs and functional residues. Next, we assess the impact of using structural alignments in pHMM performance.</p> <p>Results</p> <p>We used the SCOP database to perform our experiments. Structural alignments were obtained using the 3DCOFFEE and MAMMOTH-mult tools; sequence alignments were obtained using CLUSTALW, TCOFFEE, MAFFT and PROBCONS. We performed leave-one-family-out cross-validation over super-families. Performance was evaluated through ROC curves and paired two tailed t-test.</p> <p>Conclusion</p> <p>We observed that pHMMs derived from structural alignments performed significantly better than pHMMs derived from sequence alignment in low-identity regions, mainly below 20%. We believe this is because structural alignment tools are better at focusing on the important patterns that are more often conserved through evolution, resulting in higher quality pHMMs. On the other hand, sensitivity of these tools is still quite low for these low-identity regions. Our results suggest a number of possible directions for improvements in this area.</p

A Family of Chemoreceptors in Tribolium castaneum (Tenebrionidae: Coleoptera)

Chemoperception in invertebrates is mediated by a family of G-protein-coupled receptors (GPCR). To date nothing is known about the molecular mechanisms of chemoperception in coleopteran species. Recently the genome of Tribolium castaneum was sequenced for use as a model species for the Coleoptera. Using blast searches analyses of the T. castaneum genome with previously predicted amino acid sequences of insect chemoreceptor genes, a putative chemoreceptor family consisting of 62 gustatory receptors (Grs) and 26 olfactory receptors (Ors) was identified. The receptors have seven transmembrane domains (7TMs) and all belong to the GPCR receptor family. The expression of the T. castaneum chemoreceptor genes was investigated using quantification real- time RT-PCR and in situ whole mount RT-PCR analysis in the antennae, mouth parts, and prolegs of the adults and larvae. All of the predicted TcasGrs were expressed in the labium, maxillae, and prolegs of the adults but TcasGr13, 19, 28, 47, 62, 98, and 61 were not expressed in the prolegs. The TcasOrs were localized only in the antennae and not in any of the beetles gustatory organs with one exception; the TcasOr16 (like DmelOr83b), which was localized in the antennae, labium, and prolegs of the beetles. A group of six TcasGrs that presents a lineage with the sugar receptors subfamily in Drosophila melanogaster were localized in the lacinia of the Tribolium larvae. TcasGr1, 3, and 39, presented an ortholog to CO2 receptors in D. melanogaster and Anopheles gambiae was recorded. Low expression of almost all of the predicted chemoreceptor genes was observed in the head tissues that contain the brains and suboesophageal ganglion (SOG). These findings demonstrate the identification of a chemoreceptor family in Tribolium, which is evolutionarily related to other insect species

Amino Acid Residues Contributing to Function of the Heteromeric Insect Olfactory Receptor Complex

Olfactory receptors (Ors) convert chemical signals—the binding of odors and pheromones—to electrical signals through the depolarization of olfactory sensory neurons. Vertebrates Ors are G-protein-coupled receptors, stimulated by odors to produce intracellular second messengers that gate ion channels. Insect Ors are a heteromultimeric complex of unknown stoichiometry of two seven transmembrane domain proteins with no sequence similarity to and the opposite membrane topology of G-protein-coupled receptors. The functional insect Or comprises an odor- or pheromone-specific Or subunit and the Orco co-receptor, which is highly conserved in all insect species. The insect Or-Orco complex has been proposed to function as a novel type of ligand-gated nonselective cation channel possibly modulated by G-proteins. However, the Or-Orco proteins lack homology to any known family of ion channel and lack known functional domains. Therefore, the mechanisms by which odors activate the Or-Orco complex and how ions permeate this complex remain unknown. To begin to address the relationship between Or-Orco structure and function, we performed site-directed mutagenesis of all 83 conserved Glu, Asp, or Tyr residues in the silkmoth BmOr-1-Orco pheromone receptor complex and measured functional properties of mutant channels expressed in Xenopus oocytes. 13 of 83 mutations in BmOr-1 and BmOrco altered the reversal potential and rectification index of the BmOr-1-Orco complex. Three of the 13 amino acids (D299 and E356 in BmOr-1 and Y464 in BmOrco) altered both current-voltage relationships and K+ selectivity. We introduced the homologous Orco Y464 residue into Drosophila Orco in vivo, and observed variable effects on spontaneous and evoked action potentials in olfactory neurons that depended on the particular Or-Orco complex examined. Our results provide evidence that a subset of conserved Glu, Asp and Tyr residues in both subunits are essential for channel activity of the heteromeric insect Or-Orco complex

Computational Model of the Insect Pheromone Transduction Cascade

A biophysical model of receptor potential generation in the male moth olfactory receptor neuron is presented. It takes into account all pre-effector processes—the translocation of pheromone molecules from air to sensillum lymph, their deactivation and interaction with the receptors, and the G-protein and effector enzyme activation—and focuses on the main post-effector processes. These processes involve the production and degradation of second messengers (IP3 and DAG), the opening and closing of a series of ionic channels (IP3-gated Ca2+ channel, DAG-gated cationic channel, Ca2+-gated Cl− channel, and Ca2+- and voltage-gated K+ channel), and Ca2+ extrusion mechanisms. The whole network is regulated by modulators (protein kinase C and Ca2+-calmodulin) that exert feedback inhibition on the effector and channels. The evolution in time of these linked chemical species and currents and the resulting membrane potentials in response to single pulse stimulation of various intensities were simulated. The unknown parameter values were fitted by comparison to the amplitude and temporal characteristics (rising and falling times) of the experimentally measured receptor potential at various pheromone doses. The model obtained captures the main features of the dose–response curves: the wide dynamic range of six decades with the same amplitudes as the experimental data, the short rising time, and the long falling time. It also reproduces the second messenger kinetics. It suggests that the two main types of depolarizing ionic channels play different roles at low and high pheromone concentrations; the DAG-gated cationic channel plays the major role for depolarization at low concentrations, and the Ca2+-gated Cl− channel plays the major role for depolarization at middle and high concentrations. Several testable predictions are proposed, and future developments are discussed

The Caenorhabditis chemoreceptor gene families

Background: Chemoreceptor proteins mediate the first step in the transduction of environmental chemical stimuli, defining the breadth of detection and conferring stimulus specificity. Animal genomes contain families of genes encoding chemoreceptors that mediate taste, olfaction, and pheromone responses. The size and diversity of these families reflect the biology of chemoperception in specific species. Results: Based on manual curation and sequence comparisons among putative G-protein-coupled chemoreceptor genes in the nematode Caenorhabditis elegans, we identified approximately 1300 genes and 400 pseudogenes in the 19 largest gene families, most of which fall into larger superfamilies. In the related species C. briggsae and C. remanei, we identified most or all genes in each of the 19 families. For most families, C. elegans has the largest number of genes and C. briggsae the smallest number, suggesting changes in the importance of chemoperception among the species. Protein trees reveal family-specific and species-specific patterns of gene duplication and gene loss. The frequency of strict orthologs varies among the families, from just over 50% in two families to less than 5% in three families. Several families include large species-specific expansions, mostly in C. elegans and C. remanei. Conclusion: Chemoreceptor gene families in Caenorhabditis species are large and evolutionarily dynamic as a result of gene duplication and gene loss. These dynamics shape the chemoreceptor gene complements in Caenorhabditis species and define the receptor space available for chemosensory responses. To explain these patterns, we propose the gray pawn hypothesis: individual genes are of little significance, but the aggregate of a large number of diverse genes is required to cover a large phenotype space.JHT was supported by NIH grant RO1GM48700 and HMR by R01AI56081

Quantitative genetic variation in carbonic anhydrase isozymes from tissues of the pig-tailed macaque, Macaca nemestrina

Two isozymes of carbonic anhydrase (CA I and CA II) were quantified by a radio-immunoassay in 10 different tissues of the pig-tailed macaque. There were clearly differences in relative amounts of the two isozymes, indicating a differential regulation of these two different gene products. An inherited deficiency variant reduced red cell CA I and CA II 5000-fold and 2.7-fold, respectively. In nine other tissues, CA I was reduced from approximately twofold to 110-fold, and CA II was essentially unchanged. The CA I in deficient red cells was immunochemically and electrophoretically identical to common electrophoretic variants of CA I in the pig-tailed macaque and was enzymatically active.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/44174/1/10528_2004_Article_BF00485860.pd

Expression and Membrane Topology of Anopheles gambiae Odorant Receptors in Lepidopteran Insect Cells

A lepidopteran insect cell-based expression system has been employed to express three Anopheles gambiae odorant receptors (ORs), OR1 and OR2, which respond to components of human sweat, and OR7, the ortholog of Drosophila's OR83b, the heteromerization partner of all functional ORs in that system. With the aid of epitope tagging and specific antibodies, efficient expression of all ORs was demonstrated and intrinsic properties of the proteins were revealed. Moreover, analysis of the orientation of OR1 and OR2 on the cellular plasma membrane through the use of a novel ‘topology screen’ assay and FACS analysis demonstrates that, as was recently reported for the ORs in Drosophila melanogaster, mosquito ORs also have a topology different than their mammalian counterparts with their N-terminal ends located in the cytoplasm and their C-terminal ends facing outside the cell. These results set the stage for the production of mosquito ORs in quantities that should permit their detailed biochemical and structural characterization and the exploration of their functional properties