Determination of the affinity of each component of a composite quaternary transition-state analogue complex of creatine kinase

Molecular docking and density functional theory studies on creatine, guanidinoacetic acid, and their phosphorylated analogues binding to muscle creatine kinase

Comparative molecular docking studies on creatine and guanidinoacetic acid, as well as their phosphorylated analogues, creatine phosphate, and phosphorylated guanidinoacetic acid, are investigated. Docking and density functional theory studies are carried out for muscle creatine kinase. The changes in the geometries of the ligands before and after binding to the enzyme are investigated to explain the better binding of guanidinoacetic acid and phosphorylated guanidinoacetic acid compared to creatine and creatine phosphate.

Elevated μs-ms timescale backbone dynamics in the transition state analog form of arginine kinase

Arginine kinase catalyzes reversible phosphoryl transfer between arginine and ATP. Crystal structures of arginine kinase in an open, substrate-free form and closed, transition state analog (TSA) complex indicate that the enzyme undergoes substantial domain and loop rearrangements required for substrate binding, catalysis, and product release. Nuclear magnetic resonance (NMR) has shown that substrate-free arginine kinase is rigid on the ps-ns timescale (average S²=0.84±0.08) yet quite dynamic on the µs-ms timescale (35 residues with R_ex, 12%), and that movements of the N-terminal domain and the loop comprising residues I182-G209 are rate-limiting on catalysis. Here, NMR of the TSA-bound enzyme shows similar rigidity on the ps-ns timescale (average S²=0.91±0.05) and substantially increased μs-ms timescale dynamics (77 residues; 22%). Many of the residues displaying μs-ms dynamics in NMR Carr-Purcell-Meiboom-Gill (CPMG) ¹⁵N backbone relaxation dispersion experiments of the TSA complex are also dynamic in substrate-free enzyme. However, the presence of additional dynamic residues in the TSA-bound form suggests that dynamics extend through much of the C-terminal domain, which indicates that in the closed form, a larger fraction of the protein takes part in conformational transitions to the excited state(s). Conformational exchange rate constants (k_ex) of the TSA complex are all approximately 2500s^-1, higher than any observed in the substrate-free enzyme (800-1900s^-1). Elevated μs-ms timescale protein dynamics in the TSA-bound enzyme is more consistent with recently postulated catalytic networks involving multiple interconnected states at each step of the reaction, rather than a classical single stabilized transition state.

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.

Selectivity of kinases on the activation of tenofovir, an anti-HIV agent

Nucleoside analogues, used in HIV-therapy, need to be phosphorylated by cellular enzymes in order to become potential substrates for HIV reverse transcriptase. After incorporation into the viral DNA chain, because of lacking of their 3'-hydroxyl groups, they stop the elongation process and lead to the death of the virus. Phosphorylation of the HIV-drug derivative, tenofovir monophosphate was tested with the recombinant mammalian nucleoside diphosphate kinase (NDPK), 3-phosphoglycerate kinase (PGK), creatine kinase (CK) and pyruvate kinase (PK). Among them, only CK was found to phosphorylate tenofovir monophosphate with a reasonable rate (about 45-fold lower than with its natural substrate, ADP), while PK exhibits even lower, but still detectable activity (about 1000-fold lower compared to the value with ADP). On the other hand, neither NDPK nor PGK has any detectable activity on tenofovir monophosphate. The absence of activity with PGK is surprising, since the drug tenofovir competitively inhibits both CK and PGK towards their nucleotide substrates, with similar inhibitory constants, K(I) of 2.9 and 4.8mM, respectively. Computer modelling (docking) of tenofovir mono- or diphosphate forms to these four kinases suggests that the requirement of large-scale domain closure for functioning (as for PGK) may largely restrict their applicability for phosphorylation/activation of pro-drugs having a structure similar to tenofovir monophosphate.

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.

A method for removing effects of nonspecific binding on the distribution of binding stoichiometries: application to mass spectroscopy data

