Nucleosome structural features and intrinsic properties of the TATAAACGCC repeat sequence

An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data

Curved DNA fragments are often found near functionally important sites such as promoters and origins of replication, and hence sequence-dependent DNA curvature prediction is of great utility in genomics and bioinformatics. In light of this, an assessment of three different dinucleotide step parameters (based on gel retardation as well as crystal structure data) is carried out. These parameters (BMHT, LB and CS) are evaluated quantitatively for their ability to predict correctly the experimental results of a large set of nucleic acid sequences containing A-tracts as well as GC-rich motifs. This set contained around 40 synthetic as well as natural sequences whose solution properties have been well characterized experimentally. All three models could account reasonably well for curvature in the various DNA sequences. The CS model, where dinucleotide parameters are calculated from crystal structure data, consistently shows slightly better correlation with experimental data. Our simple analysis also indicates that presently available trinucleotide parameters fail to predict curvature in some of the well-characterized sequences. The study shows that the dinucleotide parameters with some further refinement can be used to predict sequence-dependent curvature correctly in genomic sequences.

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.

The human IL-2 gene promoter can assemble a positioned nucleosome that becomes remodeled upon T cell activation

Controlled production of the cytokine IL-2 plays a key role in the mammalian immune system. Expression from the gene is tightly regulated with no detectable expression in resting T cells and a strong induction following T cell activation. The IL-2 proximal promoter (+1 to -300) contains many well-defined transcriptional activation elements that respond to T cell stimulation. To determine the role of chromatin structure in the regulation of interleukin-2 gene transcription, nucleosome assembly across the IL-2 promoter region was examined using in vitro chromatin reconstitution assays. The IL-2 promoter assembles a nucleosome that is both translationally and rotationally positioned, spanning some of the major functional control elements. The binding of transcription factors to these elements, with the exception of the architectural protein HMGA1, was occluded by the presence of the nucleosome. Analysis of the chromatin architecture of the IL-2 gene in Jurkat T cells provided evidence for the presence of a similarly positioned nucleosome in vivo. The region encompassed by this nucleosome becomes remodeled following activation of Jurkat T cells. These observations suggest that the presence of a positioned nucleosome across the IL-2 proximal promoter may play an important role in maintaining an inactive gene in resting T cells and that remodeling of this nucleosome is important for gene activation.

The influence of the thymine C5 methyl group on spontaneous base pair breathing in DNA

Sequences of four or more AT base pairs without a 5'-TA-3' step, so-called A-tracts, influence the global properties of DNA by causing curvature of the helix axis if phased with the helical repeat and also influence nucleosome packaging. Hence it is interesting to understand this phenomenon on the molecular level, and numerous studies have been devoted to investigations of dynamical and structural features of A-tract DNA. It was early observed that anomalously slow base pair-opening kinetics were a striking physical property unique to DNA A-tracts (Leroy, J. L., Charretier, E., Kochoyan, M., and Gueron, M. (1988) Biochemistry 27, 8894-8898). Furthermore, a strong correlation between DNA curvature and anomalously slow base pair-opening dynamics was found. In the present work it is shown, using imino proton exchange measurements by NMR spectroscopy that the main contribution to the dampening of the base pair-opening fluctuations in A-tracts comes from the C5 methylation of the thymine base. Because the methyl group has been shown to have a very limited effect on the DNA curvature as well as the structure of the DNA helix, the thymine C5 methyl group stabilizes the helix directly. Empirical potential energy calculations show that methylation of the tract improves the stacking energy of a base pair with its neighbors in the tract by 3-4 kcal/mol.

Both DNA and histone fold sequences contribute to archaeal nucleosome stability

The roles and interdependence of DNA sequence and archaeal histone fold structure in determining archaeal nucleosome stability and positioning have been determined and quantitated. The presence of four tandem copies of TTTAAAGCCG in the polylinker region of pLITMUS28 resulted in a DNA molecule with increased affinity (DeltaDeltaG of approximately 700 cal mol(-1)) for the archaeal histone HMfB relative to the polylinker sequence, and the dominant, quantitative contribution of the helical repeats of the dinucleotide TA to this increased affinity has been established. The rotational and translational positioning of archaeal nucleosomes assembled on the (TTTAAAGCCG)(4) sequence and on DNA molecules selectively incorporated into archaeal nucleosomes by HMfB have been determined. Alternating A/T- and G/C-rich regions were located where the minor and major grooves, respectively, sequentially faced the archaeal nucleosome core, and identical positioning results were obtained using HMfA, a closely related archaeal histone also from Methanothermus fervidus. However, HMfA did not have similarly high affinities for the HMfB-selected DNA molecules, and domain-swap experiments have shown that this difference in affinity is determined by residue differences in the C-terminal region of alpha-helix 3 of the histone fold, a region that is not expected to directly interact with DNA. Rather this region is thought to participate in forming the histone dimer:dimer interface at the center of an archaeal nucleosome histone tetramer core. If differences in this interface do result in archaeal histone cores with different sequence preferences, then the assembly of alternative archaeal nucleosome tetramer cores could provide an unanticipated and novel structural mechanism to regulate gene expression.

