Effect of polyethylene glycol on the supercoiling free energy of DNA.

The supercoiling free energy of pUC19 DNA [2686 base pairs (bp)] was measured in various concentrations of PEG 8000 (polyethylene glycol; molecular weight 8000) by the topoisomer distribution method. The effective twist energy parameter (E(T)) that governs the supercoiling free energy declined linearly by 1.9-fold with increasing w/v % PEG from 0 to 7.5%, which lies below the threshold for intermolecular condensation. In principle, PEG could affect E(T) either via an osmotic exclusion mechanism or by altering the torsion elastic constant, bending rigidity, or self-repulsions of the DNA. Possible alterations of the DNA secondary structure and torsion elastic constant were assessed by CD spectroscopy and time-resolved fluorescence polarization anisotropy of intercalated ethidium. Up to 7.5% PEG, the secondary structure of the DNA remained largely unaltered, as evidenced by (1) the absence of any significant change in the CD spectrum, (2) an extremely small relative decrease (-0.0013) in intrinsic twist, and (3) a negligibly small change in the torsion elastic constant. The observed reduction in E(T) cannot be ascribed primarily to a decrease in torsion elastic constant, and most likely does not stem from a decrease in bending rigidity either. The decrease in medium dielectric constant due to PEG should increase the self-repulsions, and thereby increase E(T), which is opposite to the observed trend. Instead, the observed decline in E(T) is attributed to an osmotic exclusion mechanism. The change in molar volume excluded to the PEG (Delta V(ex)), when the linking difference converts from Delta l = 0 to Delta l = +/-1, was determined from the observed E(T) value and PEG osmotic pressure at each concentration. The experimental Delta V(ex) values agree well with theoretical estimates reckoned for a simple osmotic exclusion model, in which PEG is excluded by hard-core interactions from a concentric cylindrical volume around every duplex segment. The difference in volume excluded to PEG between the Delta l = 0 and the Delta l = +/-1 topoisomers is attributed entirely to the approximately 0.7 additional writhe "crossing" of two duplex strands at roughly 90 degrees, which is known to occur in the latter species. When the separation between the duplex centers at the "crossing" was adjusted so that the theoretical estimate of Delta V(ex) matched the experimental value at each PEG concentration, a value near 5.7 nm was obtained in each case. The invariance and plausible magnitude of this mean separation at the crossing provide strong support for this simple osmotic exclusion model. An alternative model, in which the PEG is excluded from the entire coil envelope of the DNA out to its radius of gyration, perhaps because it decreases the local dielectric constant, was also considered. The estimated difference in excluded volume in that case exceeds the experimental value by a factor of nearly 10(4), and could be ruled out on that basis.

The distribution of end-to-end distances of the weakly bending rod model.

The weakly bending rod model is an approximation to a worm-like chain in the limit where the ratio L(0)/P of the contour length L(0) to the persistence length P is not too large. The range of validity of the weakly bending rod model is investigated by deriving analytical expressions for its distribution of end-to-end distances P(L) and its moments and numerically comparing the results with corresponding values for the worm-like chain model. No general, closed form analytical expression for either P(L) or the average length of a worm-like chain exists, so those quantities are obtained by Monte Carlo simulations. Exact analytical expressions for and for the worm-like chain are employed in the comparison of the computed moments. Moments calculated for the approximate distributions of Daniels and Yamakawa and Stockmayer are also compared with the others. In addition, P(L) and its moments for the alternative model of Winkler et al. are compared with the others. The weakly bending rod model gives a reasonably good account of both P(L) and , m = 1,2,4, over the range L(0)/P /= 1.0. In contrast, the alternative model of Winkler et al. yields rather poor results in the rodlike domain.

Dynamic bending rigidity of a 200-bp DNA in 4 mM ionic strength: a transient polarization grating study.

