DNA Architecture,Deformability, and Nucleosome Positioning

Abstract

The positioning of DNA on nucleosomes is critical to both the organization and expression of the genetic message. Here we focus on DNA conformational signals found in the growing library of known high-resolution core-particle structures and the ways in which these features may contribute to the positioning of nucleosomes on specific DNA sequences. We survey the chemical composition of the protein-DNA assemblies and extract features along the DNA superhelical pathway—the minor-groove width and the deformations of successive base pairs—determined with reasonable accuracy in the structures. We also examine the extent to which the various nucleosome core-particle structures accommodate the observed settings of the crystallized sequences and the known positioning of the high-affinity synthetic ‘601’ sequence on DNA. We ‘thread’ these sequences on the different structural templates and estimate the cost of each setting with knowledge-based potentials that reflect the conformational properties of the DNA base-pair steps in other high-resolution protein-bound complexes.     

Sequence-directed Mapping of Nucleosome Positions

Abstract

A nucleosome DNA sequence probe is designed that combines recently derived RR/YY counter-phase and AA/TT in-phase periodical patterns. A simple nucleosome mapping procedure is introduced for prediction of the nucleosome positions in the sequence of interest, to serve as a guide for experimental studies of the chromatin structure.     

Nucleosome DNA Bendability Matrix (C. elegans)

Abstract

An original signal extraction procedure is applied to database of 146 base nucleosome core DNA sequences from C. elegans (S. M. Johnson et al. Genome Research 16, 1505–1516, 2006). The positional preferences of various dinucleotides within the 10.4 base nucleosome DNA repeat are calculated, resulting in derivation of the nucleosome DNA bendability matrix of 16x10 elements. A simplified one-line presentation of the matrix (“consensus” repeat) is…A(TTTCCGGAAA)T…. All 6 chromosomes of C. elegans conform to the bendability pattern. The strongest affinity to their respective positions is displayed by dinucleotides AT and CG, separated within the repeat by 5 bases. The derived pattern makes a basis for sequence-directed mapping of nucleosome positions in the genome of C. elegans. As the first complete matrix of bendability available the pattern may serve for iterative calculations of the species-specific matrices of bendability applicable to other genomic sequences.     

Sequence Signatures of Nucleosome Positioning in Caenorhabditis elegans

Our recent investigation in the protist Trichomonas vaginalis suggested a DNA sequence periodicity with a unit length of 120.9 nt, which represents a sequence signature for nucleosome positioning. We now extended our observation in higher eukaryotes and identified a similar periodicity of 175 nt in length in Caenorhabditis elegans. In the process of defining the sequence compositional characteristics, we found that the 10.5-nt periodicity, the sequence signature of DNA double helix, may not be sufficient for cross-nucleosome positioning but provides essential guiding rails to facilitate positioning. We further dissected nucleosome-protected sequences and identified a strong positive purine (AG) gradient from the 5′-end to the 3′-end, and also learnt that the nucleosome-enriched regions are GC-rich as compared to the nucleosome-free sequences as purine content is positively correlated with GC content. Sequence characterization allowed us to develop a hidden Markov model (HMM) algorithm for decoding nucleosome positioning computationally, and based on a set of training data from the fifth chromosome of C. elegans, our algorithm predicted 60%-70% of the well-positioned nucleosomes, which is 15%-20% higher than random positioning. We concluded that nucleosomes are not randomly positioned on DNA sequences and yet bind to different genome regions with variable stability, well-positioned nucleosomes leave sequence signatures on DNA, and statistical positioning of nucleosomes across genome can be decoded computationally based on these sequence signatures.     

The Molecular Evolution of Nucleosome Positioning Through Sequence-Dependent Deformation of the DNA Polymer

Abstract

The computational prediction of nucleosome positioning from DNA sequence now allows for in silico investigation of the molecular evolution of biophysical properties of the DNA molecule responsible for primary chromatin organization in the genome. To discern what signal components driving nucleosome positioning in the yeast genome are potentially targeted by natural selection, we compare the performance of various models predictive of nucleosome positioning within the context of a simple statistical test, the repositioned mutation test. We demonstrate that while nucleosome occupancy is driven largely by translational exclusion in response to AT content, there is also a strong signature of evolutionary conservation of regular patterns within nucleosomal DNA sequence related to the structural organization of the nucleosome core (e.g., 10-bp dinucleotide periodicity). We also use computer simulations to investigate hypothetical coding and regulatory constraints on the ability of sequence properties affecting nucleosome formation to adaptively evolve. Our results demonstrate that natural selection may act independently on different DNA sequence properties responsible for local chromatin organization. Furthermore, at least with respect to the deformation energy of the DNA molecule in the nucleosome, the presence of the genetic code has greatly restricted the ability of sequences to evolve the dynamic nucleosome organization typically observed in promoter regions.     

