Thermodynamics of DNA branching.

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.     

HU binding to bent DNA: a fluorescence resonance energy transfer and anisotropy study

HU, an architectural DNA-binding protein, either stabilizes DNA in a bent conformation or induces a bend upon binding to give other proteins access to the DNA. In this study, HU binding affinity for a bent DNA sequence relative to a linear sequence was investigated using fluorescence anisotropy measurements. A static bend was achieved by the introduction of two phased A4T4 tracts in a 20 bp duplex. Binding affinity for 20 bp duplexes containing two phased A-tracts in either a 5'-3' or 3'-5' orientation was found to be almost 10-fold higher than HU binding to a random sequence 20 bp duplex (6.1 vs 0.68 microM(-1)). The fluorescence technique of resonance energy transfer was used to quantitatively determine the static bend of the DNA duplexes and the HU-induced bend. DNA molecules were 5'-end labeled with fluorescein as the donor or rhodamine as the acceptor. From the efficiency of energy transfer, the end-to-end distance of the DNA duplexes was calculated. The end-to-end distance relative to DNA contour length (R/R(C)) yields a bend angle for the A-tract duplex of 45 +/- 7 degrees in the absence of HU and 70 +/- 3 degrees in the presence of HU. The bend angle calculated for the T4A4 tract duplex was 62 +/- 4 degrees after binding two HU dimers. Fluorescence anisotropy measurements reveal that HU binds in a 1:1 stoichiometry to the A4T4 tract duplex but a 2:1 stoichiometry to the T4A4 tract and random sequence duplex. These findings suggest that HU binding and recognition of DNA may be governed by a structural mechanism.     

Mechanism of spontaneous DNA-DNA interaction of homologous linear duplexes

Previously, we demonstrated the interaction of homologous linear duplexes with formation of four-way DNA structures on the model of five PCR products. We propose that homologous duplex interaction is initiated by the nucleation of several dissociated base pairs of the complementary ends of two fragments with Holliday junction formation, in which cross point migration occurs via spooling of DNA strands from one duplex to the other one, finally resulting in complete resolution into new or previously existing duplexes. To confirm that DNA-DNA interaction involves formation of four-way DNA structures with strand exchange at the cross point, we have demonstrated the strand exchange process between identical duplexes using homologous fragments, harboring either biotin label or (32)P-label. Incubation of the mixture resulted in the addition of (32)P-label to biotin-labeled fragments, and the intensity of (32)P-labeling of biotinylated fragments was dependent upon the incubation duration. DNA-DNA interaction is not based on surface-dependent denaturing, as Triton X-100 does not decrease the formation of complexes between DNA duplexes. The equilibrium concentration of Holliday junctions depends on the sequences of the fragment ends and the incubation temperature. The free energy of Holliday junction formation by the fragments with GC and AT ends differed by 0.6 kcal/mol. Electron microscopic analysis demonstrated that the majority of Holliday junctions harbor the cross point within a 300 base pair region of the fragment ends. This insight into the mechanism of homologous duplex interaction extends our understanding of different DNA rearrangements. Understanding of DNA-DNA interaction is of practical use for better interpretation and optimization of PCR-based analyses.     

Free energy of DNA duplex formation on short oligonucleotide microarrays

DNA/DNA duplex formation is the basic mechanism that is used in genome tiling arrays and SNP arrays manufactured by Affymetrix. However, detailed knowledge of the physical process is still lacking. In this study, we show a free energy analysis of DNA/DNA duplex formation these arrays based on the positional-dependent nearest-neighbor (PDNN) model, which was developed previously for describing DNA/RNA duplex formation on expression microarrays. Our results showed that the two ends of a probe contribute less to the stability of the duplexes and that there is a microarray surface effect on binding affinities. We also showed that free energy cost of a single mismatch depends on the bases adjacent to the mismatch site and obtained a comprehensive table of the cost of a single mismatch under all possible combination of adjacent bases. The mismatch costs were found to be correlated with those determined in aqueous solution. We further demonstrate that the DNA copy number estimated from the SNP array correlates negatively with the target length; this is presumably caused by inefficient PCR amplification for long fragments. These results provide important insights into the molecular mechanisms of microarray technology and have implications for microarray design and the interpretation of observed data.     

Structural insights into the effect of isonucleosides on B-DNA duplexes using molecular-dynamics simulations

