Effect of base composition at the center of inverted repeated DNA sequences on cruciform transitions in DNA

We have analyzed the effect of base composition at the center of symmetry of inverted repeated DNA sequences on cruciform transitions in supercoiled DNA. For this we have constructed two series of palindromic DNA sequences: one set with differing center and one set with differing center and arm sequences. The F series consists of two 96-base pair perfect inverted repeats which are identical except for the central 10 base pairs which consist of pure AT or GC base pairs. The S series was constructed such that the overall base composition of the inverted repeats was identical but in which the positioning of blocks of AT- and GC-rich sequences varied. The rate of cruciform formation for the inverted repeats in plasmid pUC8 was dramatically influenced by the 8-10 base pairs at the center of the inverted repeat. Inverted repeats with 8-10 AT base pairs in the center were kinetically much more active in cruciform formation than inverted repeats with 8-10 GC base pairs in the center. These experiments show a dominant influence of the center sequences of inverted repeats on the rate of cruciform formation.     

Counterintuitive DNA Sequence Dependence in Supercoiling-Induced DNA Melting

The metabolism of DNA in cells relies on the balance between hybridized double-stranded DNA (dsDNA) and local de-hybridized regions of ssDNA that provide access to binding proteins. Traditional melting experiments, in which short pieces of dsDNA are heated up until the point of melting into ssDNA, have determined that AT-rich sequences have a lower binding energy than GC-rich sequences. In cells, however, the double-stranded backbone of DNA is destabilized by negative supercoiling, and not by temperature. To investigate what the effect of GC content is on DNA melting induced by negative supercoiling, we studied DNA molecules with a GC content ranging from 38% to 77%, using single-molecule magnetic tweezer measurements in which the length of a single DNA molecule is measured as a function of applied stretching force and supercoiling density. At low force (<0.5pN), supercoiling results into twisting of the dsDNA backbone and loop formation (plectonemes), without inducing any DNA melting. This process was not influenced by the DNA sequence. When negative supercoiling is introduced at increasing force, local melting of DNA is introduced. We measured for the different DNA molecules a characteristic force F _char, at which negative supercoiling induces local melting of the dsDNA. Surprisingly, GC-rich sequences melt at lower forces than AT-rich sequences: F _char = 0.56pN for 77% GC but 0.73pN for 38% GC. An explanation for this counterintuitive effect is provided by the realization that supercoiling densities of a few percent only induce melting of a few percent of the base pairs. As a consequence, denaturation bubbles occur in local AT-rich regions and the sequence-dependent effect arises from an increased DNA bending/torsional energy associated with the plectonemes. This new insight indicates that an increased GC-content adjacent to AT-rich DNA regions will enhance local opening of the double-stranded DNA helix.     

A thermal denaturation study of genomic DNAs from North American minnows (Cyprinidae: Teleostei)

Base compositions and differential melting rate profiles of genomic DNAs from twenty species of North American cyprinid fishes were generated via thermal denaturation. Base pair composition expressed as % GC values ranged among the twenty species from 36.1–41.3%. This range is considerably broader than that observed at comparable taxonomic levels in other vertebrate groups. Both the range and average difference in base pair composition between species in the diverse and rapidly evolving genus Notropis were considerably greater than those between species in other North American cyprinid genera. This may indicate that genomic changes at the level of base pair composition are frequent and possibly important events in cyprinid evolution. Compositional heterogeneity and asymmetry values among the twenty species were uniform and low, respectively, suggesting that most of the species lacked DNA components in their genomes which differed substantially from their main-band DNAs in base pair composition. The melting rate profiles revealed a prominent and distinct heavy or GC-rich DNA component in the genomes of three species belonging to the subgenus Cyprinella of Notropis. These and other data suggest that the heavy melting component may reflect a large, comparatively GC-rich family of highly repeated or satellite DNA sequences common to all three genomes.     

Calculation of melting curves for DNA

