Evolution of the cytoplasmic and mitochondrial phosphagen kinases unique to annelid groups

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.

Phosphagen kinases from five groups of eukaryotic protists (Choanomonada, Alveolate, Stramenopiles, Haptophyta, and Cryptophyta): Diverse enzyme activities and phylogenetic relationship with metazoan enzymes

Among 28 groups of eukaryotes, apart from Metazoa, phosphagen kinase (PKs) is distributed in only a few protist groups, including the Choanomonada with the closest affinity to metazoans. To clarify the origin of metazoan PKs, we performed a database search and focused on 11 sequences of PK homologs from five groups of protists: the Choanomonada, Alveolata, Haptophyta, Stramenopiles, and Cryptophyta. The recombinant enzymes were prepared to determine their substrate specificity. Emiliania (Haptophyta), Anophryoides, Pseudocohnilembus, Vitrella and Chromera (Alveolata), and Monosiga (Choanomonada) all contained a gene for arginine kinase (AK). In contrast, Aphanomyces, Albugo and Ectocarpus (Stramenopiles), and Guillardia (Cryptophyta) possessed a gene for taurocyamine kinase (TK). The Guillardia TK enzyme exhibited rather strong substrate inhibition toward taurocyamine, which was analyzed using the most likely kinetic model. This was the first report of substrate inhibition in a TK. Together with the research results from other groups, the AK, TK, or creatine kinase (CK) activities have been observed sporadically in at least six groups of protists. However, it is not clear the three enzyme activities were emerged early in the evolution and divergence of protist groups, or some of enzyme activities were introduced to the protists by horizontal gene transfer. In addition, we found that seven protist enzymes examined in this study possess a myristoylation signaling sequence at the N-terminus. The amino-acid sequence around the guanidine-specificity region and the key residue at 89^th position of the protist AK and CK were homologous to those of the metazoan enzymes, but those for protist TKs were different indicating that the latter evolved independently.

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.

Mitochondrial Proteolipid Complexes of Creatine Kinase

Isoforms of creatine kinase (CK) generate and use phosphocreatine, a concentrated and highly diffusible cellular "high energy" intermediate, for the main purpose of energy buffering and transfer in order to maintain cellular energy homeostasis. The mitochondrial CK isoform (mtCK) localizes to the mitochondrial intermembrane and cristae space, where it assembles into peripherally membrane-bound, large cuboidal homooctamers. These are part of proteolipid complexes wherein mtCK directly interacts with cardiolipin and other anionic phospholipids, as well as with the VDAC channel in the outer membrane. This leads to a stabilization and cross-linking of inner and outer mitochondrial membrane, forming so-called contact sites. Also the adenine nucleotide translocator of the inner membrane can be recruited into these proteolipid complexes, probably mediated by cardiolipin. The complexes have functions mainly in energy transfer to the cytosol and stimulation of oxidative phosphorylation, but also in restraining formation of reactive oxygen species and apoptosis. In vitro evidence indicates a putative role of mtCK in mitochondrial phospholipid distribution, and most recently a role in thermogenesis has been proposed. This review summarizes the essential structural and functional data of these mtCK complexes and describes in more detail the more recent advances in phospholipid interaction, thermogenesis, cancer and evolution of mtCK.

Arginine kinase from Myzostoma cirriferum, a basal member of annelids

We assembled a phosphagen kinase gene from the Expressed Sequence Tags database of Myzostoma cirriferum, a basal member of annelids. The assembled gene sequence was synthesized using an overlap extension polymerase chain reaction method and was expressed in Escherichia coli. The recombinant enzyme (355 residues) exhibited monomeric behavior on a gel filtration column and showed strong activity only for l-arginine. Thus, the enzyme was identified as arginine kinase (AK). The two-substrate kinetic parameters were obtained and compared with other AKs. Phylogenetic analysis of amino acid sequences of phosphagen kinases indicated that the Myzostoma AK gene lineage differed from that of the polychaete Sabellastarte spectabilis AK, which is a dimer of creatine kinase (CK) origin. It is likely that the Myzostoma AK gene lineage was lost at an early stage of annelid evolution and that Sabellastarte AK evolved secondarily from the CK gene. This work contributes to our understanding of the evolution of phosphagen kinases of annelids with marked diversity.

Characterization of the arginine kinase isoforms in Caenorhabditis elegans

