Thermodynamics and premelting conformational changes of phased (dA)5 tracts

From straight to curved: A historical perspective of DNA shape

DNA is the macromolecule responsible for storing the genetic information of a cell and it has intrinsic properties such as deformability, stability and curvature. DNA Curvature plays an important role in gene transcription and, consequently, in the subsequent production of proteins, a fundamental process of cells. With recent advances in bioinformatics and theoretical biology, it became possible to analyze and understand the involvement of DNA Curvature as a discriminatory characteristic of gene-promoting regions. These regions act as sites where RNAp (ribonucleic acid-polymerase) binds to initiate transcription. This review aims to describe the formation of Curvature, as well as highlight its importance in predicting promoters. Furthermore, this article provides the potential of DNA Curvature as a distinguishing feature for promoter prediction tools, as well as outlining the calculation procedures that have been described by other researchers. This work may support further studies directed towards the enhancement of promoter prediction software.

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

Force-induced melting of DNA hairpin: Unfolding pathways and phase diagrams

Using the exact enumeration technique, we have studied the force-induced melting of a DNA hairpin on the face centered cubic lattice for two different sequences which differ in terms of loop closing base pairs. The melting profiles obtained from the exact enumeration technique is consistent with the Gaussian network model and Langevin dynamics simulations. Probability distribution analysis based on the exact density of states revealed the microscopic details of the opening of the hairpin. We showed the existence of intermediate states near the melting temperature. We further showed that different ensembles used to model single-molecule force spectroscopy setups may give different force-temperature diagrams. We delineate the possible reasons for the observed discrepancies.

Nonlinear physics opens a new paradigm for accurate transcription start site prediction

There is evidence that DNA breathing (spontaneous opening of the DNA strands) plays a relevant role in the interactions of DNA with other molecules, and in particular in the transcription process. Therefore, having physical models that can predict these openings is of interest. However, this source of information has not been used before either in transcription start sites (TSSs) or promoter prediction. In this article, one such model is used as an additional information source that, when used by a machine learning (ML) model, improves the results of current methods for the prediction of TSSs. In addition, we provide evidence on the validity of the physical model, as it is able by itself to predict TSSs with high accuracy. This opens an exciting avenue of research at the intersection of statistical mechanics and ML, where ML models in bioinformatics can be improved using physical models of DNA as feature extractors.

Origin of heat capacity increment in DNA folding: The hydration effect

BACKGROUND

Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆c_{p, F}, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion.

METHODS

Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆c_p^hyd.

RESULTS

We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol^-1 K^-1 Å^-2 and can be used to estimate the ∆c_{p, F} contribution for DNA folding.

CONCLUSIONS

We show that dehydration is one of the major sources contributing to the observed ∆c_{p, F} increment in DNA folding. Other possible sources contributing to the overall ∆c_{p, F} should be significant but appear to compensate each other to high degree. The calculated ∆c_p^hyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆c_{p, F} values. By contrast, empirical parametrizations developed for proteins result in poor ∆c_p^hyd predictions and should not be applied to DNA folding.

GENERAL SIGNIFICANCE

Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.

Sufficient minimal model for DNA denaturation: Integration of harmonic scalar elasticity and bond energies

We study DNA denaturation by integrating elasticity - as described by the Gaussian network model - with bond binding energies, distinguishing between different base pairs and stacking energies. We use exact calculation, within the model, of the Helmholtz free-energy of any partial denaturation state, which implies that the entropy of all formed "bubbles" ("loops") is accounted for. Considering base pair bond removal single events, the bond designated for opening is chosen by minimizing the free-energy difference for the process, over all remaining base pair bonds. Despite of its great simplicity, for several known DNA sequences our results are in accord with available theoretical and experimental studies. Moreover, we report free-energy profiles along the denaturation pathway, which allow to detect stable or meta-stable partial denaturation states, composed of bubble, as local free-energy minima separated by barriers. Our approach allows to study very long DNA strands with commonly available computational power, as we demonstrate for a few random sequences in the range 200-800 base-pairs. For the latter, we also elucidate the self-averaging property of the system. Implications for the well known breathing dynamics of DNA are elucidated.

Flanking A·T basepairs destabilize the B(∗) conformation of DNA A-tracts

Capillary electrophoresis has been used to characterize the interaction of monovalent cations with 26-basepair DNA oligomers containing A-tracts embedded in flanking sequences with different basepair compositions. A 26-basepair random-sequence oligomer was used as the reference; lithium and tetrabutylammonium (TBA(+)) ions were used as the probe ions. The free solution mobilities of the A-tract and random-sequence oligomers were identical in solutions containing <∼ 100 mM cation. At higher cation concentrations, the A-tract oligomers migrated faster than the reference oligomer in TBA(+) and slower than the reference in Li(+). Hence, cations of different sizes can interact very differently with DNA A-tracts. The increased mobilities observed in TBA(+) suggest that the large hydrophobic TBA(+) ions are preferentially excluded from the vicinity of the A-tract minor groove, increasing the effective net charge of the A-tract oligomers and increasing the mobility. By contrast, Li(+) ions decrease the mobility of A-tract oligomers because of the preferential localization of Li(+) ions in the narrow A-tract minor groove. Embedding the A-tracts in AT-rich flanking sequences markedly alters preferential interactions of monovalent cations with the B(∗) conformation. Hence, A-tracts embedded in genomic DNA may or may not interact preferentially with monovalent cations, depending on the relative number of A · T basepairs in the flanking sequences.

Structural waters in the minor and major grooves of DNA--a major factor governing structural adjustments of the A-T mini-helix

