Substitutions at a rheostat position in human aldolase A cause a shift in the conformational population

Substitutions at rheostat position 52 of LacI have long-range effects on the LacI conformational landscape

In proteins, amino acid changes at "rheostat" positions exhibit functional changes that vary with the substitution chosen: some substitutions enhance function, some are like wild-type, some are partially detrimental, while others abolish function. One way that substitutions might exert their complex effects is by altering protein conformational landscapes. To test this, we studied five substitutions of V52 in E. coli LacI, an experimentally-known rheostat position. For each variant, we mapped the accessible conformational landscapes by performing molecular dynamics simulations at ambient conditions and under three perturbations: increased pressure, binding to allosteric ligand "ONPF", and ONPF plus pressure. The simulated DNA binding domain landscapes were compared to published experimentally-measured parameters, and the results suggest that complex combinations of dynamic parameters and/or additional simulations in the presence of DNA are needed to predict DNA binding specificity. For the variants regulatory domains all landscapes displayed boundaries similar to wild-type, but changes within the boundaries were unique. Of these, V52A/ONPF was striking: The regulatory domains for ONPF-bound, wild-type LacI are in an "Open" conformation and, experimentally, ONPF enhances DNA binding. Four variants responded to ONPF like wild-type, but ONPF binding to V52A shifted these domains to a "Closed" conformation that is associated with diminished DNA binding for wild-type LacI. This finding predicted that ONPF's allosteric regulation of V52A would change from "anti-inducer" to "inducer", which we experimentally validated in vivo and in vitro. This supports the hypothesis that substituting rheostat positions can alter function by altering the relative populations on protein conformational landscapes.

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.

Identification of positions in human aldolase a that are neutral for apparent KM

According to evolutionary theory, many naturally-occurring amino acid substitutions are expected to be neutral or near-neutral, with little effect on protein structure or function. Accordingly, most changes observed in human exomes are also expected to be neutral. As such, accurate algorithms for identifying medically-relevant changes must discriminate rare, non-neutral substitutions against a background of neutral substitutions. However, due to historical biases in biochemical experiments, the data available to train and validate prediction algorithms mostly contains non-neutral substitutions, with few examples of neutral substitutions. Thus, available training sets have the opposite composition of the desired test sets. Towards improving a dataset of these critical negative controls, we have concentrated on identifying neutral positions - those positions for which most of the possible 19 amino acid substitutions have little effect on protein structure or function. Here, we used a strategy based on multiple sequence alignments to identify putative neutral positions in human aldolase A, followed by biochemical assays for 147 aldolase substitutions. Results showed that most variants had little effect on either the apparent Michaelis constant for substrate fructose-1,6-bisphosphate or its apparent cooperativity. Thus, these data are useful for training and validating prediction algorithms. In addition, we created a database of these and other biochemically characterized aldolase variants along with aldolase sequences and characteristics derived from sequence and structure analyses. This database is publicly available at https://github.com/liskinsk/Aldolase-variant-and-sequence-database.

Dynamics-based protein network features accurately discriminate neutral and rheostat positions

In some proteins, a unique class of nonconserved positions is characterized by their ability to generate diverse functional outcomes through single amino acid substitutions. Due to their ability to tune protein function, accurately identifying such "rheostat" positions is crucial for protein design, for understanding the impact of mutations observed in humans, and for predicting the evolution of pathogen drug resistance. However, identifying rheostat positions has been challenging, due-in part-to the absence of a clear structural relationship with binding sites. In this study, experimental data from our previous study of the Escherichia coli lactose repressor protein (LacI) was used to identify rheostat positions for which mutations tune in vivo EC50 for the allosteric ligand "IPTG." We next used the rheostat assignments to test the hypothesis that rheostat positions have unique dynamic features that will enable their identification. To that end, we integrated all-atom molecular dynamics simulations with perturbation residue response analysis. Results first revealed distinct dynamic behavior in IPTG-bound LacI compared with apo LacI, which was consistent with IPTG's role as an allosteric inducer. Next, we used a variety of dynamic features to build a classification model that discriminates experimentally characterized rheostat positions in LacI from positions with other types of substitution outcomes. In parallel, we built a second classifier model based on the 3D structural "static" network features of LacI. In comparative studies, the dynamic model better identified rheostat positions that were >8 Å from the binding site. In summary, our study provides insights into the dynamic characteristics of rheostat positions and suggests that models built on dynamic features may be useful for predicting the locations of rheostat positions in a wide range of proteins.

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.

