DNA-bending on ligand binding; change of DNA persistence length.

This paper simulates the helix-characteristic changes of apparent DNA persistence length caused by randomly distributed helix bends as induced, e.g., by DNA-bound ligand molecules. The parameters varied are the constant angle gamma of helix bending and the size alpha of the DNA drug binding site, but also the degree of DNA-ligand binding cooperativity and the helix-unwinding angle. If the size of the binding site is comparable with the helix pitch, the influence of phasing between helix bends and helix screw upon the apparent persistence length is obvious. In the accompanying paper experimental data are analyzed in terms of this theoretical background.     

Multimode Interaction of Hoechst 33258 with Eukaryotic DNA; Quantitative Analysis of the DNA Conformational Changes

Abstract

The interaction of the minor groove binding ligand Hoechst 33258 (Hoe) with natural DNA was investigated by high resolution titration rotational viscometry. Analysis of the concomitant DNA conformational changes was performed with two DNA samples of sufficiently different molar mass M, at 4°C, 22°C and 40°C, for Hoe/DNA-P ratios below r = 0.02. In this narrow r range several interaction modes could be resolved. The measured conformational changes were quantified in terms of relative changes of both apparent DNA persistence length, Δa/a, and hydrodynamically operative DNA contour length, ΔL/L. Δa/a(r) primarily is a measure of ligand-induced DNA helix stiffening, but both, Δa/a(r) and ΔL/L(r), generally depend also on ligand binding induced DNA bending or DNA unbending. The essential difference obviously is that Δa/a(r) is influenced by the randomly distributed helix bends and ΔL/L(r) by phased ones. The measurements performed at different temperatures deliver informations about existence and temperature dependent abolition of intrinsic helix curvature.

Both Hoe and netropsin (Nt) prefer binding to AT rich DNA segments, which are candidates for intrinsic DNA helix bends. But our data for Hoe interaction with calf thymus DNA (ctDNA) show characteristic differences to those for Nt-ctDNA interaction. Especially for Hoe, the mode of highest affinity is saturated already at a ligand concentration of roughly 1 nM (r = 0.0015 Hoe/DNA-P). It exhibits an unusually strong temperature dependence of the conformational DNA response. A Hoe-Nt competition experiment shows that Hoe binding to the sites of the very first Hoe mode is almost unaffected by bound Nt. But Hoe binding to the sites of the following Hoe modes does not occur due to the competition with Nt. Thus this mode of strongest Hoe-DNA interaction reflects a unique mechanism, possibly of high relevance for gene regulatory systems.     

Anisotropic flexibility of DNA depends on the base sequence. Conformation calculations of double-stranded tetranucleotides AAAA:TTTT, (AATT)2, (TTAA)2, GGGG:CCCC, (GGCC)2, (CCGG)2

The bending flexibility of six tetramers was studied in an assumption that they were extended in the both directions by regular double helices. The bends of B-DNA in different directions were considered. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less pronounced than in the perpendicular direction by the order of magnitude. Such an anisotropy is a feature of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degrees, are in agreement with the experimental value of DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove, moreover, they have an equilibrium bend of 6-12 degrees into this groove. The above inequality is caused by the stacking interaction of the bases. The bend in the central dimers is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is unessential, so that DNA remains within the limits of the B-family of forms: namely, when the helical axis is bent by 20 degrees the backbone dihedral angles vary by no more than 15 degrees. The obtained results are in accord with the X-ray structure of B-DNA dodecamer; they further substantiate our earlier model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b. p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in the equilibrium structure of certain DNA fragments. To the four "Calladine rules" two more can be added: the minor-groove steric clash of purines in the YR sequences are avoided by: (1) bending of the helix into the major groove; (2) increasing the distance between the base pairs (stretching the double helix).     

Structural perturbations induced in linear and circular DNA by the architectural protein HU from Bacillus stearothermophilus

HU is a small DNA-binding protein of eubacteria that is believed to induce or stabilize bending of the double helix and mediate nucleoid compaction in vivo. Although HU does not bind preferentially to specific DNA sequences, it is known to have high affinity for DNA sites containing structural anomalies, such as unpaired or mismatched bases, nicks, and four-way junctions. We have employed Raman spectroscopy to further investigate the structural basis of HU-DNA recognition in solution. Experiments were carried out on the homodimeric HU protein of Bacillus stearothermophilus (HUBst) and a 222-bp DNA fragment, which was isolated in linear (DNA(L222)) and circular (DNA(C222)) forms. In the absence of bound HUBst the Raman signatures of DNA(L222) and DNA(C222) are nearly superimposable, indicating that circularization produces no substantial change in the local B-DNA conformation. Conversely, the Raman signatures of DNA(L222) and DNA(C222) are perturbed significantly and specifically by HUBst binding. The HUBst-induced perturbations are markedly greater for the circularized DNA target. These results support an opportunistic molecular mechanism, in which HU binding is facilitated by intrinsic nonlinearity or flexibility in the DNA target. We propose that DNA segments which are bent or predisposed toward bending provide the high-affinity sites for HU attachment and nucleoid condensation. This model is consistent with the wide range of DNA bending angles reported in crystal structures of HU-DNA complexes.     

