A novel family of P-loop NTPases with an unusual phyletic distribution and transmembrane segments inserted within the NTPase domain

Background

Recent sequence-structure studies on P-loop-fold NTPases have substantially advanced the existing understanding of their evolution and functional diversity. These studies provide a framework for characterization of novel lineages within this fold and prediction of their functional properties.

Results

Using sequence profile searches and homology-based structure prediction, we have identified a previously uncharacterized family of P-loop NTPases, which includes the neuronal membrane protein and receptor tyrosine kinase substrate Kidins220/ARMS, which is conserved in animals, the F-plasmid PifA protein involved in phage T7 exclusion, and several uncharacterized bacterial proteins. We refer to these (predicted) NTPases as the KAP family, after Kidins220/ARMS and PifA. The KAP family NTPases are sporadically distributed across a wide phylogenetic range in bacteria but among the eukaryotes are represented only in animals. Many of the prokaryotic KAP NTPases are encoded in plasmids and tend to undergo disruption to form pseudogenes. A unique feature of all eukaryotic and certain bacterial KAP NTPases is the presence of two or four transmembrane helices inserted into the P-loop NTPase domain. These transmembrane helices anchor KAP NTPases in the membrane such that the P-loop domain is located on the intracellular side. We show that the KAP family belongs to the same major division of the P-loop NTPase fold with the AAA+, ABC, RecA-like, VirD4-like, PilT-like, and AP/NACHT-like NTPase classes. In addition to the KAP family, we identified another small family of predicted bacterial NTPases, with two transmembrane helices inserted into the P-loop domain. This family is not specifically related to the KAP NTPases, suggesting independent acquisition of the transmembrane helices.

Conclusions

We predict that KAP family NTPases function principally in the NTP-dependent dynamics of protein complexes, especially those associated with the intracellular surface of cell membranes. Animal KAP NTPases, including Kidins220/ARMS, are likely to function as NTP-dependent regulators of the assembly of membrane-associated signaling complexes involved in neurite growth and development. One possible function of the prokaryotic KAP NTPases might be in the exclusion of selfish replicons, such as viruses, from the host cells. Phylogenetic analysis and phyletic patterns suggest that the common ancestor of the animals acquired a KAP NTPase via lateral transfer from bacteria. However, an earlier transfer into eukaryotes followed by multiple losses in several eukaryotic lineages cannot be ruled out.     

Evolutionary relationship between the bacterial HPr kinase and the ubiquitous PEP-carboxykinase: expanding the P-loop nucleotidyl transferase superfamily

Similarities between protein three-dimensional structures can reveal evolutionary and functional relationships not apparent from sequence comparison alone. Here we report such a similarity between the metabolic enzymes histidine phosphocarrier protein kinase (HPrK) and phosphoenolpyruvate carboxykinase (PCK), suggesting that they are evolutionarily related. Current structure classifications place PCK and other P-loop containing nucleotidyl-transferases into different folds. Our comparison of both HPrK and PCK to other P-loop containing proteins reveals that all share a common structural motif consisting of an alphabeta segment containing the P-loop flanked by an additional beta-strand that is adjacent in space, but far apart along the sequence. Analysis also shows that HPrK/PCK differ from other P-loop containing structures no more than they differ from each other. We thus suggest that HPrK and PCK should be classified with other P-loop containing proteins, and that all probably share a common ancestor that probably contained a simple P-loop motif with different protein segments being added or lost over the course of evolution. We used the structure-based sequence alignment containing residues specific to HPrK/PCK to identify additional members of this P-loop containing family.     

Globins in nonvertebrate species: dispersal by horizontal gene transfer and evolution of the structure-function relationships

Using a new template based on an alignment of 145 nonvertebrate globins we examined several recently determined sequences of putative globins and globin-like hemeproteins. We propose that all globins have evolved from a family of ancestral, approx. 17-kDa hemeproteins, which displayed the globin fold and functioned as redox proteins. Once atmospheric O2 became available the acquisition of oxygen-binding properties was initiated, culminating in the various highly specialized functions known as present. During this evolutionary process, we suggest that (1) high oxygen affinity may have been acquired repeatedly and (2) the formation of chimeric proteins containing both a globin and a flavin binding domain was an additional and distinct evolutionary trend. Furthermore, globin-like hemeproteins encompass hemeproteins produced through convergent evolution from nonglobin ancestral proteins to carry out O2-binding functions as well as hemeproteins whose sequences exhibit the loss of some or all of the structural determinants of the globin fold. We also propose that there occurred two cases of horizontal globin gene transfer, one from an ancestor common to the ciliates Paramecium and Tetrahymena and the green alga Chlamydomonas to a cyanobacterium ancestor and the other, from a eukaryote ancestor of the yeasts Saccharomyces and Candida to a bacterial ancestor of the proteobacterial genera Escherichia, Alcaligenes, and Vitreoscilla.      

