Search

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 (K app-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 K app-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 K app-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's K app-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., (© 2024 The Protein Society.)

In Escherichia coli, the master transcription regulator catabolite repressor activator (Cra) regulates >100 genes in central metabolism. Cra binding to DNA is allosterically regulated by binding to fructose-1-phosphate (F-1-P), but the only documented source of F-1-P is from the concurrent import and phosphorylation of exogenous fructose. Thus, many have proposed that fructose-1,6-bisphosphate (F-1,6-BP) is also a physiological regulatory ligand. However, the role of F-1,6-BP has been widely debated. Here, we report that the E. coli enzyme fructose-1-kinase (FruK) can carry out its "reverse" reaction under physiological substrate concentrations to generate F-1-P from F-1,6-BP. We further show that FruK directly binds Cra with nanomolar affinity and forms higher order, heterocomplexes. Growth assays with a ΔfruK strain and fruK complementation show that FruK has a broader role in metabolism than fructose catabolism. Since fruK itself is repressed by Cra, these newly-reported events add layers to the dynamic regulation of E. coli's central metabolism that occur in response to changing nutrients. These findings might have wide-spread relevance to other γ-proteobacteria, which conserve both Cra and FruK., Competing Interests: Conflict of interests During the course of this work, Dr Bose served on the Scientific Advisory Board and was a consultant for Azitra, Inc and Merck & Co, Inc. These activities did not financially support and are unrelated to the current manuscript. The other authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2023 The Protein Society.)

The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions., Competing Interests: Declaration of Competing Interest The authors declare no competing interests., (Copyright © 2023. Published by Elsevier B.V.)

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., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and β-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + β-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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., (© The Author(s) 2023. Published by Oxford University Press.)

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., (© 2022 The Protein Society.)

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.

Some protein positions play special roles in determining the magnitude of protein function: at such "rheostat" positions, varied amino acid substitutions give rise to a continuum of functional outcomes, from wild type (or enhanced), to intermediate, to loss of function. This observed range raises interesting questions about the biophysical bases by which changes at single positions have such varied outcomes. Here, we assessed variants at position 98 in human aldolase A ("I98X"). Despite being ~17 Å from the active site and far from subunit interfaces, substitutions at position 98 have rheostatic contributions to the apparent cooperativity (n H ) associated with fructose-1,6-bisphosphate substrate binding and moderately affected binding affinity. Next, we crystallized representative I98X variants to assess structural consequences. Residues smaller than the native isoleucine (cysteine and serine) were readily accommodated, and the larger phenylalanine caused only a slight separation of the two parallel helixes. However, the diffraction quality was reduced for I98F, and further reduced for I98Y. Intriguingly, the resolutions of the I98X structures correlated with their n H values. We propose that substitution effects on both n H and crystal lattice disruption arise from changes in the population of aldolase A conformations in solution. In combination with results computed for rheostat positions in other proteins, the results from this study suggest that rheostat positions accommodate a wide range of side chains and that structural consequences manifest as shifted ensemble populations and/or dynamics changes., (© 2021 The Protein Society.)

To create bacterial transcription "circuits" for biotechnology, one approach is to recombine natural transcription factors, promoters, and operators. Additional novel functions can be engineered from existing transcription factors such as the E. coli AraC transcriptional activator, for which binding to DNA is modulated by binding L-arabinose. Here, we engineered chimeric AraC/XylS transcription activators that recognized ara DNA binding sites and responded to varied effector ligands. The first step, identifying domain boundaries in the natural homologs, was challenging because (i) no full-length, dimeric structures were available and (ii) extremely low sequence identities (≤10%) among homologs precluded traditional assemblies of sequence alignments. Thus, to identify domains, we built and aligned structural models of the natural proteins. The designed chimeric activators were assessed for function, which was then further improved by random mutagenesis. Several mutational variants were identified for an XylS•AraC chimera that responded to benzoate; two enhanced activation to near that of wild-type AraC. For an RhaR•AraC chimera, a variant with five additional substitutions enabled transcriptional activation in response to rhamnose. These five changes were dispersed across the protein structure, and combinatorial experiments testing subsets of substitutions showed significant non-additivity. Combined, the structure modeling and epistasis suggest that the common AraC/XylS structural scaffold is highly interconnected, with complex intra-protein and inter-domain communication pathways enabling allosteric regulation. At the same time, the observed epistasis and the low sequence identities of the natural homologs suggest that the structural scaffold and function of transcriptional regulation are nevertheless highly accommodating of amino acid changes., (© 2021 Wiley Periodicals LLC.)

Tetanus neurotoxin (TeNT) is an A-B toxin with three functional domains: endopeptidase, translocation (HCT), and receptor binding. Endosomal acidification triggers HCT to interact with and insert into the membrane, translocating the endopeptidase across the bilayer. Although the function of HCT is well defined, the mechanism by which it accomplishes this task is unknown. To gain insight into the HCT membrane interaction on both local and global scales, we utilized an isolated, beltless HCT variant (bHCT), which retained the ability to release potassium ions from vesicles. To examine which bHCT residues interact with the membrane, we widely sampled the surface of bHCT using 47 single-cysteine variants labeled with the environmentally sensitive fluorophore NBD. At neutral pH, no interaction was observed for any variant. In contrast, all NBD-labeled positions reported environmental change in the presence of acidic pH and membranes containing anionic lipids. We then examined the conformation of inserted bHCT using circular dichroism and intrinsic fluorescence. Upon entering the membrane, bHCT retained predominantly α-helical secondary structure, whereas the tertiary structure exhibited substantial refolding. The use of lipid-attached quenchers revealed that at least one of the three tryptophan residues penetrated deep into the hydrocarbon core of the membrane, suggesting formation of a bHCT transmembrane conformation. The possible conformational topology was further explored with the hydropathy analysis webtool MPEx, which identified a large, potential α-helical transmembrane region. Altogether, the spectroscopic evidence supports a model in which, upon acidification, the majority of TeNT bHCT entered the membrane with a concurrent change in tertiary structure., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (© 2021 The Protein Society.)

Amino acid substitutions at nonconserved protein positions can have noncanonical and "long-distance" outcomes on protein function. Such outcomes might arise from changes in the internal protein communication network, which is often accompanied by changes in structural flexibility. To test this, we calculated flexibilities and dynamic coupling for positions in the linker region of the lactose repressor protein. This region contains nonconserved positions for which substitutions alter DNA-binding affinity. We first chose to study 11 substitutions at position 52. In computations, substitutions showed long-range effects on flexibilities of DNA-binding positions, and the degree of flexibility change correlated with experimentally measured changes in DNA binding. Substitutions also altered dynamic coupling to DNA-binding positions in a manner that captured other experimentally determined functional changes. Next, we broadened calculations to consider the dynamic coupling between 17 linker positions and the DNA-binding domain. Experimentally, these linker positions exhibited a wide range of substitution outcomes: Four conserved positions tolerated hardly any substitutions ("toggle"), ten nonconserved positions showed progressive changes from a range of substitutions ("rheostat"), and three nonconserved positions tolerated almost all substitutions ("neutral"). In computations with wild-type lactose repressor protein, the dynamic couplings between the DNA-binding domain and these linker positions showed varied degrees of asymmetry that correlated with the observed toggle/rheostat/neutral substitution outcomes. Thus, we propose that long-range and noncanonical substitutions outcomes at nonconserved positions arise from rewiring long-range communication among functionally important positions. Such calculations might enable predictions for substitution outcomes at a range of nonconserved positions., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Conventionally, most amino acid substitutions at "important" protein positions are expected to abolish function. However, in several soluble-globular proteins, we identified a class of nonconserved positions for which various substitutions produced progressive functional changes; we consider these evolutionary "rheostats". Here, we report a strong rheostat position in the integral membrane protein, Na + /taurocholate (TCA) cotransporting polypeptide, at the site of a pharmacologically relevant polymorphism (S267F). Functional studies were performed for all 20 substitutions (S267X) with three substrates (TCA, estrone-3-sulfate, and rosuvastatin). The S267X set showed strong rheostatic effects on overall transport, and individual substitutions showed varied effects on transport kinetics (K m and V max ) and substrate specificity. To assess protein stability, we measured surface expression and used the Rosetta software (https://www.rosettacommons.org) suite to model structure and stability changes of S267X. Although buried near the substrate-binding site, S267X substitutions were easily accommodated in the Na + /TCA cotransporting polypeptide structure model. Across the modest range of changes, calculated stabilities correlated with surface-expression differences, but neither parameter correlated with altered transport. Thus, substitutions at rheostat position 267 had wide-ranging effects on the phenotype of this integral membrane protein. We further propose that polymorphic positions in other proteins might be locations of rheostat positions., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2020 Wiley Periodicals LLC.)

To achieve the full potential of pharmacogenomics, one must accurately predict the functional out comes 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 arange of outcomes. Finally, some nonconserved positions appear neutral, capable of accommodating all amino acid types without modifying function. This manuscript reviews the currently-known properties of rheost at 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 out comes at rheostat positions poses a challenge, utilization of this new frame of reference will further advance the application of pharmacogenomics., Competing Interests: Conflict of Interest. The authors declare that they have no conflict of interest.

Every method used to quantify biomolecular interactions has its own strengths and limitations. To quantify protein-DNA binding affinities, nitrocellulose filter binding assays with 32 P-labeled DNA quantify K d values from 10 -12 to 10 -8 M but have several technical limitations. Here, we considered the suitability of biolayer interferometry (BLI), which monitors association and dissociation of a soluble macromolecule to an immobilized species; the ratio k off /k on determines K d . However, for lactose repressor protein (LacI) and an engineered repressor protein ("LLhF") binding immobilized DNA, complicated kinetic curves precluded this analysis. Thus, we determined whether the amplitude of the BLI signal at equilibrium related linearly to the fraction of protein bound to DNA. A key question was the effective concentration of immobilized DNA. Equilibrium titration experiments with DNA concentrations below K d (equilibrium binding regime) must be analyzed differently than those with DNA near or above K d (stoichiometric binding regime). For ForteBio streptavidin tips, the most frequent effective DNA concentration was ~2 × 10 -9 M. Although variation occurred among different lots of sensor tips, binding events with K d ≥ 10 -8 M should reliably be in the equilibrium binding regime. We also observed effects from multi-valent interactions: Tetrameric LacI bound two immobilized DNAs whereas dimeric LLhF did not. We next used BLI to quantify the amount of inducer sugars required to allosterically diminish protein-DNA binding and to assess the affinity of fructose-1-kinase for the DNA-LLhF complex. Overall, when experimental design corresponded with appropriate data interpretation, BLI was convenient and reliable for monitoring equilibrium titrations and thereby quantifying a variety of binding interactions., (© 2020 The Protein Society.)

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.

To bridge biological and biochemical disciplines, the relationship between in vitro protein biochemical function and in vivo activity must be established. Such studies can (a) help determine whether properties measured in simple, dilute solutions extrapolate to the complex in vivo conditions and (b) illuminate cryptic biological factors that are new avenues for study. We have explored the in vivo-in vitro relationship for chimeras built from LacI/GalR transcription regulators. In prior studies of individual chimeras, amino acid changes that altered in vitro DNA binding affinity exhibited correlated changes in in vivo transcription repression. However, discrepancies arose when the two datasets were compared to each other: Although their DNA binding domains were identical and their in vitro binding affinities spanned the same range, their in vivo repression ranges differed by >50-fold. Here, we determined that the presence of endogenous ligand for one chimera further exacerbated the offset, but that different abilities to simultaneously bind and "loop" two DNA operators resolves the discrepancy. Indeed, results suggest that the lac operon can be looped by even weakly interacting repressor dimers. For looping-competent repressors, we measured in vitro binding to the secondary operator. Surprisingly, this was largely insensitive to amino acid changes in the repressor protein, which reflects altered specificity; this supports the emerging view that the locations of specificity determining positions can be unique to each protein homolog. In aggregate, this work illustrates how a comparative approach can enrich understanding of the in vivo-in vitro relationship and suggest unexpected avenues for future study., (© 2019 The Protein Society.)

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., (© 2018 Wiley Periodicals, Inc.)

For decades, protein engineers have endeavored to reengineer existing proteins for novel applications. Overall, protein folds and gross functions can be readily transferred from one protein to another by transplanting large blocks of sequence (i.e., domain recombination). However, predictably fine-tuning function (e.g., by adjusting ligand affinity, specificity, catalysis, and/or allosteric regulation) remains a challenge. One approach has been to use the sequences of protein families to identify amino acid positions that change during the evolution of functional variation. The rationale is that these nonconserved positions could be mutated to predictably fine-tune function. Evolutionary approaches to protein design have had some success, but the engineered proteins seldom replicate the functional performances of natural proteins. This Biophysical Perspective reviews several complexities that have been revealed by evolutionary and experimental studies of protein function. These include 1) challenges in defining computational and biological thresholds that define important amino acids; 2) the co-occurrence of many different patterns of amino acid changes in evolutionary data; 3) difficulties in mapping the patterns of amino acid changes to discrete functional parameters; 4) the nonconventional mutational outcomes that occur for a particular group of functionally important, nonconserved positions; 5) epistasis (nonadditivity) among multiple mutations; and 6) the fact that a large fraction of a protein's amino acids contribute to its overall function. To overcome these challenges, new goals are identified for future studies., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The AlloRep database (www.AlloRep.org) (Sousa et al., 2016) [1] compiles extensive sequence, mutagenesis, and structural information for the LacI/GalR family of transcription regulators. Sequence alignments are presented for >3000 proteins in 45 paralog subfamilies and as a subsampled alignment of the whole family. Phenotypic and biochemical data on almost 6000 mutants have been compiled from an exhaustive search of the literature; citations for these data are included herein. These data include information about oligomerization state, stability, DNA binding and allosteric regulation. Protein structural data for 65 proteins are presented as easily-accessible, residue-contact networks. Finally, this article includes example queries to enable the use of the AlloRep database. See the related article, "AlloRep: a repository of sequence, structural and mutagenesis data for the LacI/GalR transcription regulators" (Sousa et al., 2016) [1].

