Search

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.

A significant number of eukaryotic regulatory proteins are predicted to have disordered regions. Many of these proteins bind DNA, which may serve as a template for protein folding. Similar behavior is seen in the prokaryotic LacI/GalR family of proteins that couple hinge-helix folding with DNA binding. These hinge regions form short alpha-helices when bound to DNA but appear to be disordered in other states. An intriguing question is whether and to what degree intrinsic helix propensity contributes to the function of these proteins. In addition to its interaction with operator DNA, the LacI hinge helix interacts with the hinge helix of the homodimer partner as well as to the surface of the inducer-binding domain. To explore the hierarchy of these interactions, we made a series of substitutions in the LacI hinge helix at position 52, the only site in the helix that does not interact with DNA and/or the inducer-binding domain. The substitutions at V52 have significant effects on operator binding affinity and specificity, and several substitutions also impair functional communication with the inducer-binding domain. Results suggest that helical propensity of amino acids in the hinge region alone does not dominate function; helix-helix packing interactions appear to also contribute. Further, the data demonstrate that variation in operator sequence can overcome side chain effects on hinge-helix folding and/or hinge-hinge interactions. Thus, this system provides a direct example whereby an extrinsic interaction (DNA binding) guides internal events that influence folding and functionality.

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.

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.

Allosteric modification of ligand binding is central to LacI transcription control. Recently, the conformational change between LacI operator- and inducer-bound states was simulated with targeted molecular dynamics (TMD) [Flynn, T. C., Swint-Kruse, L., Kong, Y., Booth, C., Matthews, K. S., and Ma, J. (2003) Protein Sci., 12, 2523-2541]. Atomic-level analyses of TMD results indicate the structural importance of the core pivot region that connects the N- and C-subdomains flanking the inducer-binding site. Further, a number of LacI mutations in the core pivot have been identified recently by their altered behaviors in phenotypic screens. Biochemical characterization of three of these variants-L148F, S151P, and P320A-provides an opportunity to directly explore the role of the core pivot in repressor function. For L148F, inducer IPTG binding affinity is strengthened, whereas O(1) operator DNA binding is diminished approximately 30-fold. In contrast, O(1) binding is increased for S151P, whereas IPTG binding is decreased. UV-difference spectroscopy and urea denaturation indicate long-range effects in both variants. Interestingly, P320A binds to DNA approximately 4-fold more tightly than wild-type, yet inducer binding is unaffected. To examine linkage between the core pivot and DNA binding domains, the L148F substitution was combined with Q60G, a previously known mutant with enhanced operator affinity. The double mutant exhibits the properties of both parent proteins, resulting in near wild-type DNA binding affinity and enhanced inducer sensitivity. These features may render Q60G/L148F more cost-effective in technological applications than wild-type repressor. As a group, the behaviors of the core pivot mutants are consistent with the allosteric structural role predicted for this region by TMD and reflect the significant long-range impact that single substitutions can elicit on protein function.

Two-dimensional nuclear magnetic resonance spectroscopy has been used to monitor proton-deuterium exchange rates (kobs) for more than 30 residues in turkey ovomucoid third domain. To test whether exchange is governed by global unfolding, rates were measured over a wide range of pH and temperatures where the change in the free energy of unfolding (delta Gzerou) is known [Swint, L., & Robertson, A. D. (1993) Protein Sci. 2, 2037-2049; Swint-Kruse, L., & Robertson, A. D. (1995) Biochemistry 34, 4724-4732]. Under conditions where EX2 kinetics are observed, a subset of 6-11 residues exhibits a one-to-one correlation with global stability. These residues are all located in central regions of secondary structures. Many other sites show varied degrees of correlation with delta Gzerou, while some are slower than expected on the basis of delta Gzerou alone. Preliminary evidence suggests that the latter is due to deviation from EX2 kinetics, even though experimental conditions are relatively mild (pH* 3 and 40 degrees C) compared to those in which deviations were observed for bovine pancreatic trypsin inhibitor. These results, together with similar observations for hen egg white lysozyme and barnase, suggest that EX2 kinetics should not be assumed when interpreting exchange studies.

Thermal denaturation of turkey ovomucoid third domain (OMTKY3) has been monitored with differential scanning calorimetry (DSC) and circular dichroism in H2O and D2O, pH 1.5 to 5 and ionic strength 0.01 to 0.71. Results from DSC experiments are in good agreement with spectroscopic studies [Swint, L., & Robertson, A.D. (1993) Protein Sci. 2, 2037-2049] and fit well to the two-state unfolding model. The average ratio of calorimetric to van't Hoff enthalpies is 0.99 +/- 0.03 (n = 16), and the average value of delta Cp is 620 +/- 20 cal/(mol K) (n = 7). The free energy of unfolding (delta G(u)o) increases in the presence of salt at both pH 1.5 and 4.5. This stabilization is not due to ion binding and probably results from screened repulsive interactions between the cationic groups of OMTKY3. At very low ionic strength, the change in delta G(u)o from pH 1.5 to 4.5, delta delta GpHo, is 3.5 +/- 0.2 kcal/mol. Few interactions between ionizing groups are affected by the addition of 200 mM KCl; delta delta GpHo decreases by only 0.4 +/- 0.3 kcal/mol. Comparison of delta delta GpHo with values calculated from the pKas of all six carboxyl groups in OMTKY3 [Schaller, W.S., & Robertson, A.D. (1995) Biochemistry 34, 4714-4723] suggests that some pKas in the denatured state may be lower than those of model compounds. Moreover, calculated values of delta delta GpHo are very sensitive to modest changes in the cooperativity of proton binding.(ABSTRACT TRUNCATED AT 250 WORDS)