Structural evolution of the protein kinase-like superfamily

The protein kinase family is large and important, but it is only one family in a larger superfamily of homologous kinases that phosphorylate a variety of substrates and play important roles in all three superkingdoms of life. We used a carefully constructed structural alignment of selected kinases as the basis for a study of the structural evolution of the protein kinase-like superfamily. The comparison of structures revealed a "universal core" domain consisting only of regions required for ATP binding and the phosphotransfer reaction. Remarkably, even within the universal core some kinase structures display notable changes, while still retaining essential activity. Hence, the protein kinase-like superfamily has undergone substantial structural and sequence revision over long evolutionary timescales. We constructed a phylogenetic tree for the superfamily using a novel approach that allowed for the combination of sequence and structure information into a unified quantitative analysis. When considered against the backdrop of species distribution and other metrics, our tree provides a compelling scenario for the development of the various kinase families from a shared common ancestor. We propose that most of the so-called "atypical kinases" are not intermittently derived from protein kinases, but rather diverged early in evolution to form a distinct phyletic group. Within the atypical kinases, the aminoglycoside and choline kinase families appear to share the closest relationship. These two families in turn appear to be the most closely related to the protein kinase family. In addition, our analysis suggests that the actin-fragmin kinase, an atypical protein kinase, is more closely related to the phosphoinositide-3 kinase family than to the protein kinase family. The two most divergent families, alpha-kinases and phosphatidylinositol phosphate kinases (PIPKs), appear to have distinct evolutionary histories. While the PIPKs probably have an evolutionary relationship with the rest of the kinase superfamily, the relationship appears to be very distant (and perhaps indirect). Conversely, the alpha-kinases appear to be an exception to the scenario of early divergence for the atypical kinases: they apparently arose relatively recently in eukaryotes. We present possible scenarios for the derivation of the alpha-kinases from an extant kinase fold.     

Convergent evolution of similar enzymatic function on different protein folds: the hexokinase, ribokinase, and galactokinase families of sugar kinases.

Kinases that catalyze phosphorylation of sugars, called here sugar kinases, can be divided into at least three distinct nonhomologous families. The first is the hexokinase family, which contains many prokaryotic and eukaryotic sugar kinases with diverse specificities, including a new member, rhamnokinase from Salmonella typhimurium. The three-dimensional structure of hexokinase is known and can be used to build models of functionally important regions of other kinases in this family. The second is the ribokinase family, of unknown three-dimensional structure, and comprises pro- and eukaryotic ribokinases, bacterial fructokinases, the minor 6-phosphofructokinase 2 from Escherichia coli, 6-phosphotagatokinase, 1-phosphofructokinase, and, possibly, inosine-guanosine kinase. The third family, also of unknown three-dimensional structure, contains several bacterial and yeast galactokinases and eukaryotic mevalonate and phosphomevalonate kinases and may have a substrate binding region in common with homoserine kinases. Each of the three families of sugar kinases appears to have a distinct three-dimensional fold, since conserved sequence patterns are strikingly different for the three families. Yet each catalyzes chemically equivalent reactions on similar or identical substrates. The enzymatic function of sugar phosphorylation appears to have evolved independently on the three distinct structural frameworks, by convergent evolution. In addition, evolutionary trees reveal that (1) fructokinase specificity has evolved independently in both the hexokinase and ribokinase families and (2) glucose specificity has evolved independently in different branches of the hexokinase family. These are examples of independent Darwinian adaptation of a structure to the same substrate at different evolutionary times. The flexible combination of active sites and three-dimensional folds observed in nature can be exploited by protein engineers in designing and optimizing enzymatic function.     

Classification and evolution of P-loop GTPases and related ATPases

