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.     

Role of amino-acid residue 95 in substrate specificity of phosphagen kinases

The purpose of this study is to elucidate the mechanisms of guanidine substrate specificity in phosphagen kinases, including creatine kinase (CK), glycocyamine kinase (GK), lombricine kinase (LK), taurocyamine kinase (TK) and arginine kinase (AK). Among these enzymes, LK is unique in that it shows considerable enzyme activity for taurocyamine in addition to its original target substrate, lombricine. We earlier proposed several candidate amino acids associated with guanidine substrate recognition. Here, we focus on amino-acid residue 95, which is strictly conserved in phosphagen kinases: Arg in CK, Ile in GK, Lys in LK and Tyr in AK. This residue is not directly associated with substrate binding in CK and AK crystal structures, but it is located close to the binding site of the guanidine substrate. We replaced amino acid 95 Lys in LK isolated from earthworm Eisenia foetida with two amino acids, Arg or Tyr, expressed the modified enzymes in Escherichia coli as a fusion protein with maltose-binding protein, and determined the kinetic parameters. The K95R mutant enzyme showed a stronger affinity for both lombricine (Km=0.74 mM and kcat/Km=19.34 s(-1) mM(-1)) and taurocyamine (Km=2.67 and kcat/Km=2.81), compared with those of the wild-type enzyme (Km=5.33 and kcat/Km=3.37 for lombricine, and Km=15.31 and kcat/ Km=0.48for taurocyamine). Enzyme activity of the other mutant, K95Y, was dramatically altered. The affinity for taurocyamine (Km=1.93 and kcat/Km=6.41) was enhanced remarkably and that for lombricine (Km=14.2 and kcat/Km=0.72) was largely decreased, indicating that this mutant functions as a taurocyamine kinase. This mutant also had a lower but significant enzyme activity for the substrate arginine (Km=33.28 and kcat/Km=0.01). These results suggest that Eisenia LK is an inherently flexible enzyme and that substrate specificity is strongly controlled by the amino-acid residue at position 95.     

Mutagenesis of two acidic active site residues in human muscle creatine kinase: implications for the catalytic mechanism

Creatine kinase (CK) catalyzes the reversible phosphorylation of the guanidine substrate, creatine, by MgATP. Although several X-ray crystal structures of various isoforms of creatine kinase have been published, the detailed catalytic mechanism remains unresolved. A crystal structure of the CK homologue, arginine kinase (AK), complexed with the transition-state analogue (arginine-nitrate-ADP), has revealed two carboxylate amino acid residues (Glu225 and Glu314) within 2.8 A of the proposed transphosphorylation site. These two residues are the putative catalytic groups that may promote nucleophilic attack by the guanidine amino group on the gamma-phosphate of ATP. From primary sequence alignments of arginine kinases and creatine kinases, we have identified two homologous creatine kinase acidic amino acid residues (Glu232 and Asp326), and these were targeted for examination of their potential roles in the CK mechanism. Using site-directed mutagenesis, we have made several substitutions at these two positions. The results indicate that of these two residues the Glu232 is the likely catalytic residue while Asp326 likely performs a role in properly aligning substrates for catalysis.     

Enzymatically inactive adenylate kinase 4 interacts with mitochondrial ADP/ATP translocase

Adenylate kinase 4 (AK4) is a unique member with no enzymatic activity in vitro in the adenylate kinase (AK) family although it shares high sequence homology with other AKs. It remains unclear what physiological function AK4 might play or why it is enzymatically inactive. In this study, we showed increased AK4 protein levels in cultured cells exposed to hypoxia and in an animal model of the neurodegenerative disease amyotrophic lateral sclerosis. We also showed that short hairpin RNA (shRNA)-mediated knockdown of AK4 in HEK293 cells with high levels of endogenous AK4 resulted in reduced cell proliferation and increased cell death. Furthermore, we found that AK4 over-expression in the neuronal cell line SH-SY5Y with low endogenous levels of AK4 protected cells from H₂O₂ induced cell death. Proteomic studies revealed that the mitochondrial ADP/ATP translocases (ANTs) interacted with AK4 and higher amount of ANT was co-precipitated with AK4 when cells were exposed to H₂O₂ treatment. In addition, structural analysis revealed that, while AK4 retains the capability of binding nucleotides, AK4 has a glutamine residue instead of a key arginine residue in the active site well conserved in other AKs. Mutation of the glutamine residue to arginine (Q159R) restored the adenylate kinase activity with GTP as substrate. Collectively, these results indicate that the enzymatically inactive AK4 is a stress responsive protein critical to cell survival and proliferation. It is likely that the interaction with the mitochondrial inner membrane protein ANT is important for AK4 to exert the protective benefits to cells under stress.     

