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

Novel insights on obligate symbiont lifestyle and adaptation to chemosynthetic environment as revealed by the giant tubeworm genome

The mutualism between the giant tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has been extensively researched over the past 40 years. However, the lack of the host whole genome information has impeded the full comprehension of the genotype/phenotype interface in Riftia. Here we described the high-quality draft genome of Riftia, its complete mitogenome, and tissue-specific transcriptomic data. The Riftia genome presents signs of reductive evolution, with gene family contractions exceeding expansions. Expanded gene families are related to sulphur metabolism, detoxification, anti-oxidative stress, oxygen transport, immune system, and lysosomal digestion, reflecting evolutionary adaptations to the vent environment and endosymbiosis. Despite the derived body plan, the developmental gene repertoire in the gutless tubeworm is extremely conserved with the presence of a near intact and complete Hox cluster. Gene expression analyses establishes that the trophosome is a multi-functional organ marked by intracellular digestion of endosymbionts, storage of excretory products and haematopoietic functions. Overall, the plume and gonad tissues both in contact to the environment harbour highly expressed genes involved with cell cycle, programmed cell death, and immunity indicating a high cell turnover and defence mechanisms against pathogens. We posit that the innate immune system plays a more prominent role into the establishment of the symbiosis during the infection in the larval stage, rather than maintaining the symbiostasis in the trophosome. This genome bridges four decades of physiological research in Riftia, whilst simultaneously provides new insights into the development, whole organism functions and evolution in the giant tubeworm.

Insight into Structural Aspects of Histidine 284 of Daphnia magna Arginine Kinase

Arginine kinase (AK), a bioenergy-related enzyme, is distributed widely in invertebrates. The role of highly conserved histidines in AKs is still unascertained. In this study, the highly conserved histidine 284 (H284) in AK of Daphnia magna (DmAK) was replaced with alanine to elucidate the role of H284. We examined the alteration of catalytic activity and structural changes of H284A in DmAK. The catalytic activity of H284A was reduced dramatically compared to that in wild type (WT). Thus the crystal structure of H284A displayed several structural changes, including the alteration of D324, a hydrogen-bonding network around H284, and the disruption of π-stacking between the imidazole group of the H284 residue and the adenine ring of ATP. These findings suggest that such alterations might affect a conformational change of the specific loop consisting of G310-V322 at the antiparallel β-sheet region. Thus, we speculated that the H284 residue might play an important role in the conformational change of the specific loop when ATP binds to the substrate-binding site of DmAK.

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.

Molecular cloning and characterization of taurocyamine kinase from Clonorchis sinensis: a candidate chemotherapeutic target

Background

Adult Clonorchis sinensis lives in the bile duct and causes endemic clonorchiasis in East Asian countries. Phosphagen kinases (PK) constitute a highly conserved family of enzymes, which play a role in ATP buffering in cells, and are potential targets for chemotherapeutic agents, since variants of PK are found only in invertebrate animals, including helminthic parasites. This work is conducted to characterize a PK from C. sinensis and to address further investigation for future drug development.

Methology/Principal findings

A cDNA clone encoding a putative polypeptide of 717 amino acids was retrieved from a C. sinensis transcriptome. This polypeptide was homologous to taurocyamine kinase (TK) of the invertebrate animals and consisted of two contiguous domains. C. sinensis TK (CsTK) gene was reported and found consist of 13 exons intercalated with 12 introns. This suggested an evolutionary pathway originating from an arginine kinase gene group, and distinguished annelid TK from the general CK phylogenetic group. CsTK was found not to have a homologous counterpart in sequences analysis of its mammalian hosts from public databases. Individual domains of CsTK, as well as the whole two-domain enzyme, showed enzymatic activity and specificity toward taurocyamine substrate. Of the CsTK residues, R58, I60 and Y84 of domain 1, and H60, I63 and Y87 of domain 2 were found to participate in binding taurocyamine. CsTK expression was distributed in locomotive and reproductive organs of adult C. sinensis. Developmentally, CsTK was stably expressed in both the adult and metacercariae stages. Recombinant CsTK protein was found to have low sensitivity and specificity toward C. sinensis and platyhelminth-infected human sera on ELISA.

Conclusion

CsTK is a promising anti-C. sinensis drug target since the enzyme is found only in the C. sinensis and has a substrate specificity for taurocyamine, which is different from its mammalian counterpart, creatine.

