Search

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. Protein secondary structures–α-helices and β-sheets–are generally assumed to be fixed over evolutionary history. By leveraging sequence information and sensitive statistical techniques, this work proposes that secondary structures in naturally occurring DNA-binding proteins switched in response to stepwise mutation. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 (nH) 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 nH values. We propose that substitution effects on both nH 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. PDB Code(s): 6XMH, 6XML, 6XMM and 6XMO; [ABSTRACT FROM AUTHOR]

Substitutions within the interdomain linkers of the AraC/XyIS family proteins RhaS and RhaR were tested to determine whether side chain identity or linker structure was required for function. Neither was found crucial, suggesting that the linkers do not play a direct role in activation, but rather simply connect the two domains.

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. [ABSTRACT FROM AUTHOR]

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 β-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; π-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. [ABSTRACT FROM AUTHOR]

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•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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]