Cooperative protein-DNA interactions: effects of KCl on lambda cI binding to OR

Specifically bound lambda repressor dimers promote adjacent non-specific binding

Genetic switches frequently include DNA loops secured by proteins. Recent studies of the lambda bacteriophage repressor (CI), showed that this arrangement in which the protein links two sets of three operators separated by approximately 2.3 kbp, optimizes both the stability and dynamics of DNA loops, compared to an arrangement with just two sets of two operators. Because adjacent dimers interact pairwise, we hypothesized that the odd number of operators in each set of the lambda regulatory system might have evolved to allow for semi-specific, pair-wise interactions that add stability to the loop while maintaining it dynamic. More generally, additional CI dimers may bind non-specifically to flanking DNA sequences making the genetic switch more sensitive to CI concentration. Here, we tested this hypothesis using spectroscopic and imaging approaches to study the binding of the lambda repressor (CI) dimer protein to DNA fragments. For fragments with only one operator and a short flanking sequence, fluorescence correlation spectroscopy measurements clearly indicated the presence of two distinct DNA-CI complexes; one is thought to have a non-specifically bound CI dimer on the flanking sequence. Scanning force micrographs of CI bound to DNA with all six operators revealed wild-type or mutant proteins bound at operator positions. The number of bound, wild-type proteins increased with CI concentration and was larger than expected for strictly specific binding to operators. In contrast, a mutant that fails to oligomerize beyond a dimer, D197G, only bound to operators. These data are evidence that CI cooperativity promotes oligomerization that extends from operator sites to influence the thermodynamics and kinetics of CI-mediated looping.

Effects of electrostatic interactions on ligand dissociation kinetics

We study unbinding of multivalent cationic ligands from oppositely charged polymeric binding sites sparsely grafted on a flat neutral substrate. Our molecular dynamics simulations are suggested by single-molecule studies of protein-DNA interactions. We consider univalent salt concentrations spanning roughly a 1000-fold range, together with various concentrations of excess ligands in solution. To reveal the ionic effects on unbinding kinetics of spontaneous and facilitated dissociation mechanisms, we treat electrostatic interactions both at a Debye-Hückel (DH) (or implicit ions, i.e., use of an electrostatic potential with a prescribed decay length) level and by the more precise approach of considering all ionic species explicitly in the simulations. We find that the DH approach systematically overestimates unbinding rates, relative to the calculations where all ion pairs are present explicitly in solution, although many aspects of the two types of calculation are qualitatively similar. For facilitated dissociation (FD) (acceleration of unbinding by free ligands in solution) explicit-ion simulations lead to unbinding at lower free-ligand concentrations. Our simulations predict a variety of FD regimes as a function of free-ligand and ion concentrations; a particularly interesting regime is at intermediate concentrations of ligands where nonelectrostatic binding strength controls FD. We conclude that explicit-ion electrostatic modeling is an essential component to quantitatively tackle problems in molecular ligand dissociation, including nucleic-acid-binding proteins.

Mechanisms that Determine the Differential Stability of Stx⁺ and Stx(-) Lysogens

Phages 933W, BAA2326, 434, and λ are evolutionarily-related temperate lambdoid phages that infect Escherichia coli. Although these are highly-similar phages, BAA2326 and 933W naturally encode Shiga toxin 2 (Stx⁺), but phage 434 and λ do not (Stx⁻). Previous reports suggest that the 933W Stx⁺ prophage forms less stable lysogens in E. coli than does the Stx⁻ prophages λ, P22, and 434. The higher spontaneous induction frequency of the Stx⁺ prophage may be correlated with both virulence and dispersion of the Stx2-encoding phage. Here, we examined the hypothesis that lysogen instability is a common feature of Stx⁺ prophages. We found in both the absence and presence of prophage inducers (DNA damaging agents, salts), the Stx⁺ prophages induce at higher frequencies than do Stx⁻ prophages. The observed instability of Stx⁺ prophages does not appear to be the result of any differences in phage development properties between Stx⁺ and Stx⁻ phages. Our results indicate that differential stability of Stx⁺ and Stx⁻ prophages results from both RecA-dependent and RecA-independent effects on the intracellular concentration of the respective cI repressors.

Single-cell kinetics of a repressilator when implemented in a single-copy plasmid