The role of microhydration in structural adjustments of the AT-tract in B-DNA was studied at the B97-D/def2-SV(P) level. The (dA:dT)5 complexes with 10 water molecules in minor and 15 water molecules in major grooves were studied. The obtained network of hydrogen bonds revealed the dependence between the groove width and the types of water patterns. In the minor groove, the following patterns were observed: interstrand one-water bridges similar to that of the Dickerson "water spine" and interstrand two-water bridges. The network of structural waters in the major groove is more diverse than that in the minor groove, which agrees with crystallographic data. As the major groove is wider, it is enriched by water molecules forming two- and three-water bridges. Results suggest the nucleobase-water interactions in both grooves prevent AT-tract twisting and its "collapse" along the minor groove. Whereby, a helix structure with narrow minor and wide major grooves is formed. The structural waters affect the polynucleotide conformation so that it becomes similar to poly(dA)·poly(dT) in fibers and acquires features of the A-tracts in DNA in solution. We suggest that formation of specific water patterns in both grooves is the factor responsible for stabilization of A-tracts with a narrowed minor groove, leading in turn to their strong intrinsic bending in DNA.

Guanine radical chemistry reveals the effect of thermal fluctuations in gene promoter regions

DNA is not the static entity that structural pictures suggest. It has been longly known that it 'breathes' and fluctuates by local opening of the bases. Here we show that the effect of structural fluctuations, exhibited by AT-rich low stability regions present in some common transcription initiation regions, influences the properties of DNA in a distant range of at least 10 bp. This observation is confirmed by experiments on genuine gene promoter regions of DNA. The spatial correlations revealed by these experiments throw a new light on the physics of DNA and could have biological implications, for instance by contributing to the cooperative effects needed to assemble the molecular machinery that forms the transcription complex.

Effect of magnesium ions and temperature on the sequence-dependent curvature of DNA restriction fragments

Transient electric birefringence has been used to quantify the curvature of two DNA restriction fragments, a 199-base-pair fragment taken from the origin of replication of the M13 bacteriophage and a 207-base-pair fragment taken from the VP1 gene in the SV40 minichromosome. Stable curvature in the SV40 and M13 restriction fragments is due to a series of closely spaced A tracts, runs of 4-6 contiguous adenine residues located within 40 or 60 base pair 'curvature modules' near the center of each fragment. The M13 and SV40 restriction fragments exhibit bends of ∼ 45° in solutions containing monovalent cations and ∼ 60° in solutions containing Mg(2 +) ions. The curvature is not localized at a single site but is distributed over the various A tracts in the curvature modules. Thermal denaturation studies indicate that the curvature in the M13 and SV40 restriction fragments remains constant up to 30 °C in solutions containing monovalent cations, and up to 40 °C in solutions containing Mg(2 +) ions, before beginning to decrease slowly with increasing temperature. Hence, stable curvature in these DNA restriction fragments exists at the biologically important temperature of 37 °C.

A measure of bending in nucleic acids structures applied to A-tract DNA

A method is proposed to measure global bending in DNA and RNA structures. It relies on a properly defined averaging of base-fixed coordinate frames, computes mean frames of suitably chosen groups of bases and uses these mean frames to evaluate bending. The method is applied to DNA A-tracts, known to induce considerable bend to the double helix. We performed atomistic molecular dynamics simulations of sequences containing the A(4)T(4) and T(4)A(4) tracts, in a single copy and in two copies phased with the helical repeat. Various temperature and salt conditions were investigated. Our simulations indicate bending by roughly 10 degrees per A(4)T(4) tract into the minor groove, and an essentially straight structure containing T(4)A(4), in agreement with electrophoretic mobility data. In contrast, we show that the published NMR structures of analogous sequences containing A(4)T(4) and T(4)A(4) tracts are significantly bent into the minor groove for both sequences, although bending is less pronounced for the T(4)A(4) containing sequence. The bending magnitudes obtained by frame averaging are confirmed by the analysis of superhelices composed of repeated tract monomers.

Preparation, resonance assignment, and preliminary dynamics characterization of residue specific 13C/15N-labeled elongated DNA for the study of sequence-directed dynamics by NMR

DNA is a highly flexible molecule that undergoes functionally important structural transitions in response to external cellular stimuli. Atomic level spin relaxation NMR studies of DNA dynamics have been limited to short duplexes in which sensitivity to biologically relevant fluctuations occurring at nanosecond timescales is often inadequate. Here, we introduce a method for preparing residue-specific (13)C/(15)N-labeled elongated DNA along with a strategy for establishing resonance assignments and apply the approach to probe fast inter-helical bending motions induced by an adenine tract. Preliminary results suggest the presence of elevated A-tract independent end-fraying internal motions occurring at nanosecond timescales, which evade detection in short DNA constructs and that penetrate deep (7 bp) within the DNA helix and gradually fade away towards the helix interior.

The unique structure of A-tracts and intrinsic DNA bending

Short runs of adenines are a ubiquitous DNA element in regulatory regions of many organisms. When runs of 4–6 adenine base pairs (‘A-tracts’) are repeated with the helical periodicity, they give rise to global curvature of the DNA double helix, which can be macroscopically characterized by anomalously slow migration on polyacrylamide gels. The molecular structure of these DNA tracts is unusual and distinct from that of canonical B-DNA. We review here our current knowledge about the molecular details of A-tract structure and its interaction with sequences flanking them of either side and with the environment. Various molecular models were proposed to describe A-tract structure and how it causes global deflection of the DNA helical axis. We review old and recent findings that enable us to amalgamate the various findings to one model that conforms to the experimental data. Sequences containing phased repeats of A-tracts have from the very beginning been synonymous with global intrinsic DNA bending. In this review, we show that very often it is the unique structure of A-tracts that is at the basis of their widespread occurrence in regulatory regions of many organisms. Thus, the biological importance of A-tracts may often be residing in their distinct structure rather than in the global curvature that they induce on sequences containing them.

Properties of the Nucleic-acid Bases in Free and Watson-Crick Hydrogen-bonded States: Computational Insights into the Sequence-dependent Features of Double-helical DNA

The nucleic-acid bases carry structural and energetic signatures that contribute to the unique features of genetic sequences. Here we review the connection between the chemical structure of the constituent nucleotides and the polymeric properties of DNA. The sequence-dependent accumulation of charge on the major- and minor-groove edges of the Watson-Crick base pairs, obtained from ab initio calculations, presents unique motifs for direct sequence recognition. The optimization of base interactions generates a propellering of base-pair planes of the same handedness as that found in high-resolution double-helical structures. The optimized base pairs also deform along conformational pathways, i.e., normal modes, of the same type induced by the binding of proteins. Empirical energy computations that incorporate the properties of the base pairs account satisfactorily for general features of the next level of double-helical structure, but miss key sequence-dependent differences in dimeric structure and deformability. The latter discrepancies appear to reflect factors other than intrinsic base-pair structure.

