Hydrogen sulfide coordinates glucose metabolism switch through destabilizing tetrameric pyruvate kinase M2.

Most cancer cells reprogram their glucose metabolic pathway from oxidative phosphorylation to aerobic glycolysis for energy production. By reducing enzyme activity of pyruvate kinase M2 (PKM2), cancer cells attain a greater fraction of glycolytic metabolites for macromolecule synthesis needed for rapid proliferation. Here we demonstrate that hydrogen sulfide (H2S) destabilizes the PKM2 tetramer into monomer/dimer through sulfhydration at cysteines, notably at C326, leading to reduced PKM2 enzyme activity and increased PKM2-mediated transcriptional activation. Blocking PKM2 sulfhydration at C326 through amino acid mutation stabilizes the PKM2 tetramer and crystal structure further revealing the tetramer organization of PKM2-C326S. The PKM2-C326S mutant in cancer cells rewires glucose metabolism to mitochondrial respiration, significantly inhibiting tumor growth. In this work, we demonstrate that PKM2 sulfhydration by H2S inactivates PKM2 activity to promote tumorigenesis and inhibiting this process could be a potential therapeutic approach for targeting cancer metabolism.

Structural Basis of Nucleotide Selectivity in Pyruvate Kinase.

Nucleoside triphosphates are indispensable in numerous biological processes, with enzymes involved in their biogenesis playing pivotal roles in cell proliferation. Pyruvate kinase (PYK), commonly regarded as the terminal glycolytic enzyme that generates ATP in tandem with pyruvate, is also capable of synthesizing a wide range of nucleoside triphosphates from their diphosphate precursors. Despite their substrate promiscuity, some PYKs show preference towards specific nucleotides, suggesting an underlying mechanism for differentiating nucleotide bases. However, the thorough characterization of this mechanism has been hindered by the paucity of nucleotide-bound PYK structures. Here, we present crystal structures of Streptococcus pneumoniae PYK in complex with four different nucleotides. These structures facilitate direct comparison of the protein-nucleotide interactions and offer structural insights into its pronounced selectivity for GTP synthesis. Notably, this selectivity is dependent on a sequence motif in the nucleotide recognition site that is widely present among prokaryotic PYKs, particularly in Firmicutes species. We show that pneumococcal cell growth is significantly impaired when expressing a PYK variant with compromised GTP and UTP synthesis activity, underscoring the importance of PYK in maintaining nucleotide homeostasis. Our findings collectively advance our understanding of PYK biochemistry and prokaryotic metabolism.

Rabbit muscle pyruvate kinase activators: Synthesis, molecular docking and theoretical studies of N-substituted sulfonamide derivatives.

Pyruvate kinase (PK) activators have potential therapeutic applications in diseases such as sickle cell anemia. In this study, N-Substituted sulfonamide derivatives of 1,4-dihydropyridines were synthesized and evaluated as PK activators in vitro and using molecular docking studies. The compounds were synthesized by reacting dicarbonyl compounds with ammonium acetate, 5-nitrobenzaldehyde, and alumina sulfuric acid (ASA), followed by reduction and sulfonylation. The structures of the compounds were analyzed using spectroscopic techniques. DFT calculations provided insights into the electronic properties. Molecular docking of the compounds into the active site of PK showed favorable binding interactions. ADME evaluation indicated suitable solubility, BBB permeation, and lack of CYP450 inhibition. Overall, this study demonstrates the potential of new hybrid 1,4-dihydropyridine substituted sulfonamides as PK activators for further development. According to AC50 values, the compound (DTSF7, 0.97µM) is about 100-fold higher affective than the clinically used sulfonamide compound (AC50 = 90µM) for PK.

Evaluation of a Commercial TIMS-Q-TOF Platform for Native Mass Spectrometry.

Mass-spectrometry based assays in structural biology studies measure either intact or digested proteins. Typically, different mass spectrometers are dedicated for such measurements: those optimized for rapid analysis of peptides or those designed for high molecular weight analysis. A commercial trapped ion mobility-quadrupole-time-of-flight (TIMS-Q-TOF) platform is widely utilized for proteomics and metabolomics, with ion mobility providing a separation dimension in addition to liquid chromatography. The ability to perform high-quality native mass spectrometry of protein complexes, however, remains largely uninvestigated. Here, we evaluate a commercial TIMS-Q-TOF platform for analyzing noncovalent protein complexes by utilizing the instrument's full range of ion mobility, MS, and MS/MS (both in-source activation and collision cell CID) capabilities. The TIMS analyzer is able to be tuned gently to yield collision cross sections of native-like complexes comparable to those previously reported on various instrument platforms. In-source activation and collision cell CID were robust for both small and large complexes. TIMS-CID was performed on protein complexes streptavidin (53 kDa), avidin (68 kDa), and cholera toxin B (CTB, 58 kDa). Complexes pyruvate kinase (237 kDa) and GroEL (801 kDa) were beyond the trapping capabilities of the commercial TIMS analyzer, but TOF mass spectra could be acquired. The presented results indicate that the commercial TIMS-Q-TOF platform can be used for both omics and native mass spectrometry applications; however, modifications to the commercial RF drivers for both the TIMS analyzer and quadrupole (currently limited to m/z 3000) are necessary to mobility analyze protein complexes greater than about 60 kDa.

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.

