Energetic coupling between an oxidizable cysteine and the phosphorylatable N-terminus of human liver pyruvate kinase

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

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.

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.

Ellagic Acid and Its Metabolites as Potent and Selective Allosteric Inhibitors of Liver Pyruvate Kinase

Liver pyruvate kinase (PKL) has recently emerged as a new target for non-alcoholic fatty liver disease (NAFLD), and inhibitors of this enzyme could represent a new therapeutic option. However, this breakthrough is complicated by selectivity issues since pyruvate kinase exists in four different isoforms. In this work, we report that ellagic acid (EA) and its derivatives, present in numerous fruits and vegetables, can inhibit PKL potently and selectively. Several polyphenolic analogues of EA were synthesized and tested to identify the chemical features responsible for the desired activity. Molecular modelling studies suggested that this inhibition is related to the stabilization of the PKL inactive state. This unique inhibition mechanism could potentially herald the development of new therapeutics for NAFLD.

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 K_app-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 K_app,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.

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.

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

The pyruvate kinase (PYK) isozyme from Thermoanaerobacterium saccharolyticum (TsPYK) has previously been used in metabolic engineering for improved ethanol production. This isozyme belongs to a subclass of PYK isozymes that include an extra C-domain. Like other isozymes that include this extra C-domain, we found that TsPYK is activated by AMP and ribose-5-phosphate (R5P). Our use of sugar-phosphate analogs generated a surprising result in that IMP and GMP are allosteric inhibitors (rather than activators) of TsPYK. We believe this to be the first report of any PYK isozyme being inhibited by IMP and GMP. A truncated protein that lacks the extra C-domain is also inhibited by IMP. A screen of several other bacterial PYK enzymes (include several that have the extra-C domain) indicates that the inhibition by IMP is specific to only a subset of those isozymes.

Electrophile Modulation of Inflammation: A Two-Hit Approach

Electrophilic small molecules have gained significant attention over the last decade in the field of covalent drug discovery. Long recognized as mediators of the inflammatory process, recent evidence suggests that electrophiles may modulate the immune response through the regulation of metabolic networks. These molecules function as pleiotropic signaling mediators capable of reversibly reacting with nucleophilic biomolecules, most notably at reactive cysteines. More specifically, electrophiles target critical cysteines in redox regulatory proteins to activate protective pathways such as the nuclear factor erythroid 2-related factor 2-Kelch-like ECH-associated protein 1 (Nrf2-Keap1) antioxidant signaling pathway while also inhibiting Nuclear Factor κB (NF-κB). During inflammatory states, reactive species broadly alter cell signaling through the oxidation of lipids, amino acids, and nucleic acids, effectively propagating the inflammatory sequence. Subsequent changes in metabolic signaling inform immune cell maturation and effector function. Therapeutic strategies targeting inflammatory pathologies leverage electrophilic drug compounds, in part, because of their documented effect on the redox balance of the cell. With mounting evidence demonstrating the link between redox signaling and metabolism, electrophiles represent ideal therapeutic candidates for the treatment of inflammatory conditions. Through their pleiotropic signaling activity, electrophiles may be used strategically to both directly and indirectly target immune cell metabolism.

Distinct Hepatic PKA and CDK Signaling Pathways Control Activity-Independent Pyruvate Kinase Phosphorylation and Hepatic Glucose Production

Pyruvate kinase is an important enzyme in glycolysis and a key metabolic control point. We recently observed a pyruvate kinase liver isoform (PKL) phosphorylation site at S113 that correlates with insulin resistance in rats on a 3 day high-fat diet (HFD) and suggests additional control points for PKL activity. However, in contrast to the classical model of PKL regulation, neither authentically phosphorylated PKL at S12 nor S113 alone is sufficient to alter enzyme kinetics or structure. Instead, we show that cyclin-dependent kinases (CDKs) are activated by the HFD and responsible for PKL phosphorylation at position S113 in addition to other targets. These CDKs control PKL nuclear retention, alter cytosolic PKL activity, and ultimately influence glucose production. These results change our view of PKL regulation and highlight a previously unrecognized pathway of hepatic CDK activity and metabolic control points that may be important in insulin resistance and type 2 diabetes.

Identification of biochemically neutral positions in liver pyruvate kinase