Positioning of a Nucleosome on Mouse Satellite DNA Inserted into a Yeast Plasmid Is Determined by Its DNA Sequence and an Adjacent Nucleosome

It has earlier been shown that multiple positioning of nucleosomes on mouse satellite DNA is determined by its nucleotide sequence. To clarify whether other factors, such as boundary ones, can affect the positionings, we modified the environment of satellite DNA monomer by inserting it into a yeast plasmid between inducible GalCyc promoter and a structural region of the yeast FLP gene. We have revealed that the positions of nucleosomes on satellite DNA are identical to those detected upon reconstruction in vitro. The positioning signal (GAAAAA sequence) of satellite DNA governs nucleosome location at the adjacent nucleotide sequence as well. Upon promoter induction the nucleosome, translationally positioned on the GalCyc promoter, transfers to the satellite DNA and its location follows the positioning signal of the latter. Thus, the alternatives of positioning of a nucleosome on satellite DNA are controlled by its nucleotide sequence, though the choice of one of them is determined by the adjacent nucleosome.     

Base pair stacking in nucleosome DNA and bendability sequence pattern

DNA deformation in the nucleosome involves partial unstacking between bases and base pairs. By adjusting orientations of different base-pair stacks relative to the histone octamer surface, the optimal set of stacks and their positions is derived, resulting in a sequence pattern, theoretically best suitable for nucleosome DNA. The sequence is very much consistent with available experimental data, thus, suggesting a common eukaryotic nucleosome DNA bendability sequence pattern based exclusively on the very basics of DNA.     

Correlation between the flexibility and periodic dinucleotide patterns in yeast nucleosomal DNA sequences

Nucleosome formation and positioning, which play important roles in a number of biological processes, are thought to be related to the distinctive periodic dinucleotide patterns observed in the DNA sequence wrapped around the protein octamer. Previous research shows that flexibility is a key structural property of a nucleosomal DNA sequence. However, the relationship between the flexibility and the periodic dinucleotide patterns has received little attention in research in the past. In this study, we propose the use of three different models to measure the flexibility of yeast DNA sequences. Although the three models involve different parameters, they deliver consistent results showing that yeast nucleosomal DNA sequences are more flexible than non-nucleosomal ones. In contrast to random flexibility values along non-nucleosomal DNA sequences, the flexibility of nucleosomal DNA sequences shows a clear periodicity of 10.14 base pairs, which is consistent with the periodicity of dinucleotide distributions. We also demonstrate that there is a strong relationship between the peak positions of the flexibility and the dinucleotide frequencies. Correlation between the flexibility and the dinucleotide patterns of CA/TG, CG, GC, GG/CC, AG/CT, AC/GT and GA/TC are positive with an average value of 0.5946. The highest correlation is shown by CA/TG with a value of 0.7438 and the lowest correlation is shown by AA/TT with a value of −0.7424. The source codes and data sets are available for downloading on http://www.hy8.com/bioinformatics.htm.     

High-Resolution Mapping of Sequence-Directed Nucleosome Positioning on Genomic DNA

We have mapped in vitro nucleosome positioning on the sheep β-lactoglobulin gene using high-throughput sequencing to characterise the DNA sequences recovered from reconstituted nucleosomes. This methodology surpasses previous approaches for coverage, accuracy and resolution and, most importantly, offers a simple yet rapid and relatively inexpensive method to characterise genomic DNA sequences in terms of nucleosome positioning capacity. We demonstrate an unambiguous correspondence between in vitro and in vivo nucleosome positioning around the promoter of the gene; identify discrete, sequence-specific nucleosomal structures above the level of the canonical core particle—a feature that has implications for regulatory protein access and higher-order chromatin packing; and reveal new insights into the involvement of periodically organised dinucleotide sequence motifs of the type GG and CC and not AA and TT, as determinants of nucleosome positioning—an observation that supports the idea that the core histone octamer can exploit different patterns of sequence organisation, or structural potential, in the DNA to bring about nucleosome positioning.     

Assessing Diversity of DNA Structure-Related Sequence Features in Prokaryotic Genomes

Prokaryotic genomes are diverse in terms of their nucleotide and oligonucleotide composition as well as presence of various sequence features that can affect physical properties of the DNA molecule. We present a survey of local sequence patterns which have a potential to promote non-canonical DNA conformations (i.e. different from standard B-DNA double helix) and interpret the results in terms of relationships with organisms'' habitats, phylogenetic classifications, and other characteristics. Our present work differs from earlier similar surveys not only by investigating a wider range of sequence patterns in a large number of genomes but also by using a more realistic null model to assess significant deviations. Our results show that simple sequence repeats and Z-DNA-promoting patterns are generally suppressed in prokaryotic genomes, whereas palindromes and inverted repeats are over-represented. Representation of patterns that promote Z-DNA and intrinsic DNA curvature increases with increasing optimal growth temperature (OGT), and decreases with increasing oxygen requirement. Additionally, representations of close direct repeats, palindromes and inverted repeats exhibit clear negative trends with increasing OGT. The observed relationships with environmental characteristics, particularly OGT, suggest possible evolutionary scenarios of structural adaptation of DNA to particular environmental niches.     