Acetylation and Phosphorylation Regulate the Role of Pyruvate Kinase as a Glycolytic Enzyme or a Protein Kinase in Lamb.

Protein post-translational modifications (PTMs) play an essential role in meat quality development. However, the effect of specific PTM sites on meat proteins has not been investigated yet. The characteristics of pyruvate kinase M (PKM) were found to exhibit a close correlation with final meat quality, and thus, serine 99 (S99) and lysine 137 (K137) in PKM were mutated to study their effect on PKM function. The structural and functional properties of five lamb PKM variants, including wild-type PKM (wtPKM), PKM_S99D (S99 phosphorylation), PKM_S99A (PKM S99 dephosphorylation), PKM_K137Q (PKM K137 acetylation), and PKM_K137R (PKM K137 deacetylation), were evaluated. The results showed that the secondary structure, tertiary structure, and polymer formation were affected among different PKM variants. In addition, the glycolytic activity of PKM_K137Q was decreased because of its weakened binding with phosphoenolpyruvate. In the PKM_K137R variant, the actin phosphorylation level exhibited a decrease, suggesting a low kinase activity of PKM_K137R. The results of molecular simulation showed a 42% reduction in the interface area between PKM_K137R and actin, in contrast to wtPKM and actin. These findings are significant for revealing the mechanism of how PTMs regulate PKM function and provide a theoretical foundation for the development of precise meat quality preservation technology.

Optimizing the conversion of phosphoenolpyruvate to lactate by enzymatic channeling with mixed nanoparticle display.

Co-assembling enzymes with nanoparticles (NPs) into nanoclusters allows them to access channeling, a highly efficient form of multienzyme catalysis. Using pyruvate kinase (PykA) and lactate dehydrogenase (LDH) to convert phosphoenolpyruvic acid to lactic acid with semiconductor quantum dots (QDs) confirms how enzyme cluster formation dictates the rate of coupled catalytic flux (kflux) across a series of differentially sized/shaped QDs and 2D nanoplatelets (NPLs). Enzyme kinetics and coupled flux were used to demonstrate that by mixing different NP systems into clusters, a >10× improvement in kflux is observed relative to free enzymes, which is also ≥2× greater than enhancement on individual NPs. Cluster formation was characterized with gel electrophoresis and transmission electron microscopy (TEM) imaging. The generalizability of this mixed-NP approach to improving flux is confirmed by application to a seven-enzyme system. This represents a powerful approach for accessing channeling with almost any choice of enzymes constituting a multienzyme cascade.

Achieving unprecedented stability in lyophilized recombinase polymerase amplification with thermostable pyruvate kinase from Thermotoga maritima.

Recombinase polymerase amplification (RPA) is an isothermal DNA amplification reaction at around 41 °C using recombinase (Rec), single-stranded DNA-binding protein (SSB), strand-displacing DNA polymerase (Pol), and an ATP-regenerating enzyme. Considering the onsite use of RPA reagents, lyophilized RPA reagents with long storage stability are highly desired. In this study, as one of the approaches to solve this problem, we attempted to use a thermostable pyruvate kinase (PK). PK gene was isolated from a thermophilic bacterium Thermotoga maritima (Tma-PK). Tma-PK was expressed in Escherichia coli and purified from the cells. Tma-PK exhibited higher thermostability than human PK. The purified Tma-PK preparation was applied to RPA as an ATP-regenerating enzyme. Liquid RPA reagent with Tma-PK exhibited the same performance as that with human PK. Lyophilized RPA reagent with Tma-PK exhibited higher performance than that with human PK. Combined with our previous results of RPA reagents of thermostable Pol from a thermophilic bacterium, Aeribacillus pallidus, the results in this study suggest that thermostable enzymes are preferable to mesophilic ones as a component in lyophilized RPA reagents.

PKM2 induces mitophagy through the AMPK-mTOR pathway promoting CSFV proliferation.

