Resolving the fibrotic niche of human liver cirrhosis at single-cell level

Liver cirrhosis is a major cause of death worldwide and is characterized by extensive fibrosis. There are currently no effective antifibrotic therapies available. To obtain a better understanding of the cellular and molecular mechanisms involved in disease pathogenesis and enable the discovery of therapeutic targets, here we profile the transcriptomes of more than 100,000 single human cells, yielding molecular definitions for non-parenchymal cell types that are found in healthy and cirrhotic human liver. We identify a scar-associated TREM2+CD9+ subpopulation of macrophages, which expands in liver fibrosis, differentiates from circulating monocytes and is pro-fibrogenic. We also define ACKR1+ and PLVAP+ endothelial cells that expand in cirrhosis, are topographically restricted to the fibrotic niche and enhance the transmigration of leucocytes. Multi-lineage modelling of ligand and receptor interactions between the scar-associated macrophages, endothelial cells and PDGFRα+ collagen-producing mesenchymal cells reveals intra-scar activity of several pro-fibrogenic pathways including TNFRSF12A, PDGFR and NOTCH signalling. Our work dissects unanticipated aspects of the cellular and molecular basis of human organ fibrosis at a single-cell level, and provides a conceptual framework for the discovery of rational therapeutic targets in liver cirrhosis.

ELIXIR-UK role in bioinformatics training at the national level and across ELIXIR

ELIXIR-UK is the UK node of ELIXIR, the European infrastructure for life science data. Since its foundation in 2014, ELIXIR-UK has played a leading role in training both within the UK and in the ELIXIR Training Platform, which coordinates and delivers training across all ELIXIR members. ELIXIR-UK contributes to the Training Platform’s coordination and supports the development of training to address key skill gaps amongst UK scientists. As part of this work it acts as a conduit for nationally-important bioinformatics training resources to promote their activities to the ELIXIR community. ELIXIR-UK also leads ELIXIR’s flagship Training Portal, TeSS, which collects information about a diverse range of training and makes it easily accessible to the community. ELIXIR-UK also works with others to provide key digital skills training, partnering with the Software Sustainability Institute to provide Software Carpentry training to the ELIXIR community and to establish the Data Carpentry initiative, and taking a lead role amongst national stakeholders to deliver the StaTS project – a coordinated effort to drive engagement with training in statistics.

Phospholipase Cζ, the trigger of egg activation in mammals, is present in a non-mammalian species

The activation of the egg to begin development into an embryo is triggered by a sperm-induced increase in intracellular egg Ca2+. There has been much controversy about how the sperm induces this fundamental developmental event, but recent studies suggest that, in mammals, egg activation is triggered by a testis-specific phospholipase C: PLCζ. Since the discovery of PLCζ, it has been unclear whether its role in triggering egg activation is common to all vertebrates, or is confined to mammals. Here, we demonstrate for the first time that PLCζ is present in a non-mammalian vertebrate. Using genomic and cDNA databases, we have identified the cDNA encoding a PLCζ orthologue in the domestic chicken that, like the mammalian isoforms, is a testis-specific gene. The chicken PLCζ cDNA is 2152 bp in size and encodes an open reading frame of 639 amino acids. When injected into mouse oocytes, chicken PLCζ cRNA triggers Ca2+oscillations, indicating that it has functional properties similar to those of mammalian PLCζ. Our findings suggest that PLCζ may have a universal role in triggering egg activation in vertebrates.

Stromal cell-derived receptor 2 and cytochrome b561 are functional ferric reductases

Iron has a variety of functions in cellular organisms ranging from electron transport and DNA synthesis to adenosine triphosphate (ATP) and neurotransmitter synthesis. Failure to regulate the homeostasis of iron can lead to cognition and demyelination disorders when iron levels are deficient, and to neurodegenerative disorders when iron is in excess. In this study we show that three members of the b561 family of predicted ferric reductases, namely mouse cytochrome b561 and mouse and fly stromal cell-derived receptor 2 (SDR2), have ferric reductase activity. Given that a fourth member, duodenal cytochrome b (Dcytb), has previously been shown to be a ferric reductase, it is likely that all remaining members of this family also exhibit this activity. Furthermore, we show that the rat sdr2 message is predominantly expressed in the liver and kidney, with low expression in the duodenum. In hypotransferrinaemic (hpx) mice, sdr2 expression in the liver and kidney is reduced, suggesting that it may be regulated by iron. Moreover, we demonstrate the presence of mouse sdr2 in the choroid plexus and in the ependymal cells lining the four ventricles, through in situ hybridization analysis.