A method is reported for calculating the melting curve of a DNA molecule of random base sequence, including in the formalism the dependence of the free energy of base pair formation on the size of a denatured section. Some explicit results are shown for a “typical” base sequence, in particular the probability of helix formation at individual base pairs in several different regions of the molecule and the amount of melting from the end of the chain. Particular attention is drawn to the variation of local melting behavior from one region of the molecule to another. It is found that sections rich in AT melt at relatively low temperatures with a fairly broad transition curve, whereas regions rich in GC pairs melt at higher temperatures (as expected) with a very abrupt, local transition curve. To account qualitatively for the results one may divide melting into two kinds of processes: (a) the nucleation and growth of denatured regions, and (b) the merging together of two denatured sections at the expense of the intervening helix. The first of these processes dominates in the first stages of melting, and leads to rather broad local melting curves, whereas the second process predominates in the later stages, and occurs, in a particular part of the molecule, over a very narrow temperature range. It is estimated that the average length of a helix plus adjacent coil section at the midpoint of the transition is approximately 600 base pairs. Since transition curves which measure the local melting behavior reflect local compositions fluctuations, these curves contain information about the broad outlines of base sequence in the molecule. Some suggestions are made concerning experiments by which this potential information source could be exploited. In particular, it is pointed out that one might hope to map AT or GC rich regions at particular genetic loci in a biologically active DNA molecule. Values of the relevant parameters found earlier for the transition of homopolymers produce melting curves for a DNA of random base sequence which are in good agreement with the experimental transition curve for T2 phage DNA. Hence the present theoretical picture of the melting of polynucleotides is at least internally self-consistent.     

The organization of DNA in the mitotic and polytene chromosomes of Sciara coprophila

The organization of DNA in the mitotic metaphase and polytene chromosomes of the fungus gnat, Sciara coprophila, has been studied using base-specific DNA ligands, including anti-nucleoside antibodies. The DNA of metaphase and polytene chromosomes reacts with AT-specific probes (quinacrine, DAPI, Hoechst 33258 and anti-adenosine) and to a somewhat lesser extent with GC-specific probes (mithramycin, chromomycin A₃ and anticytidine). In virtually every band of the polytene chromosomes chromomycin A₃ fluorescence is almost totally quenched by counterstaining with the AT-specific ligand methyl green. This indicates that GC base pairs in most bands are closely interspersed with AT base pairs. The only exceptions are band IV-8A3 and the nucleolus organizer on the X. In contrast, quinacrine and DAPI fluorescence in every band is only slightly quenched by counterstaining with the GC-specific ligand actinomycin D. Thus, each band contains a moderate proportion of AT-rich DNA sequences with few interspersed GC base pairs. — The C-bands in mitotic and polytene chromosomes can be visualized by Giemsa staining after differential extraction of DNA and those in polytene chromosomes by the use of base-specific fluorochromes or antibodies without prior extraction of DNA. C-bands are located in the centromeric region of every chromosome, and the telomeric region of some. The C-bands in the polytene chromosomes contain AT-rich DNA sequences without closely interspered GC base pairs and lack relatively GC-rich sequences. However, one C-band in the centromeric region of chromosome IV contains relatively GC-rich sequences with closely interspersed AT base pairs. — C-bands make up less than 1% of polytene chromosomes compared to nearly 20% of mitotic metaphase chromosomes. The C-bands in polytene chromosomes are detectable with AT-specific or GC-specific probes while those in metaphase chromosomes are not. Thus, during polytenization there is selective replication of highly AT-rich and relatively GC-rich sequences and underreplication of the remainder of the DNA sequences in the constitutive heterochromatin.     

Non-neutral processes drive the nucleotide composition of non-coding sequences in Drosophila

The nature of the forces affecting base composition is a key question in genome evolution. There is uncertainty as to whether differences in the GC contents of non-coding sequences reflect differences in mutational bias, or in the intensity of selection or biased gene conversion. We have used a polymorphism dataset for non-coding sequences on the X chromosome of Drosophila simulans to examine this question. The proportion of GC-->AT versus AT-->GC polymorphic mutations in a locus is correlated with its GC content. This implies the action of forces that favour GC over AT base pairs, which are apparently strongest in GC-rich sequences.     

The thermal denaturation of DNA: average length and composition of denatured areas