Synthetic genetic clocks, such as the Elowitz-Leibler repressilator, will be key regulatory components of future synthetic circuits. We constructed a single-copy repressilator (SCR) by implementing the original repressilator circuit on a single-copy F-plasmid. After verifying its functionality, we studied its behaviour as a function of temperature and compared it with that of the original low-copy-number repressilator (LCR). Namely, we compared the period of oscillations, functionality (the fraction of cells exhibiting oscillations) and robustness to internal fluctuations (the fraction of expected oscillations that would occur). We found that, under optimal temperature conditions, the dynamics of the two systems differs significantly, although qualitatively they respond similarly to temperature changes. Exception to this is in the functionality, in which the SCR is higher at lower temperatures but lower at higher temperatures. Next, by adding IPTG to the medium at low and high concentrations during microscopy sessions, we showed that the functionality of the SCR is more robust to external perturbations, which indicates that the oscillatory behaviour of the LCR can be disrupted by affecting only a few of the copies in a cell. We conclude that the SCR, the first functional, synthetic, single-copy, ring-type genetic clock, is more robust to lower temperatures and to external perturbations than the original LCR. The SCR will be of use in future synthetic circuits, since it complements the array of tasks that the LCR can perform.

Determinants of bacteriophage 933W repressor DNA binding specificity

We reported previously that 933W repressor apparently does not cooperatively bind to adjacent sites on DNA and that the relative affinities of 933W repressor for its operators differ significantly from that of any other lambdoid bacteriophage. These findings indicate that the operational details of the lysis-lysogeny switch of bacteriophage 933W are unique among lambdoid bacteriophages. Since the functioning of the lysis-lysogeny switch in 933W bacteriophage uniquely and solely depends on the order of preference of 933W repressor for its operators, we examined the details of how 933W repressor recognizes its DNA sites. To identify the specificity determinants, we first created a molecular model of the 933W repressor-DNA complex and tested the predicted protein-DNA interactions. These results of these studies provide a picture of how 933W repressor recognizes its DNA sites. We also show that, opposite of what is normally observed for lambdoid phages, 933W operator sequences have evolved in such a way that the presence of the most commonly found base sequences at particular operator positions serves to decrease, rather than increase, the affinity of the protein for the site. This finding cautions against assuming that a consensus sequence derived from sequence analysis defines the optimal, highest affinity DNA binding site for a protein.

The lysis-lysogeny decision of bacteriophage 933W: a 933W repressor-mediated long-distance loop has no role in regulating 933W P(RM) activity

Our data show that unlike bacteriophage λ, repressor bound at O(L) of bacteriophage 933W has no role in regulation of 933W repressor occupancy of 933W O(R)3 or the transcriptional activity of 933W P(RM). This finding suggests that a cooperative long-range loop between repressor tetramers bound at O(R) and O(L) does not form in bacteriophage 933W. Nonetheless, 933W forms lysogens, and 933W prophage display a threshold response to UV induction similar to related lambdoid phages. Hence, the long-range loop thought to be important for constructing a threshold response in lambdoid bacteriophages is dispensable. The lack of a loop requires bacteriophage 933W to use a novel strategy in regulating its lysis-lysogeny decisions. As part of this strategy, the difference between the repressor concentrations needed to bind O(R)2 and activate 933W P(RM) transcription or bind O(R)3 and repress transcription from P(RM) is <2-fold. Consequently, P(RM) is never fully activated, reaching only ∼25% of the maximum possible level of repressor-dependent activation before repressor-mediated repression occurs. The 933W repressor also apparently does not bind cooperatively to the individual sites in O(R) and O(L). This scenario explains how, in the absence of DNA looping, bacteriophage 933W displays a threshold effect in response to DNA damage and suggests how 933W lysogens behave as "hair triggers" with spontaneous induction occurring to a greater extent in this phage than in other lambdoid phages.

Quantitative DNase footprint titration: a tool for analyzing the energetics of protein-DNA interactions

A major goal in biomedical research is to determine the mechanisms responsible for gene regulation. However, the promoters and operators that control transcription are often complex in nature, containing multiple-binding sites with which DNA-binding proteins can interact cooperatively. Quantitative DNase footprint titration is one of the few techniques capable of resolving the microscopic binding affinities responsible for the macroscopic assembly process. Here, we present a step-by-step protocol for carrying out a footprint titration experiment. We then describe how to quantify the resultant images to generate individual-site binding curves. Finally, we derive basic equations for binding at each site and present an overview of the fitting process, applying it to the anticipated results. Users should anticipate that the footprinting experiment will take 3-5 d starting from DNA template isolation to image acquisition and quantitation.

