Next Article in Journal
Salicylic Acid Is Required for Broad-Spectrum Disease Resistance in Rice
Previous Article in Journal
The Transmembrane Protease TMPRSS2 as a Therapeutic Target for COVID-19 Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The K+-Dependent and -Independent Pyruvate Kinases Acquire the Active Conformation by Different Mechanisms

by
Leticia Ramírez-Silva
1,*,
Gloria Hernández-Alcántara
1,
Carlos Guerrero-Mendiola
2,
Martin González-Andrade
1,
Adela Rodríguez-Romero
3,
Annia Rodríguez-Hernández
3,
Alan Lugo-Munguía
1,
Paul A. Gómez-Coronado
1,
Cristina Rodríguez-Méndez
1 and
Alicia Vega-Segura
1
1
Departamento de Bioquímica, Facultad de Medicina, Apartado Postal 70-159, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
2
Departamento de Genética Molecular, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
3
LANEM, Instituto de Química, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(3), 1347; https://doi.org/10.3390/ijms23031347
Submission received: 10 December 2021 / Revised: 15 January 2022 / Accepted: 17 January 2022 / Published: 25 January 2022
(This article belongs to the Section Biochemistry)

Abstract

:
Eukarya pyruvate kinases possess glutamate at position 117 (numbering of rabbit muscle enzyme), whereas bacteria have either glutamate or lysine. Those with E117 are K+-dependent, whereas those with K117 are K+-independent. In a phylogenetic tree, 80% of the sequences with E117 are occupied by T113/K114/T120 and 77% of those with K117 possess L113/Q114/(L,I,V)120. This work aims to understand these residues’ contribution to the K+-independent pyruvate kinases using the K+-dependent rabbit muscle enzyme. Residues 117 and 120 are crucial in the differences between the K+-dependent and -independent mutants. K+-independent activity increased with L113 and Q114 to K117, but L120 induced structural differences that inactivated the enzyme. T120 appears to be key in folding the protein and closure of the lid of the active site to acquire its active conformation in the K+-dependent enzymes. E117K mutant was K+-independent and the enzyme acquired the active conformation by a different mechanism. In the K+-independent apoenzyme of Mycobacterium tuberculosis, K72 (K117) flips out of the active site; in the holoenzyme, K72 faces toward the active site bridging the substrates through water molecules. The results provide evidence that two different mechanisms have evolved for the catalysis of this reaction.

1. Introduction

Rabbit muscle pyruvate kinase (WT-RMPK) was the first demonstration of an enzyme that requires monovalent cations for catalysis [1]. WT-RMPK exhibits an activity of 250 and 0.02 μmoles min−1mg−1 with and without 90 mM K+ [2], which can be partially substituted with Rb+, NH4+, and Tl+, albeit with activity 20–40% lower. In the presence of Cs+, Na+ and Li+, the enzyme exhibits only 9%, 8%, and 2%, respectively, of the maximal activity and hardly any activity (0.02%) with tris-(hydroxymethyl) aminomethane (Tris) [3] or (<0.001%) with substituted amines [4].
For many years, it was assumed that all pyruvate kinases (PKs) had an absolute requirement for K+, until 1969 when Benziman reported that the PK from Acetobacter xylinium did not require monovalent cations for activity [5]. After this first report, other examples of K+-independent enzymes arose [6,7,8,9,10,11,12,13,14]. In this respect, Laughlin and Reed, 1997, compared the sequences of the PKs from E. coli Type II, C. glutamicum and WT-RMPK [15]. The authors found that E117 (numbered according to WT-RMPK residues), close to the K+ binding site, has been substituted for Lys in the microbial enzymes. The proximity of the side chain from residue 117 to the binding site of K+ in WT-RMPK suggested that in the monovalent cation-independent enzymes, ε-NH3+ of Lys provides the inner positive charge that substitutes the activating effect of K+. Therefore, the authors constructed the mutant E117K and found it was K+-independent with ~12% of the maximal activity of the wild type enzyme with K+. Furthermore, it appeared that the activity was not stimulated by monovalent cations [15].
An extensive phylogenetic study of the PK family was conducted to establish how frequently E117 is replaced by K117 (according to the numbering of WT-RMPK) [16]. This analysis showed two clearly separated clusters of roughly the same number of sequences, one with E117 and the other with K117. All characterized enzymes with E117 exhibit K+-dependent activity, whereas those with K117 are K+-independent. Another relevant feature consists in the apparent co-evolution between the residue present at position 117 and those at positions 113, 114, and 120: 80% of the sequences that have E117 are occupied by T113/K114/T120, and 77% of those that have K117 possess L113/Q114/ (L, I, V)120. As shown in Figure 1, residues 114–118 form a short loop that connects the A and B domains. Residues 113 and 120 lie at the beginning and end of the loop, respectively. The active site sits in a cleft formed by these two domains [17].
This work focuses on understanding the individual and combined roles of highly conserved residues forming part of a ’signature’ sequence of K+-independent PKs. To this end, the residues at these positions in the K+-dependent sequence of WT-RMPK (T113/K114/E117/T120) were gradually replaced by those from K+-independent PKs (L113/Q114/K117/L120). It was expected that the quadruple mutant T113L/K114Q/E117K/T120L would confer K+-independence and would exhibit activities like those of WT-RMPK with K+. However, results indicate that K+-independent activity increases with the addition of residues L113 and Q114 to K117, obtaining half of the expected activity; but the addition of residue L120 leads to a nearly inactive enzyme. The secondary structure of mutants possessing T120L exhibited anomalous circular dichroism (CD) ellipticities compared to those of WT-RMPK and the other mutants, indicating different foldings. T120 and L120 play different roles; with the former related to the folding and closure of the lid in the K+-dependent enzymes. On the other hand, K117 or K72 (as found in Mycobacterium tuberculosis (MtbPK)) not only conferred K+-independence to PKs, but appeared key in the acquisition of the active conformation of the K+-independent PKs as demonstrated by the comparison of crystallographic data and molecular dynamic (MD) analysis of WT-RMPK (K+-dependent) and those of the MtbPK (K+-independent) [18]. Overall data provide convincing evidence that two different mechanisms have evolved for the catalysis of this reaction.

2. Results and Discussion

2.1. Effect of Li+, Na+, K+, NH4+, Rb+, Cs+, (CH3)4N+ and (CH3)3N+-CH2-CH2-OH (Choline) on the Activities of WT-RMPK and RMPK Mutants

According to previous studies [15], it was expected that mutants of RMPK with Lys in position 117 would exhibit K+-independent activity. To explore the activation of these enzymes by monovalent cations, specific effects first need to be distinguished from non-specific ionic strength effects. To this end, the effect of monovalent cations such as Li+, Na+, K+, NH4+, Rb+, Cs+, and that of bulky “inert” monovalent cations like (CH3)4NC+ and choline were studied on desalted WT-RMPK, E117K, T113L/E117K, and T113L/K114Q/E117K mutants. K114Q/E117K/T120L mutant was completely inactive and T120L and T113L/K114Q/E117K/T120L mutants were nearly inactive. Therefore, they were not included in the next experiments. (CH3)4N+ and choline are commonly used to maintain ionic strength [19] because either they do not interact with RMPK or they interact at non-specific sites on the protein [4]. An inhibitory effect due to high ionic strength (200 mM) has also been observed in RMPK with either (CH3)4N+ or activating monovalent cations [4]. Consequently, the effect of these bulky monovalent cations ((CH3)4N+ or choline) on WT-RMPK observed in Figure 2A was expected (less than 0.01% of maximal activity). Maximal activities were shown with 90 mM K+ (control) followed by Rb+ > NH4+ > Cs+ > Na+ > Li+. These activities also decreased with increasing ionic strength over 200 mM. In contrast, E117K (Figure 2B), T113L/E117K (Figure 2C), and T113L/K114Q/E117K mutants (Figure 2D) were 3, 7, and 2.7-fold activated by ~250 mM choline or (CH3)4N+, respectively. Furthermore, activating monovalent cations increased the activities of these three RMPK mutants.
The effects described herein were general for tetraalkylammonium salts, and monovalent cations and activities increased with ionic strength, except with Li+. Therefore, the effect over E117K, T113L/E117K, and T113L/K114Q/E117K mutants was due to a generalized ionic strength effect. Next, it was investigated whether the ionic strength effect of bulky and activating monovalent cations was additive. As expected, WT-RMPK (Figure 3A) expressed high activities with 90 mM K+ (white column) versus no activity without K+ and up to 100 mM (CH3)4N+ (white striped columns, not seen due to the scale). Activities were observed with 25 mM K+ and increased slightly with increasing K+ concentration but remained largely the same with different concentrations of (CH3)4N+. Hence, no additive effect was observed between K+ and bulky monovalent cation. In contrast, K+-independent activities (white columns) of E117K (Figure 3B), T113L/E117K (Figure 3C), and T113L/K114Q/E117K mutants (Figure 3D) were gradually raised by either the addition of an increasing amount of (CH3)4N+ (white striped columns) or of the bulky cation plus K+ or Rb+ (gray colored striped columns). That is to say that, the activating effect of bulky and monovalent cations were additive. Maximal activities of mutant PKs were achieved with 250 mM ionic strength, indicated with asterisks. Consequently, the next experiments were all performed in the presence of (CH3)4N+ to maintain constant ionic strength at 250 mM.
In previous studies of the E117K mutant [15,20], this effect was not observed because both used (CH3)4N+ to maintain ionic strength constant. Still, they did not explore the effect of increasing ionic strength on the activity of this mutant. Following the phylogenetic analysis [16], K+ binds to a highly conserved site. In the K+-independent PKs, three of the four residues that coordinate K+ are conserved in most enzymes. The only residue missing is T113 (numbering according to RMPK) which is Leu in K+-independent PKs. Hence, it is not surprising that these PKs still bind a monovalent cation as observed in the structure of MtbPK [18].

2.2. Effect of Activating Monovalent Cations to WT-RMPK and RMPK Mutants in the Presence of 250 mM Ionic Strength Maintained with (CH3)4N+

