Breathing and bending fluctuations in DNA modeled by an open-base-pair kink coupled to axial compression

Evaluating the role of coherent delocalized phonon-like modes in DNA cyclization

The innate flexibility of a DNA sequence is quantified by the Jacobson-Stockmayer's J-factor, which measures the propensity for DNA loop formation. Recent studies of ultra-short DNA sequences revealed a discrepancy of up to six orders of magnitude between experimentally measured and theoretically predicted J-factors. These large differences suggest that, in addition to the elastic moduli of the double helix, other factors contribute to loop formation. Here, we develop a new theoretical model that explores how coherent delocalized phonon-like modes in DNA provide single-stranded "flexible hinges" to assist in loop formation. We combine the Czapla-Swigon-Olson structural model of DNA with our extended Peyrard-Bishop-Dauxois model and, without changing any of the parameters of the two models, apply this new computational framework to 86 experimentally characterized DNA sequences. Our results demonstrate that the new computational framework can predict J-factors within an order of magnitude of experimental measurements for most ultra-short DNA sequences, while continuing to accurately describe the J-factors of longer sequences. Further, we demonstrate that our computational framework can be used to describe the cyclization of DNA sequences that contain a base pair mismatch. Overall, our results support the conclusion that coherent delocalized phonon-like modes play an important role in DNA cyclization.

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.

In the absence of writhe, DNA relieves torsional stress with localized, sequence-dependent structural failure to preserve B-form

To understand how underwinding and overwinding the DNA helix affects its structure, we simulated 19 independent DNA systems with fixed degrees of twist using molecular dynamics in a system that does not allow writhe. Underwinding DNA induced spontaneous, sequence-dependent base flipping and local denaturation, while overwinding DNA induced the formation of Pauling-like DNA (P-DNA). The winding resulted in a bimodal state simultaneously including local structural failure and B-form DNA for both underwinding and extreme overwinding. Our simulations suggest that base flipping and local denaturation may provide a landscape influencing protein recognition of DNA sequence to affect, for examples, replication, transcription and recombination. Additionally, our findings help explain results from single-molecule experiments and demonstrate that elastic rod models are strictly valid on average only for unstressed or overwound DNA up to P-DNA formation. Finally, our data support a model in which base flipping can result from torsional stress.

DNA A-tracts bending: polarization effects on electrostatic interactions across their minor groove

Bending by the DNA A-tracts constitutes a contentious issue, suggesting deficiencies in the physics employed so far. Here, we inquire as to the importance in this bending of many-body polarization effects on the electrostatic interactions across their narrow minor groove. We have done this on the basis of the findings of Jarque and Buckingham who developed a procedure based on a Monte Carlo simulation for two charges of the same sign embedded in a polarizable medium. Remarkably, the present analysis reveals that for compact DNA conformations, which result from dynamic effects, an overall attractive interaction operates between the phosphate charges; this interaction is especially strong for the narrow minor groove of the A-tracts, suggesting a tendency for DNA to bend toward this groove. This tendency is in agreement with the conclusions of electrophoretic and NMR solution studies. The present analysis is also consistent with the experimental observations that the minor groove is much more easily compressible than the major groove and the bending propensity of the A-tracts is greatly reduced at "premelting" temperatures. By contrast, the dielectric screening model predicts a repulsion between the phosphate charges and is not consistent with the aforementioned bending tendency or experimental observations.

Solvation of nucleosides in aqueous mixtures of organic solvents: relevance to DNA open basepairs