Protein families evolve functional variation by accumulating point mutations at functionally important amino acid positions. Homologs in the LacI/GalR family of transcription regulators have evolved to bind diverse DNA sequences and allosteric regulatory molecules. In addition to playing key roles in bacterial metabolism, these proteins have been widely used as a model family for benchmarking structural and functional prediction algorithms. We have collected manually curated sequence alignments for >3000 sequences, in vivo phenotypic and biochemical data for >5750 LacI/GalR mutational variants, and noncovalent residue contact networks for 65 LacI/GalR homolog structures. Using this rich data resource, we compared the noncovalent residue contact networks of the LacI/GalR subfamilies to design and experimentally validate an allosteric mutant of a synthetic LacI/GalR repressor for use in biotechnology. The AlloRep database (freely available at www.AlloRep.org) is a key resource for future evolutionary studies of LacI/GalR homologs and for benchmarking computational predictions of functional change., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

As proteins evolve, amino acid positions key to protein structure or function are subject to mutational constraints. These positions can be detected by analyzing sequence families for amino acid conservation or for coevolution between pairs of positions. Coevolutionary scores are usually rank-ordered and thresholded to reveal the top pairwise scores, but they also can be treated as weighted networks. Here, we used network analyses to bypass a major complication of coevolution studies: For a given sequence alignment, alternative algorithms usually identify different, top pairwise scores. We reconciled results from five commonly-used, mathematically divergent algorithms (ELSC, McBASC, OMES, SCA, and ZNMI), using the LacI/GalR and 1,6-bisphosphate aldolase protein families as models. Calculations used unthresholded coevolution scores from which column-specific properties such as sequence entropy and random noise were subtracted; "central" positions were identified by calculating various network centrality scores. When compared among algorithms, network centrality methods, particularly eigenvector centrality, showed markedly better agreement than comparisons of the top pairwise scores. Positions with large centrality scores occurred at key structural locations and/or were functionally sensitive to mutations. Further, the top central positions often differed from those with top pairwise coevolution scores: instead of a few strong scores, central positions often had multiple, moderate scores. We conclude that eigenvector centrality calculations reveal a robust evolutionary pattern of constraints-detectable by divergent algorithms--that occur at key protein locations. Finally, we discuss the fact that multiple patterns coexist in evolutionary data that, together, give rise to emergent protein functions., (© 2015 Wiley Periodicals, Inc.)

