Effect of ethidium on the torsion constants of linear and supercoiled DNAs

A quantitative model of a cooperative two-state equilibrium in DNA: experimental tests, insights, and predictions

Quantitative parameters for a two-state cooperative transition in duplex DNAs were finally obtained during the last 5 years. After a brief discussion of observations pertaining to the existence of the two-state equilibrium per se, the lengths, torsion, and bending elastic constants of the two states involved and the cooperativity parameter of the model are simply stated. Experimental tests of model predictions for the responses of DNA to small applied stretching, twisting, and bending stresses, and changes in temperature, ionic conditions, and sequence are described. The mechanism and significance of the large cooperativity, which enables significant DNA responses to such small perturbations, are also noted. The capacity of the model to resolve a number of long-standing and sometimes interconnected puzzles in the extant literature, including the origin of the broad pre-melting transition studied by numerous workers in the 1960s and 1970s, is demonstrated. Under certain conditions, the model predicts significant long-range attractive or repulsive interactions between hypothetical proteins with strong preferences for one or the other state that are bound to well-separated sites on the same DNA. A scenario is proposed for the activation of the ilvPG promoter on a supercoiled DNA by integration host factor.

A structural transition in duplex DNA induced by ethylene glycol

The twist energy parameter ( E T) that governs the supercoiling free energy, and the linking difference (Delta l) are measured for p30delta DNA in solutions containing 0-40 w/v % ethylene glycol (EG). A plot of E T vs -ln a w, where a w is the water activity, displays the full (reverse) sigmoidal profile of a discrete structural transition. A general theory for the effect of added osmolyte on a cooperative structural transition between two duplex states, 1 right arrow over left arrow 2, is formulated in terms of parameters applicable to individual base-pair subunits. The resulting fraction of base pairs in the 2-state ( f 2 (0)) is incorporated into expressions for the effective torsion and bending elastic constants, the effective twist energy parameter ( E T (eff)), and the change in intrinsic twist (delta l 0). Fitting the expression for E T (eff) to the measured E T values yields reasonably unambiguous estimates of E T 1 and E T 2 , the midpoint value (ln a w) 1/2, and the midpoint slope ( partial differential E T/ partial differential ln a w) 1/2, but does not yield unambiguous estimates of the equilibrium constant ( K 0), the difference in DNA-water preferential interaction coefficient (DeltaGamma), or the inverse cooperativity parameter ( J). Fitting a noncooperative model (assumed J = 1.0) to the data yields K 0 = 0.067 and DeltaGamma = -30.0 per base pair (bp). Essentially equivalent fits are provided by models with a wide range of correlated J, DeltaGamma, and K 0 values. Other results favor DeltaGamma in the range -1.0 to 0, which then requires K 0 > or = 0.914, and a cooperativity parameter, 1/ J > or = 30.0 bp. The measured delta l 0 and circular dichroism (CD) at 272 nm are found to be compatible with curves predicted using the same f 2 (0) values that best-fit the E T data. At least 7-10% of the base pairs are inferred to exist in the 2-state in 0.1 M NaCl in the complete absence of added osmolyte. Compared with the 1-state, the 2-state has a approximately 2.0- to 2.1-fold greater torsion elastic constant, a approximately 0.70-fold smaller bending elastic constant, a approximately 0.91-fold smaller E T value, a approximately 0.2% lower intrinsic twist, a somewhat lower CD near both 272 and 245 nm, and less water and/or more EG in its neighborhood. However, the relative change in preferential interaction coefficient associated with the transition is likely rather slight.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30δ DNA

Upon increasing the concentration of ethylene glycol (EG) at 37 degrees C, the twist energy parameter, E(T), which governs the supercoiling free energy, was recently found to undergo a decreasing (or reverse) sigmoidal transition with a midpoint near 20 w/v % EG. In this study, the effects of adding 20 w/v % EG on the torsion elastic constant (alpha) of linear p30delta DNA and on the hydrodynamic radius (R(H)) of a synthetic 24 bp duplex DNA were examined at both 40 and 20 degrees C. The time-resolved fluorescence intensity and fluorescence polarization anisotropy (FPA) of intercalated ethidium were measured in order to assess the effects of 20 w/v % EG on: (1) alpha; (2) R(H); (3) the lifetimes of intercalated and non-intercalated dye; (4) the amplitude of dye wobble in its binding site; and (5) the binding constant for intercalation. The effects of 20 w/v % EG on the circular dichroism (CD) spectrum of the DNA and on the emission spectrum of the free dye were also measured. At 40 degrees C, addition of 20 w/v % EG caused a substantial (1.27- to 1.35-fold) increase in alpha, a significant change in the CD spectrum, and a very small, marginally significant increase in R(H), but little or no change in the amplitude of dye wobble in its binding site or the lifetime of intercalated dye. Together with previously reported measurements of E(T), these results imply that the bending elastic constant of DNA is significantly decreased by 20 w/v % EG at 40 degrees C. At 20 degrees C, addition of 20 w/v % EG caused a marginally significant decrease in alpha and very little change in any other measured properties. Also at 20 degrees C, addition of 30 w/v % betaine caused a marginally significant increase in alpha and significant but modest change in the CD spectrum, but very little change in any other properties.

