Evolution of the Cytoplasmic and Mitochondrial Phosphagen Kinases Unique to Annelid Groups

Creatine kinase (CK) is a member of a group of phosphoryl transfer enzymes called phosphagen kinases that play a key role in cellular energy transactions in animals. Three CK isoform gene families are known—cytoplasmic CK (CK), flagellar CK (fCK), and mitochondrial CK (MiCK). Each of the isoforms has a unique gene structure (intron/exon organization). A broad array of other phosphagen kinases is present in animals. Some of these enzymes are found only in annelids and closely related groups including glyocyamine kinase (GK), lombricine kinase (LK), taurocyamine kinase (TK), and a unique arginine kinase (AK) restricted to annelids. Phylogenetic analyses of these annelid phosphagen kinases indicate that they appear to have evolved from a CK-like ancestor. To gain a greater understanding of the relationship of the CK isoforms to the annelid enzymes, we have determined the intron/exon organization of the genes for the following phosphagen kinases: Eisenia LK, Sabellastarte AK, and Arenicola mitochondrial TK (MiTK). Analysis of genomic database for the polychaete Capitella sp. yielded two putative LK genes [cytoplasmic LK and mitochondrial LK (MiLK)]. The intron/exon organization of these genes was compared with available data for cytoplasmic and mitochondrial CKs, and an annelid GK. Surprisingly, these annelid genes, irrespective of whether they are cytoplasmic (LK, AK, and GK) or mitochondrial (MiTK and MiLK), had the same 8-intron/9-exon organization and were strikingly similar to MiCK genes sharing seven of eight splice junctions. These results support the view that the MiCK gene is basal and ancestral to the phosphagen kinases unique to annelids.     

Evolution of the diverse array of phosphagen systems present in annelids

Annelids as a group express a variety of phosphagen kinases including creatine kinase (CK), glyocyamine kinase (GK), lombricine kinase (LK), taurocyamine kinase (TK) and a unique arginine kinase (AK) restricted to annelids. In prior work, we have determined and compared the intron/exon organization of the annelid genes for cytoplasmic GK, LK, AK, and mitochondrial TK and LK (MiTK and MiLK, respectively), and found that these annelid genes, irrespective of cytoplasmic or mitochondrial, have the same 8-intron/9-exon organization strikingly similar to mitochondrial CK (MiCK) genes. These results support the view that the MiCK gene is basal and ancestral to the phosphagen kinases unique to annelids. To gain a greater understanding of the evolutionary processes leading to the diversity of annelid phosphagen kinases, we determined for the first time the intron/exon organization of a cytoplasmic CK gene from a polychaete as well as that of another polychaete MiCK gene. These gene structures, coupled with a phylogenetic analyses of annelid enzymes and assessment of the fidelity of substrate specificity of some these phosphagen kinases, provide insight into the pattern of radiation of the annelid enzymes. Annelid phosphagen kinases appeared to have diverged in the following order (earliest first): (1) cytoplasmic AK, LK and TK, (2) GK, and (3) mitochondrial MiLK and MiTK. Interestingly, phylogenetic analyses showed that the above phosphagen kinases appear to be basal to all CK isoforms (mitochondrial, cytoplasmic and flagellar CKs). This somewhat paradoxical placement of CKs most likely reflects a higher rate of evolution and radiation of the annelid-specific LK, TK and GK genes than the CK isoform genes.     

The Early Evolution of the Phosphagen Kinases—Insights from Choanoflagellate and Poriferan Arginine Kinases