Understanding how each residue position contributes to protein function has been a long-standing goal in protein science. Substitution studies have historically focused on conserved protein positions. However, substitutions of nonconserved positions can also modify function. Indeed, we recently identified nonconserved positions that have large substitution effects in human liver pyruvate kinase (hLPYK), including altered allosteric coupling. To facilitate a comparison of which characteristics determine when a nonconserved position does vs does not contribute to function, the goal of the current work was to identify neutral positions in hLPYK. However, existing hLPYK data showed that three features commonly associated with neutral positions-high sequence entropy, high surface exposure, and alanine scanning-lacked the sensitivity needed to guide experimental studies. We used multiple evolutionary patterns identified in a sequence alignment of the PYK family to identify which positions were least patterned, reasoning that these were most likely to be neutral. Nine positions were tested with a total of 117 amino acid substitutions. Although exploring all potential functions is not feasible for any protein, five parameters associated with substrate/effector affinities and allosteric coupling were measured for hLPYK variants. For each position, the aggregate functional outcomes of all variants were used to quantify a "neutrality" score. Three positions showed perfect neutral scores for all five parameters. Furthermore, the nine positions showed larger neutral scores than 17 positions located near allosteric binding sites. Thus, our strategy successfully enriched the dataset for positions with neutral and modest substitutions.

Rheostat positions: A new classification of protein positions relevant to pharmacogenomics

To achieve the full potential of pharmacogenomics, one must accurately predict the functional outcomes that arise from amino acid substitutions in proteins. Classically, researchers have focused on understanding the consequences of individual substitutions. However, literature surveys have shown that most substitutions were created at evolutionarily conserved positions. Awareness of this bias leads to a shift in perspective, from considering the outcomes of individual substitutions to understanding the roles of individual protein positions. Conserved positions tend to act as “toggle” switches, with most substitutions abolishing function. However, nonconserved positions have been found equally capable of affecting protein function. Indeed, many nonconserved positions act like functional dimmer switches (“rheostat” positions): this is revealed when multiple substitutions are made at a single position. Each substitution has a different functional outcome; the set of substitutions spans a range of outcomes. Finally, some nonconserved positions appear neutral, capable of accommodating all amino acid types without modifying function. This paper reviews the currently-known properties of rheostat positions, with examples shown for pyruvate kinase, organic anion transporting polypeptide 1B1, the beta-lactamase inhibitory protein, and angiotensin-converting enzyme 2. Outcomes observed for rheostat positions have implications for the rational design of drug analogs and allosteric drugs. Furthermore, this new framework—comprising three types of protein positions—provides a new approach to interpreting disease and population-based databases of amino acid changes. In conclusion, although a full understanding of substitution outcomes at rheostat positions poses a challenge, utilization of this new frame of reference will further advance the application of pharmacogenomics.

Structural and kinetic characterization of Trypanosoma congolense pyruvate kinase

Trypanosoma are blood-borne parasites and are the causative agents of neglected tropical diseases (NTDs) affecting both humans and animals. These parasites mainly rely on glycolysis for their energy production within the mammalian host, which is why trypanosomal glycolytic enzymes have been pursued as interesting targets for the development of trypanocidal drugs. The structure-function relationships of pyruvate kinases (PYKs) from trypanosomatids (Trypanosoma and Leishmania) have been well-studied within this context. In this paper, we describe the structural and enzymatic characterization of PYK from T. congolense (TcoPYK), the main causative agent of Animal African Trypanosomosis (AAT), by employing a combination of enzymatic assays, thermal unfolding studies and X-ray crystallography.

Functional tunability from a distance: Rheostat positions influence allosteric coupling between two distant binding sites

For protein mutagenesis, a common expectation is that important positions will behave like on/off “toggle” switches (i.e., a few substitutions act like wildtype, most abolish function). However, there exists another class of important positions that manifests a wide range of functional outcomes upon substitution: “rheostat” positions. Previously, we evaluated rheostat positions located near the allosteric binding sites for inhibitor alanine (Ala) and activator fructose-1,6-bisphosphate (Fru-1,6-BP) in human liver pyruvate kinase. When substituted with multiple amino acids, many positions demonstrated moderate rheostatic effects on allosteric coupling between effector binding and phosphoenolpyruvate (PEP) binding in the active site. Nonetheless, the combined outcomes of all positions sampled the full range of possible allosteric coupling (full tunability). However, that study only evaluated allosteric tunability of “local” positions, i.e., positions were located near the binding sites of the allosteric ligand being assessed. Here, we evaluated tunability of allosteric coupling when mutated sites were distant from the allosterically-coupled binding sites. Positions near the Ala binding site had rheostatic outcomes on allosteric coupling between Fru-1,6-BP and PEP binding. In contrast, positions in the Fru-1,6-BP site exhibited modest effects on coupling between Ala and PEP binding. Analyzed in aggregate, both PEP/Ala and PEP/Fru-1,6-BP coupling were again fully tunable by amino acid substitutions at this limited set of distant positions. Furthermore, some positions exhibited rheostatic control over multiple parameters and others exhibited rheostatic effects on one parameter and toggle control over a second. These findings highlight challenges in efforts to both predict/interpret mutational outcomes and engineer functions into proteins.