Torsional rigidities of weakly strained DNAs

Measurements on unstrained linear and weakly strained large (> or =340 bp) circular DNAs yield torsional rigidities in the range C = 170-230 fJ fm. However, larger values, in the range C = 270-420 fJ fm, are typically obtained from measurements on sufficiently small (< or =247 bp) circular DNAs, and values in the range C = 300-450 fJ fm are obtained from experiments on linear DNAs under tension. A new method is proposed to estimate torsional rigidities of weakly supercoiled circular DNAs. Monte Carlo simulations of the supercoiling free energies of solution DNAs, and also of the structures of surface-confined supercoiled plasmids, were performed using different trial values of C. The results are compared with experimental measurements of the twist energy parameter, E(T), that governs the supercoiling free energy, and also with atomic force microscopy images of surface-confined plasmids. The results clearly demonstrate that C-values in the range 170-230 fJ fm are compatible with experimental observations, whereas values in the range C > or = 269 fJ fm, are incompatible with those same measurements. These results strongly suggest that the secondary structure of DNA is altered by either sufficient coherent bending strain or sufficient tension so as to enhance its torsional rigidity.

Utility and considerations of donor-donor energy migration as a fluorescence method for exploring protein structure-function

This review aims at surveying the use of electronic energy transport between chemically identical fluorophores (i.e. donors) in studies of various protein systems. Applications of intra- and interprotein energy migration are presented that make use of polarised steady-state and time-resolved fluorescence spectroscopic techniques. The donor-donor energy migration (DDEM) and the partial donor-donor energy migration (PDDEM) models for calculating distances between donor groups are exposed together with the most recent development of an extended Forster theory (EFT). Synthetic fluorescence depolarisation data that mimic time-correlated single photon counting experiments were generated using the EFT, and then further re-analysed by the different models. The results obtained were compared with the known parameters used to generate EFT data. Aspects on how to adopt the EFT in the analyses of time-correlated single photon counting experiments are also presented, as well as future aspects on using energy migration for examining protein structure.

Unbiased rotational moves for rigid-body dynamics

We introduce an unbiased protocol for performing rotational moves in rigid-body dynamics simulations. This approach-based on the analytic solution for the rotational equations of motion for an orthogonal coordinate system at constant angular velocity-removes deficiencies that have been largely ignored in Brownian dynamics simulations, namely errors for finite rotations that result from applying the noncommuting rotational matrices in an arbitrary order. Our algorithm should thus replace standard approaches to rotate local coordinate frames in Langevin and Brownian dynamics simulations.

Monte Carlo simulations of locally melted supercoiled DNAs in 20 mM ionic strength

Mesoscopic models of unmelted and locally melted supercoiled DNAs in 20 mM ionic strength are simulated over a range of linking difference from deltal = 0 to -26 turns, or superhelix density from sigma = 0 to -0.062. A domain containing m = 0, 28, or 56 melted basepairs (out of 4349 total) is modeled simply by a region of suitable length with substantially reduced torsion and bending elastic constants. Average structural properties are calculated from the saved configurations, and a reversible work protocol is used to calculate the supercoiling free energy, The cross-writhe between duplex and melted regions (defined herein) is found to be negligibly small. The total writhe, radius of gyration, and ordered elements of the diagonalized inertial tensor are found to be nearly universal functions of the residual linking difference (deltal(r)) associated with the duplex region, independent of m. However, deformability of the tertiary structure, as manifested by the variance of those same properties, is not a universal function of deltal(r)), but depends upon m.delta (SC) varies with deltal(r)) more strongly than deltal(r)) (2)due to the low ionic strength. The twist energy parameter, E (T) obtained from the simulated delta G(SC), deltal(r)), and net twisting strain of the melted region T (D), is found to be independent of m, hence also of the torsion and bending elastic constants of the melted region. However, E(T) increases linearly with -deltalr), which leads to 1). a small overestimation of E (T) for any given deltal(r)) when E(T) is determined from the observed deltal and deltal (r) by the protocol of Bauer and Benham; and 2). a significant enhancement of the apparent slope, -dE(T)/d(T), obtained via the protocol of Bauer and Benham, relative to the actual slope at fixed delta l(r). After taking these two effects into account, the theoretical and experimental values E(T) and -dE(T)/d(T) values agree rather well. For the larger deltal the melted regions are found preferentially in the linker domains between interwound arms, rather than in the apical regions at the ends of interwound arms.