Arginine kinase (AK) is a member of a large family of phosphoryl transfer enzymes called phosphagen (guanidino) kinases. AKs are present in certain protozoans, sponges, cnidarians, and both lophotrochozoan and ecdysozoan protostomes. Another phosphagen kinase, creatine kinase (CK), is found in sponges, cnidarians, and both deuterostome and protostome groups but does not appear to be present in protozoans. To probe the early evolution of phosphagen kinases, we have amplified the cDNAs for AKs from three choanoflagellates and from the hexactinellid sponge Aphrocallistes beatrix and the demosponges Suberites fuscus and Microciona prolifera. Phylogenetic analysis using maximum likelihood of these choanoflagellate and sponge AKs with other AK sequences revealed that the AK from the choanoflagellate Monosiga brevicollis clusters with the AK from the glass sponge Aphrocallistes and is part of a larger cluster containing AKs from the demosponges Suberites and Microciona as well as basal and protostome invertebrates. In contrast, AKs from Codonosiga gracilis and Monosiga ovata form a distinct cluster apart from all other AK sequences. tBLASTn searches of the recently released M. brevicollis genome database showed that this species has three unique AK genes—one virtually identical to the M. brevicollis cDNA and the other two showing great similarity to C. gracilis and M. ovata AKs. Three distinct AK genes are likely present in choanoflagellates. Two of these AKs display extensive similarity to both CKs and an AK from sponges. Previous work has shown CK evolved from an AK-like ancestor prior to the divergence of sponges. The present results provide evidence suggesting that the initial gene duplication event(s) leading to the CK lineage may have occurred before the divergence of the choanoflagellate and animal lineages.     

Characterization of creatine kinase isoforms in herring (Clupea harengus) skeletal muscle

It is known that mitochondrial creatine kinase (MtCK) in mammals is always expressed in conjunction with one of the cytosolic forms of creatine kinase (CK), either muscle-type (MM-CK) or brain-type (BB-CK) in tissues of high, sudden energy demand. The two creatine kinase (CK) isoforms were detected in herring (Clupea harengus) skeletal muscle: cytosolic CK and mitochondrial CK (MtCK) that displayed the different electrophoretic mobility. These isoforms differ in molecular weight and some biochemical properties. Isolation and purification procedures allowed to obtain purified enzymes with specific activity of the 206 μmol/min/mg for cytosolic CK and 240 μmol/min/mg for MtCK. Native M_rs of the cytosolic CK and MtCK determined by gel permeation chromatography were 86.000 and 345.000, respectively. The results indicate that one of isoforms found in herring skeletal muscle is a cytosolic dimer and the other one, is a mitochondrial octamer. Octamerization of MtCK is not an advanced feature and also exists in fish. These values correspond well with published values for MtCKs and cytosolic CK isoforms from higher vertebrate classes and even from lower invertebrates.     

Sponge paraphyly and the origin of Metazoa

In order to allow critical evaluation of the interrelationships between the three sponge classes, and to resolve the question of mono‐ or paraphyly of sponges (Porifera), we used the polymerase chain reaction (PCR) to amplify almost the entire nucleic acid sequence of the 18S rDNA from several hexactinellid, demosponge and calcareous sponge species. The amplification products were cloned, sequenced and then aligned with previously reported sequences from other sponges and nonsponge metazoans and variously distant outgroups, and trees were constructed using both neighbour‐joining and maximum parsimony methods. Our results suggest that sponges are paraphyletic, the Calcarea being more related to monophyletic Eumetazoa than to the siliceous sponges (Demospongiae, Hexactinellida). These results have important implications for our understanding of metazoan origins, because they suggest that the common ancestor of Metazoa was a sponge. They also have consequences for basal metazoan classification, implying that the phylum Porifera should be abandoned. Our results support the upgrading of the calcareous sponge class to the phylum level.     

Expression of the mitochondrial creatine kinase genes

