Structure of human muscle creatine kinase

Depletion of creatine phosphagen energetics with a covalent creatine kinase inhibitor

Creatine kinases (CKs) provide local ATP production in periods of elevated energetic demand, such as during rapid anabolism and growth. Thus, creatine energetics has emerged as a major metabolic liability in many rapidly proliferating cancers. Whether CKs can be targeted therapeutically is unknown because no potent or selective CK inhibitors have been developed. Here we leverage an active site cysteine present in all CK isoforms to develop a selective covalent inhibitor of creatine phosphagen energetics, CKi. Using deep chemoproteomics, we discover that CKi selectively engages the active site cysteine of CKs in cells. A co-crystal structure of CKi with creatine kinase B indicates active site inhibition that prevents bidirectional phosphotransfer. In cells, CKi and its analogs rapidly and selectively deplete creatine phosphate, and drive toxicity selectively in CK-dependent acute myeloid leukemia. Finally, we use CKi to uncover an essential role for CKs in the regulation of proinflammatory cytokine production in macrophages.

GPU-Accelerated All-Atom Particle-Mesh Ewald Continuous Constant pH Molecular Dynamics in Amber

Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicit-solvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replica-exchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pK_a's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pK_a shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pK_a's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pK_a upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pK_a downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.

Naturally Occurring Mutations to Muscle-Type Creatine Kinase Impact Its Canonical and Pharmacological Activities in a Substrate-Dependent Manner In Vitro

Tenofovir (TFV) is a key component of human immunodeficiency virus (HIV) pre-exposure prophylaxis (PrEP). TFV is a nucleotide analog reverse-transcriptase inhibitor prodrug that requires two separate phosphorylation reactions by intracellular kinases to form the active metabolite tenofovir-diphosphate (TFV-DP). Muscle-type creatine kinase (CKM) has previously been demonstrated to be the kinase most responsible for the phosphorylation of tenofovir-monophosphate (TFV-MP) to the active metabolite in colon tissue. Because of the importance of CKM in TFV activation, genetic variation in CKM may contribute to interindividual variability in TFV-DP levels. In the present study, we report 10 naturally occurring CKM mutations that reduced TFV-MP phosphorylation in vitro: T35I, R43Q, I92M, H97Y, R130H, R132C, F169L, Y173C, W211R, V280L, and N286I. Interestingly, of these 10, only 4—R130H, R132C, W211R, and N286I—reduced both canonical CKM activities: ADP phosphorylation and ATP dephosphorylation. Although positions 130, 132, and 286 are located in the active site, the other mutations that resulted in decreased TFV-MP phosphorylation occur elsewhere in the protein structure. Four of these eight mutations—T35I, R43Q, I92M, and W211R—were found to decrease the thermal stability of the protein. Additionally, the W211R mutation was found to impact protein structure both locally and at a distance. These data suggest a substrate-specific effect such that certain mutations are tolerated for canonical activities while being deleterious toward the pharmacological activity of TFV activation, which could influence PrEP outcomes.

Alkylation of rabbit muscle creatine kinase surface methionine residues inhibits enzyme activity in vitro

Creatine kinase (CK) catalyzes the formation of phosphocreatine from adenosine triphosphate (ATP) and creatine. The highly reactive free cysteine residue in the active site of the enzyme (Cys²⁸³) is considered essential for the enzymatic activity. In previous studies we demonstrated that Cys²⁸³ is targeted by the alkylating chemical warfare agent sulfur mustard (SM) yielding a thioether with a hydroxyethylthioethyl (HETE)-moiety. In the present study, the effect of SM on rabbit muscle CK (rmCK) activity was investigated with special focus on the alkylation of Cys²⁸³ and of reactive methionine (Met) residues. For investigation of SM-alkylated amino acids in rmCK, micro liquid chromatography-electrospray ionization high-resolution tandem-mass spectrometry measurements were performed using the Orbitrap technology. The treatment of rmCK with SM resulted in a decrease of enzyme activity. However, this decrease did only weakly correlate to the modification of Cys²⁸³ but was conclusive for the formation of Met⁷⁰-HETE and Met¹⁷⁹-HETE. In contrast, the activity of mutants of rmCK produced by side-directed mutagenesis that contained substitutions of the respective Met residues (Met⁷⁰Ala, Met¹⁷⁹Leu, and Met⁷⁰Ala/Met¹⁷⁹Leu) was highly resistant against SM. Our results point to a critical role of the surface exposed Met⁷⁰ and Met¹⁷⁹ residues for CK activity.