Cloning and developmental expression analysis of ltd-1, the Caenorhabditis elegans homologue of the mouse kyphoscoliosis (ky) gene

We have characterized the developmental expression pattern of the Caenorhabditis elegans homologue of the mouse ky gene. The Ky protein has a putative key function in muscle development and has homologues in invertebrates, fungi and a cyanobacterium. The C. elegans Ky homologue gene has been named ltd-1 for LIM and transglutaminase domains gene. The LTD-1::GFP construct is expressed in developing hypodermal cells from the twofold stage embryo through adulthood. These data define the ltd-1 gene as a novel marker for C. elegans epithelial cell development.

Genome cartography through domain annotation

The evolutionary history of eukaryotic proteins involves rapid sequence divergence, addition and deletion of domains, and fusion and fission of genes. Although the protein repertoires of distantly related species differ greatly, their domain repertoires do not. To account for the great diversity of domain contexts and an unexpected paucity of ortholog conservation, we must categorize the coding regions of completely sequenced genomes into domain families, as well as protein families.

Novel protein domains and repeats in Drosophila melanogaster: insights into structure, function, and evolution

Sequence database searching methods such as BLAST, are invaluable for predicting molecular function on the basis of sequence similarities among single regions of proteins. Searches of whole databases however, are not optimized to detect multiple homologous regions within a single polypeptide. Here we have used the prospero algorithm to perform self-comparisons of all predicted Drosophila melanogaster gene products. Predicted repeats, and their homologs from all species, were analyzed further to detect hitherto unappreciated evolutionary relationships. Results included the identification of novel tandem repeats in the human X-linked retinitis pigmentosa type-2 gene product, repeated segments in cystinosin, associated with a defect in cystine transport, and 'nested' homologous domains in dysferlin, whose gene is mutated in limb girdle muscular dystrophy. Novel signaling domain families were found that may regulate the microtubule-based cytoskeleton and ubiquitin-mediated proteolysis, respectively. Two families of glycosyl hydrolases were shown to contain internal repetitions that hint at their evolution via a piecemeal, modular approach. In addition, three examples of fruit fly genes were detected with tandem exons that appear to have arisen via internal duplication. These findings demonstrate how completely sequenced genomes can be exploited to further understand the relationships between molecular structure, function, and evolution.

A new sequence motif linking lissencephaly, Treacher Collins and oral-facial-digital type 1 syndromes, microtubule dynamics and cell migration

A previously unidentified sequence motif has been identified in the products of genes mutated in Miller-Dieker lissencephaly, Treacher Collins, oral-facial-digital type 1 and contiguous syndrome ocular albinism with late onset sensorineural deafness syndromes. An additional homologous motif was detected in a gene product fused to the fibroblast growth factor receptor type 1 in patients with an atypical stem cell myeloproliferative disorder. In total, over 100 eukaryotic intracellular proteins are shown to possess a LIS1 homology (LisH) motif, including several katanin p60 subunits, muskelin, tonneau, LEUNIG, Nopp140, aimless and numerous WD repeat-containing beta-propeller proteins. It is suggested that LisH motifs contribute to the regulation of microtubule dynamics, either by mediating dimerization, or else by binding cytoplasmic dynein heavy chain or microtubules directly. The predicted secondary structure of LisH motifs, and their occurrence in homologues of Gbeta beta-propeller subunits, suggests that they are analogues of Ggamma subunits, and might associate with the periphery of beta-propeller domains. The finding of LisH motifs in both treacle and Nopp140 reinforces previous observations of functional similarities between these nucleolar proteins. Uncharacterized LisH motif-containing proteins represent candidates for other diseases associated with aberrant microtubule dynamics and defects of cell migration, nucleokinesis or chromosome segregation.

Sequence variation and disease in the wake of the draft human genome