Effects of small neutral osmolytes on the supercoiling free energy and intrinsic twist of p30? DNA

Both theory and experiments are employed to investigate the effects of small neutral osmolytes on the average intrinsic twist (l0), the torsion and bending elastic constants, and the twist energy parameter (ET) that governs the supercoiling free energy. The experimental data for ethylene glycol and acetamide at 37 degrees C suggest, and are interpreted in terms of, a model wherein the DNA exhibits an equilibrium between two distinct conformational states that possess different numbers of bound water molecules and exhibit different intrinsic twists and torsion and bending elastic constants. Expressions are derived to relate the effective ET and l0 to the equilibrium constant, water activity (aw), and number (n) of bound water molecules released per cooperative domain undergoing the two-state transition. The variations of l0 and ET with -ln(aw) are similar for acetamide and ethylene glycol at 37 degrees C. Fitting the theory to those data yields the range n = 103-125 for ethylene glycol and n = 71-113 for acetamide, depending on the assumed value of ET for the dehydrated state. The cooperative domain size of the two-state transition is estimated to exceed about 25-30 base pairs (bp). Between 0 and 19.4 w/v % ethylene glycol, the torsion elastic constant, measured by time-resolved fluorescence polarization anisotropy (FPA), increases by 1.37-fold, whereas the measured ET decreases by 1.15-fold over that same range. The implied decrease in bending rigidity over that range is by a factor of about 0.7. The variations of l0 and ET with increasing -ln(aw) due to added ethylene glycol at 37 degrees C are far smaller than the corresponding variations observed previously at 14 and 15 degrees C. However, at 21 degrees C, upon adding either ethylene glycol or acetamide, l0 and ET initially decline steeply with increasing -ln(aw), with slopes possibly comparable to those seen at 14 and 15 degrees C, but then flatten out and follow curves similar to those at 37 degrees C. Possible origins of such mixed behavior are discussed. The effects of betaine at both 37 and 21 degrees C differ qualitatively and quantitatively in various respects from those of ethylene glycol and acetamide. Upon adding sucrose, l0 initially jumps to higher plateaus at both 37 and 21 degrees C, but its effects on ET cannot be reliably assessed, due to the limited range of -ln(aw).

Monte Carlo simulations of supercoiled DNAs confined to a plane

Recent advances in atomic force microscopy (AFM) have enabled researchers to obtain images of supercoiled DNAs deposited on mica surfaces in buffered aqueous milieux. Confining a supercoiled DNA to a plane greatly restricts its configurational freedom, and could conceivably alter certain structural properties, such as its twist and writhe. A program that was originally written to perform Monte Carlo simulations of supercoiled DNAs in solution was modified to include a surface potential. This potential flattens the DNAs to simulate the effect of deposition on a surface. We have simulated transfers of a 3760-basepair supercoiled DNA from solution to a surface in both 161 and 10 mM ionic strength. In both cases, the geometric and thermodynamic properties of the supercoiled DNAs on the surface differ significantly from the corresponding quantities in solution. At 161 mM ionic strength, the writhe/twist ratio is 1.20-1.33 times larger for DNAs on the surface than for DNAs in solution and significant differences in the radii of gyration are also observed. Simulated surface structures in 161 mM ionic strength closely resemble those observed by AFM. Simulated surface structures in 10 mM ionic strength are similar to a minority of the structures observed by AFM, but differ from the majority of such structures for unknown reasons. In 161 mM ionic strength, the internal energy (excluding the surface potential) decreases substantially as the DNA is confined to the surface. Evidently, supercoiled DNAs in solution are typically deformed farther from the minimum energy configuration than are the corresponding surface-confined DNAs. Nevertheless, the work (Delta A(int)) done on the internal coordinates, which include uniform rotations at constant configuration, during the transfer is positive and 2.6-fold larger than the decrease in internal energy. The corresponding entropy change is negative, and its contribution to Delta A(int) is positive and exceeds the decrease in internal energy by 3.6 fold. The work done on the internal coordinates during the solution-to-surface transfer is directed primarily toward reducing their entropy. Evidently, the number of configurations available to the more deformed solution DNA is vastly greater than for the less deformed surface-confined DNA.

Effect of glycerol-induced DNA conformational change on the separation of DNA fragments by capillary electrophoresis

