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. SIGNIFICANCE STATEMENT: Muscle-type creatine kinase (CKM) is important to the activation of tenofovir, a key component of HIV prophylaxis. This study demonstrates that naturally occurring CKM mutations impact enzyme function in a substrate-dependent manner such that some mutations that do not reduce canonical activities lead to reductions in the pharmacologically relevant activity. This finding at the intersection of drug metabolism and energy metabolism is important to the perspective on pharmacology of other drugs acted on by atypical drug-metabolizing enzymes.

Acetylation of muscle creatine kinase negatively impacts high-energy phosphotransfer in heart failure.

A hallmark of impaired myocardial energetics in failing hearts is the downregulation of the creatine kinase (CK) system. In heart failure patients and animal models, myocardial phosphocreatine content and the flux of the CK reaction are negatively correlated with the outcome of heart failure. While decreased CK activity is highly reproducible in failing hearts, the underlying mechanisms remains elusive. Here, we report an inverse relationship between the activity and acetylation of CK muscle form (CKM) in human and mouse failing hearts. Hyperacetylation of recombinant CKM disrupted MM homodimer formation and reduced enzymatic activity, which could be reversed by sirtuin 2 treatment. Mass spectrometry analysis identified multiple lysine residues on the MM dimer interface, which were hyperacetylated in the failing hearts. Molecular modeling of CK MM homodimer suggested that hyperacetylation prevented dimer formation through interfering salt bridges within and between the 2 monomers. Deacetylation by sirtuin 2 reduced acetylation of the critical lysine residues, improved dimer formation, and restored CKM activity from failing heart tissue. These findings reveal a potentially novel mechanism in the regulation of CK activity and provide a potential target for improving high-energy phosphoryl transfer in heart failure.

Structure-activity relationship analysis of dammarane-type natural products as muscle-type creatine kinase activators.

Muscle-type creatine kinase (CK-MM) is the target protein of ginsenosides in skeletal muscle. 20(S)-protopanaxadiol [20(S)-PPD] is an activator of CK-MM and exerts an anti-fatigue effect. In this study, twelve dammarane-type compounds were used for structure-activity relationship analysis in terms of enzyme activity, intermolecular interaction, and molecular docking. Enzyme activity analysis showed that 20(S)-PPD, 20(R)-PPD, 20(S)-protopanaxatriol [20(S)-PPT], 25-OH-PPD, 24-COOH-PPD, panaxadiol (PD), and ginsenoside Rh2 significantly increased CK-MM activity. Panaxatriol (PT), ocotillol, ginsenoside Rg1, and ginsenoside Rd had no significant influence on CK-MM activity, while jujubogenin inhibited its activity. Biolayer Interferometry (BLI) assay produced the same results as those on enzyme activity. The interaction intensity between dammarane-type compounds and CK-MM was linearly related to the compounds' maximum increment rate of enzyme activity. Molecular docking showed the following sequence of docking scores: Rd > Rg1 > Rh2 > 24-COOH-PPD > 20(S)-PPD > 20(S)-PPT > 25-OH-PPD > 20(R)-PPD > ocotillol > PT > PD > jujubogenin. We demonstrated that 20(S)-PPD was the best activator of CK-MM among the 12 dammarane-type compounds. The cyclization of the dammarane side chain, the hydroxyl group at position C6, and the glycosylation of C3, C6, and C20 reduced the ability to activate CK-MM. These findings can help in the development of enhanced CK-MM activators through structural modification.

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.

Expression of SERPINA3s in cattle: focus on bovSERPINA3-7 reveals specific involvement in skeletal muscle.