Two fused proteins combining Stichopus japonicus arginine kinase and rabbit muscle creatine kinase

Two fused proteins of dimeric arginine kinase (AK) from sea cucumber and dimeric creatine kinase (CK) from rabbit muscle, named AK-CK and CK-AK, were obtained through the expression of fused AK and CK genes. Both AK-CK and CK-AK had about 50% AK activity and about 2-fold K _m values for arginine of native AK, as well as about 50% CK activity and about 2-fold K _m values for creatine of native CK. This indicated that both AK and CK moieties are fully active in the two fused proteins. The structures of AK, CK, AK-CK, and CK-AK were compared by collecting data of far-UV circular dichroism, intrinsic fluorescence, 1-anilinonaphthalene-8-sulfonate binding fluorescence, and size-exclusion chromatography. The results indicated that dimeric AK and CK differed in the maximum emission wavelength, the exposure extent of hydrophobic surfaces, and molecular size, though they have a close evolutionary relationship. The structure and thermodynamic stability of AK, CK, AK-CK, and CK-AK were compared by guanidine hydrochloride (GdnHCl) titration. Dimeric AK was more dependent on the cooperation of two subunits than CK according to the analysis of residual AK or CK activity with GdnHCl concentration increase. Additionally, AK and CK had different denaturation curves induced by GdnHCl, but almost the same thermodynamic stability. The two fused proteins, AK-CK and CK-AK, had similar secondary structure, tertiary structure, molecular size, structure, and thermodynamic stability, which indicated that the expression order of AK and CK genes might have little effect on the characteristics of the fused proteins and might further verify the close relationship of dimeric AK and CK. Published in Russion in Biokhimiya, 2006, Vol. 71, No. 9, pp. 1208–1214.     

A new approach to secondary structure evaluation: Secondary structure prediction of porcine adenylate kinase and yeast guanylate kinase by CD spectroscopy of overlapping synthetic peptide segments

A new approach for evaluating the secondary structure of proteins by CD spectroscopy of overlapping peptide segments is applied to porcine adenylate kinase (AK1) and yeast guanylate kinase (GK3). One hundred seventy-six peptide segments of a length of 15 residues, overlapping by 13 residues and covering the complete sequences of AK1 and GK3, were synthesized in order to evaluate their secondary structure composition by CD spectroscopy. The peptides were prepared by solid phase multiple peptide synthesis method using the 9-fluorenylmethoxycarbonyl/tert-butyl strategy. The individual peptide secondary structures were studied with CD spectroscopy in a mixture of 30% trifluoroethanol in phosphate buffer (pH 7) and subsequently compared with x-ray data of AK1 and GK3. Peptide segments that cover α-helical regions of the AK1 or GK3 sequence mainly showed CD spectra with increasing and decreasing Cotton effects that were typical for appearing and disappearing α-helical structures. For segments with dominating β-sheet conformation, however, the application of this method is limited due to the stability and clustering of β-sheet segments in solution and due to the difficult interpretation of random-coiled superimposed β-sheet CD signals. Nevertheless, the results of this method especially for α-helical segments are very impressive. All α-helical and 71% of the β-sheet containing regions of the AK1 and GK3 could be identified. Moreover, it was shown that CD spectra of consecutive peptide content reveal the appearance and disappearance of α-helical secondary structure elements and help localizing them on the sequence string. © 1997 John Wiley & Sons, Inc. Biopoly 41: 213–231, 1997     

Cloning and expression of a lombricine kinase from an echiuroid worm: insights into structural correlates of substrate specificity