Mutational mimics of allosteric effectors: a genome editing design to validate allosteric drug targets

Development of drugs that allosterically regulate enzyme functions to treat disease is a costly venture. Amino acid susbstitutions that mimic allosteric effectors in vitro will identify therapeutic regulatory targets enhancing the likelihood of developing a disease treatment at a reasonable cost. We demonstrate the potential of this approach utilizing human liver pyruvate kinase (hLPYK) as a model. Inhibition of hLPYK was the first desired outcome of this study. We identified individual point mutations that: 1) mimicked allosteric inhibition by alanine, 2) mimicked inhibition by protein phosphorylation, and 3) prevented binding of fructose-1,6-bisphosphate (Fru-1,6-BP). Our second desired outcome was activation of hLPYK. We identified individual point mutations that: 1) prevented hLPYK from binding alanine, the allosteric inhibitor, 2) prevented inhibitory protein phosphorylation, or 3) mimicked allosteric activation by Fru-1,6-BP. Combining the three activating point mutations produced a constitutively activated enzyme that was unresponsive to regulators. Expression of a mutant hLPYK transgene containing these three mutations in a mouse model was not lethal. Thus, mutational mimics of allosteric effectors will be useful to confirm whether allosteric activation of hLPYK will control glycolytic flux in the diabetic liver to reduce hepatic glucose production and, in turn, reduce or prevent hyperglycemia.

Changes in the allosteric site of human liver pyruvate kinase upon activator binding include the breakage of an intersubunit cation-π bond

Human liver pyruvate kinase (hLPYK) converts phosphoenolpyruvate to pyruvate in the final step of glycolysis. hLPYK is allosterically activated by fructose-1,6-bisphosphate (Fru-1,6-BP). The allosteric site, as defined by previous structural studies, is located in domain C between the phosphate-binding loop (residues 444-449) and the allosteric loop (residues 527-533). In this study, the X-ray crystal structures of four hLPYK variants were solved to make structural correlations with existing functional data. The variants are D499N, W527H, Δ529/S531G (called GGG here) and S531E. The results revealed a conformational toggle between the open and closed positions of the allosteric loop. In the absence of Fru-1,6-BP the open position is stabilized, in part, by a cation-π bond between Trp527 and Arg538' (from an adjacent monomer). In the S531E variant glutamate binds in place of the 6'-phosphate of Fru-1,6-BP in the allosteric site, leading to partial allosteric activation. Finally, the structure of the D499N mutant does not provide structural evidence for the previously observed allosteric activation of the D499N variant.

Chokepoints in Mechanical Coupling Associated with Allosteric Proteins: The Pyruvate Kinase Example

Although the critical role of allostery in controlling enzymatic processes is well appreciated, there is a current dearth in our understanding of its underlying mechanisms, including communication between binding sites. One potential key aspect of intersite communication is the mechanical coupling between residues in a protein. Here, we introduce a graph-based computational approach to investigate the mechanical coupling between distant parts of a protein, highlighting effective pathways via which protein motion can transfer energy between sites. In this method, each residue is treated as a node on a weighted, undirected graph, in which the edges are defined by locally correlated motions of those residues and weighted by the strength of the correlation. The method was validated against experimental data on allosteric regulation in the human liver pyruvate kinase as obtained from full-protein alanine-scanning mutagenesis (systematic mutation) studies, as well as computational data on two G-protein-coupled receptors. The method provides semiquantitative information on the regulatory importance of specific structural elements. It is shown that these elements are key for the mechanical coupling between distant parts of the protein by providing effective pathways for energy transfer. It is also shown that, although there are a multitude of energy transfer pathways between distant parts of a protein, these pathways share a few common nodes that represent effective "chokepoints" for the communication.