Some structural insights into the conformations of the isonucleosides containing duplexes have been provided. Unrestrained molecular-dynamics simulations on 18-mer duplexes with isonucleosides incorporated at the 3'-end or in the center of one strand have been carried out with explicit solvent under periodic boundary conditions using the AMBER force field and the particle mesh Ewald method. The Watson–Crick hydrogen-bonding patterns of the duplexes studied remained intact throughout the simulation. For the modified duplexes, the changes observed in the inter-base pair parameters and backbone torsional angles were primarily localized at the isonucleoside-inserted area. All five structures studied remained in the B-form family. The decreased stacking abilities indicated by the large changes in inter-base pair parameters and the large changes in backbones made the modified duplexes show a minor thermal destabilization in comparison with native DNA. The MM_PBSA method for estimating binding free energies on two complementary strands was used. The results showed that the binding free energies of isonucleoside-incorporated DNA duplexes were lower than the native DNA duplex, which is in good agreement with experimental observations.      

Stabilization factors affecting duplex formation of peptide nucleic acid with DNA.

Peptide nucleic acid (PNA) is an oligonucleotide analogue in which the sugar-phosphate backbone is replaced by an N-(2-aminoethyl)glycine unit to which the nucleobases are attached. We investigated the thermodynamic behavior of PNA/DNA hybrid duplexes with identical nearest neighbors but with different sequences and chain lengths (5, 6, 7, 8, 10, 12, and 16 mers) to reveal whether the nearest-neighbor model is valid for the PNA/DNA duplex stability. CD spectra of 6, 7, and 8 mer PNA/DNA duplexes showed similar signal, while 10, 12, and 16 mer duplexes did not. The average difference in Delta G degrees (37) for short PNA/DNA duplexes with identical nearest-neighbor pairs was only 3.5%, whereas that of longer duplexes (10, 12, and 16 mers) was 16.4%. Therefore, the nearest-neighbor model seems to be useful at least for the short PNA/DNA duplexes. Thermodynamics of PNA/DNA duplexes containing 1--3 bulge residues were also studied. While the stability of the 12 mer DNA/DNA duplex decreased as the number of bulge bases increases, the number of bulge bases in PNA/DNA unchanged the duplex stability. Thus, the influence of bulge insertion in the PNA/DNA duplexes is different from that of a DNA/DNA duplex. This might be due to the different base geometry in a helix which may potentially make hydrogen bonds in a base pair and stacking interaction unfavorable compared with DNA/DNA duplexes.     

Molecular dynamics structures of peptide nucleic acid x DNA hybrid in the wild-type and mutated alleles of Ki-ras proto-oncogene--stereochemical rationale for the low affinity of PNA in the presence of an AC mismatch

The low affinity of peptide nucleic acid (PNA) to hybridize with DNA in the presence of a mismatch endows PNA with a high degree of discriminatory capacity that has been exploited in therapeutics for the selective inhibition of the expression of point-mutated genes. To obtain a structural basis for this intriguing property, molecular dynamics simulations are carried out on PNA x DNA duplexes formed at the Ki-ras proto-oncogene, comprising the point-mutated (GAT), and the corresponding wild-type (GGT) codon 12. The designed PNA forms an A...C mismatch with the wild-type sequence and a perfect A...T pair with the point mutated sequence. Results show that large movements in the pyrimidine base of the A...C mismatch cause loss of stacking, especially with its penultimate base, concomitant with a variable mismatch hydrogen bond, including its occasional absence. These, in turn, bring about dynamic water interactions in the vicinity of the mismatch. Enthalpy loss and the disproportionate entropy gain associated with these are implicated as the factors contributing to the increase in free energy and diminished stability of PNA x DNA duplex with the A...C mismatch. Absence of these in the isosequential DNA duplex, notwithstanding the A...C mismatch, is attributed to the differences in topology of PNA x DNA vis-à-vis DNA duplexes. It is speculated that similar effects might be responsible for the reduced stability observed in PNA x DNA duplexes containing other base pair mismatches, and also in mismatch containing PNA x DNA duplexes.     

Insight into why pyrrolidinyl peptide nucleic acid binding to DNA is more stable than the DNA x DNA duplex

Molecular dynamics (MD) simulations and experimental measurements of the stability of a novel pyrrolidinyl PNA binding to DNA (PNA·DNA) in both parallel and antiparallel configurations were carried out. For comparison, simulations were also performed for the DNA·DNA duplex. The conformations of the three simulated systems were found to retain well-defined base pairing and base stacking as their starting B-like structure. A large gas-phase energy repulsion of the two negatively charged sugar-phosphate backbones of the DNA strands was found to reduce the stability of the DNA·DNA duplex significantly compared with that of the PNA·DNA complexes, especially in the antiparallel binding configuration. In addition, the antiparallel PNA·DNA was observed to be less solvated than that of the other two systems. The simulated binding free energies and the experimental melting temperatures for the three investigated systems are in good agreement, indicating that the antiparallel PNA·DNA is the most stable duplex.     