Cooperative DNA binding by CI repressor is dispensable in a phage lambda variant

Complex gene regulatory circuits contain many interacting components. In principle, all of these components and interactions may be essential to the function of the circuit. Alternatively, some of them may be refinements to a simpler version of the circuit that improve its fitness. In this work, we have tested whether a particular property of a critical regulatory protein, CI, is essential to the behavior of the phage lambda regulatory circuit. In the lysogenic state, CI represses the expression of the lytic genes, allowing a stable lysogenic state, by binding cooperatively to six operators. A mutant phage lacking cooperativity because of a change in cI could not form stable lysogens; however, this defect could be suppressed by the addition of mutations that altered two cis-acting sites but did not restore cooperativity. The resulting triple mutant was able to grow lytically, form stable single lysogens, and switch to lytic growth upon prophage induction, showing a threshold response in switching similar to that of wild-type lambda. We conclude that cooperative DNA binding by CI is not essential for these properties of the lambda circuitry, provided that suppressors increase the level of CI. Unlike wild-type lysogens, mutant lysogens were somewhat unstable under certain growth conditions. We surmise that cooperativity is a refinement to a more basic circuit, and that it affords increased stability to the lysogenic state in response to environmental variations.

Indirect readout of DNA sequence by papillomavirus E2 proteins depends upon net cation uptake

The papillomavirus E2 proteins bind with high affinity to palindromic DNA sequences consisting of two highly conserved four base-pair sequences flanking a variable "spacer" of identical length (ACCG NNNN CGGT). While intimate contacts are observed between the bound proteins and conserved DNA in the available co-crystal structures, no contact is seen between the proteins and the spacer DNA. The ability of human papillomavirus strain 16 (HPV-16) E2 and bovine papillomavirus strain 1 (BPV-1) E2 to discriminate among binding sites with different spacer sequences is dependent on their sensitivity to the unique conformational and/or dynamic properties of the spacer DNA in a process termed "indirect readout". Differential sequence-specific K(+) uptake in low ionic strength solutions lacking Mg(2+) is observed upon E2 protein binding to sites containing the AATT, TTAA or ACGT spacer sequences. In contrast, the cation displacement typical of protein-DNA complex formation is observed at high K(+) concentrations or in the presence of Mg(2+). These results are interpreted to reflect the sequence-specific stabilization of bent DNA conformations by cations localized within the narrowed minor grooves of the protein-bound DNA and the intrinsic structure and flexibility of the DNA target. Mg(2+) differentially affects the binding of the HPV-16 E2 DNA binding domain (HPV16-E2/D) and the BPV-1 E2 DNA binding domain (BPV1-E2/D) to sites bearing different spacer sequences. This study suggests that monovalent and divalent cations contribute to the discrimination of DNA structure and flexibility that could in turn contribute to the specificity with which HPV16-E2/D and BPV1-E2/D mediate DNA replication and gene transcription.

Switching DNA-binding specificity by unnatural amino acid substitution

The specificity of protein–nucleic acid recognition is believed to originate largely from hydrogen bonding between protein polar atoms, primarily side-chain and polar atoms of nucleic acid bases. One way to design new nucleic acid binding proteins of novel specificity is by structure-guided alterations of the hydrogen bonding patterns of a nucleic acid–protein complex. We have used cI repressor of bacteriophage λ as a model system. In the λ-repressor–DNA complex, the ɛ-NH₂ group (hydrogen bond donor) of lysine-4 of λ-repressor forms hydrogen bonds with the amide carbonyl atom of asparagine-55 (acceptor) and the O6 (acceptor) of CG6 of operator site O_L1. Substitution of lysine-4 (two donors) by iso-steric S-(2-hydroxyethyl)-cysteine (one donor and one acceptor), by site-directed mutagenesis and chemical modification, leads to switch of binding specificity of λ-repressor from C:G to T:A at position 6 of O_L1. This suggests that unnatural amino acid substitutions could be a simple way of generating nucleic acid binding proteins of altered specificity.

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

A comparative three-dimensional model of the carboxy-terminal domain of the lambda repressor and its use to build intact repressor tetramer models bound to adjacent operator sites