Magnesium increases the curvature of duplex DNA that contains dA tracts

Distinct structural features of DNA, such as the curvature of dA tracts, are important in the recognition, packaging, and regulation of DNA. Physiologically relevant concentrations of magnesium have been found to enhance the curvature of dA tract DNAs, as monitored by solution-state NMR, indicating that the structure of DNA depends on the cations present in solution. A model is presented which accounts for the sequence-dependent effects of magnesium on DNA curvature as well as for the previously known sequence-independent effect on DNA flexibility.

Nucleosomes and the chromatin fiber

During the past year and a half, significant progress has been made in understanding the structure and dynamics of nucleosomes and the chromatin fiber, the mechanism of action of the core histone amino termini, the structure and function of histone variants, and the function of linker histones in the chromatin fiber.

Intrinsically bent DNA in the promoter regions of the yeast GAAL1-10 and GAL80 genes

Circular permutation analysis has detected fairly strong sites of intrinsic DNA bending on the promoter regions of the yeast GAL1-10 and GAL80 genes. These bends lie in functionally suggestive locations. On the promoter of the GAL1-10 structural genes, strong bends bracket nucleosome B, which lies between the UAS(G) and the GAL1 TATA. These intrinsic bends could help position nucleosome B. Nucleosome B plus two other promoter nucleosomes protect the TATA and start site elements in the inactive state of expression but are completely disrupted (removed) when GAL1-10 expression is induced. The strongest intrinsic bend ( approximately 70 degrees ) lies at the downstream edge of nucleosome B; this places it approximately 30 base pairs upstream of the GAL1 TATA, a position that could allow it to be involved in GAL1 activation in several ways, including the recruitment of a yeast HMG protein that is required for the normally robust level of GAL1 expression in the induced state (Paull, T., Carey, M., and Johnson, R. (1996) Genes Dev. 10, 2769-2781). On the regulatory gene GAL80, the single bend lies in the non-nucleosomal hypersensitive region, between a GAL80-specific far upstream promoter element and the more gene-proximal promoter elements. GAL80 promoter region nucleosomes contain no intrinsically bent DNA.

Global structure and mechanical properties of a 10-bp nucleosome positioning motif

The method of DNA cyclization kinetics reveals special properties of the TATAAACGCC sequence motif found in DNA sequences that have high affinity for core histones. Replacement of 30 bp of generic DNA by three 10-bp repeats of the motif in small cyclization constructs increases cyclization rates by two orders of magnitude. We document a 13 degrees bend in the motif and characterize the direction of curvature. The bending force constant is smaller by nearly 2-fold and there is a 35% decrease in the twist modulus, relative to generic DNA. These features are the likely source of the high affinity for bending around core histones to form nucleosomes. Our results establish a protocol for determination of the ensemble-averaged global solution structure and mechanical properties of any approximately 10-bp DNA sequence element of interest, providing information complementary to that from NMR and crystallographic structural studies.

Archaeal histone selection of nucleosome positioning sequences and the procaryotic origin of histone-dependent genome evolution

Archaeal histones and the eucaryal (eucaryotic) nucleosome core histones have almost identical histone folds. Here, we show that DNA molecules selectively incorporated by rHMfB (recombinant archaeal histone B from Methanothermus fervidus) into archaeal nucleosomes from a mixture of approximately 10(14) random sequence molecules contain sequence motifs shown previously to direct eucaryal nucleosome positioning. The dinucleotides GC, AA (=TT) and TA are repeated at approximately 10 bp intervals, with the GC harmonic displaced approximately 5 bp from the AA and TA harmonics [(GCN(3)AA or TA)(n)]. AT and CG were not strongly selected, indicating that TA not equalAT and GC not equalCG in terms of facilitating archaeal nucleosome assembly. The selected molecules have affinities for rHMfB ranging from approximately 9 to 18-fold higher than the level of affinity of the starting population, and direct the positioned assembly of archaeal nucleosomes. Fourier-transform analyses have revealed that AA dinucleotides are much enriched at approximately 10. 1 bp intervals, the helical repeat of DNA wrapped around a nucleosome, in the genomes of Eucarya and the histone-containing Euryarchaeota, but not in the genomes of Bacteria and Crenarchaeota, procaryotes that do not have histones. Facilitating histone packaging of genomic DNA has apparently therefore imposed constraints on genome sequence evolution, and since archaeal histones have no structure in addition to the histone fold, these constraints must result predominantly from histone fold-DNA contacts. Based on the three-domain universal phylogeny, histones and histone-dependent genome sequence evolution most likely evolved after the bacterial-archaeal divergence but before the archaeal-eucaryal divergence, and were subsequently lost in the Crenarchaeota. However, with lateral gene transfer, the first histone fold could alternatively have evolved after the archaeal-eucaryal divergence, early in either the euryarchaeal or eucaryal lineages.

