Zinc through the three domains of life

Zinc is one of the metal ions essential for life, as it is required for the proper functioning of a large number of proteins. Despite its importance, the annotation of zinc-binding proteins in gene banks or protein domain databases still has significant room for improvement. In the present work, we compiled a list of known zinc-binding protein domains and of known zinc-binding sequence motifs (zinc-binding patterns), and then used them jointly to analyze the proteome of 57 different organisms to obtain an overview of zinc usage by archaeal, bacterial, and eukaryotic organisms. Zinc-binding proteins are an abundant fraction of these proteomes, ranging between 4% and 10%. The number of zinc-binding proteins correlates linearly with the total number of proteins encoded by the genome of an organism, but the proportionality constant of Eukaryota (8.8%) is significantly higher than that observed in Bacteria and Archaea (from 5% to 6%). Most of this enrichment is due to the larger portfolio of regulatory proteins in Eukaryota.     

A network-based approach to classify the three domains of life

Background

Chromosomal orthologs can reveal the shared ancestral gene set and their evolutionary trends. Additionally, physico-chemical properties of encoded proteins could provide information about functional adaptation and ecological niche requirements.

Results

We analyzed 7080 genes (five groups of 1416 orthologs each) from Rhizobiales species (S. meliloti, R. etli, and M. loti, plant symbionts; A. tumefaciens, a plant pathogen; and B. melitensis, an animal pathogen). We evaluated their phylogenetic relationships and observed three main topologies. The first, with closer association of R. etli to A. tumefaciens; the second with R. etli closer to S. meliloti; and the third with A. tumefaciens and S. meliloti as the closest pair. This was not unusual, given the close relatedness of these three species. We calculated the synonymous (dS) and nonsynonymous (dN) substitution rates of these orthologs, and found that informational and metabolic functions showed relatively low dN rates; in contrast, genes from hypothetical functions and cellular processes showed high dN rates. An alternative measure of sequence variability, percentage of changes by species, was used to evaluate the most specific proportion of amino acid residues from alignments. When dN was compared with that measure a high correlation was obtained, revealing that much of evolutive information was extracted with the percentage of changes by species at the amino acid level. By analyzing the sequence variability of orthologs with a set of five properties (polarity, electrostatic charge, formation of secondary structures, molecular volume, and amino acid composition), we found that physico-chemical characteristics of proteins correlated with specific functional roles, and association of species did not follow their typical phylogeny, probably reflecting more adaptation to their life styles and niche preferences. In addition, orthologs with low dN rates had residues with more positive values of polarity, volume and electrostatic charge.

Conclusions

These findings revealed that even when orthologs perform the same function in each genomic background, their sequences reveal important evolutionary tendencies and differences related to adaptation. This article was reviewed by: Dr. Purificación López-García, Prof. Jeffrey Townsend (nominated by Dr. J. Peter Gogarten), and Ms. Olga Kamneva.     

Multi-domain proteins in the three kingdoms of life: orphan domains and other unassigned regions

Comparative studies of the proteomes from different organisms have provided valuable information about protein domain distribution in the kingdoms of life. Earlier studies have been limited by the fact that only about 50% of the proteomes could be matched to a domain. Here, we have extended these studies by including less well-defined domain definitions, Pfam-B and clustered domains, MAS, in addition to Pfam-A and SCOP domains. It was found that a significant fraction of these domain families are homologous to Pfam-A or SCOP domains. Further, we show that all regions that do not match a Pfam-A or SCOP domain contain a significantly higher fraction of disordered structure. These unstructured regions may be contained within orphan domains or function as linkers between structured domains. Using several different definitions we have re-estimated the number of multi-domain proteins in different organisms and found that several methods all predict that eukaryotes have approximately 65% multi-domain proteins, while the prokaryotes consist of approximately 40% multi-domain proteins. However, these numbers are strongly dependent on the exact choice of cut-off for domains in unassigned regions. In conclusion, all eukaryotes have similar fractions of multi-domain proteins and disorder, whereas a high fraction of repeating domain is distinguished only in multicellular eukaryotes. This implies a role for repeats in cell-cell contacts while the other two features are important for intracellular functions.     

A bioinformatic approach to the identification of bacterial proteins interacting with Toll-interleukin 1 receptor-resistance (TIR) homology domains

Members of the Toll-like receptor (TLR) family are currently under intense scrutiny for their role in the sampling and recognition of pathogens. It has already been reported that both vaccinia virus and Yersinia spp. express proteins that help them evade the TLR mediated immune response, acting through the Toll-interleukin-1 receptor-resistance (TIR) domain and leucine-rich repeat region of the host TLRs respectively. The TIR domain is involved in the dimerisation of the TLRs and their complexation with their adapter molecules. We tested here the hypothesis that bacteria have the ability to secrete proteins containing similar motifs to the intracellular TIR domains that are involved in the TIR-TIR interaction necessary for the subsequent signal transmission. Based upon their sequence homology, proteins expressing TIRs have been divided into three sub-classes, based around the TLRs, the TLR adapter proteins, and the interleukin-1 and -18 adapter proteins. The highly conserved regions from these separate sub-families were then used to identify similar bacterial proteins. The bacterial proteins identified were then included in an iterative MEME-BLAST process to broaden the search. Tollip, a known TLR antagonist and adapter protein, was included in this investigation although it does not fit into any of the three sub-classes outlined above. If suitable bacterial proteins had been identified, it would signify that certain bacteria had evolved a mechanism to aid them in avoiding detection by the innate immune system acting through the TIR domains. At this stage one has to conclude that there is no evidence currently available suggesting such a mechanism, when using the strategy applied here.     

