Magnitude and direction of DNA bending induced by screw-axis orientation: influence of sequence, mismatches and abasic sites

DNA-bending flexibility is central for its many biological functions. A new bending restraining method for use in molecular mechanics calculations and molecular dynamics simulations was developed. It is based on an average screw rotation axis definition for DNA segments and allows inducing continuous and smooth bending deformations of a DNA oligonucleotide. In addition to controlling the magnitude of induced bending it is also possible to control the bending direction so that the calculation of a complete (2-dimensional) directional DNA-bending map is now possible. The method was applied to several DNA oligonucleotides including A(adenine)-tract containing sequences known to form stable bent structures and to DNA containing mismatches or an abasic site. In case of G:A and C:C mismatches a greater variety of conformations bent in various directions compared to regular B-DNA was found. For comparison, a molecular dynamics implementation of the approach was also applied to calculate the free energy change associated with bending of A-tract containing DNA, including deformations significantly beyond the optimal curvature. Good agreement with available experimental data was obtained offering an atomic level explanation for stable bending of A-tract containing DNA molecules. The DNA-bending persistence length estimated from the explicit solvent simulations is also in good agreement with experiment whereas the adiabatic mapping calculations with a GB solvent model predict a bending rigidity roughly two times larger.     

Ring closure probabilities for DNA fragments by Monte Carlo simulation

The rate of ligation of DNA molecules into circular forms depends on the ring closure probability, commonly called the j-factor, which is a sensitive measure of the extent to which thermal fluctuations contribute to bending and twisting of DNA molecules in solution. We present a theoretical treatment of the cyclization equilibria of DNA that employs a special Monte Carlo method for generating large ensembles of model DNA chains. Using this method, the chain length dependence of the j-factor was calculated for molecules. in the size range 250 to 2000 base-pairs. The Monte Carlo results are compared with recent analytical theory and experimental data. We show that a value of 475 A for the persistence length of DNA, close to values measured by a number of other methods, is in excellent agreement with the cyclization results. Preliminary applications of the Monte Carlo method to the problem of systematically bent DNA molecules are presented. The calculated j-factor is shown to be very sensitive to the amount of bending in these fragments. This fact suggests that ligase closure measurements of systematically bent DNA molecules should be a useful method for studying sequence-directed bending in DNA.     

Higher-order structure of long repeat chromatin.

The higher-order structure of chromatin isolated from sea urchin sperm, which has a long nucleosomal DNA repeat length (approximately 240 bp), has been studied by electron microscopy and X-ray diffraction. Electron micrographs show that this chromatin forms 300 A filaments which are indistinguishable from those of chicken erythrocytes (approximately 212 bp repeat); X-ray diffraction patterns from partially oriented samples show that the edge-to-edge packing of nucleosomes in the direction of the 300 A filament axis, and the radial disposition of nucleosomes around it, are both similar to those of the chicken erythrocyte 300 A filament, which is described by the solenoid model. The invariance of the structure with increased linker DNA length is inconsistent with many other models proposed for the 300 A filament and, furthermore, means that the linker DNA must be bent. The low-angle X-ray scattering in the 300-400 A region both in vitro and in vivo differs from that of chicken erythrocyte chromatin. The nature of the difference suggests that 300 A filaments in sea urchin sperm in vivo are packed so tightly together that electron-density contrast between individual filaments is lost; this is consistent with electron micrographs of the chromatin in vitro.     

Flow birefringence of T7 phage DNA: Dependence on salt concentration

We have determined extinction angles and flow birefringence of T7 bacteriophage DNA over a wide range of shear, polymer concentration, and solvent ionic strength. From these data, information on the simple salt dependence of coil permeability to solvent and on short-range intrachain interactions (persistence length) was obtained. At all ionic strengths, our results are consistent with a partially draining coil in the Gaussian subchain dynamical theory of Rouse-Zimm-Tschoegl-Bloomfield. Salt dependence of persistence length is comparable to, although somewhat less than, that obtained previously using similar methods with a fivefold higher-molecular-weight DNA (T2 bacteriophage DNA). Possible reasons for observed discrepancies are analyzed, and the results of this work are compared in detail to other current studies of solvent ionic strength dependence in persistence length and hydrodynamic properties of DNA.     