The activation of WT-RMPK, K114Q, E117K, T113L/E117K and T113L/K114Q/E117K mutants by different concentrations of Li+, Na+, K+, NH4+, Rb+ and Cs+ were measured and compared in the presence of saturating concentrations of PEP3−, Mg-ADP complex and Mg2+free and constant ionic strength of 250 mM (Figure 4). In confirmation of previous results [3], it was observed that maximal WT-RMPK activation followed the order K+ > NH4+ > Rb+ > Cs+ > Na+ > Li+. It is important to note that the activation curves of WT-RMPK with Na+ and Li+ did saturate (see Table 1). In previous studies, WT-RMPK activation exhibited no saturation of these cations at an ionic strength of 223 mM ionic strength maintained with HEPES [21]. (CH3)4N+ may bind to non-specific sites of the enzyme favoring Na+ and Li+ occupation of the monovalent binding site. As expected K114Q mutant was K+-dependent. Ion selectivity was preserved and kinetic constants for the monovalent cations were as those of wild-type PK (Table 1). As reported for RMPK [22,23], it is expected that K114Q follow Type Ib activation by monovalent cations. In such a mechanism, M+ coordination is absolutely required for catalysis, where both kcat and kcat/Km increase hyperbolically with K+. In contrast, ion selectivity of T113L was altered; maximal activation followed the order Rb+ > NH4+ > K+ > Cs+. It was found that the polarity of residue 113 is determinant in the partition of K+ into its site [21]. E117K (Figure 4B), T113L/E117K (Figure 4C), and T113L/K114Q/E117K mutants (Figure 4D) were K+-independent. In the presence of 250 mM ionic strength maintained with (CH3)4N+, no activation was further found for these mutants in the presence of the monovalent cations, only inhibition; i.e., their activities with 90 mM Li+ was 40% of that obtained without this cation.

2.3. Kinetics of WT-RMPK and RMPK Mutants

The K+-dependent PKs signature (T113/K114/E117/T120) was gradually substituted for that of the K+-independent PKs (L113/Q114/K117/L120). It was expected that the quadruple mutant T113L/K114Q/E117K/T120L would be K+-independent and would express a similar activity to that of WT-RMPK with K+. In this respect, it was observed that those mutants containing the substitution of Thr 120 for Leu became completely or nearly completely inactive after the last purification step. Therefore, the kinetic constants of WT-RMPK and all the mutants are shown in Table 2 except those containing the mutation T120L. Single mutants with E117 showed values of kcat similar to that of WT-RMPK; whereas the gradual change of the K+-independent signature to WT-RMPK (E117K, T113L/E117K, T113L/K114Q/E117K) followed a gradual increase in kcat from 35% to ~half catalytic constant of WT-RMPK with K+. In spite of the differences in catalytic constants, affinities for PEP3− and Mg-ADP for mutants with E117 or K117 were similar to those of WT-RMPK with K+; whereas mutants with K117 exhibited Km for Mg2+free from ~10 to 50-fold greater than that of WT-RMPK (KMg2+free 0.17 ± 0.01).
It was consistently observed that K+-independent mutants exhibited very low affinities for Mg2+free compared to K+-dependent PKs in this study and elsewhere [24]. This result may be explained by the presence of water molecules through which Mg-ATP coordinates to Mg2+ and oxalate in the MtbPK structure [18], instead of direct coordination between substrates as in WT-RMPK structure [17]. The lower binding energy of Mg2+ could result from the stabilizing contribution to maintaining water molecules between the divalent cation and the other substrates.
These results indicate that K117 is responsible for the RMPK´s K+-independent activity. The addition of residues L113 and Q114 present in the K+-independent sequences, promotes an increase in the K+-independent activity and the dependence of high concentrations of Mg2+. Therefore, residues L113, Q114, and K117 can allow K+-independent activity in the RMPK, but the inclusion of L120 seems to impair the enzyme.

2.4. Circular Dichroism Spectra of WT-RMPK and RMPK Mutants

To explore the effect on WT-RMPK structure of including residues observed in K+-independent PKs, we obtained the far UV circular dichroism (CD) spectra of WT-RMPK and all mutants. CD is widely used to determine whether an expressed recombinant purified protein is folded or if a mutation affects its conformation or stability [25]. Figure 5 shows that the spectra of WT-RMPK and K114Q almost overlapped; also a small loss of CD signal of ~15% to 20% was observed in E117K, T113L/E117K and T113L/K114Q/E117K mutants and a significative loss of CD signal of ~60% to 75% was observed for T120L, K114Q/E117K/T120L and T113L/K114Q/E 117K/T120L mutants.
While most mutants are native-like, CD data also confirmed the absence of integrity of those containing the substitution of Thr 120 for Leu, which coincides with their nearly null activities. This result indicates that T120 and L120 play different roles. T120 is involved in the folding or stability of the protein; when it was substituted for Leu, less than 40% of the secondary structure and ~1% of the wild type enzyme activity remained. Therefore, this finding suggests that K+-dependent and K+-independent PKs have different folding or stability requirements.

2.5. Structural Features of the Open and Closed Conformations of WT-RMPK and of MtbPK

The skeletal muscle PK is constitutively active, and the first K+-dependent PK in which free ligand [26], partial ligand complex [27,28], and complete ligand complex [17] structural data states have been reported. Of the 97 crystal structures of PKs deposited in PDB [29], only three are K+-independent [18,30,31]. However, only the structural data of apo-, substrate-complex, effector-complex, and substrate-effector complex are available for MtbPK [18]. It is well documented how ligands induce the closure of the active site cleft in WT-RMPK [17,20], but no information has been described in this respect in the K+-independent PKs. Figure 6A shows the overlapped structural conformations of WT-RMPK with partial ligand complex (cyan/open conformation) and complete ligand complex (green/closed conformation). In contrast, Figure 6B depicts the overlapped conformations of MtbPK in the substrate-effector complex (green) and the apoenzyme (cyan). The RMPK holoenzyme structural data could not be compared with the crystal structure of the PK apoenzyme, as it is the cat´s muscle-derived [26]. Instead, the structural characteristics of the same crystal´s unit cell were compared. Figure 6A shows subunit 1 (green) with the active site completely occupied (K+, Mg2+, oxalate, and Mg-ATP) was overlapped with subunit 2 (cyan) that exhibits a partially occupied active site (K+, Mg2+, and oxalate). B domain exhibited 0° and 41° angles of rotation relative to A domain in the former and latter subunits, indicating that the closure of the active site was either totally closed or totally open, respectively [17]. The distances between T120 and S76 in both conformations are shown. In WT-RMPK, T120 moves 7.4 Å in the closed conformation, whereas in MtbPK, no movement of the B domain is induced by the substrate-effector complex. In this condition, L75 is even 1 Å farther than L75 in the apoenzyme (Figure 6B).
This finding indicates that K+-dependent PKs follow a different mechanism to achieve the active conformation than the K+-independent PKs; T120 and L75 in the former and latter groups of PKs play different roles as indicated by the nearly null activity and differences in CD data described above.

2.6. Molecular Dynamics Simulations of WT-RMPK, T120L Mutant and MtbPK in the Absence or Presence of Ligands

RMSD was used to evaluate global flexibility of structures. This analysis considers the average amount of movement of backbone atoms throughout the entire protein. RMSF estimates localized flexibility of individual amino acid residues throughout the overall structure. Three replicates of RMSD with and without ligands of WT-RMPK, T120L mutant, and MtbPK are shown in Figure 7A,C,E, respectively. The flexibility of WT-RMPK without ligands increased ~5Å from 0 ns to 50 ns and showed fluctuations of 4 to 6 Å from 50 to 100 ns; whereas RMSD of the enzyme with ligands remained without change. In contrast, RMSD of T120L mutant with or without ligands kept constant through the time-lapse studied. RMSD of MtbPK showed none and minimal changes with and without ligands, respectively, through the same period of time. On the other hand, RMSF in the presence (red) or in the absence (blue) of ligands of WT-RMPK, T120L mutant, and MtbPK are shown in Figure 7 B,D,F, respectively. In WT-RMPK, a difference of ~4Å, was observed in the residues (116 to 223) corresponding to B domain (lid) with no ligands, compared to that with them. Contrary to observations in WT-PMPK, almost no differences were observed in the residues of the “lid” in T120L mutant and MtbPK.
According to T120L mutant, although a native-like structure was simulated, it did not exhibit the flexibility of the lid characteristic of the K+-dependent PKs. Therefore, T120 is a critical residue for the closure of the lid in the catalysis of the K+-dependent PKs. In coincidence with Figure 6, the overall results of MD simulations indicate that K+-dependent and independent PKs follow different mechanisms to achieve the active conformation.

2.7. Structural Changes Induced by the Substitution of Thr for Leu 120 to the RMPK

As described before, RMPK was the first K+-dependent enzyme in which the structural data for the complete ligand complex was reported. It showed how the binding of the nucleotide brings the enzyme to a totally closed or active conformation [17]. Residues R119, K206 and D177, from the “lid”, move 6.8, 10.4 and 6.8 Å, respectively, from their position in the open conformation and establish contacts to the nucleotide in the closed conformation. K206 NZ moves 10.4 Å to a position between the 2′- and 3′-hydroxyls of the ribose moiety; the guanidinium moiety of R119 moves 6.8 Å to H-bonding distance of the β-phosphate of the nucleotide; and the carboxylate of D177 travels 6.8 Å to H-bonding contact with a water of hydration on the Mg2+ of Mg-ATP. The adenine ring fits into a pocket bounded on one side by His 77 [17].
Figure 8 compares the superposition of the molecular dynamic models of WT-RMPK and T120L mutant in the initial structures with Mg2+, PEP, Mg-ADP and K+ to understand the differences caused by this single mutation. In the T120L mutant, a structural distortion of the guanidium moiety of R119 was observed. It moved to the imidazole moiety of His77 instead of the β-phosphate of ADP. NZ of K206 moved to position of H-bonding distance of the ADP α-phosphate instead of moving between the 2′- and 3′-hydroxyls of the ribose moiety. The missing and new interactions between nucleotide and residues from the “lid” could explain the almost null activity of the T120L mutant.

2.8. Different Main Chain Orientations of Residue Lys117 in the E117K Mutant and the Corresponding Lys72 in the Apo and Holoenzyme of MtbPK