The food-borne clonorchiasis imposes public health problems on inhabitants in endemic areas. Praziquantel has been employed as an efficacious anthelminthic in large-scale campaigns as well as for individual treatment of Clonorchis sinensis human infections. Although praziquantel continues to have good efficacy, new drug development for this parasite has been recognized as a crucial issue to be investigated intensively. Clonorchis sinensis adults generate energy through glycolysis, actively utilizing exogenous glucose, and produce a large amount of eggs each day. Taurocyamine kinase (CsTK) is distributed abundantly in the locomotive and reproductive organs, and is an important enzyme in energy generation and homeostasis in adult C. sinensis. Enzymes of the glycolytic pathway are also expressed abundantly in these organs and in tegument, implying these organs play central roles which are essential for survival and reproduction of C. sinensis. The TK enzymes, including CsTK, are found only among invertebrate organisms and have substrate specificity for taurocyamine, which are significantly different from phosphagen kinases of vertebrate animals. With these molecular biological, enzymatic, and evolutionary characteristics, we propose here that CsTK could be a target for development of chemotherapeutic agents against C. sinensis and be a biomolecular model for other human-infecting trematodes.

The role of Y84 on domain 1 and Y87 on domain 2 of Paragonimus westermani taurocyamine kinase: Insights on the substrate binding mechanism of a trematode phosphagen kinase

The two-domain taurocyamine kinase (TK) from Paragonimus westermani was suggested to have a unique substrate binding mechanism. We performed site-directed mutagenesis on each domain of this TK and compared the kinetic parameters Km(Tc) and Vmax with that of the wild-type to determine putative amino acids involved in substrate recognition and binding. Replacement of Y84 on domain 1 and Y87 on domain 2 with R resulted in the loss of activity for the substrate taurocyamine. Y84E mutant has a dramatic decrease in affinity and activity for taurocyamine while Y87E has completely lost catalytic activity. Substituting H and I on the said positions also resulted in significant changes in activity. Mutation of the residues A59 on the GS region of domain 1 also caused significant decrease in affinity and activity while mutation on the equivalent position on domain 2 resulted in complete loss of activity.

Gene structure of the two-domain taurocyamine kinase from Paragonimus westermani: evidence for a distinct lineage of trematode phosphagen kinases

Taurocyamine kinase (TK) is an enzyme that catalyzes the reversible transfer of a phosphate between ATP and taurocyamine. Annelid TKs were suggested to have evolved from a CK ancestor. However, TKs from the lung fluke Paragonimus westermani comprised another lineage. Construction of phylogenetic tree and comparison of exon/intron organization showed that P. westermani TK and other trematode TKs evolved from a molluscan arginine kinase (AK) gene. Exon shuffling probably caused the changes in amino acid sequence thereby changing the affinity from AK to TK. The present study provides new insights on the evolution of phosphagen kinases found in trematodes.

A novel taurocyamine kinase found in the protist Phytophthora infestans

Phosphagen kinase (PK), which is typically in the form of creatine kinase (CK; EC 2.7.3.2) in vertebrates or arginine kinase (AK; EC 2.7.3.3) in invertebrates, plays a key role in ATP buffering systems of tissues and nerves that display high and variable rates of ATP turnover. The enzyme is also found with intermittent occurrence as AK in unicellular organisms, protist and bacteria species, suggesting an ancient origin of AK. Through a database search, we identified two novel PK genes, coding 40- and 80-kDa (contiguous dimer) enzymes in the protist Phytophthora infestans. Both enzymes showed strong activity for taurocyamine and, in addition, we detected taurocyamine in cell extracts of P. infestans. Thus, the enzyme was identified to be taurocyamine kinase (TK; EC 2.7.3.4). This was the first phosphagen kinase, other than AK, to be found in unicellular organisms. Their position on the phylogenetic tree indicates that P. infestans TKs evolved uniquely at an early stage of evolution. Occurrence of TK in protists suggests that PK enzymes show flexible substrate specificity.

Cold-adapted features of arginine kinase from the deep-sea clam Calyptogena kaikoi

The heterodont clam Calyptogena kaikoi, which inhabits depths exceeding 3,500 m where low ambient temperatures prevail, has an unusual two-domain arginine kinase (AK) with molecular mass of 80 kDa, twice that of typical AKs. The purpose of this work is to investigate the nature of the adaptations of this AK for functioning at low temperatures. Recombinant C. kaikoi AK constructs were expressed, and their two-substrate kinetic constants (k(cat), K(a), and K(ia)) were determined at 10°C and 25°C, respectively. When measured at 25°C, the K(ia) values were tenfold larger than those for corresponding K(a) values, while at 10°C, the K(ia) values decreased remarkably, but the K (a) values were almost unchanged. The Calyptogena two-domain enzyme has threefold higher catalytic efficiency, calculated by k (cat)/(K(a)(ARG)·K(ia)(ATP) ), at 10°C, than that at 25°C, reflecting adaptation for function at reduced ambient temperatures. The activation energy (E(a)) and thermodynamic parameters were determined for Calyptogena two-domain enzyme and compared with those of two-domain enzymes from mesophilic Corbicula and Anthopleura. The value for E(a) of Calyptogena enzyme were about half of those for mesophilic enzymes, and a larger decrease in entropy was observed in Calyptogena AK reaction. Although large decrease in entropy increases the ΔG(o‡) value and consequently lowers the k(cat) value, this is compensated with its lower E(a) value thereby minimizing the reduction in its k(cat) value. These thermodynamic properties, together with the kinetic ones, are also present in the separated domain 2 of the Calyptogena two-domain enzyme.

