Nucleosomal DNA Sequences Research Articles

Unveiling Nucleosome Dynamics: A Comparative Study Using All-Atom and Coarse-Grained Simulations Enhanced by Principal Component Analysis.

This study investigates nucleosome dynamics using both all-atom and coarse-grained (CG) molecular dynamics simulations, focusing on the SIRAH force field. Simulations are performed for two nucleosomal DNA sequences, ASP and Widom-601, over six microseconds at physiological salt concentrations. Comparative analysis of structural parameters, such as groove widths and base pair geometries, reveals good agreement between atomistic and CG models, though CG simulations exhibit broader conformational sampling and greater breathing motion of DNA ends. Principal component analysis (PCA) is applied to DNA structural parameters, revealing multiple free energy minima, especially in CG simulations. These findings highlight the potential of the SIRAH CG force field for studying large-scale nucleosome dynamics, offering insights into DNA repositioning and sequence-dependent behavior.

Structural basis for p53 binding to its nucleosomal target DNA sequence.

The tumor suppressor p53 functions as a pioneer transcription factor that binds a nucleosomal target DNA sequence. However, the mechanism by which p53 binds to its target DNA in the nucleosome remains elusive. Here we report the cryo-electron microscopy structures of the p53 DNA-binding domain and the full-length p53 protein complexed with a nucleosome containing the 20 base-pair target DNA sequence of p53 (p53BS). In the p53-nucleosome structures, the p53 DNA-binding domain forms a tetramer and specifically binds to the p53BS DNA, located near the entry/exit region of the nucleosome. The nucleosomal position of the p53BS DNA is within the genomic p21 promoter region. The p53 binding peels the DNA from the histone surface, and drastically changes the DNA path around the p53BS on the nucleosome. The C-terminal domain of p53 also binds to the DNA around the center and linker DNA regions of the nucleosome, as revealed by hydroxyl radical footprinting. These results provide important structural information for understanding the mechanism by which p53 binds the nucleosome and changes the chromatin structure for gene activation.

Galaxy Dnpatterntools for Computational Analysis of Nucleosome Positioning Sequence Patterns.

Nucleosomes are basic units of DNA packing in eukaryotes. Their structure is well conserved from yeast to human and consists of the histone octamer core and 147 bp DNA wrapped around it. Nucleosomes are bound to a majority of the eukaryotic genomic DNA, including its regulatory regions. Hence, they also play a major role in gene regulation. For the latter, their precise positioning on DNA is essential. In the present paper, we describe Galaxy dnpatterntools—software package for nucleosome DNA sequence analysis and mapping. This software will be useful for computational biologists practitioners to conduct more profound studies of gene regulatory mechanisms.

The biogenesis and function of nucleosome arrays

Numerous chromatin remodeling enzymes position nucleosomes in eukaryotic cells. Aside from these factors, transcription, DNA sequence, and statistical positioning of nucleosomes also shape the nucleosome landscape. The precise contributions of these processes remain unclear due to their functional redundancy in vivo. By incisive genome engineering, we radically decreased their redundancy in Saccharomyces cerevisiae. The transcriptional machinery strongly disrupts evenly spaced nucleosomes. Proper nucleosome density and DNA sequence are critical for their biogenesis. The INO80 remodeling complex helps space nucleosomes in vivo and positions the first nucleosome over genes in an H2A.Z-independent fashion. INO80 requires its Arp8 subunit but unexpectedly not the Nhp10 module for spacing. Cells with irregularly spaced nucleosomes suffer from genotoxic stress including DNA damage, recombination and transpositions. We derive a model of the biogenesis of the nucleosome landscape and suggest that it evolved not only to regulate but also to protect the genome.

Ramified rolling circle amplification for synthesis of nucleosomal DNA sequences

Nucleosomes are a crucial platform for the recruitment and assembly of protein complexes that process the DNA. Mechanistic and structural in vitro studies typically rely on recombinant nucleosomes that are reconstituted using artificial, strong-positioning DNA sequences. To facilitate such studies on native, genomic nucleosomes, there is a need for methods to produce any desired DNA sequence in an efficient manner. The current methods either do not offer much flexibility in choice of sequence or are less efficient in yield and labor. Here, we show that ramified rolling circle amplification (RCA) can be used to produce milligram amounts of a genomic nucleosomal DNA fragment in a scalable, one-pot reaction overnight. The protocol is efficient and flexible in choice of DNA sequence. It yields 10-fold more product than PCR, and rivals production using plasmids. We demonstrate the approach by producing the genomic DNA from the human LIN28B locus and show that it forms functional nucleosomes capable of binding pioneer transcription factor Oct4.

