Evolution of the deaminase fold and multiple origins of eukaryotic editing and mutagenic nucleic acid deaminases from bacterial toxin systems

The deaminase-like fold includes, in addition to nucleic acid/nucleotide deaminases, several catalytic domains such as the JAB domain, and others involved in nucleotide and ADP-ribose metabolism. Using sensitive sequence and structural comparison methods, we develop a comprehensive natural classification of the deaminase-like fold and show that its ancestral version was likely to operate on nucleotides or nucleic acids. Consequently, we present evidence that a specific group of JAB domains are likely to possess a DNA repair function, distinct from the previously known deubiquitinating peptidase activity. We also identified numerous previously unknown clades of nucleic acid deaminases. Using inference based on contextual information, we suggest that most of these clades are toxin domains of two distinct classes of bacterial toxin systems, namely polymorphic toxins implicated in bacterial interstrain competition and those that target distantly related cells. Genome context information suggests that these toxins might be delivered via diverse secretory systems, such as Type V, Type VI, PVC and a novel PrsW-like intramembrane peptidase-dependent mechanism. We propose that certain deaminase toxins might be deployed by diverse extracellular and intracellular pathogens as also endosymbionts as effectors targeting nucleic acids of host cells. Our analysis suggests that these toxin deaminases have been acquired by eukaryotes on several independent occasions and recruited as organellar or nucleo-cytoplasmic RNA modifiers, operating on tRNAs, mRNAs and short non-coding RNAs, and also as mutators of hyper-variable genes, viruses and selfish elements. This scenario potentially explains the origin of mutagenic AID/APOBEC-like deaminases, including novel versions from Caenorhabditis, Nematostella and diverse algae and a large class of fast-evolving fungal deaminases. These observations greatly expand the distribution of possible unidentified mutagenic processes catalyzed by nucleic acid deaminases.     

Polymorphic toxin systems: comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics

ABSTRACT: BACKGROUND: Proteinaceous toxins are observed across all levels of inter-organismal and intra-genomic conflicts. These include recently discovered prokaryotic polymorphic toxin systems implicated in intra-specific conflict. They are characterized by a remarkable diversity of C-terminal toxin domains generated by recombination with standalone toxin-coding cassettes. Prior analysis revealed a striking diversity of nuclease and deaminase domains among the toxin modules. We systematically investigated polymorphic toxin systems using comparative genomics, sequence and structure analysis. RESULTS: Polymorphic toxin systems are distributed across all major bacterial lineages and are delivered by at least eight distinct secretory systems. In addition to type-II, these include type-V, VI, VII (ESX), and the poorly characterized "Photorhabdus virulence cassettes (PVC)", PrsW-dependent and MuF phage-capsid-like systems. We present evidence that trafficking of these toxins is often accompanied by autoproteolytic processing catalyzed by HINT, ZU5, PrsW, caspase-like, papain-like, and a novel metallopeptidase associated with the PVC system. We identified over 150 distinct toxin domains in these systems. These span an extraordinary catalytic spectrum to include 23 distinct clades of peptidases, numerous previously unrecognized versions of nucleases and deaminases, ADP-ribosyltransferases, ADP ribosyl cyclases, RelA/SpoT-like nucleotidyltransferases, glycosyltranferases and other enzymes predicted to modify lipids and carbohydrates, and a pore-forming toxin domain. Several of these toxin domains are shared with host-directed effectors of pathogenic bacteria. Over 90 families of immunity proteins might neutralize anywhere between a single to at least 27 distinct types of toxin domains. In some organisms multiple tandem immunity genes or immunity protein domains are organized into polyimmunity loci or polyimmunity proteins. Gene-neighborhood-analysis of polymorphic toxin systems predicts the presence of novel trafficking-related components, and also the organizational logic that allows toxin diversification through recombination. Domain architecture and protein-length analysis revealed that these toxins might be deployed as secreted factors, through directed injection, or via inter-cellular contact facilitated by filamentous structures formed by RHS/YD, filamentous hemagglutinin and other repeats. Phyletic pattern and life-style analysis indicate that polymorphic toxins and polyimmunity loci participate in cooperative behavior and facultative 'cheating' in several ecosystems such as the human oral cavity and soil. Multiple domains from these systems have also been repeatedly transferred to eukaryotes and their viruses, such as the nucleo-cytoplasmic large DNA viruses. CONCLUSIONS: Along with a comprehensive inventory of toxins and immunity proteins, we present several testable predictions regarding active sites and catalytic mechanisms of toxins, their processing and trafficking and their role in intra-specific and inter-specific interactions between bacteria. These systems provide insights regarding the emergence of key systems at different points in eukaryotic evolution, such as ADP ribosylation, interaction of myosin VI with cargo proteins, mediation of apoptosis, hyphal heteroincompatibility, hedgehog signaling, arthropod toxins, cell-cell interaction molecules like teneurins and different signaling messengers.     