DNA may exhibit three different kinds of bends: 1) permanent bends; 2) slowly relaxing bends due to fluctuations in a prevailing equilibrium between differently curved secondary conformations; and 3) rapidly relaxing dynamic bends within a single potential-of-mean-force basin. The dynamic bending rigidity (kappa(d)), or equivalently the dynamic persistence length, P(d) = kappa(d)/k(B)T, governs the rapidly relaxing bends, which are responsible for the flexural dynamics of DNA on a short time scale, t < or = 10(-5) s. However, all three kinds of bends contribute to the total equilibrium persistence length, P(tot), according to 1/P(tot) congruent with 1/P(pb) + 1/P(sr) + 1/P(d), where P(pb) is the contribution of the permanent bends and P(sr) is the contribution of the slowly relaxing bends. Both P(d) and P(tot) are determined for the same 200-bp DNA in 4 mM ionic strength by measuring its optical anisotropy, r(t), from 0 to 10 micros. Time-resolved fluorescence polarization anisotropy (FPA) measurements yield r(t) for DNA/ethidium complexes (1 dye/200 bp) from 0 to 120 ns. A new transient polarization grating (TPG) experiment provides r(t) for DNA/methylene blue complexes (1 dye/100 bp) over a much longer time span, from 20 ns to 10 micros. Accurate data in the very tail of the decay enable a model-independent determination of the relaxation time (tau(R)) of the end-over-end tumbling motion, from which P(tot) = 500 A is estimated. The FPA data are used to obtain the best-fit pairs of P(d) and torsion elastic constant (alpha) values that fit those data equally well, and which are used to eliminate alpha as an independent variable. When the relevant theory is fitted to the entire TPG signal (S(t)), the end-over-end rotational diffusion coefficient is fixed at its measured value and alpha is eliminated in favor of P(d). Neither a true minimum in chi-squared nor a satisfactory fit could be obtained for P(d) anywhere in the range 500-5000 A, unless an adjustable amplitude of azimuthal wobble of the methylene blue was admitted. In that case, a well-defined global minimum and a reasonably good fit emerged at P(d) = 2000 A and (1/2) = 25 degrees. The discrimination against P(d) values <1600 A is very great. By combining the values, P(tot) = 500 A and P(d) = 2000 A with a literature estimate, P(pb) = 1370 A, a value P(sr) = 1300 A is estimated for the contribution of slowly relaxing bends. This value is analyzed in terms of a simple model in which the DNA is divided up into domains containing m bp, each of which experiences an all-or-none equilibrium between a straight and a uniformly curved conformation. With an appropriate estimate of the average bend angle per basepair of the curved conformation, a lower bound estimate, m = 55 bp, is obtained for the domain size of the coherently bent state. Previous measurements suggest that this coherent bend is not directional, or phase-locked, to the azimuthal orientation of the filament.

Manifestations of slow site exchange processes in solution NMR: a continuous Gaussian exchange model.

The effects of site exchange due to slow conformational changes in rapidly rotating molecules in solution are examined in detail. Significant gaps in the currently available theory are filled. The effects of site exchange on the lineshape, decay of a simple spin-echo, decay of the even echoes in a Carr-Purcell-Meiboom-Gill (CPMG) pulse-sequence, and decay of the transverse magnetization in a resonant spin-locking field are investigated. Both trajectory and stochastic operator approaches are formulated and shown to be completely equivalent whenever the dynamics of population transfers among the inequivalent sites is governed by either a stationary or a nonstationary Markov process. A nonstationary Markov process may result from Brownian dynamics (a stationary Markov process) in a larger conformational space that contains the subspace of inequivalent sites. A continuous Gaussian exchange model is formulated in which a nucleus undergoes continuous one-dimensional motion in a harmonic potential well that is located in a linear chemical shift gradient. The effects of this Gaussian exchange model on the lineshape, simple spin-echo decay, and decay of the even echoes of a CPMG pulse train are treated rigorously via the trajectory approach. Compact analytical expressions are obtained for the relevant correlation functions in each case. The relevant decays are found to be exponential in the very short time and long time limits, which are not necessarily experimentally significant in any given case. In the fast exchange limit the relevant decays are exponential at all times, and explicit formulas are given for their decay rates. In the long time limit, all discrete multisite models with the same intrinsic Ro2 at every site are shown to be completely equivalent to a continuous Gaussian model with appropriate relaxation time and variance of the Larmor frequency. The effects of this Gaussian exchange model on the decay of the transverse magnetization in a resonant spin-locking field are treated heuristically by a trajectory approach. The intrinsic contribution (Ro1rho) of rapid rotations and dipole-dipole interactions to relax the transverse magnetizations of two nuclei of the same kind in the presence of a (nearly) resonant spin-locking field is also derived and found to be practically the same as the intrinsic contribution, Ro2, of those same rotations to the simple and CPMG spin-echo decay rates and linewidth. Literature data for the linewidth, decay rate of the CPMG even spin-echoes, and R(1rho) decay rate for the A9-H2 protons of adenines at the central TpA step in the sequence, 5'-GCAGGTTTAAACCTCG-3', are analyzed using the Gaussian exchange model to assess the time-scale and variance of the site exchange process as well as the intrinsic Ro2 rate. Although a single Gaussian exchange process with appropriate parameters can fit these three A9-H2 data rather well, this particular "solution" cannot be reconciled with NMR relaxation data on other protons in the same DNA molecule. Rather good agreement with all of the observations is obtained by using a model of two concurrent Gaussian exchange processes, whose relaxation times, tau = 7 and 460 micros, differ in time-scale by a factor of 65. The insensitivity of R1rho in the presence of a fast site exchange process to much slower concurrent site exchange processes is explicitly demonstrated. Protocols for detecting and characterizing a second slow site exchange process are suggested.

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

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.

