Theoretical analysis of heteropolymeric transitions in superhelical DNA molecules of specified sequence

DNA superhelicity

Closing each strand of a DNA duplex upon itself fixes its linking number L. This topological condition couples together the secondary and tertiary structures of the resulting ccDNA topoisomer, a constraint that is not present in otherwise identical nicked or linear DNAs. Fixing L has a range of structural, energetic and functional consequences. Here we consider how L having different integer values (that is, different superhelicities) affects ccDNA molecules. The approaches used are primarily theoretical, and are developed from a historical perspective. In brief, processes that either relax or increase superhelicity, or repartition what is there, may either release or require free energy. The energies involved can be substantial, sufficient to influence many events, directly or indirectly. Here two examples are developed. The changes of unconstrained superhelicity that occur during nucleosome attachment and release are examined. And a simple theoretical model of superhelically driven DNA structural transitions is described that calculates equilibrium distributions for populations of identical topoisomers. This model is used to examine how these distributions change with superhelicity and other factors, and applied to analyze several situations of biological interest.

Topology and kinetics of R-loop formation

R-loops are structures containing an RNA-DNA duplex and an unpaired DNA strand. They can be formed upon "invasion" of an RNA strand into a DNA duplex, during which the RNA displaces the homologous DNA strand and binds the complementary strand. R-loops have many significant beneficial or deleterious biological effects, so it is important to understand the mechanisms for their generation and processing. We propose a model for co-transcriptional R-loop formation, in which their generation requires passage of the nascent RNA "tail" through the gap between the separated DNA strands. This passage becomes increasingly difficult with lengthening of the RNA tail. The length of the tail increases upon increasing distance between the transcription start site and the site of R-loop initiation. This causes reduced yields of R-loops with greater distance from the transcription start site. However, alternative pathways for R-loop formation are possible, involving either transient disruption of the transcription complex or the hypothetical formation of a triple-stranded structure, as a "collapsed R-loop." These alternative pathways could account for the fact that in many systems R-loops are observed very far from the transcription start site. Our model is consistent with experimental data and makes general predictions about the kinetics of R-loop formation.

DNA sequence encodes the position of DNA supercoils

The three-dimensional organization of DNA is increasingly understood to play a decisive role in vital cellular processes. Many studies focus on the role of DNA-packaging proteins, crowding, and confinement in arranging chromatin, but structural information might also be directly encoded in bare DNA itself. Here, we visualize plectonemes (extended intertwined DNA structures formed upon supercoiling) on individual DNA molecules. Remarkably, our experiments show that the DNA sequence directly encodes the structure of supercoiled DNA by pinning plectonemes at specific sequences. We develop a physical model that predicts that sequence-dependent intrinsic curvature is the key determinant of pinning strength and demonstrate this simple model provides very good agreement with the data. Analysis of several prokaryotic genomes indicates that plectonemes localize directly upstream of promoters, which we experimentally confirm for selected promotor sequences. Our findings reveal a hidden code in the genome that helps to spatially organize the chromosomal DNA.

Comparative Analysis of Mycoplasma gallisepticum vlhA Promoters

Mycoplasma gallisepticum is an intracellular parasite affecting respiratory tract of poultry that belongs to class Mollicutes. M. gallisepticum features numerous variable lipoprotein hemagglutinin genes (vlhA) that play a role in immune escape. The vlhA promoters have a set of distinct properties in comparison to promoters of the other genes. The vlhA promoters carry a variable GAA repeats region at approximately 40 nts upstream of transcription start site. The promoters have been considered active only in the presence of exactly 12 GAA repeats. The mechanisms of vlhA expression regulation and GAA number variation are not described. Here we tried to understand these mechanisms using different computational methods. We conducted a comparative analysis among several M. gallisepticum strains. Nucleotide sequences analysis showed the presence of highly conserved regions flanking repeated trinucleotides that are not linked to GAA number variation. VlhA genes with 12 GAA repeats and their orthologs in 12 M. gallisepticum strains are more conserved than other vlhA genes and have narrower GAA number distribution. We conducted comparative analysis of physicochemical profiles of M. gallisepticum vlhA and sigma-70 promoters. Stress-induced duplex destabilization (SIDD) profiles showed that sigma-70 group is characterized by the common to prokaryotic promoters sharp maxima while vlhA promoters are hardly destabilized with the region between GAA repeats and transcription start site having zero opening probability. Electrostatic potential profiles of vlhA promoters indicate the presence of the distinct patterns that appear to govern initial stages of specific DNA-protein recognition. Open state dynamics profiles of vlhA demonstrate the pattern that might facilitate transcription bubble formation. Obtained data could be the basis for experimental identification of mechanisms of phase variation in M. gallisepticum.