Mitochondrial Creatine Kinase (MtCK) is responsible for the transfer of high energy phosphate from mitochondria to the cytosolic carrier, creatine, and exists in mammals as two isoenzymes encoded by separate genes. In rats and humans, sarcomere-specific MtCK (sMtCK) is expressed only in skeletal and heart muscle, and has 87% nucleotide identity across the 1257 bp coding region. The ubiquitous isoenzyme of MtCK (uMtCK) is expressed in many tissues with highest levels in brain, gut, and kidney, and has 92% nucleotide identity between the 1254 bp coding regions of rat and human. Both genes are highly regulated developmentally in a tissue-specific manner. There is virtually no expression of sMtCK mRNA prior to birth. Unlike cytosolic muscle CK (MCK) and brain CK (BCK), there is no developmental isoenzyme switch between the MtCKs. Cell culture models representing the tissue-specific expression of either sMtCK or uMtCK are available, but there are no adequate developmental models to examine their regulation. Several animal models are available to examine the coordinate regulation of the CK gene family and include 1) Cardiac Stress by coarctation (sMtCK, BCK, and MCK), 2) Uterus and placenta during pregnancy (uMtCK and BCK), and 3) Diabetes and mitochondrial myopathy (sMtCK, BCK, and MCK). We report the details of these findings, and discuss the coordinate regulation of the genes necessary for high-energy transduction.     

Variable intron/exon structure in the oligochaete lombricine kinase gene

Lombricine kinase is an annelid enzyme that belongs to the phosphagen kinase family of which creatine kinase and arginine kinase are the typical representatives. The enzymes play important roles in the cellular energy metabolism of animals. Biochemical, physiological and molecular information with respect to lombricine kinase is limited compared to other phosphagen kinases. This study presents data on the cDNA sequences of lombricine kinase from two smaller oligochaetes, Enchytraeus sp. and Stylaria sp. The deduced amino acid sequences are analyzed and compared with other selected phosphagen kinases. The intron/exon structure of the lombricine kinase gene was determined for these two species as well as two additional oligochaetes, Lumbriculus variegatus and Tubifex tubifex, and compared with available data for annelid phosphagen kinases. The data indicate the existence of a variable organization of the proposed 8-intron/9-exon gene structure. The results provide further insights in the evolution and position of these enzymes within the phosphagen kinase family.     

Structural characterization and tissue-specific expression of the mRNAs encoding isoenzymes from two rat mitochondrial creatine kinase genes.

Creatine kinase (CK; EC 2.7.3.2) isoenzymes play prominent roles in energy transduction. Mitochondrial CK (MtCK) reversibly catalyzes the transfer of high energy phosphate to creatine and exists, in the human, as two isoenzymes encoded by separate genes. We report here the cDNA sequences of the two isoenzymes of MtCK in the rat. Rat sarcomeric MtCK has 87% nucleotide identity in the 1257 bp coding region and 82% in the 154 bp 3' untranslated region as compared with human sarcomeric MtCK. Rat ubiquitous MtCK has 92% nucleotide identity over the 1254 bp coding region with human ubiquitous MtCK and 81% identity of the 148 by 3' untranslated region. Nucleotide identity between the rat sarcomeric and ubiquitous MtCK coding regions is 70%, with no conservation of their 3' untranslated regions. Thus, MtCK sequence is conserved in a tissue-specific, rather than species-specific, manner. Conservation of the 3' untranslated regions is highly unusual and suggests a regulatory function for this region. The NH2-terminal transit peptide sequences share 82% amino acid homology between rat and human sarcomeric MtCKs and 92% homology between rat and human ubiquitous MtCKs, but have only 41% homology to each other. This tissue-specific conservation of the transit peptides suggests receptor specificity in mitochondrial uptake. Rat sarcomeric MtCK mRNA is expressed only in skeletal muscle and heart, but rat ubiquitous MtCK mRNA is expressed in many tissues, with highest levels in brain, gut and kidney. Ubiquitous MtCK mRNA levels are dramatically regulated in uterus and placenta during pregnancy. Coexpression of sarcomeric and ubiquitous MtCK with their cytosolic counterparts, MCK and BCK, respectively, supports the creatine phosphate shuttle hypothesis and suggests that expression of these genes is coordinately regulated.     

Organization of the gene for an invertebrate mitochondrial creatine kinase: comparisons with genes of higher forms and correlation of exon boundaries with functional domains