Structural and thermodynamic features of spiroiminodihydantoin damaged DNA duplexes

Oxidation of guanine or 8-oxo-7,8-dihydroguanine can produce spiroiminodihydantoin (Sp) R and S stereoisomers. Both in vitro and in vivo experiments have shown that the Sp stereoisomers are highly mutagenic, causing G --> C and G --> T transversion mutations. Therefore, they are of interest as potential endogenous cancer causing lesions. However, their structural properties in DNA duplexes remain to be elucidated. We have employed computational methods to study the Sp lesions in 11-mer DNA duplexes with A, C, G, and T partners. Molecular dynamics simulations have been carried out to obtain ensembles of structures, and the trajectories were employed to analyze the structures and compute free energies. The structural and thermodynamic analyses reveal that the Sp stereoisomers energetically favor positioning in the B-DNA major groove, with minor groove conformers also low energy in some cases, depending on the partner base. The R and S stereoisomers adopt opposite orientations with respect to the 5' to 3' direction of the modified strand. Both syn and anti glycosidic bond conformations are energetically feasible, with partner base and stereochemistry determining the preference. The lesions adversely impact base stacking and Watson-Crick hydrogen bonding interactions in the duplex, and cause groove widening. The chemical nature of the partner base determines specific hydrogen bonding and stacking properties of the damaged duplexes. The structural characteristics may relate to observed mutagenic properties of the Sp stereoisomers, including possible stereoisomer-dependent differences.     

Stability of DNA duplexes with Watson-Crick base pairs: a predicted model

The conformational stability (difference between the free energies of the folded and unfolded states, DeltaG degrees ) of a DNA duplex is considered as a function of component energy terms, hydrophobic, base stacking, hydrogen bonding, van der Waals, and electrostatic, and a trinucleotide-level helix stiffness parameter measured in terms of its Young's modulus. Hydrophobic and base stacking energy components were determined with the use of the crystal structure data of 30 DNA duplexes judicially selected within a resolution of 1.5 A, and hydrogen bonding, van der Waals and electrostatic terms were determined through an extensive review of experimental and theoretical studies. The stiffness indices for the trinucleotides were the ones realized by M. M. Gromiha [(2000) J. Biol. Phys. 26, 43-50] using the crystal structure data of 70 DNA duplexes. The unfolded state was treated in the classical way to determine its stability. Thermodynamically determined DeltaG degrees values for 111 DNA duplexes, with the number of base pairs ranging from 4 to 16, were selected in two sets, and the regression equation formed with one set was used to predict the stabilities of the other set, taking the energy components and the stiffness parameter to be independent variables. The computed energy terms indicate that the base stacking and hydrogen bonding forces are the dominant and the hydrophobic and electrostatic forces the weak partners in imparting stability to the duplexes. This model predicts DeltaG degrees values for DNA duplexes examined with a level of accuracy similar to that used for predictions made by the widely used nearest-neighbor models. The uniqueness of this model is that it combines the crystal and thermodynamic data for interpretation of conformational stability.     

Structural insight into DNA-assembled oligochromophores: crystallographic analysis of pyrene- and phenanthrene-modified DNA in complex with BpuJI endonuclease

The use of the DNA duplex as a supramolecular scaffold is an established approach for the assembly of chromophore aggregates. In the absence of detailed structural insight, the characterization of thus assembled oligochromophores is, today, largely based on solution-phase spectroscopy. Here, we describe the crystal structures of three DNA-organized chromophore aggregates. DNA hybrids containing non-nucleosidic pyrene and phenanthrene building blocks were co-crystallized with the recently described binding domain of the restriction enzyme BpuJI. Crystal structures of these complexes were determined at 2.7, 1.9 and 1.6 Å resolutions. The structures reveal aromatic stacking interactions between pyrene and/or phenanthrene units within the framework of the B-DNA duplex. In hybrids containing a single modification in each DNA strand near the end of the duplex, the two polyaromatic hydrocarbons are engaged in a face-to-face stacking orientation. Due to crystal packing and steric effects, the terminal GC base pair is disrupted in all three crystal structures, which results in a non-perfect stacking arrangement of the aromatic chromophores in two of the structures. In a hybrid containing a total of three pyrenes, crystal lattice induced end-to-end stacking of individual DNA duplexes leads to the formation of an extended aromatic π-stack containing four co-axially arranged pyrenes. The aromatic planes of the stacked pyrenes are oriented in a parallel way. The study demonstrates the value of co-crystallization of chemically modified DNA with the recombinant binding domain of the restriction enzyme BpuJI for obtaining detailed structural insight into DNA-assembled oligochromophores.     

