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.     

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.     

The structure of lombricine kinase: implications for phosphagen kinase conformational changes

Lombricine kinase is a member of the phosphagen kinase family and a homolog of creatine and arginine kinases, enzymes responsible for buffering cellular ATP levels. Structures of lombricine kinase from the marine worm Urechis caupo were determined by x-ray crystallography. One form was crystallized as a nucleotide complex, and the other was substrate-free. The two structures are similar to each other and more similar to the substrate-free forms of homologs than to the substrate-bound forms of the other phosphagen kinases. Active site specificity loop 309-317, which is disordered in substrate-free structures of homologs and is known from the NMR of arginine kinase to be inherently dynamic, is resolved in both lombricine kinase structures, providing an improved basis for understanding the loop dynamics. Phosphagen kinases undergo a segmented closing on substrate binding, but the lombricine kinase ADP complex is in the open form more typical of substrate-free homologs. Through a comparison with prior complexes of intermediate structure, a correlation was revealed between the overall enzyme conformation and the substrate interactions of His(178). Comparative modeling provides a rationale for the more relaxed specificity of these kinases, of which the natural substrates are among the largest of the phosphagen substrates.     

Origin and properties of cytoplasmic and mitochondrial isoforms of taurocyamine kinase

Taurocyamine kinase (TK) is a member of the highly conserved family of phosphagen kinases that includes creatine kinase (CK) and arginine kinase. TK is found only in certain marine annelids. In this study we used PCR to amplify two cDNAs coding for TKs from the polychaete Arenicola brasiliensis, cloned these cDNAs into the pMAL plasmid and expressed the TKs as fusion proteins with the maltose-binding protein. These are the first TK cDNA and deduced amino acid sequences to be reported. One of the two cDNA-derived amino acid sequences of TKs shows a high amino acid identity to lombricine kinase, another phosphagen kinase unique to annelids, and appears to be a cytoplasmic isoform. The other sequence appears to be a mitochondrial isoform; it has a long N-terminal extension that was judged to be a mitochondrial targeting peptide by several on-line programs and shows a higher similarity in amino acid sequence to mitochondrial creatine kinases from both vertebrates and invertebrates. The recombinant cytoplasmic TK showed activity for the substrates taurocyamine and lombricine (9% of that of taurocyamine). However, the mitochondrial TK showed activity for taurocyamine, lombricine (30% of that of taurocyamine) and glycocyamine (7% of that of taurocyamine). Neither TK catalyzed the phosphorylation of creatine. Comparison of the deduced amino acid sequences of mitochondrial CK and TK indicated that several key residues required for CK activity are lacking in the mitochondrial TK sequence. Homology models for both cytoplasmic and mitochondrial TK, constructed using CK templates, provided some insight into the structural correlation of differences in substrate specificity between the two TKs. A phylogenetic analysis using amino acid sequences from a broad spectrum of phosphagen kinases showed that annelid-specific phosphagen kinases (lombricine kinase, glycocyamine kinase and cytoplasmic and mitochondrial TKs) are grouped in one cluster, and form a sister-group with CK sequences from vertebrate and invertebrate groups. It appears that the annelid-specific phosphagen kinases, including cytoplasmic and mitochondrial TKs, evolved from a CK-like ancestor(s) early in the divergence of the protostome metazoans. Furthermore, our results suggest that the cytoplasmic and mitochondrial isoforms of TK evolved independently.     

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.     

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.     

Changing the substrate specificity of creatine kinase from creatine to glycocyamine: evidence for a highly evolved active site

Eight variants of creatine kinase were created to switch the substrate specificity from creatine to glycocyamine using a rational design approach. Changes to creatine kinase involved altering several residues on the flexible loops that fold over the bound substrates including a chimeric replacement of the guanidino specificity loop from glycocyamine kinase into creatine kinase. A maximal 2,000-fold change in substrate specificity was obtained as measured by a ratio of enzymatic efficiency (k(cat)/K(M).K(d)) for creatine vs. glycocyamine. In all cases, a change in specificity was accompanied by a large drop in enzymatic efficiency. This data, combined with evidence from other studies, indicate that substrate specificity in the phosphagen kinase family is obtained by precise alignment of substrates in the active site to maximize k(cat)/K(M).K(d) as opposed to selective molecular recognition of one guanidino substrate over another. A model for the evolution of the dimeric forms of phosphagen kinases is proposed in which these enzymes radiated from a common ancestor that may have possessed a level of catalytic promiscuity. As mutational events occurred leading to greater degrees of substrate specificity, the dimeric phosphagen kinases became evolutionary separated such that the substrate specificity could not be interchanged by a small number of mutations.     

The role of phosphagen specificity loops in arginine kinase

