Homology Induction: the use of machine learning to improve sequence similarity searches

Karwath, A. King, R. Homology induction: the use of machine learning to improve sequence similarity searches. BMC Bioinformatics. 23rd April 2002. 3:11 Additional File Describes the title organims species declaration in one string [http://www.biomedcentral.com/content/supplementary/1471- 2105-3-11-S1.doc] Sponsorship: Andreas Karwath and Ross D. King were supported by the EPSRC grant GR/L62849.

Genes involved in complex adaptive processes tend to have highly conserved upstream regions in mammalian genomes

BACKGROUND:Recent advances in genome sequencing suggest a remarkable conservation in gene content of mammalian organisms. The similarity in gene repertoire present in different organisms has increased interest in studying regulatory mechanisms of gene expression aimed at elucidating the differences in phenotypes. In particular, a proximal promoter region contains a large number of regulatory elements that control the expression of its downstream gene. Although many studies have focused on identification of these elements, a broader picture on the complexity of transcriptional regulation of different biological processes has not been addressed in mammals. The regulatory complexity may strongly correlate with gene function, as different evolutionary forces must act on the regulatory systems under different biological conditions. We investigate this hypothesis by comparing the conservation of promoters upstream of genes classified in different functional categories.RESULTS:By conducting a rank correlation analysis between functional annotation and upstream sequence alignment scores obtained by human-mouse and human-dog comparison, we found a significantly greater conservation of the upstream sequence of genes involved in development, cell communication, neural functions and signaling processes than those involved in more basic processes shared with unicellular organisms such as metabolism and ribosomal function. This observation persists after controlling for G+C content. Considering conservation as a functional signature, we hypothesize a higher density of cis-regulatory elements upstream of genes participating in complex and adaptive processes.CONCLUSION:We identified a class of functions that are associated with either high or low promoter conservation in mammals. We detected a significant tendency that points to complex and adaptive processes were associated with higher promoter conservation, despite the fact that they have emerged relatively recently during evolution. We described and contrasted several hypotheses that provide a deeper insight into how transcriptional complexity might have been emerged during evolution.

Ferrochelatase is a conserved downstream target of the blue light-sensing White collar complex in fungi.

Light is a universal signal perceived by organisms, including fungi, in which light regulates common and unique biological processes depending on the species. Previous research has established that conserved proteins, originally called White collar 1 and 2 from the ascomycete Neurospora crassa, regulate UV/blue light sensing. Homologous proteins function in distant relatives of N. crassa, including the basidiomycetes and zygomycetes, which diverged as long as a billion years ago. Here we conducted microarray experiments on the basidiomycete fungus Cryptococcus neoformans to identify light-regulated genes. Surprisingly, only a single gene was induced by light above the commonly used twofold threshold. This gene, HEM15, is predicted to encode a ferrochelatase that catalyses the final step in haem biosynthesis from highly photoreactive porphyrins. The C. neoformans gene complements a Saccharomyces cerevisiae hem15Delta strain and is essential for viability, and the Hem15 protein localizes to mitochondria, three lines of evidence that the gene encodes ferrochelatase. Regulation of HEM15 by light suggests a mechanism by which bwc1/bwc2 mutants are photosensitive and exhibit reduced virulence. We show that ferrochelatase is also light-regulated in a white collar-dependent fashion in N. crassa and the zygomycete Phycomyces blakesleeanus, indicating that ferrochelatase is an ancient target of photoregulation in the fungal kingdom.

Antibody deficiency due to a missense mutation in CD19 demonstrates the importance of the conserved tryptophan 41 in immunoglobulin superfamily domain formation.

Immunoglobulin superfamily (IgSF) domains are conserved structures present in many proteins in eukaryotes and prokaryotes. These domains are well-capable of facilitating sequence variation, which is most clearly illustrated by the variable regions in immunoglobulins (Igs) and T cell receptors (TRs). We studied an antibody-deficient patient suffering from recurrent respiratory infections and with impaired antibody responses to vaccinations. Patient's B cells showed impaired Ca(2+) influx upon stimulation with anti-IgM and lacked detectable CD19 membrane expression. CD19 sequence analysis revealed a homozygous missense mutation resulting in a tryptophan to cystein (W52C) amino acid change. The affected tryptophan is CONSERVED-TRP 41 located on the C-strand of the first extracellular IgSF domain of CD19 and was found to be highly conserved, not only in mammalian CD19 proteins, but in nearly all characterized IgSF domains. Furthermore, the tryptophan is present in all variable domains in Ig and TR and was not mutated in 117 Ig class-switched transcripts of B cells from controls, despite an overall 10% amino acid change frequency. In vitro complementation studies and CD19 western blotting of patient's B cells demonstrated that the mutated protein remained immaturely glycosylated. This first missense mutation resulting in a CD19 deficiency demonstrates the crucial role of a highly conserved tryptophan in proper folding or stability of IgSF domains.

The DUB/USP17 deubiquitinating enzymes, a multigene family within a tandemly repeated sequence

Here we report the identification of 10 human, 1 murine, and 2 rat ORFs, all of which represent additional members of the DUB/USP17 family of deubiquitinating enzymes. In addition, we demonstrate that this family constitutes part of a tandemly repeated sequence conserved throughout humans, mice, and rats. Furthermore, upon examination of the known family members we have found that the multiple genes observed, in contrast to other gene families, have arisen due to the independent expansion of an ancestral sequence within each species. This premise is further strengthened by the observation that the murine and rat genes span two exons while their human counterparts have one. These observations, in conjunction with previous work demonstrating that the DUB/USP17's are cytokine inducible and that they regulate both cell growth and survival, suggest that the DUB/USP17's are a large highly conserved family of genes that may play an important role in controlling cell fate.

ADP-Ribose-1"-Monophosphatase: a Conserved Coronavirus Enzyme That Is Dispensable for Viral Replication in Tissue Culture

Replication of the ~30-kb plus-strand RNA genome of coronaviruses and synthesis of an extensive set of subgenome-length RNAs are mediated by the replicase-transcriptase, a membrane-bound protein complex containing several cellular proteins and up to 16 viral nonstructural proteins (nsps) with multiple enzymatic activities, including protease, polymerase, helicase, methyltransferase, and RNase activities. To get further insight into the replicase gene-encoded functions, we characterized the coronavirus X domain, which is part of nsp3 and has been predicted to be an ADP-ribose-1"-monophosphate (Appr-1"-p) processing enzyme. Bacterially expressed forms of human coronavirus 229E (HCoV-229E) and severe acute respiratory syndrome-coronavirus X domains were shown to dephosphorylate Appr-1"-p, a side product of cellular tRNA splicing, to ADP-ribose in a highly specific manner. The enzyme had no detectable activity on several other nucleoside phosphates. Guided by the crystal structure of AF1521, an X domain homolog from Archaeoglobus fulgidus, potential active-site residues of the HCoV-229E X domain were targeted by site-directed mutagenesis. The data suggest that the HCoV-229E replicase polyprotein residues, Asn 1302, Asn 1305, His 1310, Gly 1312, and Gly 1313, are part of the enzyme's active site. Characterization of an Appr-1"-pase-deficient HCoV-229E mutant revealed no significant effects on viral RNA synthesis and virus titer, and no reversion to the wild-type sequence was observed when the mutant virus was passaged in cell culture. The apparent dispensability of the conserved X domain activity in vitro indicates that coronavirus replicase polyproteins have evolved to include nonessential functions. The biological significance of the novel enzymatic activity in vivo remains to be investigated.

