An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data

Curved DNA fragments are often found near functionally important sites such as promoters and origins of replication, and hence sequence-dependent DNA curvature prediction is of great utility in genomics and bioinformatics. In light of this, an assessment of three different dinucleotide step parameters (based on gel retardation as well as crystal structure data) is carried out. These parameters (BMHT, LB and CS) are evaluated quantitatively for their ability to predict correctly the experimental results of a large set of nucleic acid sequences containing A-tracts as well as GC-rich motifs. This set contained around 40 synthetic as well as natural sequences whose solution properties have been well characterized experimentally. All three models could account reasonably well for curvature in the various DNA sequences. The CS model, where dinucleotide parameters are calculated from crystal structure data, consistently shows slightly better correlation with experimental data. Our simple analysis also indicates that presently available trinucleotide parameters fail to predict curvature in some of the well-characterized sequences. The study shows that the dinucleotide parameters with some further refinement can be used to predict sequence-dependent curvature correctly in genomic sequences.     

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816–3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed ‘quasi-regularly’ throughout the genome, including both the coding and intergenic sequences. Second, there is a 10–12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in ~100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the ~100 bp long clusters of the phased A-tracts constitute the ‘structural code’ for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins.     

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.     

Attenuation of DNA-protein interactions associated with intrinsic, sequence-dependent DNA curvature.

Inherently curved DNA segments, associated with short runs of adenines, have been identified in many gene regulatory regions, yet their physiological significance remains unknown. The observations reported in this study indicate that intrinsically bent nucleic acid fragments are characterized by substantially attenuated affinities toward DNA-binding proteins involved in structural functions, such as H1 histone and protamine, as well as toward various DNA-modifying enzymes including ligases and exo- and endonucleases. Two mechanisms might be responsible for the altered binding properties. According to the first mechanism, the attenuated binding affinities and the bending represent two independent consequences of the unique structural parameters exhibited by A-tracts. Indeed, analysis of the degradation products obtained upon exposure of the curved sequences to various chemical nucleases points toward the narrowing of the DNA minor groove, a conformational modulation known to characterize A-tracts and to run along the axially-bent motifs, as a potential determinant of the observed binding attenuation. Alternatively, the conformational constraints which result from the stable bending might act to modulate the strength of DNA-protein interactions. Although the factor directly responsible for the altered binding affinities revealed by the bent sequences cannot as yet be conclusively resolved, it is proposed that a reiteration of this specific factor, being either an A-tract or a bend, in phase with the DNA helical repeat acts to amplify the modulation of the binding.(ABSTRACT TRUNCATED AT 250 WORDS)     

Condensation of oligonucleotides assembled into nicked and gapped duplexes: potential structures for oligonucleotide delivery

The condensation of nucleic acids into well-defined particles is an integral part of several approaches to artificial cellular delivery. Improvements in the efficiency of nucleic acid delivery in vivo are important for the development of DNA- and RNA-based therapeutics. Presently, most efforts to improve the condensation and delivery of nucleic acids have focused on the synthesis of novel condensing agents. However, short oligonucleotides are not as easy to condense into well-defined particles as gene-length DNA polymers and present particular challenges for discrete particle formation. We describe a novel strategy for improving the condensation and packaging of oligonucleotides that is based on the self-organization of half-sliding complementary oligonucleotides into long duplexes (ca. 2 kb). These non-covalent assemblies possess single-stranded nicks or single-stranded gaps at regular intervals along the duplex backbones. The condensation behavior of nicked- and gapped-DNA duplexes was investigated using several cationic condensing agents. Transmission electron microscopy and light-scattering studies reveal that these DNA duplexes condense much more readily than short duplex oligonucleotides (i.e. 21 bp), and more easily than a 3 kb plasmid DNA. The polymeric condensing agents, poly-l-lysine and polyethylenimine, form condensates with nicked- and gapped-DNA that are significantly smaller than condensates formed by the 3 kb plasmid DNA. These results demonstrate the ability for DNA structure and topology to alter nucleic acid condensation and suggest the potential for the use of this form of DNA in the design of vectors for oligonucleotide and gene delivery. The results presented here also provide new insights into the role of DNA flexibility in condensate formation.     

Uranyl photoprobing of nonbent A/T- and bent A-tracts. A difference of flexibility?