The toxin/immunity network of Burkholderia pseudomallei contact-dependent growth inhibition (CDI) systems

Burkholderia pseudomallei is a category B pathogen and the causative agent of melioidosis – a serious infectious disease that is typically acquired directly from environmental reservoirs. Nearly all B. pseudomallei strains sequenced to date (> 85 isolates) contain gene clusters that are related to the contact‐dependent growth inhibition (CDI) systems of γ‐proteobacteria. CDI systems from Escherichia coli and Dickeya dadantii play significant roles in bacterial competition, suggesting these systems may also contribute to the competitive fitness of B. pseudomallei. Here, we identify 10 distinct CDI systems in B. pseudomallei based on polymorphisms within the cdiA‐CT/cdiI coding regions, which are predicted to encode CdiA‐CT/CdiI toxin/immunity protein pairs. Biochemical analysis of three B. pseudomallei CdiA‐CTs revealed that each protein possesses a distinct tRNase activity capable of inhibiting cell growth. These toxin activities are blocked by cognate CdiI immunity proteins, which specifically bind the CdiA‐CT and protect cells from growth inhibition. Using Burkholderia thailandensis E264 as a model, we show that a CDI system from B. pseudomallei 1026b mediates CDI and is capable of delivering CdiA‐CT toxins derived from other B. pseudomallei strains. These results demonstrate that Burkholderia species contain functional CDI systems, which may confer a competitive advantage to these bacteria.     

MOSC domains: ancient, predicted sulfur-carrier domains, present in diverse metal-sulfur cluster biosynthesis proteins including Molybdenum cofactor sulfurases

Using computational analysis, a novel superfamily of beta-strand-rich domains was identified in the Molybdenum cofactor sulfurase and several other proteins from both prokaryotes and eukaryotes. These MOSC domains contain an absolutely conserved cysteine and occur either as stand-alone forms such as the bacterial YiiM proteins, or fused to other domains such as a NifS-like catalytic domain in Molybdenum cofactor sulfurase. The MOSC domain is predicted to be a sulfur-carrier domain that receives sulfur abstracted by the pyridoxal phosphate-dependent NifS-like enzymes, on its conserved cysteine, and delivers it for the formation of diverse sulfur-metal clusters. The identification of this domain may clarify the mechanism of biogenesis of various metallo-enzymes including Molybdenum cofactor-containing enzymes that are compromised in human type II xanthinuria.     

Delivery of CdiA Nuclease Toxins into Target Cells during Contact-Dependent Growth Inhibition