Phosphagen kinases (PKs) are well-studied enzymes involved in energy homeostasis in a wide range of animal, protozoan, and even some bacterial species. Recent genome efforts have allowed comparative work on the PKs to extend beyond the biochemistry of individual proteins to the comparative cellular physiology and examining of the role of all PK family members in an organism. The sequencing of the Caenorhabditis elegans genome and availability of sophisticated genetic tools within that system affords the opportunity to conduct a detailed physiological analysis of the PKs from a well known invertebrate for comparison with the extensive work conducted on vertebrate systems. As a first step in this effort we have carried out a detailed molecular genetic and biochemical characterization of the PKs in C. elegans. Our results reveal that C. elegans has five PK genes encoding arginine kinases that range in catalytic efficiency (kcat/KM(Arg)) from (3.1±0.6)×10(4) to (9±4)×10(5) M(-1) s(-1). This range is generally within the range seen for arginine kinases from a variety of species. Our molecular genetic and phylogenetic analysis reveals that the gene family has undergone extensive intron loss and gain within the suborder Rhabditina. In addition, within C. elegans we find evidence of gene duplication and loss. The analysis described here for the C. elegans AKs represents one of the most complete biochemical and molecular genetic analysis of a PK family within a genetically tractable invertebrate system and opens up the possibility of conducting detailed physiological comparisons with vertebrate systems using the sophisticated tools available with this model invertebrate system.

Phosphagen kinase in Schistosoma japonicum: II. Determination of amino acid residues essential for substrate catalysis using site-directed mutagenesis

Phosphagen kinases (PKs) play major roles in the regulation of energy metabolism in animals. Creatine kinase (CK) is the sole PK in vertebrates, whereas several PKs are present in invertebrates. We previously identified a contiguous dimer taurocyamine kinase (TK) from the trematode Schistosoma japonicum (Sj), a causative agent of schistosomiasis. SjTK contiguous dimer is comprised of domain 1 (D1) and domain 2 (D2). In this study, we used SjTK contiguous dimer (SjTKD1D2) or truncated single-domain constructs (SjTKD1 or SjTKD2) and employed site-directed mutagenesis to investigate the enzymatic properties of TK mutants. Mutation in SjTKD1 or SjTKD2 (D1E222G or D2E225G) caused complete loss of activity for the substrate taurocyamine. Likewise, a double mutant (D1E222GD2E225G) in the contiguous dimer (D1D2) exhibited complete loss of activity for the substrate taurocyamine. However, catalytic activity in the contiguous dimer remained in both of D1 inactive mutant (D1D2D1E222G) and D2 inactive mutant (D1D2D2E225G), suggesting that efficient catalysis of SjTKD1D2 is dependent on the activity of D1 and D2. The catalytic efficiency of the mixture of both single domains (WTD1+WTD2) showed same enzymatic properties (Km(Tauro)=0.68;Vmax/Km(Tauro)=137.04) to WTD1D2 (Km(Tauro)=0.47; Vmax/Km(Tauro)=144.30). This result suggests that the contiguous dimeric structure is not essential for the catalytic efficiencies of both domains of SjTK. Vmax/Km(Tauro) of the mixture of wild-type and inactivated domains (78.02 in WTD1+D2E225G and 128.24 in D1E222G+WTD2) were higher than the corresponding mutants (47.25 in D1D2D1E222G and 46.77 in D1D2D2E225G). To identify amino acid residues that are critical for taurocyamine binding, we performed alanine scanning mutagenesis at positions 57-63 on the guanidino specificity (GS) region of the SjTKD1, which is considered to be involved in guanidino-substrate recognition. R63A and R63Y mutants lost activity for taurocyamine, suggesting that these residues are associated with taurocyamine binding. In addition, we investigated the role of Tyr84 in D1 and found an association with substrate alignment. The Y84 residue was replaced with R, H, K, I, A, and G. Although the activities of each mutant were decreased (Vmax=2.36-67.50μmolPi/min/mgprotein), Y84 mutants possess binding affinity for taurocyamine (Km(Tauro)=3.19-10.04mM). The D1Y84R, D1Y84H, D1Y84K, and D1Y84A mutants exhibited low activity for taurocyamine, whereas the D1Y84I and D1Y84G mutants exhibited slightly decreased activity compared with the other Y84 mutants. The D1Y84K mutant lost substrate synergy between taurocyamine and ATP, suggesting that this mutation moves the position of the GS loop, similar to that of lombricine kinase (LK), and interferes with taurocyamine binding. This is the first comprehensive investigation of essential amino acid residues for substrate catalysis in trematode TK.

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.

Phosphagen kinase in Schistosoma japonicum: characterization of its enzymatic properties and determination of its gene structure