Phosphagen kinases catalyze the reversible transfer of a phosphate between ATP and guanidino substrates, a reaction that is central to cellular energy homeostasis. Members of this conserved family include creatine and arginine kinases and have similar reaction mechanisms, but they have distinct specificities for different guanidino substrates. There has not been a full structural rationalization of specificity, but two loops have been implicated repeatedly. A small domain loop is of length that complements the size of the guanidino substrate, and is located where it could mediate a lock-and-key mechanism. The second loop contacts the substrate with a valine in the methyl-substituted guanidinium of creatine, and with a glutamate in the unsubstituted arginine substrate, leading to the proposal of a discriminating hydrophobic/hydrophilic minipocket. In the present work, chimeric mutants were constructed with creatine kinase loop elements inserted into arginine kinase. Contrary to the prior rationalizations of specificity, most had measurable arginine kinase activity but no creatine kinase activity or enhanced phosphocreatine binding. Guided by structure, additional mutations were introduced in each loop, recovering arginine kinase activities as high as 15% and 64% of wild type, respectively, even though little activity would be expected in the constructs if the implicated sites had dominant roles in specificity. An atomic structure of the mismatched complex of arginine kinase with creatine and ADP indicates that specificity can also be mediated by an active site that allows substrate prealignment that is optimal for reactivity only with cognate substrates and not with close homologs that bind but do not react.     

Isoleucine 69 and valine 325 form a specificity pocket in human muscle creatine kinase

Creatine kinase (CK) catalyzes the reversible phosphorylation of creatine by ATP. From a structural perspective, the enzyme utilizes two flexible loop regions to sequester and position the substrates for catalysis. There has been debate over the specific roles of the flexible loops in substrate specificity and catalysis in CK and other related phosphagen kinases. In CK, two hydrophobic loop residues, I69 and V325, make contacts with the N-methyl group of creatine. In this study, we report the alteration of the substrate specificity of CK through the mutagenesis of V325. The V325 to glutamate mutation results in a more than 100-fold preference for glycocyamine, while mutation of V325 to alanine results in a slight preference of the enzyme for cyclocreatine (1-carboxymethyl-2-iminoimidazolidine). This study enhances our understanding of how the active sites of phosphagen kinases have evolved to recognize their respective substrates and catalyze their reactions.     

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.     

Crystal structure of brain-type creatine kinase at 1.41 A resolution

Excitable cells and tissues like muscle or brain show a highly fluctuating consumption of ATP, which is efficiently regenerated from a large pool of phosphocreatine by the enzyme creatine kinase (CK). The enzyme exists in tissue--as well as compartment-specific isoforms. Numerous pathologies are related to the CK system: CK is found to be overexpressed in a wide range of solid tumors, whereas functional impairment of CK leads to a deterioration in energy metabolism, which is phenotypic for many neurodegenerative and age-related diseases. The crystal structure of chicken cytosolic brain-type creatine kinase (BB-CK) has been solved to 1.41 A resolution by molecular replacement. It represents the most accurately determined structure in the family of guanidino kinases. Except for the N-terminal region (2-12), the structures of both monomers in the biological dimer are very similar and closely resemble those of the other known structures in the family. Specific Ca2+-mediated interactions, found between two dimers in the asymmetric unit, result in structurally independent heterodimers differing in their N-terminal conformation and secondary structure. The high-resolution structure of BB-CK presented in this work will assist in designing new experiments to reveal the molecular basis of the multiple isoform-specific properties of CK, especially regarding different subcellular locations and functional interactions with other proteins. The rather similar fold shared by all known guanidino kinase structures suggests a model for the transition state complex of BB-CK analogous to the one of arginine kinase (AK). Accordingly, we have modeled a putative conformation of CK in the transition state that requires a rigid body movement of the entire N-terminal domain by rms 4 A from the structure without substrates.     

Origin of the genes for the isoforms of creatine kinase