A model for residues 93-236 of the lambda repressor (1gfx) was predicted, based on the UmuD(') crystal structure, as part of four intact repressor molecules bound to two adjacent operator sites. The structure of region 136-230 in 1gfx was found to be nearly identical to the independently determined crystal structure of the 132-236 fragment, 1f39, released later by the PDB. Later, two more tetrameric models of the lambda repressor tetramer bound to two adjacent operator sites were constructed by us; in one of these, 1j5g, the N-domain and C-domain coordinates and hence monomer-monomer and dimer-dimer interactions are almost the same as in 1gfx, but the structure of the linker region is partly based on the linker region of the LexA dimer in 1jhe; in the other, 1lwq, the crystalline tetramer for region 140-236 has been coopted from the crystal structure deposited in 1kca, the operator DNA and N-domain coordinates of which are same as those in 1gfx and 1j5g, but the linker region is partly based on the LexA dimer structures 1jhe and 1jhh. Monomer-monomer interactions at the same operator site are stabilized by exposed hydrophobic side chains in beta-strands while cooperative interactions are mostly confined to beta(6) and some adjacent residues in both 1gfx and 1j5g. Mutational data, existence of a twofold axis relating two C-domains within a dimer, and minimization of DNA distortion between adjacent operator sites allow us to roughly position the C-domain with respect to the N-domain for both 1gfx and 1j5g. The study correlates these models with functional, biochemical, biophysical, and immunological data on the repressor in the literature. The oligomerization mode observed in the crystal structure of 132-236 may not exist in the intact repressor bound to the operator since it is shown to contradict several published biochemical data on the intact repressor.

Coupled energetics of lambda cro repressor self-assembly and site-specific DNA operator binding II: cooperative interactions of cro dimers

The bacteriophage lambda relies on interactions of the cI and cro repressors which self assemble and bind the two operators (O(R) and O(L)) of the phage genome to control the lysogenic to lytic switch. While the self assembly and O(R) binding of cI have been investigated in detail, a more complete understanding of gene regulation by phage lambda also requires detailed knowledge of the role of cro repressor as it dimerizes and binds at O(R) sites. Since dimerization and operator binding are coupled processes, a full elucidation of the regulatory energetics in this system requires that the equilibrium constants for dimerization and cooperative binding be determined. The dimerization constant for cro has been measured as a prelude to these binding studies. Here, the energetics of cro binding to O(R) are evaluated using quantitative DNaseI footprint titration techniques. Binding data for wild-type and modified O(R) site combinations have been simultaneously analyzed in concert with the dimerization energetics to obtain both the intrinsic and cooperative DNA binding energies for cro with the three O(R) sites. Binding of cro dimers is strongest to O(R)3, then O(R)1 and lastly, O(R)2. Adjacently bound repressors exhibit positive cooperativity ranging from -0.6 to -1.0 kcal/mol. Implications of these, newly resolved, energetics are discussed in the framework of a dynamic model for gene regulation. This characterization of the DNA-binding properties of cro repressor establishes the foundation on which the system can be explored for other, more complex, regulatory elements such as cI-cro cooperativity.

Mutational analysis of the thermostable arginine repressor from Bacillus stearothermophilus: dissecting residues involved in DNA binding properties

Recently the crystal structure of the DNA-unbound form of the full-length hexameric Bacillus stearothermophilus arginine repressor (ArgR) has been resolved, providing a possible explanation for the mechanism of arginine-mediated repressor-operator DNA recognition. In this study we tested some of these functional predictions by performing site-directed mutagenesis of distinct amino acid residues located in two regions, the N-terminal DNA-binding domain and the C-terminal oligomerization domain of ArgR. A total of 15 mutants were probed for their capacity to repress the expression of the reporter argC - lacZ gene fusion in Escherichia coli cells. Substitutions of highly conserved amino acid residues in the alpha2 and alpha3 helices, located in the winged helix-turn-helix DNA-binding motif, reduced repression. Loss of DNA-binding capacity was confirmed in vitro for the Ser42Pro mutant which showed the most pronounced effect in vivo. In E. coli, the wild-type B. stearothermophilus ArgR molecule behaves as a super-repressor, since recombinant E. coli host cells bearing B. stearothermophilusargR on a multicopy vector did not grow in selective minimal medium devoid of arginine and grew, albeit weakly, when l -arginine was supplied. All mutants affected in the DNA-binding domain lost this super-repressor behaviour. Replacements of conserved leucine residues at positions 87 and/or 94 in the C-terminal domain by other hydrophobic amino acid residues proved neutral or caused either derepression or stronger super-repression. Substitution of Leu87 by phenylalanine was found to increase the DNA-binding affinity and the protein solubility in the context of a double Leu87Phe/Leu94Val mutant. Structural modifications occasioned by the various amino acid substitutions were confirmed by circular dichroism analysis and structure modelling.

The initiation of DNA base excision repair of dipyrimidine photoproducts

One of the major DNA repair pathways is base excision repair, in which DNA bases that have been damaged by endogenous or exogenous agents are removed by the action of a class of enzymes known as DNA glycosylases. One subset of the known DNA glycosylases has an associated abasic lyase activity that generates a phosphodiester bond scission. The base excision pathway is completed by the sequential action of abasic endonucleases, DNA polymerases, and DNA ligases. Base excision repair of ultraviolet (UV) light-induced dipyrimidine photoproducts has been described in a variety of prokaryotic and eukaryotic organisms and phages. These enzymes vary significantly in their exact substrate specificity and in the catalytic mechanism by which repair is initiated. The prototype enzyme within this class of UV-specific DNA glycosylases is T4 endonuclease V. Endonuclease V holds the distinction of being the first glycosylase (1) to have its structure solved by X-ray diffraction of the enzyme alone as well as in complex with pyrimidine dimer-containing DNA, (2) to have its key catalytic active site residues identified, and (3) to have its mechanism of target DNA site location determined and the biological relevance of this process established. Thus, the study of endonuclease V has been critical in gaining a better understanding of the mechanisms of all DNA glycosylases.

Carboxyl-terminal domain dimer interface mutant 434 repressors have altered dimerization and DNA binding specificities

Strong dimerization of the repressor, mediated by the carboxyl (C)-terminal domain, is a prerequisite for forming a specific complex with DNA and cooperative DNA binding to form tetramers. We have generated a computer model of the C-terminal domain of the 434 repressor based on the crystal structure of the homologous UmuD' protein. This model predicts that residues in the primary sequence between 93 and 168 contribute to the dimer interface. We changed several amino acid residues located in this region. Gel filtration and crosslinking assays were used to characterize the strength and specificity of dimerization of the purified repressor C-terminal domain dimer interface mutants. These results indicate that amino acid residues K121, H139, D161 and N163 contribute to the strength and/or specificity of dimerization. The relative affinity of the bacteriophage 434 repressor for 434 operators is determined, in part, by the repressor's ability to detect sequence-dependent structural alterations in the non-contacted region at the center of an operator site. We find that the relative ability of C-terminal domain dimer interface mutant repressors to dimerize does not necessarily predict their relative abilities to bind DNA, and that these proteins are deficient in detecting non-contacted base-dependent differences in operator strength. Our results show that the structure of the DNA in complex with these mutant proteins differs from that found in wild-type repressor-operator complexes, even though the sites of these mutations lie in a separate domain from that which contacts the DNA. These observations demonstrate that the structural integrity of the C-terminal domain dimer interface is required to appropriately orient the DNA binding information contained within the DNA-contacting N-terminal domain.

Cooperative non-specific DNA binding by octamerizing lambda cI repressors: a site-specific thermodynamic analysis

Relationships between dimerization and site-specific binding have been characterized previously for wild-type and mutant cI repressors at the right operator (OR) of bacteriophage lambda DNA. However, the roles of higher-order oligomers (tetramers and octamers) that are also formed from these cI molecules have remained elusive. In this study, a clear correlation has been established between repressor oligomerization and non-specific DNA-binding activity. A modification of the quantitative DNase I footprint titration technique has been used to evaluate the degree of saturation of non-specific, OR-flanking lambda DNA by cI repressor oligomers. With the exception of one mutant, only those repressors capable of octamerizing were found to exhibit non-specific DNA-binding activity. The non-specific interaction was accurately modeled using either a one-dimensional, univalent, site-specific Ising lattice approximation, or a more traditional, multivalent lattice approach. It was found that non-specific DNA-binding by repressor oligomers is highly cooperative and energetically independent from site-specific binding at OR. Furthermore, the coupling free energy resolved for non-specific binding was similar to that of site-specific binding for each repressor, suggesting that similar structural elements may mediate the cooperative component of both binding processes. It is proposed that the state of assembly of the repressor molecule modulates its relative affinity for specific and non-specific DNA sequences. These specificities are allosterically regulated by the transmission of assembly-state information from the C-terminal domain, which mediates self-association and cooperativity, to the N-terminal domain, which primarily mediates DNA-binding. While dimers have a high affinity for their cognate sites within OR, tetramers and octamers may preferentially recognize non-specific DNA sequences. The concepts and findings developed in this study may facilitate quantitative characterization of the relationships between specific, and non-specific binding in other systems that utilize multiple modes of DNA-binding cooperativity.

Solvent perturbation of the allosteric regulation of aspartate transcarbamylase

Escherichia coli aspartate transcarbamylase (ATCase) catalyzes the first committed step in pyrimidine biosynthesis, the condensation of aspartate and carbamyl phosphate. ATCase is positively allosterically regulated by ATP and negatively regulated by CTP. We have used mild solvent perturbation to gain global molecular information about the mechanism of heterotropic allostery. The [NaCl], temperature, and osmotic pressure dependence of the enzymatic activity of ATCase has been examined in the presence and absence of allosteric effectors. The results indicate that: 1) Regulation of aspartate binding by CTP appears to involve a unique set of electrostatic interactions not involved in enzyme function in the presence of ATP or in the absence of effectors. 2) Aspartate binding is enthalpically driven in the presence and absence of allosteric effectors. 3) The apparent enthalpy and entropy of aspartate binding (delta H, delta S), and activation energy of catalysis (Ea) are substantially altered in the presence of CTP but not ATP. 4) The change in hydration of ATCase upon substrate binding is the same in the presence and absence of allosteric effectors. 5) The linkage between heterotropic and homotropic allostery is different for ATP and CTP.