Impact of single-residue mutations on protein thermal stability: The case of threonine 83 of BC2L-CN lectin

The thermal stability of trimeric lectin BC2L-CN was investigated and found to be considerably altered when mutating residue 83, originally a threonine, located at the fucose-binding loop. Mutants were analyzed using differential scanning calorimetry and isothermal microcalorimetry. Although most mutations decreased the affinity of the protein for oligosaccharide H type 1, six mutations increased the melting temperature (Tm) by >5 °C; one mutation, T83P, increased the Tm value by 18.2 °C(T83P, Tm = 96.3 °C). In molecular dynamic simulations, the investigated thermostable mutants, T83P, T83A, and T83S, had decreased fluctuations in the loop containing residue 83. In the T83S mutation, the side-chain hydroxyl group of serine formed a hydrogen bond with a nearby residue, suggesting that the restricted movement of the side-chain resulted in fewer fluctuations and enhanced thermal stability. Residue 83 is located at the interface and near the upstream end of the equivalent loop in a different protomer; therefore, fluctuations by this residue likely propagate throughout the loop. Our study of the dramatic change in thermal stability by a single amino acid mutation provides useful insights into the rational design of protein structures, especially the structures of oligomeric proteins.

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.

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.

De novo distillation of thermodynamic affinity from deep learning regulatory sequence models of in vivo protein-DNA binding

Transcription factors (TF) are proteins that bind DNA in a sequence-specific manner to regulate gene transcription. Despite their unique intrinsic sequence preferences, in vivo genomic occupancy profiles of TFs differ across cellular contexts. Hence, deciphering the sequence determinants of TF binding, both intrinsic and context-specific, is essential to understand gene regulation and the impact of regulatory, non-coding genetic variation. Biophysical models trained on in vitro TF binding assays can estimate intrinsic affinity landscapes and predict occupancy based on TF concentration and affinity. However, these models cannot adequately explain context-specific, in vivo binding profiles. Conversely, deep learning models, trained on in vivo TF binding assays, effectively predict and explain genomic occupancy profiles as a function of complex regulatory sequence syntax, albeit without a clear biophysical interpretation. To reconcile these complementary models of in vitro and in vivo TF binding, we developed Affinity Distillation (AD), a method that extracts thermodynamic affinities de-novo from deep learning models of TF chromatin immunoprecipitation (ChIP) experiments by marginalizing away the influence of genomic sequence context. Applied to neural networks modeling diverse classes of yeast and mammalian TFs, AD predicts energetic impacts of sequence variation within and surrounding motifs on TF binding as measured by diverse in vitro assays with superior dynamic range and accuracy compared to motif-based methods. Furthermore, AD can accurately discern affinities of TF paralogs. Our results highlight thermodynamic affinity as a key determinant of in vivo binding, suggest that deep learning models of in vivo binding implicitly learn high-resolution affinity landscapes, and show that these affinities can be successfully distilled using AD. This new biophysical interpretation of deep learning models enables high-throughput in silico experiments to explore the influence of sequence context and variation on both intrinsic affinity and in vivo occupancy.

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.

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.

Structural Plasticity Is a Feature of Rheostat Positions in the Human Na+/Taurocholate Cotransporting Polypeptide (NTCP)

In the Na⁺/taurocholate cotransporting polypeptide (NTCP), the clinically relevant S267F polymorphism occurs at a "rheostat position". That is, amino acid substitutions at this position ("S267X") lead to a wide range of functional outcomes. This result was particularly striking because molecular models predicted the S267X side chains are buried, and thus, usually expected to be less tolerant of substitutions. To assess whether structural tolerance to buried substitutions is widespread in NTCP, here we used Rosetta to model all 19 potential substitutions at another 13 buried positions. Again, only subtle changes in the calculated stabilities and structures were predicted. Calculations were experimentally validated for 19 variants at codon 271 ("N271X"). Results showed near wildtype expression and rheostatic modulation of substrate transport, implicating N271 as a rheostat position. Notably, each N271X substitution showed a similar effect on the transport of three different substrates and thus did not alter substrate specificity. This differs from S267X, which altered both transport kinetics and specificity. As both transport and specificity may change during protein evolution, the recognition of such rheostat positions may be important for evolutionary studies. We further propose that the presence of rheostat positions is facilitated by local plasticity within the protein structure. Finally, we note that identifying rheostat positions may advance efforts to predict new biomedically relevant missense variants in NTCP and other membrane transport proteins.