The conformations of locked nucleic acids (LNA)

We have used 2D NMR spectroscopy to study the sugar conformations of oligonucleotides containing a conformationally restricted nucleotide (LNA) with a 2'-O, 4'-C-methylene bridge. We have investigated a modified 9-mer single stranded oligonucleotide as well as three 9- and 10-mer modified oligonucleotides hybridized to unmodified DNA. The single-stranded LNA contained three modifications whereas the duplexes contained one, three and four modifications, respectively. The LNA:DNA duplexes have normal Watson-Crick base-pairing with all the nucleotides in anti-conformation. By use of selective DQF-COSY spectra we determined the ratio between the N-type (C3'-endo) and S-type (C2'-endo) sugar conformations of the nucleotides. In contrast to the corresponding single-stranded DNA (ssDNA), we found that the sugar conformations of the single-stranded LNA oligonucleotide (ssLNA) cannot be described by a major S-type conformer of all the nucleotides. The nucleotides flanking an LNA nucleotide have sugar conformations with a significant population of the N-type conformer. Similarly, the sugar conformations of the nucleotides in the LNA:DNA duplexes flanking a modification were also shown to have significant contributions from the N-type conformation. In all cases, the sugar conformations of the nucleotides in the complementary DNA strand in the duplex remain in the S-type conformation. We found that the locked conformation of the LNA nucleotides both in ssLNA and in the duplexes organize the phosphate backbone in such a way as to introduce higher population of the N-type conformation. These conformational changes are associated with an improved stacking of the nucleobases. Based on the results reported herein, we propose that the exceptional stability of the LNA modified duplexes is caused by a quenching of concerted local backbone motions (preorganization) by the LNA nucleotides in ssLNA so as to decrease the entropy loss on duplex formation combined with a more efficient stacking of the nucleobases.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

DNA binding properties of the adenovirus DNA replication priming protein pTP

The precursor terminal protein pTP is the primer for the initiation of adenovirus (Ad) DNA replication and forms a heterodimer with Ad DNA polymerase (pol). Pol can couple dCTP to pTP directed by the fourth nucleotide of the viral genome template strand in the absence of other replication proteins, which suggests that pTP/pol binding destabilizes the origin or stabilizes an unwound state. We analyzed the contribution of pTP to pTP/pol origin binding using various DNA oligonucleotides. We show that two pTP molecules bind cooperatively to short DNA duplexes, while longer DNA fragments are bound by single pTP molecules as well. Cooperative binding to short duplexes is DNA sequence independent and most likely mediated by protein/protein contacts. Furthermore, we observed that pTP binds single-stranded (ss)DNA with a minimal length of approximately 35 nt and that random ssDNA competed 25-fold more efficiently than random duplex DNA for origin binding by pTP. Remarkably, short DNA fragments with two opposing single strands supported monomeric pTP binding. pTP did not stimulate, but inhibited strand displacement by the Ad DNA binding and unwinding protein DBP. These observations suggest a mechanism in which the ssDNA affinity of pTP stabilizes Ad pol on partially unwound origin DNA.     

B-S transition in short oligonucleotides

Stretching experiments with long double-stranded DNA molecules in physiological ambient revealed a force-induced transition at a force of 65 pN. During this transition between B-DNA and highly overstretched S-DNA the DNA lengthens by a factor of 1.7 of its B-form contour length. Here, we report the occurrence of this so-called B-S transition in short duplexes consisting of 30 basepairs. We employed atomic-force-microscope-based single molecule force spectroscopy to explore the unbinding mechanism of two short duplexes containing 30 or 20 basepairs by pulling at the opposite 5' termini. For a 30-basepair-long DNA duplex the B-S transition is expected to cause a length increase of 6.3 nm and should therefore be detectable. Indeed 30% of the measured force-extension curves exhibit a region of constant force (plateau) at 65 pN, which corresponds to the B-S transition. The observed plateaus show a length between 3 and 7 nm. This plateau length distribution indicates that the dissociation of a 30-basepair duplex mainly occurs during the B-S transition. In contrast, the measured force-extension curves for a 20-basepair DNA duplex exhibited rupture forces below 65 pN and did not show any evidence of a B-S transition.     

The thermodynamic advantage of DNA oligonucleotide 'stacking hybridization' reactions: energetics of a DNA nick.