The spacing between binding sites controls the mode of cooperative DNA-protein interactions: implications for evolution of regulatory circuitry

The CI repressors of lambdoid phages bind cooperatively to adjacent binding sites. Although these binding sites are found at complex operators containing three binding sites, cooperative binding occurs only between dimers bound at two of the sites, a mode of binding termed pairwise cooperativity. At the level of regulation, pairwise cooperativity allows the proper operation of the genetic switch. At the mechanistic level, it has been proposed to result from steric distortion of the complex, such that a protein dimer bound to a central site cannot contact both flanking dimers because it "leans" towards one of the two sites. We have tested this model using the CI repressor of phage HK022. Previous work suggested that reducing the spacing between adjacent operators might allow cooperative interaction among three dimers, a mode of cooperativity we term extended. Using a set of synthetic templates, we have shown that reducing the spacing allows extended cooperativity among three binding sites. Accordingly, the mode of cooperativity changes qualitatively in response to changes in the inter-site spacing. That is, the change in spacing has major functional consequences for the properties of the complex. We suggest that such changes can play important roles in the evolution of gene regulatory circuitry and of other processes involving nucleoprotein complexes.

An aromatic stacking interaction between subunits helps mediate DNA sequence specificity: operator site discrimination by phage lambda cI repressor

Sequence specific DNA binding by regulatory proteins provides the basis for regulation of initiation of transcription. A great deal of progress has been made toward understanding sequence specific recognition by individual protein subunits. An additional level of control that needs to be understood is that due to coupling between the subunits of oligomeric regulatory proteins. An example is the bacteriophage lambda cI repressor, a dimeric protein that regulates the lysogenic to lytic genetic switch of the phage. Two levels of specificity are critical to this regulation. First, like all transcriptional regulators, dimers distinguish operator from nonspecific DNA. Direct readout of the DNA sequence by the recognition helix is considered the well understood mechanism for this. However, differential affinity for O(R)1, O(R)2 and O(R)3 is equally critical to the switch because it mediates opposing regulation of divergent promoters. Site specificity at this second level is less well understood. Conformational adaptation by both the repressor and the different operators appears to be important. To evaluate how subunit-subunit interactions are involved in this process, we investigated the effects on both dimer stability and operator binding of amino acid substitutions at the contacts between the symmetrically related helices-5 in the dimer interface. Substitutions for Tyr88 alter dimer stability and greatly perturb differential operator affinity, but generally do not affect operator versus non-operator specificity. The pattern of these effects suggests that the geometry of the face-to-face aromatic stacking interaction between symmetrically related Tyr88 in each subunit, a group in the dimer interface but far removed from the DNA binding interface, plays a critical role in operator discrimination. Conformational changes in the tertiary structure of the subunits appears to be involved. By contrast, the significant effect of I84S substitution is to greatly decrease affinity for all three operators. Presumably, the altered packing of the dimer interface causes a quarternary structural change that moves the two helix-turn-helix motifs out of register with successive DNA major grooves.