Calmodulin complexes with brain and muscle creatine kinase peptides

Calmodulin (CaM) is a ubiquitous Ca2+ sensing protein that binds to and modulates numerous target proteins and enzymes during cellular signaling processes. A large number of CaM-target complexes have been identified and structurally characterized, revealing a wide diversity of CaM-binding modes. A newly identified target is creatine kinase (CK), a central enzyme in cellular energy homeostasis. This study reports two high-resolution X-ray structures, determined to 1.24 ​Å and 1.43 ​Å resolution, of calmodulin in complex with peptides from human brain and muscle CK, respectively. Both complexes adopt a rare extended binding mode with an observed stoichiometry of 1:2 CaM:peptide, confirmed by isothermal titration calorimetry, suggesting that each CaM domain independently binds one CK peptide in a Ca2+-depended manner. While the overall binding mode is similar between the structures with muscle or brain-type CK peptides, the most significant difference is the opposite binding orientation of the peptides in the N-terminal domain. This may extrapolate into distinct binding modes and regulation of the full-length CK isoforms. The structural insights gained in this study strengthen the link between cellular energy homeostasis and Ca2+-mediated cell signaling and may shed light on ways by which cells can 'fine tune' their energy levels to match the spatial and temporal demands.

Identification of creatine kinase and alpha-1 antitrypsin as protein targets of alkylation by sulfur mustard

Sulfur mustard (SM) is a toxic chemical warfare agent deployed in several conflicts within the last 100 years and still represents a threat in terroristic attacks and warfare. SM research focuses on understanding the pathophysiology of SM and identifying novel biomarkers of exposure. SM is known to alkylate nucleophilic moieties of endogenous proteins, for example, free thiol groups of cysteine residues. The two-dimensional-thiol-differences in gel electrophoresis (2D-thiol-DIGE) technique is an initial proteomics approach to detect proteins with free cysteine residues. These amino acids are selectively labeled with infrared-maleimide dyes visualized after GE. Cysteine residues derivatized by alkylating agents are no longer accessible for the maleimide-thiol coupling resulting in the loss of the fluorescent signal of the corresponding protein. To prove the applicability of 2D-thiol-DIGE, this technology was exemplarily applied to neat human serum albumin treated with SM, to lysates from human cell culture exposed to SM as well as to human plasma exposed to CEES (chloroethyl ethyl sulfide, an SM analogue). Exemplarily, the most prominent proteins modified by SM were identified by matrix-assisted laser desorption/ionization time-of-flight (tandem) mass spectrometry, MALDI-TOF MS(/MS), as creatine kinase (CK) from human cells and as alpha-1 antitrypsin (A1AT) from plasma samples. Peptides containing the residue Cys²⁸² of CK and Cys²³² of A1AT were unambiguously identified by micro liquid chromatography-electrospray ionization high-resolution tandem-mass spectrometry (μLC-ESI MS/HR MS) as being alkylated by SM bearing the specific hydroxyethylthioethyl-(HETE)-moiety. Both peptides might represent potential biomarkers of SM exposure. This is the first report introducing these endogenous proteins as targets of SM alkylation.

Molecular docking investigation of the amantadine binding to the enzymes upregulated or downregulated in Parkinson's disease

Parkinson's disease (PD) is a progressive neurodegenerative disease. Levodopa in combination with amantadine has a demonstrated efficacy in motility impairment. An extensive investigation of some enzymes described to be upregulated or downregulated in PD was made - adenylate kinase (AK), adenine phosphoribosyltransferase (APRT), ectonucleoside triphosphate diphosphohydrolase 1 (ENTPD1), nucleoside-diphosphate kinase 3 (NDK3), purine nucleoside phosphorylase 1 (PNP1), and ecto-5'-nucleotidase (NT5E). Also, creatine kinase (CK) was included in the study because it is one of the main enzymes involved in the regulation of the nucleotide ratio in energy metabolism. To date, there is no proven link between amantadine treatment of PD and these enzymes. Because there are many AKs isoforms modified in PD, the AK was the first investigated. The molecular docking experiments allow the analysis of the selective binding of amantadine - unionized (with -NH₂ group) and ionized form (with -NH₃ ⁺ group) - to the AKs' isoforms implicated in PD. Using available X-ray 3D structures of human AKs in closed-conformation, it was demonstrated that there are notable differences between the interactions of the two forms of amantadine for the zebrafish AK1 (5XZ2), human AK2 (2C9Y), human AK5 (2BWJ), and AK from B.stearothermophilus. The cytosolic human AK1 and human AK2 mostly interact with ionized amantadine by AMP binding residues. The human AK5 interaction with ionized amantadine does not involve the residues from the catalytic site. Among other enzymes tested in the present study, APRT revealed the best results in respect of binding amantadine ionized form. The results offer a new perspective for further investigation of the connections between amantadine treatment of PD and some enzymes involved in purine metabolism.