Molecular Dynamics Simulations of a Nucleosome and Free DNA

Abstract

All atom molecular dynamics simulations (10ns) of a nucleosome and of its 146 basepairs of DNA free in solution have been conducted. DNA helical parameters (Roll, Tilt, Twist, Shift, Slide, Rise) were extracted from each trajectory to compare the conformation, effective force constants, persistence length measures, and fluctuations of nucleosomal DNA to free DNA. The conformation of DNA in the nucleosome, as determined by helical parameters, is found to be largely within the range of thermally accessible values obtained for free DNA. DNA is found to be less flexible on the nucleosome than when free in solution, however such measures are length scale dependent. A method for disassembling and reconstructing the conformation and dynamics of the nucleosome using Fourier analysis is presented. Long length variations in the conformation of nucleosomal DNA are identified other than those associated with helix repeat. These variations are required to create a proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns of the nucleosome's superhelix. Reconstruction of free DNA using selected long wavelength variations in conformation can produce either a left-handed or a right-handed superhelix. The long wavelength variations suggest 146 basepairs is a natural length of DNA to wrap around the histone core.     

Structural Mechanics of DNA Wrapping in the Nucleosome

Experimental X-ray crystal structures and a database of calculated structural parameters of DNA octamers were used in combination to analyse the mechanics of DNA bending in the nucleosome core complex. The 1kx5 X-ray crystal structure of the nucleosome core complex was used to determine the relationship between local structure at the base-step level and the global superhelical conformation observed for nucleosome-bound DNA. The superhelix is characterised by a large curvature (597°) in one plane and very little curvature (10°) in the orthogonal plane. Analysis of the curvature at the level of 10-step segments shows that there is a uniform curvature of 30° per helical turn throughout most of the structure but that there are two sharper kinks of 50° at ± 2 helical turns from the central dyad base pair. The curvature is due almost entirely to the base-step parameter roll. There are large periodic variations in roll, which are in phase with the helical twist and account for 500° of the total curvature. Although variations in the other base-step parameters perturb the local path of the DNA, they make minimal contributions to the total curvature. This implies that DNA bending in the nucleosome is achieved using the roll-slide-twist degree of freedom previously identified as the major degree of freedom in naked DNA oligomers. The energetics of bending into a nucleosome-bound conformation were therefore analysed using a database of structural parameters that we have previously developed for naked DNA oligomers. The minimum energy roll, the roll flexibility force constant and the maximum and minimum accessible roll values were obtained for each base step in the relevant octanucleotide context to account for the effects of conformational coupling that vary with sequence context. The distribution of base-step roll values and corresponding strain energy required to bend DNA into the nucleosome-bound conformation defined by the 1kx5 structure were obtained by applying a constant bending moment. When a single bending moment was applied to the entire sequence, the local details of the calculated structure did not match the experiment. However, when local 10-step bending moments were applied separately, the calculated structure showed excellent agreement with experiment. This implies that the protein applies variable bending forces along the DNA to maintain the superhelical path required for nucleosome wrapping. In particular, the 50° kinks are constraints imposed by the protein rather than a feature of the 1kx5 DNA sequence. The kinks coincide with a relatively flexible region of the sequence, and this is probably a prerequisite for high-affinity nucleosome binding, but the bending strain energy is significantly higher at these points than for the rest of the sequence. In the most rigid regions of the sequence, a higher strain energy is also required to achieve the standard 30° curvature per helical turn. We conclude that matching of the DNA sequence to the local roll periodicity required to achieve bending, together with the increased flexibility required at the kinks, determines the sequence selectivity of DNA wrapping in the nucleosome.     

Dynamics of Nucleosome Assembly and Effects of DNA Methylation