Light-scattering studies on deoxyribonucleic acid flexibility. The solution properties of a small circular deoxyribonucleic acid molecule

Previous investigations on the persistence length of DNA in solution have revealed large discrepancies between hydrodynamic results and those from light-scattering techniques which have potentially a greater resolving power. The information obtained from experiments on a small circular DNA molecule has resolved these discrepancies. The non-superhelical circular double-stranded DNA molecule from bacteriophage [unk]X174-infected cells is small enough to permit accurate light-scattering extrapolations, and its solutions have negligible anisotropy. The persistence length obtained from experimental investigations on this molecule is comparable with that obtained by hydrodynamic techniques, even with variation of the excluded-volume factor.     

Probing the Elasticity of DNA on Short Length Scales by Modeling Supercoiling under Tension

The wormlike-chain (WLC) model is widely used to describe the energetics of DNA bending. Motivated by recent experiments, alternative, so-called subelastic chain models were proposed that predict a lower elastic energy of highly bent DNA conformations. Until now, no unambiguous verification of these models has been obtained because probing the elasticity of DNA on short length scales remains challenging. Here we investigate the limits of the WLC model using coarse-grained Monte Carlo simulations to model the supercoiling of linear DNA molecules under tension. At a critical supercoiling density, the DNA extension decreases abruptly due to the sudden formation of a plectonemic structure. This buckling transition is caused by the large energy required to form the tightly bent end-loop of the plectoneme and should therefore provide a sensitive benchmark for model evaluation. Although simulations based on the WLC energetics could quantitatively reproduce the buckling measured in magnetic tweezers experiments, the buckling almost disappears for the tested linear subelastic chain model. Thus, our data support the validity of a harmonic bending potential even for small bending radii down to 3.5 nm.     

DNA multimode interaction with berenil and pentamidine; double helix stiffening, unbending and bending

The antitrypanosomal drugs berenil (Ber) and pentamidine (Pm) preferentially bind to DNA in the minor groove of A.T-rich domains. The properties of A.T clusters are essential for sequence-mediated helix bending. Groove binding drugs locally stiffen the DNA helix but may also change intrinsic helix bends or may bend straight DNA. Ligand binding to randomly distributed sites alters the apparent DNA persistence length, a. Criteria permit the distinction of the underlying mechanism(s). Helix bends, if phased with the helix screw, however, generate solenoidal superhelix components mediating an apparent change of the hydrodynamically effective DNA contour length, L. The measurement of relative changes of both, a and L, as induced by Ber or Pm is performed by titration rotational viscometry. The determination of the two quantities requires two independent measurements: the relative change of DNA intrinsic viscosity, deltay, for short (tending to rod-like) DNA molecules and for comparably long (almost coil-like) ones as a function of r, the bound drug molecules per DNA-P, and this under conditions effectively excluding intramolecular DNA-DNA crosslinking effects. At least at r< or =0.05 and < or =0.03, respectively, the two drugs virtually bind completely to a eukaryotic DNA. r ranges of different drug binding strength and, concomitantly, of different specific conformational response, could be resolved. They represent (sub)modes of different DNA sequences... Whereas the mode-specific elongation effects are fairly similar for both systems, there are pronounced quantitative differences in the relative change of DNA persistence length. The sites of highest Ber-binding strength are correlated to unbent alternating helical A.T segments followed by bent and by less bent or unbent dAn.dTn tracts straightened on Ber-binding. For Pm-DNA interaction the ligand bends the sites of highest Pm affinity. Generally, ligand induced and sequence mediated local DNA-bend removal or DNA bending, as observed for several modes of interaction with A.T rich DNA, are considered to be of gene regulatory relevance.     

Atomic force microscopy of long and short double-stranded, single-stranded and triple-stranded nucleic acids.