Bacterial contact-dependent growth inhibition (CDI) is mediated by the CdiB/CdiA family of two-partner secretion proteins. CDI systems deploy a variety of distinct toxins, which are contained within the polymorphic C-terminal region (CdiA-CT) of CdiA proteins. Several CdiA-CTs are nucleases, suggesting that the toxins are transported into the target cell cytoplasm to interact with their substrates. To analyze CdiA transfer to target bacteria, we used the CDI system of uropathogenic Escherichia coli 536 (UPEC536) as a model. Antibodies recognizing the amino- and carboxyl-termini of CdiA^UPEC536 were used to visualize transfer of CdiA from CDI^UPEC536+ inhibitor cells to target cells using fluorescence microscopy. The results indicate that the entire CdiA^UPEC536 protein is deposited onto the surface of target bacteria. CdiA^UPEC536 transfer to bamA101 mutants is reduced, consistent with low expression of the CDI receptor BamA on these cells. Notably, our results indicate that the C-terminal CdiA-CT toxin region of CdiA^UPEC536 is translocated into target cells, but the N-terminal region remains at the cell surface based on protease sensitivity. These results suggest that the CdiA-CT toxin domain is cleaved from CdiA^UPEC536 prior to translocation. Delivery of a heterologous Dickeya dadantii CdiA-CT toxin, which has DNase activity, was also visualized. Following incubation with CDI⁺ inhibitor cells targets became anucleate, showing that the D.dadantii CdiA-CT was delivered intracellularly. Together, these results demonstrate that diverse CDI toxins are efficiently translocated across target cell envelopes.     

A superfamily of proteins involved in different secretion pathways in gram-negative bacteria: modular structure and specificity of the N-terminal domain

The family of PulD proteins, which has been characterized in a wide variety of microorganisms, comprises several membrane-associated proteins essential for the transport of macromolecules across bacterial membranes. These proteins are involved in the transport of complex structures (such as phage particles, DNA) or various proteins (such as extracellular enzymes and pathogenicity determinants). Amino acid sequence analysis revealed a possible modular organisation of proteins of this superfamily, with highly conserved C-terminal domains and dissimilar N-terminal domains. In the C-terminal domain, four highly conserved regions have been found, one of them containing a remarkable common motif: (V, I)PXL(S, G)XIPXXGXLF. Structural comparisons between the N-terminal domains indicate that proteins of this superfamily can be divided into at least two subgroups, probably reflecting the existence of distinct secretion mechanisms. This implies that members of the superfamily of PulD-related proteins are independently involved in (1) the general secretory pathway, (2) a new signal-peptide-independent secretion pathway found in several bacterial pathogens, and possibly in (3) the translocation of bacteriophage particles through the bacterial cell envelope.     

Identification of novel components of NAD-utilizing metabolic pathways and prediction of their biochemical functions

Nicotinamide adenine dinucleotide (NAD) is a ubiquitous cofactor participating in numerous redox reactions. It is also a substrate for regulatory modifications of proteins and nucleic acids via the addition of ADP-ribose moieties or removal of acyl groups by transfer to ADP-ribose. In this study, we use in-depth sequence, structure and genomic context analysis to uncover new enzymes and substrate-binding proteins in NAD-utilizing metabolic and macromolecular modification systems. We predict that Escherichia coli YbiA and related families of domains from diverse bacteria, eukaryotes, large DNA viruses and single strand RNA viruses are previously unrecognized components of NAD-utilizing pathways that probably operate on ADP-ribose derivatives. Using contextual analysis we show that some of these proteins potentially act in RNA repair, where NAD is used to remove 2'-3' cyclic phosphodiester linkages. Likewise, we predict that another family of YbiA-related enzymes is likely to comprise a novel NAD-dependent ADP-ribosylation system for proteins, in conjunction with a previously unrecognized ADP-ribosyltransferase. A similar ADP-ribosyltransferase is also coupled with MACRO or ADP-ribosylglycohydrolase domain proteins in other related systems, suggesting that all these novel systems are likely to comprise pairs of ADP-ribosylation and ribosylglycohydrolase enzymes analogous to the DraG-DraT system, and a novel group of bacterial polymorphic toxins. We present evidence that some of these coupled ADP-ribosyltransferases/ribosylglycohydrolases are likely to regulate certain restriction modification enzymes in bacteria. The ADP-ribosyltransferases found in these, the bacterial polymorphic toxin and host-directed toxin systems of bacteria such as Waddlia also throw light on the evolution of this fold and the origin of eukaryotic polyADP-ribosyltransferases and NEURL4-like ARTs, which might be involved in centrosomal assembly. We also infer a novel biosynthetic pathway that might be involved in the synthesis of a nicotinate-derived compound in conjunction with an asparagine synthetase and AMPylating peptide ligase. We use the data derived from this analysis to understand the origin and early evolutionary trajectories of key NAD-utilizing enzymes and present targets for future biochemical investigations.     