The structure of the E117K mutant shows a tetramer in the asymmetric unit. Subunits A and B are the most complete ones with 519 residues out of 530, missing only the first 12 residues at the N-terminus. Due to poor density, subunit C lacks 121 amino acids, missing 11 residues at the N-terminus and from 121 to 203. Likewise, subunit D lacks 11 residues at the N-terminus, and it is missing the amino acids 118 to 215.
A superposition of subunits A and B shows almost identical structures with RMSD of 0.25 Å. However, amino acids 111 to 220 are the most different in this alignment and do not superimpose perfectly. Thus, a different conformation of the NZ atom in K117 is observed in both subunits. The NZ atom in subunit B is within H-bond distance with a water molecule, which coordinates to the Mg2+ ion near the oxalate ligand. Such an interaction does not occur in subunit A, which has much fewer water molecules in the vicinity of the oxalate.
In contrast, K117 in subunit C exhibits similar conformation to the one observed in subunit B, but shifted by about 0.5 Å. Consequently, the “lid” domain was not observed, likely due to high flexibility and a lack of complete electron density. Finally, residues K117-215 in subunit D could not be interpreted in the electron density. These couple motions of the “lid” domains could define functional transitions.
When the MtbPK structure was compared, it was evident that K72 (corresponding to the residue 117 of the RMPK mutant) is positioned towards the active site in the holo- and flips out in the apoenzyme. The partially occupied mutant E117K mutant is also facing toward the active site (Figure 9A). The distances between NZ of K72 (or K117) to O1 of the oxalate are 7.4 and 7.2 Å, respectively.
On the other hand, it is relevant to mention that T113 to T120 are localized in a loop bending region that participates in the closure or acquisition of the active catalytic conformation of K+-dependent holo-WT-RMPK (Figure 6A), whereas no bending is observed in this loop region in K+-independent holo-MtbPK (Figure 6B). In the latter, loop residues, L68-L75, exhibit the same orientation in the apo- or holoenzyme, except for K72 that flips into or out of the active site, respectively (Figure 9B). This finding indicates that a completely different mechanism is observed in the K+-independent MtbPK, where K72 switches from the inactive (apo) to the active (holo) conformation.

2.9. Different Interaction of MgATP to Oxalate, Mg2+ and K+ in the Active Sites of WT-RMPK and MtbPK

In WT-RMPK, the γ-phosphate of Mg-ATP exhibits direct interaction with oxalate, Mg2+, and K+ (Figure 10A), whereas in MtbPK these interactions occur through water molecules (Figure 10B). The latter observation suggests that water molecules may have a relevant role in the catalysis of the K+-independent enzymes.

2.10. Active Sites of the Holo and Apo MtbPK

As shown in Figure 11, K72 in the apoenzyme flips out of the active site. In contrast, in the holoenzyme, K72 faces toward the active site, interacting with a network of water molecules connecting the lid (domain B), substrates (Mg2+, oxalate, Mg-ATP, and K+), and A domain residues. Therefore, K72 is a key residue in the water network formation. A similar water network is observed in the crystal structure of K+-independent PK of Pseudomonas aeruginosa [31]. Consequently, this data confirmed that the acquisition of the active conformation of the K+-independent PKs is via K72 switching from the inactive to the active conformation.

3. Conclusions

In summary, the aim of this work was to understand the individual and combined role of changing specific residues in WT-RMPK to those observed in K+-independent PKs. K+-independent activity increased with residues L113 and Q114 to K117, with activity in the absence of K+ reaching half WT-RMPK activity with K+. In contrast, the single mutant T120L or the addition of this mutation to other mutants leads to a nearly inactive enzyme. Single mutants, T113L and K114Q, were K+-dependent with kinetic constants similar to WT-RMPK, except that ion selectivity for T113L was modified as described elsewhere [21]. Furthermore, residues in positions 117 and 120 appear decisive in the differences between the K+-dependent and -independent PKs. When Leu substituted T120, the enzyme exhibited different structural features with nearly null activity and did not close the lid to acquire its active conformation. Although simulations indicated it to be folded entirely, no movement of the residues of the “lid” was observed by MD. We conclude that T120 is involved in folding or stabilizing the protein and closing the lid of the active site in the K+-dependent PKs. When Lys was substituted for E117, the enzyme became K+-independent and acquired its active conformation by a different mechanism. The holoenzyme with E117 acquires its active conformation by closing the B domain over A domain, whereas in the K+-independent MtbPK, K72 (or K117) switches between the active and inactive conformations. In the apoenzyme, K72 flips out of the active site and, in the holoenzyme, K72 faces into the active site, connecting the substrates, residues of B domain, and A domain via water molecules. Finally, the K+ binding site is highly conserved in both branches, K+-dependent and K+-independent, of the phylogenetic tree of the PK family. In most K+-independent PKs, three of the four coordination residues of K+ are preserved. The role of K+ is not catalytically essential like in K+-dependent PKs, but may exert a weak ionic strength activation, as seen in the K+-independent mutants studied in this work.

4. Experimental Procedures

4.1. Cloning and Expression of RMPK Mutants

In all cases (Table 3), Fw and Rv mutagenic primers and external primers were used for PCR reaction. Mutagenesis was performed using the following PCR conditions: 94 °C for 5 min, 25 cycles for 30 s at 94 °C, 30 s at 62 °C, 1.5 min at 72 °C and 5 min at 72 °C. The PCR products were cloned in the pJET 1.2 vector, and mutagenesis was confirmed by automated DNA sequencing of the complete genes. The pJET 1.2 vector containing the verified sequence for each mutant was digested with EcoRI and XbaI and subcloned into the pTRC99A plasmid. All constructs with the desired mutations were transformed into competent PB25 cells for protein expression.

4.2. Cell Growth and Purification of RMPK Mutants

LB medium containing 100 μg/mL ampicillin, 50 μg/mL kanamycin, and 30 μg/mL chloramphenicol was inoculated either with PB25-E117K, or the other mutants. Expression was induced with 1 mM isopropyl 1-thio-β-D-galactopyranoside and 20 g/liter of lactose at an OD600 of about 0.8. The enzyme was purified as in [20] with some modifications. Cells were treated with lysozyme and an osmotic shock. DNAse and a protease inhibitor mixture were added to the spheroplasts, which were then broken by sonication. The suspension was centrifuged for 20 min at 20,000× g. The supernatant was precipitated with ammonium sulfate at 37%, and the resultant supernatant was collected. A second precipitation with ammonium sulfate at 65% saturation was prepared. The pellet was collected, suspended, and desalted by dialysis. The final steps involved ion exchange chromatography in DEAE and carboxymethyl-Sepharose columns previously equilibrated with HEPES pH 6.4 to 6.6 for all mutants except for the mutants containing the substitution of Thr 120 by Leu (T120L, K114Q/E117K/T120L and T113L/K114Q/E117K/T120L mutants). These mutants were poorly expressed, and about half of the proteins were in inclusion bodies. To obtain more solubilized enzymes, 500 mM NaCl was added to the buffer before lysis. The columns were equilibrated at pH 6.9, and 10% glycerol was included in all purification steps. Unlike the ~15 mg/ 4 L LB broth media recovered for other mutants, those with the T120L mutation recovered only ~2 mg/4 L. The mutants were about 95% pure, as indicated by SDS-PAGE (not shown). Purified K114Q, E117K, T113L/E117K, and T113L/K114Q/E117K mutants were stored as ammonium sulfate suspensions at 4 °C. Ammonium sulfate-free T120L exhibited no activity. Therefore, further purifications of T120L, K114Q/E117K/T120L and T113L/K114Q/E117K/T120L were stored in 20% glycerol at −20 °C in aliquots of 50 μl. However, K114Q/E117K/T120L was completely inactive and T120L and T113L/K114Q/E117K/T120L were nearly inactive. K114Q/E117K [15] and T113L [21] were not included in this study because they were studied previously.

4.3. Assays of Pyruvate Kinase Activity

Ammonium sulfate suspensions of WT-RMPK and hog muscle lactate dehydrogenase (LDH) were obtained as ammonium sulfate suspensions from SIGMA. Ammonium sulfate-free enzymes (WT-RMPK, LDH, K114Q, E117K, T113L/E117K, and T113L/K114Q/E117K) were obtained as described in [33]. Contaminating NH4+, Na+, and K+ in reaction mixtures were below the detection limit (10 μM) as indicated in [33]. Mutants containing the substitution of Thr by Leu in position 120 (T120L, K114Q/E117K/T120L and T113L/K114Q/E117K/T120L mutants) exhibited none or hardly any activities. The formation of pyruvate was measured at 25 °C in a coupled system with LDH and NADH [34]. The specific activity was not increased by the inclusion of 5-fold higher concentrations of LDH.
The reaction mixture contained 25 mM HEPES-(CH3)4NOH, pH 7.4, and the concentrations of monovalent cations (Li+, Na+, K+, NH4+, Rb+, and Cs+), Mg2+, phosphoenolypyruvate (PEP), and ADP indicated in each figure legend. The Mg-ADP complexes and free Mg2+ concentrations were calculated using the software CHELATOR [35]. The ionized PEP concentrations were calculated considering a pK value of 6.3 [36]. With the exception of Figure 2 and Figure 3, (CH3)4N+ was used to maintain at 250 mM constant ionic strength. The experiments were carried out at 25 °C, and the reaction was initiated with the addition of either RMPK or any mutant.
Protein concentrations of RMPK and of mutants were determined by measuring the absorbance at 280 nm using the extinction coefficients of 0.54 mL mg−1 cm−1 [37].

4.4. Circular Dichroism Experiments

CD measurements in the far-ultraviolet (UV) region were carried out on a Jasco J-715 spectropolarimeter. A 0.1-cm quartz cell was used and the experiments were conducted at 25 °C. The concentration of WT-RMPK and of RMPK mutants: K114Q, E117K, T120L, T113L/E117K, T113L/K114Q/E117K was K114Q/E117K/T120L, and T113L/K114Q/E117K/T120L was 100 μg/mL. Spectral scans were run from 190 to 260 nm at intervals of 1 nm, scan rates of 20 nm/min, and a time constant of 16 s. In each case, bandwidth was 1 nm, and three scans were accumulated. The spectra of blanks were subtracted from those that contained the protein. CD is expressed as θMRW (mean residual weight ellipticity).

4.5. Molecular Dynamics Simulation

4.5.1. Preparation of Initial Coordinate Files