α1-Antichymotrypsin is encoded by the unique SERPINA3 gene in humans, while it is encoded by a cluster of eight closely related genes in cattle. BovSERPINA3 proteins present a high degree of similarity and significant divergences in the reactive centre loop (RCL) domains which are responsible for the antiprotease activity. In this study, we analysed their expression patterns in a range of cattle tissues. Even if their expression is ubiquitous, we showed that the expression levels of each serpin vary in different tissues of 15-month-old Charolais bulls. Our results led us to focus on bovSERPINA3-7, one of the two most divergent members of the bovSERPINA3 family. Expression analyses showed that bovSERPINA3-7 protein presents different tissue-specific patterns with diverse degrees of N-glycosylation. Using a specific antibody raised against bovSERPINA3-7, Western blot analysis revealed a specific 96 kDa band in skeletal muscle. BovSERPINA3-7 immunoprecipitation and mass spectrometry revealed that this 96 kDa band corresponds to a complex of bovSERPINA3-7 and creatine kinase M-type. Finally, we reported that the bovSERPINA3-7 protein is present in slow-twitch skeletal myofibres. Precisely, bovSERPINA3-7 specifically colocalized with myomesin at the M-band region of sarcomeres where it could interact with other components such as creatine kinase M-type. This study opens new prospects on the bovSERPINA3-7 function in skeletal muscle and promotes opportunities for further understanding of the physiological role(s) of serpins.

Development of a Clickable Probe for Profiling of Protein Glutathionylation in the Central Cellular Metabolism of E. coli and Drosophila.

Protein glutathionylation is an important post-translational modification that regulates many cellular processes, including energy metabolism, signal transduction, and protein homeostasis. Global profiling of glutathionylated proteins (denoted as glutathionylome) is crucial for understanding redox-regulated signal transduction. Here, we developed a novel method based on click reaction and proteomics to enrich and identify the glutathionylated peptides in Escherichia coli and Drosophila lysates, in which 937 and 1,930 potential glutathionylated peptides were identified, respectively. Bioinformatics analysis showed that the cysteine residue next to negatively charged amino acid residues has a higher frequency of glutathionylation. Importantly, we found that most proteins associated with metabolic pathways were glutathionylated and that the glutathionylation sites of metabolic enzymes were highly conserved among different species. Our results indicate that the glutathione analog is a useful tool to characterize protein glutathionylation, and glutathionylation of metabolic enzymes, which play important roles in regulating cellular metabolism, is conserved.

Dynamic asymmetry and the role of the conserved active-site thiol in rabbit muscle creatine kinase.

Symmetric and asymmetric crystal structures of the apo and transition state analogue forms, respectively, of the dimeric rabbit muscle creatine kinase have invoked an "induced fit" explanation for asymmetry between the two subunits and their active sites. However, previously reported thiol reactivity studies at the dual active-site cysteine 283 residues suggest a more latent asymmetry between the two subunits. The role of that highly conserved active-site cysteine has also not been clearly determined. In this work, the S-H vibrations of Cys283 were observed in the unmodified MM isoform enzyme via Raman scattering, and then one and both Cys283 residues in the same dimeric enzyme were modified to covalently attach a cyano group that reports on the active-site environment via its infrared CN stretching absorption band while maintaining the catalytic activity of the enzyme. Unmodified and Cys283-modified enzymes were investigated in the apo and transition state analogue forms of the enzyme. The narrow and invariant S-H vibrational bands report a homogeneous environment for the unmodified active-site cysteines, indicating that their thiols are hydrogen bonded to the same H-bond acceptor in the presence and absence of the substrate. The S-H peak persists at all physiologically relevant pH's, indicating that Cys283 is protonated at all pH's relevant to enzymatic activity. Molecular dynamics simulations identify the S-H hydrogen bond acceptor as a single, long-resident water molecule and suggest that the role of the conserved yet catalytically unnecessary thiol may be to dynamically rigidify that part of the active site through specific H-bonding to water. The asymmetric and broad CN stretching bands from the CN-modified Cys283 suggest an asymmetric structure in the apo form of the enzyme in which there is a dynamic exchange between spectral subpopulations associated with water-exposed and water-excluded probe environments. Molecular dynamics simulations indicate a homogeneous orientation of the SCN probe group in the active site and thus rule out a local conformational explanation at the residue level for the multipopulation CN stretching bands. The homogeneous simulated SCN orientation suggests strongly that a more global asymmetry between the two subunits is the cause of the CN probe's broad and asymmetric infrared line shape. Together, these spectral observations localized at the active-site cysteines indicate an intrinsic, dynamic asymmetry between the two subunits that exists already in the apo form of the dimeric creatine kinase enzyme, rather than being induced by the substrate. Biochemical and methodological consequences of these conclusions are considered.