Conserved amino acid residues found in a predicted cytosolic domain of the lipopolysaccharide biosynthetic protein WecA are implicated in the recognition of UDP-N-acetylglucosamine

WecA, an integral membrane protein that belongs to a family of polyisoprenyl phosphate N-acetylhexosamine-1-phosphate transferases, is required for the biosynthesis of O-specific LPS and enterobacterial common antigen in Escherichia coli and other enteric bacteria. WecA functions as an UDP-N-acetylglucosamine (GlcNAc):undecaprenyl-phosphate GlcNAc-1-phosphate transferase. A conserved short sequence motif (His-Ile-His-His; HIHH) and a conserved arginine were identified in WecA at positions 279-282 and 265, respectively. This region is located within a predicted cytosolic segment common to all bacterial homologues of WecA. Both HIHH279-282 and the Arg265 are reminiscent of the HIGH motif (His-Ile-Gly-His) and a nearby upstream lysine, which contribute to the three-dimensional architecture of the nucleotide-binding site among various enzymes displaying nucleotidyltransferase activity. Thus, it was hypothesized that these residues may play a role in the interaction of WecA with UDP-GlcNAc. Replacement of the entire HIHH motif by site-directed mutagenesis produced a protein that, when expressed in the E. coli wecA mutant MV501, did not complement the synthesis of O7 LPS. Membrane extracts containing the mutated protein failed to transfer UDP-GlcNAc into a lipid-rich fraction and to bind the UDP-GlcNAc analogue tunicamycin. Similar results were obtained by individually replacing the first histidine (H279) of the HIHH motif as well as the Arg265 residue. The functional importance of these residues is underscored by the high level of conservation of H279 and Arg265 among bacterial WecA homologues that utilize several different UDP-N-acetylhexosamine substrates.

