The 2.4 Å structure of Zymomonas mobilis pyruvate kinase: Implications for stability and regulation.

Human liver pyruvate kinase (hlPYK) catalyzes the final step in glycolysis, the formation of pyruvate (PYR) and ATP from phosphoenolpyruvate (PEP) and ADP. Fructose 1,6-bisphosphate (FBP), a pathway intermediate of glycolysis, serves as an allosteric activator of hlPYK. Zymomonas mobilis pyruvate kinase (ZmPYK) performs the final step of the Entner-Doudoroff pathway, which is similar to glycolysis in that energy is harvested from glucose and pyruvate is generated. The Entner-Doudoroff pathway does not have FBP as a pathway intermediate, and ZmPYK is not allosterically activated. In this work, we solved the 2.4 Å X-ray crystallographic structure of ZmPYK. The protein is dimeric in solution as determined by gel filtration chromatography, but crystallizes as a tetramer. The buried surface area of the ZmPYK tetramerization interface is significantly smaller than that of hlPYK, and yet tetramerization using the standard interfaces from higher organisms provides an accessible low energy crystallization pathway. Interestingly, the ZmPYK structure showed a phosphate ion in the analogous location to the 6-phosphate binding site of FBP in hlPYK. Circular Dichroism (CD) was used to measure melting temperatures of hlPYK and ZmPYK in the absence and presence of substrates and effectors. The only significant difference was an additional phase of small amplitude for the ZmPYK melting curves. We conclude that the phosphate ion plays neither a structural or allosteric role in ZmPYK under the conditions tested. We hypothesize that ZmPYK does not have sufficient protein stability for activity to be tuned by allosteric effectors as described for rheostat positions in the allosteric homologues.

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.