Isolation, characterization, and cDNA-derived amino acid sequence of glycocyamine kinase from the tropical marine worm Namalycastis sp

The characterization of novel monomeric creatine kinases in the early branching Alveolata species, Perkinsus marinus: Implications for phosphagen kinase evolution

The genome of the unicellular molluscan parasite Perkinsus marinus contains at least five genes coding for putative creatine kinases (CK), a phosphoryl transfer enzyme which plays a key role in cellular energy transactions. Expression and kinetic analyses of three of the P. marinus CKs revealed them to be true CKs with catalytic properties in the range of typical metazoan CKs. A sequence comparison of the P. marinus CKs with a range of CK dimers and other dimeric phosphoryl transfer enzymes in this family (phosphagen kinases) showed that the P. marinus CKs lacked some of the critical residues involved in dimer stabilization, a trait all previously characterized CKs share. Size exclusion chromatography of all three expressed P. marinus CK constructs indicated they are monomeric, consistent with the observed lack of some critical dimer stabilizing residues. Phylogenetic analyses of the P. marinus CKs and putative dinoflagellate CKs with a broad range of monomeric and dimeric phosphagen kinases revealed that the Perkinsus CKs form a distinct, well-supported clade with dinoflagellate CKs which also lack the dimer stabilizing residues. Analysis of the genomic data for P. marinus showed the presence of putative genes for the two enzymes associated with creatine biosynthesis. CK in higher organisms plays a critical role in energy buffering in cell types displaying high and variable rates of ATP turnover. The presence of multiple CKs and the creatine biosynthetic pathway in P. marinus indicates that this unicellular parasite has the full complement of molecular machinery for CK-mediated energy buffering.

Diversity of phosphagen kinases in annelids: The first sequence report for a putative opheline kinase

Opheline kinase (OK) is one of the phosphagen kinases (PKs) restricted to annelids, but the amino acid sequence has not been determined yet. The OK enzyme was isolated in 1966 from the polychaete Ophelia neglecta (Opheliidae) and shown to have somewhat broader activities for the various substrates opheline, lombricine and taurocyamine. To determine the OK sequence, we analyzed the RNA sequencing data for Ophelina sp. and Thoracophelia sp., belonging to Opheliidae. Four PK sequences, namely, taurocyamine kinase (TK), creatine kinase (CK), mitochondrial CK (MiCK) and putative OK, were identified in both species, and the recombinant Ophelina enzymes were expressed in E. coli and purified. Since the substrate opheline was not commercially available, we used the partial activity toward taurocyamine to infer the enzyme specificity. The putative Ophelina OK showed lower activity to taurocyamine with a V_max/K_m nearly identical to a previously published value for an OK from a related species Ophelia neglecta. Under the same conditions, the true Ophelina TK showed much higher activity. Thus, the putative Ophelina enzyme was determined to be OK. The amino acid sequence alignment indicated that Ophelina and Thoracophelia OKs have five amino acid deletions in the GS region, like those of LKs and AKs, and the guanidino substrate specific residue was Lys, the same as LKs. In the phylogenetic tree constructed from annelid PK amino acid sequences, the OK sequences formed a distinct cluster, and it was placed near the TK and lombricine kinase (LK) clusters. This is the first report of the amino acid sequence for the OK enzyme.

Natural Products Containing 'Rare' Organophosphorus Functional Groups

Phosphorous-containing molecules are essential constituents of all living cells. While the phosphate functional group is very common in small molecule natural products, nucleic acids, and as chemical modification in protein and peptides, phosphorous can form P⁻N (phosphoramidate), P⁻S (phosphorothioate), and P⁻C (e.g., phosphonate and phosphinate) linkages. While rare, these moieties play critical roles in many processes and in all forms of life. In this review we thoroughly categorize P⁻N, P⁻S, and P⁻C natural organophosphorus compounds. Information on biological source, biological activity, and biosynthesis is included, if known. This review also summarizes the role of phosphorylation on unusual amino acids in proteins (N- and S-phosphorylation) and reviews the natural phosphorothioate (P⁻S) and phosphoramidate (P⁻N) modifications of DNA and nucleotides with an emphasis on their role in the metabolism of the cell. We challenge the commonly held notion that nonphosphate organophosphorus functional groups are an oddity of biochemistry, with no central role in the metabolism of the cell. We postulate that the extent of utilization of some phosphorus groups by life, especially those containing P⁻N bonds, is likely severely underestimated and has been largely overlooked, mainly due to the technological limitations in their detection and analysis.