Dynamical Scaling and Phase Coexistence in Topologically Constrained DNA Melting

There is a long-standing experimental observation that the melting of topologically constrained DNA, such as circular closed plasmids, is less abrupt than that of linear molecules. This finding points to an important role of topology in the physics of DNA denaturation, which is, however, poorly understood. Here, we shed light on this issue by combining large-scale Brownian dynamics simulations with an analytically solvable phenomenological Landau mean field theory. We find that the competition between melting and supercoiling leads to phase coexistence of denatured and intact phases at the single-molecule level. This coexistence occurs in a wide temperature range, thereby accounting for the broadening of the transition. Finally, our simulations show an intriguing topology-dependent scaling law governing the growth of denaturation bubbles in supercoiled plasmids, which can be understood within the proposed mean field theory.

Biosynthesis and Regulation of the Branched-Chain Amino Acids†

This review focuses on more recent studies concerning the systems biology of branched-chain amino acid biosynthesis, that is, the pathway-specific and global metabolic and genetic regulatory networks that enable the cell to adjust branched-chain amino acid synthesis rates to changing nutritional and environmental conditions. It begins with an overview of the enzymatic steps and metabolic regulatory mechanisms of the pathways and descriptions of the genetic regulatory mechanisms of the individual operons of the isoleucine-leucine-valine (ilv) regulon. This is followed by more-detailed discussions of recent evidence that global control mechanisms that coordinate the expression of the operons of this regulon with one another and the growth conditions of the cell are mediated by changes in DNA supercoiling that occur in response to changes in cellular energy charge levels that, in turn, are modulated by nutrient and environmental signals. Since the parallel pathways for isoleucine and valine biosynthesis are catalyzed by a single set of enzymes, and because the AHAS-catalyzed reaction is the first step specific for valine biosynthesis but the second step of isoleucine biosynthesis, valine inhibition of a single enzyme for this enzymatic step might compromise the cell for isoleucine or result in the accumulation of toxic intermediates. The operon-specific regulatory mechanisms of the operons of the ilv regulon are discussed in the review followed by a consideration and brief review of global regulatory proteins such as integration host factor (IHF), Lrp, and CAP (CRP) that affect the expression of these operons.

Theoretical analysis of competing conformational transitions in superhelical DNA

We develop a statistical mechanical model to analyze the competitive behavior of transitions to multiple alternate conformations in a negatively supercoiled DNA molecule of kilobase length and specified base sequence. Since DNA superhelicity topologically couples together the transition behaviors of all base pairs, a unified model is required to analyze all the transitions to which the DNA sequence is susceptible. Here we present a first model of this type. Our numerical approach generalizes the strategy of previously developed algorithms, which studied superhelical transitions to a single alternate conformation. We apply our multi-state model to study the competition between strand separation and B-Z transitions in superhelical DNA. We show this competition to be highly sensitive to temperature and to the imposed level of supercoiling. Comparison of our results with experimental data shows that, when the energetics appropriate to the experimental conditions are used, the competition between these two transitions is accurately captured by our algorithm. We analyze the superhelical competition between B-Z transitions and denaturation around the c-myc oncogene, where both transitions are known to occur when this gene is transcribing. We apply our model to explore the correlation between stress-induced transitions and transcriptional activity in various organisms. In higher eukaryotes we find a strong enhancement of Z-forming regions immediately 5′ to their transcription start sites (TSS), and a depletion of strand separating sites in a broad region around the TSS. The opposite patterns occur around transcript end locations. We also show that susceptibility to each type of transition is different in eukaryotes and prokaryotes. By analyzing a set of untranscribed pseudogenes we show that the Z-susceptibility just downstream of the TSS is not preserved, suggesting it may be under selection pressure.