Two major gene duplication events are thought to have taken place in the evolution of creatine kinases (CK) in the vertebrates - (1) the formation of distinct mitochondrial (MiCK) and cytoplasmic forms from the primordial gene and (2) subsequent formation of the sarcomeric (sar-) and ubiquitous (ubi-) isoforms of octameric MiCK and muscle (M) and brain (B) isoforms of dimeric, cytoplasmic CK. The genes of these two CK clades reflect a distant divergence as sar- and ubiMiCK genes consistently have nine protein-coding exons while M- and B-CK genes have seven protein-coding exons; these genes share only one common exon. CKs are also widely distributed in the invertebrates and it has recently been shown that MiCKs evolved well before the divergence of the major metazoan groups. In the present communication, we report the structure and topology of the gene for MiCK from the protostome marine worm Chaetopterus variopedatus. The protein-coding region of the gene for this primitive MiCK spans over 10 kb and consists of eight exons, the last five (E4-E8) have identical boundaries to the corresponding exons of sar- and ubiMiCK genes. Exon-3 of the C. variopedatus MiCK gene consists of the corresponding E3 and E4 of the vertebrate MiCKs with no intervening intron. E1 is longer and E2 is shorter in the polychaete MiCK gene than the counterpart sarcomeric and ubiquitous genes. The insertion of the intron in C. variopedatus E3 creating the two exons as well as the rearrangement of the intron between E1 and E2 must have occurred prior to or coincident with the duplication event creating the two vertebrate mitochondrial isoforms. Sarcomeric and ubiMiCKs display substantial differences from their invertebrate MiCK counterparts in properties relating to octamer stability and membrane binding. The evolutionary changes in gene topology may be a component of this functional progression.     

Isolation and characterization of the gene and cDNA encoding human mitochondrial creatine kinase

Creatine kinase (CK; EC 2.7.3.2) isoenzymes play prominent roles in energy metabolism. Nuclear genes encode three known CK subunits: cytoplasmic muscle (MCK), cytoplasmic brain (BCK), and mitochondrial (MtCK). We have isolated the gene and cDNA encoding human placental MtCK. By using a dog heart MCK cDNA-derived probe, the 7.0-kb EcoRI fragment from one cross-hybridizing genomic clone was isolated and its complete nucleotide sequence determined. A region of this clone encoded predicted amino acid sequence identical to residues 15-26 of the human heart MtCK NH2-terminal protein sequence. The human placental MtCK cDNA was isolated by hybridization to a genomic fragment encoding this region. The human placental MtCK gene contains 9 exons encoding 416 amino acids, including a 38-amino acid transit peptide, presumably essential for mitochondrial import. Residues 1-14 of human placental MtCK cDNA-derived NH2-terminal sequence differ from the human heart MtCK protein sequence, suggesting that tissue-specific MtCK mRNAs are derived from multiple MtCK genes. RNA blot analysis demonstrated abundant MtCK mRNA in adult human ventricle and skeletal muscle, low amounts in placenta and small intestine, and a dramatic increase during in vitro differentiation induced by serum-deprivation in the non-fusing mouse smooth muscle cell line, BC3H1. These findings demonstrate coordinate regulation of MtCK and cytosolic CK gene expression and support the phosphocreatine shuttle hypothesis.     

Phylogenetic Position of the Hexactinellida Within the Phylum Porifera Based on the Amino Acid Sequence of the Protein Kinase C from Rhabdocalyptus dawsoni