Characterization of a putative oomycete taurocyamine kinase: Implications for the evolution of the phosphagen kinase family

Phosphagen kinases (PKs) are known to be distributed throughout the animal kingdom, but have recently been discovered in some protozoan and bacterial species. Within animal species, these enzymes play a critical role in energy homeostasis by catalyzing the reversible transfer of a high-energy phosphoryl group from Mg⋅ATP to an acceptor molecule containing a guanidinium group. In this work, a putative PK gene was identified in the oomycete Phytophthora sojae that was predicted, based on sequence homology, to encode a multimeric hypotaurocyamine kinase. The recombinant P. sojae enzyme was purified and shown to catalyze taurocyamine phosphorylation efficiently (kcat/KM (taurocyamine) = 2 × 10(5) M(-1) s(-1)) and glycocyamine phosphorylation only weakly (kcat/KM (glycocyamine) = 2 × 10(2) M(-1) s(-1)), but lacked any observable kinase activity with the more ubiquitous guanidinium substrates, creatine or arginine. Additionally, the enzyme was observed to be dimeric but lacked cooperativity between the subunits in forming a transition state analog complex. These results suggest that protozoan PKs may exhibit more diversity in substrate specificity than was previously thought.

A diverse array of creatine kinase and arginine kinase isoform genes is present in the starlet sea anemone Nematostella vectensis, a cnidarian model system for studying developmental evolution

Phosphagen (guanidino) kinases (PK) constitute a family of homologous phosphotransferases catalyzing the reversible transfer of the high-energy phosphoryl group of ATP to naturally occurring guanidine compounds. Prior work has shown that PKs can be phylogenetically separated into two distinct groups- an arginine kinase (AK) subfamily and a creatine kinase (CK) subfamily. The latter includes three CK isoforms- cytoplasmic CK (CyCK), mitochondrial CK (MiCK) and three-domain flagellar CK (fCK). In the present study we identified six unique PK genes from the draft genome sequence of the starlet sea anemone Nematostella vectensis, a well-known model organism for understanding metazoan developmental evolution. Using reverse transcription polymerase chain reaction (RTPCR) methods, full length cDNAs were amplified for all of these PKs. These cDNAs were cloned and expressed in Escherichia coli as 6x His-tagged fusion proteins. The six PKs were identified as the three typical CK isoforms (CyCK, MiCK and fCK), two unusual AKs (a two-domain AK (2DAK) and a three-domain AK (3DAK)) and a PK which phosphorylated arginine. The latter enzyme had a very low AK activity (its apparent V(max) value being less than 0.2% that of 3DAK), lacks several key residues necessary for AK enzyme activity, and was tentatively designated as AK1. As far as we know, this constitutes the first report of an AK with the three fused AK domains. The Bayesian tree suggested that the third domain of 3DAK likely evolved from the gene for domain 2 of typical two-domain AK found widely in cnidarians. Construction of phylogenetic trees and comparison of exon-intron organizations of their respective genes indicated that the N. vectensis three-domain fCK and 3DAK evolved independently, and both enzymes are likely to be targeted to cell membranes since they have a myristoylation signal at their respective N-termini. These results complement prior work on other basal invertebrates showing that multiple CK and AK isoform genes were present at the dawn of the radiation of metazoans. The presence of isoform diversity in an organism lacking in structural complexity reflects an early imperative for targeting of PKs to particular cellular contexts such as muscle fibers, neurons, ciliated/flagellated epithelia and spermatozoa.