Quantitative nucleic acids footprinting: thermodynamic and kinetic approaches

Quantitative footprinting techniques allow a detailed analysis of the thermodynamic forces that characterize nucleic acid-ligand interactions and ligand-induced changes in nucleic acid structure by separately resolving the intrinsic and cooperative Gibbs free energy changes describing the reactions being investigated. A new implementation of the quantitative footprinting technique is the application of stopped-flow techniques to the study of kinetic reactions.

Co-operative binding of two Trp repressor dimers to alpha- or beta-centred trp operators

The alpha-centred trp operator binds one dimer of the Trp repressor, whereas the beta-centred trp operator binds two dimers of the Trp repressor (Carey et al., 1991; Haran et al., 1992). The Trp repressor with a Tyr-Gly-7 substitution binds almost as well as the wild-type Trp repressor to the alpha-centred trp operator, but it does not bind to the beta-centred trp operator. This confirms that Tyr-7 is involved in the interaction between Trp repressor dimers, as seen in the crystal structure (Lawson and Carey, 1993). Further experiments with alpha-centred trp operator variants showed that positions +/-1 of the alpha-centred trp operators play a crucial role in tetramerisation. The two innermost base pairs of the alpha-centred trp operator are not involved in contacts with the dimer of the Trp repressor binding to it. However, substitutions in these positions (T-A to G-T) effectively transform the alpha-centred trp operator into a beta-centred trp operator, and thus encourage the binding of two Trp repressor dimers to this operator. Finally, we demonstrate, with suitable heterodimers, that one subunit of each dimer suffices to bind to a beta-centred trp operator.