In this study, we have systematically compared the uranyl photocleavage of a range of bent A-tracts and nonbent TA-tracts as well as interrupted A-tracts. We demonstrate that uranyl photocleavage of A-tracts and TA-tracts is almost identical, indicating a very similar minor groove conformation. Furthermore, a 10 base pair A-tract is divided into two independent tracts by an intervening TA or GC step. Uranyl probing also clearly distinguishes the bent A4T4 and the nonbent T4A4 sequences as adopting different structures, and our interpretation of the data is consistent with a structure for the bent A4T4 sequence that resembles a continuous A-tract, whereas the nonbent T4A4 sequences are closer to two independent and opposite A-tracts that cancel each other in terms of macroscopic bending. Finally, we also note that even single TA and TAT steps are highly sensitive to uranyl photocleavage and propose that in addition to average minor groove width, uranyl also senses DNA helix flexibility/deformability. Thus, the structural difference of TA-tracts and A-tracts may to a large extent reflect a difference in flexibility, and DNA curvature may consequently require a rigid narrow minor groove conformation that creates distinct A-tract-B-DNA junctions as the predominant cause of the bending.     

Monovalent cation binding by curved DNA molecules containing variable numbers of a-tracts

Monovalent cation binding by DNA A-tracts, runs of four or more contiguous adenine or thymine residues, has been determined for two curved ∼200 basepair (bp) restriction fragments, one taken from the M13 origin of replication and the other from the VP1 gene of SV40. These two fragments have previously been shown to contain stable, centrally located bends of 44° and 46°, respectively, located within ∼60 bp “curvature modules” containing four or five irregularly spaced A-tracts. Transient electric birefringence measurements of these two fragments, sequence variants containing reduced numbers of A-tracts in the SV40 curvature module or changes in the residues flanking the A-tracts in the M13 curvature module, have been combined with the free solution electrophoretic mobilities of the same fragments using known equations to estimate the effective charge of each fragment. The effective charge is reduced, on average, by one-third charge for each A-tract in the curvature module, suggesting that each A-tract binds a monovalent cation approximately one-third of the time. Monovalent cation binding to two or more A-tracts is required to observe significant curvature of the DNA helix axis.     

Secondary conformational polymorphism of nucleic acids as a possible functional link between cellular parameters and DNA packaging processes

Circular dichroism and electron microscopy studies of various in vitro DNA packaging systems indicate that all the factors which induce and modulate the secondary conformation of DNA molecules are capable of eliciting nucleic acids condensation processes into tight, highly ordered tertiary structures as well as altering the extent of order and compactness within the resulting species. Specifically, such factors include the ionic strength, the presence of particular dehydrating agents and polyamines, as well as the pH values. It is proposed that slight alterations of these parameters induce the formation of short non-B-DNA segments that propagate as a perturbation along the B-DNA double helix. The structural fluctuations of the dsDNA molecules that result from the conformational discontinuities formed at the junction sites between the B motif and the conformationally altered segments alter the elastic response of the nucleic acids and facilitate cooperative condensation processes. Moreover, the type and frequency of the structurally modified clusters interspersed within the B conformation and determined by the environmental parameters are shown to provide a means for continuous regulation of the extent and mode of DNA packaging. The ionic strength and hydrophobic environment in the close vicinity of the DNA molecules are controlled and modulated in vivo by DNA-binding proteins such as histones and protamines; similarly, pH values and polyamine concentrations are constantly regulated in living systems. It is suggested, therefore, that the secondary structural polymorphism which characterizes the DNA molecules might display a regulatory role by acting as a functional link between cellular parameters and the extent, mode, and timing of nucleic acid packaging processes.     

A-tract polarity dominate the curvature in flanking sequences

It is well-known, but little understood, that the nucleotide sequences between phased A(4-6)-tracts (at 10-11 bp intervals) have only a slight effect on overall curvature. To explore this phenomenon, we have examined the gel-migration properties of sequences containing both A-tracts as well as G-tracts (i.e., sequences of the form G(n)C(m) or C(n)G(m), n + m > 4) in various relative positioning. We show that the composite bend of these sequences depends on their relative arrangement. When G-tracts are placed between two A-tracts, such that both tracts are repeated in phase to themselves (e.g., G(5)A(6)G(5)A(5)), or adjacent to the 3'-side of A-tracts (e.g., A(6)G(5)N(10)), they have minimal influence on the extent of bending of the composite sequence. When G-tracts are placed one helical repeat away from A-tracts (e.g., G(5)N(5)A(6)N(6)), or are adjacent only to the 5'-side of A-tracts (e.g., G(5)A(6)N(10)) their influence on the composite bend is larger. The differential behavior of AG- versus GA-tracts means that A-tracts influence their flanking sequences in a polar manner. Whereas they suppress, or make constant, the intrinsic bending characteristics of any sequence placed immediately 3' to them (and hence by definition any sequence placed between two phased A-tracts), sequences adjoining them on their 5'-side are free to modulate the overall curvature. We interpret these results as evidence for the dominant nature of the unique and nonuniform structure adopted by tracts of four adenines or more. The effects of A-tracts extend at least five base pairs into the adjoining 3' region. This is further evidence for the complexity of DNA structure and the inadequacy of simple nearest-neighbor models to explain all its manifestations.     