The sequencing phase of the human genome project will soon be over. In its wake, repertoires of sequence polymorphisms among the human population are being sampled and a battery of functional genomics projects, from gene and protein expression studies to whole proteome interaction experiments, are generating vast quantities of data. Now that the data, or the means to generate data, are available it is the application of this information in enhancing our understanding of biology that represents the next formidable challenge. Two prominent issues should be considered. First, existing data must be analysed using the best methods available. The prediction of enzymatic activity for bestrophin, whose gene is mutated in Best macular dystrophy, is described in this review. This is an example of the experimentally testable hypotheses that can result from such detailed and exhaustive analyses. Secondly, the torrents of data from high-throughput studies will need to be made more accessible to all using web-based resources that integrate and digest complementary data types. The internet sites that showcase the human genome sequence are blazing a new trail. Ultimately, the success of genome sequencing and functional genomics will be measured not by the quantity and accuracy of raw data generated, but how rapidly they can be harnessed to span the divide between genotype and phenotype.

CHROMA: consensus-based colouring of multiple alignments for publication

CHROMA annotates multiple protein sequence alignments by consensus to produce formatted and coloured text suitable for incorporation into other documents for publication. The package is designed to be flexible and reliable, and has a simple-to-use graphical user interface running under Microsoft Windows. Both the executables and source code for CHROMA running under Windows and Linux (portable command-line only) are freely available at http://www.lg.ndirect.co.uk/chroma. Software enquiries should be directed to CHROMA@lg.ndirect.co.uk.

Domain homologues of dopamine beta-hydroxylase and ferric reductase: roles for iron metabolism in neurodegenerative disorders?

One of the defining characteristics of neurodegenerative diseases, including Parkinson's, Alzheimer's and Huntington's diseases, is abnormal accumulations of iron, specifically in affected areas. Following injection of iron in rat brains, a relatively selective lesion of dopamine neurons, similar to parkinsonism, occurs. These observations indicate that Fe(II)-mediated generation of free radical species, by the Fenton reaction, might contribute to the pathoetiology of these diseases. Iron is known to possess multiple roles in the biosynthesis of catecholamines in dopaminergic neurons. These include, as Fe(II), facilitating the production of dopamine from phenylalanine by tyrosine hydroxylase, and as heme, assisting the recycling of ascorbate by cytochrome b-561 required for the generation of norepinephrine from dopamine by dopamine beta-hydroxylase. In this study, it is demonstrated that a human and mouse gene product, stromal cell-derived receptor 2, is a homologue of cytochrome b-561 and duodenal cytochrome b, and is thus predicted to be active as a ferric reductase. Moreover, this protein also contains a domain homologous to the N-terminal regulatory region of dopamine beta-hydroxylase. These findings from sequence analysis lead to a prediction that stromal cell-derived receptor 2 is a catecholamine-regulated ferric reductase active in the brain. Dysfunction of cytochrome b-561 or stromal cell-derived receptor 2, therefore, might predispose individuals to abnormal accumulation of Fe(III) and/or generation of cytotoxic free radicals as a consequence of a rapid cycling between Fe(III) and Fe(II). The hypothesis that aberrant ferric reductase activities are involved in the progression of neurodegenerative diseases should open up new avenues of research, and possibly therapy, for these devastating diseases.

On the evolution of protein folds: are similar motifs in different protein folds the result of convergence, insertion, or relics of an ancient peptide world?

This paper presents and discusses evidence suggesting how the diversity of domain folds in existence today might have evolved from peptide ancestors. We apply a structure similarity detection method to detect instances where localized regions of different protein folds contain highly similar sequences and structures. Results of performing an all-on-all comparison of known structures are described and compared with other recently published findings. The numerous instances of local sequence and structure similarities within different protein folds, together with evidence from proteins containing sequence and structure repeats, argues in favor of the evolution of modern single polypeptide domains from ancient short peptide ancestors (antecedent domain segments (ADSs)). In this model, ancient protein structures were formed by self-assembling aggregates of short polypeptides. Subsequently, and perhaps concomitantly with the evolution of higher fidelity DNA replication and repair systems, single polypeptide domains arose from the fusion of ADSs genes. Thus modern protein domains may have a polyphyletic origin.

Protein repeats: structures, functions, and evolution

Internal repetition within proteins has been a successful strategem on multiple separate occasions throughout evolution. Such protein repeats possess regular secondary structures and form multirepeat assemblies in three dimensions of diverse sizes and functions. In general, however, internal repetition affords a protein enhanced evolutionary prospects due to an enlargement of its available binding surface area. Constraints on sequence conservation appear to be relatively lax, due to binding functions ensuing from multiple, rather than, single repeats. Considerable sequence divergence as well as the short lengths of sequence repeats mean that repeat detection can be a particularly arduous task. We also consider the conundrum of how multiple repeats, which show strong structural and functional interdependencies, ever evolved from a single repeat ancestor. In this review, we illustrate each of these points by referring to six prolific repeat types (repeats in beta-propellers and beta-trefoils and tetratricopeptide, ankyrin, armadillo/HEAT, and leucine-rich repeats) and in other less-prolific but nonetheless interesting repeats.