The activity of carp muscle-specific creatine kinase at low temperature is enhanced by decreased hydrophobicity of residue 268.

Abstract The muscle-specific forms of creatine kinase in rabbit (RM-CK) and carp (M1-CK) exhibit different temperature-dependent functional properties. Replacing the glycine at residue 268 of RM-CK with asparagine increases the enzyme's activity at 10°C. In this study, we investigated how hydrophobicity of residue 268 affects the biochemical properties of RM-CK and M1-CK at low temperature. We generated three mutants of both RM-CK and M1-CK: Asp268, Lys268, and Leu268. The secondary structures of these mutants were similar, as revealed by their circular dichroism spectra. Similar to the Asn268 mutants, the Asp268 and Lys268 mutants of RM-CK and M1-CK exhibited higher specific activities at 10°C and pH 8.0. However, no such effect was observed for the RM-CK and M1-CK Leu268 mutants. While in the presence of cryoprotectant (sucrose or trehalose), the activities of wild-type RM-CK and M1-CK mutant enzymes with a hydrophobic residue at 268 were higher, and the effect was more profound at pH 8.0. It may be inferred that water molecules affect protein conformation around residue 268, thereby influencing protein stability at low temperature.

Effects of osmolytes on Pelodiscus sinensis creatine kinase: a study on thermal denaturation and aggregation.

The protective effect of osmolytes on the thermal denaturation and aggregation of Pelodiscus sinensis muscle creatine kinase (PSCK) was investigated by a combination of spectroscopic methods and thermodynamic analysis. Our results demonstrated that the addition of osmolytes, such as glycine and proline, could prevent thermal denaturation and aggregation of PSCK in a concentration-dependent manner. When the concentration of glycine and proline increased in the denatured system, the relative activation was significantly enhanced; meanwhile, the aggregation of PSCK during thermal denaturation was decreased. Spectrofluorometer results showed that glycine and proline significantly decreased the tertiary structural changes of PSCK and that thermal denaturation directly induced PSCK aggregation. In addition, we also built the 3D structure of PSCK and osmolytes by homology models. The results of computational docking simulations showed that the docking energy was relatively low and that the clustering groups were spread to the surface of PSCK, indicating that osmolytes directly protect the surface of the protein. P. sinensis are poikilothermic and quite sensitive to the change of ambient temperature; however, there were few studies on the thermal denaturation of reptile CK. Our study provides important insight into the protective effects of osmolytes on thermal denaturation and aggregation of PSCK.

Chaperone-like effects of a scFv antibody on the folding of human muscle creatine kinase.

Molecular chaperones play an essential role in assisting the folding of a myriad of nascent peptides to form different biologically active proteins. Therefore, their low substrate specificity is important for the functions of these housekeeping proteins. However, discovering chaperones which assist the folding of a particular protein can shed new light on the folding pathway of the protein, offering an interesting approach for developing specific therapeutic agents to treat protein-misfolding diseases. Screening of antibodies with chaperone-like function represents a novel strategy to meet the challenges. In this study, some single-chain variable fragment (scFv) antibodies were selected from a high-capacity phage antibody library using human muscle creatine kinase (HCK) as antigen. A scFv antibody (scFv A4) was determined to inhibit aggregation and favor recovery of the native conformation of HCK during its refolding. This antibody also increased the stability of HCK during its heat-induced unfolding process. Our findings demonstrate that scFv A4 has dual-chaperone-like activities: assisting in correct protein folding as well as protecting the native protein from unfolding. A molecular mechanism by which scFv A4 exhibits chaperone-like effects on HCK was proposed. This study demonstrates that phage antibody libraries can lead to chaperone-like proteins, and the specificity of the resulting antibody toward its antigen could provide new molecular details regarding how the chaperone interacts with the protein's unfolding and folding pathways.

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.