The effect on DNA separation of adding glycerol to the running buffer was studied using linear polyacrylamide (LPA) or poly(ethylene oxide) (PEO) as separation medium. For both LPA and PEO, it was found that the addition of 25% (v/v) glycerol to the running buffer enhanced the separation of large double-stranded DNA fragments and increased the migration time. The two buffers used, 1 x TBE (Tris-boric acid-EDTA) and 1 x TTE (Tris-N-tris(hydroxymethyl)methyl-3-aminopropanesulfonic acid (TAPS)-EDTA), showed similar improvement, but the effect on the 1 x TBE buffer was more amplified. The difference in buffer properties, such as viscosity, conductance, and pH, had little effect on the separation. We attribute the improvements made in the separation to the ability of glycerol to induce a conformational change in DNA as demonstrated by dynamic light scattering results. The presence of glycerol can increase the electrostatic interactions between the phosphate groups, decrease the hydration sphere of the polynucleotides, and compete with water to form hydrogen bonds with the side group of bases. These interactions increase the DNA contour length and reduce the effective charge over weight ratio, which can explain the experimental data. The complex formed by boric acid and glycerol had a stronger effect on the DNA conformation change than glycerol itself. This enhancement was also observed in DNA sequencing analysis.

Fluorescence resonance energy transfer as a structural tool for nucleic acids

Fluorescence resonance energy transfer is a spectroscopic method that provides distance information on macromolecules in solution in the range 20-80 A. It is particularly suited to the analysis of the global structure of nucleic acids because the long-range distance information provides constraints when modelling these important structures. The application of fluorescence resonance energy transfer to nucleic acid structure has seen a resurgence of interest in the past decade, which continues to increase. An especially exciting development is the recent extension to single-molecule studies.

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.

Restrained torsional dynamics of nuclear DNA in living proliferative mammalian cells

Physical parameters, describing the state of chromatinized DNA in living mammalian cells, were revealed by in situ fluorescence dynamic properties of ethidium in its free and intercalated states. The lifetimes and anisotropy decays of this cationic chromophore were measured within the nuclear domain, by using the ultra-sensitive time-correlated single-photon counting technique, confocal microscopy, and ultra-low probe concentrations. We found that, in living cells: 1) free ethidium molecules equilibrate between extracellular milieu and nucleus, demonstrating that the cation is naturally transported into the nucleus; 2) the intercalation of ethidium into chromatinized DNA is strongly inhibited, with relaxation of the inhibition after mild (digitonin) cell treatment; 3) intercalation sites are likely to be located in chromatin DNA; and 4) the fluorescence anisotropy relaxation of intercalated molecules is very slow. The combination of fluorescence kinetic and fluorescence anisotropy dynamics indicates that the torsional dynamics of nuclear DNA is highly restrained in living cells.

Effect of temperature on DNA secondary structure in the absence and presence of 0.5 M tetramethylammonium chloride

Changes in the average secondary structures of three different linear DNAs over the premelting region from 5 to 60 degrees C were investigated by measuring their CD spectra and also their torsion elastic constants (<alpha>) by time-resolved fluorescence polarization anisotropy. For one of these DNAs, the Haell fragment of pBR322, the apparent diffusion coefficients [Dapp(k)] at small and large scattering vectors (k) were also measured by dynamic light scattering. With increasing temperature, all three DNAs exhibited typical premelting changes in their CD spectra, and these were accompanied by 1.4- to 1.7-fold decreases in <alpha>. Also for the 1876 base pair fragment, Dapp(k) at large scattering vectors, which is sensitive to the dynamic bending rigidity, decreased by 17%, even though there was no change at small scattering vectors, where Dapp(k) = D0 is the translational diffusion coefficient of the center-of-mass. These observations demonstrate conclusively that the premelting CD changes of these DNAs are associated with a significant change in average secondary structure and mechanical properties, though not in persistence length. In the presence of 0.5 M tetramethylammonium chloride (TMA-Cl) the premelting change in CD is largely suppressed, and the corresponding changes in <alpha> and Dapp(k) at large scattering vectors are substantially diminished. These observations suggest that TMA-Cl, which binds preferentially to A.T-rich regions and stabilizes those regions (relative to G.C-rich regions) against melting, effectively stabilizes the prevailing low-temperature secondary structure sufficiently that the DNA is effectively trapped in that state over the temperature range observed.

On the origin of the temperature dependence of the supercoiling free energy