Chromatin organization and transcriptional control of gene expression in Drosophila

It is increasingly clear that the packaging of DNA in nucleosome arrays serves not only to constrain the genome within the nucleus, but also to encode information concerning the activity state of the gene. Packaging limits the accessibility of many regulatory DNA sequence elements and is functionally significant in the control of transcription, replication, repair and recombination. Here, we review studies of the heat-shock genes, illustrating the formation of a specific nucleosome array at an activatable promoter, and describe present information on the roles of DNA-binding factors and energy-dependent chromatin remodeling machines in facilitating assembly of an appropriate structure. Epigenetic maintenance of the activity state within large domains appears to be a key mechanism in regulating homeotic genes during development; recent advances indicate that chromatin structural organization is a critical parameter. The ability to utilize genetic, biochemical and cytological approaches makes Drosophila an ideal organism for studies of the role of chromatin structure in the regulation of gene expression.

Signals for TBP/TATA box recognition

The TATA box-binding protein (TBP) recognizes its target sites (TATA boxes) by indirectly reading the DNA sequence through its conformation effects (indirect readout). Here, we explore the molecular mechanisms underlying indirect readout of TATA boxes by TBP by studying the binding of TBP to adenovirus major late promoter (AdMLP) sequence variants, including alterations inside as well as in the sequences flanking the TATA box. We measure here the dissociation kinetics of complexes of TBP with AdMLP targets and, by phase-sensitive assay, the intrinsic bending in the TATA box sequences as well as the bending of the same sequence induced by TBP binding. In these experiments we observe a correlation of the kinetic stability to sequence changes within the TATA recognition elements. Comparison of the kinetic data with structural properties of TATA boxes in known crystalline TBP/TATA box complexes reveals several "signals" for TATA box recognition, which are both on the single base-pair level, as well as larger DNA tracts within the TATA recognition element. The DNA bending induced by TBP on its binding sites is not correlated to the stability of TBP/TATA box complexes. Moreover, we observe a significant influence on the kinetic stability of alteration in the region flanking the TATA box. This effect is limited however to target sites with alternating TA sequences, whereas the AdMLP target, containing an A tract, is not influenced by these changes.

DNA sequence-dependent contributions of core histone tails to nucleosome stability: differential effects of acetylation and proteolytic tail removal

Modulation of nucleosome stability in chromatin plays an important role in eukaryotic gene expression. The core histone N-terminal tail domains are believed to modulate the stability of wrapping nucleosomal DNA and the stability of the chromatin filament. We analyzed the contribution of the tail domains to the stability of nucleosomes containing selected DNA sequences that are intrinsically straight, curved, flexible, or inflexible. We find that the presence of the histone tail domains stabilizes nucleosomes containing DNA sequences that are intrinsically straight or curved. However, the tails do not significantly contribute to the free energy of nucleosome formation with flexible DNA. Interestingly, hyperacetylation of the core histone tail domains does not recapitulate the effect of tail removal by limited proteolysis with regard to nucleosome stability. We find that acetylation of the tails has the same minor effect on nucleosome stability for all the selected DNA sequences. A comparison of histone partitioning between long donor chromatin, acceptor DNA, and free histones in solution shows that the core histone tails mediate internucleosomal interactions within an H1-depleted chromatin fiber amounting to an average free energy of about 1 kcal/mol. Thus, such interactions would be significant with regard to the free energies of sequence-dependent nucleosome positioning. Last, we analyzed the contribution of the H2A/H2B dimers to nucleosome stability. We find that the intact nucleosome is stabilized by 900 cal/mol by the presence of the dimers regardless of sequence. The biological implications of these observations are discussed.