Identification of amino acid residues responsible for taurocyamine binding in mitochondrial taurocyamine kinase from Arenicola brasiliensis

In order to investigate the residues associated with binding of the substrate taurocyamine in Arenicola mitochondrial taurocyamine kinase (TK), we performed Ala-scanning of the amino acid sequence HTKTV at positions 67-71 on the GS loop, and determined apparent K(m) and V(max) (appK(m) and appV(max), respectively) of the mutant forms for the substrates taurocyamine and glycocyamine. The appK(m) values for taurocyamine of the K69A, T70A and V71A mutants were significantly increased as compared with wild-type, suggesting that these residues are associated with taurocyamine binding. Of special interest is a property of V71A mutant: its catalytic efficiency for glycocyamine was twice that for taurocyamine, indicating that the V71A mutant acts like a glycocyamine kinase, rather than a TK. The role of the amino acid residue K95 of Arenicola MiTK was also examined. K95 was replaced with R, H, Y, I, A and E. K95R, K95H and K95I have a 3-fold higher affinity for taurocyamine, and activity was largely lost in K95E. On the other hand, the K95Y mutant showed a rather unique feature; namely, an increase in substrate concentration caused a decrease in initial velocity of the reaction (substrate inhibition). This is the first report on the key amino acid residues responsible for taurocyamine binding in mitochondrial TK.

Structural basis for the mechanism and substrate specificity of glycocyamine kinase, a phosphagen kinase family member

Glycocyamine kinase (GK), a member of the phosphagen kinase family, catalyzes the Mg(2+)-dependent reversible phosphoryl group transfer of the N-phosphoryl group of phosphoglycocyamine to ADP to yield glycocyamine and ATP. This reaction helps to maintain the energy homeostasis of the cell in some multicelullar organisms that encounter high and variable energy turnover. GK from the marine worm Namalycastis sp. is heterodimeric, with two homologous polypeptide chains, alpha and beta, derived from a common pre-mRNA by mutually exclusive N-terminal alternative exons. The N-terminal exon of GKbeta encodes a peptide that is different in sequence and is 16 amino acids longer than that encoded by the N-terminal exon of GKalpha. The crystal structures of recombinant GKalphabeta and GKbetabeta from Namalycastis sp. were determined at 2.6 and 2.4 A resolution, respectively. In addition, the structure of the GKbetabeta was determined at 2.3 A resolution in complex with a transition state analogue, Mg(2+)-ADP-NO(3)(-)-glycocyamine. Consistent with the sequence homology, the GK subunits adopt the same overall fold as that of other phosphagen kinases of known structure (the homodimeric creatine kinase (CK) and the monomeric arginine kinase (AK)). As with CK, the GK N-termini mediate the dimer interface. In both heterodimeric and homodimeric GK forms, the conformations of the two N-termini are asymmetric, and the asymmetry is different than that reported previously for the homodimeric CKs from several organisms. The entire polypeptide chains of GKalphabeta are structurally defined, and the longer N-terminus of the beta subunit is anchored at the dimer interface. In GKbetabeta the 24 N-terminal residues of one subunit and 11 N-terminal residues of the second subunit are disordered. This observation is consistent with a proposal that the GKalphabeta amino acids involved in the interface formation were optimized once a heterodimer emerged as the physiological form of the enzyme. As a consequence, the homodimer interface (either solely alpha or solely beta chains) has been corrupted. In the unbound state, GK exhibits an open conformation analogous to that observed with ligand-free CK or AK. Upon binding the transition state analogue, both subunits of GK undergo the same closure motion that clasps the transition state analogue, in contrast to the transition state analogue complexes of CK, where the corresponding transition state analogue occupies only one subunit, which undergoes domain closure. The active site environments of the GK, CK, and AK at the bound states reveal the structural determinants of substrate specificity. Despite the equivalent binding in both active sites of the GK dimer, the conformational asymmetry of the N-termini is retained. Thus, the coupling between the structural asymmetry and negative cooperativity previously proposed for CK is not supported in the case of GK.

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.

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.