Predicting Reactive Cysteines with Implicit-Solvent-Based Continuous Constant pH Molecular Dynamics in Amber

Cysteines existing in the deprotonated thiolate form or having a tendency to become deprotonated are important players in enzymatic and cellular redox functions and frequently exploited in covalent drug design; however, most computational studies assume cysteines as protonated. Thus, developing an efficient tool that can make accurate and reliable predictions of cysteine protonation states is timely needed. We recently implemented a generalized Born (GB) based continuous constant pH molecular dynamics (CpHMD) method in Amber for protein pKa calculations on CPUs and GPUs. Here we benchmark the performance of GB-CpHMD for predictions of cysteine pKa's and reactivities using a data set of 24 proteins with both down- and upshifted cysteine pKa's. We found that 10 ns single-pH or 4 ns replica-exchange CpHMD titrations gave root-mean-square errors of 1.2-1.3 and correlation coefficients of 0.8-0.9 with respect to experiment. The accuracy of predicting thiolates or reactive cysteines at physiological pH with single-pH titrations is 86 or 81% with a precision of 100 or 90%, respectively. This performance well surpasses the traditional structure-based methods, particularly a widely used empirical pKa tool that gives an accuracy less than 50%. We discuss simulation convergence, dependence on starting structures, common determinants of the pKa downshifts and upshifts, and the origin of the discrepancies from the structure-based calculations. Our work suggests that CpHMD titrations can be performed on a desktop computer equipped with a single GPU card to predict cysteine protonation states for a variety of applications, from understanding biological functions to covalent drug design.

Insights into the Phosphoryl Transfer Mechanism of Human Ubiquitous Mitochondrial Creatine Kinase

Human ubiquitous mitochondrial creatine kinase (uMtCK) is responsible for the regulation of cellular energy metabolism. To investigate the phosphoryl-transfer mechanism catalyzed by human uMtCK, in this work, molecular dynamic simulations of uMtCK∙ATP-Mg2+∙creatine complex and quantum mechanism calculations were performed to make clear the puzzle. The theoretical studies hereof revealed that human uMtCK utilizes a two-step dissociative mechanism, in which the E227 residue of uMtCK acts as the catalytic base to accept the creatine guanidinium proton. This catalytic role of E227 was further confirmed by our assay on the phosphatase activity. Moreover, the roles of active site residues in phosphoryl transfer reaction were also identified by site directed mutagenesis. This study reveals the structural basis of biochemical activity of uMtCK and gets insights into its phosphoryl transfer mechanism.

Model of the Ankyrin and SOCS Box Protein, ASB9, E3 Ligase Reveals a Mechanism for Dynamic Ubiquitin Transfer

Cullin-RING E3 ligases (CRLs) are elongated and bowed protein complexes that transfer ubiquitin over 60 Å to proteins targeted for proteasome degradation. One such CRL contains the ankyrin repeat and SOCS box protein 9 (ASB9), which binds to and partially inhibits creatine kinase (CK). While current models for the ASB9-CK complex contain some known interface residues, the overall structure and precise interface of the ASB9-CK complex remains unknown. Through an integrative modeling approach, we report a third-generation model that reveals precisely the interface interactions and also fits the shape of the ASB9-CK complex as determined by small-angle X-ray scattering. We constructed an atomic model for the entire CK-targeting CRL to uncover dominant modes of motion that could permit ubiquitin transfer. Remarkably, only the correctly docked CK-containing E3 ligase and not incorrectly docked structures permitted close approach of ubiquitin to the CK substrate.

Correcting of pronated feet reduce skeletal muscle injury in young women with biomechanical abnormalities