Abolition of Intrinsically Bent DNA Structure Components in AT Clusters by Netropsin Interaction; Titration Viscometric Investigations

Abstract

It is argued that the enhancement of the apparent DNA contour length by the specifically binding non-intercalating drug netropsin (Nt) (Reinert et al., NAR 9,2335, 1981) at very low Nt/DNA-phosphate ratios essentially is the result of an abolition of periodically arranged intrinsic helix bends in A · T rich tracts of base pairs.

In the preceding paper the existence of pronounced DNA tertiary structure components has been postulated for (two species of) natural eukaryotic DNA. The resulting model suggests local apparent solenoid-related DNA tertiary structure components at high sodium ion concentration c_s, partly/totally molten out at 45/60 C. With decreasing c_s the tertiary structure components have been found to be gradually reduced, at least below c_s = 0.010 M, as titration viscometrically revealed by a gradual rise of the apparent DNA contour length (Reinert et al., JBSD 9, 537, 1991).

Hence, we performed titration viscometric analyses about Nt interaction with calf thymus DNA (ctDNA) at c_s = 0.075 M, 0.010 M and 0.004 M Na⁺. The concomitant DNA conformational changes are quantitatively described in terms of the relative changes of both DNA persistence length and hydrodynamically operative apparent DNA contour length for the three first resolved interaction modes below a Nt/DNA-P ratio of 0.03.

These experiments, together with previous respective analyses at c_s = 0.20 M Na⁺ and different temperatures (I.e.), suggest that those DNA sites binding Nt most strongly predominantly are responsible for the formation of solenoid-related DNA tertiary structure components. Most probably these are A tract-containing sequences. As the essential factor for their apparent elongation effect at low Na⁺ concentrations, a gradual alteration of the number of base pairs per helix turn seems to occur below c_s = 0.010 M Na⁺ and, concomitantly, a change in phasing between intrinsic helix bends and helix screw.     

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes.

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.     

Sequence-dependent anisotropic flexibility of B-DNA. A conformational study

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)2, d(TTAA)2, d(GGGG):d(CCCC), d(GGCC)2 and d(CCGG)2,. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degree, are in agreement with experimental value of the DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6-12 degree into this groove. The above inequality is caused by stacking interaction of the bases. The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20 degree. the backbone dihedral angles vary by no more than 15 degree. The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

HMGB binding to DNA: single and double box motifs

High mobility group (HMG) proteins are nuclear proteins believed to significantly affect DNA interactions by altering nucleic acid flexibility. Group B (HMGB) proteins contain HMG box domains known to bind to the DNA minor groove without sequence specificity, slightly intercalating base pairs and inducing a strong bend in the DNA helical axis. A dual-beam optical tweezers system is used to extend double-stranded DNA (dsDNA) in the absence as well as presence of a single box derivative of human HMGB2 [HMGB2(box A)] and a double box derivative of rat HMGB1 [HMGB1(box A+box B)]. The single box domain is observed to reduce the persistence length of the double helix, generating sharp DNA bends with an average bending angle of 99 ± 9° and, at very high concentrations, stabilizing dsDNA against denaturation. The double box protein contains two consecutive HMG box domains joined by a flexible tether. This protein also reduces the DNA persistence length, induces an average bending angle of 77 ± 7°, and stabilizes dsDNA at significantly lower concentrations. These results suggest that single and double box proteins increase DNA flexibility and stability, albeit both effects are achieved at much lower protein concentrations for the double box. In addition, at low concentrations, the single box protein can alter DNA flexibility without stabilizing dsDNA, whereas stabilization at higher concentrations is likely achieved through a cooperative binding mode.     

Selective DNA bending by a variety of bZIP proteins.