A conserved spider silk domain acts as a molecular switch that controls fibre assembly

A huge variety of proteins are able to form fibrillar structures(1), especially at high protein concentrations. Hence, it is surprising that spider silk proteins can be stored in a soluble form at high concentrations and transformed into extremely stable fibres on demand(2,3). Silk proteins are reminiscent of amphiphilic block copolymers containing stretches of polyalanine and glycine-rich polar elements forming a repetitive core flanked by highly conserved non-repetitive amino-terminal(4,5) and carboxy-terminal(6) domains. The N-terminal domain comprises a secretion signal, but further functions remain unassigned. The C-terminal domain was implicated in the control of solubility and fibre formation(7) initiated by changes in ionic composition(8,9) and mechanical stimuli known to align the repetitive sequence elements and promote beta-sheet formation(10-14). However, despite recent structural data(15), little is known about this remarkable behaviour in molecular detail. Here we present the solution structure of the C-terminal domain of a spider dragline silk protein and provide evidence that the structural state of this domain is essential for controlled switching between the storage and assembly forms of silk proteins. In addition, the C-terminal domain also has a role in the alignment of secondary structural features formed by the repetitive elements in the backbone of spider silk proteins, which is known to be important for the mechanical properties of the fibre.

Basal and antigen-induced exposure of the proline-rich sequence in CD3ε.

The CD3ε cytoplasmic tail contains a conserved proline-rich sequence (PRS) that influences TCR-CD3 expression and signaling. Although the PRS can bind the SH3.1 domain of the cytosolic adapter Nck, whether the PRS is constitutively available for Nck binding or instead represents a cryptic motif that is exposed via conformational change upon TCR-CD3 engagement (CD3Δc) is currently unresolved. Furthermore, the extent to which a cis-acting CD3ε basic amino acid-rich stretch (BRS), with its unique phosphoinositide-binding capability, might impact PRS accessibility is not clear. In this study, we found that freshly harvested primary thymocytes expressed low to moderate basal levels of Nck-accessible PRS ("open-CD3"), although most TCR-CD3 complexes were inaccessible to Nck ("closed-CD3"). Ag presentation in vivo induced open-CD3, accounting for half of the basal level found in thymocytes from MHC(+) mice. Additional stimulation with either anti-CD3 Abs or peptide-MHC ligands further elevated open-CD3 above basal levels, consistent with a model wherein antigenic engagement induces maximum PRS exposure. We also found that the open-CD3 conformation induced by APCs outlasted the time of ligand occupancy, marking receptors that had been engaged. Finally, CD3ε BRS-phosphoinositide interactions played no role in either adoption of the initial closed-CD3 conformation or induction of open-CD3 by Ab stimulation. Thus, a basal level of open-CD3 is succeeded by a higher, induced level upon TCR-CD3 engagement, involving CD3Δc and prolonged accessibility of the CD3ε PRS to Nck.