Atomic force microscopy (AFM, also called scanning force microscopy) is proving to be a useful technique for imaging DNA. Thus it is important to push the limits of AFM imaging in order to explore both what types of DNA can be reliably imaged and identified and also what substrates and methods of sample preparation are suitable. The following advances in AFM of DNA are presented here. (i) DNA molecules as short as 25 bases can be seen by AFM. The short single-stranded DNAs imaged here (25 and 50 bases long) appeared globular in the AFM, perhaps because they are all capable of intramolecular base pairing and because the DNAs were in a Mg(ll) buffer, which facilitates intramolecular cross-bridging. (ii) AFM images in air of short double-stranded DNA molecules, 100-200 bp, gave lengths consistent with A-DNA. (iii) AFM images of poly (A) show both short bent lumpy molecules with an apparent persistence length of 40 nm and long straight molecules with an apparent persistence length of 600 nm. For comparison, the apparent persistence length for double-stranded DNA from phX-174 under the same conditions was 80 nm. (iv) Structures believed to be triple- stranded DNA were seen in samples of poly(dA.poly(dT) and poly (dG).poly(dC). These structures were twice as high as double-stranded DNA and the same width. (v) Entire molecules of lambda DNA, approx. 16 micron long, were imaged clearly in overlapping scans. (vi) Plasmid DNA was imaged on oxidized silicon, although less clearly than on mica.     

Electrostatic contribution to the bending of DNA

A model is derived that accounts for the short-range electrostatic contribution to the bending of DNA molecule in solution and in complexes with proteins in terms of the non-linear Poisson-Boltzmann equation. We defined that the short-range electrostatic interactions depend on the changes of the polyion surface charge density under deformation, while the long-range interactions depend on the bending-induced changes in distances between each two points along the polyion axis. After an appropriate simplification of the Poisson-Boltzmann equation, the short-range term is calculated separately giving the lower limit for the electrostatic contribution to the DNA persistence length. The result is compared with the theoretical approaches developed earlier [M. Fixman, J. Chem. Phys. 76 (1982) 6346; M. Le Bret, J. Chem. Phys. 76 (1982) 6243] and with the experimental data. The conclusion is made that the results of Fixman-Le Bret, which took into account both types of the electrostatic interactions for a uniformly bent polyion, give the upper limit for the electrostatic persistence length at low ionic strength, and the actual behavior of the DNA persistence length lies between two theoretical limits. Only the short-range term is significant at moderate-to-high ionic strength where our results coincide with the predictions of Fixman-Le Bret. The bending of DNA on the protein surface that is accompanied by an asymmetric neutralization of the DNA charge is also analyzed. In this case, the electrostatic bending energy gives a significant favorite contribution to the total bending energy of DNA. Important implications to the mechanisms of DNA-protein interactions, particularly in the nucleosome particle, are discussed.     

Mitotic chromosome-binding activity of latency-associated nuclear antigen 1 is required for DNA replication from terminal repeat sequence of Kaposi's sarcoma-associated herpesvirus

Latency-associated nuclear antigen 1 (LANA1) of Kaposi's sarcoma-associated herpesvirus (KSHV) is implicated in the persistence of the viral genome during latent infection. It has been suggested that LANA1 tethers the viral genome to the host chromosome and also participates actively in DNA replication from the terminal repeat of KSHV. Here we show by mutational analysis that the mitotic chromosome-binding activity of LANA1 is tightly coupled to its replication activity. Thus, KSHV appears to have evolved a unique tactic for its stable maintenance.     

Theory of H1-mediated control of higher orders of structure in chromatin

It is known that the lysine-rich histone H1 induces both higher orders of folding in chromatin and donut shapes in DNA. However, these phenomena occur only on the high-salt side of a narrow transition range located at about 0.02M salt. Previous theoretical analyses of the ionic-strength dependencies of DNA persistence length and denaturation rate have provided the information that the basic rigid-rod unit in high-molecular-weight DNA is a segment 60 base pairs in length and that if the phosphate charge is neutralized, this segment will spontaneously adopt a bent conformation with radius of curvature 170 Å. On the assumption that an H1 molecule does not completely neutralize the DNA charge in its vicinity, the theory has been extended here to determine the onset of spontaneous bending as a function of salt concentration and extent of phosphate neutralization. A salt transition of the kind observed has been found for the realistic value of 82% charge neutralization, with the actual value likely to be in the neighborhood of 90%, as suggested by the measurements of Wilson and Bloomfield.¹ It is recalled that the spacer DNA length in chromatin is of about the same length as the DNA rigid-rod unit. If binding of H1 to the spacer induces, as predicted, a bent conformation of radius about 170 Å, then the observed value of about 150 Å for the outer radius of the solenoid presently thought to be the basic mode of folding for a nucleosome chain can be understood as a reflection of the inherent maximum curvature of DNA in aqueous salt solution.     

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.     

Analysis of the intrinsic bend in the M13 origin of replication by atomic force microscopy

Atomic force microscopy (AFM) has been used to image a 471-bp bent DNA restriction fragment derived from the M13 origin of replication in plasmid LITMUS 28, and a 476-bp normal, unbent fragment from plasmid pUC19. The most probable angle of curvature of the 471-bp DNA fragment is 40-50 degrees, in reasonably good agreement with the bend angle determined by transient electric birefringence, 38 degrees +/- 7 degrees. The normal 476-bp DNA fragment exhibited a Gaussian distribution of bend angles centered at 0 degrees, indicating that this fragment does not contain an intrinsic bend. The persistence length, P, was estimated to be 60 +/- 8 and 62 +/- 8 nm for the 471- and 476-bp fragments, respectively, from the observed mean-square end-to-end distances in the AFM images. Since the P-values of the normal and bent fragments are close to each other, the overall flexibility of DNA fragments of this size is only marginally affected by the presence of a stable bend. The close agreement of AFM and transient electric birefringence results validates the suitability of both methods for characterizing DNA bending and flexibility.     

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).     