A docking approach to the study of copper trafficking proteins; interaction between metallochaperones and soluble domains of copper ATPases

A structural model of the transient complex between the yeast copper chaperone Atx1 and the first soluble domain of the copper transporting ATPase Ccc2 was obtained with HADDOCK, combining NMR chemical shift mapping information with in silico docking. These two proteins are involved in copper trafficking in yeast cells. Calculations were performed starting with the copper ion either bound to Atx1 or to Ccc2 and using the experimental structures of the copper-loaded and apo forms of each protein. The copper binding motifs of the two proteins are found in close proximity. Copper tends to move from Atx1 to Ccc2, consistent with the physiological direction of transfer, with concomitant structural rearrangements, in agreement with experimental observations. The interaction is mainly of an electrostatic nature with hydrogen bonds stabilizing the complex. The structural data are relevant for a number of proteins homologous to Atx1 and Ccc2 and conserved from bacteria to humans.     

Origins of DNA replication in the three domains of life

Replication of DNA is essential for the propagation of life. It is somewhat surprising then that, despite the vital nature of this process, cellular organisms show a great deal of variety in the mechanisms that they employ to ensure appropriate genome duplication. This diversity is manifested along classical evolutionary lines, with distinct combinations of replicon architecture and replication proteins being found in the three domains of life: the Bacteria, the Eukarya and the Archaea. Furthermore, although there are mechanistic parallels, even within a given domain of life, the way origins of replication are defined shows remarkable variation.     

Estimation of phylogenetic inconsistencies in the three domains of life

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.     

Dehydrogenases from all three domains of life cleave RNA

Specific interactions of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) with RNA have been reported both in vitro and in vivo. We show that eukaryotic and bacterial GAPDH and two proteins from the hyperthermophilic archaeon Sulfolobus solfataricus, which are annotated as dehydrogenases, cleave RNA producing similar degradation patterns. RNA cleavage is most efficient at 60 degrees C, at MgCl(2) concentrations up to 5 mm, and takes place between pyrimidine and adenosine. The RNase active center of the putative aspartate semialdehyde dehydrogenase from S. solfataricus is located within the N-terminal 73 amino acids, which comprise the first mononucleotide-binding site of the predicted Rossmann fold. Thus, RNA cleavage has to be taken into account in the ongoing discussion of the possible biological function of RNA binding by dehydrogenases.     

Enhancing the mechanical stability of proteins through a cocktail approach

Rationally enhancing the mechanical stability of proteins remains a challenge in the field of single molecule force spectroscopy. Here we demonstrate that it is feasible to use a “cocktail” approach for combining more than one approach to enhance significantly the mechanical stability of proteins in an additive fashion. As a proof of principle, we show that metal chelation and protein-protein interaction can be combined to enhance the unfolding force of a protein to ∼450 pN, which is >3 times of its original value. This is also higher than the mechanical stability of most of proteins studied so far. We also extend such a cocktail concept to combine two different metal chelation sites to enhance protein mechanical stability. This approach opens new avenues to efficiently regulating the mechanical properties of proteins, and should be applicable to a wide range of elastomeric proteins.     

Prediction of food protein allergenicity: a bioinformatic learning systems approach

Food hypersensitivity is constantly increasing in Western societies with a prevalence of about 1-2% in Europe and in the USA. Among children, the incidence is even higher. Because of the introduction of foods derived from genetically modified crops on the marketplace, the scientific community, regulatory bodies and international associations have intensified discussions on risk assessment procedures to identify potential food allergenicity of the newly introduced proteins. In this work, we present a novel biocomputational methodology for the classification of amino acid sequences with regard to food allergenicity and non-allergenicity. This method relies on a computerised learning system trained using selected excerpts of amino acid sequences. One example of such a successful learning system is presented which consists of feature extraction from sequence alignments performed with the FASTA3 algorithm (employing the BLOSUM50 substitution matrix) combined with the k-Nearest-Neighbour (kNN) classification algorithm. Briefly, the two features extracted are the alignment score and the alignment length and the kNN algorithm assigns the pair of extracted features from an unknown sequence to the prevalent class among its k nearest neighbours in the training (prototype) set available. 91 food allergens from several specialised public repositories of food allergy and the SWALL database were identified, pre-processed, and stored, yielding one of the most extensively characterised repositories of allergenic sequences known today. All allergenic sequences were classified using a standard one-leave-out cross validation procedure yielding about 81% correctly classified allergens and the classification of 367 non-allergens in an independent test set resulted in about 98% correct classifications. The biocomputational approach presented should be regarded as a significant extension and refinement of earlier attempts suggested for in silico food safety assessment. Our results show that the framework described here is powerful enough to become useful as part of a multiple-procedure test scheme that also depicts other evaluation approaches such as solid phase immunoassay and tests for stability to digestions.     