Phosphagen kinases constitute a large family of enzymes catalyzing the reversible phosphorylation of guanidino acceptor compounds. These guanidino substrates differ substantially in size and chemical properties. In spite of the appearance of X-ray crystal structures for two members of this family, creatine kinase (CK) and arginine kinase (AK), the structural correlates of substrate specificity remain to be fully elucidated. We have determined the cDNA and deduced amino acid sequences for lombricine (guanidinethylphosphoserine) kinase (LK) from the echiuroid worm Urechis caupo and expressed the cDNA in Escherichia coli. The recombinant protein was purified by affinity chromatography and showed high capacity for phosphorylation of lombricine. Phosphagen kinases consist of a small, N-terminal domain and a much larger domain connected by a linker sequence. A key event in catalysis in CK and AK, and certainly all other phosphagen kinases, is a large conformational change involving involving a rotation of the two domains and the movement of two highly conserved flexible loops (one located in the small domain; the other located in the large domain of these enzymes) which clamp down on the substrates. Multiple sequence alignments of Urechis LK with the only other LK sequence available and CK, AK and glycocyamine kinase sequences, confirm the importance of the small flexible loop located in the N-terminal domain of phosphagen kinases as one component of the structural determinants of guanidine specificity. The role of the other flexible loop in the large domain in terms of substrate specificity remains questionable.     

Hypotaurocyamine kinase evolved from a gene for arginine kinase

Hypotaurocyamine kinase (HTK) is a member of the highly conserved family of phosphagen kinases that includes creatine kinase (CK) and arginine kinase (AK). HTK is found only in sipunculid worms, and it shows activities for both the substrates hypotaurocyamine and taurocyamine. Determining how HTK evolved in sipunculids is particularly insightful because all sipunculid-allied animals have AK and only some sipunculids have HTK. We determined the cDNA sequence of HTK from the sipunculid worm Siphonosoma cumanense for the first time, cloned it in pMAL plasmid and expressed it in E. coli as a fusion protein with maltose-binding protein. The cDNAderived amino acid sequence of Siphonosoma HTK showed high amino acid identity with molluscan AKs. Nevertheless, the recombinant enzyme of Siphonosoma HTK showed no activity for the substrate arginine, but showed activity for taurocyamine. Comparison of the amino acid sequences of HTK and AK indicated that the amino acid residues necessary for the binding of the substrate arginine in AK have been completely lost in Siphonosoma HTK sequence. The phylogenetic analysis indicated that the HTK amino acid sequence was placed just outside the molluscan AK cluster, which formed a sister group with the arthropod and nematode AKs. These results suggest that Siphonosoma HTK evolved from a gene for molluscan AK. Moreover, to confirm this assertion, we determined by PCR that the gene for Siphonosoma HTK has a 5-exon/4-intron structure, which is homologous with that of the molluscan AK genes. Further, the positions of splice junctions were conserved exactly between the two genes. Thus, we conclude that Siphonosoma HTK has evolved from a primordial gene for molluscan AK.     

Expression of horseshoe crab arginine kinase in Escherichia coli and site-directed mutations of the reactive cysteine peptide

Arginine kinase (AK) from the horseshoe crab Limulus polyphemus was expressed in Escherichia coli. The bulk of expressed protein resided in insoluble inclusion bodies. However, approximately 3 mg enzyme protein/l culture was present as active soluble AK. The AK-containing expression vector construct was subjected to site-directed mutagenesis via a polymerase chain reaction-based megaprimer protocol. The AK reactive cysteine peptide was engineered so that it was identical to the corresponding peptide sequence of creatine kinase, another member of the guanidino kinase enzyme family. The resulting expressed protein had a considerably reduced specific activity but was still specific for arginine/arginine phosphate. No catalytic activity was observed with other guanidine substrates (creatine, glycocyamine, taurocyamine, lombricine). The reactive cysteine peptide, characteristic of all guanidino kinases, very likely plays a minimal role in determining guanidine specificity.     

Evolution of Phosphagen Kinase VII. Isolation of Glycocyamine Kinase from the Polychaete Neanthes diversicolor and the cDNA-Derived Amino Acid Sequences of α and β Chains