Sequences and available structures were compared for all the widely distributed representatives of the P-loop GTPases and GTPase-related proteins with the aim of constructing an evolutionary classification for this superclass of proteins and reconstructing the principal events in their evolution. The GTPase superclass can be divided into two large classes, each of which has a unique set of sequence and structural signatures (synapomorphies). The first class, designated TRAFAC (after translation factors) includes enzymes involved in translation (initiation, elongation, and release factors), signal transduction (in particular, the extended Ras-like family), cell motility, and intracellular transport. The second class, designated SIMIBI (after signal recognition particle, MinD, and BioD), consists of signal recognition particle (SRP) GTPases, the assemblage of MinD-like ATPases, which are involved in protein localization, chromosome partitioning, and membrane transport, and a group of metabolic enzymes with kinase or related phosphate transferase activity. These two classes together contain over 20 distinct families that are further subdivided into 57 subfamilies (ancient lineages) on the basis of conserved sequence motifs, shared structural features, and domain architectures. Ten subfamilies show a universal phyletic distribution compatible with presence in the last universal common ancestor of the extant life forms (LUCA). These include four translation factors, two OBG-like GTPases, the YawG/YlqF-like GTPases (these two subfamilies also consist of predicted translation factors), the two signal-recognition-associated GTPases, and the MRP subfamily of MinD-like ATPases. The distribution of nucleotide specificity among the proteins of the GTPase superclass indicates that the common ancestor of the entire superclass was a GTPase and that a secondary switch to ATPase activity has occurred on several independent occasions during evolution. The functions of most GTPases that are traceable to LUCA are associated with translation. However, in contrast to other superclasses of P-loop NTPases (RecA-F1/F0, AAA+, helicases, ABC), GTPases do not participate in NTP-dependent nucleic acid unwinding and reorganizing activities. Hence, we hypothesize that the ancestral GTPase was an enzyme with a generic regulatory role in translation, with subsequent diversification resulting in acquisition of diverse functions in transport, protein trafficking, and signaling. In addition to the classification of previously known families of GTPases and related ATPases, we introduce several previously undetected families and describe new functional predictions.     

A CTP-dependent archaeal riboflavin kinase forms a bridge in the evolution of cradle-loop barrels

Proteins of the cradle-loop barrel metafold are formed by duplication of a conserved betaalphabeta-element, suggesting a common evolutionary origin from an ancestral group of nucleic acid-binding proteins. The basal fold within this metafold, the RIFT barrel, is also found in a wide range of enzymes, whose homologous relationship with the nucleic acid-binding group is unclear. We have characterized a protein family that is intermediate in sequence and structure between the basal group of cradle-loop barrels and one family of RIFT-barrel enzymes, the riboflavin kinases. We report the structure, substrate-binding mode, and catalytic activity for one of these proteins, Methanocaldococcus jannaschii Mj0056, which is an archaeal riboflavin kinase. Mj0056 is unusual in utilizing CTP rather than ATP as the donor nucleotide, and sequence conservation in the relevant residues suggests that this is a general feature of archaeal riboflavin kinases.     

Analysis of P-Loop and its Flanking Region Subsequence of Diverse NTPases Reveals Evolutionary Selected Residues

The P-loop NTPases are involved in diverse cellular functions. Members of the P-loop NTPase superfamily are characterized by presence of a highly conserved sequence pattern GxxxxGKS/T, known as Walker A motif. This motif adopts an archetypal P-loop conformation which allows accommodation of the triphosphate moiety of a bound nucleotide. Despite the presence of Walker A as a common sequence motif, P-loop NTPases exhibit extreme sequence divergence which hampers their phylogenetic or evolutionary classification. Here, we show that P-loop and its flanking region subsequence (termed as “extended-WalkerA motif”) contain distinct signatures that can be utilized to classify NTPase domain of functionally diverse proteins. We find a clearly classified group of diverse NTPases of Conserved Domain Database such as G-proteins, Ylqf, RecA like, DExDc, AAA, CPT, NK, ABC transporter and NifH proteins.     

α/β Barrel evolution and the modular assembly of enzymes: Emerging trends in the flavin oxidase/dehydrogenase family

α/β barrels have an ill-defined origin. Evidence exists which favours their divergent evolution from a common ancestral barrel and convergent evolution to a stable fold. However, recent sequence and structural information for the flavin oxidase/dehydrogenase family of barrel enzymes indicate that sub-families of α/β barrels have evolved divergently. The modular fusion of barrel domains with core structures from other gene families has also contributed to the evolution of related but catalytically distinct enzyme molecules within each sub-family of the flavin oxidases/dehydrogenases. An analysis of the structures and sequences of the flavin oxidases/dehydrogenases has now enabled studies focusing on the evolutionary origins and modular assembly of this important family of proteins to be initiated.     