The proton‐motive force is required for translocation of CDI toxins across the inner membrane of target bacteria

Contact‐dependent growth inhibition (CDI) is a mode of bacterial competition orchestrated by the CdiB/CdiA family of two‐partner secretion proteins. The CdiA effector extends from the surface of CDI⁺ inhibitor cells, binds to receptors on neighbouring bacteria and delivers a toxin domain derived from its C‐terminal region (CdiA‐CT). Here, we show that CdiA‐CT toxin translocation requires the proton‐motive force (pmf) within target bacteria. The pmf is also critical for the translocation of colicin toxins, which exploit the energized Ton and Tol systems to cross the outer membrane. However, CdiA‐CT translocation is clearly distinct from known colicin‐import pathways because ΔtolA ΔtonB target cells are fully sensitive to CDI. Moreover, we provide evidence that CdiA‐CT toxins can be transferred into the periplasm of de‐energized target bacteria, indicating that transport across the outer membrane is independent of the pmf. Remarkably, CDI toxins transferred under de‐energized conditions remain competent to enter the target‐cell cytoplasm once the pmf is restored. Collectively, these results indicate that outer‐ and inner‐membrane translocation steps can be uncoupled, and that the pmf is required for CDI toxin transport from the periplasm to the target‐cell cytoplasm.     

The Cytoplasm-Entry Domain of Antibacterial CdiA Is a Dynamic α-Helical Bundle with Disulfide-Dependent Structural Features

Many Gram‐negative bacterial species use contact-dependent growth inhibition (CDI) systems to compete with neighboring cells. CDI⁺ strains express cell-surface CdiA effector proteins, which carry a toxic C-terminal region (CdiA-CT) that is cleaved from the effector upon transfer into the periplasm of target bacteria. The released CdiA-CT consists of two domains. The C-terminal domain is typically a nuclease that inhibits cell growth, and the N-terminal “cytoplasm-entry” domain mediates toxin translocation into the target-cell cytosol. Here, we use NMR and circular dichroism spectroscopic approaches to probe the structure, stability, and dynamics of the cytoplasm-entry domain from Escherichia coli STEC_MHI813. Chemical shift analysis reveals that the CdiA-CT^MHI813 entry domain is composed of a C-terminal helical bundle and a dynamic N-terminal region containing two disulfide linkages. Disruption of the disulfides by mutagenesis or chemical reduction destabilizes secondary structure over the N-terminus, but has no effect on the C-terminal helices. Although critical for N-terminal structure, the disulfides have only modest effects on global thermodynamic stability, and the entry domain exhibits characteristics of a molten globule. We find that the disulfides form in vivo as the entry domain dwells in the periplasm of inhibitor cells prior to target-cell recognition. CdiA-CT^MHI813 variants lacking either disulfide still kill target bacteria, but disruption of both bonds abrogates growth inhibition activity. We propose that the entry domain's dynamic structural features are critical for function. In its molten globule-like state, the domain resists degradation after delivery, yet remains pliable enough to unfold for membrane translocation.     

Cell surface anchorage and ligand-binding domains of the Saccharomyces cerevisiae cell adhesion protein alpha-agglutinin, a member of the immunoglobulin superfamily.