Buffer composition mediates a switch between cooperative and independent binding of an initiator protein to DNA

The regulation of many biological processes, including DNA replication, is frequently achieved by protein-protein interactions, as well as protein-DNA interactions. Multiple protein-binding sites are often involved. For example, the replication of plasmid R6K involves binding of the initiator protein pi to seven 22-bp direct repeats (DR) in the gamma origin of replication (gamma ori). A mutant protein pi S87N has been isolated, that in Tris.borate buffer (TB) binds cooperatively to seven DR, whereas wild-type (wt) pi binds independently [Filutowicz et al., Nucleic Acids Res. 22 (1994) 4211-4215]. Surprisingly, we found that wt pi can also bind cooperatively when Tris.acetate (TA), Tris.succinate or Tris.glutamate buffers are used instead of TB. The cooperative binding of the wt pi protein was also observed in the TB buffer at high concentrations of Na2EDTA. These results suggest that pi may be able to assume two functionally distinct conformations as a result of either mutation or buffer composition. Moreover, we found that the mode of pi binding is determined not by the composition of the buffer in which the reaction was assembled, but by the composition of the electrophoresis buffer. We discuss the general implications of these findings.

A fluorescence spectroscopic study of substrate-induced conformational changes in glutaminyl-tRNA synthetase

Glutaminyl-tRNA synthetase from Escherichia coli is a member of a subgroup of aminoacyl-tRNA synthetases that do not catalyze ATP-PPi exchange in the absence of the cognate tRNA. Such behavior suggests conformational changes upon substrate binding. Two different fluorescent probes, pyrenylmaleimide and acrylodan, were used to specifically label a nonessential sulfhydryl group of GlnRS. Conformational changes induced by substrates were studied using glutaminyl-tRNA synthetase labeled with these two environment-sensitive probes. ATP was shown to cause a significant conformational change that alters the mode of binding to tRNA(Gln) to GlnRS. The alteration of the salt sensitivity pattern of tRNA(Gln) binding to GlnRS by ATP supports this. Binding of tRNA(Gln) causes a conformational change that may be different in nature for the ATP/GlnRS complex and free GlnRS. Hydrodynamic parameters deduced from fluorescence polarization studies and the use of a noncovalent probe indicate that the ATP-induced conformational change may not be global in character.