New insights into small molecules inhibitors and protein-protein interactions of VirB8 : a critical conserved component of the type IV secretion system.

Les systèmes bactériens de sécrétion de type IV (T4SS) sont constitués d’un ensemble de 8 à 12 protéines conservées. Ces dernières sont utilisées lors de la translocation de protéines, la translocation de complexes ADN-protéines mais aussi pour le transport de ces derniers au travers de la membrane cellulaire. Les T4SS, en tant que facteurs de virulence pour beaucoup de pathogènes comme Brucella suis, sont donc d’excellents modèles cibles pour le développement de médicaments d’antivirulence. Ces médicaments, en privant le pathogène de son facteur essentiel de virulence : le T4SS, constituent une alternative ou encore une amélioration des traitements antibiotiques utilisés actuellement. VirB8, un facteur d’assemblage conservé dans le T4SS, forme des dimères qui sont importants pour la fonction des T4SS dans ces pathogènes. De par ses interactions multiples, VirB8 est un excellent modèle pour l’analyse des facteurs d’assemblage mais aussi en tant que cible de médicaments qui empêcheraient son interaction avec d’autres protéines et qui, in fine, désarmeraient les bactéries en les privant de leur fonctions essentielles de virulence. À ce jour, nous savons qu’il existe un équilibre monomère-dimère et un processus d’homodimerization de VirB8 dont l’importance est vitale pour la fonctionnement biologique des T4SSs. En se basant sur des essais quantitatifs d’interaction, nous avons identifié (i) des sites potentiels d’interaction avec d’autres protéines VirB du T4SS mais aussi (ii) isolé des petites molécules inhibitrices afin de tester la fonction protéique de VirB8. Afin de déterminer les acides aminés importants pour l’hétérodimérization de VirB8 avec VirB10, nous avons effectué des expériences de mutagenèse aléatoire, de phage display et d’arrimage moléculaire in silico. Ces expériences ont démontré l’importance de trois acides aminés localisés sur le feuillet β : R160, S162, T164 et I165. Ces derniers seraient importants pour l’association de VirB8 avec VirB10 étant donné que leur mutagenèse entraine une diminution de la formation du complexe VirB8-VirB10. L’objectif actuel de notre projet de recherche est de pouvoir mieux comprendre mais aussi d’évaluer le rôle de VirB8 dans l’assemblage du T4SS. Grace à un méthode de criblage adaptée à partir de la structure de VirB8, nous avons pu identifié une petite molécule inhibitrice BAR-068, qui aurait un rôle prometteur dans l’inhibition du T4SS. Nous avons utilisé la spectroscopie par fluorescence, l’essai à deux hybrides, le cross-linking et la cristallographie afin de déterminer le mécanisme d'interaction existant entre VirB8 et BAR-068. Ces travaux pourraient permettre de nombreuses avancées, notamment en termes de compréhension des mécanismes d’inhibition du T4SS. Notre objectif ultime est de pouvoir caractériser la séquence d’évènements essentiels à l’assemblage et au fonctionnement du T4SS. De manière globale, notre projet de recherche permettrait de révéler les grands principes d’assemblage des protéines membranaires, les processus de sécrétion de protéines chez les bactéries mais aussi de proposer une nouvelle stratégie lors du développement de drogues antimicrobiennes.

Phylogeny, function, and evolution of the cupins, a structurally conserved, functionally diverse superfamily of proteins