Monte Carlo simulations using temperature-invariant torsional and bending rigidities fail to predict the rather steep decline of the experimental supercoiling free energy with increasing temperature, and consequently fail to predict the correct sign and magnitude of the supercoiling entropy. To illustrate this problem, values of the twist energy parameter (E(T)), which governs the supercoiling free energy, were simulated using temperature-invariant torsion and bending potentials and compared to experimental data on pBR322 over a range of temperatures. The slope, -dE(T)/dT, of the simulated values is also compared to the slope derived from previous calorimetric data. The possibility that the discrepancies arise from some hitherto undetected temperature dependence of the torsional rigidity was investigated. The torsion elastic constant of an 1876-bp restriction fragment of pBR322 was measured by time-resolved fluorescence polarization anisotropy of intercalated ethidium over the range 278-323 K, and found to decline substantially over that interval. Simulations of a 4349-bp model DNA were performed using these measured temperature-dependent torsional rigidities. The slope, -dE(T)/dT, of the simulated data agrees satisfactorily with the slope derived from previous calorimetric measurements, but still lies substantially below that of Duguet's data. Models that involve an equilibrium between different secondary structure states with different intrinsic twists and torsion constants provide the most likely explanation for the variation of the torsion constant with T and other pertinent observations.

Extension of torsionally stressed DNA by external force

Metropolis Monte Carlo simulation was used to study the elasticity of torsionally stressed double-helical DNA. Equilibrium distributions of DNA conformations for different values of linking deficit, external force, and ionic conditions were simulated using the discrete wormlike chain model. Ionic conditions were specified in terms of DNA effective diameter, i.e., hard-core radius of the model chain. The simulations show that entropic elasticity of the double helix depends on how much it is twisted. For low amounts of twisting (less than about one turn per twist persistence length) the force versus extension is nearly the same as in the completely torsionally relaxed case. For more twisting than this, the molecule starts to supercoil, and there is an increase in the force needed to realize a given extension. For sufficiently large amounts of twist, the entire chain is plectonemically supercoiled at low extensions; a finite force must be applied to obtain any extension at all in this regime. The simulation results agree well with the results of recent micromanipulation experiments.

Fluorescence anisotropy of DNA/DAPI complex: torsional dynamics and geometry of the complex

Fluorescence depolarization of synthetic polydeoxynucleotide/4'-6-diamidino-2-phenylindole dihydrochloride complexes has been investigated as a function of dye/polymer coverage. At low coverage, fluorescence depolarization is due to local torsional motions of the DNA segment where the dye resides. At relatively high coverage, fluorescence depolarization is dominated by energy transfer to other dye molecules along the DNA. The extent of the observed depolarization due to torsional motion depends on the angle the dye molecule forms with the DNA helical axis. A large torsional motion and a small angle produce the same depolarization as a small torsional motion and a large projection angle. Furthermore, the extent of transfer critically depends on the relative orientation of dye molecules along the DNA. The effect of multiple transfer is examined using a Monte Carlo approach. The measurement of depolarization with transfer, at high coverage, allows determination of the dye orientation about the DNA helical axis. The value of the torsional spring constant is then determined, at very low coverage, for few selected polydeoxynucleotides.

Conformation of short DNA fragments by modulated fluorescence polarization anisotropy

The technique of fluorescence polarization anisotropy (FPA) decay of intercalated ethidium has been used to study DNA conformation and dynamics, which are being recognized as primary determinants in transcription control and other cellular processes. Frequency modulated FPA when applied to two DNA molecules, a "straight" 50 base-pairs duplex fragment, and a bent fragment of similar length, has yielded different rotational diffusion coefficients for the two fragments. The data have been processed with an analytical model and with Brownian dynamics simulations, obtaining a good fit and a quantitative agreement between the two models. Both analyses have confirmed that one fragment can be described as a straight cylinder, while the other fragment is bent, with an angle estimated to be 45 degrees +/- 3 degrees. FPA has proved to be very powerful in determining simple conformations of short DNA duplexes and also particularly apt to probe the dynamical features of DNA fragments where conventional methods are either too cumbersome or fail to give quantitative results. In addition, the ligand no longer behaves ideally due to its complex structure and charge distribution. Thus for the protein the slope is no longer related simply to the net ligand charge, and the PB model gives a much larger slope than the LL model.(ABSTRACT TRUNCATED AT 400 WORDS)

Monte Carlo simulations of supercoiling free energies for unknotted and trefoil knotted DNAs