There is often an interest in knowing, for a given ligand concentration, how many protein molecules have one, two, three, etc. ligands bound in a specific manner. This is a question that cannot be addressed using conventional ensemble techniques. Here, a mathematical method is presented for separating specific from nonspecific binding in nonensemble studies. The method provides a way to determine the distribution of specific binding stoichiometries at any ligand concentration when using nonensemble (e.g., single-molecule) methods. The applicability of the method is demonstrated for ADP binding to creatine kinase using mass spectroscopy data. A major advantage of our method, which can be applied to any protein-ligand system, is that no previous information regarding the mechanism of ligand interaction is required.

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.

Deconstructing the native state: energy landscapes, function, and dynamics of globular proteins

Proteins are highly complex molecules with features exquisitely selected by nature to carry out essential biological functions. Physical chemistry and polymer physics provide us with the tools needed to make sense of this complexity. Upon translation, many proteins fold to a thermodynamically stable form known as the native state. The native state is not static, but consists of a hierarchy of conformations, that are continuously explored through dynamics. In this review we provide a brief introduction to some of the core concepts required in the discussion of the protein native dynamics using energy landscapes ideas. We first discuss recent works which have challenged the structure-function paradigm by demonstrating function in disordered proteins. Next we examine the hierarchical organization in the energy landscapes using atomistic molecular dynamics simulations and principal component analysis. In particular, the role of direct and water-mediated contacts in sculpting the landscape is elaborated. Another approach to studying the native state ensemble is based on choosing high-resolution order parameters for computing one- or two-dimensional free energy surfaces. We demonstrate that 2D free energy surfaces provide rich thermodynamic and kinetic information about the native state ensemble. Brownian dynamics simulations on such a surface indicate that protein conformational dynamics is weakly activated. Finally, we briefly discuss implicit and coarse-grained protein models and emphasize the solvent role in determining native state structure and dynamics.

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.

Stereospecificity, substrate, and inhibitory properties of nucleoside diphosphate analogs for creatine and pyruvate kinases

Antiviral alpha-P-borano substituted NTPs are promising chain terminators targeting HIV reverse transcriptase (RT). Activation of antiviral nucleoside diphosphates (NDPs) to NTPs may be carried out by pyruvate kinase (PK) and creatine kinase (CK). Herein, are presented the effects of nucleobase, ribose, and alpha-phosphate substitutions on substrate specificities of CK and PK. Both enzymes showed two binding modes and negative cooperativity with respect to substrate binding. The stereospecificity and inhibition of ADP phosphorylation by alpha-P-borano substituted NDP (NDPalphaB) stereoisomers were also investigated. The Sp-ADPalphaB isomer was a 70-fold better substrate for CK than the Rp isomer, whereas PK preferred the Rp isomer of NDPalphaBs. For CK, the Sp-ADPalphaB isomer was a competitive inhibitor; for PK, the Rp-ADPalphaB isomer was a poor competitive inhibitor and the Sp-ADPalphaB isomer was a poor non-competitive inhibitor. Taken together, these data suggest that, although the Rp-NDPalphaB isomer would be minimally phosphorylated by CK or PK, it should not inhibit either enzyme.

Kinetic analysis of two purified forms of arginine kinase: absence of cooperativity in substrate binding of dimeric phosphagen kinase

Arginine kinase from sea urchin eggs and sea cucumber muscle are dimeric enzymes, unlike the more widely distributed monomeric enzyme found in other invertebrates. Both purified enzymes exhibited features characteristic of the monomeric arginine kinases including pH optima, formation of a catalytic dead-end complex (enzyme-MgADP-arginine) and stabilization of this complex by monovalent anions. A complete analysis of initial velocity data, in both directions for each substrate, indicated that substrate binding cooperativity was either minimal or non-existent. Unlike many other multi-subunit enzymes, the significance of the dimeric state of the phosphagen kinases remains unclear. These present results would suggest that (a) cooperativity, or so-called synergism in substrate binding is not a characteristic of the dimeric state of the protein and (b) the functional significance of the dimeric state is not related to the ability of some of these enzymes to undergo cooperativity in substrate binding. The significance of the dimeric state for the creatine kinases and arginine kinases remains to be established.

A deconvolution method for the separation of specific versus nonspecific interactions in noncovalent protein-ligand complexes analyzed by ESI-FT-ICR mass spectrometry