Toward the goal of understanding how open basepairs in DNA interact with their heterogeneous environment, we have studied the steady-state intrinsic fluorescence properties of the purine and pyrimidine deoxynucleosides in organic solvents in the presence of small amounts of water. The organic solvents used in the present study were: n-butanol, acetonitrile, methanol, n-propanol, isopropanol, and isobutanol. For n-butanol and acetonitrile, which have a high degree of amphiphilicity and weak hydrogen bonding ability, respectively, the fluorescence spectral properties of the purines are found to depend on the sequence of steps in which the aqueous mixtures were formed. By contrast, no such dependence was observed in the mixtures with any of the other solvents used in the present study. Moreover, no such dependence was observed for the pyrimidines. These findings suggest that the final solvation network around the purines is dependent on the nature of the environment to which they were initially exposed. This would tend to present an impediment to the closing of AT or GC basepairs in DNA that become open as a result of structural fluctuations, DNA bending, or protein-DNA interactions.

The trapping of a spontaneously "flipped-out" base from double helical nucleic acids by host-guest complexation with beta-cyclodextrin: the intrinsic base-flipping rate constant for DNA and RNA

Beta-cyclodextrin, which forms stable host-guest complexes with purine bases, induces the melting of RNA and DNA duplexes below their normal melting temperatures. Alpha-cyclodextrin, which does not form stable complexes, has no effect on either RNA or DNA. Gamma-cyclodextrin, which forms weaker complexes, has no effect on RNA and a smaller effect than beta-cyclodextrin on DNA. The rate of melting is kinetically first-order in duplex and, above about 20 mM beta-cyclodextrin, is independent of the beta-cyclodextrin concentration with a first-order rate constant, common to both RNA and DNA, of (3.5 +/- 0.5) x 10(-3) s(-1) at 61 degrees C (DNA) and at 50 degrees C (RNA). This is taken to be the rate constant for spontaneous "flipping out" of a base from within the duplex structure of the nucleic acids, the exposed base being rapidly trapped by beta-cyclodextrin. Like beta-cyclodextrin, nucleic acid methyltransferases bind the target base for methylation in a site that requires it to have flipped out of its normal position in the duplex. The spontaneous flip-out rate constant of around 10(-3) s(-1) is near the value of k(cat) for the methyltransferases (ca. 10(-3) to 10(-1) s(-1)). In principle, the enzymes, therefore, need effect little or no catalysis of the flipping-out reaction. Nevertheless, the flip-out rate in enzyme/DNA complexes is much faster. This observation suggests that the in vivo circumstances may differ from in vitro models or that factors other than a simple drive toward higher catalytic power have been influential in the evolution of these enzymes.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

Studies of DNA bending and flexibility using gel electrophoresis

Gel electrophoretic methods have become established as primary tools in the study and elucidation of sequence-directed curvature both in free DNA and in the operator DNA of several site-specific nucleoprotein complexes. Results using them have been generally consistent with physical methods sensitive to DNA structure and conformation in those instances where direct comparisons can be made, and in a number of cases, gel methods have provided unique information not presently available from other techniques. Two basic strategies have been used: one based upon anomalous gel mobility effects; and a second based upon cyclization properties of curved DNA. Within each of these categories, various approaches have been used, some of which can lead, in favorable cases, to quantitative estimation of bending angles. In this review, the various gel-based methods that have been used to date are critically discussed and the qualitative and quantitative information that can be obtained from them is evaluated. A number of possible structural models for DNA curvature are described and a distinction is drawn between static or fixed bending and bending due to anisotropic flexibility at specific sequence loci. The importance and roles of gel electrophoretic methods in providing experimental approaches to this question are discussed. It is suggested that both static curvature and anisotropic flexibility in operator DNA may provide much of the basis for indirect readout of sequence information by specific site-binding regulatory proteins.

Thermodynamic characterization of daunomycin-DNA interactions: comparison of complete binding profiles for a series of DNA host duplexes