Phosphagen kinase of the giant tubeworm Riftia pachyptila. Cloning and expression of cytoplasmic and mitochondrial isoforms of taurocyamine kinase

The giant tubeworm Riftia pachyptila lives at deep-sea hydrothermal vents along the East Pacific Rise and the Galapagos Rift. The large size and high growth rate of R. pachyptila is supported by an endosymbiotic relationship with a chemosynthetic bacterium. Elucidation of the regulation of energy metabolism of the giant tubeworm remains an interesting problem. The purpose of this study is to determine the cDNA sequence of phosphagen kinase, one of the most important enzymes in energy metabolism, and to characterize its function. Two phosphagen kinase cDNA sequences amplified from the cDNA library of R. pachyptila showed high derived amino acid sequence identity (74%) with those of cytoplasmic taurocyamine kinase (TK) and mitochondrial TK from an annelid Arenicola brasiliensis. The cytoplasmic form of the Riftia recombinant enzyme showed stronger activity for the substrates taurocyamine and also considerable activity for lombricine (21% that of taurocyamine). The mitochondrial form, which was structurally similar to mitochondrial creatine kinase, showed stronger activity for taurocyamine, and a broader activity for various guanidine compounds: glycocyamine (35% that of taurocyamine), lombricine (31%) and arginine (3%). Both forms showed no activity for creatine. The difference in substrate specificities between the cytoplasmic and mitochondrial forms might be attributable to the large difference in the amino acid sequence of the GS region and/or several key amino acid residues for establishing guanidine substrate specificity. Based on these results, we conclude that Riftia contains at least two forms of TK as phosphagen kinase. We also report the kinetic parameters, Km and kcat, of Arenicola and Riftia TKs for the first time. The Km values for taurocyamine of Arenicola and Riftia TKs ranged from 0.9 to 4.0 mM and appear to be comparable to those of other annelid-specific enzymes, lombricine kinase and glycocyamine kinase, but are significantly lower than those of Neanthes cytoplasmic and mitochondrial creatine kinases. Comparison of kcat/Km value in various annelid phosphagen kinases indicates that Arenicola mitochondrial TK has the highest catalytic efficiency (16.2 s-1 mM-1). In Arenicola TKs, the mitochondrial form has seven-fold higher efficiency than the cytoplasmic form.

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.

Amino acid residues 62 and 193 play the key role in regulating the synergism of substrate binding in oyster arginine kinase

The purpose of this study is to clarify the amino acid residues responsible for the synergism in substrate binding of arginine kinase (AK), a key enzyme in invertebrate energy metabolism. AKs contain a pair of highly conserved amino acids (D62 and R193) that form an ion pair, and replacement of these residues can cause a pronounced loss of activity. Interestingly, in the oyster Crassostrea AK, these residues are replaced by an N and a K, respectively. Despite this replacement, the enzyme retains high activity and moderate synergism in substrate binding (Kd/Km=2.3). We replaced the N62 by G or D and the K193 by G or R in Crassostrea AK, and also constructed the double mutants of N62G/K193G and N62D/K193R. All of the mutants retained 50-90% of the wild-type activity. In N62G and N62D mutants, the Kmarg for arginine binding was comparable to that of wild-type enzyme, but the Kdarg was increased 2-5-fold, resulting in a strong synergism (Kd/Km=4.9-11.3). On the other hand, in K193G and K193R mutants, the Kmarg was increased 4-fold, and synergism was lost almost completely (Kd/Km=1.0-1.4). The N62G/K193G double mutant showed similar characteristics to the K193G and K193R mutants. Another double mutant, N62D/K193R, similar to the amino acid pair in the wild-type enzyme, had characteristics similar to those of the wild-type enzyme. These results indicate that the amino acid residues 62 and 193 play the key role in mediating the synergism in substrate binding of oyster arginine kinase.