Biomechanical abnormalities of pronated feet accompanied by functional leg length disparity may increase the risk of skeletal muscle injury. Objective of the study is to prove that correction of pronated feet by the foot orthoses will reduce the creatine kinase-MM (CK-MM) concentrations as the muscle injury indicator. The design study was double blind randomized clinical trials with control. Research subjects were divided into two groups, group 1 used the foot orthoses while group 2 did not used the foot orthoses. The whole subject examined the concentrations of the CK-MM enzyme before, and 24–72 hours after the walking test. The walking test was conducted 15 minutes with maximum speed. The concentration of the CK-MM enzyme before walking test on treatment group was 70.07±15.33 International Unit (IU), similar with the control group was 69.85±17.03 IU (P=0.971). The increased in CK-MM enzyme concentrations 45 hours after the walking test was lower in the treatment group (7.8±9 IU) than the control group (22.0±11.5 IU) (P=0.001). The CK-MM enzyme concentrations continued to decline in the treatment group after the second walking test (77.21±17.47 IU), and after the third walking test (69.86±11.88 IU) (P=0.018). The foot orthoses for correcting the pronated feet on the young women with biomechanical abnormalities is able to reduce the degree of the skeletal muscle injury after walking activity.

A single residue substitution accounts for the significant difference in thermostability between two isoforms of human cytosolic creatine kinase

Creatine kinase (CK) helps maintain homeostasis of intracellular ATP level by catalyzing the reversible phosphotransfer between ATP and phosphocreatine. In humans, there are two cytosolic CK isoforms, the muscle-type (M) and the brain-type (B), which frequently function as homodimers (hMMCK and hBBCK). Interestingly, these isoenzymes exhibit significantly different thermostabilities, despite high similarity in amino acid sequences and tertiary structures. In order to investigate the mechanism of this phenomenon, in this work, we first used domain swapping and site-directed mutagenesis to search for the key residues responsible for the isoenzyme-specific thermostability. Strikingly, the difference in thermostability was found to principally arise from one single residue substitution at position 36 (Pro in hBBCK vs. Leu in hMMCK). We then engaged the molecular dynamics simulations to study the molecular mechanism. The calculations imply that the P36L substitution introduces additional local interactions around residue 36 and thus further stabilizes the dimer interface through a complex interaction network, which rationalizes the observation that hMMCK is more resistant to thermal inactivation than hBBCK. We finally confirmed this molecular explanation through thermal inactivation assays on Asp36 mutants that were proposed to devastate the local interactions and thus the dimer associations in both isoenzymes.

The substrate-free and -bound crystal structures of the duplicated taurocyamine kinase from the human parasite Schistosoma mansoni

The taurocyamine kinase from the blood fluke Schistosoma mansoni (SmTK) belongs to the phosphagen kinase (PK) family and catalyzes the reversible Mg(2+)-dependent transfer of a phosphoryl group between ATP and taurocyamine. SmTK is derived from gene duplication, as are all known trematode TKs. Our crystallographic study of SmTK reveals the first atomic structure of both a TK and a PK with a bilobal structure. The two unliganded lobes present a canonical open conformation and interact via their respective C- and N-terminal domains at a helix-mediated interface. This spatial arrangement differs from that observed in true dimeric PKs, in which both N-terminal domains make contact. Our structures of SmTK complexed with taurocyamine or l-arginine compounds explain the mechanism by which an arginine residue of the phosphagen specificity loop is crucial for substrate specificity. An SmTK crystal was soaked with the dead end transition state analog (TSA) components taurocyamine-NO3 (2-)-MgADP. One SmTK monomer was observed with two bound TSAs and an asymmetric conformation, with the first lobe semiclosed and the second closed. However, isothermal titration calorimetry and enzyme kinetics experiments showed that the two lobes function independently. A small angle x-ray scattering model of SmTK-TSA in solution with two closed active sites was generated.

Effects of the Fc-III tag on activity and stability of green fluorescent protein and human muscle creatine kinase

The Fc-III tag is a newly developed fusion tag that can be applied to protein purification and detection. In the present work, we use the Fc-III-tagged green fluorescent protein (GFP) and human muscle creatine kinase (CK) as model systems to investigate effects of the Fc-III tag on activities and stabilities of the expressed multicysteine-containing proteins. Our results show the Fc-III tag has no adverse effects on the fluorescence of GFP and reduces the occurrence of GFP misfolding due to incorrect Cys oxidation compared with the His-tagged protein. The activity and stability of the Fc-III-tagged CK is slightly lower than that of the tag-free CK, but is higher than that of the His-tagged CK as determined by the ratio of the oxidized versus reduced CK. A major portion of His-tagged CK is in its oxidized form, while that of the Fc-III-tagged CK is in its reduced form. A folding model of CK with different tags was proposed, which may provide insights into the effect of the Fc-III tag on the conformations of disulfide-bridged proteins.