Disrupting of E79 and K138 interaction is responsible for human muscle creatine kinase deficiency diseases.

Creatine kinase (CK) is a key enzyme for cellular energy metabolism, catalyzing the reversible phosphoryl transfer from phosphocreatine to ADP in vertebrates. Due to its important physiological functions, the acquired CK somatic mutations are closely correlated to diseases. In this study, the E79G point mutation was identified in two acute myocardial infarction patients with muscle CK activity deficiency. The crystal structure of CK indicates that the E79 and K138 interaction plays key roles in sustaining the recognition between N-terminal and C-terminal domains. Mutations of these residues caused pronounced loss of activity, conformational changes and distinct substrate synergism alteration. Moreover, spectroscopic spectra experiments suggested that mutations disrupting this hydrogen bond impaired the secondary and tertiary structure of CK. Meanwhile, protein folding experiments implied that mutations lead them to the partially unfolded state which made them easier to be inactivated and unfolded under environmental stresses. Furthermore, this partially unfolded state upon environmental stresses might gradually decrease the CK level in the patients. Thus, these results might provide clues in the mechanism of CK deficiency diseases.

Chaperone-like effect of the linker on the isolated C-terminal domain of rabbit muscle creatine kinase.

Intramolecular chaperones (IMCs), which are specific domains/segments encoded in the primary structure of proteins, exhibit chaperone-like activity against the aggregation of the other domains in the same molecule. In this research, we found that the truncation of the linker greatly promoted the thermal aggregation of the isolated C-terminal domain (CTD) of rabbit muscle creatine kinase (RMCK). Either the existence of the linker covalently linked to CTD or the supply of the synthetic linker peptide additionally could successfully protect the CTD of RMCK against aggregation in a concentration-dependent manner. Truncated fragments of the linker also behaved as a chaperone-like effect with lower efficiency, revealing the importance of its C-terminal half in the IMC function of the linker. The aggregation sites in the CTD of RMCK were identified by molecular dynamics simulations. Mutational analysis of the three key hydrophobic residues resulted in opposing effects on the thermal aggregation between the CTD with intact or partial linker, confirming the role of linker as a lid to protect the hydrophobic residues against exposure to solvent. These observations suggested that the linkers in multidomain proteins could act as IMCs to facilitate the correct folding of the aggregation-prone domains. Furthermore, the intactness of the IMC linker after proteolysis modulates the production of off-pathway aggregates, which may be important to the onset of some diseases caused by the toxic effects of aggregated proteolytic fragments.

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.

Folding studies on muscle type of creatine kinase from Pelodiscus sinensis.

A folding study of creatine kinase from Pelodiscus sinensis has not yet been reported. To gain more insight into structural and folding mechanisms of P. sinensis CK (PSCK), denaturants such as SDS, guanidine HCl, and urea were applied in this study. We purified PSCK from the muscle of P. sinensis and conducted inhibition kinetics with structural unfolding studies under various conditions. The results revealed that PSCK was completely inactivated at 1.8 mM SDS, 1.05 M guanidine HCl, and 7.5 M urea. The kinetics via time-interval measurements showed that the inactivation by SDS, guanidine HCl, and urea were all first-order reactions with kinetic processes shifting from monophase to biphase at increasing concentrations. With respect to tertiary structural changes, PSCK was unfolded in different ways; SDS increased the hydrophobicity but retained the most tertiary structural conformation, while guanidine HCl and urea induced conspicuous changes in tertiary structures and initiated kinetic unfolding mechanisms. Our study provides information regarding PSCK and enhances our knowledge of the reptile-derived enzyme folding.

