Correlations between scaffold/matrix attachment region (S/MAR) binding activity and DNA duplex destabilization energy

Scaffold or matrix-attachment regions (S/MARs) are thought to be involved in the organization of eukaryotic chromosomes and in the regulation of several DNA functions. Their characteristics are conserved between plants and humans, and a variety of biological activities have been associated with them. The identification of S/MARs within genomic sequences has proved to be unexpectedly difficult, as they do not appear to have consensus sequences or sequence motifs associated with them. We have shown that S/MARs do share a characteristic structural property, they have a markedly high predicted propensity to undergo strand separation when placed under negative superhelical tension. This result agrees with experimental observations, that S/MARs contain base-unpairing regions (BURs). Here, we perform a quantitative evaluation of the association between the ease of stress-induced DNA duplex destabilization (SIDD) and S/MAR binding activity. We first use synthetic oligomers to investigate how the arrangement of localized unpairing elements within a base-unpairing region affects S/MAR binding. The organizational properties found in this way are applied to the investigation of correlations between specific measures of stress-induced duplex destabilization and the binding properties of naturally occurring S/MARs. For this purpose, we analyze S/MAR and non-S/MAR elements that have been derived from the human genome or from the tobacco genome. We find that S/MARs exhibit long regions of extensive destabilization. Moreover, quantitative measures of the SIDD attributes of these fragments calculated under uniform conditions are found to correlate very highly (r2>0.8) with their experimentally measured S/MAR-binding strengths. These results suggest that duplex destabilization may be involved in the mechanisms by which S/MARs function. They suggest also that SIDD properties may be incorporated into an improved computational strategy to search genomic DNA sequences for sites having the necessary attributes to function as S/MARs, and even to estimate their relative binding strengths.     

Susceptibility to superhelically driven DNA duplex destabilization: a highly conserved property of yeast replication origins

Strand separation is obligatory for several DNA functions, including replication. However, local DNA properties such as A+T content or thermodynamic stability alone do not determine the susceptibility to this transition in vivo. Rather, superhelical stresses provide long-range coupling among the transition behaviors of all base pairs within a topologically constrained domain. We have developed methods to analyze superhelically induced duplex destabilization (SIDD) in genomic DNA that take into account both this long-range stress-induced coupling and sequence-dependent local thermodynamic stability. Here we apply this approach to examine the SIDD properties of 39 experimentally well-characterized autonomously replicating DNA sequences (ARS elements), which function as replication origins in the yeast Saccharomyces cerevisiae. We find that these ARS elements have a strikingly increased susceptibility to SIDD relative to their surrounding sequences. On average, these ARS elements require 4.78 kcal/mol less free energy to separate than do their immediately surrounding sequences, making them more than 2,000 times easier to open. Statistical analysis shows that the probability of this strong an association between SIDD sites and ARS elements arising by chance is approximately 4 × 10⁻¹⁰. This local enhancement of the propensity to separate to single strands under superhelical stress has obvious implications for origin function. SIDD properties also could be used, in conjunction with other known origin attributes, to identify putative replication origins in yeast, and possibly in other metazoan genomes.     

Computational and in vitro analysis of destabilized DNA regions in the interferon gene cluster: potential of predicting functional gene domains

Recent evidence adds support to a traditional concept according to which the eukaryotic nucleus is organized into functional domains by scaffold or matrix attachment regions (S/MARs). These regions have previously been predicted to have a high potential for stress-induced duplex destabilization (SIDD). Here we report the parallel results of binding (reassociation) and computational SIDD analyses for regions within the human interferon gene cluster on the short arm of chromosome 9 (9p22). To verify and further refine the biomathematical methods, we focus on a 10 kb region in the cluster with the pseudogene IFNWP18 and the interferon alpha genes IFNA10 and IFNA7. In a series of S/MAR binding assays, we investigate the promoter and termination regions and additional attachment sequences that were detected in the SIDD profile. The promoters of the IFNA10 and the IFNA7 genes have a moderate approximately 20% binding affinity to the nuclear matrix; the termination sequences show stronger association (70-80%) under our standardized conditions. No comparable destabilized elements were detected flanking the IFNWP18 pseudogene, suggesting that selective pressure acts on the physicochemical properties detected here. In extended, noncoding regions a striking periodicity is found of rather restricted SIDD minima with scaffold binding potential. By various criteria, the underlying sequences represent a new class of S/MARs, thought to be involved in a higher level organization of the genome. Together, these data emphasize the relevance of SIDD calculations as a valid approach for the localization of structural, regulatory, and coding regions in the eukaryotic genome.     

