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

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.

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 arginine kinase gene family

Arginine kinase (AK), catalyzing the reversible transfer of phosphate from MgATP to arginine yielding phosphoarginine and MgADP, is widely distributed throughout the invertebrates and is also present in certain protozoa. Typically, these proteins are found as monomers targeted to the cytoplasm, but true dimeric and contiguous dimeric AKs as well as mitochondrial AK activities have been observed. In the present study, we have obtained the sequences of the genes for AKs from two distantly related molluscs-the cephalopod Nautilus pompilius and the bivalve Crassostrea gigas. These new data were combined with available gene structure data (exon/intron organization) extracted from EST and genome sequencing project databases. These data, comprised of 23 sequences and gene structures from Protozoa, Cnidaria, Platyhelminthes, Mollusca, Arthropoda and Nematoda, provide great insight into the evolution and divergence of the AK family. Sequence and phylogenetic analyses clearly show that the AKs are homologous having arisen from some common ancestor. However, AK gene organization is highly divergent and variable. Molluscan AK genes typically have a highly conserved six-exon/five-intron organization, a structure that is very similar to that of the platyhelminth Schistosoma mansoni Arthropod and nematode AK genes have fewer introns, while the cnidarian and protozoan genes each display unique exon/intron organization when compared to the other AK genes. The non-conservative nature of the AK genes is in sharp contrast to the relatively high degree of conservation of intron positions seen in a homologous enzyme creatine kinase (CK). The present results also show that gene duplication and subsequent fusion events forming unusual two-domain AKs occurred independently at least four times as these contiguous dimers are present in Protozoa, Cnidaria, Platyhelminthes and Mollusca. Detailed analyses of the amino acid sequences indicate that two AKs (one each from Drosophila and Caenorhabditis) have what appear to be N-terminal mitochondrial targeting sequences, providing the first evidence for true mitochondrial AK genes. The AK gene family is ancient and the lineage has undergone considerable divergence as well as multiple duplication and fusion events.

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.

Alternative splicing produces transcripts coding for alpha and beta chains of a hetero-dimeric phosphagen kinase

Glycocyamine kinase (GK) catalyzes the reversible phosphorylation of glycocyamine (guanidinoacetate), a reaction central to cellular energy homeostasis in certain animals. GK is a member of the phosphagen kinase enzyme family and appears to have evolved from creatine kinase (CK) early in the evolution of multi-cellular animals. Prior work has shown that GK from the polychaete Neanthes (Nereis) diversicolor exits as a hetero-dimer in vivo and that the two polypeptide chains (termed alpha and beta) are coded for by unique transcripts. In the present study, we demonstrate that the GK from a congener Nereis virens is also hetero-dimeric and is coded for by alpha and beta transcripts, which are virtually identical to the corresponding forms in N. diversicolor. The GK gene from N. diversicolor was amplified by PCR. Sequencing of the PCR products showed that the alpha and beta chains are the result of alternative splicing of the GK primary mRNA transcript. These results also strongly suggest that this gene underwent an early tandem exon duplication event. Full-length cDNAs for N. virens GKalpha and GKbeta were individually ligated into expression vectors and the resulting constructs used to transform Escherichia coli expression hosts. Regardless of expression conditions, minimal GK activity was observed in both GKalpha and GKbeta constructs. Inclusion bodies for both were harvested, unfolded in urea and alpha chains, beta chains and mixtures of alpha and beta chains were refolded by sequential dialysis. Only modest amounts of GK activity were observed when alpha and beta were refolded individually. In contrast, when refolded the alpha and beta mixture yielded highly active hetero-dimers, as validated by size exclusion chromatography, electrophoresis and mass spectrometry, with a specific activity comparable to that of natural GK. The above evidence suggests that there is a preference for hetero-dimer formation in the GKs from these two polychaetes. The evolution of the alternate splicing and an additional exon in these GKs, producing alpha and beta transcripts, can be viewed as a possible compensation for a mutation(s) in the original gene, which most likely coded for a homo-dimeric protein.

C-terminal Lysines Determine Phospholipid Interaction of Sarcomeric Mitochondrial Creatine Kinase

High affinity interaction between octameric mitochondrial creatine kinase (MtCK) and the phospholipid cardiolipin in the inner mitochondrial membrane plays an important role in metabolite channeling between MtCK and inner membrane adenylate translocator, which itself is tightly bound to cardiolipin. Three C-terminal basic residues revealed as putative cardiolipin anchors in the x-ray structures of MtCK and corresponding to lysines in human sarcomeric MtCK (sMtCK) were exchanged by in vitro mutagenesis (K369A/E, K379Q/A/E, K380Q/A/E) to yield double and triple mutants. sMtCK proteins were bacterially expressed, purified to homogeneity, and verified for structural integrity by enzymatic activity, gel filtration chromatography, and CD spectroscopy. Interaction with cardiolipin and other acidic phospholipids was quantitatively analyzed by light scattering, surface plasmon resonance, and fluorescence spectroscopy. All mutant sMtCKs showed a strong decrease in vesicle cross-linking, membrane affinity, binding capacity, membrane ordering capability, and binding-induced changes in protein structure as compared with wild type. These effects did not depend on the nature of the replacing amino acid but on the number of exchanged lysines. They were moderate for Lys-379/Lys-380 double mutants but pronounced for triple mutants, with a 30-fold lower membrane affinity and an entire lack of alterations in protein structure compared with wild-type sMtCK. However, even triple mutants partially maintained an increased order of cardiolipin-containing membranes. Thus, the three C-terminal lysines determine high affinity sMtCK/cardiolipin interaction and its effects on MtCK structure, whereas low level binding and some effect on membrane fluidity depend on other structural components. These results are discussed in regard to MtCK microcompartments and evolution.

Kinetic properties and structural characteristics of an unusual two-domain arginine kinase of the clam Corbicula japonica

Arginine kinase (AK) from the clam Corbicula japonica is a unique enzyme in that it has an unusual two-domain structure with molecular mass of 80 kDa. It lacks two functionally important amino acid residues, Asp-62 and Arg-193, which are conserved in other 40 kDa AKs and are assumed to be key residues for stabilizing the substrate-bound structure. K m arg and Vmax values for the recombinant two-domain AK were determined. These values were close to those of usual 40 kDa AKs, although Corbicula AK lacks the functionally important Asp-62 and Arg-193. Domain 2 of Corbicula AK was separated from the two-domain enzyme and was expressed in Escherichia coli. Domain 2 still exhibited activity. However, kinetic parameters for domain 2 appeared to be slightly, but significantly, different from those of two-domain AK. Thus, it is likely that the formation of the contiguous dimer alters the kinetic properties of its constituent domains significantly. Comparison of K d arg and K m arg for two-domain AK and its domain 2 showed that the affinity of the enzyme for arginine is greater in the presence of substrate ATP than in its absence. Presumably this difference is correlated with the large structural differences in the enzyme in the presence or absence of substrate, namely open and closed structures. We expressed three mutants of Corbicula AK domain 2 (His-60 to Gly or Arg, Asp-197 to Gly), and determined their K m arg and Vmax values. The affinity for the substrate arginine in mutant enzymes was reduced considerably, accompanied by a decrease in Vmax. These results suggest that His-60 and Asp-197 affect the substrate binding system, and are consistent with the hypothesis that a hydrogen bond is formed between His-60 and Asp-197 in Corbicula AK as a substitute for the Asp-62 and Arg-193 bond in normal AKs.