alpha-Agglutinin is a cell adhesion glycoprotein expressed on the cell wall of Saccharomyces cerevisiae alpha cells. Binding of alpha-agglutinin to its ligand a-agglutinin, expressed by a cells, mediates cell-cell contact during mating. Analysis of truncations of the 650-amino-acid alpha-agglutinin structural gene AG alpha 1 delineated functional domains of alpha-agglutinin. Removal of the C-terminal hydrophobic sequence allowed efficient secretion of the protein and loss of cell surface attachment. This cell surface anchorage domain was necessary for linkage to a glycosyl phosphatidylinositol anchor. A construct expressing the N-terminal 350 amino acid residues retained full a-agglutinin-binding activity, localizing the binding domain to the N-terminal portion of alpha-agglutinin. A 278-residue N-terminal peptide was inactive; therefore, the binding domain includes residues between 278 and 350. The segment of alpha-agglutinin between amino acid residues 217 and 308 showed significant structural and sequence similarity to a consensus sequence for immunoglobulin superfamily variable-type domains. The similarity of the alpha-agglutinin-binding domain to mammalian cell adhesion proteins suggests that this structure is a highly conserved feature of adhesion proteins in diverse eukaryotes.     

Coordinated delivery and function of bacterial MARTX toxin effectors

Bacteria often coordinate virulence factors to fine‐tune the host response during infection. These coordinated events can include toxins counteracting or amplifying effects of another toxin or though regulating the stability of virulence factors to remove their function once it is no longer needed. Multifunctional autoprocessing repeats‐in toxin (MARTX) toxins are effector delivery toxins that form a pore into the plasma membrane of a eukaryotic cell to deliver multiple effector proteins into the cytosol of the target cell. The function of these proteins includes manipulating actin cytoskeletal dynamics, regulating signal transduction pathways and inhibiting host secretory pathways. Investigations into the molecular mechanisms of these effector domains are providing insight into how the function of some effectors overlap and regulate one another during infection. Coordinated crosstalk of effector function suggests that MARTX toxins are not simply a sum of all their parts. Instead, modulation of cell function by effector domains may depend on which other effector domain are co‐delivered. Future studies will elucidate how these effectors interact with each other to modulate the bacterial host interaction.     

Prokaryotic structural maintenance of chromosomes (SMC) proteins: distribution, phylogeny, and comparison with MukBs and additional prokaryotic and eukaryotic coiled-coil proteins.

Structural maintenance of chromosomes (SMC) proteins are known to be essential for chromosome segregation in some prokaryotes and in eukaryotes. A systematic search for the distribution of SMC proteins in prokaryotes with fully or partially sequenced genomes showed that they form a larger family than previously anticipated and raised the number of known prokaryotic homologs to 54. Secondary structure predictions revealed that the length of the globular N-terminal and C-terminal domains is extremely well conserved in contrast to the hinge domain and coiled-coil domains which are considerably shorter in several bacterial species. SMC proteins are present in all gram-positive bacteria and in nearly all archaea while they were found in less than half of the gram-negative bacteria. Phylogenetic analyses indicate that the SMC tree roughly resembles the 16S rRNA tree, but that cyanobacteria and Aquifex aeolicus obtained smc genes by lateral transfer from archaea. Fourteen out of 22 smc genes located in fully sequenced genomes seem to be co-transcribed with a second gene out of six different gene families, indicating that the deduced gene products might be involved in similar functions. The SMC proteins were compared with other prokaryotic proteins with long coiled-coil domains. The lengths of different protein domains and signature sequences allowed to differentiate SMCs, MukBs, which were found to be confined to gamma proteobacteria, and two subfamilies of COG 0419 including the SbcC nuclease from E. coli. A phylogenetic analysis was performed including the prokaryotic coiled-coil proteins as well as SMCs and Rad18 proteins from selected eukaryotes.     

Functional and phylogenetic properties of the pore-forming beta-barrel transporters of the Omp85 family