Molecular modeling of closed circular DNA thermodynamic ensembles

Many modeling studies of supercoiled DNA are based on equilibrium structures from theoretical calculations or energy minimization. Since closed circular DNAs are flexible, it is possible that errors are introduced by calculating properties from a single minimum energy structure, rather than from a complete thermodynamic ensemble. We have investigated this question using molecular dynamics simulations on a low resolution molecular mechanics model in which each base pair is represented by three points (a plane). This allows the inclusion of sequence-dependent variations of tip, inclination, and twist. Three kinds of sequences were tested: (1) homogeneous DNA, in which all base pairs have the helicoidal parameters of an ideal, average B-DNA; (2) random sequence DNA; and (3) curved DNA. We examined the rate of convergence of various structural parameters. Convergence for most of these is slowest for homogeneous sequences, more rapid for random sequences, and most rapid for curved sequences. The most slowly converging parameter is the antipodes profile. In a plasmid with N base pairs (bp), the antipodes distance is the distance d _ij from base pair i to base pair j halfway around the plasmid, j = i + N/2. The antipodes profile at time t is a plot of d _ij over the range i = 1, N/2. In a homogeneous plasmid, convergence requires that the antipodes profile averaged over time must be flat. Even in the small plasmids examined here, the average properties of the ensembles were found to differ from those of static equilibrium structures. These effects will be even more dramatic for larger plasmids. Further, average and dynamic properties are affected by both plasmid size and sequence. © 1996 John Wiley & Sons, Inc.     

The positive aspects of stress: strain initiates domain decondensation (SIDD).

The conventional string-based bioinformatic methods of genomic sequence analysis are often insufficient to identify DNA regulatory elements, since many of these do not have a recognizable motif. Even in case a sequence pattern is known to be associated with an element it may only partially mediate its function. This suggests that properties not correlated with the details of base sequence contribute to regulation. One of these attributes is the DNA strand-separation potential, known as SIDD (stress-induced duplex destabilization) which facilitates the access of tracking proteins and the formation of local secondary structures. Using the type 1 interferon gene cluster as a paradigm, we demonstrate that the imprints in a SIDD profile coincide with chromatin domain borders and with DNAse I hypersensitive sites to which regulatory potential could be assigned. The approach permits the computer-guided identification of yet unknown, mostly remote sites and the design of artificial elements with predictable properties for multiple applications.     

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments’ propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between ‘omics’ and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles’ local geometry consists of straight segments connected by sharp bends compressing the DNA’s inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles’ inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes’ nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.     

A theoretical study of formation of DNA noncanonical structures under negative superhelical stress

The development of statistical mechanical models of the formation of noncanonical structures in circular DNA and the finding of the energy parameters for these models made it possible to predict the appearance of such structures in a DNA with any given sequence. It does not seem feasible, however, to perform such calculations for DNA sequences of considerable length by allowing for all the possible states. We propose a special algorithm for calculating the thermodynamic characteristics of various conformational rearrangements in DNA that occur under negative supercoilings, allowing for several possible states of each base pair in the chain. Calculations have been performed for a number of natural DNAs. According to these calculations, the most likely noncanonical structures in DNA under normal conditions are cruciform structures and the Z form. The results of the calculations are compared with the experimental data reported in the literature. State diagrams have been computed for a number of inserts in circular DNA that can adopt both the cruciform conformation and the left-handed helical Z form.     

A Theoretical Study of Formation of DNA Noncanonical Structures Under Negative Superhelical Stress

Abstract

The development of statistical mechanical models of the formation of noncanonical structures in circular DNA and the finding of the energy parameters for these models made it possible to predict the appearance of such structures in a DNA with any given sequence. It does not seem feasible, however, to perform such calculations for DNA sequences of considerable length by allowing for all the possible states. We propose a special algorithm for calculating the thermodynamic characteristics of various conformational rearrangements in DNA that occur under negative supercoilings, allowing for several possible states of each base pair in the chain. Calculations have been performed for a number of natural DNAs. According to these calculations, the most likely noncanonical structures in DNA under normal conditions are cruciform structures and the Z form. The results of the calculations are compared with the experimental data reported in the literature. State diagrams have been computed for a number of inserts in circular DNA that can adopt both the cruciform conformation and the left-handed helical Z form.     