Secondary and Supersecondary Structure of Proteins in Light of the Structure of Hydrophobic Cores

The traditional classification of protein structures (with regard to their supersecondary and tertiary conformation) is based on an assessment of conformational similarities between various polypeptide chains and particularly on the presence of specific secondary structural motifs. Mutual relations between secondary folds determine the overall shape of the protein and may be used to assign proteins to specific families (such as the immunoglobulin-like family). An alternative means of conducting structural assessment focuses on the structure of the protein's hydrophobic core. In this case, the protein is treated as a quasi-micelle, which exposes hydrophilic residues on its surface while internalizing hydrophobic residues. The accordance between the actual distribution of hydrophobicity in a protein and its corresponding theoretical ("idealized") distribution can be determined quantitatively, which, in turn, enables comparative analysis of structures regarded as geometrically similar (as well as geometrically divergent structures which are nevertheless regarded as similar in the sense of the fuzzy oil drop model). In this scope, the protein may be compared to an "intelligent micelle," where local disorder is often intentional and related to biological function-unlike traditional surfactant micelles which remain highly symmetrical throughout and do not carry any encoded information.

RheoScale: A tool to aggregate and quantify experimentally determined substitution outcomes for multiple variants at individual protein positions

Human mutations often cause amino acid changes (variants) that can alter protein function or stability. Some variants fall at protein positions that experimentally exhibit "rheostatic" mutation outcomes (different amino acid substitutions lead to a range of functional outcomes). In ongoing studies of rheostat positions, we encountered the need to aggregate experimental results from multiple variants, to describe the overall roles of individual positions. Here, we present "RheoScale" which generates quantitative scores to discriminate rheostat positions from those with "toggle" (most substitutions abolish function) or "neutral" (most substitutions have wild-type function) outcomes. RheoScale scores facilitate correlations of experimental data (such as binding affinity or stability) with structural and bioinformatic analyses. The RheoScale calculator is encoded into a Microsoft Excel workbook and an R script. Example analyses are shown for three model protein systems, including one assessed via deep mutational scanning. The RheoScale calculator quickly and efficiently provided quantitative descriptions that were in good agreement with prior qualitative observations. As an example application, scores were compared to the example proteins' structures; strong rheostat positions tended to occur in dynamic locations. In the future, RheoScale scores can be easily integrated into computational studies to facilitate improved algorithms for predicting outcomes of human variants.

Pathways crossing mammalian and plant sulfenomic landscapes

Reactive oxygen species (ROS) and especially hydrogen peroxide, are potent signaling molecules that activate cellular defense responses. Hydrogen peroxide can provoke reversible and irreversible oxidative posttranslational modifications on cysteine residues of proteins that act in diverse signaling circuits. The initial oxidation product of cysteine, sulfenic acid, has emerged as a biologically relevant posttranslational modification, because it is the primary sulfur oxygen modification that precedes divergent series of additional adaptations. In this review, we focus on the functional consequences of sulfenylation for both mammalian and plant proteins. Furthermore, we created compendia of sulfenylated proteins in human and plants based on mass spectrometry experiments, thereby defining the current plant and human sulfenomes. To assess the evolutionary conservation of sulfenylation, the sulfenomes of human and plants were compared based on protein homology. In total, 185 human sulfenylated proteins showed homology to sulfenylated plant proteins and the conserved sulfenylation targets participated in specific biological pathways and metabolic processes. Comprehensive functional studies of sulfenylation remains a future challenge, with multiple candidates suggested by mass spectrometry awaiting scrutinization.

The self-inhibitory nature of metabolic networks and its alleviation through compartmentalization

Metabolites can inhibit the enzymes that generate them. To explore the general nature of metabolic self-inhibition, we surveyed enzymological data accrued from a century of experimentation and generated a genome-scale enzyme-inhibition network. Enzyme inhibition is often driven by essential metabolites, affects the majority of biochemical processes, and is executed by a structured network whose topological organization is reflecting chemical similarities that exist between metabolites. Most inhibitory interactions are competitive, emerge in the close neighbourhood of the inhibited enzymes, and result from structural similarities between substrate and inhibitors. Structural constraints also explain one-third of allosteric inhibitors, a finding rationalized by crystallographic analysis of allosterically inhibited L-lactate dehydrogenase. Our findings suggest that the primary cause of metabolic enzyme inhibition is not the evolution of regulatory metabolite–enzyme interactions, but a finite structural diversity prevalent within the metabolome. In eukaryotes, compartmentalization minimizes inevitable enzyme inhibition and alleviates constraints that self-inhibition places on metabolism.