Using a combination of spectroscopic and calorimetric techniques, we have determined complete thermodynamic binding profiles (delta G degree, delta H degree, and delta S degree) for the complexation of daunomycin to a series of 10 polymeric DNA duplexes. We find the resulting drug binding data to be sensitive to the base composition and sequence of the host duplex, with the binding free energies ranging from -7.5 to -10.8 kcal/mol of bound drug and the binding enthalpies ranging from +4.11 to -10.76 kcal/mol of bound drug at 25 degrees C. The smaller range in the free energy term reflects the impact of large enthalpy-entropy compensations. We observe that the three synthetic duplexes which exhibit the highest daunomycin binding affinities all contain GC (or IC) base pairs as part of alternating purine/pyrimidine sequence motifs, with these high binding affinities being strongly enthalpy driven at 25 degrees C. Specific comparisons between the binding profiles for daunomycin complexation with select pairs of host duplexes lead to the following observations: (1) The presence or absence of a major-groove methyl group does not alter daunomycin binding thermodynamics. (2) The presence or absence of a minor-groove amino group does alter daunomycin binding thermodynamics. (3) Duplexes with different base compositions but identical minor-groove functionality exhibit similar daunomycin binding thermodynamics. (4) Homopolymeric duplexes composed of either AT or AU base pairs, but not GC base pairs, exhibit large enthalpy-entropy compensations in their daunomycin binding profiles. We propose interpretations of these and other features of our thermodynamic data in terms of specific daunomycin-DNA interactions deduced from available structural data.

Atomic structure of the DNA repair [4Fe-4S] enzyme endonuclease III

The crystal structure of the DNA repair enzyme endonuclease III, which recognizes and cleaves DNA at damaged bases, has been solved to 2.0 angstrom resolution with an R factor of 0.185. This iron-sulfur [4Fe-4S] enzyme is elongated and bilobal with a deep cleft separating two similarly sized domains: a novel, sequence-continuous, six-helix domain (residues 22 to 132) and a Greek-key, four-helix domain formed by the amino-terminal and three carboxyl-terminal helices (residues 1 to 21 and 133 to 211) together with the [4Fe-4S] cluster. The cluster is bound entirely within the carboxyl-terminal loop with a ligation pattern (Cys-X6-Cys-X2-Cys-X5-Cys) distinct from all other known [4Fe-4S] proteins. Sequence conservation and the positive electrostatic potential of conserved regions identify a surface suitable for binding duplex B-DNA across the long axis of the enzyme, matching a 46 angstrom length of protected DNA. The primary role of the [4Fe-4S] cluster appears to involve positioning conserved basic residues for interaction with the DNA phosphate backbone. The crystallographically identified inhibitor binding region, which recognizes the damaged base thymine glycol, is a seven-residue beta-hairpin (residues 113 to 119). Location and side chain orientation at the base of the inhibitor binding site implicate Glu112 in the N-glycosylase mechanism and Lys120 in the beta-elimination mechanism. Overall, the structure reveals an unusual fold and a new biological function for [4Fe-4S] clusters and provides a structural basis for studying recognition of damaged DNA and the N-glycosylase and apurinic/apyrimidinic-lyase mechanisms.

Effects of Na+ on the persistence length and excluded volume of T7 bacteriophage DNA

Total intensity, Rayleigh light scattering has been used to measure the rms radius, second virial coefficient, persistence length, and excluded volume of homogeneous T7 bacteriophage DNA as a function of Na+ concentration (0.005 to 3.0 M). All parameters decrease sharply as [Na+] increases, and tend to level off at high Na+. The variation of persistence length with [Na+] is consistent with predictions from counterion condensation theory.

The elastic resilience of DNA can induce all-or-none structural transitions in the nucleosome core particle

DNA on the surface of the histone octamer in the native nucleosome core particle is modeled as a circumferentially wound elastic line on the surface of a cylinder. In a model for the radial transition, the line is allowed to straighten, and thus lose energy, by swinging off the surface, but it is impeded in such an excursion by a radial force field representing the attractive interaction between DNA and histone octamer. In a model for the axial transition, the line may straighten by becoming more parallel to a generator of the cylinder while remaining on the surface. In this mode of straightening, dimer-tetramer or tetramer-tetramer interfaces are disrupted, and the resulting energy gain impedes the transition. Both radial and axial transitions are predicted to occur in all-or-none fashion. We propose that these models are related to the abrupt transitions actually observed in the nucleosome core particle.

