Inhibition of Pyruvate Kinase From Thermoanaerobacterium saccharolyticum by IMP Is Independent of the Extra-C Domain

Are Allosteric Mechanisms Conserved Among Homologues?

Conservation of allosteric mechanisms among homologues is often assumed but seldom tested. This assumption underpins key concepts like coevolution of residues involved in allosteric mechanisms and the comparison of structures of two different homologues to gain insights into allosteric mechanisms. As an initial assessment of whether allosteric mechanisms are conserved among homologues, this work reviews what is known about the allosteric mechanisms of liver pyruvate kinase (LPYK) vs. skeletal muscle pyruvate kinase (M1PYK), framed within a two-ligand allosteric energy cycle description of allosteric regulation. Selective observations from other PYK homologues are included when relevant. The primary focus of this review is on functional data, while expressing caution regarding the interpretation of allosteric mechanisms based solely on available X-ray crystallographic structures. Additionally, this review considers types of data that are currently lacking for these two PYK homologues, highlighting potential techniques that could be valuable for evaluating the conservation of allosteric mechanisms among homologues. In particular, a hybrid tetramer technique that has been used to study bacterial phosphofructokinases 1 is summarized. Interestingly, despite a high degree of similarity (66.5% sequence identity) between the LPYK and rM1PYK proteins, the available functional comparisons do not provide strong evidence for conserved allosteric mechanisms. Lastly, we consider whether insights into native allosteric mechanisms are relevant to allosteric mechanisms associated with allosteric drug designs.

Rheostatic contributions to protein stability can obscure a position's functional role

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a "stability rheostat." The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

Divalent cations in human liver pyruvate kinase exemplify the combined effects of complex-equilibrium and allosteric regulation

There is growing recognition that the functional outcome of binding of an allosteric regulator to a protein/enzyme is influenced by the presence of other ligands. Here, this complexity is exemplified in the allosteric regulation of human liver pyruvate kinase (hLPYK) that is influenced by the presence of a range of divalent cation types and concentrations. For this system, fructose-1,6-bisphosphate (activator) and alanine (inhibitor) both influence the protein's affinity for the substrate, phosphoenolpyruvate (PEP). Mg2+, Mn2+, Ni2+, and Co2+ were the primary divalent cations evaluated, although Zn2+, Cd2+, V2+, Pb2+, Fe2+, and Cu2+also supported activity. Allosteric coupling between Fru-1,6-BP and PEP and between Ala and PEP varied depending on divalent cation type and concentration. Due to complicating interactions among small molecules, we did not attempt the fitting of response trends and instead we discuss a range of potential mechanisms that may explain those observed trends. Specifically, observed "substrate inhibition" may result from substrate A in one active site acting as an allosteric regulator for the affinity for substrate B in a second active site of a multimer. We also discuss apparent changes in allosteric coupling that can result from a sub-saturating concentration of a third allosteric ligand.

Odd one out? Functional tuning of Zymomonas mobilis pyruvate kinase is narrower than its allosteric, human counterpart

Various protein properties are often illuminated using sequence comparisons of protein homologs. For example, in analyses of the pyruvate kinase multiple sequence alignment, the set of positions that changed during speciation ("phylogenetic" positions) were enriched for "rheostat" positions in human liver pyruvate kinase (hLPYK). (Rheostat positions are those which, when substituted with various amino acids, yield a range of functional outcomes). However, the correlation was moderate, which could result from multiple biophysical constraints acting on the same position during evolution and/or various sources of noise. To further examine this correlation, we here tested Zymomonas mobilis PYK (ZmPYK), which has <65% sequence identity to any other PYK sequence. Twenty-six ZmPYK positions were selected based on their phylogenetic scores, substituted with multiple amino acids, and assessed for changes in Kapp-PEP . Although we expected to identify multiple, strong rheostat positions, only one moderate rheostat position was detected. Instead, nearly half of the 271 ZmPYK variants were inactive and most others showed near wild-type function. Indeed, for the active ZmPYK variants, the total range of Kapp,PEP values ("tunability") was 40-fold less than that observed for hLPYK variants. The combined functional studies and sequence comparisons suggest that ZmPYK has evolved functional and/or structural attributes that differ from the rest of the family. We hypothesize that including such "orphan" sequences in MSA analyses obscures the correlations used to predict rheostat positions. Finally, results raise the intriguing biophysical question as to how the same protein fold can support rheostat positions in one homolog but not another.