The nucleosome is the fundamental packing unit of the eukaryotic genome, and CpG methylation is an epigenetic modification associated with gene repression and silencing. We investigated nucleosome assembly mediated by histone chaperone Nap1 and the effects of CpG methylation based on three-color single molecule FRET measurements, which enabled direct monitoring of histone binding in the context of DNA wrapping. According to our observation, (H3-H4)₂ tetramer incorporation must precede H2A-H2B dimer binding, which is independent of DNA termini wrapping. Upon CpG methylation, (H3-H4)₂ tetramer incorporation and DNA termini wrapping are facilitated, whereas proper incorporation of H2A-H2B dimers is inhibited. We suggest that these changes are due to rigidified DNA and increased random binding of histones to DNA. According to the results, CpG methylation expedites nucleosome assembly in the presence of abundant DNA and histones, which may help facilitate gene packaging in chromatin. The results also indicate that the slowest steps in nucleosome assembly are DNA termini wrapping and tetramer positioning, both of which are affected heavily by changes in the physical properties of DNA.     

The DNA Sequence-dependence of Nucleosome Positioning in vivo and in vitro

Abstract

The contribution of histone-DNA interactions to nucleosome positioning in vivo is currently a matter of debate. We argue here that certain nucleosome positions, often in promoter regions, in yeast may be, at least in part, specified by the DNA sequence. In contrast other positions may be poorly specified. Positioning thus has both statistical and DNA-determined components. We further argue that the relative affinity of the octamer for different DNA sequences can vary and therefore the interaction of histones with the DNA is a ‘tunable’ property.     

Sequence Structure of Hidden 10.4-base Repeat in the Nucleosomes of C. elegans

Abstract

By measuring prevailing distances between YY, YR, RR, and RY dinucleotides in the large database of the nucleosome DNA fragments from C. elegans, the consensus sequence structure of the nucleosome DNA repeat of C. elegans was reconstructed: (YYYYYRRRRR)_n. An actual period was estimated to be 10.4 bases. The pattern is fully consistent with the nucleosome DNA patterns of other eukaryotes, as established earlier, and, thus, the YYYYYRRRRR repeat can be considered as consensus nucleosome DNA sequence repeat across eukaryotic species. Similar distance analysis for [A, T] dinucleotides suggested the related pattern (TTTYTARAAA)_n where the TT and AA dinucleotides display rather out of phase behavior, contrary to the “AA or TT” in-phase periodicity, considered in some publications. A weak 5-base periodicity in the distribution of TA dinucleotides was detected.     

Plasticity and Epigenetic Inheritance of Centromere-specific Histone H3 (CENP-A)-containing Nucleosome Positioning in the Fission Yeast

Nucleosomes containing the specific histone H3 variant CENP-A mark the centromere locus on each chromatin and initiate kinetochore assembly. For the common type of regional centromeres, little is known in molecular detail of centromeric chromatin organization, its propagation through cell division, and how distinct organization patterns may facilitate kinetochore assembly. Here, we show that in the fission yeast S. pombe, a relatively small number of CENP-A/Cnp1 nucleosomes are found within the centromeric core and that their positioning relative to underlying DNA varies among genetically homogenous cells. Consistent with the flexible positioning of Cnp1 nucleosomes, a large portion of the endogenous centromere is dispensable for its essential activity in mediating chromosome segregation. We present biochemical evidence that Cnp1 occupancy directly correlates with silencing of the underlying reporter genes. Furthermore, using a newly developed pedigree analysis assay, we demonstrated the epigenetic inheritance of Cnp1 positioning and quantified the rate of occasional repositioning of Cnp1 nucleosomes throughout cell generations. Together, our results reveal the plasticity and the epigenetically inheritable nature of centromeric chromatin organization.     

DNA Curvature of the Tobacco GRS Repetitive Sequence Family and its Relation to Nucleosome Positioning

Abstract

Recently, a highly repetitive DNA sequence family (GRS) from tobacco was described in our laboratory. These sequences were found to be localized predominantly in the pericen-tromeric heterochromatin of tobacco chromosomes. To test the hypothesis that these sequences play an important role in the formation of heterochromatin, we investigated the DNA curvature of the GRS sequences and its possible impact to the chromatin structure at these loci. Application of the nearest-neighbour wedge model of intrinsic DNA curvature for the GRS1 family member predicted two loci of curvature: a major bend at the 5′ end of the sequence and a minor bend of opposite direction at the centre of the GRS1. The presence of the major and the minor loci of DNA curvature was studied experimentally using permutation analysis and site-directed mutagenesis. The experimental results were consistent with the computer predictions. We gave evidence that the described DNA curvature is also present in the entire GRS family. Genomic statistical sequencing showed the conservation of the major bend sequence determinants in the members of the GRS family. To investigate the chromatin structure at the GRS sequences, we determined the nucleosome positioning in vivo at these sequences using thermal cycle primer extension. A relation between the curvature pattern and the histone octamer position was observed: the major bend is excluded from the nucleosome surface to the linker region, while the minor bend is distributed along the core DNA The suggestion is made that the sequences in the minor locus of curvature define the rotational setting of the nucleosome, and a possible role of the major bend as a factor, which defines the translational setting, is discussed.