Plagiarized bacterial genes in the human book of life

The initial analysis of the human genome draft sequence reveals that our 'book of life' is multi-authored. A small but significant proportion of our genes owes their heritage not to antecedent eukaryotes but instead to bacteria. The publicly funded Human Genome Project study indicates that about 0.5% of all human genes were copied into the genome from bacterial sources. Detailed sequence analyses point to these 'horizontal gene transfer' events having occurred relatively recently. So how did the human 'book of life' evolve to be a chimaera, part animal and part bacterium? And what was the probable evolutionary impact of such gene plagiarism?

Syncoilin, a novel member of the intermediate filament superfamily that interacts with alpha-dystrobrevin in skeletal muscle

Dystrophin coordinates the assembly of a complex of structural and signaling proteins that are required for normal muscle function. A key component of the dystrophin protein complex is alpha-dystrobrevin, a dystrophin-associated protein whose absence results in neuromuscular junction defects and muscular dystrophy. To gain further insights into the role of alpha-dystrobrevin in skeletal muscle, we used the yeast two-hybrid system to identify a novel alpha-dystrobrevin-binding partner called syncoilin. Syncoilin is a new member of the intermediate filament superfamily and is highly expressed in skeletal and cardiac muscle. In normal skeletal muscle, syncoilin is concentrated at the neuromuscular junction, where it colocalizes and coimmunoprecipitates with alpha-dystrobrevin-1. Expression studies in mammalian cells demonstrate that, while alpha-dystrobrevin and syncoilin associate directly, overexpression of syncoilin does not result in the self-assembly of intermediate filaments. Finally, unlike many components of the dystrophin protein complex, we show that syncoilin expression is up-regulated in dystrophin-deficient muscle. These data suggest that alpha-dystrobrevin provides a link between the dystrophin protein complex and the intermediate filament network at the neuromuscular junction, which may be important for the maintenance and maturation of the synapse.

Issues in predicting protein function from sequence

Identifying homologues, defined as genes that arose from a common evolutionary ancestor, is often a relatively straightforward task, thanks to recent advances made in estimating the statistical significance of sequence similarities found from database searches. The extent by which homologues possess similarities in function, however, is less amenable to statistical analysis. Consequently, predicting function by homology is a qualitative, rather than quantitative, process and requires particular care to be taken. This review focuses on the various approaches that have been developed to predict function from the scale of the atom to that of the organism. Similarities in homologues' functions differ considerably at each of these different scales and also vary for different domain families. It is argued that due attention should be paid to all available clues to function, including orthologue identification, conservation of particular residue types, and the co-occurrence of domains in proteins. Pitfalls in database searching methods arising from amino acid compositional bias and database size effects are also discussed.

Initial sequencing and analysis of the human genome

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

Sialidase-like Asp-boxes: sequence-similar structures within different protein folds

Sequence similarity is the most common measure currently used to infer homology between proteins. Typically, homologous protein domains show sequence similarity over their entire lengths. Here we identify Asp box motifs, initially found as repeats in sialidases and neuraminidases, in new structural and sequence contexts. These motifs represent significantly similar sequences, localized to beta hairpins within proteins that are otherwise different in sequence and three-dimensional structure. By performing a combined sequence- and structure-based analysis we detect Asp boxes in more than nine protein families, including bacterial ribonucleases, sulfite oxidases, reelin, netrins, some lipoprotein receptors, and a variety of glycosyl hydrolases. Although the function common to each of these proteins, if any, remains unclear, we discuss possible functions of Asp boxes on the basis of previously determined experimental results and discuss different evolutionary scenarios for the origin of Asp-box containing proteins.

JmjC: cupin metalloenzyme-like domains in jumonji, hairless and phospholipase A2beta

On the basis of significant sequence similarity, we have identified JmjC domains in more than 100 eukaryotic and bacterial sequences. These include human hairless, mutated in individuals with alopecia universalis, retinoblastoma-binding protein 2 and several putative chromatin-associated proteins. JmjC domains are predicted to be metalloenzymes that adopt the cupin fold, and are candidates for enzymes that regulate chromatin remodelling.