Metabolites act as enzyme inhibitors, but their global impact on metabolism has scarcely been considered. Here, the authors generate a human genome-wide metabolite-enzyme inhibition network, and find that inhibition occurs largely due to limited structural diversity of metabolites, leading to a global constraint on metabolism which subcellular compartmentalization minimizes.

Whole-protein alanine-scanning mutagenesis of allostery: A large percentage of a protein can contribute to mechanism

Many studies of allosteric mechanisms use limited numbers of mutations to test whether residues play "key" roles. However, if a large percentage of the protein contributes to allosteric function, mutating any residue would have a high probability of modifying allostery. Thus, a predicted mechanism that is dependent on only a few residues could erroneously appear to be supported. We used whole-protein alanine-scanning mutagenesis to determine which amino acid sidechains of human liver pyruvate kinase (hL-PYK; approved symbol PKLR) contribute to regulation by fructose-1,6-bisphosphate (Fru-1,6-BP; activator) and alanine (inhibitor). Each nonalanine/nonglycine residue of hL-PYK was mutated to alanine to generate 431 mutant proteins. Allosteric functions in active proteins were quantified by following substrate affinity over a concentration range of effectors. Results show that different residues contribute to the two allosteric functions. Only a small fraction of mutated residues perturbed inhibition by alanine. In contrast, a large percentage of mutated residues influenced activation by Fru-1,6-BP; inhibition by alanine is not simply the reverse of activation by Fru-1,6-BP. Moreover, the results show that Fru-1,6-BP activation would be extremely difficult to elucidate using a limited number of mutations. Additionally, this large mutational data set will be useful to train and test computational algorithms aiming to predict allosteric mechanisms.

Exploring the limits of the usefulness of mutagenesis in studies of allosteric mechanisms

The outcome of structure-guided mutational analyses is often used in support of postulated mechanisms of protein allostery. However, the limits of how informative mutations can be in understanding allosteric mechanisms are not completely clear. Here, we report an exercise to evaluate whether mutational data can support a simplistic mechanistic model, developed with minimal data inputs. Due to the lack of a mechanism to explain how alanine allosterically modifies the affinity of human liver pyruvate kinase (approved symbol PKLR) for its substrate, phosphoenolpyruvate, we proposed a speculative allosteric mechanism for this system. Within the allosteric amino-acid-binding site (something in the effector site must, of necessity, contribute to the allosteric mechanism), we implemented multiple mutational strategies: (1) site-directed random mutagenesis at positions that contact bound alanine and (2) mutations to probe specific questions. Despite acknowledged inadequacies used to formulate the speculative mechanism, many mutations modified the allosteric coupling constant (Q_ax ) consistent with that mechanism. The observed support for this speculative mechanism leaves us to ponder the best use of mutational data in structure-function studies of allosteric mechanisms. The mutational databank derived from this exercise has an independent value for training and testing algorithms specific to allostery.

Benchmarking predictions of allostery in liver pyruvate kinase in CAGI4

The Critical Assessment of Genome Interpretation (CAGI) is a global community experiment to objectively assess computational methods for predicting phenotypic impacts of genomic variation. One of the 2015-2016 competitions focused on predicting the influence of mutations on the allosteric regulation of human liver pyruvate kinase. More than 30 different researchers accessed the challenge data. However, only four groups accepted the challenge. Features used for predictions ranged from evolutionary constraints, mutant site locations relative to active and effector binding sites, and computational docking outputs. Despite the range of expertise and strategies used by predictors, the best predictions were marginally greater than random for modified allostery resulting from mutations. In contrast, several groups successfully predicted which mutations severely reduced enzymatic activity. Nonetheless, poor predictions of allostery stands in stark contrast to the impression left by more than 700 PubMed entries identified using the identifiers "computational + allosteric." This contrast highlights a specialized need for new computational tools and utilization of benchmarks that focus on allosteric regulation.

Are all regions of folded proteins that undergo ligand-dependent order-disorder transitions targets for allosteric peptide mimetics?