A method to separate specific and nonspecific noncovalent interactions observed in ESI mass spectra between a protein and its ligands is presented. Assuming noncooperative binding, the specific ligand binding is modeled as a statistical distribution on identical binding sites. For the nonspecific fraction we assume a statistical distribution on a large number of "nonspecific" interacting sites. The model was successfully applied to the noncovalent interaction between the protein creatine kinase (CK) and its ligands adenosine diphosphate (ADP) and adenosine triphosphate (ATP) that both exhibit nonspecific binding in the mass spectrum. The two sequential dissociation constants obtained by applying our method are K(1,diss) = 11.8 +/- 1.5 microM and K(2,diss) = 48 +/- 6 microM for ADP. For ATP, the constants are K(1,diss) = 27 +/- 7 microM and K(2,diss) = 114 +/- 27 microM. All constants are in good correlation with reported literature values. The model should be valuable for systems with a large dissociation constant that require high ligand concentrations and thus have increased potential of forming nonspecific adducts.

Synthesis of alpha-P-modified nucleoside diphosphates with ethylenediamine

This report describes a one-pot synthesis of alpha-P-borano-, alpha-P-thio-, and alpha-P-seleno-modified nucleoside diphosphate analogues that are otherwise difficult to obtain. The key step involves the intramolecular nucleophilic attack by an amino group in 5 to remove the gamma-phosphate. The absolute configurations of P-diastereomers were confirmed by analysis of their 1H NMR. Affinity studies revealed that the nucleoside boranodiphosphates are potentially useful in antiviral research.

Regulation of NDR protein kinase by hydrophobic motif phosphorylation mediated by the mammalian Ste20-like kinase MST3

NDR protein kinases are involved in the regulation of cell cycle progression and morphology. NDR1/NDR2 protein kinase is activated by phosphorylation on the activation loop phosphorylation site Ser281/Ser282 and the hydrophobic motif phosphorylation site Thr444/Thr442. Autophosphorylation of NDR is responsible for phosphorylation on Ser281/Ser282, whereas Thr444/Thr442 is targeted by an upstream kinase. Here we show that MST3, a mammalian Ste20-like protein kinase, is able to phosphorylate NDR protein kinase at Thr444/Thr442. In vitro, MST3 selectively phosphorylated Thr442 of NDR2, resulting in a 10-fold stimulation of NDR activity. MOB1A (Mps one binder 1A) protein further increased the activity, leading to a fully active kinase. In vivo, Thr442 phosphorylation after okadaic acid stimulation was potently inhibited by MST3KR, a kinase-dead mutant of MST3. Knockdown of MST3 using short hairpin constructs abolished Thr442 hydrophobic motif phosphorylation of NDR in HEK293F cells. We conclude that activation of NDR is a multistep process involving phosphorylation of the hydrophobic motif site Thr444/2 by MST3, autophosphorylation of Ser281/2, and binding of MOB1A.

Transition state stabilization by six arginines clustered in the active site of creatine kinase

Six fully conserved arginine residues (R129, R131, R235, R291, R319, and R340) closely grouped in the nucleotide binding site of rabbit muscle creatine kinase (rmCK) were mutated; four to alanine and all six to lysine. Kinetic analyses in the direction of phosphocreatine formation showed that all four alanine mutants led to substantial losses of activity with three (R129A, R131A, and R235A) having no detectable activity. All six lysine mutants retained variable degrees of reduced enzymatic activity. Static quenching of intrinsic tryptophan fluorescence was used to measure the binding constants for MgADP and MgATP. Nucleotide binding was at most only modestly affected by mutation of the arginine residues. Thus, the cluster of arginines seem to be primarily responsible for transition state stabilization which is further supported by the observation that none of the inactive mutants demonstrated the ability to form a transition analogue complex of MgADP.nitrate.creatine as determined by fluorescence quenching assays. As a whole, the results suggest that the most important role these residues play is to properly align the substrates for stabilization of the phosphoryl transfer reaction.

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.

The role of phosphagen specificity loops in arginine kinase

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

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.