Effects of Na+ and Mg2+ on the structures of supercoiled DNAs: comparison of simulations with experiments.

Recent cryo-electron microscopy (cryo-EM) results suggest that sufficient NaCl concentration (> or approximately 0.1 M) and superhelix density (> or approximately-0.05) cause circular DNAs to adopt highly extended, tightly interwound configurations, in which the strands are laterally contiguous along almost their entire length. Millimolar levels of MgCl2 reportedly act synergistically with NaCl to produce similar conformations. However, Monte Carlo simulations with purely repulsive interduplex forces failed to reproduce such structures. In the present work, solution measurements of particular physical properties were performed both to characterize the effects of Na+ and Mg2+ on DNA structure and to provide quantitative tests of Monte Carlo simulations of circular DNAs. Supercoiled p30 delta DNAs in 10 mM Tris plus 0, 0.122, and 0.1 M NaCl, and 0.1 M NaCl plus 4 mM Mg2+ were examined by static and dynamic light scattering (LS and DLS), time-resolved fluorescence polarization anisotropy (FPA) of intercalated ethidium, and circular dichroism (CD) spectroscopy. Upon addition of 0.122 M NaCl, the radius of gyration (Rg) decreased substantially, which indicates that p30 delta adopts a more compact structure. This contradicts the cryo-EM studies, where molecular extension and Rg both increase upon adding 0.1 M NaCl. In 0.1 M NaCl, the torsion constant measured by FPA is practically invariant to superhelix density, and the plateau diffusion coefficient at large scattering vector (Dplat) is likewise nearly the same at both relaxed and native superhelix densities. Such invariance is difficult to reconcile with any transition from relaxed circles to tightly interwound structures with laterally contiguous strands. Metropolis Monte Carlo simulations were performed to generate canonically distributed sets of structures, from which average Do values and scattered intensity ratios, [symbol: see text]I (zero) [symbol: see text]/[symbol: see text] l(k) [symbol: see text], were calculated. Agreement between simulations and experiments in regard to [symbol: see text] I(O) [symbol: see text] /[symbol: see text] I(k) [symbol: see text], D(zero) and the supercoiling free energy, delta Gsc (delta l), is remarkably good for the most extensively studied p30 delta samples. The simulated structures exhibit no sign of very tight interwinding with extensive lateral contacts, but instead exhibit most probable superhelix diameters of 85 to 90 A. When 4 mM Mg2+ was added to native supercoiled p30 delta in 0.1 M NaCl, Rg decreased, D(zero) increased, and the longest internal relaxation rate (1/tau 2(zero)) increased, all of which indicate a further overall contraction of the molecular envelope. The torsion constant exhibited a slight increase that is hardly statistically significant. In this case, agreement between the simulations and experiments was only semi-quantitative for most samples investigated, although the predicted contraction was exhibited by all five samples of p30 delta and one of pBR322 DNA. The simulated structures in 0.1 M NaCl plus 4 mM Mg2+ again showed no sign of extensive lateral contacts. A plausible explanation is proposed for the highly extended, tightly interwound structures seen in cryo-EM, and explicitly tested by Monte Carlo simulations of a 1000 bp circular DNA at +25 and -50 degrees C. Structures identical to those seen in cryo-EM are in fact the equilibrium structures in the simulations at -50 degrees C, and the estimated time for equilibration (2.3 x 10(-6) second) is much smaller than the estimated time for vitrification (1 x 10(-4) second).

Effect of bending strain on the torsion elastic constant of DNA.