Poly(dA:dT) tracts: major determinants of nucleosome organization

Homopolymeric stretches of deoxyadenosine nucleotides (A's) on one strand of double-stranded DNA, referred to as poly(dA:dT) tracts or A-tracts, are overabundant in eukaryotic genomes. They have unusual structural, dynamic, and mechanical properties, and may resist sharp bending. Such unusual material properties, together with their overabundance in eukaryotes, raised the possibility that poly(dA:dT) tracts might function in eukaryotes to influence the organization of nucleosomes at many genomic regions. Recent genome-wide studies strongly confirm these ideas and suggest that these tracts play major roles in chromatin organization and genome function. Here we review what is known about poly(dA:dT) tracts and how they work.

Experimental and theoretical studies of sequence effects on the fluctuation and melting of short DNA molecules

Understanding the melting of short DNA sequences probes DNA at the scale of the genetic code and raises questions which are very different from those posed by very long sequences, which have been extensively studied. We investigate this problem by combining experiments and theory. A new experimental method allows us to make a mapping of the opening of the guanines along the sequence as a function of temperature. The results indicate that non-local effects may be important in DNA because an AT-rich region is able to influence the opening of a base pair which is about 10 base pairs away. An earlier mesoscopic model of DNA is modified to correctly describe the timescales associated with the opening of individual base pairs well below melting, and to properly take into account the sequence. Using this model to analyze some characteristic sequences for which detailed experimental data on the melting is available (Montrichok et al 2003 Europhys. Lett. 62 452), we show that we have to introduce non-local effects of AT-rich regions to get acceptable results. This brings a second indication that the influence of these highly fluctuating regions of DNA on their neighborhood can extend to some distance.

Unique Physical Signature of DNA Curvature and Its Implications for Structure and Dynamics

A particularly sensitive birefringence technique is used to analyze a curved DNA fragment with 118 bp and a standard DNA with 119 bp. At salt concentrations from 0.5 to 10 mM, both fragments show the usual negative stationary birefringence and monotonic transients - differences are relatively small. At 100 mM salt the curved DNA shows a positive stationary birefringence and non-monotonic transients with processes having amplitudes of opposite sign, whereas signals of the standard DNA remain as usual. Transients induced by reversal of the field vector indicate the existence of a permanent dipole for the curved DNA. 2-MHz-ac pulses induce a negative stationary birefringence in both DNAs. These results are consistent with calculations on models for curved DNA predicting a quasi-permanent dipole and a positive dichroism/birefringence. The quasi-permanent dipole results from the loss of symmetry in the charge distribution of the curved polyelectrolyte. The appearance of the unique signature of curvature at high salt is mainly due to a strong decrease of the polarizability by about 2 orders of magnitude. The special mode of orientation resulting from the quasi-permanent dipole is expected to contribute to the gel migration anomaly. The time constants of birefringence decay for the curved fragment are shorter than those of the 119 bp fragment by a factor of approximately 1.10 at 0.6 mM salt, whereas this factor is approximately 1.20 at 100 mM Na+. If both fragments were normal DNA with 3.4 A rise per base pair, the factor would be approximately 1.02. At high salt and high electric field strengths the factor increases up to 1.37. The implications for the bending dynamics and the potential to distinguish static from dynamic persistence by field reversal experiments are discussed. The dependence of the curvature on the salt concentration indicated by the time constants is consistent with a clear decrease of the electrophoretic anomaly at decreasing salt concentration.

Studies of DNA dumbbells VIII. Melting analysis of DNA dumbbells with dinucleotide repeat stem sequences

Melting curves and circular dichroism spectra were measured for a number of DNA dumbbell and linear molecules containing dinucleotide repeat sequences of different lengths. To study effects of different sequences on the melting and spectroscopic properties, six DNA dumbbells whose stems contain the central sequences (AA)(10), (AC)(10), (AG)(10), (AT)(10), (GC)(10), and (GG)(10) were prepared. These represent the minimal set of 10 possible dinucleotide repeats. To study effects of dinucleotide repeat length, dumbbells with the central sequences (AG)(n), n = 5 and 20, were prepared. Control molecules, dumbbells with a random central sequence, (RN)(n), n = 5, 10, and 20, were also prepared. The central sequence of each dumbbell was flanked on both sides by the same 12 base pairs and T(4) end-loops. Melting curves were measured by optical absorbance and differential scanning calorimetry in solvents containing 25, 55, 85, and 115 mM Na(+). CD spectra were collected from 20 to 45 degrees C and [Na(+)] from 25 to 115 mM. The spectral database did not reveal any apparent temperature dependence in the pretransition region. Analysis of the melting thermodynamics evaluated as a function of Na(+) provided a means for quantitatively estimating the counterion release with melting for the different sequences. Results show a very definite sequence dependence, indicating the salt-dependent properties of duplex DNA are also sequence dependent. Linear DNA molecules containing the (AG)(n) and (RN)(n), sequences, n = 5, 10, 20, and 30, were also prepared and studied. The linear DNA molecules had the exact sequences of the dumbbell stems. That is, the central repeat sequence in each linear duplex was flanked on both sides by the same 12-bp sequence. Melting and CD studies were also performed on the linear DNA molecules. Comparison of results obtained for the same sequences in dumbbell and linear molecular environments reveals several interesting features of the interplay between sequence-dependent structural variability, sequence length, and the unconstrained (linear) or constrained (dumbbell) molecular environments.

Examination of the premelting transition of DNA A-tracts using a fluorescent adenosine analogue

The fluorescent adenosine analogue 4-amino-8-(2-deoxy-beta-d-ribofuranosyl)-5'-O-dimethoxytrityl-6-methyl-7(8H)-pteridone (6MAP) has been used to perform residue specific analyses of DNA A-tracts during the premelting transition. DNA A-tracts, which exhibit sequence-induced curvature, adopt a B-DNA conformation as a function of increasing temperature. Fluorescence melting curves indicate that 6MAP is a more sensitive reporter of the premelting transition than UV absorption spectroscopy. Further, residue specific fluorescence analyses of A-tract and control duplexes reveal that some of the conformational changes associated with the premelting transition occur within A-tract regions. Analyses of the energetics of the premelting transition indicate that ApA steps make a larger enthalpic contribution to the premelting transition than ApT steps. To explore the effect of cations on the premelting transition, fluorescence melts were performed in the presence of NH(4)(+), Mg(2+), and low (0.05 M) and high (0.5 M) concentrations of Na(+). These studies show that the fluorescence intensity changes associated with the premelting transition are sensitive to cation type and concentration and are larger and more pronounced in the presence of 0.5 M Na(+), NH(4)(+), and Mg(2+). Incorporation of 6MAP into longer duplexes containing phased A-tracts shows that the local environment of adenosines in phased A-tracts is similar to that of individual A-tracts. Fluorescence quenching results indicate that ApA and ApT steps within A-tracts are less solvent exposed than their counterparts in control sequence isomers, possibly because of the narrowed minor groove of A-tract sequences.

Transcriptional organization of the temperature-sensitive transfer system from the IncHI1 plasmid R27

One of the characteristic features of IncHI1 plasmids is a thermosensitive process of conjugation, which is optimal between 22 °C and 30 °C but inhibited at 37 °C. R27, the prototypical IncHI1 plasmid, contains transfer genes clustered in two regions of the plasmid, Tra1 and Tra2. In the present study, transcriptional analyses of thetragenes were undertaken at both 30 °C and 37 °C. Screening of 38tragenes showed thattragenes are transcriptionally linked in six operons, three in each Tra region. RT-PCR analysis showed that gene expression was reduced at 37 °C relative to that observed at 30 °C. The transcription start sites of the six transcripts were identified, promoters and upstream regions were cloned, and transcription was tested at both temperatures. In cells grown at 37 °C, in the presence of R27, the promoters were inhibited, except for promoters of the H operon and AN operon. Conditions that influenced DNA topology, such as osmolarity, anaerobiosis, quorum sensing and acidity, showed no significant influence on transfer frequency. These results should facilitate future understanding of the basis of temperature-sensitive transfer in this large conjugative plasmid.

Structure and axis curvature in two dA6 x dT6 DNA oligonucleotides: comparison of molecular dynamics simulations with results from crystallography and NMR spectroscopy

Molecular dynamics (MD) simulations have been performed on the A6 containing DNA dodecamers d(GGCAAAAAACGG) solved by NMR and d(CGCAAAAAAGCG) solved by crystallography. The experimental structures differ in the direction of axis bending and in other small but important aspects relevant to the DNA curvature problem. Five nanosecond MD simulations of each sequence have been performed, beginning with both the NMR and crystal forms as well as canonical B-form DNA. The results show that all simulations converge to a common form in close proximity to the observed NMR structure, indicating that the structure obtained in the crystal is likely a strained form due to packing effects. A-tracts in the MD model are essentially straight. The origin of axis curvature is found at pyrimidine-purine steps in the flanking sequences.

Centromere parC of plasmid R1 is curved

The centromere sequence parC of Escherichia coli low-copy-number plasmid R1 consists of two sets of 11 bp iterated sequences. Here we analysed the intrinsic sequence-directed curvature of parC by its migration anomaly in polyacrylamide gels. The 159 bp long parC is strongly curved with anomaly values (k-factors) close to 2. The properties of the parC curvature agree with those of other curved DNA sequences. parC contains two regions of 5-fold repeated iterons separated by 39 bp. We modified 4 bp within this intermediate sequence so that we could analyse the two 5-fold repeated regions independently. The analysis shows that the two repeat regions are not independently curved parts of parC but that the overall curvature is a property of the whole fragment. Since the centromere sequence of an E.coli plasmid as well as eukaryotic centromere sequences show DNA curvature, we speculate that curvature might be a general property of centromeres.

Topological measurement of an A-tract bend angle: effect of magnesium

Sequences of four to six adenine residues, termed A-tracts, have been shown to produce curvature in the DNA double helix. A-tracts have been used extensively as reference standards to quantify bending induced by other sequences as well as by DNA binding proteins when they bind to their sites. However, the ability of an A-tract to serve as such a standard is hampered by the wide variation of values reported for the amount of bend conferred by an A-tract. One experimental condition that differs in these studies is the presence of divalent cation. To evaluate this effect, a new application of a topological method, termed rotational variant analysis, is used here to measure for the first time the effect of the presence of magnesium ion on the bend angle conferred by an A-tract. This method, which has the unique ability to measure a bend angle in the presence or absence of magnesium ion, demonstrates that magnesium ion markedly increases the bend angle. For example, when measured in a commonly used gel electrophoretic buffer, the bend angle conferred by a tract of six adenine residues increases from about 7 degrees in the absence of magnesium ion to 19 degrees in the presence of 3.9 mM magnesium ion. This quantitative demonstration of substantial magnesium ion dependence has several important implications. First, it explains discrepancies among bend values reported in various previous studies, particularly those employing gel electrophoretic versus other solution methods. In addition, these findings necessitate substantial revisions of the conclusions in a large number of studies that have used A-tract DNA as the bend angle reference standard in comparison measurements. Finally, any such future studies employing this comparison methodology will need to use the same sequence analyzed in the original measurements as well as replicate the original measurement conditions (e.g. ionic composition and temperature).

Bubble nucleation and cooperativity in DNA melting

Bubbles in DNA are related to fundamental processes such as duplication and transcription. Using a new ensemble technique to trap intermediate states, we present direct measurements of the average length of the denaturation bubble and the statistical weights of the bubble states in the temperature-driven melting of DNA oligomers. For a bubble flanked by double-stranded regions, we find a nucleation size of approximately 20 bases, and a broad distribution of bubble sizes. However, for bubbles opening at the ends of the molecule there is no nucleation threshold. The measured statistical weights of different conformations agree with the predictions of the thermodynamic models in the case of unzipping from the ends; however, internal bubble states are not completely described by the models. The measurements further show that, due to end effects, the melting transition becomes a two-state process only in the limit of a molecule length L approximately 1 bp.

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.

Analysis of DNA bending by transient electric birefringence

Transient electric birefringence has been used to analyze DNA bending in six restriction fragments containing 171, 174, 207, 263, 289, and 471 bp in three different low ionic strength buffers. The target fragments contain sequences corresponding to the apparent bend centers in pUC19 and Litmus 28, previously identified by the circular permutation assay (Strutz, K.; Stellwagen, N. C. Electrophoresis 1996, 17, 989-995). The target fragments migrate anomalously slowly in polyacrylamide gels and exhibit birefringence relaxation times that are shorter than those of restriction fragments of the same size, taken from nonbent regions of the same plasmids. Apparent bend angles ranging from 30 degrees to 41 degrees were calculated for the target fragments by tau-ratio method. The bend angles of four of the target fragments were independent of temperature from 4 degrees C to 20 degrees C, but decreased when the temperature was increased to 37 degrees C. The bend angles of the other two target fragments were independent of temperature over the entire range examined, 4 degrees -37 degrees C. Hence, the thermal stability of sequence-dependent bends in random-sequence DNA is variable. The bend angles of five of the six target fragments were independent of the presence or absence of Mg2+ ions in the solution, indicating most of the target fragments were stably bent or curved, rather than anisometrically flexible. Restriction fragments containing 219 and 224 bp, with sequences somewhat offset from the sequence of the 207 bp fragment, were also studied. Comparison of the tau-ratios of these overlapping fragments allowed both the bend angle and bend position to be independently determined. These methods should be useful for analyzing sequence-dependent bending in other random-sequence DNAs.

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 temperature-dependent interaction of neutral red with calf thymus DNA

Neutral red (NR) is used as a probe to study the temperature and concentration dependent interaction of a cationic dye with nucleic acid. A temperature-dependent interaction of NR with calf thymus DNA (CT DNA) has been studied by differential pulse voltammetry (DPV), UV-Visible absorption, circular dichroism (CD) and fluorescence spectroscopy. The experimental results of increasing peak current, changes in the UV-Visible absorption and fluorescence spectra of NR and decreasing the induced circular dichroism (ICD) intensity show that (i) the binding mode of NR molecules is changed from intercalating into DNA base pairs to aggregating along the DNA double helix and (ii) the orientation of NR chromophore in DNA double helix is also changed with the temperature.

RpoD promoters in Campylobacter jejuni exhibit a strong periodic signal instead of a -35 box

We have used a hidden Markov model (HMM) to identify the consensus sequence of the RpoD promoters in the genome of Campylobacter jejuni. The identified promoter consensus sequence is unusual compared to other bacteria, in that the region upstream of the TATA-box does not contain a conserved -35 region, but shows a very strong periodic variation in the AT-content and semi-conserved T-stretches, with a period of 10-11 nucleotides. The TATA-box is in some, but not all cases, preceded by a TGx, similar to an extended -10 promoter. We predicted a total of 764 presumed RpoD promoters in the C.jejuni genome, of which 654 were located upstream of annotated genes. A similar promoter was identified in Helicobacter pylori, a close phylogenetic relative of Campylobacter, but not in Escherichia coli, Vibrio cholerae, or six other Proteobacterial genomes, or in Staphylococcus aureus. We used upstream regions of high confidence genes as training data (n=529, for the C.jejuni genome). We found it necessary to limit the training set to genes that are preceded by an intergenic region of >100bp or by a gene oriented in the opposite direction to be able to identify a conserved sequence motif, and ended up with a training set of 175 genes. This leads to the conclusion that the remaining genes (354) are more rarely preceded by a (RpoD) promoter, and consequently that operon structure may be more widespread in C.jejuni than has been assumed by others. Structural predictions of the regions upstream of the TATA-box indicates a region of highly curved DNA, and we assume that this facilitates the wrapping of the DNA around the RNA polymerase holoenzyme, and offsets the absence of a conserved -35 binding motif.

Topological measurement of an A-tract bend angle: comparison of the bent and straightened states

It is well established that an A-tract imparts curvature to the DNA double helix. Constructs of such A-tracts have been used as bend standards in a large number of both structural and functional studies, and A-tracts can confer significant activation in transcription. An accurate value for the bend angle induced by an A-tract is centrally important to all such studies, but the estimates reported for the bend angle of an A-tract differ by greater than threefold. To address this problem, we have used the rotational variant method to measure the angle of DNA curvature conferred by a tract of six adenine bases (A6 tract). The original version of the method measured a protein-induced bend angle independent of external standards. It compared the effect of bent and straight forms of the sequence on the topology of a DNA plasmid in which the sequence is cloned as a series of tandem repeats. To adapt the approach to the measurement of an intrinsic bend, high temperature was used to generate the straightened reference state, with the required topological relaxation being performed by a hyperthermophile topoisomerase. Appropriate plasmids containing tandem repeats of A-tracts were constructed and topologically analyzed in this manner. The bend value measured at 4 degrees C was 26(+/-2), and decreased linearly to 17(+/-2) at 37 degrees C. The relationship to other estimates and the application of these values are discussed.

Determination of base and backbone contributions to the thermodynamics of premelting and melting transitions in B DNA

In previous papers of this series the temperature-dependent Raman spectra of poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) were used to characterize structurally the melting and premelting transitions in DNAs containing consecutive A.T and alternating A.T/T.A base pairs. Here, we describe procedures for obtaining thermodynamic parameters from the Raman data. The method exploits base-specific and backbone-specific Raman markers to determine separate thermodynamic contributions of A, T and deoxyribosyl-phosphate moieties to premelting and melting transitions. Key findings include the following: (i) Both poly(dA).poly(dT) and poly(dA-dT). poly(dA-dT) exhibit robust premelting transitions, due predominantly to backbone conformational changes. (ii) The significant van't Hoff premelting enthalpies of poly(dA).poly(dT) [DeltaH(vH)(pm) = 18.0 +/- 1.6 kcal x mol(-1) (kilocalories per mole cooperative unit)] and poly(dA-dT).poly(dA-dT) (DeltaH(vH)(pm) = 13.4 +/- 2.5 kcal x mol(-1)) differ by an amount (approximately 4.6 kcal x mol(-1)) estimated as the contribution from three-centered inter-base hydrogen bonding in (dA)(n).(dT)(n) tracts. (iii) The overall stacking free energy of poly(dA). poly(dT) [-6.88 kcal x mol(bp)(-1) (kilocalories per mole base pair)] is greater than that of poly(dA-dT). poly(dA-dT) (-6.31 kcal x mol(bp)(-1)). (iv) The difference between stacking free energies of A and T is significant in poly(dA).poly(dT) (DeltaDeltaG(st) = 0.8 +/- 0.3 kcal. mol(bp)(-1)), but marginal in poly(dA-dT).poly(dA-dT) (DeltaDeltaG(st) = 0.3 +/- 0.3 kcal x mol(bp)(-1)). (v) In poly(dA). poly(dT), the van't Hoff parameters for melting of A (DeltaH(vH)(A) = 407 +/- 23 kcal.mol(-1), DeltaS(vH)(A) = 1166 +/- 67 cal. degrees K(-1) x mol(-1), DeltaG(vH(25 degrees C))(A) = 60.0 +/- 3.2 kcal x mol(-1)) are clearly distinguished from those of T (DeltaH(vH)(T) = 185 +/- 38 kcal x mol(-1), DeltaS(vH)(T) = 516 +/- 109 cal. degrees K(-1) x mol(-1), DeltaG(vH(25 degrees C))(T) = 27.1 +/- 5.5 kcal x mol(-1)). (vi) Similar relative differences are observed in poly(dA-dT). poly(dA-dT) (DeltaH(vH)(A) = 333 +/- 54 kcal x mol(-1), DeltaS(vH)(A) = 961 +/- 157 cal. degrees K(-1) x mol(-1), DeltaG(vH(25 degrees C))(A) = 45.0 +/- 7.6 kcal x mol(-1); DeltaH(vH)(T) = 213 +/- 30 kcal x mol(-1), DeltaS(vH)(T) = 617 +/- 86 cal. degrees K(-1) x mol(-1), DeltaG(vH(25 degrees C))(T) = 29.3 +/- 4.9 kcal x mol(-1)). The methodology employed here distinguishes thermodynamic contributions of base stacking, base pairing and backbone conformational ordering in the molecular mechanism of double-helical B DNA formation.

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