A study of the mechanism of the chaperone-like function of an scFv of human creatine kinase by computer simulation

A new application of antibodies is to use them as macromolecular chaperones. Protein antigens usually have multiple epitopes, thus, there may be a plurality of antibodies binding to one antigen. However, not all antibodies that bind to one antigen could act as a chaperone. Experiments show that some screened anti-human creatine kinase single chain antibodies (scFV) could assist in the folding and stabilizing of the enzyme, while others could not. We built the model of the single chain antibody (scFv-A4) that increased the stability of human creatine kinase (HCK) by the homology modeling method. Epitopes of human creatine kinase were predicted by computer and then the binding of scFv-A4 and HCK was modeled with computer. The calculation results were further combined with the peptide array membrane experiment results to obtain reliable models for the scFv-A4-HCK complex. Based on the above study we gave an explanation about how scFv-A4 could act as a macromolecular chaperone assisting the folding of HCK. This study provides an approach for predicting antigen-antibody binding mode and also a useful theoretical guidance for the study of antibodies' chaperone-like function.

Dynamical properties of the loop 320s of substrate-free and substrate-bound muscle creatine kinase by NMR: evidence for independent subunits

Muscle creatine kinase (MCK; EC2.7.3.2) is a 86 kDa homodimer that belongs to the family of guanidino kinases. MCK has been intensively studied for several decades, but it is still not known why it is a dimer because this quaternary structure does not translate into obvious structural or functional advantages over the homologous monomeric arginine kinase. In particular, it remains to be demonstrated whether MCK subunits are independent. Here, we describe NMR chemical-shift perturbation and relaxation experiments designed to study the active site 320s flexible loop of this enzyme. The analysis was performed with the enzyme in its ligand-free and MgADP-complexed forms, as well as with the transition-state analogue abortive complex (MCK-Mg-ADP-creatine-nitrate ion). Our data indicate that each subunit can bind substrates independently.

Lestaurtinib inhibits histone phosphorylation and androgen-dependent gene expression in prostate cancer cells

Background

Epigenetics is defined as heritable changes in gene expression that are not based on changes in the DNA sequence. Posttranslational modification of histone proteins is a major mechanism of epigenetic regulation. The kinase PRK1 (protein kinase C related kinase 1, also known as PKN1) phosphorylates histone H3 at threonine 11 and is involved in the regulation of androgen receptor signalling. Thus, it has been identified as a novel drug target but little is known about PRK1 inhibitors and consequences of its inhibition.

Methodology/Principal Finding

Using a focused library screening approach, we identified the clinical candidate lestaurtinib (also known as CEP-701) as a new inhibitor of PRK1. Based on a generated 3D model of the PRK1 kinase using the homolog PKC-theta (protein kinase c theta) protein as a template, the key interaction of lestaurtinib with PRK1 was analyzed by means of molecular docking studies. Furthermore, the effects on histone H3 threonine phosphorylation and androgen-dependent gene expression was evaluated in prostate cancer cells.

Conclusions/Significance

Lestaurtinib inhibits PRK1 very potently in vitro and in vivo. Applied to cell culture it inhibits histone H3 threonine phosphorylation and androgen-dependent gene expression, a feature that has not been known yet. Thus our findings have implication both for understanding of the clinical activity of lestaurtinib as well as for future PRK1 inhibitors.

Dissimilarity in the folding of human cytosolic creatine kinase isoenzymes

Creatine kinase (CK, EC 2.7.3.2) plays a key role in the energy homeostasis of excitable cells. The cytosolic human CK isoenzymes exist as homodimers (HMCK and HBCK) or a heterodimer (MBCK) formed by the muscle CK subunit (M) and/or brain CK subunit (B) with highly conserved three-dimensional structures composed of a small N-terminal domain (NTD) and a large C-terminal domain (CTD). The isoforms of CK provide a novel system to investigate the sequence/structural determinants of multimeric/multidomain protein folding. In this research, the role of NTD and CTD as well as the domain interactions in CK folding was investigated by comparing the equilibrium and kinetic folding parameters of HMCK, HBCK, MBCK and two domain-swapped chimeric forms (BnMc and MnBc). Spectroscopic results indicated that the five proteins had distinct structural features depending on the domain organizations. MBCK BnMc had the smallest CD signals and the lowest stability against guanidine chloride-induced denaturation. During the biphasic kinetic refolding, three proteins (HMCK, BnMc and MnBc), which contained either the NTD or CTD of the M subunit and similar microenvironments of the Trp fluorophores, refolded about 10-fold faster than HBCK for both the fast and slow phase. The fast folding of these three proteins led to an accumulation of the aggregation-prone intermediate and slowed down the reactivation rate thereby during the kinetic refolding. Our results suggested that the intra- and inter-subunit domain interactions modified the behavior of kinetic refolding. The alternation of domain interactions based on isoenzymes also provides a valuable strategy to improve the properties of multidomain enzymes in biotechnology.

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

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