Individual-site binding data and the energetics of protein-DNA interactions

Individual-site isotherms for the binding of bacteriophage lambda repressor to the left and right lambda operators have been determined [D. F. Senear, M. Brenowitz, M. A. Shea, and G. K. Ackers (1986) Biochemistry, Vol. 25, pp. 7344-7354.] using the DNAse protection technique [footprinting; D. J. Galas and A. Schmitz (1978) Nucleic Acids Research, Vol. 5, pp. 3157-3170]. These extensive data have been interpreted with a quantitative model that emphasized cooperative interactions between adjacently bound ligands [occupied <==> occupied interactions; G. K. Ackers, A. D. Johnson, and M. A. Shea (1982) Proceedings of the National Academy of Science, USA, Vol. 79, pp. 1129-1133]. Overlooked in this model are the effects of cooperative interactions between a site containing a bound ligand and its neighboring unoccupied site (occupied <==> unoccupied interactions). This paper reinterprets the existing data with a model that considers occupied <==> unoccupied as well as occupied <==> occupied interactions. The results yield parameters that differ substantially from those already reported. A discussion on the advisability of ignoring occupied <==> unoccupied interactions is included.

Effects of anomalous migration and DNA to protein ratios on resolution of equilibrium constants from gel mobility-shift assays

Numerical resolution of binding constants from data generated by the gel mobility-shift assay is dependent on the experimental resolution of the different ligation states of the DNA. Previously we showed that the populations of the intermediate ligation states at partial saturation with the protein ligand are extremely sensitive to cooperativity (Senear, D. F. and Brenowitz, M. J. Biol. Chem. 1991, 266, 13661-13671). This makes accurate gel mobility-shift data extremely useful to the demonstration of cooperativity. However, the accuracy with which the intermediate ligation state populations are resolved has been questioned. Thus, two additional and related questions are now considered. First, what information is available if the intermediate ligation state populations are not used in the analysis of binding constants. Second, is accurate information obtained from those states under conditions of high DNA concentration. These questions are addressed by using the interactions of the lambda cI repressor protein with the three site operator, OR, and the interaction of the E. coli GalR protein with the single site operator, OE. Both simulated and experimental data are analyzed. The results point to two conclusions. First, precise resolution of all macroscopic constants for binding of proteins to DNA is critically dependent on the intermediate ligation state populations; resolution is limited to at most two DNA sites if these states are not used in the analysis. Second, when the DNA and protein concentrations used in the titrations are comparable, the resolution of binding constants is extremely sensitive to experimental uncertainty in the macromolecule concentrations.

Site-specific enthalpic regulation of DNA transcription at bacteriophage lambda OR

Binding of cI repressor to DNA fragments containing the three specific binding sites of the right operator (OR) of bacteriophage lambda was studied in vitro over the temperature range 5-37 degrees C by quantitative footprint titration. The individual-site isotherms, obtained for binding repressor dimers to each site of wild-type OR and to appropriate mutant operator templates, were analyzed for the Gibbs energies of intrinsic binding and pairwise cooperative interactions. It is found that dimer affinity for each of the three sites varies inversely with temperature, i.e., the binding reactions are enthalpy driven, unlike many protein-DNA reactions. By contrast, the magnitude of the pairwise cooperativity terms describing interaction between adjacently site-bound repressor dimers is quite small. This result in combination with the recent finding that repressor monomer-dimer assembly is highly enthalpy driven (with delta H degrees = -16 kcal mol-1) [Koblan, K. S., & Ackers, G. K. (1991) Biochemistry 30, 7817-7821] indicates that the associative contacts between site-bound repressors that mediate cooperativity are unlikely to be the same as those responsible for dimerization. The intrinsic binding enthalpies for all three sites are negative (exothermic) and nearly temperature-invariant, indicating no heat capacity changes on the scale of those inferred in other protein-DNA systems. However, the three operator sites are affected differentially by temperature: the intrinsic binding free energies for sites OR1 and OR3 change in parallel over the entire range, delta H0OR1 = -23.3 +/- 4.0 kcal mol-1 and delta H0OR3 = -22.7 +/- 1.2 kcal mol-1.(ABSTRACT TRUNCATED AT 250 WORDS)