Melting of the solvent structure around a RNA duplex: a molecular dynamics simulation study

From three 2.4-ns molecular dynamics simulations of the r(CpG)(12) duplex conducted at 5, 25 and 37 degrees C, a strong temperature dependence of the dynamics of the water molecules and ions located in the first nucleic acid coordination shell is observed. At 5 degrees C, the highest residence times of bound water molecules exceed 1 ns while, at 37 degrees C, they decrease to 0.5 ns in agreement with available NMR data. Similar temperature dependencies are observed for the potassium ions bound to the duplex. In this temperature range, the structure of the RNA helix remains essentially unchanged. Thus, the observed alterations correspond to a 'premelting' of the solvent structure around the duplex. It is proposed that, before the nucleic acid structure melts, the entropy of the solvent increases to a point where it is no longer compensated by the enthalpic contribution of solute-solute and solute-solvent interactions. At this stage, the weakest structural elements start to melt. In other terms, the experimentally observed melting processes are preceded by a melting of the more labile solvent structure.

Temperature dependence of the Raman spectrum of DNA. II. Raman signatures of premelting and melting transitions of poly(dA).poly(dT) and comparison with poly(dA-dT).poly(dA-dT)

The temperature dependence of the Raman spectrum of poly(dA).poly(dT) (dA: deoxyadenosine; dT: thymidine), a model for DNA containing consecutive adenine.thymine (A.T) pairs, has been analyzed using a spectrometer of high spectral precision and sensitivity. Three temperature intervals are distinguished: (a) premelting (10 < t < 70 degrees C), in which the native double helix is structurally altered but not dissociated into single strands; (b) melting (70 < t < 80 degrees C), in which the duplex is dissociated into single strands; and (c) postmelting (80 < t degrees C), in which no significant structural change can be detected. The distinctive Raman difference signatures observed between 10 and 70 degrees C and between 70 and 80 degrees C are interpreted in terms of the structural changes specific to premelting and melting transitions, respectively. Premelting alters the low-temperature conformation of the deoxyribose-phosphate backbone and eliminates base hydrogen bonding that is distinct from canonical Watson-Crick hydrogen bonding; these premelting perturbations occur without disruption of base stacking. Conversely, melting eliminates canonical Watson-Crick pairing and base stacking. The results are compared with those reported previously on poly(dA-dT).poly(dA-dT), the DNA structure consisting of alternating A.T and T.A pairs (L. Movileanu, J. M. Benevides, and G. J. Thomas, Jr. Journal of Raman Spectroscopy, 1999, Vol. 30, pp. 637-649). Poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) exhibit strikingly dissimilar temperature-dependent Raman profiles prior to the onset of melting. However, the two duplexes exhibit very similar melting transitions, including the same Raman indicators of ruptured Watson-Crick pairing, base unstacking and collapse of backbone order. A detailed analysis of the data provides a comprehensive Raman assignment scheme for adenosine and thymidine residues of B-DNA, delineates Raman markers diagnostic of consecutive A.T and alternating A.T/T.A tracts of DNA, and identifies the distinct Raman difference signatures for premelting and melting transitions in the two types of sequences.