Viruses exploit the host cell's energy metabolism system to support their replication. Mitochondria, known as the powerhouse of the cell, play a critical role in regulating cell survival and virus replication. Our prior research indicated that the classical swine fever virus (CSFV) alters mitochondrial dynamics and triggers glycolytic metabolic reprogramming. However, the role and mechanism of PKM2, a key regulatory enzyme of glycolytic metabolism, in CSFV replication remain unclear. In this study, we discovered that CSFV enhances PKM2 expression and utilizes PKM2 to inhibit pyruvate production. Using an affinity purification coupled mass spectrometry system, we successfully identified PKM as a novel interaction partner of the CSFV non-structural protein NS4A. Furthermore, we validated the interaction between PKM2 and both CSFV NS4A and NS5A through co-immunoprecipitation and confocal analysis. PKM2 was found to promote the expression of both NS4A and NS5A. Moreover, we observed that PKM2 induces mitophagy by activating the AMPK-mTOR signaling pathway, thereby facilitating CSFV proliferation. In summary, our data reveal a novel mechanism whereby PKM2, a metabolic enzyme, promotes CSFV proliferation by inducing mitophagy. These findings offer a new avenue for developing antiviral strategies. IMPORTANCE: Viruses rely on the host cell's material-energy metabolic system for replication, inducing host metabolic disorders and subsequent immunosuppression-a major contributor to persistent viral infections. Classical swine fever virus (CSFV) is no exception. Classical swine fever is a severe acute infectious disease caused by CSFV, resulting in significant economic losses to the global pig industry. While the role of the metabolic enzyme PKM2 (pyruvate dehydrogenase) in the glycolytic pathway of tumor cells has been extensively studied, its involvement in viral infection remains relatively unknown. Our data unveil a new mechanism by which the metabolic enzyme PKM2 mediates CSFV infection, offering novel avenues for the development of antiviral strategies.

The role of PKM2 in cancer progression and its structural and biological basis.

Pyruvate kinase M2 (PKM2), a subtype of pyruvate kinase (PK), has been shown to play an important role in the development of cancer. It regulates the last step of glycolytic pathway. PKM2 has both pyruvate kinase and protein kinase activity, and the conversion of these two functions of PKM2 depends on the mutual change of dimer and tetramer. The dimerization of PKM2 can promote the proliferation and growth of tumor cells, so inhibiting the dimerization of PKM2 is essential to curing cancer. The aggregation of PKM2 is regulated by both endogenous and exogenous cofactors as well as post-translational modification (PTM). Although there are many studies on the different aggregation of PKM2 in the process of tumor development, there are few summaries in recent years. In this review, we first introduce the role of PKM2 in various biological processes of tumor growth. Then, we summarize the aggregation regulation mechanism of PKM2 by various endogenous cofactors such as Fructose-1, 6-diphosphate (FBP), various amino acids, and post-translational modification (PTMs). Finally, the related inhibitors and agonists of PKM2 are summarized to provide reference for regulating PKM2 aggregation in the treatment of cancer in the future.

PYK-SubstitutionOME: an integrated database containing allosteric coupling, ligand affinity and mutational, structural, pathological, bioinformatic and computational information about pyruvate kinase isozymes.

Interpreting changes in patient genomes, understanding how viruses evolve and engineering novel protein function all depend on accurately predicting the functional outcomes that arise from amino acid substitutions. To that end, the development of first-generation prediction algorithms was guided by historic experimental datasets. However, these datasets were heavily biased toward substitutions at positions that have not changed much throughout evolution (i.e. conserved). Although newer datasets include substitutions at positions that span a range of evolutionary conservation scores, these data are largely derived from assays that agglomerate multiple aspects of function. To facilitate predictions from the foundational chemical properties of proteins, large substitution databases with biochemical characterizations of function are needed. We report here a database derived from mutational, biochemical, bioinformatic, structural, pathological and computational studies of a highly studied protein family-pyruvate kinase (PYK). A centerpiece of this database is the biochemical characterization-including quantitative evaluation of allosteric regulation-of the changes that accompany substitutions at positions that sample the full conservation range observed in the PYK family. We have used these data to facilitate critical advances in the foundational studies of allosteric regulation and protein evolution and as rigorous benchmarks for testing protein predictions. We trust that the collected dataset will be useful for the broader scientific community in the further development of prediction algorithms. Database URL https://github.com/djparente/PYK-DB.

Roadmap to Pyruvate Kinase M2 Modulation - A Computational Chronicle.

