A new method for identification of protein (sub)families in a set of proteins based on hydropathy distribution in proteins

Structural similarity among proteins is reflected in the distribution of hydropathicity along the amino acids in the protein sequence. Similarities in the hydropathy distributions are obvious for homologous proteins within a protein family. They also were observed for proteins with related structures, even when sequence similarities were undetectable. Here we present a novel method that employs the hydropathy distribution in proteins for identification of (sub)families in a set of (homologous) proteins. We represent proteins as points in a generalized hydropathy space, represented by vectors of specifically defined features. The features are derived from hydropathy of the individual amino acids. Projection of this space onto principal axes reveals groups of proteins with related hydropathy distributions. The groups identified correspond well to families of structurally and functionally related proteins. We found that this method accurately identifies protein families in a set of proteins, or subfamilies in a set of homologous proteins. Our results show that protein families can be identified by the analysis of hydropathy distribution, without the need for sequence alignment. (C) 2005 Wiley-Liss, Inc.

Functional recycling of C2 domains throughout evolution: A comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level,the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various subfamilies. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces' our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics. (C) 2003 Elsevier Ltd. All rights reserved.

Probing the S100 protein family through genomic and functional analysis

The EF-hand superfamily of calcium binding proteins includes the S100, calcium binding protein, and troponin subfamilies. This study represents a genome, structure, and expression analysis of the S100 protein family, in mouse, human, and rat. We confirm the high level of conservation between mammalian sequences but show that four members, including S100A12, are present only in the human genome. We describe three new members of the S100 family in the three species and their locations within the S100 genomic clusters and propose a revised nomenclature and phylogenetic relationship between members of the EF-hand superfamily. Two of the three new genes were induced in bone-marrow-derived macrophages activated with bacterial lipopolysaccharide, suggesting a role in inflammation. Normal human and murine tissue distribution profiles indicate that some members of the family are expressed in a specific manner, whereas others are more ubiquitous. Structure-function analysis of the chemotactic properties of murine S100A8 and human S100A12, particularly within the active hinge domain, suggests that the human protein is the functional homolog of the murine protein. Strong similarities between the promoter regions of human S100A12 and murine S100A8 support this possibility. This study provides insights into the possible processes of evolution of the EF-hand protein superfamily. Evolution of the S100 proteins appears to have occurred in a modular fashion, also seen in other protein families such as the C2H2-type zinc-finger family. (C) 2004 Elsevier Inc. All rights reserved.

Systematic characterization of the zinc-finger-containing proteins in the mouse transcriptome

