DNA bending by bHLH charge variants

We wish to understand the role of electrostatics in DNA stiffness and bending. The DNA charge collapse model suggests that mutual electrostatic repulsions between neighboring phosphates significantly contribute to DNA stiffness. According to this model, placement of fixed charges near the negatively charged DNA surface should induce bending through asymmetric reduction or enhancement of these inter-phosphate repulsive forces. We have reported previously that charged variants of the elongated basic-leucine zipper (bZIP) domain of Gcn4p bend DNA in a manner consistent with this charge collapse model. To extend this result to a more globular protein, we present an investigation of the dimeric basic-helix–loop–helix (bHLH) domain of Pho4p. The 62 amino acid bHLH domain has been modified to position charged amino acid residues near one face of the DNA double helix. As observed for bZIP charge variants, DNA bending toward appended cations (away from the protein:DNA interface) is observed. However, unlike bZIP proteins, DNA is not bent away from bHLH anionic charges. This finding can be explained by the structure of the more globular bHLH domain which, in contrast to bZIP proteins, makes extensive DNA contacts along the binding face.


INTRODUCTION
The DNA double helix is a relatively stiff biopolymer with a persistence length of $150 bp. Over this length the average deflection of the helix axis is one radian under physiological salt conditions and in the absence of proteins (1)(2)(3). In contrast, very long DNA molecules spontaneously collapse to conformations whose end-to-end distances are proportional to the square root of the contour length (4). These characteristics imply that DNA is locally stiff yet globally flexible. Fundamental cellular processes such as nucleosome-mediated compaction (5), assembly and regulation of transcriptional machinery (6)(7)(8), and genetic recombination (9,10) require energetically unfavorable DNA deformation. Protein binding must provide the energetic driving force required for sharp DNA bending. We are interested in how asymmetric electrostatic interactions between proteins and DNA contribute to the generation of DNA-bending forces (3,(11)(12)(13)(14)(15)(16). Proteinmediated DNA bending is a crucial element in many essential biological processes, yet the contributions of electrostatic effects remain poorly understood.
Previous results from our laboratory support the idea that simple DNA-binding proteins can be engineered to induce DNA bending in opposite directions by positioning positive or negative charges asymmetrically on one face of the DNA duplex (13). Considerable efforts are still focused on understanding how neighboring cations and the phosphate diester backbone affect the stiffness of DNA (17)(18)(19)(20)(21)(22)(23). These findings support the proposals of Rich, Manning and others who theorized that charge-charge repulsions between neighboring phosphate diesters are at least partially responsible for the intrinsic rigidity of DNA (24)(25)(26).
The basic-leucine zipper (bZIP) domain of the Saccharomyces cerevisiae transcription factor Gcn4p ( Figure 1A) has been studied extensively as a biophysical model because of its simple secondary structure and orthology to the mammalian oncogenes Fos and Jun (27)(28)(29)(30)(31)(32)(33). The bZIP domain of Gcn4p is a 57 amino acid peptide sufficient for homo-dimerization and sequence-specific DNA binding (34). In the absence of DNA, the basic region of free bZIP peptides remains relatively disordered (35,36), whereas at high-protein concentration the leucine zipper region is natively folded into an alpha-helical coiledcoil. DNA binding potentiates folding, and crystal structures suggest that the basic region adopts a continuous alphahelical structure when bound to DNA (34,37). Previous efforts in our laboratory have shown that the bZIP domain is adaptable as a scaffold for charge variations, while preserving the ability of the modified bZIP peptides to bind DNA (13,15,16).
Electrophoretic assays have been widely used to detect DNA bending because mobility in native gels is related to DNA shape (38)(39)(40). Our own electrophoretic studies using DNA phasing probes were influenced by the work *To whom correspondence should be addressed. Tel: +1 507 284 9041; Fax: +1 507 284 2053; Email: maher@mayo.edu Ó 2006 The Author(s). This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/ by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
of Schepartz (30,41) and Kerppola (29,33,42) and have suggested a role for electrostatics in DNA bending. For example, the bZIP domain of transcription factor Gcn4p was used as a scaffold to position charged amino acids near one face of the DNA duplex to explore if local charges induce DNA bending ( Figure 1A). When bound to DNA, charged residues on Gcn4p have the potential to create laterally asymmetric charge distributions that appear to modulate inter-phosphate repulsions in a manner consistent with predictions (3,11,24,25,43,44). Specifically, DNA bends toward cationic charges that asymmetrically reduce phosphate repulsions, and away from anionic charges that enhance these repulsions.
DNA bending interpretations in electrophoretic phasing assays involving bZIP proteins can be controversial due to concerns that observed mobility changes could reflect changes in the shape, flexibility and/or gel-matrix interactions of the elongated peptide structure (28,33,45,46). In particular, DNA cyclization kinetics experiments have suggested that certain bZIP proteins do not bend DNA (46,47). Our previous work on the effect of appending globular domains onto bZIP peptides showed that apparent charge-dependent DNA bending was preserved in these more symmetrical complexes but was not detected by minicircle binding experiments (15,16). We reconciled the results by proposing that the bending potential for bZIP-DNA complexes is quite flat, so that different average geometries are observed with different methods (15). In the present work, we study a relatively globular protein, the basic helix-loop-helix (bHLH) domain of yeast Pho4p ( Figure 1B and C), as an alternative scaffold for positioning charges asymmetrically near DNA. Charge variants of Pho4p provide an independent opportunity for exploration of DNA bending by laterally asymmetric charge interactions and of protein-DNA complex flexibility.
Pho4p, a 312 amino acid transcription factor in S.cerevisae, belongs to the PHO gene family, which encodes >30 proteins responsible for regulation of exogenous inorganic phosphate uptake (48,49). At high exogenous phosphate, Pho4p is rendered inactive through phosphorylation by CDK complex orthologs, Pho80p-Pho85p (50). In phosphate-poor media Pho4p is hypophosphorylated and imported into the nucleus, where it acts together with Pho2p to activate genes required for phosphate scavenging (51) by specific recognition of the E-box sequence, 5 0 -CACGTG ( Figure 1D) (52,53). The C-terminal bHLH domain of Pho4p, the domain responsible for both selective DNA binding and homodimerization, includes a basic residue-rich region followed by a HLH fold (54)(55)(56)(57). The X-ray co-crystal structure of the 62 amino acid Pho4p bHLH bound to DNA has been solved at 2.8 s (58). The resulting model ( Figure 1B) shows that the loops of the bHLH dimers play a critical role in conferring stability and specificity to Pho4p DNA binding through extensive protein-phosphodiester backbone contacts.
Here we show that, as seen for the bZIP domain of Gcn4p, charge variants of the bHLH domain of Pho4p appear to induce DNA bending toward appended cations. Unlike bZIP peptides, anionic charge variants of the bHLH peptide do not cause DNA bending in the opposite direction. We propose that DNA bending toward the bound protein is prevented by DNA contacts made by the loop region of the bHLH dimer.

Oligonucleotides
Unmodified oligonucleotides were synthesized using standard phosphoramidite chemistry. Oligonucleotides were purified by denaturing PAGE (20% acrylamide; 29:1 acrylamide/ bisacrylamide, 8 M urea), extracted from gel fragments (10 mM Tris-HCl, 1 mM EDTA and 300 mM NaOAc, pH 5.2), and desalted using C 18 reverse phase cartridges (Sep-Pak; Waters Corporation). Oligonucleotide concentrations were determined at 260 nm using nearest neighbor molar extinction coefficients as described previously (11). Oligonucleotides were annealed at $50 mM in 250 mM NaCl, by heating to 95 C and cooling to room temperature overnight.
bHLH:DNA-binding affinity studies DNA probes for binding affinity assays containing either an E-box (5 0 -CGCGCGCACGTGCCG; E-box underlined) or CRE (5 0 -CGATGCATGACGTCATGC; CRE underlined) sequence were prepared by annealing a self-complementary oligonucleotide to generate duplexes with 3 0 recessed termini. The annealed probe (400 pmole) was end-labeled using 20 U of the Klenow fragment of E.coli DNA polymerase I and 10 mCi each of [a-32 P]dGTP and dCTP (3000 Ci/mmol) in 70 mM Tris-HCl, pH 7.6, 10 mM MgCl 2 and 5 mM DTT for 45 min at 37 C, followed by incubation with 5 mM each of dGTP and dCTP for 15 min to polymerize all termini. Radiolabeled probes were purified from unincorporated nucleotides using spin columns (Chromaspin TE-10; Clontech) and diluted for binding affinity assays. bHLH samples were dissolved in water and concentrations were determined by the method of Pace et al. (60) followed by dilution into binding buffer [1· phosphate-buffered saline containing 5% glycerol and 0.025% NP-40]. Binding reactions were performed by incubation of peptide dilutions (10 pM to 1 mM) with annealed radiolabeled DNA probes (100 pM) in a total volume of 20 ml of binding buffer supplemented with 100 mg/ml BSA. Binding reactions were incubated on ice for 30 min before analysis on native polyacrylamide gels (12% acrylamide; 29:1 acrylamide/bisacrylamide) in 0.5· TBE buffer run for 750 V-h at 22 C. Dried gels were analyzed using a STORM 840 scanner (Amersham). Equilibrium dissociation constant (K d ) values were estimated by measuring the fractional occupancy (q) of the radiolabeled duplex [D] as a function of the bHLH concentration [P]. The K d was estimated with non-linear least squares fitting to Equation 1, using KaleidagraphÔ software (Synergy Software).
1 Circular dichroism spectroscopy CD spectra were recorded using a Jasco 810 spectropolarimeter. Far UV-CD spectra were obtained in the continuous mode, collecting measurements every 1 nm with an averaging time of 5 s at 295 K between 200-260 nm. bHLH samples were 10 mM in 0.2 cm path length cells. DNA probes containing the E-box sequence (underlined) were generated by annealing the single self-complementary oligonucleotide (5 0 -CGCGCGCACGTGCCGCGCG) at 500 mM. Spectra were acquired in binding buffer with either 10 mM free bHLH peptides or bHLH complexes with E-box DNA at a molar ratio of 1:1.1. All baseline buffer and nucleic acid contributions to ellipticity were subtracted after spectral collection. Mean residual ellipticity (Q MRE , deg cm 2 /dmol) was calculated using the relation: where Q mdeg is the ellipticity in mdeg, l is the path length in cm, c is the molar peptide concentration and n is the number of residues in the peptide. Alpha-helical estimates of bHLH peptides were based upon the previously published equation (61): where Q MRE is the mean residual ellipticity at 222 nm. For Q 0 and Q 100 representing 0 and 100% helicity, values of À2340 and À30 300 (deg cm 2 /dmol) were used.
Electrophoretic phasing assays of DNA bending DNA probes for phasing analyses were generated and radiolabeled by PCR from derivatives of plasmids pDP-AP-1-21, pDP-AP-1-23, pDP-AP-1-25, pDP-AP-1-26, pDP-AP-1-28 and pDP-AP-1-30 (30) that had been modified through sitedirected mutagenesis of the original AP-1 binding site to generate E-box-binding sites suitable for bHLH binding. DNA probes were incubated with sufficient peptide for $50% mobility shift in 20 ml reactions containing binding buffer and 100 mg/ml BSA. Reactions were incubated on ice for 30 min before native PAGE (6% acrylamide; 29:1 acrylamide/bisacrylamide) in 0.5· TBE buffer for 2000 V-h at 22 C. Dried gels were analyzed as above.

Phasing analysis calculations
Phasing analysis was performed essentially as described previously (28,62,63). Migration distances for free and bound DNA probes were measured and normalized to the average mobility of each set of DNA probes according to: where m rel is the normalized mobility, m is the measured mobility of each probe and m avg is the average mobility of a group of phasing probes sharing a common bound peptide. Values of m rel were plotted against the spacer distance (bp spacing between the center of the E-box site and the center of the proximal A 5 -tract) for fitting to the phasing function: where m rel is the normalized mobility of the species, S is the spacer length (bp), P PH is the phasing period (10.5 bp/turn), S T is the trans spacer length (predicted length in bp where protein-induced and A 5 -tract bends are in trans), and A PH , the amplitude of the phasing function. The value of A PH was approximated by least squares curve fitting using Kalei-daGraph. The bend angle, a B , was then estimated from the empirical relation: where a C is the reference bend angle of the phased array of A 5 tracts (54 ), and k is an empirically determined constant ($1.09 in this case).

DNA bending energy predictions
The DDG contributions to the overall process of protein-DNA complex formation, defined as DG (variant) À DG (SAA), were estimated from experimental data or else were proposed as reasonable guesses that rationalize experimental observations. The individual DDG contributions that come from the experiment were estimated as follows: For protein folding, the folding reaction is considered to be two-state, with the equilibrium constant calculated as follows: where f native is the fraction of folded protein. The relative free energy for folding is calculated from the ratio of equilibrium constants: The free energy cost of DNA bending is calculated as described previously (64). The equation assumes a persistence length of 150 bp: where the bend angle Q, in degrees, is evenly distributed over L base pairs. For example, the free energy cost of 600 cal/mol shown in Figure 5 corresponds to a bend of 24.5 over 14 bp.
The total of all the DDG contributions associated with protein-DNA complex formation is given by the observed dissociation constants.

Experimental design
The native bHLH domain of Pho4p ( Figure 1B and C), representing residues 248-312 from the wild-type sequence (59), was modified to laterally position asymmetric charges near one face of the DNA. bHLH SHORT peptides were generated by insertion of three additional residues (XXX, Figure 1B and C) N-terminal to the basic domain to position charge modifications near the phosphodiester backbone of DNA. This strategy placed charged residues at positions similar to those studied previously using the bZIP domain of Gcn4p ( Figure 1A and B) (13). Comparison of the crystal structures shows that the bHLH peptide, despite its longer primary sequence, folds into a more globular structure ( Figure 1A and B). The bHLH LONG peptide group, while positioning charges in a manner identical to bHLH SHORT , was predicted from studies of bZIP peptides (unpublished data) to stabilize folding by addition of 12 N-terminal residues ( Figure 1C). The sequences of the CRE and E-box binding sites for Gcn4p and Pho4p, respectively, are shown in Figure 1D.

Binding affinities of bHLH charge variants
Before assessing their abilities to bend DNA, the DNAbinding affinities of Pho4p bHLH charge variants were determined by protein titrations against a constant concentration of an 18 bp DNA duplex containing the E-box-binding site ( Figure 1D), and were assayed by electrophoretic mobility shift ( Figure 2A). The fits shown in Figure 2B (bHLH SHORT ) and 2C (bHLH LONG ) were used to estimate apparent dissociation constants (K d ) for each of three charge variants for both peptide lengths ( Table 1). The data for both short and long bHLH peptides show that the neutral charge variant has the highest affinity for the E-box DNA, followed by the anionic form (Table 1). For both peptides, the cationic bHLH displayed 3-to 4-fold lower affinities for DNA ( Table 1). The weaker binding affinities of cationic (KKK) charge variants were initially puzzling, as positive charges near the DNA surface might be expected to promote protein-DNA interactions. It is noteworthy that the appendage of the 12 N-terminal residues in bHLH LONG peptides reduces the differences in DNA binding affinities among the charge variants ( Figure 2B and C and Table 1). Control experiments verified that all peptides retained proper DNA sequence specificity in binding assays (data not shown).

Folding of free and DNA-bound bHLH charge variants
To test if the observed binding affinity changes were due to changes in the coupling between protein folding and DNA binding, we characterized the folding of bHLH charge variants using circular dichroism (CD) spectroscopy (65). Changes in protein conformation induced upon DNA binding can be detected as spectral changes in the far UV CD region between 210 and 230 nm, a region largely unaffected by DNA, but useful in observing protein secondary structure. Under conditions comparable to those used in binding studies, far-UV CD data were collected for both free peptide and DNA-peptide complexes to determine the structural effects of charge insertion and DNA binding.
The far-UV CD spectra of the free bHLH peptides (10 mM) are shown in Figure 3. Based on the molar ellipticity values at   Table 2). The results for the neutral charge variant (SAA) are consistent with the expected helical content of the peptides (free bHLH SHORT : 43%; free bHLH LONG : 37%) based upon the number of bHLH residues expected to be helical when not bound to DNA. For free bHLH SHORT peptides, insertion of cationic residues (KKK) destabilizes folding whereas addition of anionic residues (EEE) stabilizes folding ( Figure 3A and Table 2). These effects are significantly reduced in bHLH LONG peptides, demonstrating the strong stabilizing effect the 12 amino acid leader segment. We have observed similar effects for corresponding charge variants of the bZIP domain of Gcn4p (data not shown). The ratio of molar ellipticity at 205 nm to the ellipticity at 222 nm, diagnostic for alpha helix content, confirms that unbound bHLH charge variants remain partially unfolded. Partial unfolding is consistent with a previous report that noted a stabilizing effect of buffer cation composition (66). Upon addition of 1.1 molar equivalents of cognate E-box DNA to the bHLH charge variant peptides, alpha-helical content was increased in all cases ( Figure 3 and Table 2). For the bHLH SHORT peptides, comparison of the results for the cationic variant (KKK) to the neutral (SAA) and anionic (EEE) variants suggests that the cationic variant remains incompletely alpha-helical, despite a clear conformational change upon DNA binding ( Figure 3A and Table 2). In contrast, bHLH LONG peptides displayed similar folding enhancements in the presence of DNA ( Figure 3B and Table 2). In all cases, the random-coil minimum at 207 nm shifted (by 5 nm) to 212 nm and was less prominent relative to the alpha-helical signature with a Q 212 :Q 222 ratio of $0.9 (66). This change is consistent with an increase in alpha-helical structure and a reduction in the random-coil folded state. The calculated helical content of the neutral bHLH:DNA complex is again consistent with predicted results in both long and short bHLH peptides. These findings confirm reports of far-UV CD spectral changes upon bHLH binding to DNA (66).
Our data highlight the importance of electrostatic effects in the folding of the bHLH domain. Placement of additional negative charges near the peptide N-terminus enhances the intrinsic helical content of the bHLH fold whereas additional positive charge reduces helical content. These results are consistent with the expected electrostatic interactions between charged residues near the terminus of an alpha helix and the intrinsic helix macrodipole (67,68), which are favorable for negatively charged residues near the N-terminus and unfavorable for positively charged residues. Our CD results show that cationic bHLH charge variant samples the unfolded state much more often than the other variants, explaining the $20fold decrease in binding affinity for the E-box DNA in the case of the bHLH SHORT cationic variant. The reduction in the charge-dependent variation of folding stabilities of the bHLH LONG peptide group further supports this notion, as the inserted charges are placed further from the amino terminus, and so the electrostatic-macrodipole interaction is greatly reduced leading to more uniform folding stabilities and presumably binding affinities. The reduction in binding affinity seen for the SAA variant and the lack of the expected increase in affinity upon improved folding of the KKK variant may be due to steric interference between the extended helix and the DNA or to loss of contacts formed previously by less structured DNA tails.