Regulation of sodium-calcium exchanger activity by creatine kinase under energy-compromised conditions

Na(+)/Ca(2+) exchanger (NCX) is one of the major mechanisms for removing Ca(2+) from the cytosol especially in cardiac myocytes and neurons, where their physiological activities are triggered by an influx of Ca(2+). NCX contains a large intracellular loop (NCXIL) that is responsible for regulating NCX activity. Recent evidence has shown that proteins, including kinases and phosphatases, associate with NCX1IL to form a NCX1 macromolecular complex. To search for the molecules that interact with NCX1IL and regulate NCX1 activity, we used the yeast two-hybrid method to screen a human heart cDNA library and found that the C-terminal region of sarcomeric mitochondrial creatine kinase (sMiCK) interacted with NCX1IL. Moreover, both sMiCK and the muscle-type creatine kinase (CKM) coimmunoprecipitated with NCX1 using lysates of cardiacmyocytes and HEK293T cells that transiently expressed NCX1 and various creatine kinases. Both sMiCK and CKM were able to produce a recovery in the decreased NCX1 activity that was lost under energy-compromised conditions. This regulation is mediated through a putative PKC phosphorylation site of sMiCK and CKM. The autophosphorylation and the catalytic activity of sMiCK and CKM are not required for their regulation of NCX1 activity. Our results suggest a novel mechanism for the regulation of NCX1 activity.

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.

Structural studies of human brain-type creatine kinase complexed with the ADP-Mg2+-NO3- -creatine transition-state analogue complex

Creatine kinase is a member of the phosphagen kinase family, which catalyzes the reversible phosphoryl transfer reaction that occurs between ATP and creatine to produce ADP and phosphocreatine. Here, three structural aspects of human-brain-type-creatine-kinase (hBB-CK) were identified by X-ray crystallography: the ligand-free-form at 2.2A; the ADP-Mg2+, nitrate, and creatine complex (transition-state-analogue complex; TSAC); and the ADP-Mg2+-complex at 2.0A. The structures of ligand-bound hBB-CK revealed two different monomeric states in a single homodimer. One monomer is a closed form, either bound to TSAC or the ADP-Mg2+-complex, and the second monomer is an unliganded open form. These structural studies provide a detailed mechanism indicating that the binding of ADP-Mg2+ alone may trigger conformational changes in hBB-CK that were not observed with muscle-type-CK.

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.

Detection of local polarity and conformational changes at the active site of rabbit muscle creatine kinase with a new arginine-specific fluorescent probe

A new polarity-sensitive fluorescent probe, 3-(4-chloro-6-p-glyoxal-phenoxy-1,3,5-triazinylamino)-7-(dimethylamino)-2-methylphenazine (CGTDP), is synthesized for selective labeling of active-site arginine residues. The probe comprises a neutral red moiety as a polarity-sensitive fluorophore and a phenylglyoxal unit as an arginine-specific labeling group. The probe exhibits a sensitive response of shift of fluorescence maximum emission wavelength to solvent polarity only instead of pH or temperature, which leads to the use of the probe in detecting the local polarity and conformational changes of the active site of rabbit muscle creatine kinase (CK) denatured by pH or temperature. The polarity of the active site domain has been first found to correspond to a dielectric constant of about 44, and the conformational change of the active site directly revealed by CGTDP occurs far before that of CK as a whole disclosed by the intrinsic tryptophan fluorescence during acid or thermal denaturation. The present strategy may provide a useful method to detect the local polarity and conformational changes of the active sites of many enzymes that employ arginine residues as anion recognition sites under different denaturation conditions.

Effects of the single point genetic mutation D54G on muscle creatine kinase activity, structure and stability