Sequence-dependent nucleosome sliding in rotation-coupled and uncoupled modes revealed by molecular simulations.

While nucleosome positioning on eukaryotic genome play important roles for genetic regulation, molecular mechanisms of nucleosome positioning and sliding along DNA are not well understood. Here we investigated thermally-activated spontaneous nucleosome sliding mechanisms developing and applying a coarse-grained molecular simulation method that incorporates both long-range electrostatic and short-range hydrogen-bond interactions between histone octamer and DNA. The simulations revealed two distinct sliding modes depending on the nucleosomal DNA sequence. A uniform DNA sequence showed frequent sliding with one base pair step in a rotation-coupled manner, akin to screw-like motions. On the contrary, a strong positioning sequence, the so-called 601 sequence, exhibits rare, abrupt transitions of five and ten base pair steps without rotation. Moreover, we evaluated the importance of hydrogen bond interactions on the sliding mode, finding that strong and weak bonds favor respectively the rotation-coupled and -uncoupled sliding movements.

Hypomethylated domain-enriched DNA motifs prepattern the accessible nucleosome organization in teleosts

BackgroundGene promoters in vertebrate genomes show distinct chromatin features such as stably positioned nucleosome array and DNA hypomethylation. The nucleosomes are known to have certain sequence preferences, and the prediction of nucleosome positioning from DNA sequence has been successful in some organisms such as yeast. However, at gene promoters where nucleosomes are much more stably positioned than in other regions, the sequence-based model has failed to work well, and sequence-independent mechanisms have been proposed.ResultsUsing DNase I-seq in medaka embryos, we demonstrated that hypomethylated domains (HMDs) specifically possess accessible nucleosome organization with longer linkers, and we reassessed the DNA sequence preference for nucleosome positioning in these specific regions. Remarkably, we found with a supervised machine learning algorithm, k-mer SVM, that nucleosome positioning in HMDs is accurately predictable from DNA sequence alone. Specific short sequences (6-mers) that contribute to the prediction are specifically enriched in HMDs and distribute periodically with approximately 200-bp intervals which prepattern the position of accessible linkers. Surprisingly, the sequence preference of the nucleosome and linker in HMDs is opposite from that reported previously. Furthermore, the periodicity of specific motifs at hypomethylated promoters was conserved in zebrafish.ConclusionThis study reveals strong link between nucleosome positioning and DNA sequence at vertebrate promoters, and we propose hypomethylated DNA-specific regulation of nucleosome positioning.

Dynamic unwrapping of nucleosomes by HsRAD51 that includes sliding and rotational motion of histone octamers.

Wrapping of genomic DNA into nucleosomes poses thermodynamic and kinetic barriers to biological processes such as replication, transcription, repair and recombination. Previous biochemical studies have demonstrated that in the presence of adenosine triphosphate (ATP) the human RAD51 (HsRAD51) recombinase can form a nucleoprotein filament (NPF) on double-stranded DNA (dsDNA) that is capable of unwrapping the nucleosomal DNA from the histone octamer (HO). Here, we have used single molecule Förster Resonance Energy Transfer (smFRET) to examine the real time nucleosome dynamics in the presence of the HsRAD51 NPF. We show that oligomerization of HsRAD51 leads to stepwise, but stochastic unwrapping of the DNA from the HO in the presence of ATP. The highly reversible dynamics observed in single-molecule trajectories suggests an antagonistic mechanism between HsRAD51 binding and rewrapping of the DNA around the HO. These stochastic dynamics were independent of the nucleosomal DNA sequence or the asymmetry created by the presence of a linker DNA. We also observed sliding and rotational oscillations of the HO with respect to the nucleosomal DNA. These studies underline the dynamic nature of even tightly associated protein–DNA complexes such as nucleosomes.

Conformational selection and dynamic adaptation upon linker histone binding to the nucleosome.

Linker histones are essential for DNA compaction in chromatin. They bind to nucleosomes in a 1:1 ratio forming chromatosomes. Alternative configurations have been proposed in which the globular domain of the linker histone H5 (gH5) is positioned either on- or off-dyad between the nucleosomal and linker DNAs. However, the dynamic pathways of chromatosome assembly remain elusive. Here, we studied the conformational plasticity of gH5 in unbound and off-dyad nucleosome-bound forms with classical and accelerated molecular dynamics simulations. We find that the unbound gH5 converts between open and closed conformations, preferring the closed form. However, the open gH5 contributes to a more rigid chromatosome and restricts the motion of the nearby linker DNA through hydrophobic interactions with thymidines. Moreover, the closed gH5 opens and reorients in accelerated simulations of the chromatosome. Brownian dynamics simulations of chromatosome assembly, accounting for a range of amplitudes of nucleosome opening and different nucleosome DNA sequences, support the existence of both on- and off-dyad binding modes of gH5 and reveal alternative, sequence and conformation-dependent chromatosome configurations. Taken together, these findings suggest that the conformational dynamics of linker histones and nucleosomes facilitate alternative chromatosome configurations through an interplay between induced fit and conformational selection.