To modulate transcription, a variety of input signals must be sensed by genetic regulatory proteins. In these proteins, flexibility and disorder are emerging as common themes. Prokaryotic regulators generally have short, flexible segments, whereas eukaryotic regulators have extended regions that lack predicted secondary structure (intrinsic disorder). Two examples illustrate the impact of flexibility and disorder on gene regulation: the prokaryotic LacI/GalR family, with detailed information from studies on LacI, and the eukaryotic family of Hox proteins, with specific insights from investigations of Ultrabithorax (Ubx). The widespread importance of structural disorder in gene regulatory proteins may derive from the need for flexibility in signal response and, particularly in eukaryotes, in protein partner selection., (© 2015 by The American Society for Biochemistry and Molecular Biology, Inc.)

In prokaryotes, the construction of synthetic, multi-input promoters is constrained by the number of transcription factors that can simultaneously regulate a single promoter. This fundamental engineering constraint is an obstacle to synthetic biologists because it limits the computational capacity of engineered gene circuits. Here, we demonstrate that complex multi-input transcriptional logic gating can be achieved through the use of ligand-inducible chimeric transcription factors assembled from the LacI/GalR family. These modular chimeras each contain a ligand-binding domain and a DNA-binding domain, both of which are chosen from a library of possibilities. When two or more chimeras have the same DNA-binding domain, they independently and simultaneously regulate any promoter containing the appropriate operator site. In this manner, simple transcriptional AND gating is possible through the combination of two chimeras, and multiple-input AND gating is possible with the simultaneous use of three or even four chimeras. Furthermore, we demonstrate that orthogonal DNA-binding domains and their cognate operators allow the coexpression of multiple, orthogonal AND gates. Altogether, this work provides synthetic biologists with novel, ligand-inducible logic gates and greatly expands the possibilities for engineering complex synthetic gene circuits.