beta-Barrel-shaped channels of the Omp85 family are involved in the translocation or assembly of proteins of bacterial, mitochondrial, and plastidic outer membranes. We have compared these proteins to understand the evolutionary development of the translocators. We have demonstrated that the proteins from proteobacteria and mitochondria have a pore diameter that is at least five times smaller than found for the Omp85 in cyanobacteria and plastids. This finding can explain why Omp85 from cyanobacteria (but not the homologous protein from proteobacteria) was remodeled to become the protein translocation pore after endosymbiosis. Further, the pore-forming region of the Omp85 proteins is restricted to the C terminus. Based on a phylogenetic analysis we have shown that the pore-forming domain displays a different evolutionary relationship than the N-terminal domain. In line with this, the affinity of the N-terminal domain to the C-terminal region of the Omp85 from plastids and cyanobacteria differs, even though the N-terminal domain is involved in gating of the pore in both groups. We have further shown that the N-terminal domain of nOmp85 takes part in homo-oligomerization. Thereby, the differences in the phylogeny of the two domains are explained by different functional constraints acting on the regions. The pore-forming domain, however, is further divided into two functional regions, where the distal C terminus itself forms a dimeric pore. Based on functional and phylogenetic analysis, we suggest an evolutionary scenario that explains the origin of the contemporary translocon.     

Autoprocessing of the Vibrio cholerae RTX toxin by the cysteine protease domain

Vibrio cholerae RTX is a large multifunctional bacterial toxin that causes actin crosslinking. Due to its size, it was predicted to undergo proteolytic cleavage during translocation into host cells to deliver activity domains to the cytosol. In this study, we identified a domain within the RTX toxin that is conserved in large clostridial glucosylating toxins TcdB, TcdA, TcnA, and TcsL; putative toxins from V. vulnificus, Yersinia sp., Photorhabdus sp., and Xenorhabdus sp.; and a filamentous/hemagglutinin-like protein FhaL from Bordetella sp. In vivo transfection studies and in vitro characterization of purified recombinant protein revealed that this domain from the V. cholerae RTX toxin is an autoprocessing cysteine protease whose activity is stimulated by the intracellular environment. A cysteine point mutation within the RTX holotoxin attenuated actin crosslinking activity suggesting that processing of the toxin is an important step in toxin translocation. Overall, we have uncovered a new mechanism by which large bacterial toxins and proteins deliver catalytic activities to the eukaryotic cell cytosol by autoprocessing after translocation.     

Secretory and endocytic pathways converge in a dynamic endosomal system in a primitive protozoan

Leishmania are a group of primitive eukaryotic trypanosomatid protozoa that are apically polarized with a flagellum at their anterior end. Surrounding the base of the flagellum is the flagellar reservoir that constitutes the site for endocytosis and exocytosis in these organisms. In the present study, we define a novel multivesicular tubular compartment involved in the intracellular trafficking of macromolecules in Leishmania . This dynamic structure appears to subtend the flagellar reservoir and extends towards the posterior end of the cell. Functional domains of several surface-expressed proteins, such as the gp63 glycosyl phosphatidyl inositol anchor and the 3'nucleotidase/nuclease transmembrane domain were fused to green fluorescent protein. These chimeric proteins were found to traffic through the secretory pathway and, while reaching their intended destinations, also accumulated within the intracellular tubular compartment. Using various compounds that are efficient fluid-phase markers used to track endocytosis in higher eukaryotes, we showed that this tubular compartment constitutes an important station in the endocytic pathway of these cells. Based on our functional observations of its role in the trafficking of expressed proteins and endocytosed markers, this compartment appears to have properties similar to endosomes of higher eukaryotes.     

Structural conservation of the PIN domain active site across all domains of life