Recent analyses of genes encoding proteins typical for multicellularity, especially adhesion molecules and receptors, favor the conclusion that all metazoan phyla, including the phylum Porifera (sponges), are of monophyletic origin. However, none of these data includes cDNA encoding a protein from the sponge class Hexactinellida. We have now isolated and characterized the cDNA encoding a protein kinase C, belonging to the C subfamily (cPKC), from the hexactinellid sponge Rhabdocalyptus dawsoni. The two conserved regions, the regulatory part with the pseudosubstrate site, the two zinc fingers, and the C2 domain, as well as the catalytic domain were used for phylogenetic analyses. Sequence alignment and construction of a phylogenetic tree from the catalytic domains revealed that the yeast Saccharomyces cerevisiae and the protozoan Trypanosoma brucei are at the base of the tree, while the hexactinellid R. dawsoni branches off first among the metazoan sequences; the other two classes of the Porifera, the Calcarea (the sequence from Sycon raphanus was used) and the Demospongiae (sequences from Geodia cydonium and Suberites domuncula were used), branch off later. The statistically robust tree also shows that the two cPKC sequences from the higher invertebrates Drosophila melanogaster and Lytechinus pictus are most closely related to the calcareous sponge. This finding was also confirmed by comparing the regulatory part of the kinase gene. We suggest, that (i) within the phylum Porifera, the class Hexactinellida diverged first from a common ancestor to the Calcarea and the Demospongiae, which both appeared later, and (ii) the higher invertebrates are more closely related to the calcareous sponges. Received: 6 August 1997 / Accepted: 24 October 1997     

The role of an absolutely conserved tryptophan residue in octamer formation and stability in mitochondrial creatine kinases

In most organisms, mitochondrial creatine kinase (MtCK) is present as dimers and octamers with the latter predominating under physiological conditions. An absolutely conserved tryptophan residue (Trp-264 in chicken sarcomeric MtCK) appears to play a key role in octamer stability. Recently, it has been shown that the sponge Tethya aurantia, a member of the most ancient group of living multi-cellular animals, expresses an obligate, dimeric MtCK that lacks this absolutely conserved tryptophan residue, instead possessing a tyrosine in this position. In the present study we confirm that the absolutely conserved tryptophan residue is lacking in other sponge MtCKs where it is instead substituted by histidine or asparagine. Site directed mutations of the Trp-264 in expression constructs of chicken sarcomeric MtCK and the octameric MtCK from the marine worm Chaetopterus destabilized the octameric quaternary structure producing only dimers. A Tyr-->Trp mutation in an expression construct of the Tethya MtCK construct failed to produce octamerization; Tyr-->His and Tyr-->Asn mutations also yielded dimers. These results, in conjunction with analysis of homology models of Chaetopterus and Tethya MtCKs, strongly support the view that while the absolutely conserved tryptophan residue is important in octamer stability, octamer formation involves a complex suite of interactions between a variety of residues.     

Evolution of Metazoan Cell Junction Proteins: The Scaffold Protein MAGI and the Transmembrane Receptor Tetraspanin in the Demosponge Suberites domuncula

Until recently the positioning of the sponges (phylum Porifera) within the metazoan systematics was hampered by the lack of molecular evidence for the existence of junctional structures in the surface cell layers. In this study two genes related to the tight junctions are characterized from the demosponge Suberites domuncula: tetraspanin (SDTM4SF), a cell surface receptor, and MAGI (SDMAGI), a MAGUK (membrane-associated guanylate kinase homologue) protein. Especially the MAGI protein is known in other metazoan animal phyla to exist exclusively in tight junctions. The characteristic domains of MAGI proteins (six PDZ domains, two WW domains, and a truncated guanylate kinase motif) are conserved in the sponge protein. The functional analysis of SDMAGI done by in situ hybridization shows its expression in the surface epithelial layers (exopinacoderm and endopinacoderm). Northern blot studies reveal that expression of SDMAGI and SDTM4SF increases after formation of the pinacoderm layer in the animals as well as in primmorphs. These results support earlier notions that sponges contain junctional structures. We conclude that sponges contain epithelia whose cells are organized by cell junctions.The sequence from Suberites domuncula reported here, the protein membrane-associated guanylate kinase with an inverted arrangement (MAGI), is deposited in the EMBL/GenBank database under accession number AJ580406.     

Evolution and Divergence of the Genes for Cytoplasmic,Mitochondrial, and Flagellar Creatine Kinases