The stresses imposed on DNA within organisms can drive the molecule from its standard B-form double-helical structure into other conformations at susceptible sites within the sequence. We present a theoretical method to calculate this transition behavior due to stresses induced by supercoiling. We also develop a numerical algorithm that calculates the transformation probability of each base pair in a user-specified DNA sequence under stress. We apply this method to analyze the competition between transitions to strand separated and left-handed Z-form structures. We find that these two conformations are both competitive under physiological environmental conditions, and that this competition is especially sensitive to temperature. By comparing its results to experimental data we also show that the algorithm properly describes the competition between melting and Z-DNA formation. Analysis of large gene sets from various organisms shows a correlation between sites of stress-induced transitions and locations that are involved in regulating gene expression.

Theoretical analysis of the stress induced B-Z transition in superhelical DNA

We present a method to calculate the propensities of regions within a DNA molecule to transition from B-form to Z-form under negative superhelical stresses. We use statistical mechanics to analyze the competition that occurs among all susceptible Z-forming regions at thermodynamic equilibrium in a superhelically stressed DNA of specified sequence. This method, which we call SIBZ, is similar to the SIDD algorithm that was previously developed to analyze superhelical duplex destabilization. A state of the system is determined by assigning to each base pair either the B- or the Z-conformation, accounting for the dinucleotide repeat unit of Z-DNA. The free energy of a state is comprised of the nucleation energy, the sequence-dependent B-Z transition energy, and the energy associated with the residual superhelicity remaining after the change of twist due to transition. Using this information, SIBZ calculates the equilibrium B-Z transition probability of each base pair in the sequence. This can be done at any physiologically reasonable level of negative superhelicity. We use SIBZ to analyze a variety of representative genomic DNA sequences. We show that the dominant Z-DNA forming regions in a sequence can compete in highly complex ways as the superhelicity level changes. Despite having no tunable parameters, the predictions of SIBZ agree precisely with experimental results, both for the onset of transition in plasmids containing introduced Z-forming sequences and for the locations of Z-forming regions in genomic sequences. We calculate the transition profiles of 5 kb regions taken from each of 12,841 mouse genes and centered on the transcription start site (TSS). We find a substantial increase in the frequency of Z-forming regions immediately upstream from the TSS. The approach developed here has the potential to illuminate the occurrence of Z-form regions in vivo, and the possible roles this transition may play in biological processes.

Strand separation in negatively supercoiled DNA

We consider Benham's model for strand separation in negatively supercoiled circular DNA, and study denaturation as function of the linking difference density kappa<0. We propose a statistical version of this model, based on bayesian segmentation methods of current use in bioinformatics; this leads to new algorithms with priors adapted to supercoiled DNA, taking into account the random nature of the free energies needed to denature base pairs.

The analysis of stress-induced duplex destabilization in long genomic DNA sequences

We present a method for calculating predicted locations and extents of stress-induced DNA duplex destabilization (SIDD) as functions of base sequence and stress level in long DNA molecules. The base pair denaturation energies are assigned individually, so the influences of near neighbors, methylated bases, adducts, or lesions can be included. Sample calculations indicate that copolymeric energetics give results that are close to those derived when full near-neighbor energetics are used; small but potentially informative differences occur only in the calculated SIDD properties of moderately destabilized regions. The method presented here for analyzing long sequences calculates the destabilization properties within windows of fixed length N, with successive windows displaced by an offset distance d(o). The final values of the relevant destabilization parameters for each base pair are calculated as weighted averages of the values computed for each window in which that base pair appears. This approach implicitly assumes that the strength of the direct coupling between remote base pairs that is induced by the imposed stress attenuates with their separation distance. This strategy enables calculations of the destabilization properties of DNA sequences of any length, up to and including complete chromosomes. We illustrate its utility by calculating the destabilization properties of the entire E. coli genomic DNA sequence. A preliminary analysis of the results shows that promoters are associated with SIDD regions in a highly statistically significant manner, suggesting that SIDD attributes may prove useful in the computational prediction of promoter locations in prokaryotes.

Stress-induced DNA duplex destabilization (SIDD) in the E. coli genome: SIDD sites are closely associated with promoters