Protein families might evolve paralogous functions on their common tertiary scaffold in two ways. First, the locations of functionally-important sites might be "hard-wired" into the structure, with novel functions evolved by altering the amino acid (e.g. Ala vs Ser) at these positions. Alternatively, the tertiary scaffold might be adaptable, accommodating a unique set of functionally important sites for each paralogous function. To discriminate between these possibilities, we compared the set of functionally important sites in the six largest paralogous subfamilies of the LacI/GalR transcription repressor family. LacI/GalR paralogs share a common tertiary structure, but have low sequence identity (≤ 30%), and regulate a variety of metabolic processes. Functionally important positions were identified by conservation and co-evolutionary sequence analyses. Results showed that conserved positions use a mixture of the "hard-wired" and "accommodating" scaffold frameworks, but that the co-evolution networks were highly dissimilar between any pair of subfamilies. Therefore, the tertiary structure can accommodate multiple networks of functionally important positions. This possibility should be included when designing and interpreting sequence analyses of other protein families. Software implementing conservation and co-evolution analyses is available at https://sourceforge.net/projects/coevolutils/.

The millions of protein sequences generated by genomics are expected to transform protein engineering and personalized medicine. To achieve these goals, tools for predicting outcomes of amino acid changes must be improved. Currently, advances are hampered by insufficient experimental data about nonconserved amino acid positions. Since the property "nonconserved" is identified using a sequence alignment, we designed experiments to recapitulate that context: Mutagenesis and functional characterization was carried out in 15 LacI/GalR homologs (rows) at 12 nonconserved positions (columns). Multiple substitutions were made at each position, to reveal how various amino acids of a nonconserved column were tolerated in each protein row. Results showed that amino acid preferences of nonconserved positions were highly context-dependent, had few correlations with physico-chemical similarities, and were not predictable from their occurrence in natural LacI/GalR sequences. Further, unlike the "toggle switch" behaviors of conserved positions, substitutions at nonconserved positions could be rank-ordered to show a "rheostatic", progressive effect on function that spanned several orders of magnitude. Comparisons to various sequence analyses suggested that conserved and strongly co-evolving positions act as functional toggles, whereas other important, nonconserved positions serve as rheostats for modifying protein function. Both the presence of rheostat positions and the sequence analysis strategy appear to be generalizable to other protein families and should be considered when engineering protein modifications or predicting the impact of protein polymorphisms.