Glycocyamine kinase (GK) was isolated from the marine polychaete Neanthes diversicolor by gel filtration, DEAE-cellulose chromatography, butyl-Toyopearl hydrophobic chromatography, and chromatofocusing. The GK was eluted as a single peak on the latter three chromatographies, and the molecular mass for the native GK was estimated to be about 80 kDa. The SDS–PAGE showed that the isolated GK consists of two distinct subunits in equal proportion, α and β chains, with molecular masses of 42.2 and 43.8 kDa, respectively. The present results suggest that the Neanthes GK has a heterodimeric structure. The cDNAs for α and β chains of Neanthes GK were amplified by PCR and their cDNA-derived amino acid sequences were determined. The α and β chains are composed of 374 and 390 amino acids, and the molecular masses were calculated to be 42,392 and 43,966 Da, respectively, in good agreement with the apparent masses on SDS–PAGE. The β chain has a characteristic N-terminal extension of 15 amino acids, and all of the sequence differences between α and β chains were restricted in the N-terminal region of 50 residues. The overall sequence identity was 92%. The occurrence of heterodimeric nature in Neanthes GK is of great interest from the evolutionary point of view, because the heterodimeric structure is only known for creatine kinase MB-isozyme specific for mammalian heart muscle among phosphagen kinases.     

A new type of metal-binding site in cobalt- and zinc-containing adenylate kinases isolated from sulfate-reducers Desulfovibrio gigas and Desulfovibrio desulfuricans ATCC 27774

Adenylate kinase (AK) mediates the reversible transfer of phosphate groups between the adenylate nucleotides and contributes to the maintenance of their constant cellular level, necessary for energy metabolism and nucleic acid synthesis. The AK were purified from crude extracts of two sulfate-reducing bacteria (SRB), Desulfovibrio (D.) gigas NCIB 9332 and Desulfovibrio desulfuricans ATCC 27774, and biochemically and spectroscopically characterised in the native and fully cobalt- or zinc-substituted forms. These are the first reported adenylate kinases that bind either zinc or cobalt and are related to the subgroup of metal-containing AK found, in most cases, in Gram-positive bacteria. The electronic absorption spectrum is consistent with tetrahedral coordinated cobalt, predominantly via sulfur ligands, and is supported by EPR. The involvement of three cysteines in cobalt or zinc coordination was confirmed by chemical methods. Extended X-ray absorption fine structure (EXAFS) indicate that cobalt or zinc are bound by three cysteine residues and one histidine in the metal-binding site of the “LID” domain. The sequence ¹²⁹Cys-X₅-His-X₁₅-Cys-X₂-Cys of the AK from D. gigas is involved in metal coordination and represents a new type of binding motif that differs from other known zinc-binding sites of AK. Cobalt and zinc play a structural role in stabilizing the LID domain.     

Distinct organization of energy metabolism in HL-1 cardiac cell line and cardiomyocytes

Expression and function of creatine kinase (CK), adenylate kinase (AK) and hexokinase (HK) isoforms in relation to their roles in regulation of oxidative phosphorylation (OXPHOS) and intracellular energy transfer were assessed in beating (B) and non-beating (NB) cardiac HL-l cell lines and adult rat cardiomyocytes or myocardium. In both types of HL-1 cells, the AK2, CKB, HK1 and HK2 genes were expressed at higher levels than the CKM, CKMT2 and AK1 genes. Contrary to the saponin-permeabilized cardiomyocytes the OXPHOS was coupled to mitochondrial AK and HK but not to mitochondrial CK, and neither direct transfer of adenine nucleotides between CaMgATPases and mitochondria nor functional coupling between CK-MM and CaMgATPases was observed in permeabilized HL-1 cells. The HL-1 cells also exhibited deficient complex I of the respiratory chain. In conclusion, contrary to cardiomyocytes where mitochondria and CaMgATPases are organized into tight complexes which ensure effective energy transfer and feedback signaling between these structures via specialized pathways mediated by CK and AK isoforms and direct adenine nucleotide channeling, these complexes do not exist in HL-1 cells due to less organized energy metabolism.     

Mechanism of inhibition of arginine kinase by flavonoids consistent with thermodynamics of docking simulation