UV resonance Raman and circular dichroism studies of a DNA duplex containing an A(3)T(3) tract: evidence for a premelting transition and three-centered H-bonds

The presence of A(n) and A(n)T(n) tracts in double-helical sequences perturbs the structural properties of DNA molecules, resulting in the formation of an alternate conformation to standard B-DNA known as B'-DNA. Evidence for a transition occurring prior to duplex melting in molecules containing A(n) tracts was previously detected by circular dichroism (CD) and calorimetric studies. This premelting transition was attributed to a conformational change from B'- to B-DNA. Structural features of A(n) and A(n)T(n) tracts revealed by X-ray crystallography include a large degree of propeller twisting of adenine bases, narrowed minor grooves, and the formation of three-centered H-bonds between dA and dT bases. We report UV resonance Raman (UVRR) and CD spectroscopic studies of two related DNA dodecamer duplexes, d(CGCAAATTTGCG)(2) (A(3)T(3)) and d(CGCATATATGCG)(2) [(AT)(3)]. These studies address the presence of three-centered H-bonds in the B' conformation and gauge the impact of these putative H-bonds on the structural and thermodynamic properties of the A(3)T(3) duplex. UVRR and CD spectra reveal that the premelting transition is only observed for the A(3)T(3) duplex, is primarily localized to the dA and dT bases, and is associated with base stacking interactions. Spectroscopic changes associated with the premelting transition are not readily detectable for the sugar-phosphate backbone or the cytosine and guanosine bases. The temperature-dependent concerted frequency shifts of dA exocyclic NH(2) and dT C4=O vibrational modes suggest that the A(3)T(3) duplex forms three-centered hydrogen bonds at low temperatures, while the (AT)(3) duplex does not. The enthalpy of this H-bond, estimated from the thermally induced frequency shift of the dT C4=O vibrational mode, is approximately 1.9 kJ/mol or 0.46 kcal/mol.

Raman spectroscopy of an O(2)-Co(II)bleomycin-calf thymus DNA adduct: alternate polymer conformations

Bleomycin (Blm) is an antitumor agent which binds to specific sequences of DNA and as HO(2)-Fe(III)Blm causes single and double strand cleavage. In the present investigation, binding of O(2)-Co(II)Blm to a native DNA polymer, calf thymus DNA, was examined using conventional Raman spectroscopy. O(2)-Co(II)Blm is a model for O(2)-Fe(II)Blm, the direct precursor of HO(2)-Fe(III)Blm. Although the DNA polymer retained a predominant B-form structure, Raman spectral evidence was obtained for localized structural changes to A, C and Z-DNA forms. The presence of these alternate DNA forms within B-DNA implied the presence of B/A, B/C and B/Z junctions. The observed changes in DNA secondary structure were attributed to perturbation of structural water resulting from binding of O(2)-Co(II)Blm within the minor groove.

Analysis and interpretation of ligand-DNA binding isotherms

Binding studies provide information of fundamental and central importance for the complete understanding of ligand-DNA interactions. Studies of ligand binding to long natural DNA samples, to synthetic deoxypolynucleotides of simple repeating sequence, and to oligonucleotides of defined sequence are all needed to begin to understand the interaction in detail. Binding studies provide entry into the thermodynamics of the DNA interactions, which in turn provides great insight into the molecular forces that drive the binding process. This chapter summarizes both model-dependent and -independent approaches for the analysis and interpretation of binding isotherms, and should serve as a concise guide for handling experimental data.

Magnesium increases the curvature of duplex DNA that contains dA tracts

Distinct structural features of DNA, such as the curvature of dA tracts, are important in the recognition, packaging, and regulation of DNA. Physiologically relevant concentrations of magnesium have been found to enhance the curvature of dA tract DNAs, as monitored by solution-state NMR, indicating that the structure of DNA depends on the cations present in solution. A model is presented which accounts for the sequence-dependent effects of magnesium on DNA curvature as well as for the previously known sequence-independent effect on DNA flexibility.

The presence of two modes of binding to calf thymus DNA by metal-free bleomycin: a low frequency Raman study

Double-stranded DNA is targeted by bleomycin in cancer cells and ambiguity exists as to its mode of DNA binding. A conventional Raman study was performed on drug/DNA complexes in which the low frequency spectral region (560-930 cm(-1)) was examined at two temperatures (19 and 30 degrees C). At 30 degrees C, a global Raman hypochromism was observed consistent with partial intercalation of the bithiazole moiety. At 19 degrees C, Raman hypochromism (increased base pair stacking) was detected for bands associated with GC base pairs while Raman hyperchromism (base pair destacking) was evident for bands associated with AT base pairs. These results suggest that intercalation of the bithiazole moiety occurs with greater disruption of the more efficiently stacked AT base pairs at the lower temperature. Evidence for minor groove binding was indicated by an increase in the population of bands corresponding to C3' endo sugar conformations resulting from drug induced local desolvation of the DNA polymer.

A systematic approach toward the analysis of drug-DNA interactions using Raman spectroscopy: the binding of metal-free bleomycins A(2) and B(2) to calf thymus DNA

Bleomycins A(2) and B(2) are the two active components in the antineoplastic drug Blenoxane. DNA is targeted by this drug in cancer cells and the mode of action of this drug involves DNA binding. Ambiguity exists as to the way in which bleomycin binds to DNA. Raman spectroscopy was used to examine both calf thymus DNA and a bleomycin/DNA complex at two temperatures. A curvefitting technique was applied to these spectra for a spectral region obscured by many overlapping bands associated with the nucleotide bases in order to derive information about frequencies, bandwidths, and intensities of the vibrational modes in this region. This allowed identification and analysis of bands associated with specific assigned nucleotide base residues. Upon binding of bleomycin, several significant changes in bandwidth, intensities, and frequencies relative to uncomplexed DNA were observed consistently at both higher (30 degrees C) and lower (19 degrees C) temperature. The data presented here support at least a partial intercalation mode of binding for bleomycin that is temperature dependent and more pronounced at the more physiologically relevant temperature of 30 degrees C.

Genome organisation and chromatin structure in Escherichia coli

We have analysed the complete sequence of the Escherichia coli K12 isolate MG1655 genome for chromatin-associated protein binding sites, and compared the predicted location of predicted sites with experimental expression data from 'DNA chip' experiments. Of the dozen proteins associated with chromatin in E. coli, only three have been shown to have significant binding preferences: integration host factor (IHF) has the strongest binding site preference, and FIS sites show a weak consensus, and there is no clear consensus site for binding of the H-NS protein. Using hidden Markov models (HMMs), we predict the location of 608 IHF sites, scattered throughout the genome. A subset of the IHF sites associated with repeats tends to be clustered around the origin of replication. We estimate there could be roughly 6000 FIS sites in E. coli, and the sites tend to be localised in two regions flanking the replication termini. We also show that the regions upstream of genes regulated by H-NS are more curved and have a higher AT content than regions upstream of other genes. These regions in general would also be localised near the replication terminus.

A statistical mechanical model for predicting B-DNA curvature and flexibility

A statistical mechanical model taking into account the symmetric twisting, tilting, sliding fluctuations and asymmetric rolling fluctuations has been proposed to predict the macroscopic curvature and flexibility of B-DNA. Based on the statistical data of structural parameters of double helix in nucleic acid database and the related theoretical analysis, the equilibrium angular parameters (Omega, rho and tau) describing the orientation of successive base-pair planes, the translation parameters (D(y)) along the long axis of neighboring base-pair step and the corresponding force constants are arranged for ten dimers appropriately. Under the assumption of independent angular parameters, independent base-pair steps and a simple energy function, we can calculate the macroscopic curvature and the flexibility of DNA sequences through the transformation matrix <R> and the Boltzmann ensemble average. The predictions on curvature and flexibility of DNA have been compared with the corresponding experimental data. The agreement is remarkably good. It is demonstrated that the lowering of the temperature does increase the DNA curvature.

The curvature of dA tracts is temperature dependent

The curvature of dA tracts has been proposed to be important in the recognition, packaging, and regulation of DNA. The effects of dA tracts on the gel mobility, rate of cyclization, and other properties of DNA have been extensively studied. The consensus value for the curvature induced by a single dA tract is about 18 degrees. There are two main competing models for the origin of the curvature of dA tracts. One model assigns the central role to sequence-dependent steric clashes and the other to sequence dependent interactions with cations. The temperature dependence of the shape functions, the molecule specific part of the diffusion coefficients, of a set of six DNAs has been examined here. The set contains DNAs with dA tracts in or out of phase with respect to the helical repeat as well as those with scrambled dA-dT regions. The results show that the curvature of dA tracts is highly temperature dependent and that the curvature is largely melted out by 40 degrees C. The curvature melts out before there is significant premelting, or breathing of the dA tracts or the scrambled dA-dT regions. The curvature does not appear to reach a plateau value at low temperatures. A qualitative model for the melting of the curvature of dA tracts is proposed.

Ecologic genomics of DNA: upstream bending in prokaryotic promoters

After our analysis of the distribution of predicted intrinsic curvature along all available complete prokaryotic genomes, the genomes were divided into two groups. Curvature distribution in all prokaryotes of the first group indicated a substantial fraction of promoters characterized by intrinsic DNA curvature located within or upstream of the promoter region. We did not find this peculiar DNA curvature distribution in prokaryotes in the second group. Remarkably, all bacteria of the first group were mesophilic, whereas many prokaryotes of the second group were hyperthermophilic. We hypothesize that DNA curvature plays a biologic role in gene regulation in mesophilic as opposed to hyperthermophilic prokaryotes, i.e., DNA curvature presumably has a functional adaptive significance determined by temperature selection.

Conformational isomers of curved DNA molecules can be observed by polyacrylamide gel electrophoresis

Dimers, trimers and higher multimers of two 147-base pair restriction fragments called 12 A and 12B, obtained from the MspI digest of plasmid pBR322, migrate as sharp bands in agarose and dilute polyacrylamide gels, indicating that they are homogeneous in molecular weight. However, the electrophoretic bands corresponding to multimers of the curved fragment 12A are split into sharp sub-bands in more concentrated polyacrylamide gels. The relative intensities and spacing of the sub-bands depend on the number of monomers in the multimer, the pH of the buffer, and the presence or absence of divalent cations in the solution. Since band splitting is not observed for the normal 12B multimers under any gel-running conditions, the sub-bands observed for multimers of the curved fragment 12A must be attributed to conformational isomers which are in slow exchange on the electrophoretic time scale. Band splitting is also observed for multimers of a curved DNA fragment containing the kinetoplast bending locus and for plasmid pUC19 linearized by digestion with certain restriction enzymes. Plasmid pUC19 contains two nearly equidistant regions of intrinsic curvature (Strutz, K., Stellwagen, N. C., Electrophoresis 1996, 17, 989-995). Hence, DNA molecules containing two or more regions of curvature exist as discrete subpopulations of conformational isomers which can be observed as separate bands migrating in polyacrylamide gels.

The question of long-range allosteric transitions in DNA

The question of long-range allosteric transitions of DNA secondary structure and their possible involvement in transcriptional activation is discussed in the light of new results. A variety of recent evidence strongly supports a fluctuating long-range description of DNA secondary structure. Balanced equilibria between two or more different secondary structures, and the occurrence of very large domain sizes, have been documented in several instances. Long-range allosteric effects stemming from changes in sequence or secondary structure over a small region of the DNA have been observed to extend over distances up to hundreds of base pairs in some cases. The discovery that coherent bending strain beyond a threshold level in small (N < or = 250 base pairs (bp)] circular DNAs significantly alters the DNA secondary structure has important implications, especially for transcriptional activators that either bend the DNA directly or are involved in the formation of DNA loops of sufficiently small size (N < or = 250 bp). Whether the RNA polymerase is activated primarily via protein: protein contacts, as is widely believed, or instead via a bend-induced allosteric transition of the DNA in such a small loop, is now an open question. Binding of the transcriptional activator Sp1 to linear DNA induces a remarkably long-range change in its secondary structure, and catabolite activator protein binding to a supercoiled DNA behaves similarly, though possibly for different reasons. Compelling evidence for a bend-induced long-range structural transmission effect of the transcriptional activator integration host factor on RNA polymerase activity was recently reported. These results may augur a new paradigm in which allosteric transitions of duplex DNA, as well as of the proteins, are involved in the regulation of transcription.