Aberrant folding of important proteins caused by genetic mutations is closely correlated to many diseases. Due to the important physiological role in excitable cells, the activity and level of creatine kinase (CK) play a crucial role in maintaining body functions. Muscle CK deficiency disease was identified by an unusual CK activity decrease in an acute myocardial infarction patient caused by the single point mutation D54G. In this research, it was found that the D54G mutant had substantially decreased activity, substrate binding affinity and stability. Spectroscopic experiments indicated that the mutation impaired the structure of CK, which resulted in a partially unfolded state with more hydrophobic exposure and exposed Trp residues. The inability to fold to the functional compact state made the mutant be prone to aggregate upon microenvironmental stresses, and might gradually decrease the CK level of the patient.

Very fast empirical prediction and rationalization of protein pKa values

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.

Conformational change in the C-terminal domain is responsible for the initiation of creatine kinase thermal aggregation

Protein conformational changes may be associated with particular properties such as its function, transportation, assembly, tendency to aggregate, and potential cytotoxicity. In this research, the conformational change that is responsible for the fast destabilization and aggregation of rabbit muscle creatine kinase (EC 2.7.3.2) induced by heat was studied by intrinsic fluorescence and infrared spectroscopy. A pretransitional change of the tryptophan microenvironments was found from the intrinsic fluorescence spectra. A further analysis of the infrared spectra using quantitative second-derivative and two-dimensional correlation analysis indicated that the changes of the beta-sheet structures in the C-terminal domain and the loops occurred before the formation of intermolecular cross-beta-sheet structures and the unfolding of alpha-helices. These results suggested that the pretransitional conformational changes in the active site and the C-terminal domain might result in the modification of the domain-domain interactions and the formation of an inactive dimeric form that was prone to aggregate. Our results highlighted the fact that some minor conformational changes, which were usually negligible or undetectable by normal methods, might play a crucial role in protein stability and aggregation. Our results also suggested that the changes in domain-domain interactions, but not the dissociation of the dimer, might play a crucial role in the thermal denaturation and aggregation of this dimeric two-domain protein.

Loop movement and catalysis in creatine kinase

Recently the crystal structure of creatine kinase from Torpedocalifornica was determined to 2.1 A. The dimeric structure revealed two different forms in the unit cell: one monomer was bound to a substrate, MgADP, and the other monomer was bound to a transition-state analogue complex composed of MgADP, nitrate and creatine. The most striking difference between the structures is the movement of two loops (comprising residues 60-70 and residues 323-333) into the active site in the transition state structure. This loop movement effectively occludes the active site from solvent, and the loops appear to be locked into place by a salt bridge formed between His66 and Asp326. His66 is of particular interest as it is located within a PGHP motif conserved in all creatine kinases but not found in other guanidino kinases. We have carried out alanine-scanning mutagenesis of each of the residues in the PGHP motif and determined that only the His66 plays a significant role in the creatine kinase reaction. Although neither residue interacts directly with the substrate, the interaction His66 and Asp326 appears to be important in providing the precise alignment of substrates necessary for phosphoryl group transfer. Finally, it is clear that neither His66 nor Asp326 are responsible for the pKs observed in the pH-rate profile for HMCK.

Relating structure to mechanism in creatine kinase

Found in all vertebrates, creatine kinase catalyzes the reversible reaction of creatine and ATP forming phosphocreatine and ADP. Phosphocreatine may be viewed as a reservoir of "high-energy phosphate" which is able to supply ATP, the primary energy source in bioenergetics, on demand. Consequently, creatine kinase plays a significant role in energy homeostasis of cells with intermittently high energy requirements. The enzyme is of clinical importance and its levels are routinely used as an indicator of myocardial and skeletal muscle disorders and for the diagnosis of acute myocardial infarction. First identified in 1928, the enzyme has undergone intensive investigation for over 75 years. There are four major isozymes, two cytosolic and two mitochondrial, which form dimers and octamers, respectively. Depending on the pH, the enzyme operates by a random or an ordered bimolecular mechanism, with the equilibrium lying towards phosphocreatine production. Evidence suggests that conversion of creatine to phosphocreatine occurs via the in-line transfer of a phosphoryl group from ATP. A recent X-ray structure of creatine kinase bound to a transition state analog complex confirmed many of the predictions based on kinetic, spectroscopic, and mutagenesis studies. This review summarizes and correlates the more significant mechanistic and structural studies on creatine kinase.

Local dynamics measured by hydrogen/deuterium exchange and mass spectrometry of creatine kinase digested by two proteases