Arginine kinase plays a vital role in invertebrate homeostasis by buffering ATP concentrations. Arginine kinase might serve as a target for environmentally friendly insect-selective pesticides, because it differs notably from its counterpart in vertebrates, creatine kinase. In this study, two members of the flavonoid family, quercetin (QU) and luteolin (LU), were identified as novel noncompetitive inhibitors of locust arginine kinase. They were found to have inhibition parameters (K_i) of 11.2 and 23.9 μM, respectively. By comparing changes in the activity and intrinsic fluorescence of AK, the inhibition mechanisms of these flavonoids were found to involve binding to Trp residues in the active site. This was determined by examination of the static quenching parameter K_sv. The main binding forces between flavonoids and AK were found to be hydrophobic based on the thermodynamic parameters of changes in enthalpy (ΔH), entropy (ΔS), and free energy (ΔG) and on docking simulation results. Molecular docking analyses also suggested that flavonoids could bind to the active site of AK and were close to the Trp 221 in active site. Molecular simulation results mimic the experimental results, indicated that QU had a lower binding energy and a stronger inhibitory effect on AK than LU, suggesting that the extra hydroxyl group in QU might increase binding ability.     

Molecular characterization and kinetic properties of a novel two-domain taurocyamine kinase from the lung fluke Paragonimus westermani

Taurocyamine kinase (TK) was previously reported to be restricted to certain marine annelids; however, the present study has proven otherwise. The lung fluke Paragonimuswestermani has a contiguous two-domain TK with a mass of 80 216 Da consisting of 713 amino acid residues sharing higher sequence identity with molluscan arginine kinase (AK). Both domains of P. westermani TK have significant activity for the substrate taurocyamine and exhibited synergism during substrate binding. Since TK plays a key role in energy metabolism and is not present in mammals, inhibitors against P. westermani TK could be effective novel chemotherapeutic agents and could be utilized for the development of specific diagnostic tools for the detection of paragonimiasis.     

Expression and localization of epitope-tagged protein kinase CK2

Protein kinase CK2, formerly known as casein kinase II, is a ubiquitous protein serine/threonine kinase. The enzyme exists in tetrameric complexes composed of two catalytic (CK2α and/or CK2α′) subunits and two subunits (CK2β) that appear to have a role in modulating the activity of the catalytic subunits. With the exception of their unrelated carboxy-terminal domains, the two isozymic forms of mammalian CK2 display extensive sequence identity. Furthermore, CK2α and CK2α′ exhibit remarkable conservation between species, suggesting that they may have unique functions. In the present study, the cDNAs encoding CK2α and CK2α′ were modified by addition of the hemagglutinin tag of the influenza virus at the amino terminus of the respective proteins. The epitope-tagged proteins were transfected into Cos-7 cells and the localization of the expressed proteins determined by indirect immunofluorescence using monoclonal antibodies specific for the epitope tag. The use of transfection favors the formation of homotetrameric complexes (i.e., α₂β₂, α′₂β₂) instead of heterotetrameric complexes (i.e., αα′β₂) that are present in many cells. Epitope-tagged CK2α and CK2α′ displayed kinase activity and the ability to form complexes with CK2β. The results of these studies also indicate definitively that CK2α and CK2α′ are both localized predominantly within the nucleus. Mutation of conserved lysine residues within the ATP binding domains of CK2α and CK2α′ resulted in loss of kinase activity. However, examination of these mutants indicates that kinase activity is not essential for formation of complexes between subunits of CK2 and is not required for nuclear localization of CK2. J. Cell. Biochem. 64: 525–537. © 1997 Wiley-Liss, Inc.     

Molecular cloning and nucleotide sequence of tuna growth hormone cDNA

cDNA for mRNA of tuna growth hormone (GH) was cloned by screening a cDNA library constructed from tuna pituitary gland poly(A)⁺ RNA. The nucleotide sequence of cDNA (911 bases) revealed an open reading frame of 615 nucleotides, including a sequence (51 bases) for a possible secretory protein leader peptide. Noncoding regions were found in the nucleotide sequences up- (5′-terminal: 65 bases) and down- (3′-terminal: 231 bases) stream of the open reading frame. An amino-acid sequence deduced from the nucleotide sequence of the cDNA was identical with that determined in the purified tuna GH. Tuna GH was composed of 187 amino acids, and had a calculated molecular weight of 21275. Amino-acid sequencing showed that there was one possible N-glycosylation site at Asn (Asn-Cys-Thr). Tuna GH showed amino-acid sequence homologies with chum salmon (67%), yellow tail (90%) and with human (32%) growth hormones.