A spectral study of melting curves of DNA ranging from 73 to 32% AT indicates that the base ratio of sequences melting within DNA are a linear function of temperature. A study of partially denatured DNA by electron microscopy, reversible renaturation and fractionation on hydroxylapatite suggests that the melting curve of DNA represents the melting of sequences which average 3-4 million daltons in length. These sequences appear to be a combination of two areas, one which is high in AT and denatures in the first three-quarters of the melting curve, and one which is high in GC and denatures in the final quarter. The length of these sequences appears to vary between 1.5-6 million daltons.     

The influence of DNA base ratios on intragenic mutational spectra in prokaryotes

Because of the degeneracy of the Genetic Code, a given amino acid sequence can be written in terms of either GC- or AT-rich DNA-specified mRNA codons. It is shown that the genetic consequences in terms of detectable single base pair substitution mutations are very different in such GC- or AT-rich sequences. Furthermore the occurrences of potential frameshift mutation hot-spots, in runs of repeated or reiterated base pairs, and sites of mutationally important base methylation, also differ markedly between AT- and GC-rich sequences. Thus the evolutionary pathways which can be followed by using single or multiple changes in the amino acid sequence of a given polypeptide will be very different in prokaryotes with GC-rich or AT-rich DNA.     

Modification of the melting properties of duplex DNA by attachment of a GC-rich DNA sequence as determined by denaturing gradient gel electrophoresis.

The melting behavior of a DNA fragment carrying the mouse beta maj-globin promoter was investigated as a means of establishing procedures for separating DNA fragments differing by any single base substitution using the denaturing gradient gel electrophoresis procedure of Fischer and Lerman (1,2). We find that attachment of a 300 base pair GC-rich DNA sequence, termed a GC-clamp, to a 135 bp DNA fragment carrying the mouse beta-globin promoter significantly alters the pattern of DNA melting within the promoter. When the promoter is attached to the clamp, the promoter sequences melt without undergoing strand dissociation. The calculated distribution of melting domains within the promoter differs markedly according to the relative orientation of the clamp and promoter sequences. We find that the behavior of DNA fragments containing the promoter and clamp sequences on denaturing gradient polyacrylamide gels is in close agreement with the theoretical melting calculations. These studies provide the basis for critical evaluation of the parameters for DNA melting calculations, and they establish conditions for determining whether all single base substitutions within the promoter can be separated on denaturing gradient gels.     

Conformation and reactivity of DNA in the complex with proteins. III. Helix-coil transition and conformational studies of model complexes of DNA''s with poly-L-histidine.

Differences in the interaction of poly-L-histidine with DNA of various base composition have been demonstrated using melting and CD measurements. The two types of complexes formed with DNA at pH values below the pK of 5.9 and in the region of pH 6.5 are very different in their CD spectral properties. The binding effects with highly protonated poly-L-histidine are AT-dependent as reflected by large negative CD spectra indicating the formation of psi-DNA as a condensed state of the double helix. GC-rich DNA may, however, also form psi-DNA structures with poly-L-histidine under certain conditions. At pH 6.5 complex formation with the weakly protonated polypeptide is GC-dependent. From the results it is concluded that protonated poly-L-histidine interacts more specifically at AT base pairs, prabably along the small groove while the weakly protonated poly-L-histidine tends to interact preferentially with GC regions which seems to occur rather in the large groove.     

Heterogeneity in E. coli RecBCD Helicase-DNA Binding and Base Pair Melting