Preparation of Recombinant Centromeric Nucleosomes and Formation of Complexes with Nonhistone Centromere Proteins.

Centromeres are present on each chromosome to direct proper segregation during cell division. The understanding of how the histone H3 variant, CENP-A, epigenetically marks the location of the centromere on the chromosome has been advanced, in part, through the study of histone complexes, nucleosomes, and nucleosomal complexes with nonhistone centromere proteins. In this chapter, we describe the preparation of recombinant versions of these complexes. The methodology is firmly rooted in classic nucleosome reconstitution methods, but we highlight the aspects of the preparations that diverge from those used for the methods established with canonical histones. We also provide a method for producing PCR-amplified nucleosomal DNA sequences in milligram quantities that is particularly useful for studies where multiple sequences and/or chemical modifications are desired. Lastly, we describe our approach to assemble and analyze a complex between the recombinant human CENP-A nucleosome and one of its binding partners, CENP-C.

Sequence structure of Lowary/Widom clones forming strong nucleosomes

Lowary and Widom selected from random sequences those which form exceptionally stable nucleosomes, including clone 601, the current champion of strong nucleosome (SN) sequences. This unique sequence database (LW sequences) carries sequence elements which confer stability on the nucleosomes formed on the sequences, and, thus, may serve as source of information on the structure of “ideal” or close to ideal nucleosome DNA sequence. An important clue is also provided by crystallographic study of Vasudevan and coauthors on clone 601 nucleosomes. It demonstrated that YR·YR dinucleotide stacks (primarily TA·TA) follow one another at distances 10 or 11 bases or multiples thereof, such that they all are located on the interface between DNA and histone octamer. Combining this important information with alignment of the YR-containing 10-mers and 11-mers from LW sequences, the bendability matrices of the stable nucleosome DNA are derived. The matrices suggest that the periodically repeated TA (YR), RR, and YY dinucleotides are the main sequence features of the SNs. This consensus coincides with the one for recently discovered SNs with visibly periodic DNA sequences. Thus, the experimentally observed stable LW nucleosomes and SNs derived computationally appear to represent the same entity – exceptionally stable SNs.

Review fifteen years of search for strong nucleosomes

Don Crothers, Mikael Kubista, Jon Widom, and their teams have been first to look for strong nucleosomes, in a bid to reveal the nucleosome positioning pattern(s) carried by the nucleosome DNA sequences. They were first to demonstrate that the nucleosome stability correlates with 10-11 base sequence periodicity, and that the strong nucleosomes localize preferentially in centromeres. This review describes these findings and their connection to recent discovery of the strong nucleosomes (SNs) with visibly periodic nucleosome DNA sequences.

Computational discovery of feature patterns in nucleosomal DNA sequences

The identification of important factors that affect nucleosome formation is critical to clarify nucleosome-forming mechanisms and the role of the nucleosome in gene regulation. Various features reported in the literature led to our hypothesis that multiple features can together contribute to nucleosome formation. Therefore, we compiled 779 features and developed a pattern discovery and scoring algorithm FFN (Finding Features for Nucleosomes) to identify feature patterns that are differentially enriched in nucleosome-forming sequences and nucleosome-depletion sequences. Applying FFN to genome-wide nucleosome occupancy data in yeast and human, we identified statistically significant feature patterns that may influence nucleosome formation, many of which are common to the two species. We found that both sequence and structural features are important in nucleosome occupancy prediction. We discovered that, even for the same feature combinations, variations in feature values may lead to differences in predictive power. We demonstrated that the identified feature patterns could be used to assist nucleosomal sequence prediction.

Nucleosomes shape DNA polymorphism and divergence.

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.

DNA Physical Parameters Modulate Nucleosome Positioning in the Saccharomyces cerevisiae Genome

Nucleosome positioning plays essential roles in various cellular processes. Although many efforts have been made in this area, the rules defining nucleosome positioning is still elusive. In the present study, DNA physical parameters derived from atomistic molecular dynamic simulations were introduced to analyze nucleosomal and linker DNA sequences. The distinct structural patterns between nucleosomal and linker sequences indicate that DNA physical parameters are suitable to describe nucleosomal DNA sequences and to reveal physical mechanisms of nucleosome positioning. Further analysis of DNA flexibility around regulatory regions indicates that nucleosome positioning is closely correlated with sequence flexibility. These results demonstrate that DNA physical parameters are useful for the in silico nucleosome positioning prediction.