We have investigated DNA bending by bZIP family proteins that can bind to the AP-1 site. DNA bending is widespread, although not universal, among members of this family. Different bZIP protein dimers induced distinct DNA bends. The DNA bend angles ranged from virtually 0 to greater than 40 degrees as measured by phasing analysis and were oriented toward both the major and the minor grooves at the center of the AP-1 site. The DNA bends induced by the various heterodimeric complexes suggested that each component of the complex induced an independent DNA bend as previously shown for Fos and Jun. The Fos-related proteins Fra1 and Fra2 bent DNA in the same orientation as Fos but induced smaller DNA bend angles. ATF2 also bent DNA toward the minor groove in heterodimers formed with Fos, Fra2, and Jun. CREB and ATF1, which favor binding to the CRE site, did not induce significant DNA bending. Zta, which is a divergent member of the bZIP family, bent DNA toward the major groove. A variety of DNA structures can therefore be induced at the AP-1 site through combinatorial interactions between different bZIP family proteins. This diversity of DNA structures may contribute to regulatory specificity among the plethora of proteins that can bind to the AP-1 site.     

MPD and DNA Bending in Crystals and in Solution

Bending of 15 to 24° is observed within crystal structures ofB-DNA duplexes, is strongly sequence-dependent, and exhibits no correlation with the concentration of MPD (2-methyl-2,4-pentanediol) in the crystallizing solution. Two types of bends are observed: facultative bends or flexible hinges at junctions between regions of G·C and A·T base-pairs, and a persistent and almost obligatory bend at the center of the sequence R-G-C-Y. Only A-tracts are characteristically straight and unbent in every crystal structure examined to date. A detailed examination of normal vector plots for individual strands of a double helix provides an explanation, in terms of the stacking properties of guanine and adenine bases. The effect of high MPD concentrations, in both solution and crystal, is to decrease local bending somewhat without removing it altogether. MPD gel retardation experiments provide no basis for choosing among the three models that seek to explain macroscopic curvature of DNA by means of microscopic bending: junction bending, bent A-tracts, or bent general- sequence DNA. Crystallographic data on the straightness of A-tracts, the bendability of non-A sequences, and the identity of inclination angles in A-tract and non-A-tractB-DNA support only the general-sequence bending model. The pre-melting transition observed in A-tract DNA probably represents a relaxation of stiff adenine stacks to a flexible conformation more typical of general-sequence DNA.     

Modeling helix-turn-helix protein-induced DNA bending with knowledge-based distance restraints.

A crucial element of many gene functions is protein-induced DNA bending. Computer-generated models of such bending have generally been derived by using a presumed bending angle for DNA. Here we describe a knowledge-based docking strategy for modeling the structure of bent DNA recognized by a major groove-inserting alpha-helix of proteins with a helix-turn-helix (HTH) motif. The method encompasses a series of molecular mechanics and dynamics simulations and incorporates two experimentally derived distance restraints: one between the recognition helix and DNA, the other between respective sites of protein and DNA involved in chemical modification-enabled nuclease scissions. During simulation, a DNA initially placed at a distance was "steered" by these restraints to dock with the binding protein and bends. Three prototype systems of dimerized HTH DNA binding were examined: the catabolite gene activator protein (CAP), the phage 434 repressor (Rep), and the factor for inversion stimulation (Fis). For CAP-DNA and Rep-DNA, the root mean square differences between model and x-ray structures in nonhydrogen atoms of the DNA core domain were 2.5 A and 1.6 A, respectively. An experimental structure of Fis-DNA is not yet available, but the predicted asymmetrical bending and the bending angle agree with results from a recent biochemical analysis.     

Analysis of birefringence decay profiles for nucleic acid helices possessing bends: the tau-ratio approach.

For nucleic acid helices in the 100-200-bp range, a central bend or point of flexibility increases the rate of rotational diffusion. In a transient electric birefringence (TEB) experiment, this increase is manifest as a reduction in the terminal (slowest) birefringence decay time. Previous experimental and theoretical work has demonstrated that the ratio of the decay times for a bent/flexible molecule and its fully duplex (linear) counterpart represents a sensitive, quantifiable measure of the apparent bend angle (tau-ratio approach). In the current work, we have examined the influence of helix parameters (e.g., persistence length, helix rise, diameter) on the tau-ratio for a given bend. The tau-ratio is found to be remarkably insensitive to variations and/or uncertainties in the helix parameters, provided that one employs bent and control molecules with the same sequence and length (apart from the bend itself). Although a single tau-ratio determination normally does not enable one to distinguish between fixed and flexible bends, such a distinction can be made from a set of tau-ratios for molecules possessing two variably phased bends. A number of additional uncertainties are examined, including errors in the estimation of the dimensions of nonhelix elements that are responsible for bends; such errors can, in principle, be estimated by performing a series of measurements for molecules of varying length.     

Sequence-Dependent Anisotropic Flexibility of B-DNA A Conformational Study

Abstract

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)₂, d(TTAA) ₂, d(GGGG):d(CCCC), d(GGCC) ₂ and d(CCGG) ₂. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5–7°, are in agreement with experimental value of the DNA persistence length.

Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6–12° into this groove. The above inequality is caused by stacking interaction of the bases.