E. coli RecBCD, a helicase/nuclease involved in double stranded (ds) DNA break repair, binds to a dsDNA end and melts out several DNA base pairs (bp) using only its binding free energy. We examined RecBCD-DNA initiation complexes using thermodynamic and structural approaches. Measurements of enthalpy changes for RecBCD binding to DNA ends possessing pre-melted ssDNA tails of increasing length suggest that RecBCD interacts with ssDNA as long as 17–18 nucleotides and can melt at least 10–11 bp upon binding a blunt DNA end. Cryo-EM structures of RecBCD alone and in complex with a blunt-ended dsDNA show significant conformational heterogeneities associated with the RecB nuclease domain (RecB^Nuc) and the RecD subunit. In the absence of DNA, 56% of RecBCD molecules show no density for the RecB nuclease domain, RecB^Nuc, and all RecBCD molecules show only partial density for RecD. DNA binding reduces these conformational heterogeneities, with 63% of the molecules showing density for both RecD and RecB^Nuc. This suggests that the RecB^Nuc domain is dynamic and influenced by DNA binding. The major RecBCD-DNA structural class in which RecB^Nuc is docked onto RecC shows melting of at least 11 bp from a blunt DNA end, much larger than previously observed. A second structural class in which RecB^Nuc is not docked shows only four bp melted suggesting that RecBCD complexes transition between states with different extents of DNA melting and that the extent of melting regulates initiation of helicase activity.     

Segments of mitochondrial DNA of yeast carrying antibiotic resistances

Summary Mitochondrial DNA was isolated from an oligomycin-resistant petite mutant of yeast, Saccharomyces cerevisiae. It had repeated sequences of 3600 base pairs. This segment was about one twentieth of the whole mtDNA of wild type yeast, which had a size of 74 kilo base pairs.This segment of mtDNA had one cleavage site for a restriction endonuclease, Hind II, which was more resistant to cleavage than the other Hind II sites in wild type mtDNA. It had two cleavage sites for Hha I and gave two Hha fragments, which were arranged alternatively. Digestion with Hae III gave four fragments and these fragments were mapped.Mitochondrial DNA of this mutant showed a loss of heterogeneity in a melting profile. It melted within a narrow range of temperature, which was similar to that of poly dA·poly dT. Its differential melting curve was significantly different from that of wild type mtDNA.Mapping of mtDNA of a wild type yeast was carried out with restriction endonucleases. Fragments of mtDNA, which were isolated from petites carrying oligomycin-erythromycin-chloramphenicol-resistance and erythromycin-chloramphenicol resistance were also mapped. Loci of oligomycin-resistance, erythromycin-resistance and chloramphenicol-resistance were investigated based on the maps of Eco R I fragments and Hind II fragments.     

Thermodynamics of DNA Containing Very Stable Chemically Modified Base Pairs

Abstract

DNA chemical modifications caused by the binding of some antitumor drugs give rise to a very strong local stabilization of the double helix. These sites melt at a temperature that is well above the melting temperatures of ordinary AT and GC base pairs. In this work we have examined the melting behavior of DNA containing very stable sites. Analytical expressions were derived and used to evaluate the thermodynamic properties of homopolymers DNA with several different distributions of stable sites. The results were extended to DNA with a heterogeneous sequence of AT and GC base pairs. The results were compared to the melting properties of DNA with ordinary covalent interstrand cross-links. It was found that, as with an ordinary interstrand cross-link, a single strongly stabilized site makes a DNA's melting temperature (T_m ) independent of strand concentration. However in contrast to a DNA with an interstrand cross-link, a strongly stabilized site makes the DNA's T_m independent of DNA length and equal to T _∞, the melting temperature of an infinite length DNA with the same GC-content and without a stabilized site. Moreover, at a temperature where more than 80% of base pairs are melted, the number of ordinary (non-modified) helical base pairs (n) is independent of both the DNA length and the location of the stabilized sites. For this condition, n(T) = (2ω-a) S (1- S ) and S = exp[ΔS(T∞-T)/(RT)] where ω is the number of strongly stabilized sites in the DNA chain, a is the number of DNA ends that contain a stabilized site, and ΔS, T, and R are the base pair entropy change, the temperature, and the universal gas constant per mole. The above expression is valid for a temperature interval that corresponds to n<0.2N for ω=1, and n<0.1N for ω>1, where N is the number of ordinary base pairs in the DNA chain.     

Nanomechanics of Diaminopurine-Substituted DNA