Thermodynamics of cation-induced DNA condensation

An estimation of the various free energy contributions to DNA collapse into toroidal particles is made, considering DNA bending and segment mobility, electrostatic repulsions between DNA chains, and attractive forces resulting from correlated counterion fluctuations. It is shown that the process of DNA condensation becomes spontaneous in the presence of divalent cations in methanol, and in the presence of tri- or tetravalent cations in water media. This is a consequence of the large decrease in the electrostatic repulsion between charged DNA segments, allowing the attractive force resulting from correlated fluctuations of bound counterions to become dominant. Our calculations indicate that short DNA fragments would condense into multimolecular particles in order to maximize the attractive force due to counterion fluctuations.

Internal dynamics of linear and superhelical DNA as studied by photon correlation spectroscopy

Photon correlation spectroscopy has been used to study the translational (D0) and internal dynamics of monodisperse linear and polydisperse superhelical plasmid-DNAs. Scattering vector dependent correlation functions were measured and analyzed with the inverse Laplace transform CONTIN written by S. Provencher. For scattering vectors (q) lower than 1.3 X 10(5) cm-1, D0 can be separated from internal dynamics. Using the measured D0 value, the q dependence of internal modes was determined. Extrapolation of the internal relaxation times for q----0 yields the longest internal relaxation time t1. This time can be connected to the persistence length in terms of the Berg-Soda model, which describes the molecule as a semiflexible circular polymer with hydrodynamic interactions. The calculated length of 76 nm for DNA, I = 0.15 mol/L, is a little higher than the one obtained from static light scattering data without excluded volume corrections. A comparison of experimental and simulated correlation functions for the Berg-Soda model shows that the model gives a fairly good description of the dynamics of the linear molecule, whereas large discrepancies between model and experimental functions are observed for the superhelical DNA. Small differences between model and experimental functions are mainly attributed to the neglect of the torsional modes that may be coupled to bending and flexing modes. For the superhelical DNA the agreement is improved if the calculation is carried out with a linear molecule, with shorter contour length and increased diameter. Both quantities can be derived from the known superhelix tilt angle.

Sequence-dependent kinks induced in curved DNA

In certain curved DNA fragments without AA dinucleotides, the gel retardation anomaly associated with curvature passes through a maximum with fragment length, indicating length (and electric field) dependent structural transitions in the DNA. We suggest that thermally induced stereochemical kinks in DNA are stabilized in the gel, thus relieving the effects of curvature. These kinks are shown to occur specifically at CA/TG and TA/TA stacks. Other physical and biological evidence points to frequent structural dislocations at CA and TA steps. These reversible sequence dependent kinks may therefore represent a novel class of structural protein-DNA recognition elements.

Computer simulation of hydrodynamic properties of semiflexible macromolecules: randomly broken chains, wormlike chains, and analysis of properties of DNA

The translational and rotational diffusion coefficients and the intrinsic viscosity of semiflexible, randomly broken, and wormlike chains have been obtained by Monte Carlo simulation in the context of the rigid-body treatment. Both approximate and rigorous rigid-body hydrodynamics are used, so that the error introduced by the approximate methods can be evaluated. A randomly broken chain and a wormlike chain having the same contour length and persistence length have the same radius of gyration but different values for any of the hydrodynamic properties. The two types of chains are compared in this regard. Considering that the cross section of the chain is represented by a cylinder better than by a string of spheres, we devise a cylindrical correction to be applied to the results simulated for chains of beads. Application is made to the analysis of experimental data for the translational and rotational coefficients of DNA fragments with up to 10(3) base pairs, obtaining the persistence length for each model. The values for the wormlike chain agree well with model-independent values obtained from radii of gyration and with other literature data at varying ionic strength. The randomly broken chain is equally able to reproduce the experimental length dependence of the properties, but the resulting persistence length may be too high.