Long range structural communication between sequences in supercoiled DNA. Sequence dependence of contextual influence on cruciform extrusion mechanism

Sequence context may profoundly alter the character of structural transitions in supercoiled DNA (Sullivan, K. M., and Lilley, D. M. J. (1986) Cell 47, 817-827). The A + T-rich sequences of ColE1, which flank the inverted repeat, are responsible for cruciform extrusion following a mechanistic pathway which proceeds via a relatively large denatured region. This C-type mechanism results in kinetic properties which are very different from those of the S-type pathway, the normal mechanism of cruciform extrusion in the absence of the ColE1 flanking sequences. We have analyzed the sequence requirements for the induction of the C-type pathway. The 100-base pair left side sequence of ColE1 (colL) was subjected to systematic deletion using Bal31 exonucleolysis, showing that removal of 30 base pairs from its right end abolished extrusion by the C-type process. A cloned oligonucleotide of the same 30-base pair sequence was sufficient to confer C-type cruciform extrusion on an adjacent inverted repeat. An A + T-rich sequence from Drosophila was found to act like the ColE1 sequences. We have studied the effects of introducing sequences between the A + T-rich colL, and the inverted repeat on which it acts. A range of such fragments was found, from those which augment the effect of colL to those which block it completely. In general, it appears that the ability of a sequence to block the effect of colL depends on both the length and G + C content of the fragment. The sequences which are responsible for the extrusion by the C-type pathway are termed C-type inducing sequences, while sequences which are interposed between the inducing sequence and the inverted repeat, and which may either augment or attenuate the effect, but which cannot function as inducing sequences in isolation, are termed transmitting sequences. The results of these studies are most readily consistent with long range destabilization of DNA structure via telestability effects.     

Evidence for long poly(dA).poly(dT) tracts in D. discoideum DNA at high frequencies and their preferential avoidance of nucleosomal DNA core regions

The eukaryote, Dictyostelium discoideum, has one of the most (A+T) rich genomes studied to date. Isolated nuclear D. discoideum DNA (AX3 strain) was used to qualitatively determine the frequency and length distribution of long (dA).(dT) homopolymer tracts in this genome, in comparison to the less (A+T) rich calf thymus and Schistosoma mansoni DNAs that had few observable long tracts. These experimental data accurately reflect the significantly elevated frequencies of long tracts found computationally within the D. discoideum intron and flanking sequences, but not exons. PCR amplification of long (dA).(dT) homopolymer tract containing sequences was carried out. Then experimental biotinylated (dT)18 probe hybridization to the PCR amplified DNA showed that the long (dA).(dT) homopolymer tracts were enriched in D. discoideum sequences only hundreds of base pair in length, under conditions where no equivalent hybridization was observed to S. mansoni DNA or calf DNA sequences. Similar probe hybridization to DNA isolated following micrococcal nuclease digestion of D. discoideum chromatin demonstrated that long (dA).(dT) homopolymer tracts were more highly enriched in nucleosomal DNA lengths that included the internucleosomal linker as compared to shorter linker free mononucleosomal lengths. This observation is in agreement with the frequency of tract spacing results calculated from GenBank sequence data. These frequency data indicate that adjacent long tracts plus the intervening spacer DNA are found at peak lengths (average 42 bp), exactly characteristic of the internucleosomal spacer region of D. discoideum chromatin and are in sufficient number to be found in nearly half of all nucleosomes. Compared to shuffled tract sequence controls, these lengths of adjacent long tracts plus the intervening spacer DNA were found to be significantly enriched. Lesser enrichments are observed at lengths corresponding to adjacent tracts being separated by nucleosomal core length DNA sequences (145-185 bp). These data strongly suggest that adjacent long tracts occur spaced at selected lengths so as to avoid the central core regions of nucleosomes and instead are found localized within internucleosomal DNA linker and core edge regions in D. discoideum chromatin.     

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA.

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.     

Length Mutations in Human Mitochondrial DNA