The possibility of intrinsic local curvature in DNA toroids

A variety of solution conditions are known to induce collapse of linear DNA into a compact configuration without dramatic change of local structure. When visualized, these compact forms frequently have a toroidal appearance. We ask whether the molecular basis of the toroidal shape can be a stable curvature of isolated rodlike DNA segments. Application of the classical Euler-Lagrange theory for the buckling of elastic rods provides us with an affirmative answer. Specifically, we see that, in principle, sufficient addition of inert polymer to DNA solutions can induce buckling of DNA segments. However, no attempt is made to quantitate the Euler-Lagrange condition for sufficiency in terms of added polymer concentration. We find also that complete neutralization of the phosphate charge is more than sufficient to induce buckling of DNA segments of length comparable to a Kuhn segment. The quantitative argument involves comparison of buckling forces provided by polyelectrolyte theory with the Euler-Lagrange criterion. Knowledge of the ionic-strength dependence of DNA bending stiffness (persistence length) is not required.     

Use of isogenic lines and simultaneous probing to identify DNA markers tightly linked to the tm-2a gene in tomato

The Tm-2a gene of tomato confers resistance to the viral pathogen, tobacco mosaic virus. Like many economically important plant genes, Tm-2a has been characterized phenotypically and by classical linkage analysis, yet nothing is known about its gene product. We report here the isolation of two DNA clones which are very tightly linked to the Tm-2a gene. These clones were identified by testing 122 genomic clones as hybridization probes against Southern blots consisting of DNA from pairs of nearly isogenic lines with or without the Tm-2a gene. Screening such a large number of clones in a short period of time was facilitated by co-labeling and simultaneous probing of sets of up to 10 random genomic clones. Tightly linked clones were distinguished by the fact that they exhibited one or more restriction fragment length polymorphisms between the nearly isogenic lines. Tight linkage of the clones with Tm-2a was verified in a segregating F(2) population. Both mapped to the same locus 0.4 +/- 0.4 centimorgans away from Tm-2a and may provide starting points for a genomic ;;walk' to this gene. Due to the availability of isogenic lines in many plant species, the strategy outlined in this paper should be widely applicable for selecting DNA clones tightly linked to genes of interest.