Phosphagen kinases (PKs) play a major role in the regulation of energy metabolism in animals. Creatine kinase (CK) is the sole PK in vertebrates, whereas several PKs are present in invertebrates. Here, we report the enzymatic properties and gene structure of PK in the trematode Schistosoma japonicum (Sj). SjPK has a unique contiguous dimeric structure comprising domain 1 (D1) and domain 2 (D2). The three states of the recombinant SjPK (D1, D2, and D1D2) show a specific activity for the substrate taurocyamine. The comparison of the two domains of SjPK revealed that D1 had a high turnover rate (kcat=52.91) and D2 exhibited a high affinity for taurocyamine (Km(Tauro) =0.53±0.06). The full-length protein exhibited higher affinity for taurocyamine (Km(Tauro) =0.47±0.03) than the truncated domains (D1=1.30±0.10, D2=0.53±0.06). D1D2 also exhibited higher catalytic efficiency (kcat/Km(Tauro) =82.98) than D1 (40.70) and D2 (29.04). These results demonstrated that both domains of SjTKD1D2 interacted efficiently and remained functional. The three-dimensional structure of SjPKD1 was constructed by the homology modeling based on the transition state analog complex state of Limulus AK. This protein model of SjPKD1 suggests that the overall structure is almost conserve between SjPKD1 and Limulus AK except for the flexible loops, that is, particularly guanidino-specificity (GS) region, which is associated with the recognition of the corresponding guanidino substrate. The constructed NJ tree and the comparison of exon/intron organization suggest that SjTK has evolved from an arginine kinase (AK) gene. SjTK has potential as a novel antihelminthic drug target as it is absent in mammals and its strong activity may imply a significant role for this protein in the energy metabolism of the parasite.

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.

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.

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.

Insights into shell deposition in the Antarctic bivalve Laternula elliptica: gene discovery in the mantle transcriptome using 454 pyrosequencing

Background

The Antarctic clam, Laternula elliptica, is an infaunal stenothermal bivalve mollusc with a circumpolar distribution. It plays a significant role in bentho-pelagic coupling and hence has been proposed as a sentinel species for climate change monitoring. Previous studies have shown that this mollusc displays a high level of plasticity with regard to shell deposition and damage repair against a background of genetic homogeneity. The Southern Ocean has amongst the lowest present-day CaCO₃ saturation rate of any ocean region, and is predicted to be among the first to become undersaturated under current ocean acidification scenarios. Hence, this species presents as an ideal candidate for studies into the processes of calcium regulation and shell deposition in our changing ocean environments.

Results

454 sequencing of L. elliptica mantle tissue generated 18,290 contigs with an average size of 535 bp (ranging between 142 bp-5.591 kb). BLAST sequence similarity searching assigned putative function to 17% of the data set, with a significant proportion of these transcripts being involved in binding and potentially of a secretory nature, as defined by GO molecular function and biological process classifications. These results indicated that the mantle is a transcriptionally active tissue which is actively proliferating. All transcripts were screened against an in-house database of genes shown to be involved in extracellular matrix formation and calcium homeostasis in metazoans. Putative identifications were made for a number of classical shell deposition genes, such as tyrosinase, carbonic anhydrase and metalloprotease 1, along with novel members of the family 2 G-Protein Coupled Receptors (GPCRs). A membrane transport protein (SEC61) was also characterised and this demonstrated the utility of the clam sequence data as a resource for examining cold adapted amino acid substitutions. The sequence data contained 46,235 microsatellites and 13,084 Single Nucleotide Polymorphisms(SNPs/INDELS), providing a resource for population and also gene function studies.

Conclusions

This is the first 454 data from an Antarctic marine invertebrate. Sequencing of mantle tissue from this non-model species has considerably increased resources for the investigation of the processes of shell deposition and repair in molluscs in a changing environment. A number of promising candidate genes were identified for functional analyses, which will be the subject of further investigation in this species and also used in model-hopping experiments in more tractable and economically important model aquaculture species, such as Crassostrea gigas and Mytilus edulis.

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.

cDNA identification, comparison and phylogenetic aspects of lombricine kinase from two oligochaete species

Creatine kinase and arginine kinase are the typical representatives of an eight-member phosphagen kinase family, which play important roles in the cellular energy metabolism of animals. The phylum Annelida underwent a series of evolutionary processes that resulted in rapid divergence and radiation of these enzymes, producing the greatest diversity of the phosphagen kinases within this phylum. Lombricine kinase (EC 2.7.3.5) is one of such enzymes and sequence information is rather limited compared to other phosphagen kinases. This study presents data on the cDNA sequences of lombricine kinase from two oligochaete species, the California blackworm (Lumbriculus variegatus) and the sludge worm (Tubifex tubifex). The deduced amino acid sequences are analyzed and compared with other selected phosphagen kinases, including two additional lombricine kinase sequences extracted from DNA databases and provide further insights in the evolution and position of these enzymes within the phosphagen kinase family. The data confirms the presence of a deleted region within the flexible loop (the GS region) of all six examined lombricine kinases. A phylogenetic analysis of these six lombricine kinases clearly positions the enzymes together in a small subcluster within the larger creatine kinase (EC 2.7.3.2) clade.