We present the first analysis of stress-induced DNA duplex destabilization (SIDD) in a complete chromosome, the Escherichia coli K12 genome. We used a newly developed method to calculate the locations and extents of stress-induced destabilization to single-base resolution at superhelix density sigma = -0.06. We find that SIDD sites in this genome show a statistically highly significant tendency to avoid coding regions. And among intergenic regions, those that either contain documented promoters or occur between divergently transcribing coding regions, and hence may be inferred to contain promoters, are associated with strong SIDD sites in a statistically highly significant manner. Intergenic regions located between convergently transcribing genes, which are inferred not to contain promoters, are not significantly enriched for destabilized sites. Statistical analysis shows that a strongly destabilized intergenic region has an 80% chance of containing a promoter, whereas an intergenic region that does not contain a strong SIDD site has only a 24% chance. We describe how these observations may illuminate specific mechanisms of regulation, and assist in the computational identification of promoter locations in prokaryotes.

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.

DNA topology-mediated control of global gene expression in Escherichia coli

Because the level of DNA superhelicity varies with the cellular energy charge, it can change rapidly in response to a wide variety of altered nutritional and environmental conditions. This is a global alteration, affecting the entire chromosome and the expression levels of all operons whose promoters are sensitive to superhelicity. In this way, the global pattern of gene expression may be dynamically tuned to changing needs of the cell under a wide variety of circumstances. In this article, we propose a model in which chromosomal superhelicity serves as a global regulator of gene expression in Escherichia coli, tuning expression patterns across multiple operons, regulons, and stimulons to suit the growth state of the cell. This model is illustrated by the DNA supercoiling-dependent mechanisms that coordinate basal expression levels of operons of the ilv regulon both with one another and with cellular growth conditions.

Single and multiple topologically driven structural transitions in DNA

We derive some exact general results concerning the behavior of topological absorbers (i.e., sequences undergoing topologically driven structural transitions) in closed circular DNA molecules. Starting from the formal physical framework that covers all known structural transitions, like those from standard B-DNA to nonstandard conformers Z-DNA, H-DNA, cruciform-DNA, melt-DNA or others, we develop a reduced state space description that leads to an analytically simplified "black box" view of absorbers. The latter contains only a single state variable-the total sequence unwinding u describing the topological state of the absorber. We show that the statistical mechanics of u is determined by the (one-dimensional) absorption free energy function G(abs) and find explicit expressions for G(abs) and for moments <u(n)> in terms of the standard experimental observable-the absorption function alpha:=<u>. The reduced state space method is then applied to systems consisting of several interacting topologically coupled absorbers and a formula predicting their collective behavior (superposition) in terms of their individual absorptions is derived. Using these results we formulate and discuss solution methods for two basic types of inverse problems that turn out to be fundamental for future absorber construction.

An initiation element in the yeast CUP1 promoter is recognized by RNA polymerase II in the absence of TATA box-binding protein if the DNA is negatively supercoiled

Purified RNA polymerase II initiated transcription from the yeast CUP1 promoter fused to a C-less cassette if the DNA was negatively supercoiled. Relaxed plasmid was not transcribed. Transcription did not require addition of any other transcription factors. TATA box-binding protein (TBP) was not detectable in the polymerase preparation and the TATA box was not required. Deletion analysis of the CUP1 promoter revealed that a 25-bp element containing the initiation region was sufficient for recognition by polymerase. Two transcription start sites were mapped, one of which is identical to one of the two major start sites observed in vivo. Our observations can be accounted for by using a theoretical analysis of the probability of DNA melting within the plasmid as a function of superhelix density: the CUP1 initiation element is intrinsically unstable to superhelical stress, permitting entry of the polymerase, which then scans the DNA to locate the start site. In support of this analysis, the CUP1 promoter was sensitive to mung bean nuclease. These observations and a previous theoretical analysis of yeast genes support the idea that promoters are stress points within the DNA superhelix. The role of transcription factors might be to mark the promoter and to regulate specific melting of promoter DNA.

Free energy and information contents of conformons in proteins and DNA