LacI/GalR transcription regulators have extensive, non-conserved interfaces between their regulatory domains and the 18 amino acids that serve as 'linkers' to their DNA-binding domains. These non-conserved interfaces might contribute to functional differences between paralogs. Previously, two chimeras created by domain recombination displayed novel functional properties. Here, we present a synthetic protein family, which was created by joining the LacI DNA-binding domain/linker to seven additional regulatory domains. Despite 'mismatched' interfaces, chimeras maintained allosteric response to their cognate effectors. Therefore, allostery in many LacI/GalR proteins does not require interfaces with precisely matched interactions. Nevertheless, the chimeric interfaces were not silent to mutagenesis, and preliminary comparisons suggest that the chimeras provide an ideal context for systematically exploring functional contributions of non-conserved positions. DNA looping experiments revealed higher order (dimer-dimer) oligomerization in several chimeras, which might be possible for the natural paralogs. Finally, the biological significance of repression differences was determined by measuring bacterial growth rates on lactose minimal media. Unexpectedly, moderate and strong repressors showed an apparent induction phase, even though inducers were not provided; therefore, an unknown mechanism might contribute to regulation of the lac operon. Nevertheless, altered growth correlated with altered repression, which indicates that observed functional modifications are significant.

One emphasis of the Gibbs Conference on Biothermodynamics is the value of thermodynamic measurements for understanding behaviors of biological systems. In this study, the correlation between thermodynamic measurements of in vitro DNA binding affinity with in vivo transcription repression was investigated for two transcription repressors. In the first system, which comprised an engineered LacI/GalR homolog, mutational changes altered the equilibrium constant for binding DNA. Changes correlated with altered repression, but estimates of in vivo repressor concentration suggest a ≥25-fold discrepancy with in vitro conditions. In the second system, changes in ligand binding to BirA altered dimerization and subsequent DNA occupancy. Again, these changes correlate with altered in vivo repression, but comparison with in vitro measurements reveals a ~10-fold discrepancy. Further analysis of each system suggests that the observed discrepancies between in vitro and in vivo results reflect the contributions of additional equilibria to the transcription repression process., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Concomitant with the genomic era, many bioinformatics programs have been developed to identify functionally important positions from sequence alignments of protein families. To evaluate these analyses, many have used the LacI/GalR family and determined whether positions predicted to be "important" are validated by published experiments. However, we previously noted that predictions do not identify all of the experimentally important positions present in the linker regions of these homologs. In an attempt to reconcile these differences, we corrected and expanded the LacI/GalR sequence set commonly used in sequence/function analyses. Next, a variety of analyses were carried out (1) for the entire LacI/GalR sequence set and (2) for a subset of homologs with functionally-important "YxPxxxAxxL" motifs in their linkers. This strategy was devised to determine whether predictions could be improved by knowledge-based sequence sorting and-for some analyses-did increase the number of linker positions identified. However, two functionally important linker positions were not reliably identified by any analysis. Finally, we compared the new predictions to all known experimental data for E. coli LacI and three homologous linkers. From these, we estimate that >50% of positions are important to the functions of the LacI/GalR homologs. In corollary, neutral positions might occur less frequently and might be easier to detect in sequence analyses. Although analyses have successfully guided mutations that partially exchange protein functions, a better experimental understanding of the sequence/function relationships in protein families would be helpful for uncovering the remaining rules used by nature to evolve new protein functions., (Copyright © 2011 Wiley-Liss, Inc.)