Creatine kinase (CK) is a member of a family of phosphoryl transfer enzymes called phosphagen (guanidino) kinases which play a central role in cellular energy homeostasis. There are three CK isoform gene groups, each coding for proteins targeted to different intracellular compartments--cytoplasmic (CytCK), mitochondrial (MtCK) and flagellar (FlgCK). The former two CKs are either dimeric or octameric while FlgCKs are contiguous trimers consisting of three fused, complete CK domains. Conventional wisdom supports the view that CKs evolved from a cytoplasmic, monomeric ancestral protein closely related to a phosphagen kinase homologue, arginine kinase (AK). Recently, it has been shown that a demosponge (Phylum Porifera) expresses a true MtCK and two dimeric, protoflagellar CKs (protoflgCK) with great similarity to FlgCKs. To further probe the early evolution of CK, we have obtained additional sequences for Mt- and protoflgCKs from two more demosponges and from three hexactinellid (glass) sponges as well as an MtCK sequence from a basal metazoan cnidarian. Phylogenetic analyses using Maximum Likelihood (ML) of these new CK sequences with other CKs and phosphagen kinases yielded a consensus tree containing an assemblage of MtCKs and a supercluster consisting of protoflg-, Flg- and CytCKs. The MtCKs appear basal in the tree topology consistent with prior results. Within the protoflg-, Flg- and CytCK supercluster, the protoflgCKs appear to be allied to the domains of the FlgCKs, although the support is not robust. PCR amplification of genomic DNA and sequencing of the genes for Mt- and protoflgCK from the demosponge Suberites fuscus showed that the sponge MtCK shares four-five common intron:exon boundaries with invertebrate, protochordate and vertebrate MtCKs supporting a common ancestry and the extreme conservation of intron:exon organization in MtCK genes. The protoflgCK gene organization was highly divergent in relation to other CK genes but shares a common intron:exon boundary with domain 2 of the gene for the FlgCK from the tunicate Ciona intestinalis, providing support for the linkage of the protoflgCKs with the FlgCKs. Our results show that the two, major CK gene lineages are present in arguably the oldest, extant metazoan group, the hexactinellid sponges, indicating that these two genes are ancient and confirming prior work that the MtCK gene is likely basal and ancestral.     

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.     

Induced fit in guanidino kinases--comparison of substrate-free and transition state analog structures of arginine kinase

Arginine kinase (AK) is a member of the guanidino kinase family that plays an important role in buffering ATP concentration in cells with high and fluctuating energy demands. The AK specifically catalyzes the reversible phosphoryl transfer between ATP and arginine. We have determined the crystal structure of AK from the horseshoe crab (Limulus polyphemus) in its open (substrate-free) form. The final model has been refined at 2.35 A with a final R of 22.3% (R(free) = 23.7%). The structure of the open form is compared to the previously determined structure of the transition state analog complex in the closed form. Classically, the protein would be considered two domain, but dynamic domain (DynDom) analysis shows that most of the differences between the two structures can be considered as the motion between four rigid groups of nonsequential residues. ATP binds near a cluster of positively charged residues of a fixed dynamic domain. The other three dynamic domains close the active site with separate hinge rotations relative to the fixed domain. Several residues of key importance for the induced motion are conserved within the phosphagen kinase family, including creatine kinase. Substantial conformational changes are induced in different parts of the enzyme as intimate interactions are formed with both substrates. Thus, although induced fit occurs in a number of phosphoryl transfer enzymes, the conformational changes in phosphagen kinases appear to be more complicated than in prior examples.     

Expression, purification from inclusion bodies, and crystal characterization of a transition state analog complex of arginine kinase: a model for studying phosphagen kinases.

Phosphagen kinases catalyze the reversible transfer of a phosphoryl group between guanidino phosphate compounds and ADP, thereby regenerating ATP during bursts of cellular activity. Large quantities of highly pure arginine kinase (EC 2.7.3.3), the phosphagen kinase present in arthropods, have been isolated from E. coli, into which the cDNA for the horseshoe crab enzyme had been cloned. Purification involves size exclusion and anion exchange chromatographies applied in the denatured and refolded states. The recombinant enzyme has been crystallized as a transition state analog complex. Near complete native diffraction data have been collected to 1.86 A resolution. Substitution of a recombinant source for a natural one, improvement in the purification, and data collection at cryo temperatures have all yielded significant improvements in diffraction.     

Characterization of a novel bacterial arginine kinase from Desulfotalea psychrophila

Phosphagen kinases are found throughout the animal kingdom and catalyze the transfer of a high-energy gamma phosphoryl-group from ATP to a guanidino group on a suitable acceptor molecule such as creatine or arginine. Recent genome sequencing efforts in several proteobacteria, including Desulfotalea psychrophila LSv54, Myxococcus xanthus, Sulfurovum sp. NBC37-1, and Moritella sp. PE36 have revealed what appears to be a phosphagen kinase homolog present in their genomes. Based on sequence comparisons these putative homologs bear a strong resemblance to arginine kinases found in many invertebrates and some protozoa. We describe here a biochemical characterization of one of these homologs from D. psychrophila expressed in E. coli that confirms its ability to reversibly catalyze phosphoryl transfer from ATP to arginine. A phylogenetic analysis suggests that these bacteria homologs are not widely distributed in proteobacteria species. They appear more related to protozoan arginine kinases than to similar proteins seen in some Gram-positive bacteria that share key catalytic residues but encode protein tyrosine kinases. This raises the possibility of horizontal gene transfer as a likely origin of the bacterial arginine kinases.