The cupin superfamily is a group of functionally diverse proteins that are found in all three kingdoms of life, Archaea, Eubacteria, and Eukaryota. These proteins have a characteristic signature domain comprising two histidine- containing motifs separated by an intermotif region of variable length. This domain consists of six beta strands within a conserved beta barrel structure. Most cupins, such as microbial phosphomannose isomerases (PMIs), AraC- type transcriptional regulators, and cereal oxalate oxidases (OXOs), contain only a single domain, whereas others, such as seed storage proteins and oxalate decarboxylases (OXDCs), are bi-cupins with two pairs of motifs. Although some cupins have known functions and have been characterized at the biochemical level, the majority are known only from gene cloning or sequencing projects. In this study, phylogenetic analyses were conducted on the conserved domain to investigate the evolution and structure/function relationships of cupins, with an emphasis on single- domain plant germin-like proteins (GLPs). An unrooted phylogeny of cupins from a wide spectrum of evolutionary lineages identified three main clusters, microbial PMIs, OXDCs, and plant GLPs. The sister group to the plant GLPs in the global analysis was then used to root a phylogeny of all available plant GLPs. The resulting phylogeny contained three main clades, classifying the GLPs into distinct subfamilies. It is suggested that these subfamilies correlate with functional categories, one of which contains the bifunctional barley germin that has both OXO and superoxide dismutase (SOD) activity. It is proposed that GLPs function primarily as SODs, enzymes that protect plants from the effects of oxidative stress. Closer inspection of the DNA sequence encoding the intermotif region in plant GLPs showed global conservation of thymine in the second codon position, a character associated with hydrophobic residues. Since many of these proteins are multimeric and enzymatically inactive in their monomeric state, this conservation of hydrophobicity is thought to be associated with the need to maintain the various monomer- monomer interactions. The type of structure-based predictive analysis presented in this paper is an important approach for understanding gene function and evolution in an era when genomes from a wide range of organisms are being sequenced at a rapid rate.

Leafy, terminal flower1 and agamous are functionally conserved but do not regulate terminal flowering and floral determinacy in Impatiens balsamina

In Impatiens balsamina a lack of commitment of the meristem during floral development leads to the continuous requirement for a leaf-derived floral signal. In the absence of this signal the meristem reverts to leaf production. Current models for Arabidopsis state that LEAFY (LFY) is central to the integration of floral signals and regulates flowering partly via interactions with TERMINAL FLOWER1 (TFL1) and AGAMOUS (AG). Here we describe Impatiens homologues of LFY, TFL1 and AG (IbLFY, IbTFL1 and IbAG) that are highly conserved at a sequence level and demonstrate homologous functions when expressed ectopically in transgenic Arabidopsis. We relate the expression patterns of IbTFL1 and IbAG to the control of terminal flowering and floral determinacy in Impatiens. IbTFL1 is involved in controlling the phase of the axillary meristems and is expressed in axillary shoots and axillary meristems which produce inflorescences, but not in axillary flowers. It is not involved in maintaining the terminal meristem in either an inflorescence or indeterminate state. Terminal flowering in Impatiens appears therefore to be controlled by a pathway that uses a different integration system than that regulating the development of axillary flowers and branches. The pattern of ovule production in Impatiens requires the meristem to be maintained after the production of carpels. Consistent with this morphological feature IbAG appears to specify stamen and carpel identity, but is not sufficient to specify meristem determinacy in Impatiens.

Relation of pretreatment sequence diversity in NS5A region of HCV genotype 1 with immune response between pegylated-INF/ribavirin therapy outcomes

Hepatitis C virus (HCV) infection frequently persists despite substantial virus-specific immune responses and the combination of pegylated interferon (INF)-alpha and ribavirin therapy. Major histocompatibility complex class I restricted CD8+ T cells are responsible for the control of viraemia in HCV infection, and several studies suggest protection against viral infection associated with specific HLAs. The reason for low rates of sustained viral response (SVR) in HCV patients remains unknown. Escape mutations in response to cytotoxic T lymphocyte are widely described; however, its influence in the treatment outcome is ill understood. Here, we investigate the differences in CD8 epitopes frequencies from the Los Alamos database between groups of patients that showed distinct response to pegylated alpha-INF with ribavirin therapy and test evidence of natural selection on the virus in those who failed treatment, using five maximum likelihood evolutionary models from PAML package. The group of sustained virological responders showed three epitopes with frequencies higher than Non-responders group, all had statistical support, and we observed evidence of selection pressure in the last group. No escape mutation was observed. Interestingly, the epitope VLSDFKTWL was 100% conserved in SVR group. These results suggest that the response to treatment can be explained by the increase in immune pressure, induced by interferon therapy, and the presence of those epitopes may represent an important factor in determining the outcome of therapy.