Universal full-length nucleosome mapping sequence probe

For the computational sequence-directed mapping of the nucleosomes, the knowledge of the nucleosome positioning motifs – 10–11 base long sequences – and respective matrices of bendability, is not sufficient, since there is no justified way to fuse these motifs in one continuous nucleosome DNA sequence. Discovery of the strong nucleosome (SN) DNA sequences, with visible sequence periodicity allows derivation of the full-length nucleosome DNA bendability pattern as matrix or consensus sequence. The SN sequences of three species (A. thaliana, C. elegans, and H. sapiens) are aligned (512 sequences for each species), and long (115 dinucleotides) matrices of bendability derived for the species. The matrices have strong common property – alternation of runs of purine–purine (RR) and pyrimidine–pyrimidine (YY) dinucleotides, with average period 10.4 bases. On this basis the universal [R,Y] consensus of the nucleosome DNA sequence is derived, with exactly defined positions of respective penta- and hexamers RRRRR, RRRRRR, YYYYY, and YYYYYY.

Strong nucleosomes of A. thaliana concentrate in centromere regions

Earlier identified strongest nucleosome DNA sequences of A. thaliana, those with visible 10–11 base sequence periodicity, are mapped along chromosomes. Resulting positional distributions reveal distinct maxima, one per chromosome, located in the centromere regions. Sequence-directed nucleosome mapping demonstrates that the strong nucleosomes (SNs) make tight arrays, several ‘parallel’ nucleosomes each, suggesting a columnar chromatin structure. The SNs represent a new class of centromeric nucleosomes, presumably, participating in synapsis of chromatids and securing the centromere architecture.

Visible periodicity of strong nucleosome DNA sequences

Fifteen years ago, Lowary and Widom assembled nucleosomes on synthetic random sequence DNA molecules, selected the strongest nucleosomes and discovered that the TA dinucleotides in these strong nucleosome sequences often appear at 10–11 bases from one another or at distances which are multiples of this period. We repeated this experiment computationally, on large ensembles of natural genomic sequences, by selecting the strongest nucleosomes – i.e. those with such distances between like-named dinucleotides, multiples of 10.4 bases, the structural and sequence period of nucleosome DNA. The analysis confirmed the periodicity of TA dinucleotides in the strong nucleosomes, and revealed as well other periodic sequence elements, notably classical AA and TT dinucleotides. The matrices of DNA bendability and their simple linear forms – nucleosome positioning motifs – are calculated from the strong nucleosome DNA sequences. The motifs are in full accord with nucleosome positioning sequences derived earlier, thus confirming that the new technique, indeed, detects strong nucleosomes. Species- and isochore-specific variations of the matrices and of the positioning motifs are demonstrated. The strong nucleosome DNA sequences manifest the highest hitherto nucleosome positioning sequence signals, showing the dinucleotide periodicities in directly observable rather than in hidden form.

Nucleosomal TATA-switch: Competing orientations of TATA on the nucleosome

Transcription is known to be affected by the rotational setting of the transcription response elements within nucleosomes. We studied the rotational positioning of the TATA box, the most universal promoter motif. We applied a bioinformatic high-resolution nucleosome mapping technique to eukaryotic promoters. Our results show that the nucleosome DNA sequence harboring the TATA box encodes alternative rotational positions for the same piece of DNA. This may serve for switching the gene activity on and off.

The Effect of Nucleosome Stability on Gene Expression Level

Nucleosomes, which are the basic packaging units of chromatin, are stably positioned in promoters upstream of most stress-induced genes. These promoter nucleosomes are generally thought to repress gene expression due to exclusion; they prevent transcription factors from accessing their binding sites on the DNA. Some promoter nucleosomes, however, do not directly occlude transcription factor binding sites, and therefore their role in gene expression remains unknown. To understand the non-exclusive role of nucleosomes in gene expression, we designed a model promoter in budding yeast where a nucleosome intervenes between a transcription factor binding site in nucleosome depleted region and the transcription start site. To vary nucleosome stability, we constructed promoter variants with several different GC contents in the nucleosomal DNA sequence. We then measured the downstream gene expression level from these promoter variants at different induction levels using a fluorescent protein reporter. Our preliminary data show that relatively high nucleosome stability does not always correlate with relatively low gene expression level. This result suggests that nucleosomes might contribute to gene expression in a cooperative manner either by bringing otherwise distant DNA sites close together or by anchoring the SWI/SNF chromatin remodeling complex for efficient removal of adjacent nucleosomes.