Creatine kinase (CK) plays a central role in energy homeostasis in cells that display high and variable rates of energy turnover. A number of CK genes exist, each being targeted to particular intracellular compartments. In the vertebrates, two genes code for proteins which form homo- and heterodimers targeted to the cytoplasm, while two additional genes code for primarily octameric proteins targeted to the mitochondrial intermembrane space. Yet another gene is present in certain groups which codes for three fused, complete CK domains and is typically targeted to the flagellar membrane of primitive-type spermatozoa. CK is widely distributed in protochordates and both protostome and deuterostome invertebrate groups. The evolutionary relationships of these CK genes have not been fully elucidated. The present communication reports new cDNA-derived deduced amino acid sequences for four cytoplasmic and three mitochondrial CKs and one flagellar CK from lophotrochozoan, protostome invertebrates as well as a new cytoplasmic CK sequence from a protochordate tunicate. These new sequences, coupled with available sequences in the databases and sequences extracted from genome sequencing projects, provide revealing insights into the evolution and divergence of CK genes. Phylogenetic analyses showed that single cytoplasmic, mitochondrial, and flagellar CK genes were present prior to the divergence of the protostomes and deuterostomes. The flagellar CK gene may have evolved within the cytoplasmic gene clade, although the evidence is somewhat equivocal. The two cytoplasmic genes in the vertebrates, and most likely the two mitochondrial genes, evolved after the divergence of the craniates from the protochordates. Comparison of the structure of the genes for selected cytoplasmic, mitochondrial, and flagellar CKs revealed two identical intron boundaries, further reinforcing the notion of a common evolutionary origin, but also showed patterns of changes in structure consistent with each gene type. These studies show that the cytoplasmic, mitochondrial, and flagellar CK genes are rather ancient and that there has been a systematic pattern of duplication and divergence consistent with changing nature of energy demands and physicochemical environment in the cells where they are expressed.[Reviewing Editor: Martin Kreitman]     

Over-expression, purification and characterization of the oligomerization dynamics of an invertebrate mitochondrial creatine kinase

Mitochondrial creatine kinase (MtCK) plays a central role in energy homeostasis within cells that display high and variable rates of ATP turnover. Vertebrate MtCKs exist primarily as octamers but readily dissociate into constituent dimers under a variety of circumstances. MtCK is an ancient protein that is also found in invertebrates including sponges, the most primitive of all multi-cellular animals. We have cloned, expressed, and purified one of these invertebrate MtCKs from a marine polychaete worm, Chaetopterus variopedatus (CVMtCK). Size exclusion chromatography and dynamic light scattering (DLS) were used to characterize oligomeric state in comparison with that of octameric chicken sarcomeric isoform (SarMtCK). At protein concentrations >1 mg/ml, CVMtCK was predominantly octameric (>90%). When diluted to 0.1 mg/ml, CVMtCK dissociated into dimers much more rapidly than SarMtCK when observed under identical conditions. The rate of dissociation for both MtCKs increased as temperature rose from 2 to 28 degrees C, and in CVMtCK, fell at higher incubation temperatures. The fraction of octameric CVMtCK at equilibrium increased with temperature and then fell. Temperature transition studies showed that octamers and dimers were rapidly interconvertible on a similar time scale. Importantly, when CVMtCK was converted to the transition state analog complex (TSAC), both size exclusion chromatography and DLS showed that there was minimal dissociation of octamers into dimers while SarMtCK octamers were highly unstable as the TSAC. These results clearly show distinct differences in octamer stability between CVMtCK and SarMtCK, which could impact function under physiological circumstances. Furthermore, the large yield of recombinant protein and high stability of CVMtCK in the TSAC suggest that this protein might be a good target for crystallization efforts.     

Polarity factor 'Frizzled' in the demosponge Suberites domuncula: identification, expression and localization of the receptor in the epithelium/pinacoderm(1)