A new Monte Carlo (MC) algorithm is proposed for simulating inextensible circular chains with finite twisting and bending rigidity. This new algorithm samples the relevant Riemann volume elements in a uniform manner, when the constraining potential vanishes. Simulations are performed for filaments comprising 170 subunits, each containing approximately 28 bp, which corresponds to a DNA length of 4770 bp. The bending rigidity is chosen to yield a persistence length, P = 500 A, and the intersubunit potential is taken to be a hard-cylinder potential with diameter d = 50 A. This value of d yields the same second virial coefficient as the electrostatic potential obtained by numerical solution of the Poisson-Boltzmann equation for 150 mM salt. Simulations are performed for unknotted circles and also for trefoil knotted circles using two different values of the torsional rigidity, C = (2.0 and 3.0) x 10(-19) dyne cm2. In the case of unknotted circles, the simulated supercoiling free energy varies practically quadratically with linking difference delta l. The simulated twist energy parameter ET, its slope dET/dT, and the mean reduced writhe <w>/delta l for C = 3 x 10(-19) dyne cm2 all agree well with recent simulations for unknotted circles using the polygon-folding algorithm with identical P, d, and C. The simulated ET vs. delta l data for C = 2.0 x 10(-19) dyne cm2 agree rather well with recent experimental data for p30 delta DNA (4752 bp), for which the torsional rigidity, C = 2.07 x 10(-19) dyne cm2, was independently measured. The experimental data for p30 delta are enormously more likely to have arisen from C = 2.0 x 10(-19) than from C = 3.0 x 10(-19) dyne cm2. Serious problems with the reported experimental assessments of ET for pBR322 and their comparison with simulated data are noted. In the case of a trefoil knotted DNA, the simulated value, (ET)tre, exceeds that of the unknotted DNA, (ET)unk, by approximately equal to 1.40-fold at magnitude of delta l = 1.0, but declines to a plateau about 1.09-fold larger than (ET)unk when magnitude of delta l > or = 15. Although the predicted ratio, (ET)tre/(ET)unk approximately equal to 1.40, agrees fairly well with recent experimental measurements on a 5600-bp DNA, the individual measured ET values, like some of those reported for pBR322, are so large that they cannot be simulated using P = 500 A, d = 50 A, and any previous experimental estimate of C.

Effect of ethidium binding and superhelix density on the supercoiling free energy and torsion and bending constants of p30 delta DNA

Topoisomer distributions created by the action of topoisomerase I on p30 delta DNA in the presence of various concentrations of ethidium are measured and analyzed using recently developed theory to obtain the twist energy parameter (ET) that governs the free energy of supercoiling in each case. Competitive dialysis experiments to investigate the relative affinity of ethidium for linear and supercoiled DNAs at different binding ratios are assayed fluorometrically and the results are analyzed using related theory. The topoisomer distributions and fluorescence intensity ratios agree well with the theory, which is based on the assumption that the supercoiling free energy varies quadratically with the effective linking difference, regardless of ethidium binding or superhelix density. The topoisomer distribution experiments alone yield an average best-fit value, ET = 950 +/- 80, independent of ethidium binding ratio from r = 0 to 0.082, while the combined topoisomer distribution and ethidium binding experiments yield an average best-fit value, ET = 1030 +/- 90, which is essentially independent of ethidium binding ratio from r = 0 to 0.082 and superhelix density from sigma = 0 to (-)0.053. One may conclude that the supercoiling free-energy-varies quadratically with effective linking difference over the entire range of observed ethidium binding ratios and superhelix densities. The independently measured torsion constant (alpha) of p30 delta DNA is likewise essentially independent of superhelix density and ethidium binding ratio. The observed invariance of ET and alpha implies that the bending constant kappa beta is similarly invariant to superhelix density and ethidium binding ratio. The apparently ideal behavior displayed by p30 delta DNA is not exhibited by pBR322 DNA, which is discussed in the following companion paper.

A model for the binding of E. coli single-strand binding protein to supercoiled DNA

A model is proposed for the binding of E. coli single strand binding protein (SSB) to supercoiled DNA. The basic tetrameric binding units of SSB are assumed to bind in pairs to the complementary single strands of a locally melted region. The cooperativity of the binding includes contributions from both protein-protein and base-pair stacking interactions. Each bound SSB tetramer is assumed to unwind l = 34 bp, which implies an unwinding angle of 3.27 turns. The resulting loss of superhelical strain is the essential driving force for binding SSB to supercoiled DNAs. All molecular parameters entering into the theory are estimated from available data, except for the composite binding constant (Ka), which is adjusted to best-fit the theory to the fluorescence quenching (FQ) and diffusion coefficient (D0) data of Langowski et al. Very good fits are obtained with optimum values of Ka that are consistent with estimates from other data. This binding model predicts several noteworthy features. (1) SSB binds essentially always in a single contiguous stack on a supercoiled plasmid, and relative fluctuations in stack length are quite small, in agreement with results of electron microscopy studies. (2) The progressive loss of superhelical strain with increasing bound ligand decreases the affinity of the DNA for SSB. This anti-cooperativity offsets the cooperativity of the binding and causes apparent saturation of the binding at rather low binding ratios. Consequently, over the limited span of the measurements, the FQ data can also be satisfactorily fitted by a non-cooperative model comprising a small number of independent sites. (3) When SSB binds to a population of different topoisomers, the distribution of linking differences of the resulting complexes is extremely narrow. Thus, SSB acts to level any differences in superhelical strain in a population of topoisomers. Finally, the effects of restricting binding to a region comprising only part of the plasmid are assessed.