DNA bending induced by covalently bound drugs. Gel electrophoresis and chemical probe studies

Modification of nucleotide residues arising from the covalent binding of a drug or as a result of irradiation with ultraviolet light can induce distortion of the DNA double helix. The purpose of this review is to show that, from investigation of the electrophoretic mobility of the modified DNA fragments, one can deduce whether the distortions behave more as the centers of directed bends or as hinge joints. It is also demonstrated that chemical probes are a complementary tool for the analysis of distortions at the nucleotide level.

Base pair opening pathways in B-DNA

Molecular modeling is used to study the opening pathways of bases within a B-DNA oligomer. It is demonstrated that many open states are possible for a single base pair, although a preference for opening towards the major groove of the double helix is found. In addition we show that opening is strongly influenced by the nature of the base involved and is also coupled in many cases to DNA bending.

The DNA bending by acetylaminofluorene residues and by apurinic sites

We have studied the distortions induced in double-stranded oligonucleotides by covalently bound acetylaminofluorene residues and by apurinic sites. Within the acetylaminofluorene-modified oligonucleotide three base-pairs are unpaired as detected by the chemical probes chloroacetaldehyde and osmium tetroxide. These two probes reveal that the bases adjacent to the apurinic site are paired. In both the modified double-stranded oligonucleotides, the backbone on the 5' side of the modification is more reactive with 1,10-phenanthroline copper than the backbone on the 3' side. On polyacrylamide gels, the ligated multimers of acetylaminofluorene or apurinic site-modified oligonucleotides migrate slower than the multimers of the unmodified oligonucleotides. It is suggested that the acetylaminofluorene-modified guanine residues and the apurinic sites behave more as hinge joints than as the centres of directed bends.

Interplay of hydrophobic and electrostatic interactions in biopolymer chromatography. Effect of salts on the retention of proteins

The effect of salt on the retention behavior of proteins in electrostatic and hydrophobic interaction chromatography is described by a three-parameter equation, log k' = A - B log ms + Cms, where k' is the retention factor and ms is the molality of the salt in the eluent. Parameter B, termed the electrostatic interaction parameter, depends on the characteristic charge of the protein and the salt counterion and governs the change of retention with the salt concentration in ion-exchange chromatography. According to the model the magnitude of the hydrophobic interaction parameter C is determined by the hydrophobic contact area upon protein binding at the stationary phase surface and the properties of the salt as measured by its molal surface tension increment. Retention data measured at different salt concentrations in the eluent on a variety of ion exchangers can be fitted to the above equation which yields U-shaped plots of log k' against log ms. The limiting slopes of the appropriate plots at sufficiently low and high salt concentrations can be used to evaluate the electrostatic and hydrophobic interaction parameters, respectively. The approach, which is based on a combination of established treatments of electrostatic and hydrophobic interactions offers a convenient framework for analyzing retention data in biopolymer high-performance liquid chromatography and for the characterization of stationary phases. Furthermore, it may facilitate some characterization of protein molecules on the basis of their retention behavior as a function of the concentration and nature of the salt in the eluent. In the treatment of electrostatic interactions use is made of the counterion condensation theory that is believed to make possible a more comprehensive analysis than the traditional stoichiometric ion-exchange model which assumes binding of the proteins by coulombic interactions at discrete sites. The treatment of hydrophobic interactions is based on an adaptation of the solvophobic theory which predicts that the hydrophobic portion of the free energy of binding is proportional to the hydrophobic contact area and the microthermodynamic surface tension of the aqueous salt solution. Despite its simplicity the theory was successful in explaining the observed effect of the nature and concentration of salt in the eluent, the pH and the effect of the density of fixed charges at the surface of the stationary phase in the absence of specific salt effects.

Different binding modes of spermine to A-T and G-C base pairs modulate the bending and stiffening of the DNA double helix