Curvature of mouse satellite DNA and condensation of heterochromatin

Cloned, sequenced mouse satellite DNA exhibits properties characteristic of molecules that possess a stable curvature. Circularly permuted fragments containing the region predicted to bend were used to map the curvature relative to DNA sequence. The altered mobility of these fragments in polyacrylamide gels is reversed when gels are run in the presence of distamycin A, a drug that binds preferentially to AT-rich DNA. Treatment of living mouse cells with this drug dramatically reduces the condensation of centromeric heterochromatin, the exclusive location of satellite sequences. In situ hybridization of satellite probes to extended chromosomes at the electron microscope level shows that satellite does not comprise a single block but is distributed throughout the centromere region. Based on these experiments, we hypothesize that the structure of mouse satellite DNA is an important feature of centromeric heterochromatin condensation.     

DNA chip replication for a personalized DNA chip

We report the replication technology of DNA chip using by sequence specific localization of nucleic acids via hybridization and electric transfer of the nucleic acids onto a new substrate without losing their array information. The denatured DNA fragments are first spotted and UV-cross-linked on a nylon membrane. The membrane is then immersed and hybridized in a DNA mixture solution that contains all complementary sequences of the nucleic acids to be hybridized with the DNA fragments on the membrane. The hybridized DNA fragments are transferred to another membrane at the denatured condition. After separating two membranes, the transferred membrane contains a complementary array of DNA fragments. This method can be used for the replication of the same copy of DNA chip repeatedly and moreover could be applied for a personalized DNA chip fabrication, where specific information of each spot of DNA chip is originated from the genetic information of a personal sample.     

Role of DNA-DNA interactions on the structure and thermodynamics of bacteriophages Lambda and P4

Electrostatic interactions play an important role in both packaging of DNA inside bacteriophages and its release into bacterial cells. While at physiological conditions DNA strands repel each other, the presence of polyvalent cations such as spermine and spermidine in solutions leads to the formation of DNA condensates. In this study, we discuss packaging of DNA into bacteriophages P4 and Lambda under repulsive and attractive conditions using a coarse-grained model of DNA and capsids. Packaging under repulsive conditions leads to the appearance of the coaxial spooling conformations; DNA occupies all available space inside the capsid. Under the attractive potential both packed systems reveal toroidal conformations, leaving the central part of the capsids empty. We also present a detailed thermodynamic analysis of packaging and show that the forces required to pack the genomes in the presence of polyamines are significantly lower than those observed under repulsive conditions. The analysis reveals that in both the repulsive and attractive regimes the entropic penalty of DNA confinement has a significant non-negligible contribution into the total energy of packaging. Additionally we report the results of simulations of DNA condensation inside partially packed Lambda. We found that at low densities DNA behaves as free unconfined polymer and condenses into the toroidal structures; at higher densities rearrangement of the genome into toroids becomes hindered, and condensation results in the formation of non-equilibrium structures. In all cases packaging in a specific conformation occurs as a result of interplay between bending stresses experienced by the confined polymer and interactions between the strands.     

Rotational dynamics of curved DNA fragments studied by fluorescence polarization anisotropy

The rotational dynamics of short DNA fragments with or without intrinsic curvature were studied using time-resolved phase fluorimetry of intercalated ethidium with detection of the anisotropy. Parameters determined were the spinning diffusion coefficient of the DNA fragments about the long axis and the zero-time ethidium fluorescence anisotropy. We find a significant decrease in the spinning diffusion coefficient for all curved fragments compared to the straight controls. This decrease is likewise evident in rotational diffusion coefficients computed from DNA structures obtained by a curvature prediction program for these sequences. Using a hinged-cylinder model, we can identify the change in rotational diffusion coefficient with a permanent bend of 13-16 degrees per helix turn for the sequences studied. Moreover, for some of the curved fragments an increased flexibility has to be assumed in addition to the permanent bend in order to explain the data.     

DNA packaging induced by micellar aggregates: a novel in vitro DNA condensation system.

