|Home | About | Journals | Submit | Contact Us | Français|
Lactose repressor protein (LacI) utilizes an allosteric mechanism to regulate transcription in E. coli, and the transition between inducer- and operator-bound states has been simulated by targeted molecular dynamics (TMD). The side chains of amino acids 149 and 193 interact and were predicted by TMD simulation to play a critical role in the early stages of the LacI conformational change. D149 contacts IPTG directly, and variations at this site provide the opportunity to dissect its role in inducer binding and signal transduction. Single mutants at D149 or S193 exhibit minimal change in operator binding, and alterations in inducer binding parallel changes in operator release, indicating normal allosteric response. The observation that the double mutant D149A/S193A exhibits wild-type properties excludes the requirement for inter-residue hydrogen bond formation in the allosteric response. The double mutant D149C/S193C purified from cell extracts shows decreased sensitivity to inducer binding, while retaining wild-type binding affinities and kinetic constants for both operator and inducer. By manipulating cysteine oxidation, we show that the more reduced state of D149C/S193C responds to inducer more similarly to wild-type protein, whereas the more oxidized state displays diminished inducer sensitivity. These features of D149C/S193C indicate that the novel disulfide bond formed in this mutant impedes the allosteric transition, consistent with the role of this region predicted by TMD simulation. Together, these results establish the requirement for flexibility of spatial relationship between D149 and S193 rather than a specific D149-S193 interaction in the LacI allosteric response to inducer.
Allosteric regulation is a major mechanism of control in many biological processes, including cell signaling as well as genetic and metabolic regulation (1, 2). Although allosteric behavior is a common feature for many regulatory proteins, a detailed atomic-level description of this mechanism remains largely elusive. Lactose repressor protein (LacI) is a prototypic genetic regulatory protein that regulates transcription of the lac operon genes (3), which encode enzymes responsible for metabolism of lactose in Escherichia coli (3, 4). High affinity binding of repressor to operator sequences within the lac operon (5) precludes transcription of the lac metabolic enzymes by RNA polymerase (4). Diminished affinity for the target operator DNA and a consequent increase in lac mRNA transcription result from a conformational change in response to binding a metabolite of lactose (1,6-allolactose) to a site within LacI distant from the DNA site (4).
LacI is a well-characterized tetrameric protein comprised of identical monomers. The structure of the monomer can be dissected into three major elements (Figure 1A and 1B): an N-terminal DNA binding domain (amino acids 1–60); a core domain (amino acids 61–330) that binds inducer and forms the dimer interface; and a C-terminal domain (amino acids 331–360) that generates the tetramer interface (6–9). The core domain of LacI is structurally homologous with the venus fly-trap (VFT) family that includes bacterial periplasmic binding proteins (PBPs) involved in transport and the extracellular domain of G-protein coupled receptors (GPCRs) as well as other receptors that transduce extracellular stimuli into intracellular signals and exhibit subdomain movements associated with function (10–12). In this larger functional context, significant efforts have been made to study allosteric behavior in LacI, which brings enormous biochemical and structural information to facilitate more detailed analysis (4, 8, 13–17).
The three-dimensional structures of the core domain and full-length LacI in both operator- and IPTG- bound conformations are available (7, 8, 16). The shift in LacI structure between operator- and inducer-bound conformations requires both rotation and translation of the N-subdomain within the core domain (Figure 1B) (8, 16). Amino acids D149 and S193 are located near the IPTG binding pocket at the interface of the N- and C- subdomains (Figure 1A). Shifts in the side chains at D149 and S193 as well as substantial movement of the backbone of D149 are observed between DNA- and IPTG-bound LacI structures (Figure 1C). D149 occupies a key position at the start of a flexible loop that comprises residues from 149 to 156. In the DNA-bound structure, a single bond is formed between the backbone of S193 and the side chain of D149; however, in the IPTG-bound structure, the D149 side chain forms two hydrogen bonds to the S193 side chain in each monomer (Figure 1C) (8, 16). This alteration suggests D149 and S193 might be participants in determining the distinct conformational states of the repressor.
Recently, a TMD simulation was utilized to simulate the conformational transition pathways in LacI (18). This prediction for the operator-bound to IPTG-bound forms shows significant movement in the flexible loop containing D149 at early stages of the conformational change. Mobility in this flexible loop allows the side chain of residue 149 to move ~4Å, forming hydrogen bonds with the side chain at 193. These hydrogen bonds may secure the flexible loop in place (18). This event appears critical in the initial stages of propagating the message of inducer-binding through the N-subdomain of LacI to the DNA binding domain. In the intermediate state of TMD simulation (Figure 1D) (18), the side chain of residue 149 points to a position neither close to the IPTG-bound nor to the DNA-bound state. This observation further supports the hypothesis that the region around residue 149 contributes to transmitting the allosteric signal. Consistent with their predicted importance, alterations at positions 149 and 193 resulted in diminished response to IPTG binding in phenotypic screens (14). Further, because D149 is among the residues that directly contact inducer IPTG, exploring this region provides an opportunity to examine whether IPTG binding and signal transduction capacities can be uncoupled. To dissect the structural and functional influence of this specific region, multiple single substitutions at positions 149 or 193 were examined. Double substitutions, D149A/S193A and D149C/S193C, were constructed and their biochemical and biophysical characteristics examined to explore the role of specific side chain interaction and flexibility of the region in LacI allosteric response.
Plasmid pJC1 (9) containing the complete LacI gene was used as the expression vector in E. coli BLIM cells (19). All substitutions were derived by site-directed mutagenesis (Quickchange, Stratagene Inc.) using the corresponding oligonucleotides. Full sequencing of each LacI gene (SeqWright Inc.) verified the presence of only the expected mutation. For protein expression, cells transformed with specific plasmids were grown in 2 × YT liquid media in a shaking incubator at 37 for ~20 hours. The cells were centrifuged and resuspended in breaking buffer (0.2 M Tris-HCl, pH 7.6, 0.2 M KCl, 0.01 M Mg acetate, 5% glucose, 0.3 mM DTT and 0.3 µM PMSF) containing lysozyme and frozen at −20. The frozen cells were slowly thawed and treated with 80 µl of 10 mg/ml DNaseI and MgCl2 (final concentration of 10 mM). Following centrifugation, the supernatant was mixed with ammonium sulfate to final 37% saturation, centrifuged, and the resuspended pellet was dialyzed overnight against 0.09 M KP buffer (0.09 M potassium phosphate, pH 7.6, 0.3 mM DTT, and 5% glucose) at 4. After removing any precipitate by centrifugation, the supernatant was applied to a phosphocellulose column equilibrated in 0.09 M KP buffer and eluted with a gradient from 0.12–0.30 M KP buffer (0.12–0.30 M potassium phosphate, pH 7.6, 0.3 mM DTT, and 5% glucose). Protein purity was assessed by using SDS-PAGE and was ≥ 90%. Stoichiometric DNA binding activity for each repressor variant was determined by filter binding assay described in the following section. Activities for purified proteins were ≥ 90%.
DNA binding experiments were performed by nitrocellulose filter binding assay (20, 21) with 40-bp natural operator (5' – TGTTGTGTGGAATTGTGAGCGGATAACAATTTCACACAGG – 3') (Biosource International). This assay was carried out in FB buffer (0.01 M Tris-HCl, pH 7.4, 0.15 M KCl, 0.1 mM DTT, 0.1 mM EDTA, and 5% DMSO) with 100 µg/mL of bovine serum albumin. Protein at various dilutions (1 × 10−13 to 5 × 10−9 M) was first mixed with 32P-labeled DNA. To determine DNA activity under stoichiometric conditions, the percentage of active repressor was measured with DNA concentration at least 10-fold above the Kd. To measure DNA binding affinity, DNA concentration was set at least 10-fold below the Kd. After a short incubation, the mixtures were filtered through nitrocellulose using a 96-well plate and exposed to a Fuji phosphorimaging plate overnight. The plate was read by a Fuji phosphorimager, and the radioactivity was quantified and analyzed with Igor Pro (Wavemetrics) to determine the binding affinity with the following equation:
where Yobs is the observed level of retained radioactivity at a given protein concentration, Ymax is the level of radioactivity measured at saturation, c is the background radioactivity detected when no repressor was present, and Kd is the equilibrium dissociation constant. The value of the Hill coefficient, n, is generally ~1 for LacI.
IPTG binding to repressor proteins was monitored by the fluorescence emission shift for intrinsic tryptophan that has been shown to accompany IPTG binding (22). The protein solution was placed in a 1-cm2 quartz cuvette at room temperature, and the fluorescence emission intensity change at wavelength > 350 nm was monitored at varying IPTG concentration (22). Protein concentration was fixed at ~1.5 × 10−7 M monomer in fluorescence buffer containing 0.01 M Tris-HCl, 0.15 M KCl, pH 7.4, and IPTG concentration was in the range of 1 × 10−8 to 1 × 10−3 M. The experiment was performed using SLM-Aminco AB2 spectrofluorometer with an excitation wavelength of 285 nm. The emission spectra were recorded between 300 and 380 nm. The binding affinity for the protein was determined by fitting the data in Igor Pro using equation 2,
where Kd is the equilibrium dissociation constant; Yobs is the measured fluoresence at a given IPTG concentration; Ymax is the maximum change in fluorescence signal between zero and saturating inducer; c is a constant background coefficient; and n is the Hill coefficient, which is ~1 in absence of operator.
The operator release assay measures the effect of inducer sugar on DNA binding affinity (23). IPTG ranging from 10−8 to 10−3 M was added following mixture of the protein and labeled operator DNA in FB buffer with 100 µg/mL of bovine serum albumin. After incubation for ~30 min, nitrocellulose filter binding was used to monitor the level of DNA bound in the presence of varied inducer concentration. The concentration of DNA was at least 10-fold below the Kd obtained from the operator binding assay (~10−12 M), and the concentration of protein was set to ~80% of maximal DNA binding for the specific repressor protein. Data were analyzed by Igor Pro using a variation of equation 2. In this case, Kd is substituted by [IPTG]mid, corresponding to the concentration of IPTG required to release 50% of the protein-DNA complex. Yobs is the measured radioactivity at a specific inducer concentration; Ymax is the maximum radioactivity differential observed, and c is the radioactivity with no protein present.
The association rates for IPTG binding to LacI were measured using a stopped flow system equipped with fluorescence detection (Applied Photophysics Ltd.) (24). The association processes were monitored as the fluorescence decrease due to IPTG binding. The excitation wavelength was 285 nm, and fluorescence emission was monitored using a 350 nm cut-off filter. Time courses were measured for the association of repressor with IPTG under pseudo-first-order conditions. Final LacI concentration was 1.0 × 10−6 M, and IPTG concentrations were varied from 0.05 mM to 0.2 mM in buffer containing 0.01 M Tris-HCl, pH 7.4, 0.15 M KCl. The association rate constants were determined by fitting observed rate constant (kobs) using the equation:
where kassoc and kdissoc are the association and dissociation rate constants, respectively.
Operator dissociation rate analysis utilized the nitrocellulose binding assay (25). LacI (1.8 × 10−11 M) and 32P-labeled DNA (3.6 × 10−11 M) in a total volume of 10 mL were equilibrated for 15 min at room temperature in FB buffer containing 0.01 M Tris-HCl, pH 7.4, 0.15 M KCl, 0.1 mM EDTA, and 5% DMSO with 100 µg/mL of bovine serum albumin. An excess of unlabeled DNA over the initial DNA concentration (90-fold excess) was added to the equilibrated solution, immediately followed by inverting the tube 3 times gently. Duplicate aliquots (100 µl) were withdrawn and filtered on a nitrocellulose membrane at the desired times. Reactions were monitored until equilibrium was reached at ~ 2 hr. The equilibrium value, corresponding to background retention of DNA, was subtracted from each time point. The dissociation rate constant was determined by fitting the data to a single exponential model.
The oxidized state of the double-cysteine mutant was obtained by utilizing a Hitrap desalting column (GE Healthcare Corp.) to remove DTT from the protein solution. After elution from the column, the solution was dialyzed in air-purged fluorescence buffer (0.01 M Tris-HCl, 0.15 M KCl, pH 7.4) for ~1 hr to further enhance the oxidized state. The corresponding reduced state measurement was achieved by exposure of the protein to a high concentration of DTT and 6 M urea followed by rapid elution through a Hitrap desalting column using fluorescence buffer to remove excess DTT.
Ellman’s reagent, 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB), was used to determine the free sulfhydryl content in LacI (26, 27). A volume of 100 µl of the denatured protein sample (with a final monomer concentration in the range ~1–10 µM in the presence of 6 M urea) was added to 50 µl of reaction buffer (0.1 M sodium phosphate, pH 7.6, 1 mM EDTA) and 20 µl of Ellman’s reagent. To measure the background, 100 µl of reaction buffer was utilized instead of protein solution. After incubation at room temperature to complete the reaction (~15 min), the absorbance at 412 nm was determined against buffer background to quantify the amount of product released from the reaction between an active thiolate of protein cysteines and DTNB using an extinction coeffient ε = 14290 cm−1M−1 .
A series of single substitutions for D149 was generated to explore the functional impact of alterations in electrostatics, polarity, size, hydrogen-bonding capacity, and hydrophobic character. For operator binding affinity, D149 and S193 single mutants display values within ~2-fold of wild-type LacI in the absence of IPTG (Figure 2; Table 1). In the presence of 1 mM IPTG, all the single mutants have sensitivity to IPTG binding similar to wild-type LacI, with the possible exception of D149C (Figure 2; Table 1). For IPTG binding, substitutions with N, S and A at residue D149 decreased binding affinity only slightly. The affinities of D149Q and D149C for IPTG were decreased by ~4-fold, whereas those for D149T and D149E were decreased by ~13- and ~23-fold, respectively. Both S193A and S193C displayed comparable IPTG binding affinities to wild-type LacI (Figure 3; Table 1).
To examine allosteric communication between DNA- and IPTG-binding sites, we quantified the release of DNA from the LacI•DNA complex in response to inducer. This assay determines the concentration of IPTG required to release 50% of the operator, which reflects the combined contributions of IPTG binding and the concomitant conformational change. For position 149, N, S or A substitutions yielded essentially wild-type operator release values, whereas substitution with E, Q, T and C shifted the mid-point of operator release. For D149E, which has most severe effect, [IPTG]mid for operator release was increased ~50-fold. D149T had a ~20-fold higher value, and values for both D149Q and D149 C increased ~7-fold. For S193 single mutants, both S193A and S193C displayed operator release values comparable to wild-type LacI (Figure 3; Table 1). Therefore, a subset of the D149 substitutions results in diminished response to inducer binding with normal DNA affinity.
The alteration in IPTG release can derive from impaired IPTG binding or an impediment in the allosteric transition. The induction ratio, measured by the ratio of the operator release mid-point ([IPTG]mid) divided by the measured Kd for inducer binding, allows comparison between different LacI variants. Using this value, effects on allosteric communication can be distinguished from decreased IPTG binding affinity. For wild-type LacI, this ratio is ~1.5. The ratios for all mutants are within ~2-fold of wild-type LacI. Notice that although operator release [IPTG]mid values for D149E, D149Q, D149T and D149C are larger than that for wild-type LacI, their IPTG binding affinities are weaker than wild-type LacI, resulting in induction ratios comparable to wild-type (Table 1). These data suggest that substitutions at residue 149 do not significantly impact LacI allostery, even when inducer binding affinity is affected.
Since single mutants at D149 or S193 did not significantly disrupt LacI allostery, two double mutants, D149A/S193A and D149C/S193C, were generated. The TMD simulation suggested that hydrogen bond formation between the side chains of D149 and S193 might be important to function (Figure 1C, 1D). Therefore, D149A/S193A was constructed to abolish the potential for hydrogen bond formation. Further, because the flexible loop that begins with D149 undergoes high mobility in the simulation, D149C/S193C was designed to constrain movement with a disulfide bond. In both operator-and IPTG- bound LacI structures, the distances between the Cα atoms of residues 149 and 193 are within a range favorable for disulfide formation as found for other systems (28).
The D149A/S193A double mutant displays operator and inducer binding affinities, as well as an [IPTG]mid value for operator release, within ~2-fold of wild-type values (Figure 4; Table 2). The inducibility for D149A/S193A is close to wild-type, demonstrating that the ability to form hydrogen bonds between residue 149 and 193 is not critical for allosteric signal transduction. Interestingly, the double mutant D149C/S193C has operator binding affinity comparable to wild-type LacI in the absence of IPTG; whereas in the presence of 1 mM IPTG, this mutant still binds to operator, reflecting decreased sensitivity to inducer binding (Figure 4; Table 2). This double-cysteine mutant exhibited IPTG binding affinity within ~2-fold of wild-type LacI. Resistance of D149C/S193C to release of operator by inducer binding was confirmed by the operator release assay, with an over 12-fold increased [IPTG]mid value; further, the induction ratio is ~6-fold higher than that for wild-type LacI (Table 2).
To explore whether reciprocal behavior could be measured for operator impact on IPTG binding, [14C]-IPTG binding was measured in the assay buffer for operator binding. The retention of IPTG at protein concentrations just above the protein•DNA Kd and stoichiometric operator DNA was diminished by ~2-fold for D149C/S193C compared to wild-type LacI (data not shown), consistent with expectations from DNA binding results. These distinct properties of D149C/S193C meet the criteria for an allosteric LacI mutant, and flexibility for movement in the region encompassed by D149 and/or S193 appears critical for allosteric regulation.
To ensure that the differences for D149C/S193C were not kinetic in origin, wild-type and D149C/S193C association with IPTG and dissociation from operator were measured. The IPTG association rate constant for wild-type protein is 1.1 × 105 M−1s−1 at pH 7.4, consistent with data from previous work (24, 29, 30). The value for D149C/S193C is within ~2-fold of that for wild-type LacI (Figure 5; Table 3). Dissociation rate constants were derived from the plot of kassoc versus concentration of IPTG, and the D149C/S193C value was ~3-fold that for wild-type LacI. Thus, the disulfide bond between residues 149 and 193 does not substantially impact the rate constants for D149C/S193C binding to IPTG. The equilibrium constant for each reaction computed from the corresponding rate constants was comparable within ~2-fold to the value for the Kd determined directly from fluorescence titration (Table 3).
The operator dissociation rate constants were measured by adding unlabeled operator DNA to trap released repressor (25). The dissociation rate constants for wild-type LacI and D149C/S193C measured under same conditions were essentially identical (Figure 6; Table 3). These data demonstrate that introducing a disulfide bond between residues 149 and 193 does not significantly influence kinetics for repressor-operator dissociation.
To assess directly the formation of the disulfide bond in D149C/S193C, the number of free cysteines per monomer in this double mutant was determined by Ellman’s reaction in 6 M urea to provide a quantitative measure of available thiols (26, 27). The results demonstrate that the ratio of thionitrobenzoate anion produced over [monomer] is 2.9 ± 0.1 for wild-type LacI, indicating that each monomer exhibits the expected 3 free cysteines (residues 107, 140 and 281) (31). For the double mutant D149C/S193C, 4.0 ± 0.1 free cysteines are detected in the native monomeric state, suggesting there is a mixture of the oxidized and reduced forms of the protein present in solution. Following treatment with DTT and rapid removal by a desalting column, 4.6 ± 0.1 free cysteines are detected, indicating more reduced D149C/S193C is obtained. Also, the D149C/S193C mutant can be pushed toward the oxidized state by removing DTT through a desalting column and dialyzing in air-bubbled buffer. Under this condition, 3.3 ± 0.1 free cysteines per monomer are detected.
Since D149C/S193C with greater and lesser oxidation can be obtained experimentally, the impact of varied extent of disulfide bond formation was further examined. As controls, the thermodynamic properties of wild-type protein produced under the same conditions for both oxidization and reduction were also examined. The operator binding affinities for oxidized- and reduced-state D149C/S193C are decreased compared to untreated protein, consistent with results for wild-type LacI under the same conditions (Figure 7; Table 4). In the presence of a high concentration of 1 mM IPTG, D149C/S193C in both states can still bind to operator DNA with higher affinity than wild-type LacI (Figure 7; Table 4), with the oxidized protein binding ~9-fold more tightly. IPTG binding affinities for oxidized and reduced D149C/S193C are comparable to the untreated double mutant and wild-type LacI. Importantly, the operator release data demonstrate that the oxidized and reduced forms of D149C/S193C have distinct sensitivity to IPTG, reflected in the ~10-fold difference between their [IPTG]mid values (Table 4). Therefore, the induction ratios for oxidized and reduced D149C/S193C are > 10-fold different. The difference in response to IPTG binding between the oxidized and reduced forms of D149C/S193C confirms that the disulfide bond introduced impedes allosteric signal transduction within the LacI structure.
Despite broad exploration of LacI as a regulatory protein, neither D149 nor S193 have been targeted for biochemical analysis of substitutions. Previous phenotypic analysis demonstrates that alterations at these sites result in proteins insensitive to IPTG binding (14). However, since D149 is among the amino acids that have direct contact with inducer IPTG, this phenotype could be caused by either diminished IPTG binding affinity or impaired allosteric transition. In this study, a series of D149 and S193 point mutations were generated to examine the effects of different amino acid side chains on function. The induction ratios for the D149 single mutants are within ~2-fold that for wild-type protein, suggesting that coupling between inducer binding and the allosteric conformational change is not altered. In addition, both S193 single mutants, S193A and S193C, display wild-type properties. The data for single substitutions at either D149 or S193 demonstrate that substitutions at these sites either exert minimal effect on LacI binding properties or alter both inducer binding affinity and operator release, without significantly affecting the allosteric transition (i.e., there is no impact on the induction ratio).
In different ligand-bound structures, distinct interactions are observed between D149 or S193 and other surrounding amino acids. For example, the IPTG-bound LacI structure evinces at least one hydrogen bond between the R197 polypeptide backbone and the S193 backbone or side chain (8). In the operator-bound conformation, backbone hydrogen bonds are present between D149 and L128/Y126 (16). The double mutant D149A/S193A was designed to explore the importance of hydrogen bond formation by precluding the side-chain contribution at these positions. Interestingly, the double-alanine substitution does not affect binding properties of LacI, with operator binding and release properties similar to the wild-type protein. Loss of hydrogen bond capacity may be compensated by the proximal alanine-alanine hydrophobic interaction to promote allosteric signaling.
D149 is at the start of a loop (residues 149–156), for which the flexibility was postulated to be a key factor in the structural transition (18). The TMD results for LacI indicate that mobility in this loop accompanies the early steps in propagating the inducer-binding signal from the binding pocket to the DNA-binding site (18). To reduce flexibility within this region, residues that could form a disulfide bond were introduced by converting residues 149 and 193 to cysteines. This double-cysteine mutant indeed forms a disulfide bond between residues 149 and 193, and this variant binds with wild-type affinity to both operator and inducer. However, the D149C/S193C protein exhibits diminished impact of inducer binding on operator affinity compared to wild-type LacI and has an induction ratio ~9, over 6-fold greater than wild-type LacI and other mutants, including the single substitutions D149C and S193C.
To confirm that disulfide bond formation was the basis for the change in induction ratio, D149C/S193C was exposed to conditions designed to “push” towards either the oxidized- or reduced-state. Oxidized-D149C/S193C displays an even higher induction ratio (~16) than the protein as purified from bacterial cells (~9), whereas the reduced-state double mutant displays an induction ratio near that for wild-type LacI. This substantial difference between oxidized- and reduced-D149C/S193C is consistent with the hypothesis that flexibility in the 149–156 region is necessary for signal transmission at the initial stage of the allosteric transition in response to inducer binding.
Limited motions at either position 149 or 193 in the double-cysteine mutant D149C/S193C may influence various interactions with other residues. The TMD simulation indicates that D149 also forms a hydrogen bond with the backbone of F161, just before the flexible loop is stabilized by interaction with S193 (18). F161 is located in the center of a hydrophobic cluster that also experiences a high degree of mobility during the simulation. If interaction with F161 is impaired, the function of the entire hydrophobic cluster during the allosteric transition may be compromised. Further, the intermediate state of TMD simulation reveals the side chain of D149 transiently flips towards residue N125, which provides the opportunity to form a hydrogen bond between side chains of D149 and N125 during the process of transmitting the IPTG binding signal to the DNA-binding domain (18). However, an interaction between resides 149 and 125 is not present in the extreme DNA- and IPTG-bound conformations (8, 16).
For S193, a hydrogen bond between side chains of residues S193 and R197 is observed in the DNA-bound structure (16). During TMD simulation, this side-chain bond breaks, and new hydrogen bonds are generated between the backbone of R197 and side chain or backbone of S193 (18). However, the interaction between these side chains is not resumed in the inducer-bound structure, indicating its importance only in the DNA-bound conformation. Interactions between D149 or S193 and other amino acids may contribute to allosteric regulation in LacI, but the double-cysteine mutant D149C/S193C limits the movement of these residues, which may consequently affect the conformational change. Which of these interactions individually or in combination is crucial for the structural shifts that are responsible for allostery is not discernible from these data. However, the importance of flexibility in this region is clearly evident.
Protein flexibility appears to play a general role in propagating a signal from a localized perturbation (e.g., ligand binding or point mutation), allowing long-range communication that is fundamental to many protein functions (e.g., 32–36). In an analysis of more than 50 pairs of protein structures, each in two different states, structural motion carries information over long distances in both sequence and space (37). Interestingly, the motions observed favor residues located in the less constrained loops and at the protein surface, and ~2-fold greater motion is observed in allosteric proteins (37). A very recent model that describes allosteric motion for multiple proteins suggests that contact rearrangement is key to the quaternary changes that occur in LacI (38), consistent with the results reported in this study.
In LacI, substitution of other residues in this solvent-exposed flexible loop that encompasses amino acids 149–156 results in altered inducer and operator binding affinities (23). For example, despite induction ratios similar to wild-type LacI, the two mutations L148F and S151P are an “opposite” pair: L148F exhibits lowered affinity for operator and increased affinity for IPTG, whereas S151P has enhanced operator binding and diminished IPTG binding affinity (23). These results are consistent with a role for this flexible loop in regulating LacI allosteric equilibrium. Thus, residues in this flexible loop appear to make a key contribution to allosteric regulation, perhaps by impacting the range and nature of structural states available to the protein.
Overall, the mutants described herein provide specific insight into the role of D149 and S193 in LacI allosteric response. The experimental results demonstrate that impeding mobility in the D149/S193 region by disulfide bond formation interrupts allosteric signal transmission from the IPTG-binding pocket to the DNA-binding domain without significantly affecting IPTG and operator binding. Integration of biochemical examination of LacI variants with structural analysis and computational simulation provides an important avenue to complete our understanding of allosteric communication within LacI. This work provides a direct demonstration of the general utility of coupling experimental work with computational simulation for identifying regions of a protein crucial for function.
We thank members of the Matthews laboratory for helpful discussion and suggestions on experiments and Dr. Liskin Swint-Kruse for sage advice and critical reading of the manuscript.
†This work was supported by grants from the Robert A. Welch Foundation (C-576) and the National Institutes of Health (GM22441) to Kathleen S. Matthews.