Hydrogen/deuterium exchange coupled to mass spectrometry has been used to investigate the structure and dynamics of native dimeric cytosolic muscle creatine kinase. The protein was incubated in D2O for various time. After H/D exchange and rapid quenching of the reaction, the partially deuterated protein was cleaved in parallel by two different proteases (pepsin or type XIII protease from Aspergillus saitoi) to increase the sequence coverage and spatial resolution of deuterium incorporation. The resulting peptides were analyzed by liquid chromatography coupled to mass spectrometry. In comparison with the 3D structure of MM-CK, the analysis of the two independent proteolysis deuteration patterns allowed us to get new insights into CK local dynamics as compared to a previous study using pepsin [Mazon et al. Protein Science 13 (2004) 476-486]. In particular, we obtained more information on the kinetics and extent of deuterium exchange in the N- and C-terminal extremities represented by the 1-22 and 362-380 pepsin peptides. Indeed, we observed a very different behaviour of the 1-12 and 13-22 type XIII protease peptides, and similarly for the 362-373 and 374-380 peptides. Moreover, comparison of the deuteration patterns of type XIII protease segments of the large 90-126 pepsin peptide led us to identify a small relatively dynamic region (108-114).

Hydrogen/deuterium exchange studies of native rabbit MM-CK dynamics

Creatine kinase (CK) isoenzymes catalyse the reversible transfer of a phosphoryl group from ATP onto creatine. This reaction plays a very important role in the regulation of intracellular ATP concentrations in excitable tissues. CK isoenzymes are highly resistant to proteases in native conditions. To appreciate localized backbone dynamics, kinetics of amide hydrogen exchange with deuterium was measured by pulse-labeling the dimeric cytosolic muscle CK isoenzyme. Upon exchange, the protein was digested with pepsin, and the deuterium content of the resulting peptides was determined by liquid chromatography coupled to mass spectrometry (MS). The deuteration kinetics of 47 peptides identified by MS/MS and covering 96% of the CK backbone were analyzed. Four deuteration patterns have been recognized: The less deuterated peptides are located in the saddle-shaped core of CK, whereas most of the highly deuterated peptides are close to the surface and located around the entrance to the active site. Their exchange kinetics are discussed by comparison with the known secondary and tertiary structures of CK with the goal to reveal the conformational dynamics of the protein. Some of the observed dynamic motions may be linked to the conformational changes associated with substrate binding and catalytic mechanism.

Determinants of cysteine pKa values in creatine kinase and α1-antitrypsin

The structural determinants of the unusually low pK(a) values of Cys282 in human creatine kinase and Cys232 in alpha1-antitrypsin were studied computationally. We have demonstrated that hydrogen bonding to the cysteine residue is the prime determinant for both proteins. In the case of creatine kinase, the hydrogen bond donors are a serine side chain and an amide NH-group, while in alpha1-antitrypsin the donor is an amide NH. Each hydrogen bond lowers the pK(a) by between 0.8 and 1.5 pH units. The 1.1-unit lowering due to the Ser284-Cys282 hydrogen bond is in good agreement with the 1.2-unit difference between the Cys282 pK(a) value of wild-type and the S284A mutant of creatine kinase.

Asparagine 285 plays a key role in transition state stabilization in rabbit muscle creatine kinase

To explore the possibility that asparagine 285 plays a key role in transition state stabilization in phosphagen kinase catalysis, the N285Q, N285D, and N285A site-directed mutants of recombinant rabbit muscle creatine kinase (rmCK) were prepared and characterized. Kinetic analysis of phosphocreatine formation showed that the catalytic efficiency of each N285 mutant was reduced by approximately four orders of magnitude, with the major cause of activity loss being a reduction in k(cat) in comparison to the recombinant native CK. The data for N285Q still fit a random-order, rapid-equilibrium mechanism, with either MgATP or creatine binding first with affinities very nearly equal to those for native CK. However, the affinity for the binding of the second substrate is reduced approximately 10-fold, suggesting that addition of a single methylene group at position 285 disrupts the symphony of substrate binding. The data for the N285A mutant only fit an ordered binding mechanism, with MgATP binding first. Isosteric replacement to form the N285D mutant has almost no effect on the K(M) values for either creatine or MgATP, thus the decrease in activity is due almost entirely to a 5000-fold reduction in k(cat). Using the quenching of the intrinsic CK tryptophan fluorescence by added MgADP (Borders et al. 2002), it was found that, unlike native CK, none of the mutants have the ability to form a quaternary TSAC. We use these data to propose that asparagine 285 indeed plays a key role in transition state stabilization in the reaction catalyzed by creatine kinase and other phosphagen kinases.

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

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