Although the classical view of how proteins function relied on well folded structures, it is now recognized that the functions of many proteins are dependent on being intrinsically disordered. The primary consideration in this work is the intermediate group of proteins that are overall well folded, but which contain small regions that undergo order-disorder transitions. In particular, the current focus is on those order-disorder transitions that are energetically coupled to ligand binding. As exemplified by the case of human liver pyruvate kinase (hL-PYK), peptides that mimic the sequence of the order-disorder region can be used as allosteric regulators of the enzyme. On the basis of this example and others reported in the literature, we propose that a similar use of peptides that mimic protein regions that experience ligand-dependent order-disorder transitions can be a generalized initiation point for the development of allosteric drugs.

1-Deoxynojirimycin Alleviates Liver Injury and Improves Hepatic Glucose Metabolism in db/db Mice

The present study investigated the effect of 1-Deoxynojirimycin (DNJ) on liver injury and hepatic glucose metabolism in db/db mice. Mice were divided into five groups: normal control, db/db control, DNJ-20 (DNJ 20 mg·kg⁻¹·day⁻¹), DNJ-40 (DNJ 40 mg·kg⁻¹·day⁻¹) and DNJ-80 (DNJ 80 mg·kg⁻¹·day⁻¹). All doses were treated intravenously by tail vein for four weeks. DNJ was observed to significantly reduce the levels of serum triglyceride (TG), total cholesterol (TC), low density lipoprotein cholesterol (LDL-C) and liver TG, as well as activities of serum alanine aminotransferase (ALT), and aspartate transaminase (AST); DNJ also alleviated macrovesicular steatosis and decreased tumor necrosis factor α (TNF-α), interleukin-1 (IL-1), interleukin-6 (IL-6) levels in liver tissue. Furthermore, DNJ treatment significantly increased hepatic glycogen content, the activities of hexokinase (HK), pyruvate kinase (PK) in liver tissue, and decreased the activities of glucose-6-phosphatase (G6Pase), glycogen phosphorylase (GP), and phosphoenolpyruvate carboxykinase (PEPCK). Moreover, DNJ increased the phosphorylation of phosphatidylinositol 3 kinase (PI3K) on p85, protein kinase B (PKB) on Ser473, glycogen synthase kinase 3β (GSK-3β) on Ser9, and inhibited phosphorylation of glycogen synthase (GS) on Ser645 in liver tissue of db/db mice. These results demonstrate that DNJ can increase hepatic insulin sensitivity via strengthening of the insulin-stimulated PKB/GSK-3β signal pathway and by modulating glucose metabolic enzymes in db/db mice. Moreover, DNJ also can improve lipid homeostasis and attenuate hepatic steatosis in db/db mice.

Distinguishing the interactions in the fructose 1,6-bisphosphate binding site of human liver pyruvate kinase that contribute to allostery

In the study of allosteric proteins, understanding which effector-protein interactions contribute to allosteric activation is important both for designing allosteric drugs and for understanding allosteric mechanisms. The antihyperglycemic target, human liver pyruvate kinase (hL-PYK), binds its allosteric activator, fructose 1,6-bisphosphate (Fru-1,6-BP), such that the 1'-phosphate interacts with side chains of Arg501 and Trp494 and the 6'-phosphate interacts with Thr444, Thr446, Ser449 (i.e., the 444-449 loop), and Ser531. Additionally, backbone atoms from the 527-533 loop interact with a sugar ring hydroxyl and the two effector phosphate moieties. An effector analogue series indicates that only one phosphate on the sugar is required for activation. However, singly phosphorylated sugars, including Fru-1-P and Fru-6-P, bind with a Kix in the range of 0.07-1 mM. The second phosphate of Fru-1,6-BP causes tight effector binding, because this native effector has a Kix of 0.061 μM. Glucose 1,6-bisphosphate and ribulose 1,5-bisphosphate bind in the 0.07-1 mM range. The contrast with a higher Fru-1,6-BP binding indicates specificity for the fructose sugar conformation. Site-directed random mutagenesis at each residue that contacts bound Fru-1,6-BP showed that a negative charge introduced at position 531 mimics allosteric activation, even in the absence of Fru-1,6-BP. Collectively, analogue and mutagenesis studies are consistent with the 527-533 loop playing a key role in allosteric function. Deletion mutations that shortened the 527-533 loop were expected to prevent formation of hydrogen bonds between backbone atoms on the loop and Fru-1,6-BP. Indeed, Fru-1,6-BP did not activate these loop-shortened mutant proteins. Previous structural comparisons of M1-PYK and M2-PYK indicate that the 527-533 loop makes interactions across a subunit interface when an activator is not present. Mutating the hL-PYK subunit interface interactions among Trp527, Arg528, and Asp499 mimics allosteric activation. Considered with published structures, these results are consistent with (1) the two phosphates of Fru-1,6-BP docking to Arg501/Trp494 and the 444-449 loop, respectively, and (2) the formation of hydrogen bonds among Fru-1,6-BP and backbone atoms of the 527-533 loop pulling this loop away from the subunit interface, which results in breaking of the Trp527-Arg528-Asp499 interactions to elicit an allosteric response.