The torsion constants of both circular and linear forms of the same 181 bp DNA were investigated by time-resolved fluorescence polarization anisotropy (FPA) of intercalated ethidium. The ratio of intrinsic ethidium binding constants of the circular and linear species was determined from the relative fluorescence intensities of intercalated and non-intercalated dye in each case. Possible changes in secondary structure were also probed by circular dichroism (CD) spectroscopy. Upon circularization, the torsion constant increased by a factor of 1.42, the intrinsic binding constant for ethidium increased by about fourfold, and the CD spectrum underwent a significant change. These effects are attributed to an altered secondary structure induced by the bending strain. Quantitative agreement between torsion constants obtained from the present FPA studies and previous topoisomer distribution measurements on circular DNAs containing 205 to 217 bp removes a long-standing apparent discrepancy between those two methods. After storage at 4 degrees C for eight months, the torsion constant of the circular DNA increased by about 1.25-fold, whereas that of the linear DNA remained unchanged. For these aged circles, both the torsion constant and intrinsic binding constant ratio lie close to the corresponding values obtained previously for a 247 bp DNA by analyzing topoisomer distributions created in the presence of various amounts of ethidium. The available evidence strongly implies that torsion constants measured for small circular DNAs with less than 250 bp are specific to the altered secondary structure(s) therein, and are not applicable to linear and much larger circular DNAs with lower mean bending strains.

Position-dependent internal motions and effective correlation times for magnetization transfer in DNA.

An effective correlation time that accounts for position dependence of the combined local angular motions and collective twisting and bending deformations, as well as for anisotropic uniform rotation, is defined in terms of the magnetization-transfer rate and expressed in terms of molecular parameters. Application to measured magnetization-transfer rates from H6 to H5 of cytosine in cases where all other relevant data, including the uniform rotational diffusion coefficients, are known, suggests that the amplitude of local angular motion, as well as that due to collective deformations, is significantly greater for a penultimate base pair than for base pairs near the center of the molecule, and that such amplitudes might be approximately transferable from one molecule to another. Protocols are suggested for using estimated ratios of effective correlation times in the initial calibrations of internuclear distances and in the subsequent structure-refinement process.

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.

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)

A test of the model-free formulas. Effects of anisotropic rotational diffusion and dimerization.

The effects of anisotropy of rotational diffusion and extent of dimerization on the performance of the simple and extended model-free formulas are investigated. Numerically exact 15N NMR R1, R2, and NOE data are simulated for cylindrically symmetric species and also for mixtures of spherical monomers and dimers with the same internal dynamics. The relevant internuclear vectors are assumed to move in isotropic deflection potentials fixed at different orientations in the molecule and to exhibit a single relaxation time in their internal correlation functions. The simple model-free formula is fitted to these simulated data in order to obtain the best-fit order parameter and internal relaxation time for each nucleus. Fitting is accomplished by the standard data-analysis protocol, in which a single common global correlation time is adjusted, and also by an alternative protocol, in which the global correlation time is adjusted separately for each nucleus. The extended model-free formula is likewise fitted to these same data. With noise-free data, the simple model-free formula and standard protocol yield remarkably good best-fit internal motion parameters up to moderate anisotropies (r = D parallel/D perpendicular = 2.0), but some or many of the NMR relaxation data are not fitted well even for quite modest anisotropies (r = 1.3). The extended model-free formula yields an improved fit to the NMR data, but predicts substantial amplitudes of nonexistent slow internal motions (tau approximately greater than 0.2 ns) for many of the nuclei for all r > or = 1.3. The simple model-free formula with the alternate protocol yields even better internal motion parameters than the standard protocol and also an excellent fit to the NMR relaxation data. The best-fit global correlation time for each nucleus corresponds very closely to the theoretical correlation time defined herein. Knowledge of these times would allow one not only to estimate the anisotropy of diffusion but also in favorable cases to infer the existence of slow internal motions. Inclusion of typical statistical errors in R1, R2, and NOE, or modest exchange contributions to R2, seriously degrades the performance of the simple model-free formula with either protocol, especially in regard to the internal relaxation times, which can exhibit very large deviations from their input value even for spherical diffusors. When the fraction of monomers existing as dimers lies in the range 0.1 < or = fd < or = 0.8, none of the three model-free approaches tested yields reliable internal motion parameters.

Effects of different cations on the hydrodynamic radius of DNA.