Sequence-specific conformational strains (SSCS) of biopolymers that carry free energy and genetic information have been called conformons, a term coined independently by two groups over two and a half decades ago [Green, D.E., Ji, S., 1972. The electromechanochemical model of mitochondrial structure and function. In: Schultz, J., Cameron, B.F. (Eds.), Molecular Basis of Electron Transport. Academic Press, New York, pp. 1-44; Volkenstein, M.V., 1972. The Conformon. J. Theor. Biol. 34, 193-195]. Conformons provide the molecular mechanisms necessary and sufficient to account for all biological processes in the living cell on the molecular level in principle--including the origin of life, enzymic catalysis, control of gene expression, oxidative phosphorylation, active transport, and muscle contraction. A clear example of SSCS is provided by SIDD (strain-induced duplex destabilization) in DNA recently reported by Benham [Benham, C.J., 1996a. Duplex destabilization in superhelical DNA is predicted to occur at specific transcriptional regulatory regions. J. Mol. Biol. 255, 425-434; Benham, C.J., 1996b. Computation of DNA structural variability--a new predictor of DNA regulatory regions. CABIOS 12(5), 375-381]. Experimental as well as theoretical evidence indicates that conformons in proteins carry 8-16 kcal/mol of free energy and 40-200 bits of information, while those in DNA contain 500-2500 kcal/mol of free energy and 200-600 bits of information. The similarities and differences between conformons and solitons have been analyzed on the basis of the generalized Franck-Condon principle [Ji, S., 1974a. A general theory of ATP synthesis and utilization. Ann. N.Y. Acad. Sci. 227, 211-226; Ji, S., 1974b. Energy and negentropy in enzymic catalysis. Ann. N.Y. Acad. Sci. 227, 419-437]. To illustrate a practical application, the conformon theory was applied to the molecular-clamp model of DNA gyrase proposed by Berger and Wang [Berger, J.M., Wang, J.C., 1996. Recent developments in DNA topoisomerases II structure and mechanism. Curr. Opin. Struct. Biol. 6(1), 84-90], leading to the proposal of an eight-step molecular mechanism for the action of the enzyme. Finally, a set of experimentally testable predictions has been formulated on the basis of the conformon theory.

Electrostatic-undulatory theory of plectonemically supercoiled DNA

We present an analytical calculation of the electrostatic interaction in a plectonemic supercoil within the Poisson-Boltzmann approximation. Undulations of the supercoil strands arising from thermal motion couple nonlinearly with the electrostatic interaction, giving rise to a strong enhancement of the bare interaction. In the limit of fairly tight winding, the free energy of a plectonemic supercoil may be split into an elastic contribution containing the bending and torsional energies and an electrostatic-undulatory free energy. The total free energy of the supercoil is minimized according to an iterative scheme, which utilizes the special symmetry inherent in the usual elastic free energy of the plectoneme. The superhelical radius, opening angle, and undulation amplitudes in the radius and pitch are obtained as a function of the specific linking difference and the concentration of monovalent salt. Our results compare favorably with the experimental values for these parameters of Boles et al. (1990. J. Mol. Biol. 213:931-951). In particular, we confirm the experimental observation that the writhe is a virtually constant fraction of the excess linking number over a wide range of superhelical densities. Another important prediction is the ionic strength dependence of the plectonemic parameters, which is in reasonable agreement with the results from computer simulations.

Denaturation of supercoiled DNA: a Monte Carlo study

A theoretical investigation of the denaturation characteristics of a supercoiled DNA has been presented employing a Metropolis Monte Carlo algorithm to examine the overall melting profiles of a supercoiled plasmid as the temperature is varied. We show that in contrast to a previously presented algorithm, this much simpler method is sufficient to explain almost all the overall denaturation characteristics and it also correctly calculates the detailed denaturation probabilities of each base pair at various degrees of supercoiling. We also present for the first time a theoretical investigation of the alkaline denaturation of a supercoiled plasmid. Although one can qualitatively reproduce the denaturation profiles using the present Monte Carlo algorithm, the agreement with experiment is not as good as in the case of thermal denaturation. The possible sources of discrepancy between theory and experiment have been discussed.

Activation of gene expression by a novel DNA structural transmission mechanism that requires supercoiling-induced DNA duplex destabilization in an upstream activating sequence

We have previously demonstrated that integration host factor (IHF)-mediated activation of transcription from the ilvPG promoter of Escherichia coli requires a supercoiled DNA template and occurs in the absence of specific interactions between IHF and RNA polymerase. In this report, we describe a novel, supercoiling-dependent, DNA structural transmission mechanism for this activation. We provide theoretical evidence for a supercoiling-induced DNA duplex destabilized (SIDD) structure in the A + T-rich, ilvPG regulatory region between base pair positions +1 and -160. We show that the region of this SIDD sequence immediately upstream of an IHF binding site centered at base pair position -92 is, in fact, destabilized by superhelical stress and that this duplex destabilization is inhibited by IHF binding. Thus, in the presence of IHF, the negative superhelical twist normally absorbed by this DNA structure in the promoter distal half of the SIDD sequence is transferred to the downstream portion of the SIDD sequence containing the ilvPG promoter site. This IHF-mediated translocation of superhelical energy facilitates duplex destabilization in the -10 region of the downstream ilvPG promoter and activates transcription by increasing the rate of open complex formation.