The kyphoscoliosis (ky) mouse is deficient in hypertrophic responses and is caused by a mutation in a novel muscle-specific protein

The ky mouse mutant exhibits a primary degenerative myopathy preceding chronic thoraco-lumbar kyphoscoliosis. The histopathology of the ky mutant suggests that Ky protein activity is crucial for normal muscle growth and function as well as the maturation and stabilization of the neuromuscular junction. Muscle hypertrophy in response to increasing demand is deficient in the ky mutant, whereas adaptive fibre type shifts take place. The ky locus has previously been localized to a small region of mouse chromosome 9 and we have now identified the gene and the mutation underlying the kyphoscoliotic mouse. The ky transcript encodes a novel protein that is detected only in skeletal muscle and heart. The identification of the ky gene will allow detailed analysis of the impact of primary myopathy on idiopathic scoliosis in mice and man.

Alternative splicing of dystrobrevin regulates the stoichiometry of syntrophin binding to the dystrophin protein complex

Dystrophin coordinates the assembly of a complex of structural and signalling proteins that is required for normal muscle function. A key component of the dystrophin-associated protein complex (DPC) is alpha-dystrobrevin, a dystrophin-related and -associated protein whose absence results in muscular dystrophy and neuromuscular junction defects [1,2]. The current model of the DPC predicts that dystrophin and dystrobrevin each bind a single syntrophin molecule [3]. The syntrophins are PDZ-domain-containing proteins that facilitate the recruitment of signalling proteins such as nNOS (neuronal nitric oxide synthase) to the DPC [4]. Here we show, using yeast two-hybrid analysis and biochemical binding studies, that alpha-dystrobrevin in fact contains two independent syntrophin-binding sites in tandem. The previously undescribed binding site is situated within an alternatively spliced exon of alpha-dystrobrevin, termed the variable region-3 (vr3) sequence, which is specifically expressed in skeletal and cardiac muscle [5,6]. Analysis of the syntrophin-binding region of dystrobrevin reveals a tandem pair of predicted alpha helices with significant sequence similarity. These alpha helices, each termed a syntrophin-binding motif, are also highly conserved in dystrophin and utrophin. Together these data show that there are four potential syntrophin-binding sites per dystrophin complex in skeletal muscle: two on dystrobrevin and two on dystrophin or utrophin. Furthermore, alternative splicing of dystrobrevin provides a mechanism for regulating the stoichiometry of syntrophin association with the DPC. This is likely to have important consequences for the recruitment of specific signalling molecules to the DPC and ultimately for its function.

Proteins of the endoplasmic-reticulum-associated degradation pathway: domain detection and function prediction

Sequence database searches, using iterative-profile and Hidden-Markov-model approaches, were used to detect hitherto-undetected homologues of proteins that regulate the endoplasmic reticulum (ER)-associated degradation pathway. The translocon-associated subunit Sec63p (Sec=secretory) was shown to contain a domain of unknown function found twice in several Brr2p-like RNA helicases (Brr2=bad response to refrigeration 2). Additionally, Cue1p (Cue=coupling of ubiquitin conjugation to ER degradation), a yeast protein that recruits the ubiquitin-conjugating (UBC) enzyme Ubc7p to an ER-associated complex, was found to be one of a large family of putative scaffolding-domain-containing proteins that include the autocrine motility factor receptor and fungal Vps9p (Vps=vacuolar protein sorting). Two other yeast translocon-associated molecules, Sec72p and Hrd3p (Hrd=3-hydroxy-3-methylglutaryl-CoA reductase degradation), were shown to contain multiple tetratricopeptide-repeat-like sequences. From this observation it is suggested that Sec72p associates with a heat-shock protein, Hsp70, in a manner analogous to that known for Hop (Hsp70/Hsp90 organizing protein). Finally, the luminal portion of Ire1p (Ire=high inositol-requiring), thought to convey the sensing function of this transmembrane kinase and endoribonuclease, was shown to contain repeats similar to those in beta-propeller proteins. This finding hints at the mechanism by which Ire1p may sense extended unfolded proteins at the expense of compact folded molecules.

Identification of distant homologues of fibroblast growth factors suggests a common ancestor for all beta-trefoil proteins