Effect of ethidium binding and superhelix density on the apparent supercoiling free energy and torsion constant of pBR322 DNA

The value of the twist energy parameter (ET) of pBR322 is determined near zero superhelix density from topoisomer distributions created under various conditions. The resulting value, ET = 1155 +/- 65, at 37 degrees C is essentially unaffected by adding 10 mM Mg2+, or by changing the kind of Topo I from chicken-red-cell to calf-thymus. This value significantly exceeds that (ET = 950 +/- 80) measured for p30 delta DNA under identical conditions by the same method in the preceding paper. Decreasing the temperature from 37 to 21 degrees C yields a slightly larger value, ET = 1340 +/- 130, but the statistical significance of the increase is marginal. Attempts to determine reliable ET values for pBR322 at higher superhelix densities by ethidium binding were frustrated by the fact that good fits of the equilibrium dialysis results could not be achieved using a single value of ET. Moreover, the curves of apparent ET versus binding ratio r vary considerably from one preparation to another, and for a given preparation vary with time after cell lysis up to about seven weeks, after which they settle in to nearly reproducible behavior. The apparent ET values obtained from competitive dialysis experiments are typically rather low (ET approximately 700) for small r and nearly native superhelix density, and rise up to 1300 to 1500 with increasing binding ratio (up to r = 0.055) and decreasing negative superhelix density.(ABSTRACT TRUNCATED AT 250 WORDS)

Circularization of small DNAs in the presence of ethidium: a theoretical analysis

A rigorous theory is developed for ethidium binding to linear and circular DNAs and for the ratios of topoisomers produced upon ligation of an equilibrium population of noncovalently closed circles in the presence of ethidium. Assuming an unwinding angle theta E = 26 degrees for intercalated ethidium, optimum values of the intrinsic binding constant, KE = 7.16 x 10(4) M-1, the intrinsic twist, l0 = 23.746 turns, and twist energy parameter, Et = 5250, are obtained by fitting the present theory to the data of Shore and Baldwin [(1993) Journal of Molecular Biology, Vol. 170, pp. 983-1007] for a 247 base pair DNA. A very good fit is achieved with these optimum values, but a poor fit results when the parameters estimated by Shore and Baldwin are employed in the same theory. Three assumptions employed in the analysis of Shore and Baldwin are found to be not strictly valid. Adoption of the present substantially larger Et value as representative of their short DNAs would allow the Et vs N data of Shore and Baldwin to conform to the shape predicted by Shimada and Yamakawa [(1985) Journal of Molecular Biology, Vol. 184, pp. 319-329] and Frank-Kamenetskii et al. [(1985) Journal of Biomolecular Structure and Dynamics, Vol. 2, pp. 1005-1012], and would imply that all of their DNAs exist in a substantially stiffer than normal state.

The amplitude of local angular motion of purines in DNA in solution

Nuclear magnetic resonance and optical experiments are combined to determine the rms amplitude of local angular motion of purines in DNA in solution. A 12 base-pair duplex DNA with the sequence d(CGCGAATTCGCG)2 is deuterated at the H8 positions of adenine and guanine by exchange with solvent at 55 degrees C. The deuterium nmr spectrum of this DNA is measured at 30 mg/mL at 30 degrees C in an 11.76 Tesla magnet (76.75 MHz). The time-resolved fluorescence polarization anisotropies (FPA) of this same sample and also a greatly diluted sample (0.215 mg/mL) were measured after addition of ethidium. FPA measurements of the dilute sample yield the hydrodynamic radius, RH = 9.94 +/- 0.2 A, while those at the nmr concentration are employed to characterize the collective motions in terms of either an enhanced viscosity or dimer formation. The rms amplitude of local angular motion was determined by analyzing the 2H-nmr spectrum, in particular the line width, using recently developed theory for the transverse relaxation rate (RQ2) together with essential information about the collective motions from these and other optical studies. When the principal-axis frame of the electric field gradient tensor is assumed to undergo overdamped libration around each of its three body-fixed axes in an isotropic deflection potential, then the rms amplitude of local angular motion around any single axis is found to lie in the range 10 degrees-11 degrees, provided the high DNA concentration acts to enhance the viscosity, and is about 9 degrees-11 degrees, if it acts to produce end-to-end dimers. The proton nmr relaxation data of Eimer et al. are reanalyzed and shown to yield an rms amplitude of angular motion of the cytosine H5-H6 internuclear vector of 9 degrees-10 degrees, depending upon its orientation with respect to the helix axis. In all of these analyses, full account is taken of the collective twisting and bending deformations, which have a small but significant effect on the results. It is shown that the rms amplitudes of local angular motion do not depend strongly on the model (potential), provided that isotropic rotation around the same number of axes is allowed and that one compares rms angles of the same dimensionality. The rms amplitudes of local angular motion in solution are comparable to those observed for the same sequence at low levels of hydration in the solid state.