The explosion of protein sequences deduced from genetic code has led to both a problem and a potential resource: Efficient data use requires interpreting the functional impact of sequence change without experimentally characterizing each protein variant. Several groups have hypothesized that interpretation could be aided by analyzing the sequences of naturally occurring homologues. To that end, myriad sequence/function analyses have been developed to predict which conserved, semi-conserved, and nonconserved positions are functionally important. These positions must be discriminated from the nonconserved positions that are functionally silent. However, the assumptions that underlie sequence analyses are based on experimental results that are sparse and usually designed to address different questions. Here, we use three homologues from a test family common to bioinformatics-the LacI/GalR transcription repressors-to test a common assumption: If a position is functionally important for one family member, it has similar importance in all homologues. We generated experimental sequence/function information for each nonconserved position in the 18 amino acids that link the DNA-binding and regulatory domains of three LacI/GalR homologues. We find that the functional importance of each position is preserved among the three linkers, albeit to different degrees. We also find that every linker position contributes to function, which has twofold implications. (1) Since the linker positions range from highly conserved to semi-conserved to nonconserved and contribute to affinity, selectivity, and allosteric response, we assert that sequence/function analyses must identify positions in the LacI/GalR linkers to be qualified as "successful". Many analyses overlook this region since most of the residues do not directly contact ligand. (2) No position in the LacI/GalR linker is functionally silent. This finding is inconsistent with another underlying principle of many analyses: Using sequence sets to discriminate important from non-contributing positions obligates silent positions, which denotes that most homologues tolerate a variety of amino acid substitutions at the position without functional change. Instead, additional combinatorial mutants in the LacI/GalR linkers show that particular substitutions can be silent in a context-dependent manner. Thus, specific permutations of sequence change (rather than change at silent positions) would facilitate neutral drift during evolution. Finally, the combinatorial mutants also reveal functional synergy between semi- and nonconserved positions. Such functional relationships would be missed by analyses that rely primarily upon co-evolution.