Activity and function of rabbit muscle-specific creatine kinase at low temperature by mutation at gly268 to asn268.

Carp muscle-specific creatine kinase M1 isoenzyme (M1-CK) seems to have evolved to adapt to synchronized changes in body temperature and intracellular pH. When gly(268) in rabbit muscle-specific creatine kinase was substituted with asn(268) as found in carp M1-CK, the rabbit muscle-specific CK G286N mutant specific activity at pH 8.0 and 10°C was more than 2-fold higher than that in the wild-type rabbit enzyme. Kinetic studies showed that K(m) values of the rabbit CK G268N mutant were similar to those of the wild-type rabbit enzyme, yet circular dichroism spectra showed that the overall secondary structures of the mutant enzyme, at pH 8.0 and 5°C, were almost identical to the carp M1-CK enzyme. The X-ray diffraction pattern of the mutant enzyme crystal revealed that amino acid residues involved in substrate binding are closer to one another than in the rabbit enzyme, and the cysteine283 active site of the mutant enzyme points away from the ADP binding site. At pH 7.4-8.0 and 35-10°C, with a smaller substrate, dADP, specific activities of the mutant enzyme were consistently higher than the wild-type rabbit enzyme and more similar to the carp M1-CK enzyme. Thus, the smaller active site of the RM-CK G268N mutant may be one of the reasons for its improved activity at low temperature.

Characterization of caged compounds binding to proteins by NMR spectroscopy.

Photolysable caged ligands are used to investigate protein function and activity. Here, we investigate the binding properties of caged nucleotides and their photo released products to well established but evolutionary and structurally unrelated nucleotide-binding proteins, rabbit muscle creatine kinase (RMCK) and human annexin A6 (hAnxA6), using saturation transfer difference NMR spectroscopy. We detect the binding of the caged nucleotides and discuss the general implications on interpreting data collected with photolysable caged ligands using different techniques. Strategies to avoid non-specific binding of caged compound to certain proteins are also suggested.

Effect of mixed crowding on refolding of human muscle creatine kinase.

The effects of mixed crowding agents containing both sucrose and dextran 70 on refolding process of human muscle creatine kinase (HCK) were studied by enzyme activity assay and aggregation measurements. The results showed that sucrose and dextran have opposite effect on parameters of HCK during refolding: reactivation yield, refolding rates and amount of aggregation, as they were both used in the mixed crowding agents. The exclusion volume effect of dextran and osmophobic effect of sucrose on HCK refolding can be counteracted by each other: sucrose bated the aggregation induced by dextran and increased the final reactivation yield and refolding rate of the slow track, while dextran inhibited the effect of sucrose to prevent aggregation and help correct folding. The effects of human cyclophilin 18 (hCyp18) and casein on folding of HCK were also studied in crowding conditions, and it was found that the chaperone function of hCyp18 was additive with sucrose but blocked by dextran, and that the aggregation-improving effect of casein was additive with dextran 70 but impaired by sucrose to certain extent. This study indicates osmolytes and macromolecule crowding agents play different roles in the physiological conditions and could lead to better understanding of protein folding in the intracellular environment.

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.

Dissecting the key residues crucial for the species-specific thermostability of muscle-type creatine kinase.

Species-specific protein thermal stability is closely correlated to the living conditions of the organism, especially to its body temperature. In this research, human and zebrafish muscle-type creatine kinases (MMCKs) were taken as model proteins to investigate the molecular adaptation of proteins in poikilothermal and homoiothermal animals. Both the optimal temperature for catalysis and the thermal stability of human MMCK was much higher than those of zebrafish MMCK. Sequence alignment identified 9 amino acid variations conserved in either the teleost MMCKs or the mammal and electric ray MMCKs. Bidirectional mutations were performed to find the residues with beneficial mutations. The results showed that two residues close to the dimer interface of MMCK, the 46th and 146th residue, were crucial for species-specific thermal stability.