The influence of base composition (and sequence) on the process of interaction between synthetic polynucleotides and spermine, has been investigated using ultraviolet (including second derivative) spectroscopy, and electric dichroism. Different binding modes of spermine to poly(dG-dC) as compared to A-T containing polynucleotides, were evidenced. An interaction with the N7 and O6 of guanine is probably partially involved in the former case while simple electrostatic interaction with the phosphate groups would dominate in the latter. In the intermediate binding range (spermine over DNA phosphate molar ratios Sp/P of the order of 0.1 to 0.2), the complexes with poly(dA).poly(dT) and those with poly(dA-dT) displayed an important contribution of a permanent dipole moment to the orientation mechanism, as detected by the application of bipolar pulses in electric dichroism experiments. Just prior to precipitation (at Sp/P slightly larger than 0.3), these polynucleotides show electric dichroism and relaxation times characteristics corresponding to toroidal particles formation resulting from a bending of their chains. This implies asymmetric binding to phosphate sites on A-T containing polynucleotides. At low Sp/P ratios, spermine induced a stiffening of poly(dG-dC). No influence of spermine on the orientation mechanism of this polynucleotide was detected for Sp/P values ranging from zero to 0.35. The spermine-induced bending of A-T rich regions thus appears to be essential for DNA condensation into toroidal particles.

Energetic coupling between DNA bending and base pair opening

The pathway for base pair opening within a B-DNA duplex is investigated by theoretical molecular modeling. The results show that the disruption of a single base pair is energetically compatible with the deductions made from hydrogen exchange measurements. In addition, it is found that the opening process is greatly facilitated by DNA bending and that, conversely, once a base pair is disrupted, DNA can bend very easily. It appears that the energetic coupling between these two processes may play an important role in many biological reactions involving nucleic acid distortion.

Low-frequency collective motion in biomacromolecules and its biological functions

The occurrence of low-frequency motion in biomacromolecules, which had long been speculated upon with a great deal of skepticism, is now a clearly established phenomenon and has been convincingly demonstrated. The next stage in the process of its elucidation appropriately concerns the determination of its origin, the development of a feasible and effective model for the calculations involved and, more importantly, the extending of investigations on its biological roles in order to gain insights into the various interesting mechanisms underlying the dynamic processes occurring in biomacromolecules. Confronted with such a task, this review has been written with the aim of stimulating further, through a systematic and comprehensive description, developments to open up this exciting frontier of molecular biology, especially from the viewpoint of biological functions.

Principles of sequence-dependent flexure of DNA

The curvature of a bent rod may be defined in several different, but equivalent ways. The best way of describing the curvature of double-helical DNA is by an angle of turning per base step. Curvature comes mainly from the angle of roll between successive base-pairs, and this is defined as positive when the angle opens up on the minor groove side of the bases. DNA forms a plane curve if the roll angle values along the molecule alternate periodically between positive and negative, with a complete period equal to the helical repeat. It is known from studies of crystallized oligomers that the roll angles for particular dinucleotide steps have preferred values, or lie in preferred ranges of values. Therefore the formation of a plane curve will be easier with some base sequences of DNA than with others. We set up a computer algorithm for determining the ease with which DNA of given sequence will adopt a curved form. The algorithm has two different sets of constants: in model 1 the base step parameters come from an inspection of crystallized oligomers, and in model 2 data from a statistical survey of the incidence of dinucleotide steps in a large number of samples of chicken erythrocyte core DNA is incorporated. Both forms of the algorithm successfully locate the dyad of the nucleosome sequence (modulo 10) in a frog gene, and suggest strongly that sequence-dependent flexural properties of DNA play a part in the recognition of binding sites by nucleosome cores.

Origin of low-frequency motions in biological macromolecules. A view of recent progress in the quasi-continuity model