Zinc-finger-containing proteins can be classified into evolutionary and functionally divergent protein families that share one or more domains in which a zinc ion is tetrahedrally coordinated by cysteines and histidines. The zinc finger domain defines one of the largest protein superfamilies in mammalian genomes; 46 different conserved zinc finger domains are listed in InterPro (http://www.ebi.ac.uk/InterPro). Zinc finger proteins can bind to DNA, RNA, other proteins, or lipids as a modular domain in combination with other conserved structures. Owing to this combinatorial diversity, different members of zinc finger superfamilies contribute to many distinct cellular processes, including transcriptional regulation, mRNA stability and processing, and protein turnover. Accordingly, mutations of zinc finger genes lead to aberrations in a broad spectrum of biological processes such as development, differentiation, apoptosis, and immunological responses. This study provides the first comprehensive classification of zinc finger proteins in a mammalian transcriptome. Specific detailed analysis of the SP/Kruppel-like factors and the E3 ubiquitin-ligase RING-H2 families illustrates the importance of such an analysis for a more comprehensive functional classification of large protein families. We describe the characterization of a new family of C2H2 zinc-finger-containing proteins and a new conserved domain characteristic of this family, the identification and characterization of Sp8, a new member of the Sp family of transcriptional regulators, and the identification of five new RING-H2 proteins.

Domains of the TGN: Coats, tethers and G proteins

The trans-Golgi network is the major sorting compartment of the secretory pathway for protein, lipid and membrane traffic. There is a constant flow of membrane and cargo to and from this compartment. Evidence is emerging that the trans-Golgi network has multiple biochemically and functionally distinct subdomains, each of which contributes to the combined sorting and transport requirements of this dynamic compartment. The recruitment of distinct arrays of protein complexes to trans-Golgi network membranes is likely to produce the diversity of structure and biochemistry observed amongst subdomains that serve to generate different carriers or maintain resident trans-Golgi network components. This review discusses how these subdomains may be formed and examines the molecular players involved, including G proteins, clathrin adaptors and golgin tethers. Diversity within these protein families is highlighted and shown to be critical for the functionality of the trans-Golgi network, as a mediator of protein sorting and membrane transport, and for the maintenance of Golgi structure.

AAA ATPases regulate membrane association of yeast oxysterol binding proteins and sterol metabolism

The yeast genome encodes seven oxysterol binding protein homologs, Osh1p-Osh7p, which have been implicated in regulating intracellular lipid and vesicular transport. Here, we show that both Osh6p and Osh7p interact with Vps4p, a member of the AAA ( ATPases associated with a variety of cellular activities) family. The coiled-coil domain of Osh7p was found to interact with Vps4p in a yeast two-hybrid screen and the interaction between Osh7p and Vps4p appears to be regulated by ergosterol. Deletion of VPS4 induced a dramatic increase in the membrane-associated pools of Osh6p and Osh7p and also caused a decrease in sterol esterification, which was suppressed by overexpression of OSH7. Lastly, overexpression of the coiled-coil domain of Osh7p (Osh7pCC) resulted in a multi-vesicular body sorting defect, suggesting a dominant negative role of Osh7pCC possibly through inhibiting Vps4p function. Our data suggest that a common mechanism may exist for AAA proteins to regulate the membrane association of yeast OSBP proteins and that these two protein families may function together to control subcellular lipid transport.

Nuclear translocation of cell-surface receptors: Lessons from fibroblast growth factor

The nuclear localization of a number of growth factors, cytokine ligands and their receptors has been reported in various cell lines and tissues. These include members of the fibroblast growth factor (FGF), epidermal growth factor and growth hormone families. Accordingly, a number of nuclear functions have begun to emerge for these protein families. The demonstration of functional interactions of these proteins with the nuclear import machinery has further supported their functions as nuclear signal transducers. Here, we review the membrane- trafficking machinery and pathways demonstrated to regulate this cell surface to nucleus-trafficking event and highlight the many remaining unanswered questions. We focus on the FGF family, which is providing many of the clues as to the process of this unusual phenomenon.

Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a metachain to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely. [Bayesian phylogenetic inference; heating parameter; Markov chain Monte Carlo; replicated chains.]

A rapidly evolving secretome builds and patterns a sea shell

Background: Instructions to fabricate mineralized structures with distinct nanoscale architectures, such as seashells and coral and vertebrate skeletons, are encoded in the genomes of a wide variety of animals. In mollusks, the mantle is responsible for the extracellular production of the shell, directing the ordered biomineralization of CaCO3 and the deposition of architectural and color patterns. The evolutionary origins of the ability to synthesize calcified structures across various metazoan taxa remain obscure, with only a small number of protein families identified from molluskan shells. The recent sequencing of a wide range of metazoan genomes coupled with the analysis of gene expression in non-model animals has allowed us to investigate the evolution and process of biomineralization in gastropod mollusks. Results: Here we show that over 25% of the genes expressed in the mantle of the vetigastropod Haliotis asinina encode secreted proteins, indicating that hundreds of proteins are likely to be contributing to shell fabrication and patterning. Almost 85% of the secretome encodes novel proteins; remarkably, only 19% of these have identifiable homologues in the full genome of the patellogastropod Lottia scutum. The spatial expression profiles of mantle genes that belong to the secretome is restricted to discrete mantle zones, with each zone responsible for the fabrication of one of the structural layers of the shell. Patterned expression of a subset of genes along the length of the mantle is indicative of roles in shell ornamentation. For example, Has-sometsuke maps precisely to pigmentation patterns in the shell, providing the first case of a gene product to be involved in molluskan shell pigmentation. We also describe the expression of two novel genes involved in nacre (mother of pearl) deposition. Conclusion: The unexpected complexity and evolvability of this secretome and the modular design of the molluskan mantle enables diversification of shell strength and design, and as such must contribute to the variety of adaptive architectures and colors found in mollusk shells. The composition of this novel mantle-specific secretome suggests that there are significant molecular differences in the ways in which gastropods synthesize their shells.

Phosphoregulators: Protein kinases and protein phosphatases of mouse

With the completion of the human and mouse genome sequences, the task now turns to identifying their encoded transcripts and assigning gene function. In this study, we have undertaken a computational approach to identify and classify all of the protein kinases and phosphatases present in the mouse gene complement. A nonredundant set of these sequences was produced by mining Ensembl gene predictions and publicly available cDNA sequences with a panel of InterPro domains. This approach identified 561 candidate protein kinases and 162 candidate protein phosphatases. This cohort was then analyzed using TribeMCL protein sequence similarity clustering followed by CLUSTALV alignment and hierarchical tree generation. This approach allowed us to (1) distinguish between true members of the protein kinase and phosphatase families and enzymes of related biochemistry, (2) determine the structure of the families, and (3) suggest functions for previously uncharacterized members. The classifications obtained by this approach were in good agreement with previous schemes and allowed us to demonstrate domain associations with a number of clusters. Finally, we comment on the complementary nature of cDNA and genome-based gene detection and the impact of the FANTOM2 transcriptome project.

Low frequency of CHEK2 1100delC allele in Australian multiple-case breast cancer families: functional analysis in heterozygous individuals

A protein-truncating variant of CHEK2, 1100delC, is associated with a moderate increase in breast cancer risk. We have determined the prevalence of this allele in index cases from 300 Australian multiple-case breast cancer families, 95% of which had been found to be negative for mutations in BRCA1 and BRCA2. Only two (0.6%) index cases heterozygous for the CHEK2 mutation were identified. All available relatives in these two families were genotyped, but there was no evidence of co-segregation between the CHEK2 variant and breast cancer. Lymphoblastoid cell lines established from a heterozygous carrier contained approximately 20% of the CHEK2 1100delC mRNA relative to wild-type CHEK2 transcript. However, no truncated CHK2 protein was detectable. Analyses of expression and phosphorylation of wild-type CHK2 suggest that the variant is likely to act by haploinsufficiency. Analysis of CDC25A degradation, a downstream target of CHK2, suggests that some compensation occurs to allow normal degradation of CDC25A. Such compensation of the 1100delC defect in CHEK2 might explain the rather low breast cancer risk associated with the CHEK2 variant, compared to that associated with truncating mutations in BRCA1 or BRCA2.

Isolation, solution structure, and insecticidal activity of Kalata B2, a circular protein with a twist: Do Mobius strips exist in nature?

A large number of macrocyclic miniproteins with diverse biological activities have been isolated from the Rubiaceae, Violaceae, and Cucurbitaceae plant families in recent years. Here we report the three-dimensional structure determined using H-1 NMR spectroscopy and demonstrate potent insecticidal activity for one of these peptides, kalata B2. This peptide is one of the major components of an extract from the leaves of the plant Oldenlandia affinis. The structure consists of a distorted triple-stranded beta-sheet and a cystine knot arrangement of the disulfide bonds and is similar to those described for other members of the cyclotide family. The unique cyclic and knotted nature of these molecules makes them a fascinating example of topologically complex proteins. Examination of the sequences reveals that they can be separated into two subfamilies, one of which contains a larger number of positively charged residues and has a bracelet-like circularization of the backbone. The second subfamily contains a backbone twist due to a cis-peptidyl-proline bond and may conceptually be regarded as a molecular Mobius strip. Kalata B2 is the second putative member of the Mobius cyclotide family to be structurally characterized and has a cis-peptidyl-proline bond, thus validating the suggested name for this subfamily of cyclotides. The observation that kalata B2 inhibits the growth and development of Helicoverpa armigera larvae suggests a role for the cyclotides in plant defense. A comparison of the sequences and structures of kalata B1 and B2 provides insight into the biological activity of these peptides.

Discovery of cyclotide-like protein sequences in graminaceous crop plants: Ancestral precursors of circular proteins?

Cyclotides are peptides from plants of the Rubiaceae and Violaceae families that have the unusual characteristic of a macrocylic backbone. They are further characterized by their incorporation of a cystine knot in which two disulfides, along with the intervening backbone residues, form a ring through which a third disulfide is threaded. The cyclotides have been found in every Violaceae species screened to date but are apparently present in only a few Rubiaceae species. The selective distribution reported so far raises questions about the evolution of the cyclotides within the plant kingdom. In this study, we use a combined bioinformatics and expression analysis approach to elucidate the evolution and distribution of the cyclotides in the plant kingdom and report the discovery of related sequences widespread in the Poaceae family, including crop plants such as rice ( Oryza sativa), maize ( Zea mays), and wheat ( Triticum aestivum), which carry considerable economic and social importance. The presence of cyclotide-like sequences within these plants suggests that the cyclotides may be derived from an ancestral gene of great antiquity. Quantitative RT-PCR was used to show that two of the discovered cyclotide-like genes from rice and barley ( Hordeum vulgare) have tissue-specific expression patterns.

Localization of neogenin protein during morphogenesis in the mouse embryo

Neogenin, a close relative of the axon guidance receptor DCC, has been shown to be a receptor for members of the Netrin and Repulsive Guidance Molecule families. Recent studies have begun to uncover a role for Neogenin in organogenesis. Here we examine the localization of Neogenin protein in the developing mouse embryo (embryonic day 14.5) when organogenesis is progressing rapidly. We observe that Neogenin protein is restricted to distinct tissue layers within a given organ. In some embryonic epithelia such as the gut and pancreas, Neogenin protein is predominantly polarized to the basal surfaces of the epithelial cells. In contrast, Neogenin is restricted to mesenchymal cells within the lung and kidney. Neogenin is also seen in differentiating skeletal muscle and condensing cartilage throughout the embryo, and in the trigeminal and dorsal root ganglia of the peripheral nervous system. This study supports the emerging role for Neogenin as a key receptor in the establishment of the morphological architecture in many developing organ systems.