Linking hydrothermal geochemistry to organismal physiology: physiological versatility in Riftia pachyptila from sedimented and basalt-hosted vents

Much of what is known regarding Riftia pachyptila physiology is based on the wealth of studies of tubeworms living at diffuse flows along the fast-spreading, basalt-hosted East Pacific Rise (EPR). These studies have collectively suggested that Riftia pachyptila and its chemoautotrophic symbionts are physiologically specialized, highly productive associations relying on hydrogen sulfide and oxygen to generate energy for carbon fixation, and the symbiont's nitrate reduction to ammonia for energy and biosynthesis. However, Riftia also flourish in sediment-hosted vents, which are markedly different in geochemistry than basalt-hosted systems. Here we present data from shipboard physiological studies and global quantitative proteomic analyses of Riftia pachyptila trophosome tissue recovered from tubeworms residing in the EPR and the Guaymas basin, a sedimented, hydrothermal vent field. We observed marked differences in symbiont nitrogen metabolism in both the respirometric and proteomic data. The proteomic data further suggest that Riftia associations in Guaymas may utilize different sulfur compounds for energy generation, may have an increased capacity for energy storage, and may play a role in degrading exogenous organic carbon. Together these data reveal that Riftia symbionts are far more physiologically plastic than previously considered, and that -contrary to previous assertions- Riftia do assimilate reduced nitrogen in some habitats. These observations raise new hypotheses regarding adaptations to the geochemical diversity of habitats occupied by Riftia, and the degree to which the environment influences symbiont physiology and evolution.

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.

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.

Crystallization and X-ray analysis of the Schistosoma mansoni guanidino kinase

The 716-amino-acid guanidino kinase from the parasitic flatworm Schistosoma mansoni results from the fusion of two guanidino kinase subunits. Crystals of this 80 kDa protein have been obtained in the monoclinic space group P2(1), with unit-cell parameters a = 52.7, b = 122.1, c = 63.2 A, beta = 108.5 degrees . Synchrotron data were collected to 2.8 A resolution on ESRF beamline ID29. The structure was solved by the molecular-replacement method, using the 357-amino-acid structure of the arginine kinase from Trypanosoma cruzi as the search model.

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.

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.

Comparison of kinetic constants of creatine kinase isoforms

The purpose of this study was to elucidate the functional differences between the CK isoforms by cloning the cDNAs of 12 CK isoforms: the M and B cytoplasmic forms and uMiCK from mouse, the M1, M2 and B cytoplasmic forms from Danio rerio, M1 and M2 cytoplasmic forms from the lower vertebrate Lampetra japonica, a cytoplasmic CK and a MiCK from the marine worm Neanthes diversicolor, and a cytoplasmic CK and a MiCK from the soft coral Dendronephthya gigantea. These were expressed in Escherichia coli as a fusion protein with maltose-binding protein, and kinetic constants (K(m), K(d) and k(cat)) of all the recombinant enzymes, except for the unstable Dendronephthya cytoplasmic CK, were determined for the forward reaction. The kinetic constants of the M- and B-forms of the mouse and Danio cytoplasmic CKs differed significantly, with the K(m) for creatine (K(m)Cr) of M-CK being three- to nine-fold higher than that of B-CK, possibly reflecting differences in the concentration of creatine in muscle and brain cells. The mouse uMiCK had the lowest K(m)Cr value among the CK isoforms. In addition, it also exhibited a strong synergism for substrate binding (K(d)/K(m)=11.8). These results indicate that uMiCK has unique characteristics compared with other CK isoforms. Two subisoforms of M-CK were found in the lower vertebrate L. japonica, and the kinetic constants of recombinant M1- and M2-CKs differed significantly. The M1- and M2-CKs were expressed in skeletal muscle with a ratio of 7:3, while M1-CK was the predominant subisoform in the testis. The kinetic constants of cytoplasmic CK from the marine worm Neanthes were significantly different from those of Neanthes MiCK, possibly indicating that functional differences among CK isoforms occurred at least before the divergence of annelids from other protostome invertebrates.

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.