The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20°, the backbone dihedral angles vary by no more than 15°.

The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of “mini-kinks” separated by a half-pitch of the double helix, i.e. by 5–6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimentional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

Bent DNA structures associated with several origins of replication are recognized by a unique enzyme from trypanosomatids.

Sequence-directed bending of the DNA double helix is a conformational variation found in both prokaryotic and eukaryotic organisms. The utilization of bent DNA structures from various sources as specific signals recognized by an enzyme is demonstrated here using a unique endonuclease purified from trypanosomatid cells. Crithidia fasciculata nicking enzyme was previously shown to recognize specifically the bent structure found in kinetoplast DNA minicircles. The binding constant measured for this specific interaction is of two orders of magnitude higher than that measured for the binding of the enzyme to a non-curved sequence. As determined by binding competition and mobility shift electrophoresis analyses, this enzyme recognizes the sequence-directed bends associated with the origins of replication of bacteriophage lambda and simian virus 40 (SV40), as well as that located within the autonomously replicating sequence (ARS1) region of the yeast S. cerevisiae.     

Bending of synthetic bacteriophage 434 operators by bacteriophage 434 proteins.

The extent of DNA bending induced by 434 repressor, its amino terminal DNA binding domain (R1-69), and 434 Cro was studied by gel shift assay. The results show that 434 repressor and R1-69 bend DNA to the same extent. 434 Cro-induced DNA bends are similar to those seen with the 434 repressor proteins. On approximately 265 base pair fragments, the cyclic AMP receptor protein of Escherichia coli (CRP) produces larger mobility shifts than does 434 repressor. This indicates that the 434 proteins bend DNA to a much smaller extent than does CRP. The effects of central operator sequence on intrinsic and 434 protein-induced DNA bending was also examined by gel shift assay. Two 434 operators having different central sequences and affinities for 434 proteins display no static bending. The amount of gel shift induced by 434 repressor on these operators is identical, showing that the 434 repressor bends operators with different central sequences to the same extent. Hence, mutations in the central region of the operator do not influence the bent structure of the unbound or bound operator.     

Interaction of netropsin and distamycin with deoxyribonucleic acid: electric dichroism study

We report dichroism and equilibrium binding studies of netropsin (Net) and distamycin A3 (Dist) binding to deoxyribonucleic acid (DNA). We show that at low degrees of binding (r) to calf thymus DNA, Net induces a considerable increase in the apparent DNA length (14 A/drug molecule bound), closely analogous to the results reported earlier for Dist. In addition, we show that chicken erythrocyte DNA shows length changes similar to those of calf thymus DNA upon distamycin binding. DNA length reaches a maximum at 1 bound drug/20-30 base pairs and then decreases to its initial value by r = 0.1. This effect is not seen for two other DNAs with nearly identical A + T base pair content and may therefore arise from the details of base sequence or base modification in eukaryotic DNA. We also show that Dist binding to calf thymus DNA at low r values is positively cooperative and shows a DNA affinity which is primarily nonionic. We demonstrate that independent of the DNA to which they are bound, the Net and Dist transition moments are inclined by 43 +/- 3 degrees from the helix axis, consistent with the idea that both drugs bind inside and parallel to the DNA small groove. From dichroism measurements, we show that the conformational change induced in calf thymus DNA by Dist does not kink or bend the helix and does not substantially alter the average inclination of the bases. Finally, we outline a statistical mechanical theory for calculation of binding isotherms when binding is coupled to a DNA structural change.     

Contribution of the intrinsic curvature to measured DNA persistence length

The persistence length of DNA, a, depends both on the intrinsic curvature of the double helix and on the thermal fluctuations of the angles between adjacent base-pairs. We have evaluated two contributions to the value of a by comparing measured values of a for DNA containing a generic sequence and for an "intrinsically straight" DNA. In each 10 bp segment of the intrinsically straight DNA an initial sequence of five bases is repeated in the sequence of the second five bases, so any bends in the first half of the segment are compensated by bends in the opposite direction in the second half. The value of a for the latter DNA depends, to a good approximation, on thermal fluctuations only; there is no intrinsic curvature. The values of a were obtained from measurements of the cyclization efficiency for short DNA fragments, about 200 bp in length. This method determines the persistence length of DNA with exceptional accuracy, due to the very strong dependence of the cyclization efficiency of short fragments on the value of a. We find that the values of a for the two types of DNA fragment are very close and conclude that the contribution of the intrinsic curvature to a is at least 20 times smaller than the contribution of thermal fluctuations. The relationship between this result and the angles between adjacent base-pairs, which specify the intrinsic curvature, is analyzed.