Until recently, it was assumed that polarity and axis formation have evolved only in metazoan phyla higher than Cnidaria. One key molecule involved in the signal transduction causing tissue polarity is Frizzled, a seven-transmembrane receptor that is activated by the Wnt family of secreted proteins. We report the isolation and characterization of a Frizzled gene from the demosponge Suberites domuncula (Sd-Fz). The deduced polypeptide comprises all characteristic domains known from Frizzled receptors of higher metazoans. In situ hybridization studies show that Sd-Fz is expressed in cells close to the surface of the sponges and in the pinacocytes of some canals. Northern blot analysis demonstrates its upregulation during the formation of three-dimensional sponge cell aggregates in culture. These data provide for the first time experimental evidence that already in the lowest metazoan phylum (Porifera) genes are present which are very likely involved in tissue polarity.     

The role of Arg-96 in Danio rerio creatine kinase in substrate recognition and active center configuration

In creatine kinases (CKs), the amino acid residue-96 is a strictly conserved arginine. This residue is not directly associated with substrate binding, but it is located close to the binding site of the substrate creatine. On the other hand, the residue-96 is known to be involved in expression in the substrate specificity of various other phosphagen (guanidino) kinases, since each enzyme has a specific residue at this position: arginine kinase (Tyr), glycocyamine kinase (Ile), taurocyamine kinase (His) and lombricine kinase (Lys). To gain a greater understanding of the role of residue-96 in CKs, we replaced this residue in zebra fish Danio rerio cytoplasmic CK with other 19 amino acids, and expressed these constructs in Escherichia coli. All the twenty recombinant enzymes, including the wild-type, were obtained as soluble form, and their activities were determined in the forward direction. Compared with the activity of wild-type, the R96K mutant showed significant activity (8.3% to the wild-type), but 10 mutants (R96Y, A, S, E, H, T, F, C, V and N) showed a weak activity (0.056–1.0%). In the remaining mutants (R96Q, G, M, P, L, W, D and I), the activity was less than 0.05%. Our mutagenesis studies indicated that Arg-96 in Danio CK can be substituted for partially by Lys, but other replacements caused remarkable loss of activity. From careful inspection of the crystal structures (transition state analog complex (TSAC) and open state) of Torpedo cytoplasmic CK, we found that the side chain of R96 forms hydrogen bonds with A339 and D340 only in the TSAC structure. Based on the assumption that CKs consist of four dynamic domains (domains 1–3, and fixed domain), the above hydrogen bonds act to link putative domains 1 and 3 in TSAC structure. We suggest that residue-96 in CK and equivalent residues in other phosphagen kinases, which are structurally similar, have dual roles: (1) one involves in distinguishing guanidino substrates, and (2) the other plays a key role in organizing the hydrogen-bond network around residue-96 which offers an appropriate active center for the high catalytic turnover. The mode of development of the network appears to be unique each phosphagen kinase, reflecting evolution of each enzyme.     

Molecular Evolution of TEPP Protein Genes in Metazoans

TEPP is a gene expressed in human reproductive organs such as testis, prostate, and placenta. Here, identification and molecular evolutionary analysis of TEPP proteins in various metazoan animals including deuterostomes (chordates, hemichordates, and echinoderms), lophotrochozoans (mollusks and annelids), and cnidarians (sea anemone and coral) are reported. A multiple sequence alignment revealed two highly conserved regions in TEPP proteins that had no similarity to any other known domains or proteins. Genomic sequence analysis showed frequent shifting of the splice sites of intron 1 in mammalian TEPP genes. A comparison of the intron positions in the coding region showed that the exon/intron structure of the TEPP gene was established in an early metazoan ancestor and that independent loss of a single intron occurred in echinoderms and in vertebrates. The urochordate tunicate TEPP genes are intronless, possibly due to replacement of the original gene by a retrogene. No homolog was detected in birds, insects, nematodes, and teleost fishes despite the extensive sequence data of these species, implying that the TEPP gene might be lost in these lineages.