Modulation of DNA supercoiling by interaction with netropsin and other minor groove binders

The assay of DNA unwinding by ethidium, followed by sedimentation velocity techniques, was applied to complexes of supercoiled plasmid DNA with different non-intercalating drugs which strongly and sequence-specifically bind to DNA. Compared with the behaviour of naked DNA, most of the complexes exhibit an increase in the critical EB/nucleotide binding ratio associated with the principal minimum in the sedimentation profile. Using netropsin (Nt) as the paradigm of the minor groove binders investigated, the drug-induced alterations in various structural parameters of both the relaxed and supercoiled form of DNA are described. Whereas winding number, helical repeat (both being defined with reference to a surface normal), and linking number of the superhelical DNA remain constant in our experiments, its twist number, surface twist, number of superhelical turns as well as the absolute values of linking number difference, superhelix density, and writhing number increase on binding of Nt. Correspondingly, compared with the naked relaxed DNA a higher linking number (or twist number, or winding number), a higher average duplex winding angle and a lower helical repeat have to be assigned to the relaxed Nt-DNA complex. The various minor groove binders investigated were found to differ considerably in their efficiency to alter the structure of supercoiled DNA.

Dynamics and structures of DNA: long-range effects of a 16 base-pair (CG)8 sequence on secondary structure

The effects of inserting 16 base pair (bp) of alternating CG [(CG)8] near the middle of a much longer restriction fragment (1097 bp) are investigated by measuring various properties that are sensitive to secondary and tertiary structure. Results for this fragment are compared with those for a control fragment (1089 bp) with the identical sequence except at the insert. Another fragment (1382 bp), which contains a 296-bp extension at the 5'-end of the 1089-bp control fragment, is also used as a secondary control in some experiments. When the 1097-bp (CG)8 insert fragment is compared with the control fragments in 0.1 M NaCl buffer, the (CG)8 insert is found to induce disproportionately large relative changes in the molar ellipticity at 273 nm ([theta]273), the torsion constant (alpha) measured by fluorescence polarization anisotropy, the optical melting profile, and the susceptibility to S1 nuclease. Estimates of the minimum distance over which the (CG)8 insert alters the secondary structure range from 330 to 550 bp. With increasing NaCl concentration, the 1097-bp insert fragment undergoes a structural transition between 2.0 and 2.5 M that is manifested in the apparent diffusion coefficient (Dplat) from dynamic light scattering at large scattering vector. This transition, which is not exhibited by the control DNAs, is presumed to involve formation of Z-helix at the insert. However, the observed decrease in (Dplat) is attributed to an increase in bending rigidity, which perforce must be globally distributed far beyond the (CG)8 insert per se. In 4.25 M NaCl (but not in 0.1 M NaCl), the addition of 1 ethidium dye per 300 bp induces an extensive structural transition in the 1097 bp (CG)8 insert fragment. This transition, which also is not exhibited by the control DNAs, significantly decreases the bending rigidity, doubles [theta]273, and takes place on a time scale of a few days. Removal of ethidium and salt by dialysis vs 0.1 M NaCl buffer restores the original properties of the 1097-bp (CG)8 insert fragment. The present results are consistent with a (fluctuating, long-range) description of the secondary structure in which a given short sequence transiently fluctuates among two or more distinct secondary structures that extend over much larger domains of variable position and size, and whose relative stabilities depend on distant as well as close-lying base pairs.

Twist constraints on linker DNA in the 30-nm chromatin fiber: implications for nucleosome phasing

Previous work has shown that nucleosome repeat lengths, and hence linker DNA lengths, are preferentially quantized to a set of values differing by integral multiples of the helical twist of DNA. An explanation was proposed in which this preferential quantitation is due to twist constraints on linker DNA arising from nucleosome-nucleosome interactions in folded chromatin. Here we report the results of a study, using ethidium intercalation, designed to test whether twist constraints do indeed exist. Electron microscopy reveals that ethidium intercalation causes decondensation of dinucleosomes. Direct measurement of the free energy of intercalation by fluorescence spectroscopy reveals competition between chromatin folding and ethidium intercalation. Results from other laboratories establish that these effects of ethidium are due to ethidium-induced changes in the twist of linker DNA, and not to a variety of other possible effects. We conclude that twist constraints on linker DNA do exist. These may explain the observation of preferentially quantized linker DNA lengths. Implications of these results for mechanisms of nucleosome phasing and the mechanisms of drug action are discussed.