'Stacking hybridization reactions' wherein two or more short DNA oligomers hybridize in a contiguous tandem orientation onto a longer complementary DNA single strand have been employed to enhance a variety of analytical oligonucleotide hybridization schemes. If the short oligomers anneal in perfect head-to-tail register the resulting duplex contains a nick at every boundary between hybridized oligomers. Alternatively, if the short oligomers do not hybridize precisely in register, i.e. single strand regions on the longer strand are left unbound, gaps are formed between regions where short oligomers bind. The resulting gapped DNA duplexes are considerably less stable than their nicked duplex analogs. Formation of base pair stacking interactions between neighboring oligomers at the nicks that do not occur in gapped duplexes has been proposed as the source of the observed added stability. However, quantitative evidence supporting this hypothesis for DNA has not been reported. Until now, a direct comparison of the thermodynamics of DNA nicks versus DNA gaps has not been performed. In this communication we report such a comparison. Analysis of optical melting experiments in a well defined molecular context enabled quantitative evaluations of the relative thermodynamic difference between nicked and gapped DNA duplexes. Results of the analysis reveal that a nick may be energetically favored over a gap by at least 1.4 kcal/mol and perhaps as much as 2.4 kcal/mol. The presence of a 5'phosphate at a nick or gap fails to significantly affect their stabilities.     

Cytarabine-induced destabilization of a model Okazaki fragment.

Cytarabine is a potent anticancer drug that interferes with elongation of the lagging strand at the replication fork during DNA synthesis. The effects of cytarabine substitution on the structural and thermodynamic properties of a model Okazaki fragment were investigated using UV hyperchromicity and 1H NMR spectroscopy to determine how cytarabine alters the physicochemical properties of Okazaki fragments that are intermediates during DNA replication. Two model Okazaki fragments were prepared corresponding to a primary initiation site for DNA replication in the SV40 viral genome. One model Okazaki fragment consisted of five ribo- and seven deoxyribonucleotides on the hybrid strand, together with its complementary (DNA) strand. The second model Okazaki fragment was identical to the first with the exception of cytarabine substitution for deoxycytidine at the third DNA nucleotide of the hybrid strand. Thermodynamic parameters for the duplex to single strand transition for each model Okazaki fragment were calculated from the concentration dependence of the T m at 260 nm. Cytarabine significantly decreased the stability of this model Okazaki fragment, decreasing the melting temperature from 46.8 to 42.4 degrees C at a concentration of 1.33 x 10(-5) M. The free energy for the duplex to single strand transition was 1.2 kcal/mol less favorable for the cytarabine-substituted Okazaki fragment relative to the control at 37 degrees C. Analysis of the temperature dependence of the imino1H resonances for the two duplexes demonstrated that cytarabine specifically destabilized the DNA:DNA duplex portion of the model Okazaki fragment. These results are consistent with inhibition of lagging strand DNA synthesis by cytarabine substitution resulting from destabilization of the DNA:DNA duplex portion of Okazaki fragments in vivo .     

HU binding to DNA: evidence for multiple complex formation and DNA bending

HU, a nonspecific histone-like DNA binding protein, participates in a number of genomic events as an accessory protein and forms multiple complexes with DNA. The HU-DNA binding interaction was characterized by fluorescence, generated with the guanosine analogue 3-methyl-8-(2-deoxy-beta-D-ribofuranosyl)isoxanthopterin (3-MI) directly incorporated into DNA duplexes. The stoichiometry and equilibrium binding constants of complexes formed between HU and 13 and 34 bp DNA duplexes were determined using fluorescence anisotropy and analytical ultracentrifugation. These measurements reveal that three HU molecules bind to the 34 bp duplexes, while two HU molecules bind to the 13 bp duplex. The data are well described by an independent binding site model, and the association constants for the first binding event for both duplexes are similar (approximately 1 x 10(6) M(-1)), indicating that HU binding affinity is independent of duplex length. Further analysis of the binding curves in terms of a nonspecific binding model is indicative that HU binding to DNA exhibits little to no cooperativity. The fluorescence intensity also increases upon HU binding, consistent with decreased base stacking and increased solvent exposure of the 3-MI fluorescence probe. These results are suggestive of a local bending or unwinding of the DNA. On the basis of these results we propose a model in which bending of DNA accompanies HU binding. Up to five complex bands are observed in gel mobility shift assays of HU binding to the 34 bp duplexes. We suggest that protein-induced bending of the DNA leads to the observation of complexes in the gel, which have the same molecular weight but different relative mobilities.