2,6-diaminopurine (DAP) is a nucleobase analog of adenine. When incorporated into double-stranded DNA (dsDNA), it forms three hydrogen bonds with thymine. Rare in nature, DAP substitution alters the physical characteristics of a DNA molecule without sacrificing sequence specificity. Here, we show that in addition to stabilizing double-strand hybridization, DAP substitution also changes the mechanical and conformational properties of dsDNA. Thermal melting experiments reveal that DAP substitution raises melting temperatures without diminishing sequence-dependent effects. Using a combination of atomic force microscopy (AFM), magnetic tweezer (MT) nanomechanical assays, and circular dichroism spectroscopy, we demonstrate that DAP substitution increases the flexural rigidity of dsDNA yet also facilitates conformational shifts, which manifest as changes in molecule length. DAP substitution increases both the static and dynamic persistence length of DNA (measured by AFM and MT, respectively). In the static case (AFM), in which tension is not applied to the molecule, the contour length of DAP-DNA appears shorter than wild-type (WT)-DNA; under tension (MT), they have similar dynamic contour lengths. At tensions above 60 pN, WT-DNA undergoes characteristic overstretching because of strand separation (tension-induced melting) and spontaneous adoption of a conformation termed S-DNA. Cyclic overstretching and relaxation of WT-DNA at near-zero loading rates typically yields hysteresis, indicative of tension-induced melting; conversely, cyclic stretching of DAP-DNA showed little or no hysteresis, consistent with the adoption of the S-form, similar to what has been reported for GC-rich sequences. However, DAP-DNA overstretching is distinct from GC-rich overstretching in that it happens at a significantly lower tension. In physiological salt conditions, evenly mixed AT/GC DNA typically overstretches around 60 pN. GC-rich sequences overstretch at similar if not slightly higher tensions. Here, we show that DAP-DNA overstretches at 52 pN. In summary, DAP substitution decreases the overall stability of the B-form double helix, biasing toward non-B-form DNA helix conformations at zero tension and facilitating the B-to-S transition at high tension.     

Interactions of DNA with divalent metal ions. IV. Competitive studies of Mn2+ binding to AT- and GC-rich DNAs

It has been inferred from previous studies that Mn²⁺ ions bind preferentially to G·C base pairs in DNA, and it has even been suggested that this preference for G·C pairs might be responsible for some of the Mn²⁺ specific effects observed in various biochemical reactions. In this paper we investigate the AT/GC preference of Mn²⁺ by direct competition studies in which AT-rich DNA was dialyzed against GC-rich DNA in the presence of varying amounts of Mn²⁺. Analysis of these results demonstrates that over a wide range of Mn²⁺/DNA(P) molar ratios, Mn²⁺ binds to A·T and to G·C base pairs with virtually identical affinity, although in a somewhat different mode. Both the present and previous nmr, uv, CD, and melting studies are discussed in terms of the different modes of binding of Mn²⁺ to single- and double-stranded DNA.     

Vanishing GC-rich isochores in mammalian genomes

To understand the origin and evolution of isochores-the peculiar spatial distribution of GC content within mammalian genomes-we analyzed the synonymous substitution pattern in coding sequences from closely related species in different mammalian orders. In primate and cetartiodactyls, GC-rich genes are undergoing a large excess of GC --> AT substitutions over AT --> GC substitutions: GC-rich isochores are slowly disappearing from the genome of these two mammalian orders. In rodents, our analyses suggest both a decrease in GC content of GC-rich isochores and an increase in GC-poor isochores, but more data will be necessary to assess the significance of this pattern. These observations question the conclusions of previous works that assumed that base composition was at equilibrium. Analysis of allele frequency in human polymorphism data, however, confirmed that in the GC-rich parts of the genome, GC alleles have a higher probability of fixation than AT alleles. This fixation bias appears not strong enough to overcome the large excess of GC --> AT mutations. Thus, whatever the evolutionary force (neutral or selective) at the origin of GC-rich isochores, this force is no longer effective in mammals. We propose a model based on the biased gene conversion hypothesis that accounts for the origin of GC-rich isochores in the ancestral amniote genome and for their decline in present-day mammals.     

Thermodynamics of DNA containing very stable chemically modified base pairs

