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.)

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.)

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.)

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.)

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.)

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.)

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.)

The crystal structures of lactose repressor protein (LacI) provide static endpoint views of the allosteric transition between DNA- and IPTG-bound states. To obtain an atom-by-atom description of the pathway between these two conformations, motions were simulated with targeted molecular dynamics (TMD). Strikingly, this homodimer exhibited asymmetric dynamics. All asymmetries observed in this simulation are reproducible and can begin on either of the two monomers. Asymmetry in the simulation originates around D149 and was traced back to the pre-TMD equilibrations of both conformations. In particular, hydrogen bonds between D149 and S193 adopt a variety of configurations during repetitions of this process. Changes in this region propagate through the structure via noncovalent interactions of three interconnected pathways. The changes of pathway 1 occur first on one monomer. Alterations move from the inducer-binding pocket, through the N-subdomain beta-sheet, to a hydrophobic cluster at the top of this region and then to the same cluster on the second monomer. These motions result in changes at (1) side chains that form an interface with the DNA-binding domains and (2) K84 and K84', which participate in the monomer-monomer interface. Pathway 2 reflects consequent reorganization across this subunit interface, most notably formation of a H74-H74rsquo; pi-stacking intermediate. Pathway 3 extends from the rear of the inducer-binding pocket, across a hydrogen-bond network at the bottom of the pocket, and transverses the monomer-monomer interface via changes in H74 and H74rsquo;. In general, intermediates detected in this study are not apparent in the crystal structures. Observations from the simulations are in good agreement with biochemical data and provide a spatial and sequential framework for interpreting existing genetic data.

LacI and PurR are highly homologous proteins. Their functional units are homodimers, with an N-terminal DNA binding domain that comprises the helix-turn-helix (HTH), N-linker, and hinge regions from both monomers. Hinge structural changes are known to occur upon DNA dissociation but are difficult to monitor experimentally. The initial steps of hinge unfolding were therefore examined using molecular dynamics simulations, utilizing a truncated, chimeric protein comprising the LacI HTH/N-linker and PurR hinge. A terminal Gly-Cys-Gly was added to allow "dimerization" through disulfide bond formation. Simulations indicate that differences in LacI and PurR hinge primary sequence affect the quaternary structure of the hinge x hinge' interface. However, these alternate hinge orientations would be sterically restricted by the core domain. These results prompted detailed comparison of recently available DNA-bound structures for LacI and truncated LacI(1-62) with the PurR structure. Examination revealed that different N-linker and hinge contacts to the core domain of the partner monomer (which binds effector molecule) affect the juxtapositions of the HTH, N-linker, and hinge regions in the DNA binding domain. In addition, the two full-length repressors exhibit significant differences in the interactions between the core and the C-linker connection to the DNA binding domain. Both linkers and the hinge have been implicated in the allosteric response of these repressors. Intriguingly, one functional difference between these two proteins is that they exhibit opposite allosteric response to effector. Simulations and observed structural distinctions are correlated with mutational analysis and sequence information from the LacI/GalR family to formulate a mechanism for fine-tuning individual repressor function.

The repressor proteins of the LacI/GalR family exhibit significant similarity in their secondary and tertiary structures despite less than 35% identity in their primary sequences. Furthermore, the core domains of these oligomeric repressors, which mediate dimerization, are homologous with the monomeric periplasmic binding proteins, extending the issue of plasticity to quaternary structure. To elucidate the determinants of assembly, a structure-based alignment has been created for three repressors and four periplasmic binding proteins. Contact maps have also been constructed for the three repressor interfaces to distinguish any conserved interactions. These analyses show few strict requirements for assembly of the core N-subdomain interface. The interfaces of repressor core C-subdomains are well conserved at the structural level, and their primary sequences differ significantly from the monomeric periplasmic binding proteins at positions equivalent to LacI 281 and 282. However, previous biochemical and phenotypic analyses indicate that LacI tolerates many mutations at 281. Mutations at LacI 282 were shown to abrogate assembly, but for Y282D this could be compensated by a second-site mutation in the core N-subdomain at K84 to L or A. Using the link between LacI assembly and function, we have further identified 22 second-site mutations that compensate the Y282D dimerization defect in vivo. The sites of these mutations fall into several structural regions, each of which may influence assembly by a different mechanism. Thus, the 360-amino acid scaffold of LacI allows plasticity of its quaternary structure. The periplasmic binding proteins may require only minimal changes to facilitate oligomerization similar to the repressor proteins.