Functionally conserved cis-regulatory elements of COL18A1 identified through zebrafish transgenesis

Type XVIII collagen is a component of basement membranes, and expressed prominently in the eye, blood vessels, liver, and the central nervous system. Homozygous mutations in COL18A1 lead to Knobloch Syndrome, characterized by ocular defects and occipital encephalocele. However, relatively little has been described on the role of type XVIII collagen in development, and nothing is known about the regulation of its tissue-specific expression pattern. We have used zebrafish transgenesis to identify and characterize cis-regulatory sequences controlling expression of the human gene. Candidate enhancers were selected from non-coding sequence associated with COL18A1 based on sequence conservation among mammals. Although these displayed no overt conservation with orthologous zebrafish sequences, four regions nonetheless acted as tissue-specific transcriptional enhancers in the zebrafish embryo, and together recapitulated the major aspects of col18a1 expression. Additional post-hoc computational analysis on positive enhancer sequences revealed alignments between mammalian and teleost sequences, which we hypothesize predict the corresponding zebrafish enhancers; for one of these, we demonstrate functional overlap with the orthologous human enhancer sequence. Our results provide important insight into the biological function and regulation of COL18A1, and point to additional sequences that may contribute to complex diseases involving COL18A1. More generally, we show that combining functional data with targeted analyses for phylogenetic conservation can reveal conserved cis-regulatory elements in the large number of cases where computational alignment alone falls short. (C) 2009 Elsevier Inc. All rights reserved.

Sequence and function of lysosomal and digestive cathepsin D-like proteinases of Musca domestica midgut

Musca domestica larvae display in anterior and middle midgut contents, a proteolytic activity with pH optimum of 3.0-3.5 and kinetic properties like cathepsin D. Three cDNAs coding for preprocathepsin D-like proteinases (ppCAD 1, ppCAD 2, ppCAD 3) were cloned from a M. domestica midgut cDNA library. The coded protein sequences included the signal peptide, propeptide and mature enzyme that has all conserved catalytic and substrate binding residues found in bovine lysosomal cathepsin D. Nevertheless, ppCAD 2 and ppCAD 3 lack the characteristic proline loop and glycosylation sites. A comparison among the sequences of cathepsin D-like enzymes from some vertebrates and those found in M. domestica and in the genomes of Aedes aegypti, Drosophila melanogaster, Tribolium castaneum, and Bombyx mori showed that only flies have enzymes lacking the proline loop (as defined by the motif: DxPxPx(G/A)P), thus resembling vertebrate pepsin. ppCAD 3 should correspond to the digestive cathepsin D-like proteinase (CAD) found in enzyme assays because: (1) it seems to be the most expressed CAD, based on the frequency of ESTs found. (2) The mRNA for CAD 3 is expressed only in the anterior and proximal middle midgut. (3) Recombinant procathepsin D-like proteinase (pCAD 3), after auto-activation has a pH optimum of 2.5-3.0 that is close to the luminal pH of M. domestica midgut. (4) Immunoblots of proteins from different tissues revealed with anti-pCAD 3 serum were positive only in samples of anterior and middle midgut tissue and contents. (5) CAD 3 is localized with immunogold inside secretory vesicles and around microvilli in anterior and middle midguit cells. The data support the view that on adapting to deal with a bacteria-rich food in an acid midgut region, M. domestica digestive CAD resulted from the same archetypical gene as the intracellular cathepsin D, paralleling what happened with vertebrates. The lack of the proline loop may be somehow associated with the extracellular role of both pepsin and digestive CAD 3. (C) 2009 Elsevier Ltd. All rights reserved.