Cloning, expression, characterization and phylogenetic analysis of arginine kinase from greasyback shrimp (Metapenaeus ensis)

Arginine kinase (AK) plays an important role in cellular energy metabolism in invertebrate. The encoding AK gene from Shrimp Metapenaeus ensis (M. ensis) was cloned in prokaryotic expression plasmid pET-28a, and it was then expressed in Escherichia coil in dissoluble form. The recombinant protein was purified by following three chromatography steps in turn: CM-Cellulose cation-exchange, Sephacryl S-100HR gel filtrate and DEAE-Sepharose anion-exchange. The purified AK's apparent K(m) was 2.33+/-0.1 and 1.59+/-0.2 mM for ATP and l-arginine, respectively, while its optimum pH and temperature was 8.5 and 30 degrees C in the process of forward reaction, respectively. Phylogenetic analysis of cDNA-derived amino acid sequences for the AKs indicated a close affinity of M. ensis and another shrimp (Litopenaeus vannamei).

A novel arginine kinase from the shrimp Neocaridina denticulata: the fourth arginine kinase gene lineage

Arginine kinases (AK) evolved independently at least three times: first at an early stage of phosphagen kinase evolution (typical AK), second from the cytoplasmic creatine kinase (CK) gene later in metazoan evolution (Stichopus AK) and third from the mitochondrial CK gene in the course of annelid radiation (Sabellastarte AK). Here, we present a possible fourth AK lineage. We amplified cDNA encoding three AKs (AK1, AK2 and AK3) from the shrimp Neocaridina denticulata, and determined their amino acid sequences (355-356 residues). Each cDNA sequence was cloned in a pET plasmid and the corresponding recombinant kinase was expressed in E. coli. The AKs showed monomeric nature similar to typical AKs on an analytical gel filtration column. While the amino acid sequence of AK2 corresponded to that of typical AK, containing the conserved key residues established in Limulus AK for the substrate binding site, those of AK1 and AK3 lacked some of these key residues, indicating a similar evolution to Stichopus and Sabellastarte AKs. Moreover, phylogenetic analysis of phosphagen kinases indicated that Neocaridina AK1 and AK3 diverged at the deepest branching point close to the root of the tree and formed a distinct cluster outside the typical AK cluster, which included Neocaridina AK2. Kinetic constants of Neocaridina AKs were similar to those of other AKs. However, activation energy (E(a)) for the transition state of AK1 and AK3 was about 1.5-fold larger than that of AK2. The DeltaH(o++) values for AK1 and AK3 were also about 1.5-fold larger than that for AK2, but all three DeltaG(o++) values were very similar (71-72 kJ/mol); this results in similar reaction velocities for the three AK reactions due to a lower decrease in entropy in the activated complexes of AK1 and AK3 reactions compared to that of AK2. These findings suggest that Neocaridina AK1 and AK3 are unique from the known three AK gene lineages and therefore comprises a fourth AK lineage.

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.

Two-domain arginine kinase from the deep-sea clam Calyptogena kaikoi--evidence of two active domains

The cDNA and deduced amino acid sequences for arginine kinase (AK) from the deep-sea clam Calyptogena kaikoi have been determined revealing an unusual two-domain (2D) structure with molecular mass of 80 kDa, twice that of normal AK. The amino acid sequences of both domains contain most of the residues thought to be required for substrate binding found in the horseshoe crab Limulus polyphemus AK, a well studied system for which several X-ray crystal structures exist. However, two highly conserved residues, D62 and R193, that form a salt bridge thereby stabilizing the substrate-bound structure have been replaced by G and N in domain 1, and G and P in domain 2, respectively. The present effort probes whether both domains of Calyptogena AK are catalytically competent. Recombinant constructs of the wild-type enzyme of both single domains, and of selected mutants of the Calyptogena AK have been expressed as fusion proteins with the maltose-binding protein. The wild-type two-domain enzyme (2D[WT]) had high AK activity (k(cat)=23 s(- 1), average value of the two domains), and the single domain 2 (D2[WT]) showed 1.5-times higher activity (k(cat)=38 s(- 1)) than the wild-type 2D[WT]. Interestingly, the single domain 1 (D1[WT]) showed only a very low activity (k(cat) approximately 0.016 s(- 1)). Introduction of a Y68A mutation in both domains virtually abolished catalytic activity. On the other hand, significant residual activity was observed (k(cat)=2.8 s(- 1)), when the Y68A mutation was introduced only into domain 2 of the two-domain enzyme. A similar mutation in domain 1 of the two-domain enzyme reduced activity to a much lower extent (k(cat)=11.1 s(- 1)). Although the domains of this "contiguous" dimeric AK each have catalytic capabilities, the presence of domain 2 strongly influences the stability and activity of domain 1.

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.