The coordinates corresponding to the structure of RMPK and MtbPK were obtained from the Protein Data Bank (PDB, http://www.rcsb.org) (accessed on 13 January 2022). The structure used was 1A49.pdb refined at 2.1 Å [17] for RMPK and PDB entry 5ws9 for MtbPK (refined at 1.9 Å). The structures had the hydrogens removed before the initial preparation of the molecular dynamics (MD). The T120L mutant was constructed from the same structure (PDB entry 1A49.pdb) using the PyMOL program [38], and subsequently refined with Rosetta relax 3.4 [39]. The structures of the ligands ADP, ATP, oxalate, and PEP were obtained from the PDB co-crystallized structure when available or constructed using Spartan’16 Parallel Suite for Linux software. All structures were minimized using Gaussian 09, revision A.02 at DTP B3LYP/3-21G level of theory. Partial charges and force fields of the ligands were generated using the antechamber program in AmberTools19.

4.5.2. Molecular Dynamics Simulation

All the structural complexes were verified, cleaned, and ordered with the pdb4amber scrip before starting the preparation in order to generate suitable topologies from the LEaP module of AmberTools19 [40,41,42]. Each structure and complex were subjected to the following protocol: hydrogen and other missing atoms were added using the LEaP module with the leaprc.protein.ff19SB parameter set, Cl counterions were added to neutralize the system; the complexes were then solvated in an octahedral box of explicit TIP3P model water molecules localizing the box limits at 12 Å from the protein surface. MD simulations were performed at 1 atm and 298.15 K, maintained with the Berendsen barostat and thermostat, using periodic boundary conditions and particle mesh Ewald sums (grid spacing of 1 Å) for treating long-range electrostatic interactions with a 10 Å cutoff for computing direct interactions. The SHAKE algorithm was used to satisfy bond constraints, allowing the employment of a two fs time step for the integration of Newton’s equations as recommended in the Amber package [41,42]. Amber leaprc.protein.ff19SB force field [43] parameters were used for all residues. All calculations were made using graphics processing units (GPU) accelerated MD engine in AMBER (pmemd.cuda), a program package that runs entirely on CUDA-enabled GPUs [44]. The protocol consisted of performing a minimization of the initial structure, followed by 2500 heating steps (50 ps) and depressurization (50 ps) at 298.15 K and 1.0 atm pressure, respectively. Finally, the system is equilibrated with 250,000 steps (500 ps) before starting the production of MD. The production of the MD consisted of 100 ns (50,000,000 steps) for each complex. Frames were saved at ten ps intervals for subsequent analysis. Three replicates of the molecular dynamics simulations were carried out. Each of them starting from the minimization step. All data from molecular dynamics studies are available at the following link http://bq.unam.mx/ijms-23-01347/ (accessed on 13 January 2022).

4.5.3. Trajectory Analysis

All analyses were done using CPPTRAJ [45] part of AmberTools19 utilities and OriginPro 9.1. The calculations of RMSD and Root Mean Square Fluctuations (RMSF) were made, considering the C, CA, and N; for the distances, only the CA was used. The charts were built with Origin 9.0, and the trends were adjusted with the function processing smooth (method Lowess). VMD and PyMOL were used to visualize and create the images from the MD.
All calculations were made using a system HP Cluster Platform 3000 SL supercomputer “MIZTLI” with a processing capacity of 118 TFlop/s. It has 5312 Intel E5-2670 processing cores, 16 NVIDIA m2090 cards, a total RAM of 15,000 Gbytes, and a mass storage system of 750 Terabytes (http://www.super.unam.mx/) (accessed on 13 January 2022).

4.6. Crystallization of the Mutant E117K

The pure E117K mutant was dialyzed into a buffer containing 50 mM Tris pH 7.0. The enzyme was concentrated in the same buffer, and 1.0 mM ATP, 1.0 mM MgCl2, and 0.2 mM oxalate were added to the protein solution, which ended with a final concentration of 7 mg/mL. The vapor diffusion method in a hanging drop mode was used to set up the crystallization drops, using 24 well plates. Wells contained 500 μL of crystallization buffer, and drops were 4.0 μL with a 1:1 ratio of enzyme to crystallization buffer. After the first set of crystallization experiments, the best-looking crystals appeared in condition No. 46 of the Hampton Research Crystal Screen I (Journey Aliso Viejo, CA, USA). This solution contains 0.2 M calcium acetate, 0.1 M sodium cacodylate pH 6.5, and 18% w/v Polyethylene Glycol 8000. However, diffracting-quality crystals grew from drops containing 0.2 M magnesium acetate instead of calcium acetate. For this experiment, the well´s crystallization solution was covered with a 1:1 mixture of paraffin and silicone oils from Hampton Research (Journey, Aliso Viejo, CA, USA).

4.7. Data Collection and Reduction

X-ray diffraction data for the best crystals were collected on beamline X6A at the National Synchrotron Light Source (Upton, NY, USA), at 100 K using 25% glycerol as cryoprotectant. A Quantum 210 CCD detector was used to collect 180 images with an oscillation range of 1.0° and 230 mm from the detector. Frames were integrated using XDS [46], with an initial space group P222, which was confirmed as P212121 using pointless in CCP4 [47]. Data were scaled with SCALA from CCP4 [47], and the mtz file was analyzed using Xtriage from Phenix [48], which suggested no problems with the data or space group selected. The Matthew’s coefficient of 2.2 and 0.45 fraction of solvent suggested four copies of the E117K pyruvate kinase mutant in the asymmetric unit. Data collection statistics are shown in Table 4.

4.8. Molecular Replacement and Model Refinement

Initial phases for the orthorhombic P212121 were determined by molecular replacement using chain A of the pyruvate kinase´s atomic coordinates from rabbit muscle (PDB 1AQF) with the program Phaser from Phenix [49]. A tetramer of the E117K mutant was obtained as the solution. The initial tetrameric model was refined using rigid-body refinement, and non-crystallographic (NCS) restraints were then applied. The resulting model underwent iterative rounds of refinement with Phenix.refine and manual rebuilding with Coot [50,51]. After the initial stages of refinement, oxalate, acetate, solvent, and ions were included, and the CCP4 cif files for all molecules except water were used during refinement. Refinement was finished when the best possible Rfree and Rwork statistics were obtained, and further adjustments no longer improved the model. Validation of the Mg2+ metal ions was done using the CheckMyMetal server [52].

Author Contributions

L.R.-S. designed, purified, performed kinetic of WT-RMPK and K114Q mutant, evaluated the data, performed CD experiments of all mutants, prepared the figures and wrote the manuscript with input of C.G.-M., G.H.-A. build up double, triple and quadruple mutants and purified T120L mutant and performed kinetic experiments. A.L.-M., purified and performed kinetic experiments of E117K and E117K/T113L/K114Q mutants. P.A.G.-C., purified and performed kinetic experiments of T113L/E117K mutant. C.R.-M., build up mutants K114Q/E117K/T120L and T113L/K114Q/E117K/T120L, and purified them. A.V.-S. and L.R.-S. performed ionic strength experiments of WT-RMPK and all active mutants. C.G.-M. performed the analyses of structural data of MtbPK and WT-RMPK and made the structural figures. A.R.-R. and A.R.-H. obtained the crystals of mutant E117K, collected X-ray data and resolved the structure of the protein. M.G.-A. performed the molecular dynamics of WT-RMPK, T120L mutant, and MtbPK. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Dirección General de Apoyo al Personal Académico-Universidad Nacional Autónoma de México, Grant IN201021 (LR-S), Grant IN205821 (GH-A), Grant IA204116 (MG-A) and the postdoctoral fellowship of (CG-M). Super computer Grant LANCAD UNAM DGTIC-313 (MG-A). AR-R thanks (FORDECYT-PRONACES/87163).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank Vivian Stojanoff at the X6A beamline for assistance with X-ray data collection. The X6A beamline and the National Synchrotron Light Source, Brookhaven National Laboratory, are funded by NIH–NIGMS under agreement GM-0080 and the US Department of Energy under contract No. DE-AC02-98CH10886, respectively. They also thank LANEM-IQ for characterization of the E117K mutant crystals. AR-R thanks (FORDECYT-PRONACES/87163). The authors also thank Nallely Cabrera for technical assistance, Juan Pablo Pardo, José García-Trejo, and Marietta Tuena de Gómez-Puyou for valuable suggestions and Willian Newell and Andrés Ramírez-Silva for careful revision of the manuscript.

Conflicts of Interest

The authors declare that they have no conflict of interest with the contents of this article.

Abbreviations

Pyruvate kinasePK
Wild type rabbit muscle pyruvate kinaseWT-RMPK
Choline(CH3)3N+-CH2-CH2-OH
Circular dichroismCD
Molecular dynamicsMD
Mycobacterium tuberculosis pyruvate kinaseMtbPK
N-(2-Hydroxyethyl)piperazine-N′-(2-ethanesulfonic acid)HEPES
phosphoenolpyruvatePEP
Luria BertaniLB
hog muscle lactate dehydrogenaseLDH

References

  1. Kachmar, J.F.; Boyer, P.D. Kinetic analysis of enzyme reactions: II. The potassium activation and calcium inhibition of pyruvic phosphoferase. J. Biol. Chem. 1953, 200, 669–682. [Google Scholar] [CrossRef]
  2. Ramírez-Silva, L.; Ferreira, S.T.; Nowak, T.; Tuena de Gómez-Puyou, M.; Gómez-Puyou, A. Dimethyl sulfoxide promotes the K+-independent activity of pyruvate kinase and the acquisition of the active catalytic conformation. Eur. J. Biochem. 2001, 268, 3267–3274. [Google Scholar] [CrossRef] [PubMed]
  3. Kayne, F.J. Thallium (I) activation of pyruvate kinase. Arch. Biochem. Biophys. 1971, 143, 232–239. [Google Scholar] [CrossRef]
  4. Nowak, T. Conformational changes required for pyruvate kinase activity as modulated by monovalent cations. J. Biol. Chem. 1976, 251, 73–78. [Google Scholar] [CrossRef]
  5. Benziman, M. Factors affecting the activity of pyruvate kinase of Acetobacter xylinum. Biochem. J. 1969, 112, 631–636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Liao, C.L.; Atkinson, D.E. Regulation at the phosphoenolpyruvate branchpoint in Azotobacter vinelandii: Pyruvate kinase. J. Bacteriol. 1971, 106, 37–44. [Google Scholar] [CrossRef] [Green Version]
  7. Smart, J.B.; Pritchard, G.G. Regulation of pyruvate kinase from Propionibacterium shermanii. Arch. Microbiol. 1979, 122, 281–288. [Google Scholar] [CrossRef]
  8. Kapoor, R.; Venkitasubramanian, T.A. Glucose 6-phosphate activation of pyruvate kinase from Mycobacterium smegmatis. Biochem. J. 1981, 193, 435–440. [Google Scholar] [CrossRef] [Green Version]
  9. Waygood, E.B.; Rayman, M.K.; Sanwal, B.D. The control of pyruvate kinase of Escherichia coli. II. Effectors and regulatory properties of the enzyme activated by ribose 5-phosphate. Can. J. Biochem. 1975, 53, 444–454. [Google Scholar] [CrossRef]
  10. Busto, F.; Del Valle, P.; Soler, J. Some kinetic properties of pyruvate kinase from Phycomyces blakesleeanus. Biochem. Cell Biol. 1988, 66, 148–157. [Google Scholar] [CrossRef]
  11. Jetten, M.S.M.; Gubler, M.E.; Lee, S.H.; Sinskey, A.J. Structural and functional analysis of pyruvate kinase from Corynebacterium glutamicum. Appl. Environ. Microbiol. 1994, 60, 2501–2507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Steiner, P.; Fussenegger, M.; Bailey, J.E.; Sauer, U. Cloning and expression of the Zymomonas mobilis pyruvate kinase gene in Escherichia coli. Gene 1988, 220, 31–38. [Google Scholar] [CrossRef]
  13. Schramm, A.; Siebers, B.; Tjaden, B.; Brinkmann, H.; Hensel, R. Pyruvate kinase of the hyperthermophilic Crenarchaeote Thermoproteus tenax: Physiological role and phylogenetic aspects. J. Bacteriol. 2000, 182, 2001–2009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Johnsen, U.; Hansen, T.; Schönheit, P. Comparative analysis of pyruvate kinases from the hyperthermophilic archaea Archaeoglobus fulgidus, Aeropyrum pernix, and Pyrobaculum aerophilum. J. Biol. Chem. 2003, 278, 25417–25427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Laughlin, L.T.; Reed, G. The monovalent cation requirement of rabbit muscle pyruvate kinase is eliminated by substitution of lysine for glutamate 117. Arch. Biochem. Biophys. 1997, 348, 262–267. [Google Scholar] [CrossRef]
  16. Oria-Hernández, J.; Riveros-Rosas, H.; Ramírez-Silva, L. Dichotomic phylogenetic tree of the pyruvate kinase family. K+-dependent and –independent enzymes. J. Biol. Chem. 2006, 281, 30717–30724. [Google Scholar] [CrossRef] [Green Version]
  17. Larsen, T.M.; Benning, M.M.; Rayment, I.; Reed, G.H. Structure of the bis(Mg2+)-ATP-oxalate complex of the rabbit muscle pyruvate kinase at 2.1 Å resolution: ATP binding over a barrel. Biochemistry 1998, 37, 6247–6255. [Google Scholar] [CrossRef]
  18. Zhong, W.; Cui, L.; Goh, B.C.; Cai, O.; Ho, P.; Chionh, Y.H.; Yuan, M.; El Sahili, A.; Fothergill-Gilmore, L.; Walkinshaw, M.D.; et al. Allosteric pyruvate kinase-based “logic gate” synergistically senses energy and sugar levels in Mycobacterium tuberculosis. Nat. Commun. 2017, 8, 1986. [Google Scholar] [CrossRef]
  19. Suelter, C.H. Monovalent cations in enzyme-catalyzed reactions. In Metal Ions in Biological Systems; Sigel, H., Ed.; Dekker: New York, NY, USA, 1974; Volume 3, pp. 201–251. [Google Scholar]
  20. Oria-Hernández, J.; Cabrera, N.; Pérez-Montfort, R.; Ramírez-Silva, L. Pyruvate kinase revisited. The activating efffect of K+. J. Biol. Chem. 2005, 280, 37924–37929. [Google Scholar] [CrossRef] [Green Version]
  21. Ramírez-Silva, L.; Guerrero-Mendiola, C.; Cabrera, N. The importance of polarity in the evolution of the K+-binding site of pyruvate kinase. Int. J. Mol. Sci. 2014, 15, 22214–22226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Di Cera, E. A structural perspective on enzymes activated by monovalent cations. J. Biol. Chem. 2006, 281, 1305–1308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Page, M.J.; Di Cera, E. Role of Na+ and K+ in enzyme function. Physiol. Rev. 2006, 86, 1049–1092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Guerrero-Mendiola, C.; García-Trejo, J.J.; Encalada, R.; Saavedra, E.; Ramírez-Silva, L. The contribution of two isozymes to the pyruvate kinase activity of Vibrio cholera: One K+-dependent constitutively active and another K+-independent with essential allosteric activation. PLoS ONE 2017, 12, e0178673. [Google Scholar] [CrossRef] [PubMed]
  25. Geenfield, N. Using circular dichroism spectra to estimate protein secondary structure. Nat. Protoc. 2006, 1, 2876–2890. [Google Scholar] [CrossRef] [PubMed]
  26. Allen, S.C.; Muirhead, H. Refined three-dimensional structure of cat muscle (M1) pyruvate kinase at a resolution of 2.6 Å. Acta Crystallogr. Sect. D Biol. Crystallogr. 1996, 52, 499–504. [Google Scholar] [CrossRef]
  27. Larsen, T.M.; Laughlin, L.T.; Holden, H.M.; Rayment, I.; Reed, G.H. Structure of Rabbit Muscle Pyruvate Kinase Complexed with Mn2+, K+, and Pyruvate. Biochemistry 1994, 33, 6301–6309. [Google Scholar] [CrossRef]
  28. Larsen, T.M.; Benning, M.M.; Wesenberg, G.E.; Rayment, I.; Reed, G.H. Ligand-Induced Domain Movement in Pyruvate Kinase: Structure of the Enzyme from Rabbit Muscle with Mg2+, K+, and L-Phospholactate at 2.7 Å Resolution. Arch. Biochem. Biophys. 1997, 354, 199–206. [Google Scholar] [CrossRef]
  29. Schormann, N.; Hayden, K.L.; Lee, P.; Banerjee, S.; Chattopadhyay, D. An overview of structure, function, and regulation of pyruvate kinases. Protein Sci. 2019, 28, 1771–1784. [Google Scholar] [CrossRef]
  30. Solomons, T.G.; Johnsen, U.; Schönheit, P.; Davies, C. 3-Phosphoglycerate is an Allosteric Activator of Pyruvate Kinase from the Hyperthermophilic Archaeon Pyrobaculum aerophilum. Biochemistry 2013, 52, 5865–5875. [Google Scholar] [CrossRef] [PubMed]
  31. Abdelhamid, Y.; Brear, P.; Greenhalgh, J.; Chee, X.; Rahman, T.; Welch, M. Evolutionary plasticity in the allosteric regulator binding site of pyruvate kinase isoform PykA from Pseudomonas aeruginosa. J. Biol. Chem. 2019, 294, 15505–15528. [Google Scholar] [CrossRef]
  32. Roe, D.R.; Cheatham, T.E., III. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef] [PubMed]
  33. Ramírez-Silva, L.; Tuena de Gómez-Puyou, M.; Gómez-Puyou, A. Water-Induced transitions in the K+ requirements for the activity of Pyruvate Kinase entrapped in Reverse Micelles. Biochemistry 1993, 32, 5332–5338. [Google Scholar] [CrossRef] [PubMed]
  34. Büchner, T.; Pleiderer, G. Pyruvate kinase from muscle. In Methods in Enzymology; Colowick, S., Kaplan, N., Eds.; Academic Press: New York, NY, USA, 1955; Volume 1, pp. 435–440. [Google Scholar]
  35. Shoemakers, T.J.M.; Visser, G.J.; Flik, G.; Theuvent, P.R. CHELATOR: An improved method for computing metal ion concentrations in physiological solutions. Biotechniques 1992, 12, 870–879. [Google Scholar]
  36. Doughtery, T.M.; Cleland, W.W. pH studies of the chemical mechanism of rabbit muscle pyruvate kinase. 2. Physiological substrates and phosphoenol-αbutyrate. Biochemistry 1985, 24, 5875–5880. [Google Scholar]
  37. Cottam, G.L.; Hollenberg, P.F.; Coon, M.J. Subunit structure of rabbit muscle pyruvate kinase. J. Biol. Chem. 1969, 244, 1481–1486. [Google Scholar] [CrossRef]
  38. DeLano, W.L. Use of PYMOL as a communications tool for molecular science. Am. Chem. Soc. 2004, 228, U313–U314. [Google Scholar]
  39. Ferrari, A.J.R.; Gozzo, F.C.; Martinez, L. Statistical force-field for structural modeling using chemical cross-linking/mass spectrometry distance constraints. Bioinformatics 2019, 35, 3005–3012. [Google Scholar] [CrossRef]
  40. David, A.C.; Josh, B.; Betz, R.M.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; et al. AMBER 2015. 2015. Available online: https://ambermd.org/doc12/Amber15.pdf (accessed on 13 January 2022).
  41. Case, D.A.; Darden, T.A.; Cheatham, T.E., III; Simmerling, C.L.; Wang, J.; Duke, R.E.; Luo, R.; Walker, R.C.; Zhang, W.; Merz, K.M.; et al. AMBER 12. 2012. Available online: https://ambermd.org/doc12/Amber12.pdf (accessed on 13 January 2022).
  42. Case, D.A.; Cheatham, T.E., III; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [Green Version]
  43. Walker, R.C.; Crowley, M.F.; Case, D.A. The implementation of a fast and accurate QM/MM potential method in Amber. J. Comput. Chem. 2008, 29, 1019–1031. [Google Scholar] [CrossRef] [Green Version]
  44. Tian, C.; Kasavajhala, K.; Belfon, K.A.A.; Raguette, L.; Huang, H.; Migues, A.N.; Bickel, J.; Wang, Y.; Pincay, J.; Wu, Q.; et al. ff19SB: Amino-Acid-Specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution. J. Chem. Theory Comput. 2020, 16, 528–552. [Google Scholar] [CrossRef]
  45. Salomon-Ferrer, R.; Gotz, A.W.; Poole, D.; Le Grand, S.; Walker, R.C. Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. J. Chem. Theory Comput. 2013, 9, 3878–3888. [Google Scholar] [CrossRef] [PubMed]
  46. Kabsch, W. Xds. Acta. Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 125–132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Evans, P.R. Scaling and assessment of data quality. Acta. Crystallogr. Sect. D Biol. Crystallogr. 2006, 62, 72–82. [Google Scholar] [CrossRef] [PubMed]
  48. Adams, P.D.; Afonine, P.V.; Bunkóczi, G.; Chen, V.B.; Davis, I.W.; Echols, N.; Headd, J.J.; Hung, L.-W.; Kapral, G.J.; Grosse-Kunstleve, R.W.; et al. A comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 213–221. [Google Scholar] [CrossRef] [Green Version]
  49. McCoy, A.J.; Grosse-Kunstleve, R.W.; Adams, P.D.; Winn, M.D.; Storoni, L.C.; Read, R.J. Phaser crystallographic software. J. Appl. Crystallogr. 2007, 40, 658–674. [Google Scholar] [CrossRef] [Green Version]
  50. Afonine, P.V.; Grosse-Kunstleve, R.W.; Echols, N.; Headd, J.J.; Moriarty, N.W.; Mustyakimov, M.; Terwilliger, T.C.; Urzhumtsev, A.; Zwart, P.H.; Adams, P.D. Towards automated crystallographic structure refinement with phenix refine. Acta. Crystallogr. D. Biol. Crystallogr. 2012, 68, 352–367. [Google Scholar] [CrossRef] [Green Version]
  51. Emsley, P.; Cowtan, K. Coot: Model-building tools for molecular graphics. Acta Crystallogr. Sect. D Biol. Crystallogr. 2004, 60, 2126–2132. [Google Scholar] [CrossRef] [Green Version]
  52. Zheng, H.; Cooper, D.R.; Porebski, P.J.; Shabalin, I.G.; Handing, K.B.; Minor, W. CheckMyMetal: A macromolecular metal-binding validation tool. Acta Crystallogr. Sect. D Biol. Crystallogr. 2017, 73, 223–233. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Ribbon representation of the boundary region between domains A (gray) and B (blue) of WT-RMPK. According to Oria-Hernández et al. (2006) [16], short loop residues T113, K114, E117, and T120 co-evolve in a consensus sequence in the family of the K+-dependent PKs in 80% of the 121 sequences; whereas L113, Q114, K117 and (L, V, I) 120 co-evolve in the K+-independent PKs in 77% of the 106 sequences. Residues of RMPK are illustrated in sticks and K+ in purple. The figure was constructed with PDB 1A49, subunit A.
Figure 1. Ribbon representation of the boundary region between domains A (gray) and B (blue) of WT-RMPK. According to Oria-Hernández et al. (2006) [16], short loop residues T113, K114, E117, and T120 co-evolve in a consensus sequence in the family of the K+-dependent PKs in 80% of the 121 sequences; whereas L113, Q114, K117 and (L, V, I) 120 co-evolve in the K+-independent PKs in 77% of the 106 sequences. Residues of RMPK are illustrated in sticks and K+ in purple. The figure was constructed with PDB 1A49, subunit A.
Ijms 23 01347 g001
Figure 2. Effect of ionic strength on the activities of WT-RMPK (A), E117K (B), T113L/E117K (C), and T113L/K114Q/E117K (D) mutants induced by the indicated monovalent cations. In the absence of monovalent cations, the reaction mixtures contained ionic strengths of 99.5, 50, and 126.8 mM in (BD), respectively. The activities in these conditions for (BD) were 14 ± 3, 10 ± 1 and 34 ± 4, respectively. In the presence of 90 mM K+, the activity for (A) was 263 ± 15 and ionic strength 143 mM. These activities were normalized to 100%, which corresponded to the control activity for each PK. Besides 25 mM HEPES-(CH3)4NOH pH 7.4, 0.2 mM NADH, 8 μg/ml LDH and a reaction mixture with saturating concentrations of substrates, Li+ (), Na+ (▲), K+ (∎), NH4+ (), Rb+ (▼), Cs+ (), choline (●) or (CH3)4N+ () were added to complete the ionic strengths indicated in the figure. Reaction mixtures contained: (A) 0.64 mM PEP3−, 2.5 mM Mg-ADP, 1.69 mM Mg2+free; (B) 0.63 mM PEP3−, 2 mM Mg-ADP and 19 mM Mg2+free; (C) 0.54 mM PEP3−, 0.7 mM Mg-ADP and 3 mM Mg2+free and (D) 0.67 mM PEP3−, 1.84 mM Mg-ADP and 24 mM Mg2+free. Assays were performed at 25 °C, and the reaction was started by the addition of PK. Amounts of PK ranged from 0.1 to 0.2 μg/mL. Standard deviation bars of three experiments are shown.
Figure 2. Effect of ionic strength on the activities of WT-RMPK (A), E117K (B), T113L/E117K (C), and T113L/K114Q/E117K (D) mutants induced by the indicated monovalent cations. In the absence of monovalent cations, the reaction mixtures contained ionic strengths of 99.5, 50, and 126.8 mM in (BD), respectively. The activities in these conditions for (BD) were 14 ± 3, 10 ± 1 and 34 ± 4, respectively. In the presence of 90 mM K+, the activity for (A) was 263 ± 15 and ionic strength 143 mM. These activities were normalized to 100%, which corresponded to the control activity for each PK. Besides 25 mM HEPES-(CH3)4NOH pH 7.4, 0.2 mM NADH, 8 μg/ml LDH and a reaction mixture with saturating concentrations of substrates, Li+ (), Na+ (▲), K+ (∎), NH4+ (), Rb+ (▼), Cs+ (), choline (●) or (CH3)4N+ () were added to complete the ionic strengths indicated in the figure. Reaction mixtures contained: (A) 0.64 mM PEP3−, 2.5 mM Mg-ADP, 1.69 mM Mg2+free; (B) 0.63 mM PEP3−, 2 mM Mg-ADP and 19 mM Mg2+free; (C) 0.54 mM PEP3−, 0.7 mM Mg-ADP and 3 mM Mg2+free and (D) 0.67 mM PEP3−, 1.84 mM Mg-ADP and 24 mM Mg2+free. Assays were performed at 25 °C, and the reaction was started by the addition of PK. Amounts of PK ranged from 0.1 to 0.2 μg/mL. Standard deviation bars of three experiments are shown.
Ijms 23 01347 g002
Figure 3. Effect of the addition of various concentrations of K+ (or Rb+) and (CH3)4N+ on the activities of WT-RMPK (A), E117K (B), T113L/E117K (C) and T113L/K114Q/E117K (D) mutants. Ionic strength and activities with 90 mM K+ (A) and without monovalent cations (BD) were those indicated in the legend of Figure 2. The normalized control activities were shown in white columns in (AD). White striped columns were in the presence of 25, 50, 75, and 100 mM (CH3)4N+ concentrations without K+ (or Rb+). In addition to the indicated (CH3)4N+ concentration, light gray, gray, dark gray, and black striped columns were in the presence of the 25, 50, 75, and 100 mM K+ (or Rb+) concentrations, respectively. The highest activities of RMPK mutants were achieved with ionic strength of 250 mM pointed out with *. Standard deviation bars of three experiments are shown.
Figure 3. Effect of the addition of various concentrations of K+ (or Rb+) and (CH3)4N+ on the activities of WT-RMPK (A), E117K (B), T113L/E117K (C) and T113L/K114Q/E117K (D) mutants. Ionic strength and activities with 90 mM K+ (A) and without monovalent cations (BD) were those indicated in the legend of Figure 2. The normalized control activities were shown in white columns in (AD). White striped columns were in the presence of 25, 50, 75, and 100 mM (CH3)4N+ concentrations without K+ (or Rb+). In addition to the indicated (CH3)4N+ concentration, light gray, gray, dark gray, and black striped columns were in the presence of the 25, 50, 75, and 100 mM K+ (or Rb+) concentrations, respectively. The highest activities of RMPK mutants were achieved with ionic strength of 250 mM pointed out with *. Standard deviation bars of three experiments are shown.
Ijms 23 01347 g003
Figure 4. Effect of monovalent cations on the activities of WT-RMPK (A), K114Q (inset of A), E117K (B), T113L/E117K (C), and T113L/K114Q/E117K (D) mutants. The activities of WT-RMPK (A) and K114Q mutant (inset of A) and of the normalized activities of E117K (B), T113L/E117K (C) and T113L/K114Q/E117K (D) mutants in the presence of the indicated concentrations of Li+ (), Na+ (▲), K+ (■), NH4+ (), Rb+ (▼) or Cs+ () are shown. In the absence of monovalent cations and ionic strength of 250 mM, the activities of (BD) were 36 ± 6, 51 ± 1 and 82 ± 7, respectively. The experimental conditions were as in Figure 2, except that the ionic strength was maintained constant at 250 mM with (CH3)4N+. The reaction mixture for K114Q mutant, contained 1.5 mM PEP3−, 5.4 mM Mg-ADP and 3 mM Mg2+free and 135 mM of ionic strength. Standard deviation bars of three experiments are shown.
Figure 4. Effect of monovalent cations on the activities of WT-RMPK (A), K114Q (inset of A), E117K (B), T113L/E117K (C), and T113L/K114Q/E117K (D) mutants. The activities of WT-RMPK (A) and K114Q mutant (inset of A) and of the normalized activities of E117K (B), T113L/E117K (C) and T113L/K114Q/E117K (D) mutants in the presence of the indicated concentrations of Li+ (), Na+ (▲), K+ (■), NH4+ (), Rb+ (▼) or Cs+ () are shown. In the absence of monovalent cations and ionic strength of 250 mM, the activities of (BD) were 36 ± 6, 51 ± 1 and 82 ± 7, respectively. The experimental conditions were as in Figure 2, except that the ionic strength was maintained constant at 250 mM with (CH3)4N+. The reaction mixture for K114Q mutant, contained 1.5 mM PEP3−, 5.4 mM Mg-ADP and 3 mM Mg2+free and 135 mM of ionic strength. Standard deviation bars of three experiments are shown.
Ijms 23 01347 g004
Figure 5. Far-UV CD spectra of WT-RMPK (○) and of the RMPK mutants: E117K (□), K114Q (△), T113L/E117K (●), T113L/K114Q/E117K (■), T120L (▽), K114Q/E117K/T120L (▲) and T113L/K114Q/E117K/T120L (▼). The spectra were obtained in mixtures that contained each pyruvate kinase at a concentration of 100 μg/mL in 25 mM phosphate buffer pH 7.0 at 25 °C in cell pathway of 0.1 cm. CD is expressed as θMRW (mean residual weight ellipticity).
Figure 5. Far-UV CD spectra of WT-RMPK (○) and of the RMPK mutants: E117K (□), K114Q (△), T113L/E117K (●), T113L/K114Q/E117K (■), T120L (▽), K114Q/E117K/T120L (▲) and T113L/K114Q/E117K/T120L (▼). The spectra were obtained in mixtures that contained each pyruvate kinase at a concentration of 100 μg/mL in 25 mM phosphate buffer pH 7.0 at 25 °C in cell pathway of 0.1 cm. CD is expressed as θMRW (mean residual weight ellipticity).
Ijms 23 01347 g005
Figure 6. Ribbon representation of open and closed conformations of WT-RMPK (A) and of MtbPK (B). According to the movement of the upper B domain over the central A domain of WT-RMPK, panels show the superposition of open/inactive (cyan) and closed/active (green) conformers. In both conformers, the distances between O of T120 to O S76 or O of L75 to S37 are shown. Panel A was constructed from the coordinates of PDB 1A49 subunit A (green) with its active site completely occupied (oxalate, Mg2+, -Mg-ATP, and K+) (closed conformation) and subunit B (cyan) with its active site partially occupied (oxalate, Mg2+and K+) (open conformation) [17]. Panel B was constructed from the coordinates of PDB 5WS9 subunit A (green) with its active site completely occupied (oxalate, Mg2+ Mg-ATP, and K+) and of PDB 5WRP subunit A (cyan) with its active site empty [18]. This figure was built using the program PYMOL v 2.3.1.
Figure 6. Ribbon representation of open and closed conformations of WT-RMPK (A) and of MtbPK (B). According to the movement of the upper B domain over the central A domain of WT-RMPK, panels show the superposition of open/inactive (cyan) and closed/active (green) conformers. In both conformers, the distances between O of T120 to O S76 or O of L75 to S37 are shown. Panel A was constructed from the coordinates of PDB 1A49 subunit A (green) with its active site completely occupied (oxalate, Mg2+, -Mg-ATP, and K+) (closed conformation) and subunit B (cyan) with its active site partially occupied (oxalate, Mg2+and K+) (open conformation) [17]. Panel B was constructed from the coordinates of PDB 5WS9 subunit A (green) with its active site completely occupied (oxalate, Mg2+ Mg-ATP, and K+) and of PDB 5WRP subunit A (cyan) with its active site empty [18]. This figure was built using the program PYMOL v 2.3.1.
Ijms 23 01347 g006
Figure 7. RMSD and RMSF from molecular trajectories of WT-RMPK (A,B), T120L mutant (C,D), and MtbPK (E,F). Initial structures were in the presence of Mg2+, PEP, Mg-ADP and K+ for the K+-dependent WT-RMPK and T120L mutant and without K+ for the K+-independent MtbPK. The presence or absence of ligands are indicated by different colors in (A,C,E), whereas red and blue lines indicated with and without ligands in (B,D,F). Triplicates of molecular dynamics were carried out and the analyzes were made with CPPTRAJ [32]. For the particular case of MtbPK, runs 1 and 2 overlapped, either with or without ligands (E).
Figure 7. RMSD and RMSF from molecular trajectories of WT-RMPK (A,B), T120L mutant (C,D), and MtbPK (E,F). Initial structures were in the presence of Mg2+, PEP, Mg-ADP and K+ for the K+-dependent WT-RMPK and T120L mutant and without K+ for the K+-independent MtbPK. The presence or absence of ligands are indicated by different colors in (A,C,E), whereas red and blue lines indicated with and without ligands in (B,D,F). Triplicates of molecular dynamics were carried out and the analyzes were made with CPPTRAJ [32]. For the particular case of MtbPK, runs 1 and 2 overlapped, either with or without ligands (E).
Ijms 23 01347 g007
Figure 8. Comparison of the interactions of ADP (magenta and yellow) in the active site of the WT-RMPK (green) and T120L mutant (cyan). In the WT-RMPK, NH2 of R119 is at H-bond distance (3.13 and 3.11 Å) of O3B and O1B of ADP, respectively. In contrast, in the T120L mutant, NH2 of R119 is at 2.92 and 3.16 Å of ND1 and CE1 of H77, respectively. In WT-RMPK, NZ of K206 is at H-bond distance (2.57 and 3.19 Å) of O3´ and O2´of ribose, respectively; whereas in T120L mutant, NZ of K206 is at H-bond distance (2.96 Å) of O1A of ADP. This figure was prepared using the molecular dynamic models of both enzymes in the initial structures described in Figure 7 from PDB 149A.
Figure 8. Comparison of the interactions of ADP (magenta and yellow) in the active site of the WT-RMPK (green) and T120L mutant (cyan). In the WT-RMPK, NH2 of R119 is at H-bond distance (3.13 and 3.11 Å) of O3B and O1B of ADP, respectively. In contrast, in the T120L mutant, NH2 of R119 is at 2.92 and 3.16 Å of ND1 and CE1 of H77, respectively. In WT-RMPK, NZ of K206 is at H-bond distance (2.57 and 3.19 Å) of O3´ and O2´of ribose, respectively; whereas in T120L mutant, NZ of K206 is at H-bond distance (2.96 Å) of O1A of ADP. This figure was prepared using the molecular dynamic models of both enzymes in the initial structures described in Figure 7 from PDB 149A.
Ijms 23 01347 g008
Figure 9. Comparison between residues K117 in E117K mutant (magenta), K72 in the apo (cyan), and holoenzyme (green) in MtbPK (A). Close up of residues L68 to L75 of MtbPK forming part of the short loop that connects A domain to B domain in the apo (cyan) and holoenzyme (green) on MtbPK (B). In (A), the superposition of the E117K mutant coordinates (PDB 7R6Y subunit A; magenta), the holo-MtbPK (PDB 5WS9 subunit A; green), and the apo-MtbPK (PDB 5WRP subunit A; cyan) (18) is shown. Residue K72 in MtbPK corresponds to residue 117 in E117K mutant. The holoenzyme (green) contained oxalate, Mg2+, Mg-ATP, and K+; however, the nucleotide was omitted. The oxalate and Mg2+ present in the RMPK mutant are not shown. In (B), the figure shows that apo (cyan) or holoenzyme residues overlap or have the same orientation in the loop except for K72. This figure was prepared using PYMOL v 2.3.1.
Figure 9. Comparison between residues K117 in E117K mutant (magenta), K72 in the apo (cyan), and holoenzyme (green) in MtbPK (A). Close up of residues L68 to L75 of MtbPK forming part of the short loop that connects A domain to B domain in the apo (cyan) and holoenzyme (green) on MtbPK (B). In (A), the superposition of the E117K mutant coordinates (PDB 7R6Y subunit A; magenta), the holo-MtbPK (PDB 5WS9 subunit A; green), and the apo-MtbPK (PDB 5WRP subunit A; cyan) (18) is shown. Residue K72 in MtbPK corresponds to residue 117 in E117K mutant. The holoenzyme (green) contained oxalate, Mg2+, Mg-ATP, and K+; however, the nucleotide was omitted. The oxalate and Mg2+ present in the RMPK mutant are not shown. In (B), the figure shows that apo (cyan) or holoenzyme residues overlap or have the same orientation in the loop except for K72. This figure was prepared using PYMOL v 2.3.1.
Ijms 23 01347 g009
Figure 10. Differential coordination of Mg-ATP in the active site of a K+-dependent PK (WT-RMPK) (A) and in a K+-independent PK (MtbPK) (B). Panel A and Panel B show the coordination of oxalate, Mg2+, Mg-ATP, and K+ in the active sites of RMPK (PDB 1A49 subunit A; green) and of MtbPK (PDB 5WS9 subunit A; cyan), respectively. No residues of the enzymes are shown.
Figure 10. Differential coordination of Mg-ATP in the active site of a K+-dependent PK (WT-RMPK) (A) and in a K+-independent PK (MtbPK) (B). Panel A and Panel B show the coordination of oxalate, Mg2+, Mg-ATP, and K+ in the active sites of RMPK (PDB 1A49 subunit A; green) and of MtbPK (PDB 5WS9 subunit A; cyan), respectively. No residues of the enzymes are shown.
Ijms 23 01347 g010
Figure 11. Ribbon representation of the active sites of MtbPK in the holo (green) and the apoenzyme (cyan); and MD. Superposition of the coordinates of the holo- (PDB 5WS9 subunit A; green) to the apoenzyme (PDB 5WRP subunit A; cyan) are shown. As observed in Figure 9, K72 flips out or into the active site in the apo- and holoenzyme, respectively. ε-amino group of K72 interacts with Mg2+, Oxalate, Mg-ATP, K+, residues of A domain and of the lid (B domain) via a water network. Numbering is according to MtbPK. The bottom part shows the structural models of the MD of MtbPK every 20 ns, observing the same displacement of the K72 (green sticks). Ligands are represented as follows: ATP (blue sticks) Oxalate (yellow sticks), Mg2+ (orange spheres) and K+ (magenta spheres).
Figure 11. Ribbon representation of the active sites of MtbPK in the holo (green) and the apoenzyme (cyan); and MD. Superposition of the coordinates of the holo- (PDB 5WS9 subunit A; green) to the apoenzyme (PDB 5WRP subunit A; cyan) are shown. As observed in Figure 9, K72 flips out or into the active site in the apo- and holoenzyme, respectively. ε-amino group of K72 interacts with Mg2+, Oxalate, Mg-ATP, K+, residues of A domain and of the lid (B domain) via a water network. Numbering is according to MtbPK. The bottom part shows the structural models of the MD of MtbPK every 20 ns, observing the same displacement of the K72 (green sticks). Ligands are represented as follows: ATP (blue sticks) Oxalate (yellow sticks), Mg2+ (orange spheres) and K+ (magenta spheres).
Ijms 23 01347 g011
Table 1. Kinetics of the activation of WT-RMPK and K114Q mutant by monovalent cations at saturating concentrations of PEP3−, Mg-ADP complex, and Mg2+free.The data from Figure 3A were fitted to the Michaelis–Menten equation (v = Vmax S/Km + S) and to the Hill equation (v = Vmax Sn/K0.5n+ Sn) (Origin version 7.0). The mean and standard deviation from three experiments are shown. kcat/aK values are expressed in log form.
Table 1. Kinetics of the activation of WT-RMPK and K114Q mutant by monovalent cations at saturating concentrations of PEP3−, Mg-ADP complex, and Mg2+free.The data from Figure 3A were fitted to the Michaelis–Menten equation (v = Vmax S/Km + S) and to the Hill equation (v = Vmax Sn/K0.5n+ Sn) (Origin version 7.0). The mean and standard deviation from three experiments are shown. kcat/aK values are expressed in log form.
M+.Ionic radius (pm)WT-RMPKK114Q
kcat (s−1)K (mM)nlog E (s−1 M−1)kcat (s−1)K (mM)nlog E (s−1 M−1)
Li+6916 ± 353 ± 25___0.47N.D.N.D.N.D.N.D.
Na+102138 ± 848 ± 41.7 ± 0.13.46N.D.N.D.N.D.N.D.
K+138857 ± 165.8 ± 0.21.5 ± 0.15.17970 ± 2411 ± 0.61.4 ± 0.14.95
NH4+147596 ± 242.5 ± 0.21.8 ± 0.25.38742 ± 253.5 ± 0.21.6 ± 0.15.33
Rb+148462 ± 205.4 ± 0.41.4 ± 0.24.93691 ± 247.0 ± 0.51.5 ± 0.14.99
Cs+170198 ± 1220 ± 31.5 ± 0.13.99223 ± 1315 ± 1.81.4 ± 0.24.17
aK represents the K0.5 and Km for the data fitted to the Hill and Michaelis–Menten equations, respectively.
Table 2. Kinetic constants for PEP3−, MgADP complex and Mg2+free of WT-RMPK, K114Q-RMPK, E117K-RMPK, T113L/E117K-RMPK and T113L/K114Q/E117K-RMPK. The experiments for PEP3− were carried out in the presence of saturating MgADP complex and saturating Mg2+free, whereas those for MgADP and Mg2+free were performed in the presence of saturating concentrations of the other substrates. The experimental conditions were as in Figure 2 except for the absence of monovalent cations for pyruvate kinases mutants with K117 (K+-independent enzymes) and the inclusion of (CH3)4NCl to maintain ionic strength at 250 mM. Kinetic data (not shown) were fitted to the Hill equation (v = Vmax Sn/K0.5n + Sn) and to the Michaelis–Menten equation (v = Vmax S/Km + S) (Origin version 7.0). The mean and standard deviation from three experiments are shown. kcat/aK values are expressed in log form. Data of aT113L and bT113L were taken from Table S1 and Table S2 of Supplementary material, respectively, [21].
Table 2. Kinetic constants for PEP3−, MgADP complex and Mg2+free of WT-RMPK, K114Q-RMPK, E117K-RMPK, T113L/E117K-RMPK and T113L/K114Q/E117K-RMPK. The experiments for PEP3− were carried out in the presence of saturating MgADP complex and saturating Mg2+free, whereas those for MgADP and Mg2+free were performed in the presence of saturating concentrations of the other substrates. The experimental conditions were as in Figure 2 except for the absence of monovalent cations for pyruvate kinases mutants with K117 (K+-independent enzymes) and the inclusion of (CH3)4NCl to maintain ionic strength at 250 mM. Kinetic data (not shown) were fitted to the Hill equation (v = Vmax Sn/K0.5n + Sn) and to the Michaelis–Menten equation (v = Vmax S/Km + S) (Origin version 7.0). The mean and standard deviation from three experiments are shown. kcat/aK values are expressed in log form. Data of aT113L and bT113L were taken from Table S1 and Table S2 of Supplementary material, respectively, [21].
RMPK & RMPK MutantsWith E117RMPK MutantsWith K117
kcat (s−1)K (mM)nLog kcat/K (s−1 M−1)kcat (s−1)K (mM)nLog kcat/K (s−1 M−1)
PEPPEP
WT703 ± 280.063 ± 0.008__7.05E117K249 ± 160.063 ± 0.00__6.60
aT113L725 ± 490.108 ± 0.018__6.83T113L/E117K280 ± 80.048 ± 0.002__6.77
K114Q805 ± 260.15 ± 0.014__6.74T113L/K114Q/E117K363 ± 40.067 ± 0.001.82 ± 0.176.73
MgADP MgADP
WT691 ± 160.23 ± 0.010__6.48E117K261 ± 120.360 ± 0.0351.23 ± 0.095.86
bT113L831 ± 1201.07 ± 0.26__5.89T113L/E117K363 ± 280.100 ± 0.030__6.56
K14Q690 ± 190.54 ± 0.034__6.10T113L/K114Q/E117K419 ± 80.640 ± 0.0301.82 ± 0.175.82
Mg2+free Mg2+free
WT865 ± 120.17 ± 0.0101.46 ± 0.16.71E117K284 ± 44.9 ± 0.3__4.76
K114QN.D.N.D.N.D.N.D.T113L/E117K458 ± 751.6 ± 0.51.39 ± 0.35.46
T113LN.D.N.D.N.D.N.D.T113L/K114Q/E117K371 ± 168.7 ± 0.71.44 ± 0.14.63
K represents the K0.5 and Km for the data fitted to the Hill and Michaelis–Menten equations, respectively.
Table 3. List of DNA templates and primer sequences used in PCR.
Table 3. List of DNA templates and primer sequences used in PCR.
NameTemplate DNAPrimer Sequence 5′→3′
K114QWT-RMPKFW_CTGGACACTCAGGGACCCGAG
Rv_CTCGGGTCCCTGAGTGTCCAG
T120LWT-RMPKFw_ATCCGGCTAGGCCTCATCAAGGGC
Rv_GCCCTTGATGAGGCCTAGCCGGAT
T113L/E117KE117K mutant 1Fw_CTGGACCTTAAGGGACCCAAGATCCGGACG
Rv_AGTCCGGATCTTGGGTCCCTTAAGGTCCAG
T113L/K114Q/E117KE117K mutant 1Fw_CTCTGGACCTGCAGGGACCC
Rv_GGGTCCCTGCAGGTCCAGAG
K114Q/E117K/T120LK114Q mutantFw_GGACCCAAGATCCGGCTAGGCCTCATCAAGGGC
Rv_GCCCTTGATGAGGCCTAGCCGGATCTTGGGTCC
T113L/K114K/E117K/
T120L
T113L/K114Q/E117K
mutant
External (pTRC99A)Fw_TAATCATCCGGCTCGTATAATGTG
Rv_GGCTGAAAATCTTCTCTCATCCGC
1 Reference [20].
Table 4. Data collection and refinement statistics.
Table 4. Data collection and refinement statistics.
Data Collection
Wavelength [Å]0.98
Space groupP212121
Unit cell dimensions a, b, c [Å]108.34, 121.83, 161.37
Resolution range [Å] a48.18–2.25 (2.37–2.25)
Rmerge (%) b0.093 (0.57)
I/σI17.5 (3.0)
Observations
Total
Unique
733102 (89732)
101836 (14698)
Completeness99.9 (100.0)
Redundancy7.2 (6.1)
CC1/20.99 (0.81)
Refinement Statistics
Rwork/Rfree (%) c0.20/0.25
No. of atoms:15009
Protein14307
Ligands60
Metal Ions4
Solvent638
B value (Å2) 34.0
RMSD bond lengths (Å)0.006
RMSD bond angles (°)0.898
Ramachandran favored (%)98.73
Ramachandran allowed (%)1.11
Ramachandran outliers (%)0.16
PDB code 7R6Y
a Statistics for the highest resolution shell is shown in parenthesis. b Rmerge = ΣhklΣi |Ii(hkl)–<I(hkl)> | /ΣhklΣi Ii(hkl) c Rfree was calculated from a subset of reflections (4.9%), which were not used during refinement.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ramírez-Silva, L.; Hernández-Alcántara, G.; Guerrero-Mendiola, C.; González-Andrade, M.; Rodríguez-Romero, A.; Rodríguez-Hernández, A.; Lugo-Munguía, A.; Gómez-Coronado, P.A.; Rodríguez-Méndez, C.; Vega-Segura, A. The K+-Dependent and -Independent Pyruvate Kinases Acquire the Active Conformation by Different Mechanisms. Int. J. Mol. Sci. 2022, 23, 1347. https://doi.org/10.3390/ijms23031347

AMA Style

Ramírez-Silva L, Hernández-Alcántara G, Guerrero-Mendiola C, González-Andrade M, Rodríguez-Romero A, Rodríguez-Hernández A, Lugo-Munguía A, Gómez-Coronado PA, Rodríguez-Méndez C, Vega-Segura A. The K+-Dependent and -Independent Pyruvate Kinases Acquire the Active Conformation by Different Mechanisms. International Journal of Molecular Sciences. 2022; 23(3):1347. https://doi.org/10.3390/ijms23031347

Chicago/Turabian Style

Ramírez-Silva, Leticia, Gloria Hernández-Alcántara, Carlos Guerrero-Mendiola, Martin González-Andrade, Adela Rodríguez-Romero, Annia Rodríguez-Hernández, Alan Lugo-Munguía, Paul A. Gómez-Coronado, Cristina Rodríguez-Méndez, and Alicia Vega-Segura. 2022. "The K+-Dependent and -Independent Pyruvate Kinases Acquire the Active Conformation by Different Mechanisms" International Journal of Molecular Sciences 23, no. 3: 1347. https://doi.org/10.3390/ijms23031347

APA Style

Ramírez-Silva, L., Hernández-Alcántara, G., Guerrero-Mendiola, C., González-Andrade, M., Rodríguez-Romero, A., Rodríguez-Hernández, A., Lugo-Munguía, A., Gómez-Coronado, P. A., Rodríguez-Méndez, C., & Vega-Segura, A. (2022). The K+-Dependent and -Independent Pyruvate Kinases Acquire the Active Conformation by Different Mechanisms. International Journal of Molecular Sciences, 23(3), 1347. https://doi.org/10.3390/ijms23031347

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop