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

DNA A-tracts are not curved in solutions containing high concentrations of monovalent cations

The intrinsic curvature of seven 98 bp DNA molecules containing up to four centrally located A6-tracts has been measured by gel and capillary electrophoresis as a function of the number and arrangement of the A-tracts. At low cation concentrations, the electrophoretic mobility observed in polyacrylamide gels and in free solution decreases progressively with the increasing number of phased A-tracts, as expected for DNA molecules with increasingly curved backbone structures. Anomalously slow electrophoretic mobilities are also observed for DNA molecules containing two pairs of phased A-tracts that are out of phase with each other, suggesting that out-of-phase distortions of the helix backbone do not cancel each other out. The mobility decreases observed for the A-tract samples are due to curvature, not cation binding in the A-tract minor groove, because identical free solution mobilities are observed for a molecule with four out-of-phase A-tracts and one with no A-tracts. Surprisingly, the curvature of DNA A-tracts is gradually lost when the monovalent cation concentration is increased to ∼200 mM, regardless of whether the cation is a hydrophilic ion like Na+, NH4+, or Tris+ or a hydrophobic ion like tetrabutylammonium. The decrease in A-tract curvature with increasing ionic strength, along with the known decrease in A-tract curvature with increasing temperature, suggests that DNA A-tracts are not significantly curved under physiological conditions.

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.

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes -2.3 kcal mol(-1) to -0.9 kcal mol(-1) (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of +1.3 kcal mol(-1). Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.

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.

Escherichia coli DNA adenine methyltransferase: the structural basis of processive catalysis and indirect read-out

We have investigated the structural basis of processive GATC methylation by the Escherichia coli DNA adenine methyltransferase, which is critical in chromosome replication and mismatch repair. We determined the contribution of the orthologically conserved phosphate interactions involving residues Arg(95), Asn(126), Asn(132), Arg(116), and Lys(139), which directly contact the DNA outside the cognate recognition site (GATC) to processive catalysis, and that of residue Arg(137), which is not conserved and contacts the DNA backbone within the GATC sequence. Alanine substitutions at the conserved positions have large impacts on processivity yet do not impact k(cat)/K(m)(DNA) or DNA affinity (K(D)(DNA)). However, these mutants cause large preferences for GATC sites varying in flanking sequences when considering the pre-steady state efficiency constant k(chem)/K(D)(DNA). These changes occur mainly at the level of the methylation rate constant, which results in the observed decreases in processive catalysis. Thus, processivity and catalytic efficiency (k(cat)/K(m)(DNA)) are uncoupled in these mutants. These results reveal that the binding energy involved in DNA recognition contributes to the assembly of the active site rather than tight binding. Furthermore, the conserved residues (Arg(95), Asn(126), Asn(132), and Arg(116)) repress the modulation of the response of the enzyme to flanking sequence effects. Processivity impacted mutants do not show substrate-induced dimerization as is observed for the wild type enzyme. This study describes the structural means by which an enzyme that does not completely enclose its substrate has evolved to achieve processive catalysis, and how interactions with DNA flanking the recognition site alter this processivity.

Topological analysis of plasmid chromatin from yeast and mammalian cells

Yeast has proven to be a powerful system for investigation of chromatin structure. However, the extent to which yeast chromatin can serve as a model for mammalian chromatin is limited by the significant number of differences that have been reported. To further investigate the structural relationship between the two chromatins, we have performed a DNA topological analysis of pRSSVO, a 5889 base-pair plasmid that can replicate in either yeast or mammalian cells. When grown in mammalian cells, pRSSVO contains an average of 33 negative supercoils, consistent with one nucleosome per 181 bp. This is close to the measured nucleosome repeat length of 190 bp. However, when grown in yeast cells, pRSSVO contains an average of only 23 negative supercoils, which is indicative of only one nucleosome per 256 bp. This is dramatically different from the measured nucleosome repeat length of 165 bp. To account for these observations, we suggest that yeast chromatin is composed of relatively short ordered arrays of nucleosomes with a repeat of 165 bp, separated by substantial gaps, possibly corresponding to regulatory regions.

Effect of organic cosolvents on the free solution mobility of curved and normal DNA molecules

The free solution mobilities of curved and normal 199-bp DNA fragments have been measured in buffer solutions containing various quantities of the organic cosolvents methanol, ethanol, 2-propanol, 2-methyl-2,4-pentanediol (MPD), ethylene glycol, and ACN, using CE. The curved fragment, taken from the VP1 gene of SV40, contains five unevenly spaced A- and T-tracts in a centrally located "curvature module"; the A- and T-tracts have been mutated to other sequences in the normal 199-bp fragment. The free solution mobility of the curved 199-bp fragment is significantly lower than that of its normal counterpart in aqueous solutions [Stellwagen, E., Lu, Y. J., Stellwagen, N. C., Nucleic Acids Res. 2005, 33, 4425-4432]. The mobilities of both the curved and normal fragments decrease with increasing cosolvent concentration, due to the effect of the cosolvent on the viscosity and dielectric constant of the solution. The mobility differences between the curved and normal 199-bp fragments and the mobility ratios decrease approximately linearly with the increasing mole fraction of cosolvent in the solution. Hence, MPD and other organic cosolvents affect DNA electrophoretic mobility by a common mechanism, most likely the preferential hydration of the DNA surface that occurs in aqueous cosolvents. The gradual loss of the anomalously slow mobility of the curved 199-bp fragment with increasing cosolvent concentration, combined with other data in the literature, suggests that preferential hydration gradually widens the narrow A-tract minor groove, releasing site-bound counterions in the minor groove and shifting the conformation toward that of normal DNA.

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

Purification and Characterization of the Simian Virus 40 Transcription Elongation Complex

The transcriptional regulatory region of the simian virus 40 minichromosome that is being transcribed in the cell is nucleosome-free, while that of the non-transcribed minichromosome is nucleosome covered. Although additional studies have shown that the two structures are otherwise similar, the precision of these indirect studies has not been sufficient to determine if the transition between the two involves nucleosome displacement or nucleosome sliding. In order to address this question directly, we have developed a new function-based affinity isolation method that is capable of purifying the native transcription elongation complex of a single gene from mammalian cells. The simian virus 40 transcription elongation complex was purified by this method and the topological linking number of its DNA was compared directly to that of the bulk, non-transcribed minichromosome. The results show that the two types of minichromosome contain the same number of nucleosomes as well as nucleosomal structure. These findings indicate that interconversion between the non-transcribing and transcribing states is accomplished by a remodeling event involving nucleosome sliding rather than nucleosome displacement.

Gapped DNA and cyclization of short DNA fragments

We use the cyclization of small DNA molecules, approximately 200 bp in length, to study conformational properties of DNA fragments with single-stranded gaps. The approach is extremely sensitive to DNA conformational properties and, being complemented by computations, allows a very accurate determination of the fragment's conformational parameters. Sequence-specific nicking endonucleases are used to create the 4-nt-long gap. We determined the bending rigidity of the single-stranded region in the gapped DNA. We found that the gap of 4 nt in length makes all torsional orientations of DNA ends equally probable. Our results also show that the gap has isotropic bending rigidity. This makes it very attractive to use gapped DNA in the cyclization experiments to determine DNA conformational properties, since the gap eliminates oscillations of the cyclization efficiency with the DNA length. As a result, the number of measurements is greatly reduced in the approach, and the analysis of the data is greatly simplified. We have verified our approach on DNA fragments containing well-characterized intrinsic bends caused by A-tracts. The obtained experimental results and theoretical analysis demonstrate that gapped-DNA cyclization is an exceedingly sensitive and accurate approach for the determination of DNA bending.

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.

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

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.

Temperature sensing by the dsrA promoter

Synthesis of the small regulatory RNA DsrA is under temperature control. The minimal dsrA promoter of 36 bp contains sufficient information to ensure such regulation. In vivo, we have analyzed the critical elements responsible for the temperature control of dsrA by using a collection of chimeric promoters combining various elements of the dsrA promoter and the lacUV5 promoter, which does not respond to temperature. Our results favor an RNA polymerase-DNA interaction model instead of a trans-acting factor for temperature regulation. While all of the elements of the dsrA promoter contribute to temperature-sensitive expression, the sequence of the -10 box and the spacer region are the essential elements for the thermal response of the dsrA promoter. The proper context for these promoter elements, including at least one of the flanking elements, the -35 region or the start site region, is also required. Point mutations demonstrate that the sequence of the -10 box imposes constraints on the length and the sequence of the spacer and/or its AT richness, even at low temperature. These results show a complex interdependence of different regions in the promoter for temperature regulation.