The PIN (PilT N‐terminus) domain is a compact RNA‐binding protein domain present in all domains of life. This 120‐residue domain consists of a central and parallel β sheet surrounded by α helices, which together organize 4–5 acidic residues in an active site that binds one or more divalent metal ions and in many cases has endoribonuclease activity. In bacteria and archaea, the PIN domain is primarily associated with toxin–antitoxin loci, consisting of a toxin (the PIN domain nuclease) and an antitoxin that inhibits the function of the toxin under normal growth conditions. During nutritional or antibiotic stress, the antitoxin is proteolytically degraded causing activation of the PIN domain toxin leading to a dramatic reprogramming of cellular metabolism to cope with the new situation. In eukaryotes, PIN domains are commonly found as parts of larger proteins and are involved in a range of processes involving RNA cleavage, including ribosomal RNA biogenesis and nonsense‐mediated mRNA decay. In this review, we provide a comprehensive overview of the structural characteristics of the PIN domain and compare PIN domains from all domains of life in terms of structure, active site architecture, and activity.     

The monomeric GIY-YIG homing endonuclease I-BmoI uses a molecular anchor and a flexible tether to sequentially nick DNA

The GIY-YIG nuclease domain is found within protein scaffolds that participate in diverse cellular pathways and contains a single active site that hydrolyzes DNA by a one-metal ion mechanism. GIY-YIG homing endonucleases (GIY-HEs) are two-domain proteins with N-terminal GIY-YIG nuclease domains connected to C-terminal DNA-binding and they are thought to function as monomers. Using I-BmoI as a model GIY-HE, we test mechanisms by which the single active site is used to generate a double-strand break. We show that I-BmoI is partially disordered in the absence of substrate, and that the GIY-YIG domain alone has weak affinity for DNA. Significantly, we show that I-BmoI functions as a monomer at all steps of the reaction pathway and does not transiently dimerize or use sequential transesterification reactions to cleave substrate. Our results are consistent with the I-BmoI DNA-binding domain acting as a molecular anchor to tether the GIY-YIG domain to substrate, permitting rotation of the GIY-YIG domain to sequentially nick each DNA strand. These data highlight the mechanistic differences between monomeric GIY-HEs and dimeric or tetrameric GIY-YIG restriction enzymes, and they have implications for the use of the GIY-YIG domain in genome-editing applications.     

The emergence of catalytic and structural diversity within the beta-clip fold

The beta-clip fold includes a diverse group of protein domains that are unified by the presence of two characteristic waist-like constrictions, which bound a central extended region. Members of this fold include enzymes like deoxyuridine triphosphatase and the SET methylase, carbohydrate-binding domains like the fish antifreeze proteins/Sialate synthase C-terminal domains, and functionally enigmatic accessory subunits of urease and molybdopterin biosynthesis protein MoeA. In this study, we reconstruct the evolutionary history of this fold using sensitive sequence and structure comparisons methods. Using sequence profile searches, we identified novel versions of the beta-clip fold in the bacterial flagellar chaperone FlgA and the related pilus protein CpaB, the StrU-like dehydrogenases, and the UxaA/GarD-like hexuronate dehydratases (SAF superfamily). We present evidence that these versions of the beta-clip domain, like the related type III anti-freeze proteins and C-terminal domains of sialic acid synthases, are involved in interactions with carbohydrates. We propose that the FlgA and CpaB-like proteins mediate the assembly of bacterial flagella and Flp pili by means of their interactions with the carbohydrate moieties of peptidoglycan. The N-terminal beta-clip domain of the hexuronate dehydratases appears to have evolved a novel metal-binding site, while their C-terminal domain is likely to adopt a metal-binding TIM barrel-like fold. Using structural comparisons, we show that the beta-clip fold can be further classified into two major groups, one that includes the SAF, SET, dUTPase superfamilies, and the other that includes the phage lambda head decoration protein, the beta subunit of urease and the C-terminal domain of the molybdenum cofactor biosynthesis protein MoeA. Structural comparisons also suggest the beta-clip fold was assembled through the duplication of a three-stranded unit. Though the three-stranded units are likely to have had a common origin, we present evidence that complete beta-clip domains were assembled through such duplications, independently on multiple occasions. There is also evidence for circular permutation of the basic three-stranded unit on different occasions in the evolution of the beta-clip unit. We also describe how assembly of this fold from a basic three-stranded unit has been utilized to accommodate a variety of activities in its different versions.