The effects of different cations on the hydrodynamic radius (RH) of a 48-bp synthetic DNA are measured by time-resolved fluorescence polarization anisotropy of intercalated ethidium. Relative statistical errors in RH are only approximately 1%. With increasing cation concentration, Na+ causes a small decrease in RH, Cs+ causes a somewhat larger decrease by up to 0.5 A at 100 mM, and (CH3CH2)4N+ causes an increase in RH by approximately 0.5 A at 100 mM. The qualitatively different effects of these monovalent cations indicates that the changes in RH with cation concentration do not arise primarily from electrolyte friction. Divalent cations cause much larger increases in RH with increasing cation concentration. Mg2+ causes an increase in RH by up to 1.0 A at 24.4 mM, and Mn2+ causes an increase in RH by up to 1.6 A at 24.4 mM. These effects are independent of DNA concentration. There is some positive correlation between the order of effects of the different cations on RH and the order of their effects on interhelical hydration forces. It is suggested that these different ions affect RH either by altering the hydration layer or possibly by some effect on DNA structure, such as stabilizing bends.

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.

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.

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

The torsion elastic constants (alpha) of linear pBR322 (4363 bp) and pUC8 (2717 bp) DNAs and supercoiled pBR322 and pJMSII (4375 bp) DNAs are measured in 0.1 M NaCl as a function of added ethidium/base-pair (EB/BP) ratio by studying the fluorescence polarization anisotropy (FPA) of the intercalated ethidium. The time-resolved FPA is measured by using a picosecond dye laser for excitation and time-correlated single photon counting detection. Previously developed theory for the emission anisotropy is generalized to incorporate rotations of the transition dipole due to excitation transfer. The excitation transfers are simulated by a Monte Carlo procedure (Genest et al., Biophys. Chem. 1 (1974) 266-278) and the consequent rotations of the transition dipole are superposed on the Brownian rotations. After accounting for excitation transfer, the torsion constants of the linear DNAs are found to be essentially independent of intercalated ethidium up to a binding ratio r = 0.10 dye/bp. Dynamic light scattering measurements on linear pUC8 DNA confirm that the torsion constant is independent of binding ratio up to r = 0.20 dye/bp. If alpha d denotes the torsion constant between ethidium and a base-pair, and alpha 0 that between two base-pairs, then our data imply that alpha d/alpha 0 lies in the range 0.65 to 1.64 with a most probable value of 1.0. The torsion constants of supercoiled DNAs decrease substantially with increasing binding ratio even after accounting for excitation transfer. At the binding ratio r* = 0.064, where the superhelix density vanishes and superhelical strain is completely relaxed, the torsion constant of the supercoiled pBR322 DNA/dye complex lies below that of the corresponding linear DNA/dye complex by about 30%. This contradicts the conventional view according to which linear, nicked circular, and supercoiled DNA/dye complexes with r = r* should coexist with the same concentration of free dye, display the same distribution of bound dye, and exhibit identical secondary structures, twisting and bending rigidities, and FPA dynamics. These and other observations imply the existence of metastable secondary structure in freshly relaxed supercoiled DNAs. A tentative explanation is presented for these and other unexpected observations on supercoiled DNAs.

Evidence for allosteric transitions in secondary structure induced by superhelical stress.

Previous studies suggest that the global secondary structures of native supercoiled and equilibrium linear DNAs may differ somewhat. Recent evidence also indicates that metastable secondary structure commonly persists following complete relaxation of the superhelical stress by intercalating dyes or by the action of topoisomerase I. In this work, the torsion constants (alpha) of pBR322, pUC8 and M13mp7 (replicative form) DNAs are determined by time-resolved fluorescence polarization anisotropy at various times subsequent to linearization. In all three cases, the torsion constants are relatively low immediately after linearization, and evolve for eight to ten weeks before reaching their apparent equilibrium values. It is shown in detail how the persistence of metastable secondary structure, subsequent to relaxation of superhelical stress, necessarily implies that one or more transitions in equilibrium secondary structure are induced as the superhelix density is varied from zero to native, or vice versa. Samples of pUC8 dimer (5434 base-pairs) with different superhelix densities are prepared by the action of topoisomerase I in the presence of various amounts of ethidium. Their median linking number differences are determined by standard band counting methods. The translational diffusion coefficient (Do) and the plateau diffusion coefficient (Dplat) characterizing internal motions over short distances (225 A) are determined by dynamic light-scattering. The torsion constant (alpha) between base-pairs and the circular dichroism spectrum are also measured for each sample. Curves of Dplat, Do, alpha and molar ellipticity ([theta]) (at the minimum near 250 nm) versus superhelix density (sigma) are constructed. The curve of Do versus sigma is very similar to that for sedimentation coefficient versus sigma for simian virus 40 (SV40) and polyoma DNAs. The curves of Dplat, Do, alpha and [theta] versus sigma show that, with increasing negative superhelix density, a structural transition occurs near sigma = -0.020 to an intermediate state with low torsion constant, and a second structural transition occurs near sigma = -0.035 to a state that exhibits more normal properties by sigma = -0.048. These data are consistent with the hypothesis that supercoiling induces two successive allosteric transitions to alternative global secondary structures. The data are much less consistent with the hypothesis that supercoiling induces some radical secondary structure at one or a few sites of small extent at sigma = -0.020, and at other sites at sigma = -0.035, or with hypotheses based on changes in tertiary structure alone.(ABSTRACT TRUNCATED AT 400 WORDS)