Structural Evolution of the Protein Kinase–Like Superfamily

The protein kinase family is large and important, but it is only one family in a larger superfamily of homologous kinases that phosphorylate a variety of substrates and play important roles in all three superkingdoms of life. We used a carefully constructed structural alignment of selected kinases as the basis for a study of the structural evolution of the protein kinase–like superfamily. The comparison of structures revealed a “universal core” domain consisting only of regions required for ATP binding and the phosphotransfer reaction. Remarkably, even within the universal core some kinase structures display notable changes, while still retaining essential activity. Hence, the protein kinase–like superfamily has undergone substantial structural and sequence revision over long evolutionary timescales. We constructed a phylogenetic tree for the superfamily using a novel approach that allowed for the combination of sequence and structure information into a unified quantitative analysis. When considered against the backdrop of species distribution and other metrics, our tree provides a compelling scenario for the development of the various kinase families from a shared common ancestor. We propose that most of the so-called “atypical kinases” are not intermittently derived from protein kinases, but rather diverged early in evolution to form a distinct phyletic group. Within the atypical kinases, the aminoglycoside and choline kinase families appear to share the closest relationship. These two families in turn appear to be the most closely related to the protein kinase family. In addition, our analysis suggests that the actin-fragmin kinase, an atypical protein kinase, is more closely related to the phosphoinositide-3 kinase family than to the protein kinase family. The two most divergent families, α-kinases and phosphatidylinositol phosphate kinases (PIPKs), appear to have distinct evolutionary histories. While the PIPKs probably have an evolutionary relationship with the rest of the kinase superfamily, the relationship appears to be very distant (and perhaps indirect). Conversely, the α-kinases appear to be an exception to the scenario of early divergence for the atypical kinases: they apparently arose relatively recently in eukaryotes. We present possible scenarios for the derivation of the α-kinases from an extant kinase fold.     

The origins of modern proteomes

Distributions of phylogenetically related protein domains (fold superfamilies), or FSFs, among the three Superkingdoms (trichotomy) are assessed. Very nearly 900 of the 1200 FSFs of the trichotomy are shared by two or three Superkingdoms. Parsimony analysis of FSF distributions suggests that the FSF complement of the last common ancestor to the trichotomy was more like that of modern eukaryotes than that of archaea and bacteria. Studies of length distributions among members of orthologous families of proteins present in all three Superkingdoms reveal that such lengths are significantly longer among eukaryotes than among bacteria and archaea. The data also reveal that proteins lengths of eukaryotes are more broadly distributed than they are within archaeal and bacterial members of the same orthologous families. Accordingly, selective pressure for a minimal size is significantly greater for orthologous protein lengths in archaea and bacteria than in eukaryotes. Alignments of orthologous proteins of archaea, bacteria and eukaryotes are characterized by greater sequence variation at their N-terminal and C-terminal domains, than in their central cores. Length variations tend to be localized in the terminal sequences; the conserved sequences of orthologous families are localized in a central core. These data are consistent with the interpretation that the genomes of the last common ancestor (LUCA) encoded a cohort of FSFs not very different from that of modern eukaryotes. Divergence of bacterial and archaeal genomes from that common ancestor may have been accompanied by more intensive reductive evolution of proteomes than that expressed in eukaryotes. Dollo's Law suggests that the evolution of novel FSFs is a very slow process, while laboratory experiments suggests that novel protein genesis from preexisting FSFs can be relatively rapid. Reassortment of FSFs to create novel proteins may have been mediated by genetic recombination before the advent of more efficient splicing mechanisms.     

Sequence and structure classification of kinases

Kinases are a ubiquitous group of enzymes that catalyze the phosphoryl transfer reaction from a phosphate donor (usually ATP) to a receptor substrate. Although all kinases catalyze essentially the same phosphoryl transfer reaction, they display remarkable diversity in their substrate specificity, structure, and the pathways in which they participate. In order to learn the relationship between structural fold and functional specificities in kinases, we have done a comprehensive survey of all available kinase sequences (>17,000) and classified them into 30 distinct families based on sequence similarities. Of these families, 19, covering nearly 98% of all sequences, fall into seven general structural folds for which three-dimensional structures are known. These fold groups include some of the most widespread protein folds, such as Rossmann fold, ferredoxin fold, ribonuclease H fold, and TIM beta/alpha-barrel. On the basis of this classification system, we examined the shared substrate binding and catalytic mechanisms as well as variations of these mechanisms in the same fold groups. Cases of convergent evolution of identical kinase activities occurring in different folds are discussed.     