DNA chemical modifications caused by the binding of some antitumor drugs give rise to a very strong local stabilization of the double helix. These sites melt at a temperature that is well above the melting temperatures of ordinary AT and GC base pairs. In this work we have examined the melting behavior of DNA containing very stable sites. Analytical expressions were derived and used to evaluate the thermodynamic properties of homopolymer DNA with several different distributions of stable sites. The results were extended to DNA with a heterogeneous sequence of AT and GC base pairs. The results were compared to the melting properties of DNA with ordinary covalent interstrand cross-links. It was found that, as with an ordinary interstrand cross-link, a single strongly stabilized site makes a DNA's melting temperature (T(m)) independent of strand concentration. However in contrast to a DNA with an interstrand cross-link, a strongly stabilized site makes the DNA's T(m) independent of DNA length and equal to T(infinity), the melting temperature of an infinite length DNA with the same GC-content and without a stabilized site. Moreover, at a temperature where more than 80% of base pairs are melted, the number of ordinary (non-modified) helical base pairs (n) is independent of both the DNA length and the location of the stabilized sites. For this condition, n(T) = (2 omega-a)S/(1-S) and S = exp[DeltaS(T(infinity)-T)/(RT)] where omega is the number of strongly stabilized sites in the DNA chain, a is the number of DNA ends that contain a stabilized site, and DeltaS, T, and R are the base pair entropy change, the temperature, and the universal gas constant per mole. The above expression is valid for a temperature interval that corresponds to n<0.2N for omega=1, and n<0.1N for omega>1, where N is the number of ordinary base pairs in the DNA chain.     

Long-range structural effects in supercoiled DNA: statistical thermodynamics reveals a correlation between calculated cooperative melting and contextual influence on cruciform extrusion

We have previously described [K. M. Sullivan and D. M. J. Lilley (1986) Cell 47, 817-827] a set of sequences, called C-type inducing sequences, which cause cruciform extrusion by adjacent inverted repeats to occur by an abnormal kinetic pathway involving a large denatured region of DNA. In this paper we apply statistical thermodynamic DNA helix melting theory to these sequences. We find a marked correlation between the ability of sequences to confer C-type cruciform character experimentally and their calculated propensity to undergo cooperative melting, and no exceptions have been found. The correlations are both qualitative and quantitative. Thus the ColE1 flanking sequences behave as single melting units, while the DNA of the S-type plasmid pIRbke8 exhibits no propensity to melt in the region of the bke cruciform. The results of the calculations are also fully consistent with the following experimental observations: 1. the ability of the isolated colL and colR fragments of the ColE1 flanking sequences, as well as the short sequence col30, to confer C-type character; 2. C-type induction by an A + T rich Drosophila sequence; 3. low-temperature cruciform extrusion by an (AT)34 sequence; 4. the effect of changing sequences at a site 90 base pairs (bp) removed from the inverted repeat; 5. the effects of systematic deletion of the colL sequence; and 6. the effects of insertion of various sequences in between the colL sequence and the xke inverted repeat. These studies show that telestability effects on thermal denaturation as predicted from equilibrium helix melting theory of linear DNA molecules may explain all the features that are revealed by studying the extrusion of cruciforms in circular DNA molecules subjected to superhelical stress.     

Study of DNA melting in the region of the inversion of relative stability of AT and GC pairs]

Systematic data on the dependence of the melting curve parameters of DNA from different organisms on the concentration of salt (C2H5)5NBr have been obtained. The melting curves were studied by spectrophotometric as well as by microcalorimetric methods. The DNA melting range width is shown to pass through the minimum value delta0T = 0.6 +/- 0.1 degrees at the point of inversion of relative stability of AT and GC pairs that corresponds to the concentration of (C2H5)4NBr equal to 2.9 +/- 0.1 M. This concentration, as well as the value of delta0T, are the same for different DNA's of common chemical structure. The T2 and T4 DNA containing hydroxymethylated and glucosylated cytosine residues show an anomalous behaviour. The enthalpy of melting falls very slowly as the salt concentration increases. The possible causes of the observed value of delta0T are discussed. A conclusion is drawn that the main factor which governs the DNA melting process in the region of inversion of the relative stability of AT and GC pairs is the heterogeneity of stacking interaction between different base pairs.