The recent progress in the quasi-continuity model and its applications in studying the low-frequency internal motions of biological macromolecules have been surveyed. Emphasis is placed on revealing the origin of this kind of internal collective motion, which involves many atoms and has significant biological functions. In light of such a line, the low-frequency motions in alpha-helix structure, beta-structure (including beta-sheet and beta-barrel), and DNA double-helix structure, the three most fundamental component elements in biological macromolecules, are discussed, and the corresponding physical pictures described. It turns out that the low-frequency motion in biological macromolecules originates from their two common intrinsic characteristics, i.e., they possess a series of weak bonds, such as hydrogen bonds and salt bridges, and a substantial mass distributed over the region containing those weak bonds.

Micrococcal nuclease as a DNA structural probe: its recognition sequences, their genomic distribution and correlation with DNA structure determinants

We have analyzed micrococcal nuclease (MNase) DNA cleavage patterns at the sequence level by examining 2.3 X 10(3) base-pairs of data derived from the Drosophila melanogaster 44D larval cuticle locus. Within this region, MNase preferentially cleaved 140 sites. Clusters of these sites appear to generate the preferential MNase eukaryotic DNA cleavage sites seen on agarose gels at roughly 100 to 300 base-pair intervals. These clusters of preferential cleavage sites rarely occur within gene coding regions. The analysis revealed that duplex DNA sequences preferentially cleaved by MNase are generally determined by a single strand sequence: d(A-T)n, where n greater than or equal to 1, flanked by a 5' dC or dG. Cleavage of the other strand is generally staggered 5' by several nucleotides and occurs even if such sequences are absent on that strand. An empirical predictive DNA cleavage model derived from a statistical analysis of the sequence level data was applied to seven eukaryotic gene loci of known sequence. The predicted patterns were in good general agreement with the previously observed eukaryotic gene/spacer cleavage pattern. Statistical analysis also revealed that sites of predicted preferential DNA cleavage occur less frequently in protein coding regions than for randomized sequences of the same length and nucleotide content. Comparison of the MNase cleavage patterns to the sequence-dependent pattern of binding energies between duplex DNA strands indicates that MNase preferentially cleaves sequences with low helix stability.

The twist, writhe, and linking number distributions in closed circular DNA

We discuss the predictions which follow from the assumption of statistically independent twist and writhe distributions of given variances in circular DNA with single-strand nicks. The nature of the topoisomer distribution produced upon covalent closure of the nicks is described, as well as the nature of the twist and writhe distributions in the fully-closed molecules. In particular, we show how the distributions depend on the magnitudes of the given variances, and how the relative magnitudes of the variances can be deduced from experiment. One additional consequence of the theory is the prediction of a necessary difference between the temperature coefficient of the twist in nicked versus fully-closed circular DNA. The ratio of the two twist coefficients turns out to depend only on the ratio of the twist and writhe variances in nicked DNA.

The trajectory of a stiff rod in a curved potential energy trough. An initial result for short nucleosomal rods

The equilibrium trajectory of the axis of a rod subject to an externally imposed curved potential energy trough tends to conform to the shape of the curved trough, but also tends to be straight because of elastic resistance to bending. The actual path of the axis is a balance between the two extremes. We consider a potential energy trough centered along a circular arc of radius R. For a rod of small length compared to R, we show that the axis at equilibrium forms an arc of a circle of radius greater than R. The value of the radius of the axial path depends on the relative values of the Hooke's Law bending constant for the rod and the depth and width of the trough. Motivation for the calculation is provided by nucleosomal DNA, which conforms to the surface of a roughly cylindrical histone core at physiological ionic strength, but is observed to unwind into a partially extended conformation at very low ionic strength. We suggest that the rigidity to bending of short DNA segments becomes sufficiently great at low ionic strength to overcome attractive interactions with the histone surface. Alternately, of course, if during the cell cycle mutually attractive forces between DNA and histone core are weakened at constant ionic strength, the same type of unfolding would be expected to occur as the strength of the DNA-histone contacts drops below the level required to overcome elastic resistance to bending of the DNA rod.