SURVEY AND SUMMARY: holliday junction resolvases and related nucleases: identification of new families, phyletic distribution and evolutionary trajectories

Holliday junction resolvases (HJRs) are key enzymes of DNA recombination. A detailed computer analysis of the structural and evolutionary relationships of HJRs and related nucleases suggests that the HJR function has evolved independently from at least four distinct structural folds, namely RNase H, endonuclease, endonuclease VII-colicin E and RusA. The endonuclease fold, whose structural prototypes are the phage lambda exonuclease, the very short patch repair nuclease (Vsr) and type II restriction enzymes, is shown to encompass by far a greater diversity of nucleases than previously suspected. This fold unifies archaeal HJRs, repair nucleases such as RecB and Vsr, restriction enzymes and a variety of predicted nucleases whose specific activities remain to be determined. Within the RNase H fold a new family of predicted HJRs, which is nearly ubiquitous in bacteria, was discovered, in addition to the previously characterized RuvC family. The proteins of this family, typified by Escherichia coli YqgF, are likely to function as an alternative to RuvC in most bacteria, but could be the principal HJRs in low-GC Gram-positive bacteria and AQUIFEX: Endonuclease VII of phage T4 is shown to serve as a structural template for many nucleases, including MCR:A and other type II restriction enzymes. Together with colicin E7, endonuclease VII defines a distinct metal-dependent nuclease fold. As a result of this analysis, the principal HJRs are now known or confidently predicted for all bacteria and archaea whose genomes have been completely sequenced, with many species encoding multiple potential HJRs. Horizontal gene transfer, lineage-specific gene loss and gene family expansion, and non-orthologous gene displacement seem to have been major forces in the evolution of HJRs and related nucleases. A remarkable case of displacement is seen in the Lyme disease spirochete Borrelia burgdorferi, which does not possess any of the typical HJRs, but instead encodes, in its chromosome and each of the linear plasmids, members of the lambda exonuclease family predicted to function as HJRs. The diversity of HJRs and related nucleases in bacteria and archaea contrasts with their near absence in eukaryotes. The few detected eukaryotic representatives of the endonuclease fold and the RNase H fold have probably been acquired from bacteria via horizontal gene transfer. The identity of the principal HJR(s) involved in recombination in eukaryotes remains uncertain; this function could be performed by topoisomerase IB or by a novel, so far undetected, class of enzymes. Likely HJRs and related nucleases were identified in the genomes of numerous bacterial and eukaryotic DNA viruses. Gene flow between viral and cellular genomes has probably played a major role in the evolution of this class of enzymes. This analysis resulted in the prediction of numerous previously unnoticed nucleases, some of which are likely to be new restriction enzymes.     

A study of the evolution of inverted-topology repeats from LeuT-fold transporters using AlignMe

X-ray crystal structures have revealed that numerous secondary transporter proteins originally categorized into different sequence families share similar structures, namely, the LeuT fold. The core of this fold consists of two units of five transmembrane helices, whose conformations have been proposed to exchange to form the two alternate states required for transport. That these two units are related implies that LeuT-like transporters evolved from gene-duplication and fusion events. Thus, the origins of this structural repeat may be relevant to the evolution of transport function. However, the lack of significant sequence similarity requires sensitive sequence search methods for analyzing their evolution. To this end, we developed a software application called AlignMe, which can use various types of input information, such as residue hydrophobicity, to perform pairwise alignments of sequences and/or of hydropathy profiles of (membrane) proteins. We used AlignMe to analyze the evolutionary relationships between repeats of the LeuT fold. In addition, we identified proteins from the so-called DedA family that potentially share a common ancestor with these repeats. DedA domains have been implicated in, e.g., selenite uptake; they are found widely distributed across all kingdoms of life; two or more DedA domains are typically found per genome, and some may adopt dual topologies. These results suggest that DedA proteins existed in ancient organisms and may function as dimers, as required for a would-be ancestor of the LeuT fold. In conclusion, we provide novel insights into the evolution of this important structural motif and thus potentially into the alternating-access mechanism of transport itself.     

Evolutionary Ancestry of Eukaryotic Protein Kinases and Choline Kinases

