Using networks to identify fine structural differences between functionally distinct protein states

The vast increase in available data from the "-omics" revolution has enabled the fields of structural proteomics and structure prediction to make great progress in assigning realistic three-dimensional structures to each protein molecule. The challenge now lies in determining the fine structural details that endow unique functions to sequences that assume a common fold. Similar problems are encountered in understanding how distinct conformations contribute to different phases of a single protein's dynamic function. However, efforts are hampered by the complexity of these large, three-dimensional molecules. To overcome this limitation, structural data have been recast as two-dimensional networks. This analysis greatly reduces visual complexity but retains information about individual residues. Such diagrams are very useful for comparing multiple structures, including (1) homologous proteins, (2) time points throughout a dynamics simulation, and (3) functionally different conformations of a given protein. Enhanced structural examination results in new functional hypotheses to test experimentally. Here, network representations were key to discerning a difference between unliganded and inducer-bound lactose repressor protein (LacI), which were previously presumed to be identical structures. Further, the interface of unliganded LacI was surprisingly similar to that of the K84L variant and various structures generated by molecular dynamics simulations. Apo-LacI appears to be poised to adopt the conformation of either the DNA- or inducer-bound structures, and the K84L mutation appears to freeze the structure partway through the conformational transition. Additional examination of the effector binding pocket results in specific hypotheses about how inducer, anti-inducer, and neutral sugars exert their effects on repressor function.     

Rheostat functional outcomes occur when substitutions are introduced at nonconserved positions that diverge with speciation

When amino acids vary during evolution, the outcome can be functionally neutral or biologically‐important. We previously found that substituting a subset of nonconserved positions, “rheostat” positions, can have surprising effects on protein function. Since changes at rheostat positions can facilitate functional evolution or cause disease, more examples are needed to understand their unique biophysical characteristics. Here, we explored whether “phylogenetic” patterns of change in multiple sequence alignments (such as positions with subfamily specific conservation) predict the locations of functional rheostat positions. To that end, we experimentally tested eight phylogenetic positions in human liver pyruvate kinase (hLPYK), using 10–15 substitutions per position and biochemical assays that yielded five functional parameters. Five positions were strongly rheostatic and three were non‐neutral. To test the corollary that positions with low phylogenetic scores were not rheostat positions, we combined these phylogenetic positions with previously‐identified hLPYK rheostat, “toggle” (most substitution abolished function), and “neutral” (all substitutions were like wild‐type) positions. Despite representing 428 variants, this set of 33 positions was poorly statistically powered. Thus, we turned to the in vivo phenotypic dataset for E. coli lactose repressor protein (LacI), which comprised 12–13 substitutions at 329 positions and could be used to identify rheostat, toggle, and neutral positions. Combined hLPYK and LacI results show that positions with strong phylogenetic patterns of change are more likely to exhibit rheostat substitution outcomes than neutral or toggle outcomes. Furthermore, phylogenetic patterns were more successful at identifying rheostat positions than were co‐evolutionary or eigenvector centrality measures of evolutionary change.     

Subdividing repressor function: DNA binding affinity, selectivity, and allostery can be altered by amino acid substitution of nonconserved residues in a LacI/GalR homologue

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.     

Functionally important positions can comprise the majority of a protein's architecture

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.     

Integrated insights from simulation, experiment, and mutational analysis yield new details of LacI function

Protein structural change underlies many signal transduction processes. Although end-state structures are known for various allosteric proteins, intermediates are difficult to observe. Recently, targeted molecular dynamics simulation (TMD) was used to examine the conformational transition and predict relevant intermediates for wild-type lactose repressor (LacI). A catalog of involved residues suggests that the transition of this homodimer is asymmetric and that K84 is a prominent participant in the dynamic N-subdomain interface. Previous experiments indicated that hydrophobic substitutions at position 84 engender slowed, biphasic inducer binding kinetics, which might reflect the same phenomena observed in TMD. Here, we report biochemical confirmation that DNA and inducer binding remain allosterically linked in K84A and K84L, albeit with a differential smaller than that found in wild-type LacI. Other features of these mutant proteins are consistent with an allosteric conformational shift that approximates that of the wild type. As a consequence, these repressors can be utilized to explore an unanswered question about LacI function: How many inducers (one or two per dimer) are required to diminish operator affinity? The biphasic natures of the K84L and K84A inducer association rates allow direct correlation between the two distinct inducer binding events and operator release. Indeed, the kinetics of operator release for the K84A and K84L closely parallel those for the second inducer binding event. Together with implications from previous equilibrium results for wild-type and mutant proteins, these kinetic data demonstrate that binding of two inducers per dimeric DNA binding unit is required to release the operator in these variant LacI proteins.