Packaged DNA. An elastic model

We review and deepen a theory of elastic bending of DNA on a persistence length scale. In a regime of extensive charge neutralization the axis of the double helix is elastically unstable when straight. Its stable bent conformation allows nucleation of DNA toruses and in principle could direct the supercoiled (solenoid) form of a polynucleosome. The Euler theory of elastic instability of macroscopic rods gives a partial description of the intrinsic ability of DNA to form locally stable bends. A different, quasi-Eulerian theory can be based on what is probably the dominant bending mechanism of DNA in solution-flexible kinking at the sites of open base pairs. This predictive theory is in quantitative agreement with the observed value (about 16 nm) for the minimum radius of torus holes. Stability of the inner torus ring is achieved when DNA phosphate groups are about 90% neutralized by trivalent cations, another prediction that is consistent with the observed formation of toruses in these conditions. The predicted stable radius of curvature of charge-neutralized DNA is also equal to the radial dimension of a maximally contracted polynucleosome supercoil as measured by neutron scattering (17 nm), but further experimental investigation of the geometrical disposition of the spacer DNA regions in the solenoid will be necessary to rule out the possibility of accidental agreement for this complex system. We stress again the experimental reality and probable importance of open base pairs in the equilibrium solution conformation of DNA.

Torsional and bending rigidity of the double helix from data on small DNA rings

We have calculated the variance of equilibrium distribution of a circular wormlike polymer chain over the writhing number, [Wr)2), as a function of the number of Kuhn statistical segments, n. For large n these data splice well with our earlier results obtained for a circular freely jointed polymer chain. Assuming that [delta Lk)2) = [delta Tw)2) we have compared our results with experimental data on the chain length dependence of the [delta Lk)2) value recently obtained by Horowitz and Wang for small DNA rings. This comparison has shown an excellent agreement between theory and experiment and yielded a reliable estimate of the torsional and bending rigidity parameters. Namely, the torsional rigidity constant is C = 3.0.10(-19) erg cm, and the bending rigidity as expressed in terms of the DNA persistence length is a = 500 A. The obtained value of C agrees well with earlier estimates by Shore and Baldwin as well as by Horowitz and Wang whereas the a value is in accord with the data of Hagerman. We have found the data of Shore and Baldwin on the chain length dependence of the [delta Lk)2) value to be entirely inconsistent with our theorectical results.

The flexibility of alternating dA-dT sequences

The flexibility of alternating poly (dA-dT) has been investigated by the technique of transient electric dichroism. Rotational relaxation times, which are very sensitive to changes in the end-to-end length of flexible polymers, are determined from the field free dichroism decay curves of four, well defined fragments of poly (dA-dT) ranging in size from 136 to 270 base pairs. Persistence lengths, calculated from the results of Hagerman and Zimm (Biopolymers (1981) 29, 1481-1502), are in the range 200-250 A. This makes alternating dA-dT sequences about twice as flexible as naturally occurring, "random" sequence DNA. Considering a bend around a nucleosome, for example, this difference in persistence length translates to an energy difference between poly (dA-dT) and random sequence DNA of 0.17 kT/base pair or 1 kcal per 10 base pair stretch. This energy difference is sufficiently large to suggest that dA-dT sequences could serve as markers in DNA packaging, for example, at sites where DNA must tightly bend to accommodate structures.

Low-frequency vibrations of DNA molecules

A model for calculating the low-frequency modes in DNA molecules is presented. The present model is associated with the 'breathing' of a DNA molecule as well as its complementary hydrogen bonds. The calculated results show excellent agreement with the observed low-frequency wavenumber (30 cm-1). Consequently, such an internal motion as reflected in the proposed model might be the origin of the observed low-frequency vibration in DNA molecules. This is helpful for investigating the relevant biological functions, which so far have been discussed by many scientists.