Pyruvate kinase M2 (PKM2) has surfaced as a potential target for anti-cancer therapy. PKM2 is known to be overexpressed in the tumor cells and is a critical metabolic conduit in supplying the augmented bioenergetic demands of the recalcitrant cancer cells. The presence of PKM2 in structurally diverse tetrameric as well as dimeric forms has opened new avenues to design novel modulators. It is also a truism to state that drug discovery has advanced significantly from various computational techniques like molecular docking, virtual screening, molecular dynamics, and pharmacophore mapping. The present review focuses on the role of computational tools in exploring novel modulators of PKM2. The structural features of various isoforms of PKM2 have been discussed along with reported modulators. An extensive analysis of the structure-based and ligand- based in silico methods aimed at PKM2 modulation has been conducted with an in-depth review of the literature. The role of advanced tools like QSAR and quantum mechanics has been established with a brief discussion of future perspectives.

Identification of residues involved in allosteric signal transmission from amino acid binding site of pyruvate kinase muscle isoform 2.

PKM2 is a rate-limiting enzyme in the glycolytic process and is involved in regulating tumor proliferation. Several amino acids (AAs) such as Asn, Asp, Val, and Cys have been shown to bind to the AA binding pocket of PKM2 and modulate its oligomeric state, substrate binding affinity, and activity. Although previous studies have attributed that the main chain and side chain of bound AAs are responsible for initiating signal to regulate PKM2, the signal transduction pathway remains elusive. To identify the residues involved in signal transfer process, N70 and N75 located at two ends of a ß strand connecting the active site and AA binding pocket were altered. Biochemical studies of these variants with various AA ligands (Asn, Asp, Val, and Cys), illustrate that N70 and N75, along with ß1 connecting these residues are part of the signal transduction pathway between the AA binding pocket and the active site. The results demonstrate that mutation of N70 to D prevents the transfer of the inhibitory signal mediated by Val and Cys, whereas N75 to L alteration blocks the activating signal initiated by Asn and Asp. Taken together, this study confirms that N70 is one of the residues responsible for transmitting the inhibitory signal and N75 is involved in the activation signal flow.

Serendipitous Identification of a Covalent Activator of Liver Pyruvate Kinase.

Enzymes are effective biological catalysts that accelerate almost all metabolic reactions in living organisms. Synthetic modulators of enzymes are useful tools for the study of enzymatic reactions and can provide starting points for the design of new drugs. Here, we report on the discovery of a class of biologically active compounds that covalently modifies lysine residues in human liver pyruvate kinase (PKL), leading to allosteric activation of the enzyme (EC50 =0.29 µM). Surprisingly, the allosteric activation control point resides on the lysine residue K282 present in the catalytic site of PKL. These findings were confirmed by structural data, MS/MS experiments, and molecular modelling studies. Altogether, our study provides a molecular basis for the activation mechanism and establishes a framework for further development of human liver pyruvate kinase covalent activators.

Identification of Natural Compounds as Inhibitors of Pyruvate Kinase M2 for Cancer Treatment.

The reliance of tumor cells on aerobic glycolysis is one of the emerging hallmarks of cancer. Pyruvate kinase M2 (PKM2), an important enzyme of glycolytic pathway, is highly expressed in a number of cancer cells. Tumor cells heavily depend on PKM2 to fulfill their divergent energetic and biosynthetic requirements, suggesting it as novel drug target for cancer therapies. Based on this context, we performed enzymatic-assay-based screening of the in-house phenolic compounds library for the identification of PKM2 inhibitors. This screening identified silibinin, curcumin, resveratrol, and ellagic acid as potential inhibitors of PKM2 with IC50 values of 0.91 µM, 1.12 µM, 3.07 µM, and 4.20 µM respectively. For the determination of Ki constants and the inhibition type of hit compounds, Lineweaver-Burk graphs were plotted. Silibinin and ellagic acid performed the competitive inhibition of PKM2 with Ki constants of 0.61 µM and 5.06 µM, while curcumin and resveratrol were identified as non-competitive inhibitors of PKM2 with Ki constants of 1.20 µM and 7.34 µM. The in silico screening of phenolic compounds against three binding sites of PKM2 provided insight into the binding pattern and functionally important amino residues of PKM2. Further, the evaluation of cytotoxicity via MTT assay demonstrated ellagic acid as potent inhibitor of cancer cell growth (IC50 = 20 µM). These results present ellagic acid, silibinin, curcumin, and resveratrol as inhibitors of PKM2 to interrogate metabolic reprogramming in cancer cells. This study has also provided the foundation for further research to validate the potential of identified bioactive entities for PKM2 targeted-cancer therapies.

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.

A magneto-controlled biocatalytic cascade with a fluorescent output.