The Nrf2 regulatory network provides an interface between redox and intermediary metabolism

Nuclear factor-erythroid 2 p45-related factor 2 (Nrf2, also called Nfe2l2) is a transcription factor that regulates the cellular redox status. Nrf2 is controlled through a complex transcriptional/epigenetic and post-translational network that ensures its activity increases during redox perturbation, inflammation, growth factor stimulation and nutrient/energy fluxes, thereby enabling the factor to orchestrate adaptive responses to diverse forms of stress. Besides mediating stress-stimulated induction of antioxidant and detoxification genes, Nrf2 contributes to adaptation by upregulating the repair and degradation of damaged macromolecules, and by modulating intermediary metabolism. In the latter case, Nrf2 inhibits lipogenesis, supports β-oxidation of fatty acids, facilitates flux through the pentose phosphate pathway, and increases NADPH regeneration and purine biosynthesis; these observations suggest Nrf2 directs metabolic reprogramming during stress.

Liraglutide ameliorates glycometabolism and insulin resistance through the upregulation of GLUT4 in diabetic KKAy mice

Liraglutide, a long-lasting glucagon‑like peptide‑1 analogue, has been used for the treatment of patients with type 2 diabetes mellitus since 2009. In this study, we investigated the anti-diabetic effects and mechanisms of action of liraglutide in a spontaneous diabetic animal model, using KK/Upj-Ay/J (KKAy) mice. The KKAy mice were divided into 2 groups, the liraglutide group (mice were treated with 250 µg/kg/day liraglutide) and the model group (treated with an equivalent amount of normal saline). C57BL/6J mice were used as the controls (treated with an equivalent amount of normal saline). The treatment period lasted 6 weeks. During this treatment period, fasting blood glucose (FBG) levels and the body weight of the mice were measured on a weekly basis. Our results revealed that liraglutide significantly decreased FBG levels, the area under the curve following a oral glucose tolerance test and insulin tolerance test, increased serum insulin levels, reduced homeostasis model assessment of insulin resistance and increased the insulin sensitivity index. Furthermore, liraglutide ameliorated glycometabolism dysfunction by increasing glycolysis via hexokinase and glycogenesis via pyruvate kinase activation. An ultrastructural examination of the pancreas revealed that liraglutide improved the damaged state of islet β cells and increased the number of insulin secretory granules. The real-time PCR results revealed that the gene expression of glucose transporter 4 (GLUT4) increased following treatment with liraglutide. Liraglutide also upregulated the protein expression of GLUT4 in liver tissue and skeletal muscle. Our results suggest that liraglutide ameliorates glycometabolism and insulin resistance in diabetic KKAy mice by stimulating insulin secretion, increasing glycogenesis and glycolysis and upregulating the expression of GLUT4.

The redox biochemistry of protein sulfenylation and sulfinylation

Controlled generation of reactive oxygen species orchestrates numerous physiological signaling events (Finkel, T. (2011) Signal transduction by reactive oxygen species. J. Cell Biol. 194, 7-15). A major cellular target of reactive oxygen species is the thiol side chain (RSH) of Cys, which may assume a wide range of oxidation states (i.e. -2 to +4). Within this context, Cys sulfenic (Cys-SOH) and sulfinic (Cys-SO2H) acids have emerged as important mechanisms for regulation of protein function. Although this area has been under investigation for over a decade, the scope and biological role of sulfenic/sulfinic acid modifications have been recently expanded with the introduction of new tools for monitoring cysteine oxidation in vitro and directly in cells. This minireview discusses selected recent examples of protein sulfenylation and sulfinylation from the literature, highlighting the role of these post-translational modifications in cell signaling.