Determination of the structures of fibroblast growth factors and interleukin-1s has previously revealed that they both adopt a beta-trefoil fold, similar to those found in Kunitz soybean trypsin inhibitors, ricin-like toxins, plant agglutinins and hisactophilin. These families possess distinct functions and occur in different subcellular localisations, and they appear to lack significant similarities in their sequences, ligands and modes of ligand binding. We have analysed the significance of sequence identities observed after structure alignment and provide statistical evidence that these beta-trefoil proteins are all homologues, having arisen from a common ancestor. In addition, we have explored the sequence space of all beta-trefoil proteins and have determined that the actin-binding proteins fascins, and other proteins of unknown function, are beta-trefoil family homologues. Unlike other beta-trefoil proteins, the triplicated repeats in each of the four beta-trefoil domains of fascins are significantly similar in sequence. This hints at how the beta-trefoil fold arose from the duplication of an ancestral gene encoding a homotrimeric single-repeat protein. The combined analysis of structure and sequence databases for detecting significant similarities is suggested as a highly sensitive approach to determining the common ancestry of extremely divergent homologues.

Regulation of chromatin structure by site-specific histone H3 methyltransferases

The organization of chromatin into higher-order structures influences chromosome function and epigenetic gene regulation. Higher-order chromatin has been proposed to be nucleated by the covalent modification of histone tails and the subsequent establishment of chromosomal subdomains by non-histone modifier factors. Here we show that human SUV39H1 and murine Suv39h1--mammalian homologues of Drosophila Su(var)3-9 and of Schizosaccharomyces pombe clr4--encode histone H3-specific methyltransferases that selectively methylate lysine 9 of the amino terminus of histone H3 in vitro. We mapped the catalytic motif to the evolutionarily conserved SET domain, which requires adjacent cysteine-rich regions to confer histone methyltransferase activity. Methylation of lysine 9 interferes with phosphorylation of serine 10, but is also influenced by pre-existing modifications in the amino terminus of H3. In vivo, deregulated SUV39H1 or disrupted Suv39h activity modulate H3 serine 10 phosphorylation in native chromatin and induce aberrant mitotic divisions. Our data reveal a functional interdependence of site-specific H3 tail modifications and suggest a dynamic mechanism for the regulation of higher-order chromatin.

More than 1,000 putative new human signalling proteins revealed by EST data mining

Cloning procedures aided by homology searches of EST databases have accelerated the pace of discovery of new genes, but EST database searching remains an involved and onerous task. More than 1.6 million human EST sequences have been deposited in public databases, making it difficult to identify ESTs that represent new genes. Compounding the problems of scale are difficulties in detection associated with a high sequencing error rate and low sequence similarity between distant homologues. We have developed a new method, coupling BLAST-based searches with a domain identification protocol, that filters candidate homologues. Application of this method in a large-scale analysis of 100 signalling domain families has led to the identification of ESTs representing more than 1,000 novel human signalling genes. The 4,206 publicly available ESTs representing these genes are a valuable resource for rapid cloning of novel human signalling proteins. For example, we were able to identify ESTs of at least 106 new small GTPases, of which 6 are likely to belong to new subfamilies. In some cases, further analyses of genomic DNA led to the discovery of previously unidentified full-length protein sequences. This is exemplified by the in silico cloning (prediction of a gene product sequence using only genomic and EST sequence data) of a new type of GTPase with two catalytic domains.

Homology-based method for identification of protein repeats using statistical significance estimates

Short protein repeats, frequently with a length between 20 and 40 residues, represent a significant fraction of known proteins. Many repeats appear to possess high amino acid substitution rates and thus recognition of repeat homologues is highly problematic. Even if the presence of a certain repeat family is known, the exact locations and the number of repetitive units often cannot be determined using current methods. We have devised an iterative algorithm based on optimal and sub-optimal score distributions from profile analysis that estimates the significance of all repeats that are detected in a single sequence. This procedure allows the identification of homologues at alignment scores lower than the highest optimal alignment score for non-homologous sequences. The method has been used to investigate the occurrence of eleven families of repeats in Saccharomyces cerevisiae, Caenorhabditis elegans and Homo sapiens accounting for 1055, 2205 and 2320 repeats, respectively. For these examples, the method is both more sensitive and more selective than conventional homology search procedures. The method allowed the detection in the SwissProt database of more than 2000 previously unrecognised repeats belonging to the 11 families. In addition, the method was used to merge several repeat families that previously were supposed to be distinct, indicating common phylogenetic origins for these families.