Evidence for a conceptually novel DNA packaging process is presented. X-ray scattering, electron microscopy, and circular dichroism measurements indicate that in the presence of positively charged micellar aggregates and flexible anionic polymers, such as negatively charged polypeptides or single-stranded RNA species, a complex is formed in which DNA molecules are partially embedded within a micellar scaffold and partially condensed into highly packed chiral structures. Based on studies of micelle-DNA and micelle-flexible anionic polymer systems, as well as on the known effects of a high charge density upon the micellar organization, a DNA packaging model is proposed. According to this model, the DNA induces the elongation of the micelles into rodlike aggregates, forming a closely packed matrix in which the DNA molecules are immobilized. In contrast, the flexible anionic polymers stabilize clusters of spherical micelles which are proposed to effect a capping of the rodlike micelles, thus arresting their elongation and creating surfactant-free segments of the DNA that are able to converge and collapse. Thus, unlike other in vitro DNA packaging systems, in which condensation follows encounters between charge-neutralized DNA molecules, a prepackaging phase where the DNA is immobilized within a matrix is proposed in this case. Cellular and nuclear membranes have been implicated in DNA packaging processes in vivo, and negatively charged polyelectrolytes were shown to be involved in the processes. These observations, combined with the basic tenets of the DNA condensation system described here, allow for the progression to the study of more elaborate model systems and thus might lead to insights into the nature and roles of the intricate in vivo DNA-membrane complexes.     

Increased temperature and 2-methyl-2,4-pentanediol change the DNA structure of both curved and uncurved adenine/thymine-rich sequences

DNA curvature is affected by elevated temperature and dehydrating agents such as 2-methyl-2,4-pentanediol (MPD) (used in crystallization). This effect of MPD has been ascribed to a specific distortion of the structure of adenine tracts (A-tracts), probably through a deformation of the characteristic narrow minor groove. Uranyl photoprobing indicates that a narrowed minor groove is present in all A/T regions containing four or more A/T base pairs. Consequently, this technique may be employed to study conformational changes in other A/T-rich sequences than pure A-tracts. In this study we use uranyl photoprobing to demonstrate that the effect of elevated temperature and MPD is analogous on both "normal" and curve-inducing A/T-rich sequences. The results therefore indicate that under these conditions the minor groove is widened in all A/T sequences and not only in pure A-tracts as previously suggested. Thus, the rather subtle structural difference of AT regions and A-tracts in nonbent DNA versus A-tracts in bent DNA may be quantitative rather than qualitative; i.e., the structure is more persistent and/or rigid in bent DNA.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

The influence of the thymine C5 methyl group on spontaneous base pair breathing in DNA

Sequences of four or more AT base pairs without a 5'-TA-3' step, so-called A-tracts, influence the global properties of DNA by causing curvature of the helix axis if phased with the helical repeat and also influence nucleosome packaging. Hence it is interesting to understand this phenomenon on the molecular level, and numerous studies have been devoted to investigations of dynamical and structural features of A-tract DNA. It was early observed that anomalously slow base pair-opening kinetics were a striking physical property unique to DNA A-tracts (Leroy, J. L., Charretier, E., Kochoyan, M., and Gueron, M. (1988) Biochemistry 27, 8894-8898). Furthermore, a strong correlation between DNA curvature and anomalously slow base pair-opening dynamics was found. In the present work it is shown, using imino proton exchange measurements by NMR spectroscopy that the main contribution to the dampening of the base pair-opening fluctuations in A-tracts comes from the C5 methylation of the thymine base. Because the methyl group has been shown to have a very limited effect on the DNA curvature as well as the structure of the DNA helix, the thymine C5 methyl group stabilizes the helix directly. Empirical potential energy calculations show that methylation of the tract improves the stacking energy of a base pair with its neighbors in the tract by 3-4 kcal/mol.     

A refined prediction method for gel retardation of DNA oligonucleotides from dinucleotide step parameters: reconciliation of DNA bending models with crystal structure data

The development and assessment of a prediction method for gel retardation and sequence dependent curvature of DNA based on dinulcleotide step parameters are described. The method is formulated using the Babcock-Olson equations for base pair step geometry (1) and employs Monte Carlo simulated annealing for parameter optimization against experimental data. The refined base pair step parameters define a stuctural construct which, when the width of observed parameter distributions is taken into account, is consistent with the results of DNA oligonucleotide crystal structures. The predictive power of the method is demonstrated and tested via comparisons with DNA bending data on sets of sequences not included in the training set, including A-tracts with and without periodic helix phasing, phased A4T4 and T4A4 motifs, a sequence with a phased GGGCCC motif, some "unconventional" helix phasing sequences, and three short fragments of kinetoplast DNA from Crithidia fasiculata that exhibit significantly different behavior on non-denaturing polyacrylamide gels. The nature of the structural construct produced by the methodology is discussed with respect to static and dynamic models of structure and representations of bending and bendability. An independent theoretical account of sequence dependent chemical footprinting results is provided. Detailed analysis of sequences with A-tract induced axis bending forms the basis for a critical discussion of the applicability of wedge models,junction models and non A-tract, general sequence models for understanding the origin of DNA curvature at the molecular level.     

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.