By high-resolution, restriction mapping of mitochondrial DNAs purified from 112 human individuals, we have identified 14 length variants caused by small additions and deletions (from about 6 to 14 base pairs in length). Three of the 14 length differences are due to mutations at two locations within the D loop, whereas the remaining 11 occur at seven sites that are probably within other noncoding sequences and at junctions between coding sequences. In five of the nine regions of length polymorphism, there is a sequence of five cytosines in a row, this sequence being comparatively rare in coding DNA. Phylogenetic analysis indicates that, in most of the polymorphic regions, a given length mutation has arisen several times independently in different human lineages. The average rate at which length mutations have been arising and surviving in the human species is estimated to be many times higher for noncoding mtDNA than for noncoding nuclear DNA. The mystery of why vertebrate mtDNA is more prone than nuclear DNA to evolve by point mutation is now compounded by the discovery of a similar bias toward rapid evolution by length mutation.     

Derivation of nearest-neighbor properties from data on nucleic acid oligomers. I. Simple sets of independent sequences and the influence of absent nearest neighbors

The constraints on combinations of nearest neighbors in nucleic acid sequences and the numbers of independent sequences needed to describe nearest-neighbor properties of oligomers and polymers are derived and summarized. It has been pointed out in previous work [D. M. Gray and I. Tinoco, Jr. (1970) Biopolymers, Vol. 9, pp. 223-244; R. F. Goldstein and A. S. Benight (1992) Biopolymers, Vol. 32, pp. 1679-1693] that these constraints restrict the information available from measurements of properties of sequence combinations. The emphasis in this paper is on the properties of oligomer sequences that vary in length, where each nucleotide or base pair at the end of the sequence makes a significant contribution to the measured property by interacting with its boundary of fixed sequence or solvent. In such cases it is not be possible to determine values of properties of individual nearest neighbors, except for the like neighbors [e.g., d(A-A), d(G-G), d(T-T), and d(C-C) nucleotide neighbors in single-stranded DNA or d(A-A)/d(T-T) and d(G-G)/d(C-C) base pair neighbors in double-stranded DNA], solely from measurements of properties of different sequences. Even values for properties of the like neighbors cannot be determined from such oligomeric sequences if the sequences are all of the same length. Nearest-neighbor properties of oligomer sequences that vary in length can be summarized in terms of the values for independent sets of sequences that are nearest neighbors and monomers all with boundaries of the fixed sequence or solvent. Straightforward combinations of the values for the independent sequences will give the values of the property for any dependent sequence, without explicit knowledge of the individual nearest-neighbor values. These considerations have important consequences for the derivation of widely used thermodynamic parameters, as discussed in the following paper.     

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.     

Overlapping redundant septuplets identical with regulatory elements of HIV-1 and SV40.

Overlapping redundant short oligomers in DNA sequences of retroviruses and papovaviruses have been identified. For each sequence, a search procedure determines the 5% short oligomers of the same length with the highest ratios of observed to expected occurrences based on singlet composition of the sequence. These short oligomers are referred to as compositionally-assessed redundant sequence elements (COARSEs). A pair of COARSEs overlapping by at least one base is considered to be a COARSE overlap. Most COARSE overlaps of the 7th order (overlapping septuplets) are found in long terminal repeats of retroviruses and in the regulatory control regions of papovaviruses SV40, BK and JC. Many of the 7th order COARSE overlaps in HIV-1 and SV40 are identical with regulatory elements determined experimentally. On the contrary, very few of the most frequently occurring oligomer overlaps, which are defined differently from COARSE overlaps, are present in the regulatory regions of retroviruses and papovaviruses. Examining DNA sequences of other genomes by the COARSE overlap method may identify putative regulatory regions.     

Base pair opening in three DNA-unwinding elements

DNA-unwinding elements are specific base sequences that are located in the origin of DNA replication where they provide the start point for strand separation and unwinding of the DNA double helix. In the present work we have obtained the first characterization of the opening of individual base pairs in DNA-unwinding elements. The three DNA molecules investigated reproduce the 13-mer DNA-unwinding elements present in the Escherichia coli chromosome. The base sequences of the three 13-mers are conserved in the origins of replication of enteric bacterial chromosomes. The exchange of imino protons with solvent protons was measured for each DNA as a function of the concentration of exchange catalyst using nuclear magnetic resonance spectroscopy. The exchange rates provided the rates and the equilibrium constants for opening of individual base pairs in each DNA at 20 degrees C. The results reveal that the kinetics and energetics of the opening reactions for AT/TA base pairs are different in the three DNA-unwinding elements due to long range effects of the base sequence. These differences encompass the AT/TA base pairs that are conserved in various bacterial genomes. Furthermore, a qualitative correlation is observed between the kinetics and energetics of opening of AT/TA base pairs and the location of the corresponding DNA-unwinding element in the origin of DNA replication.