Phospholipase D regulation and localisation is dependent upon a phosphatidylinositol 4,5-biphosphate-specific PH domain

The signalling pathway leading, for example, to actin cytoskeletal reorganisation, secretion or superoxide generation involves phospholipase D (PLD)-catalysed hydrolysis of phosphatidylcholine to generate phosphatidic acid, which appears to mediate the messenger functions of this pathway. Two PLD genes (PLD1 and PLD2) with similar domain structures have been doned and progress has been made in identifying the protein regulators of PLD1 activation, for example Arf and Rho family members. The activities of both PLD isoforms are dependent on phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) and our sequence analysis suggested the presence of a pleckstrin homology (PH) domain in PLD1, although its absence has also been daimed. Investigation of the inositide dependence showed that a bis-phosphorylated lipid with a vicinal pair of phosphates was required for PLD1 activity. Furthermore, PLD1 bound specifically and with high affinity to lipid surfaces containing PI(4,5)P2 independently of the substrate phosphatidylcholine, suggesting a key role for the PH domain in PLD function. Importantly, a glutathione-S-transferase (GST) fusion protein comprising GST and the PH domain of PLD1 (GST-PLD1-PH) also bound specifically to supported lipid monolayers containing PI(4,5)P2. Point mutations within the PLD1 PH domain inhibited enzyme activity, whereas deletion of the domain both inhibited enzyme activity and disrupted normal PLD1 localisation. Thus, the functional PH domain regulates PLD by mediating its interaction with polyphosphoinositide-containing membranes; this might also induce a conformational change, thereby regulating catalytic activity.

SMART: a web-based tool for the study of genetically mobile domains