A biocatalytic cascade based on concerted operation of pyruvate kinase and luciferase with a bioluminescent output was switched reversibly between low and high activity by applying an external magnetic field at different positions or removing it. The enzymes participating in the reaction cascade were bound to magnetic nanoparticles to allow their translocation or aggregation/dispersion to be controlled by the magnetic field. The reaction intensity, measured as the bioluminescent output, was dependent on the effective distances between the enzymes transported on the magnetic nanoparticles controlled by the magnets.

Reversible amyloids of pyruvate kinase couple cell metabolism and stress granule disassembly.

Cells respond to stress by blocking translation, rewiring metabolism and forming transient messenger ribonucleoprotein assemblies called stress granules (SGs). After stress release, re-establishing homeostasis and disassembling SGs requires ATP-consuming processes. However, the molecular mechanisms whereby cells restore ATP production and disassemble SGs after stress remain poorly understood. Here we show that upon stress, the ATP-producing enzyme Cdc19 forms inactive amyloids, and that their rapid re-solubilization is essential to restore ATP production and disassemble SGs in glucose-containing media. Cdc19 re-solubilization is initiated by the glycolytic metabolite fructose-1,6-bisphosphate, which directly binds Cdc19 amyloids, allowing Hsp104 and Ssa2 chaperone recruitment and aggregate re-solubilization. Fructose-1,6-bisphosphate then promotes Cdc19 tetramerization, which boosts its activity to further enhance ATP production and SG disassembly. Together, these results describe a molecular mechanism that is critical for stress recovery and directly couples cellular metabolism with SG dynamics via the regulation of reversible Cdc19 amyloids.

In Silico Study of the Structure and Ligand Preference of Pyruvate Kinases from Cyanobacterium Synechocystis sp. PCC 6803.

Finding reliable cheap sources for producing chemicals and materials is always challenging. During recent decades, photosynthetic organisms such as cyanobacteria, which used CO2 as a carbon source for making products, have attracted a great deal of attention. Among cyanobacteria, Synechocystis sp. PCC 6803 has been considered as a model strain and has some desirable features that make it suitable for use as an industrial strain. Pyruvate kinase (PK) catalyzes the transformation of phosphoenolpyruvate (PEP) to pyruvate in the last step of glycolysis that is an essential enzyme to produce adenosine triphosphate (ATP) in all organisms. Therefore, it plays a critical role in regulating cell metabolism. However, active and allosteric sites of PK and allosteric mechanisms governing PK activity are poorly understood in many bacteria. This study was aimed to provide more insight into PKs of Synechocystis sp. PCC 6803, using in silico methods. The results indicated that predicted structures of PKs from Synechocystis sp. PCC 6803 are reliable and can be considered for further studies. Molecular docking studies suggested that for predicted structures of sll0587 and sll1275, respectively, there are three and two possible active or allosteric sites. Furthermore, molecular interaction analysis of modeled structures proposes that sll0587 is strongly inhibited by ATP and when ATP concentration is low, this isoenzyme is active.

Rheostat functional outcomes occur when substitutions are introduced at nonconserved positions that diverge with speciation.

When amino acids vary during evolution, the outcome can be functionally neutral or biologically-important. We previously found that substituting a subset of nonconserved positions, "rheostat" positions, can have surprising effects on protein function. Since changes at rheostat positions can facilitate functional evolution or cause disease, more examples are needed to understand their unique biophysical characteristics. Here, we explored whether "phylogenetic" patterns of change in multiple sequence alignments (such as positions with subfamily specific conservation) predict the locations of functional rheostat positions. To that end, we experimentally tested eight phylogenetic positions in human liver pyruvate kinase (hLPYK), using 10-15 substitutions per position and biochemical assays that yielded five functional parameters. Five positions were strongly rheostatic and three were non-neutral. To test the corollary that positions with low phylogenetic scores were not rheostat positions, we combined these phylogenetic positions with previously-identified hLPYK rheostat, "toggle" (most substitution abolished function), and "neutral" (all substitutions were like wild-type) positions. Despite representing 428 variants, this set of 33 positions was poorly statistically powered. Thus, we turned to the in vivo phenotypic dataset for E. coli lactose repressor protein (LacI), which comprised 12-13 substitutions at 329 positions and could be used to identify rheostat, toggle, and neutral positions. Combined hLPYK and LacI results show that positions with strong phylogenetic patterns of change are more likely to exhibit rheostat substitution outcomes than neutral or toggle outcomes. Furthermore, phylogenetic patterns were more successful at identifying rheostat positions than were co-evolutionary or eigenvector centrality measures of evolutionary change.