DNA bending by Pho4p charge variants
We employed a well-established electrophoretic DNA phasing assay that permits quantitative estimation of both the magnitude and direction of DNA bending. Our goal was to determine how asymmetric charge interactions at the Pho4p bHLH binding site (E-box) can lead to changes in observed DNA bending. Figure 4A illustrates the phasing probe design based on constructs originally developed by Schepartz (30) and Crothers (40). These $300 bp DNA duplexes contain a central binding site for Pho4p (E-box) variably spaced from an array of three phased A 5 -tracts that serves as a curved reference. This spacer distance is varied in the five phasing probes by a total of 10 bp (in 2-3 bp increments) in order to generate unique helical phasings between the phased A 5 tract and the E-box. DNA probes in which the phased A 5 tract and protein-induced bend occur on the same helical face (cis) will migrate more slowly due to a large composite helical deflection, whereas DNA probes in which the deformations lie on opposite helical faces (trans) will migrate more rapidly as the two deformations tend to cancel.
Results of electrophoretic phasing experiments for bHLH peptides are shown in Figure 4B and C. Probe migration rates were normalized, fit to a phasing function and plotted against spacer length ( Figure 4D). The calculated bend angles (Table 3) show that the free E-box site has little intrinsic bending (estimated bend angle of 1 -2 toward the major groove). bHLH binding bends the E-box site slightly, with Helical content, expressed as a fraction of the total number of residues, was determined as described in Materials and Methods (Equation 3). In the case of this bHLH SHORT peptide, the maximum helical content of the free peptide (Helix 1 + Helix 2 + loop helix) is predicted to be 28 of 65 residues ($43%), in the case of bHLH LONG peptide, the maximum helical content of the free peptide (Helix 1 + Helix 2 + loop helix) is predicted to be 28 of 75 residues ($37%). The maximum helical content of the bHLH:DNA peptide (Helix 1 + Helix 2 + loop helix + basic region) is predicted to be 43-47 of 65 residues ($66-72%) for bHLH SHORT and 48-52 of 75 residues ($64-69%).  Table 3). These findings are consistent with previous crystallographic studies that show very little DNA bending in the bHLH:DNA complex (58). In contrast to the neutral variants, the cationic charge variants (KKK) showed a strong apparent bend, estimated at 26 -29 , toward the minor groove at the center of the E-box. This amounts to $24 of net bending relative to the SAA:DNA complex. These data are consistent with previous reports for cationic bZIP charge variants based on electrophoretic phasing assays (13,15,63). Both the anionic bHLH charge variants (EEE) caused little DNA bending, yielding minor-groove directed bends of 3 -4 (Table 3). This result contrasts with previous observations for bZIP charge variants, where approximately equal but opposite DNA bending (anionic: major groove; cationic: minor groove) was observed for opposite charges. Why does the bHLH scaffold permit DNA bending in only one direction? Since the DNA-binding affinities and specificities of the neutral (SAA) and anionic (EEE) bHLH variants are comparable, the lack of DNA bending by bHLH (EEE) is presumably not due to misfolding or loss of DNA-binding activity. Inspection of the crystal structures of both bZIP and bHLH suggests that the inability of the anionic charge variant of bHLH to bend DNA toward the major groove is explained by the bulky loop residues of the bHLH, which are known to be in contact with DNA and involved in enhancing DNA binding and specificity (58). The Pho4p bHLH loop directly contacts DNA on the major groove face of the E-box. Unlike the case for bZIP peptides, this presents a considerable steric obstacle to DNA bending in that direction.
To explicate the linkages among charged bHLH peptide folding, DNA bending and changes in protein-DNA contacts, we sketch a plausible free energy balance model in Figure 5. Figure 5A compares DNA bending by bHLH charge variants (top) and bZIP charge variants (bottom), illustrating the steric hindrance that prevents DNA from binding toward the bHLH protein. Figure 5B rationalizes the bending, folding and binding affinity results coherently in terms of DDG contributions (relative free energies for protein-DNA complex formation) in the EEE and KKK complexes relative to the SAA complex. The total of all the DDG contributions relative to the SAA state must equal the observed DDG for binding, which is about +1.8 kcal/mol for bHLH SHORT and +1.0 kcal/mole for bHLH LONG . The description below uses measured or calculated DDG contributions where possible, as described in Materials and Methods. Otherwise, the values are chosen to be reasonable and to rationalize experimental observations. The point of the analysis is not to claim that these values are correct but simply to list some of the contributions that must be considered and to show that small free energy changes can rationalize DNA bending by the KKK variant, lack of bending by the EEE variant, and the experimentally determined DDG for bHLH:DNA complex formation.
The elements of the model are as follows: (i) The KKK variant disrupts folding, which decreases the available binding free energy relative to SAA because more of the binding free energy must be coupled to folding and/or because contacts made by the N-terminal tail are lost. This free energy cost is larger for the bHLH SHORT variant because it does not fold as well. The free energy cost for folding can be estimated from the CD data by assuming that the EEE peptide is 95% folded and calculating the equilibrium constants for folding. The free energy cost shown for disrupting Helix 1 contacts is not known from experiment. (ii) The free energy needed to bend DNA by $25 over $15 bp is small, only $0.6 kcal/mol [assuming a persistence length of 150 bp; (7)], and the bending is paid for by improved contacts with the N-terminal positive charges, again chosen arbitrarily to rationalize the existence of bending while being consistent with measured energetics. (iii) The improved Helix 1 contacts in the KKK variants cause some distortion of contacts on the loop side, which costs binding free energy. The net effect is that as long as the sum of DNA bending energy and the loss of loop contacts is less than the free energy gained from electrostatic contacts, then the bent state will be a free energy minimum. (iv) For the EEE variants, the small negative DDG expected from enhanced protein folding is opposed by unfavorable steric contacts made by the N-terminal helix and/or unfavorable EEE-DNA electrostatic interaction, leading to a slight decrease in binding affinity that is again smaller for the bHLH LONG variants. As discussed above, DNA bending is not observed because the free energy cost of compressing the protein-DNA interface is greater than the free energy gain from decreasing electrostatic repulsion by bending away from Helix 1.
Attempts to test some of these hypotheses directly using bHLH variants with either extensions of helix 1 or substituting residues with less bulky side-chains in the loop region were unsuccessful because these peptides failed to bind DNA (data not shown). These results are consistent with the model in Figure 5, although not conclusive. The essential result is that it is understandable how DNA bending can occur only with the KKK charge variant, even though the KKK peptide binds DNA less tightly than the SAA and EEE peptides: Contacts between the loop region and the DNA are disrupted by the apparent bending away from the KKK variants, and the disruption is reflected in the decreased binding affinity of the KKK variants. Electrostatic repulsion or structural changes destabilize the binding of EEE variants but cannot induce bending. Bend angle estimates are generated from Equations 2-4 as described in Materials and Methods. Mean and standard deviations were obtained from three independent measurements. b Bend direction from Equations 4-6 refers to the DNA groove that is narrowed by bending, relative to a reference frame at the center of the E-box sequence.