The reversible phosphorylation of proteins catalyzed by protein kinases in eukaryotes supports an important role for eukaryotic protein kinases (ePKs) in the emergence of nucleated cells in the third superkingdom of life. Choline kinases (ChKs) could also be critical in the early evolution of eukaryotes, because of their function in the biosynthesis of phosphatidylcholine, which is unique to eukaryotic membranes. However, the genomic origins of ePKs and ChKs are unclear. The high degeneracy of protein sequences and broad expansion of ePK families have made this fundamental question difficult to answer. In this study, we identified two class-I aminoacyl-tRNA synthetases with high similarities to consensus amino acid sequences of human protein-serine/threonine kinases. Comparisons of primary and tertiary structures supported that ePKs and ChKs evolved from a common ancestor related to glutaminyl aminoacyl-tRNA synthetases, which may have been one of the key factors in the successful of emergence of ancient eukaryotic cells from bacterial colonies.     

The Structure of the NTPase That Powers DNA Packaging into Sulfolobus Turreted Icosahedral Virus 2

Biochemical reactions powered by ATP hydrolysis are fundamental for the movement of molecules and cellular structures. One such reaction is the encapsidation of the double-stranded DNA (dsDNA) genome of an icosahedrally symmetric virus into a preformed procapsid with the help of a genome-translocating NTPase. Such NTPases have been characterized in detail from both RNA and tailed DNA viruses. We present four crystal structures and the biochemical activity of a thermophilic NTPase, B204, from the nontailed, membrane-containing, hyperthermoacidophilic archaeal dsDNA virus Sulfolobus turreted icosahedral virus 2. These are the first structures of a genome-packaging NTPase from a nontailed, dsDNA virus with an archaeal host. The four structures highlight the catalytic cycle of B204, pinpointing the molecular movement between substrate-bound (open) and empty (closed) active sites. The protein is shown to bind both single-stranded and double-stranded nucleic acids and to have an optimum activity at 80°C and pH 4.5. The overall fold of B204 places it in the FtsK-HerA superfamily of P-loop ATPases, whose cellular and viral members have been suggested to share a DNA-translocating mechanism.     

Human OLA1 defines an ATPase subfamily in the Obg family of GTP-binding proteins

Purine nucleotide-binding proteins build the large family of P-loop GTPases and related ATPases, which perform essential functions in all kingdoms of life. The Obg family comprises a group of ancient GTPases belonging to the TRAFAC (for translation factors) class and can be subdivided into several distinct protein subfamilies. The founding member of one of these subfamilies is the bacterial P-loop NTPase YchF, which had so far been assumed to act as GTPase. We have biochemically characterized the human homologue of YchF and found that it binds and hydrolyzes ATP more efficiently than GTP. For this reason, we have termed the protein hOLA1, for human Obg-like ATPase 1. Further biochemical characterization of YchF proteins from different species revealed that ATPase activity is a general but previously missed feature of the YchF subfamily of Obg-like GTPases. To explain ATP specificity of hOLA1, we have solved the x-ray structure of hOLA1 bound to the nonhydrolyzable ATP analogue AMPPCP. Our structural data help to explain the altered nucleotide specificity of YchF homologues and identify the Ola1/YchF subfamily of the Obg-related NTPases as an exceptional example of a single protein subfamily, which has evolved altered nucleotide specificity within a distinct protein family of GTPases.     

Structural and functional diversity of the microbial kinome

The eukaryotic protein kinase (ePK) domain mediates the majority of signaling and coordination of complex events in eukaryotes. By contrast, most bacterial signaling is thought to occur through structurally unrelated histidine kinases, though some ePK-like kinases (ELKs) and small molecule kinases are known in bacteria. Our analysis of the Global Ocean Sampling (GOS) dataset reveals that ELKs are as prevalent as histidine kinases and may play an equally important role in prokaryotic behavior. By combining GOS and public databases, we show that the ePK is just one subset of a diverse superfamily of enzymes built on a common protein kinase-like (PKL) fold. We explored this huge phylogenetic and functional space to cast light on the ancient evolution of this superfamily, its mechanistic core, and the structural basis for its observed diversity. We cataloged 27,677 ePKs and 18,699 ELKs, and classified them into 20 highly distinct families whose known members suggest regulatory functions. GOS data more than tripled the count of ELK sequences and enabled the discovery of novel families and classification and analysis of all ELKs. Comparison between and within families revealed ten key residues that are highly conserved across families. However, all but one of the ten residues has been eliminated in one family or another, indicating great functional plasticity. We show that loss of a catalytic lysine in two families is compensated by distinct mechanisms both involving other key motifs. This diverse superfamily serves as a model for further structural and functional analysis of enzyme evolution.