Cladogenesis, coalescence and the evolution of the three domains of life

In this article, we explore the large-scale structure of the tree of life by using a simple model with a constant number of species and rates of speciation that equal the rates of extinction. In addition, we discuss the consequences of horizontal gene transfer for the concept of a most recent common ancestor of all living organisms (cenancestor). A simple null hypothesis based on coalescence theory explains some features of the observed topologies of the tree of life. Simulations of genes and organismal lineages suggest that there was no single common ancestor that contained all the genes ancestral to those shared among the three domains of life. Each contemporary molecule has its own history that traces back to an individual molecular cenancestor. However, these molecular ancestors were likely to be present in different organisms and at different times.     

Noncellulosomal cohesin- and dockerin-like modules in the three domains of life

The high-affinity cohesin–dockerin interaction was originally discovered as modular components, which mediate the assembly of the various subunits of the multienzyme cellulosome complex that characterizes some cellulolytic bacteria. Until recently, the presence of cohesins and dockerins within a bacterial proteome was considered a definitive signature of a cellulosome-producing bacterium. Widespread genome sequencing has since revealed a wealth of putative cohesin- and dockerin-containing proteins in Bacteria, Archaea, and in primitive eukaryotes. The newly identified modules appear to serve diverse functions that are clearly distinct from the classical cellulosome archetype, and the vast majority of parent proteins are not predicted glycoside hydrolases. In most cases, only a few such genes have been identified in a given microorganism, which encode proteins containing but a single cohesin and/or dockerin. In some cases, one or the other module appears to be missing from a given species, and in other cases both modules occur within the same protein. This review provides a bioinformatics-based survey of the current status of cohesin- and dockerin-like sequences in species from the Bacteria, Archaea, and Eukarya. Surprisingly, many identified modules and their parent proteins are clearly unrelated to cellulosomes. The cellulosome paradigm may thus be the exception rather than the rule for bacterial, archaeal, and eukaryotic employment of cohesin and dockerin modules.     

Transitional forms between the three domains of life and evolutionary implications

The question as to the origin and relationship between the three domains of life is lodged in a phylogenetic impasse. The dominant paradigm is to see the three domains as separated. However, the recently characterized bacterial species have suggested continuity between the three domains. Here, we review the evidence in support of this hypothesis and evaluate the implications for and against the models of the origin of the three domains of life. The existence of intermediate steps between the three domains discards the need for fusion to explain eukaryogenesis and suggests that the last universal common ancestor was complex. We propose a scenario in which the ancestor of the current bacterial Planctomycetes, Verrucomicrobiae and Chlamydiae superphylum was related to the last archaeal and eukaryotic common ancestor, thus providing a way out of the phylogenetic impasse.     

The early evolution of lipid membranes and the three domains of life

All cell membranes are composed of glycerol phosphate phospholipids, and this commonality argues for the presence of such phospholipids in the last common ancestor, or cenancestor. However, phospholipid biosynthesis is very different between bacteria and archaea, leading to the suggestion that the cenancestor was devoid of phospholipid membranes. Recent phylogenomic studies challenge this view, suggesting that the cenancestor did possess complex phospholipid membranes. Here, we discuss the implications of these recent findings for membrane evolution in archaea and bacteria, and for the origin of the eukaryotic cell.     

Rigid domains in proteins: An algorithmic approach to their identification

A rigid domain, defined here as a tertiary structure common to two or more different protein conformations, can be identified numerically from atomic coordinates by finding sets of residues, one in each conformation, such that the distance between any two residues within the set belonging to one conformation is the same as the distance between the two structurally equivalent residues within the set belonging to any other conformation. The distance between two residues is taken to be the distance between their respective α carbon atoms. With the methods of this paper we have found in the deoxy and oxy conformations of the human hemoglobin α₁β₁ dimer a rigid domain closely related to that previously identified by Baldwin and Chothia (J. Mol. Biol. 129:175–220,1979). We provide two algorithms, both using the difference-distance matrix, with which to search for rigid domains directly from atomic coordinates. The first finds all rigid domains in a protein but has storage and processing demands that become prohibitively large with increasing protein size. The second, although not necessarily finding every rigid domain, is computationally tractable for proteins of any size. Because of its efficiency we are able to search protein conformations recursively for groups of non-intersecting domains. Different protein conformations, when aligned by superimposing their respective domain structures; can be examined for structural differences in regions complementing a rigid domain. © 1995 Wiley-Liss, Inc.     

The replication clamp-loading machine at work in the three domains of life

Sliding clamps are ring-shaped proteins that tether DNA polymerases to DNA, which enables the rapid and processive synthesis of both leading and lagging strands at the replication fork. The clamp-loading machinery must repeatedly load sliding-clamp factors onto primed sites at the replication fork. Recent structural and biochemical analyses provide unique insights into how these clamp-loading ATPase machines function to load clamps onto the DNA. Moreover, these studies highlight the evolutionary conservation of the clamp-loading process in the three domains of life.