Summary
The work presented here supports the interpretations of previous experiments suggesting that electrostatic effects can cause DNA bending. Gel-based DNA bending assays have been compared previously to other DNA bending assays such as DNA cyclization and minicircle binding, with variable results (15,45,46). There has been concern that discrepancies might be explained by the elongated shape of bZIP peptides bound to DNA. The use of the more globular bHLH domain of Pho4p complements our previous elongated bZIP phasing assay results as it demonstrates that electrostatic effects in DNA bending are observed in the context of a less elongated, more globular peptide complex. We also observed electrostatic effects on DNA bending when we previously appended globular domains to the N-termini of bZIP peptides (15). Taken together, these results argue against the concern that electrophoretic phasing assay results are dominated by peptide shape. The difference in bending behavior between bHLH and bZIP charge variants likely results from a difference in how peptides pack against DNA, as the simple helical structure and lack of bulky domains on the bZIP peptide permits DNA bending in both directions, whereas the larger bHLH domain only permits bending away from the protein.
Our findings support the electrostatic collapse model of DNA bending and are consistent with (but do not prove) the somewhat controversial notion that DNA phosphate diester charges contribute significantly to DNA stiffness (3,11,25,26,44,69,70). Complementary approaches to these gel-based assays are currently under investigation. The red DNA resembles the DNA in the original structures, proposed to be the DNA shape for neutral appended residues. Steric interactions due to bulky loop residues of Pho4p allow appended cations (KKK charge variant) to cause DNA bending toward the minor groove (green), but do not allow substantial anion-mediated DNA bending (EEE charge variant) toward the major groove (blue). The bZIP charge variants have less steric hindrance and allow DNA bending in either direction. (B) Speculative energy balance for Pho4 EEE and KKK variants. DDG values relative to the SAA variant are shown. The values for protein folding, DNAbending energy and the total DDG are calculated as described in Materials and Methods. The other DDG values, for the disruption of H1 contacts upon substitution, the formation of H1 contacts upon DNA bending and the disruption of loop contacts upon DNA bending, are reasonable guesses that rationalize bending only for the KKK variant and sum to the correct total DDG .