Stress-induced duplex DNA destabilization in scaffold/matrix attachment regions

S/MARs are DNA elements 300 to several thousand base-pairs long, which are operationally defined by their affinity for the nuclear scaffold or matrix. S/MARs occur exclusively in eukaryotic genomes, where they mediate several functions. Because S/MARs do not have a clearcut consensus sequence, the characteristics that define their activity are thought to be structural. Ubiquitous S/MAR binding proteins have been identified, but to date no unique binding sequence or structural motif has been found. Here we show by computational analysis that S/MARs conform to a specific design whose essential attribute is the presence of stress-induced base-unpairing regions (BURs). Stress-induced destabilization (SIDD) profiles are calculated using a previously developed statistical mechanical procedure in which the superhelical deformation is partitioned between strand separation, twisting within denatured regions, and residual superhelicity. The results of these calculations show that BURs exhibit a succession of evenly spaced destabilized sites that would render part or all of the S/MAR sequence single stranded at sufficient superhelicity. These analyses are performed for a range of sequenced S/MAR elements from the borders of eukaryotic gene domains, from centromeres, and from positions where S/MARs are known to support the action of an enhancer. The results reported here are in excellent agreement with earlier in vitro chemical reactivity studies. This approach demonstrates the potential for computational analysis to predict the points of division of the eukaryotic genome into functional units (domains), and also to locate certain cis-regulatory sequences.

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.

Energetics of the strand separation transition in superhelical DNA

In this paper the values of three free energy parameters governing the superhelical strand separation transition are determined by analysis of available experimental data. These are the free energy, a, needed to initiate a run of separation, the torsional stiffness, C, associated with interstrand winding of the two single strands comprising a separated site and the coefficient, K, of the quadratic free energy associated to residual linking. The experimental data used in this analysis are the locations and relative amounts of strand separation occurring in the pBR322 DNA molecule and the measured residual linking, both evaluated over a range of negative linking differences. The analytic method used treats strand separation as a heteropolymeric, co-operative, two-state transition to a torsionally deformable alternative conformation, which takes place in a circular DNA molecule constrained by the constancy of its linking number. The values determined for these parameters under the experimental conditions (T = 310 K, pH = 7.0, monovalent cation concentration = 0.01 M) are a = 10.84(+/- 0.2) kcal/mol, C = 2.5(+/- 0.3) x 10(-13) erg/rad2 and K = 2350(+/- 80) RT/N, where N is the molecular length in base-pairs. In order to assess the accuracy of the author's theoretical methods, these free energy parameters are incorporated into the analysis of superhelical strand separation in different molecules and under other conditions than those used in their evaluation. First, the temperature dependence of transition is treated, then superhelical strand separation is analyzed in a series of DNA molecules having systematic sequence modifications, and the results of these theoretical analyses are compared with those from experiments. In all molecules, transition is predicted in the range of linking differences where it is seen experimentally. Moreover, it occurs at the specific sequence locations that the analysis predicts, and with approximately the predicted relative amounts of transition at each location. The known sensitivities of this transition to changes of temperature and to small sequence modifications are predicted in a quantitatively precise manner by the theoretical results. The demonstrated high-level precision of these theoretical methods provides a tool for the screening of DNA sequences for sites susceptible to superhelical strand separation, some of which may have regulatory or other biological significance.

Negative supercoiling increases the sensitivity of plasmid DNA to single-strand break induction by X-rays

Negatively supercoiled topoisomers of the plasmid pIBI 30 were irradiated with 250 kV X-rays and assayed for strand scission by agarose gel electrophoresis. The survival of supercoiled molecules (Form I) decreased exponentially with increasing X-ray exposure and the dose required to reduce the fraction of DNA in Form I to 37% of its value in unirradiated controls (D37) decreased with increasing negative superhelicity. This enhanced radiation sensitivity of underwound DNA is tentatively attributed to the transient denaturation of the double helix that increases the susceptibility of individual strands to free radical attack.