SMART (a Simple Modular Architecture Research Tool) allows the identification and annotation of genetically mobile domains and the analysis of domain architectures (http://SMART.embl-heidelberg.de ). More than 400 domain families found in signalling, extra-cellular and chromatin-associated proteins are detectable. These domains are extensively annotated with respect to phyletic distributions, functional class, tertiary structures and functionally important residues. Each domain found in a non-redundant protein database as well as search parameters and taxonomic information are stored in a relational database system. User interfaces to this database allow searches for proteins containing specific combinations of domains in defined taxa.

IMPALA: matching a protein sequence against a collection of PSI-BLAST-constructed position-specific score matrices

MOTIVATION

Many studies have shown that database searches using position-specific score matrices (PSSMs) or profiles as queries are more effective at identifying distant protein relationships than are searches that use simple sequences as queries. One popular program for constructing a PSSM and comparing it with a database of sequences is Position-Specific Iterated BLAST (PSI-BLAST).

RESULTS

This paper describes a new software package, IMPALA, designed for the complementary procedure of comparing a single query sequence with a database of PSI-BLAST-generated PSSMs. We illustrate the use of IMPALA to search a database of PSSMs for protein folds, and one for protein domains involved in signal transduction. IMPALA's sensitivity to distant biological relationships is very similar to that of PSI-BLAST. However, IMPALA employs a more refined analysis of statistical significance and, unlike PSI-BLAST, guarantees the output of the optimal local alignment by using the rigorous Smith-Waterman algorithm. Also, it is considerably faster when run with a large database of PSSMs than is BLAST or PSI-BLAST when run against the complete non-redundant protein database.

beta-propeller repeats and a PDZ domain in the tricorn protease: predicted self-compartmentalisation and C-terminal polypeptide-binding strategies of substrate selection

Prokaryotic proteases demonstrate a variety of substrate-selection strategies that prevent uncontrolled protein degradation. Proteasomes and ClpXP-like proteases form oligomeric structures that exclude large substrates from central solvated chambers containing their active sites. Monomeric prolyl oligopeptidases have been shown to contain beta-propeller structures that similarly reduce access to their catalytic residues. By contrast, Tsp-like enzymes contain PDZ domains that are thought to specifically target C-terminal polypeptides. We have investigated the sequence of Thermoplasma acidophilum tricorn protease using recently-developed database search methods. The tricorn protease is known to associate into a 20 hexamer capsid enclosing an extremely large cavity that is 37 nm in diameter. It is unknown, however, how this enzyme selects its small oligopeptide substrates. Our results demonstrate the presence in tricorn protease of a PDZ domain and two predicted six-bladed beta-propeller domains. We suggest that the PDZ domain is involved in targeting non-polar C-terminal peptides, similar to those generated by the T. acidophilum proteasome, whereas the beta-propeller domains serve to exclude large substrates from the tricorn protease active site in a similar manner to that previously indicated for prolyl oligopeptidase.

Raf-like Ras/Rap-binding domains in RGS12- and still-life-like signalling proteins

Ras proteins play critical roles in regulating cell growth and differentiation, and mutated Ras genes are expressed in a variety of human cancers. Consequently, much interest has centered on the binding partners of Ras, including the Ras-binding domain (RBD) of Raf kinase. Here evidence is presented that domains homologous to the Raf RBD are present in tandem in RGS12, RGS14 and LOCO, and singly in molecules similar to mouse Tiam-1. In addition, RGS12, RGS14 and LOCO are shown to contain single "LGN motifs" that are guanine nucleotide exchange factors specific for the alpha-subunit of G proteins. These findings indicate "cross-talk" interactions between signalling pathways involving Ras and Rap and pathways involving Rho, Rac and G alpha GTPases.

Analysis of mutations in the tudor domain of the survival motor neuron protein SMN

Autosomal recessive childhood onset spinal muscular atrophy (SMA) is a leading cause of infant mortality caused by mutations in the survival motor neuron (SMN) gene. The SMN protein is involved in RNA processing and is localised in structures called GEMs in the nucleus. Nothing is yet understood about why mutations in SMN gene result in the selective motor neuron loss observed in patients. The SMN protein domains conserved across several species may indicate functionally significant regions. Exon 3 of SMN contains homology to a tudor domain, where a Type I SMA patient has been reported to harbour a missense mutation. We have generated missense mutants in this region of SMN and have tested their ability to form GEMs when transfected into HeLa cells. Our results show such mutant SMN proteins still localise to GEMs. Furthermore, exon 7 deleted SMN protein appears to exert a dominant negative effect on localisation of endogenous SMN protein. However, exon 3 mutant protein and exon 5 deleted protein exert no such effect.

The cytoplasmic helical linker domain of receptor histidine kinase and methyl-accepting proteins is common to many prokaryotic signalling proteins

Mutations in the cytoplasmic linker regions of receptor histidine kinase and chemoreceptor proteins have been shown previously to significantly impair receptor functions. Here we demonstrate significant sequence similarities between these regions in numerous histidine kinases, methyl-accepting proteins, adenylyl cyclases and other prokaryotic signalling proteins. It is suggested that these 'HAMP domains' possess roles of regulating the phosphorylation or methylation of homodimeric receptors by transmitting the conformational changes in periplasmic ligand-binding domains to cytoplasmic signalling kinase and methyl-acceptor domains.

Eukaryotic signalling domain homologues in archaea and bacteria. Ancient ancestry and horizontal gene transfer

Phyletic distributions of eukaryotic signalling domains were studied using recently developed sensitive methods for protein sequence analysis, with an emphasis on the detection and accurate enumeration of homologues in bacteria and archaea. A major difference was found between the distributions of enzyme families that are typically found in all three divisions of cellular life and non-enzymatic domain families that are usually eukaryote-specific. Previously undetected bacterial homologues were identified for# plant pathogenesis-related proteins, Pad1, von Willebrand factor type A, src homology 3 and YWTD repeat-containing domains. Comparisons of the domain distributions in eukaryotes and prokaryotes enabled distinctions to be made between the domains originating prior to the last common ancestor of all known life forms and those apparently originating as consequences of horizontal gene transfer events. A number of transfers of signalling domains from eukaryotes to bacteria were confidently identified, in contrast to only a single case of apparent transfer from eukaryotes to archaea.

Protein families in multicellular organisms

The complete sequence of the nematode worm Caenorhabditis elegans contains the genetic machinery that is required to undertake the core biological processes of single cells. However, the genome also encodes proteins that are associated with multicellularity, as well as others that are lineage-specific expansions of phylogenetically widespread families and yet more that are absent in non-nematodes. Ongoing analysis is beginning to illuminate the similarities and differences among human proteins and proteins that are encoded by the genomes of the multicellular worm and the unicellular yeast, and will be essential in determining the reliability of transferring experimental data among phylogenetically distant species.