The lactose repressor protein (LacI) was among the very first genetic regulatory proteins discovered, and more than 1000 members of the bacterial LacI/GalR family are now identified. LacI has been the prototype for understanding how transcription is controlled using small metabolites to modulate protein association with specific DNA sites. This understanding has been greatly expanded by the study of other LacI/GalR homologues. A general picture emerges in which the conserved fold provides a scaffold for multiple types of interactions - including oligomerization, small molecule binding, and protein-protein binding - that in turn influence target DNA binding and thereby regulate mRNA production. Although many different functions have evolved from this basic scaffold, each homologue retains functional flexibility: For the same protein, different small molecules can have disparate impact on DNA binding and hence transcriptional outcome. In turn, binding to alternative DNA sequences may impact the degree of allosteric response. Thus, this family exhibits a symphony of variations by which transcriptional control is achieved.

In protein families, conserved residues often contribute to a common general function, such as DNA-binding. However, unique attributes for each homolog (e.g. recognition of alternative DNA sequences) must arise from variation in other functionally-important positions. The locations of these "specificity determinant" positions are obscured amongst the background of varied residues that do not make significant contributions to either structure or function. To isolate specificity determinants, a number of bioinformatics algorithms have been developed. When applied to the LacI/GalR family of transcription regulators, several specificity determinants are predicted in the 18 amino acids that link the DNA-binding and regulatory domains. However, results from alternative algorithms are only in partial agreement with each other. Here, we experimentally evaluate these predictions using an engineered repressor comprising the LacI DNA-binding domain, the LacI linker, and the GalR regulatory domain (LLhG). "Wild-type" LLhG has altered DNA specificity and weaker lacO(1) repression compared to LacI or a similar LacI:PurR chimera. Next, predictions of linker specificity determinants were tested, using amino acid substitution and in vivo repression assays to assess functional change. In LLhG, all predicted sites are specificity determinants, as well as three sites not predicted by any algorithm. Strategies are suggested for diminishing the number of false negative predictions. Finally, individual substitutions at LLhG specificity determinants exhibited a broad range of functional changes that are not predicted by bioinformatics algorithms. Results suggest that some variants have altered affinity for DNA, some have altered allosteric response, and some appear to have changed specificity for alternative DNA ligands.

Many mutations that impact protein function occur at residues that do not directly contact ligand. To understand the functional contributions from the sequence that links the DNA-binding and regulatory domains of the LacI/GalR homologues, we have created a chimeric protein (LLhP), which comprises the LacI DNA-binding domain, the LacI linker, and the PurR regulatory domain. Although DNA binding site residues are identical in LLhP and LacI, thermodynamic measurements of DNA binding affinity show that LLhP does not discriminate between alternative DNA ligands as well as LacI. In addition, small-angle scattering experiments show that LLhP is more compact than LacI. When DNA is released, LacI shows a 20 A increase in length that was previously attributed to unfolding of the linker. This change is not seen in apo-LLhP, even though the linker sequences of the two proteins are identical. Together, results indicate that long-range functional and structural changes are propagated across the interface that forms between the linker and regulatory domain. These changes could be mediated via the side chains of several linker residues that contact the regulatory domains of the naturally occurring proteins, LacI and PurR. Substitution of these residues in LLhP leads to a range of functional effects. Four variants exhibit altered affinity for DNA, with no changes in selectivity or allosteric response. Another two result in proteins that bind operator DNA with very low affinity and no allosteric response, similar to LacI binding nonspecific DNA sequences. Two more substitutions simultaneously diminish affinity, enhance allostery, and profoundly alter DNA ligand selectivity. Thus, positions within the linker can be varied to modulate different aspects of repressor function.