Dependence of the torsional rigidity of DNA on base composition.

The Escherichia coli phage 434 repressor binds as a dimer to the operator of the DNA helix. Although the centre of the operator is not in contact with protein, the repressor binding affinity can be reduced at least 50-fold by changing the sequence there: operators with A.T base pairs near their centre bind the repressor more strongly than do operators with G.C base pairs at the same positions. To explain these observations, it has been proposed that the base composition at the centre of the operator affects the affinity of the operator for repressor by altering the ease with which operator DNA can undergo the torsional deformation necessary for complex formation. In this model, the variation in binding affinity would require the torsion constant to have specific values and to change in a sequence-dependent manner. We have now measured torsion constants for DNAs with widely different base compositions. Our results indicate that the torsion constants depend only slightly on the overall composition, and firmly delimit the range of values for each. Even the upper-limit values are much too small to account for the observed changes in affinity of the 434 repressor. These results rule out simple models that rely on substantial generic differences in torsion constant between A.T-rich sequences and G.C-rich sequences, although they do not rule out the possibility of particular sequences having abnormal torsion constants.

Interaction of chloroquine with linear and supercoiled DNAs. Effect on the torsional dynamics, rigidity, and twist energy parameter.

The magnitude and uniformity of the torsion elastic constant (alpha) of linear pBR322 DNA and supercoiled pBR322 DNAs with high-twist (sigma = -0.083) and normal-twist (sigma = -0.48) are measured in 0.1 M NaCl as a function of added chloroquine/base-pair ratio (chl/bp) by studying the fluorescence polarization anisotrophy (FPA) of intercalated ethidium dye. The time-resolved FPA is measured by using a picosecond dye laser for excitation and time-correlated single-photon counting detection. A general theory is developed for the binding of ligands that unwind superhelical DNAs, and the simultaneous binding of two different intercalators is treated in detail. The equilibrium constant (K) for binding chloroquine to linear pBR322 DNA and the number (r) of bound chloroquines per base pair are determined from the relative amplitude ratio of the slow (normally intercalated) and fast (free) components in the decay of the (probe) ethidium fluorescence intensity as a function of chl/bp. For chloroquine binding to supercoiled pBR322 DNAs, the intrinsic binding constant is assumed to be the same as for the linear DNA, but the twist energy parameter ET (N times the free energy to change the linking number from 0 to 1 in units of kBT) is regarded as adjustable. Using the best-fit ET, the binding ratios r are calculated for each chl/bp ratio. Twist energy parameters are also determined for ethidium binding to these supercoiled DNAs by competitive dialysis. For chloroquine binding, we obtain ET = 360 and 460 respectively for the normal-twist and high-twist supercoiled DNAs. For ethidium binding the corresponding values are ET = 280 +/- 70 and 347 +/- 50. Like other dye-binding values, these are substantially lower than those obtained by ligation methods. In the absence of chloroquine, the torsion constants of all three DNAs are virtually identical, alpha = (5.0 +/- 0.4) x 10(-12) dyn.cm. For linear pBR322 DNA, the magnitude and uniformity of alpha remain unaltered by intercalated chloroquine up to r = 0.19. This finding argues that the FPA is not significantly relaxed by diffusion of any kinks or solitons. If alpha d denotes the torsion constant between a dye and a base pair and alpha 0 that between two base pairs, then our data imply that alpha d/alpha 0 lies in the range 0.65-1.64, with a most probable value of 1.0.(ABSTRACT TRUNCATED AT 400 WORDS)