Distinctive sequence patterns in metazoan and yeast nucleosomes: implications for linker histone binding to AT-rich and methylated DNA

Interpretable deep residual network uncovers nucleosome positioning and associated features

Nucleosomes represent elementary building units of eukaryotic chromosomes and consist of DNA wrapped around a histone octamer flanked by linker DNA segments. Nucleosomes are central in epigenetic pathways and their genomic positioning is associated with regulation of gene expression, DNA replication, DNA methylation and DNA repair, among other functions. Building on prior discoveries that DNA sequences noticeably affect nucleosome positioning, our objective is to identify nucleosome positions and related features across entire genome. Here, we introduce an interpretable framework based on the concepts of deep residual networks (NuPoSe). Trained on high-coverage human experimental MNase-seq data, NuPoSe is able to learn sequence and structural patterns associated with nucleosome organization in human genome. NuPoSe can be also applied to unseen data from different organisms and cell types. Our findings point to 43 informative features, most of them constitute tri-nucleotides, di-nucleotides and one tetra-nucleotide. Most features are significantly associated with the nucleosomal structural characteristics, namely, periodicity of nucleosomal DNA and its location with respect to a histone octamer. Importantly, we show that features derived from the 27 bp linker DNA flanking nucleosomes contribute up to 10% to the quality of the prediction model. This, along with the comprehensive training sets, deep-learning architecture, and feature selection method, may contribute to the NuPoSe's 80-89% classification accuracy on different independent datasets.

Histone N-tails modulate sequence-specific positioning of nucleosomes

The precise mechanisms governing sequence-dependent positioning of nucleosomes on DNA remain unknown in detail. Existing algorithms, taking into account the sequence-dependent deformability of DNA and its interactions with the histone globular domains, predict rotational setting of only 65% of human nucleosomes mapped in vivo. To uncover novel factors responsible for the nucleosome positioning, we analyzed potential involvement of the histone N-tails in this process. To this aim, we reconstituted the H2A/H4 N-tailless nucleosomes on human BRCA1 DNA (~100 kb) and compared their positions and sequences with those of the wild-type nucleosomes. In the case of H2A tailless nucleosomes, the AT content of DNA sequences is changed locally at superhelical location (SHL) ±4, while maintaining the same rotational setting as their wild-type counterparts. Conversely, the H4 tailless nucleosomes display widespread changes of the AT content near SHL ±1 and SHL ±2, where the H4 N-tails interact with DNA. Furthermore, a substantial number of H4 tailless nucleosomes exhibit rotational setting opposite to that of the wild-type nucleosomes. Thus, our findings strongly suggest that the histone N-tails are operative in selection of nucleosome positions, which may have wide-ranging implications for epigenetic modulation of chromatin states.

DNA Sequence-Dependent Binding of Linker Histone gH1 Regulates Nucleosome Conformations

Sequence-dependent binding between DNA and proteins in chromatin is an essential part of gene expression. Linker histone H1 is an important protein in the regulation of chromatin compartmentalization and compaction, and its binding with the nucleosome is sensitive to the DNA sequence. Although the interactions of H1 and DNA have been widely investigated, the mechanism of nucleosome conformation changes induced by the DNA-sequence-dependent binding with gH1 (globular H1.0) remains largely unclear at the atomic level. In the present molecular dynamics simulations, both linker and dyad DNAs were mutated to investigate the conformational changes of the nucleosome induced by the sequence-dependent binding of gH1 based on the on-dyad binding mode. Our results indicate that gH1 is insensitive to the DNA sequence of the dyad DNA but presents an apparent preference to linker DNA with an AT-rich sequence. Moreover, this specific binding induces the entry/exit region of a nucleosome to a tight conformation and regulates the accessibility of core histones. Considering that the entry/exit region of the nucleosome is a crucial binding site for many functional proteins related to gene expression, the conformational change at this region could represent an important gene regulation signal.

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.

A critical role for linker DNA in higher-order folding of chromatin fibers

Nucleosome-nucleosome interactions drive the folding of nucleosomal arrays into dense chromatin fibers. A better physical account of the folding of chromatin fibers is necessary to understand the role of chromatin in regulating DNA transactions. Here, we studied the unfolding pathway of regular chromatin fibers as a function of single base pair increments in linker length, using both rigid base-pair Monte Carlo simulations and single-molecule force spectroscopy. Both computational and experimental results reveal a periodic variation of the folding energies due to the limited flexibility of the linker DNA. We show that twist is more restrictive for nucleosome stacking than bend, and find the most stable stacking interactions for linker lengths of multiples of 10 bp. We analyzed nucleosomes stacking in both 1- and 2-start topologies and show that stacking preferences are determined by the length of the linker DNA. Moreover, we present evidence that the sequence of the linker DNA also modulates nucleosome stacking and that the effect of the deletion of the H4 tail depends on the linker length. Importantly, these results imply that nucleosome positioning in vivo not only affects the phasing of nucleosomes relative to DNA but also directs the higher-order structure of chromatin.

Minimal Cylinder Analysis Reveals the Mechanical Properties of Oncogenic Nucleosomes

Histone variants regulate replication, transcription, DNA damage repair, and chromosome segregation. Though widely accepted as a paradigm, it has not been rigorously demonstrated that histone variants encode unique mechanical properties. Here, we present a new theoretical approach called minimal cylinder analysis that uses strain fluctuations to determine the Young's modulus of nucleosomes from all-atom molecular dynamics simulations. Recently, we validated this computational tool against in vitro single-molecule nanoindentation of histone variant nucleosomes. In this report, we further extend minimal cylinder analysis to study the biophysical properties of hybrid nucleosomes that are known to exist in human cancer cells and contain H3 histone variants CENP-A and H3.3. Here, we report that the heterotypic nucleosome has an intermediate elasticity (8.5 ± 0.5 MPa) compared to CENP-A (6.2 ± 0.4 MPa) and H3 (9.8 ± 0.7 MPa) and that the dynamics of both canonical and CENP-A nucleosomes are preserved and partitioned across the nucleosome pseudodyad. Furthermore, we investigate the mechanism by which the elasticity of these heterotypic nucleosomes augments cryptic binding surfaces. From these analyses, we predict that the heterotypic nucleosome is permissive to the binding of one copy of the kinetochore protein CENP-C while still retaining a closed DNA end configuration required for linker histone H1 to bind. We discuss that the ectopic deposition of CENP-A in cancer by H3.3 chaperones HIRA and DAXX may fortuitously result in hybrid nucleosome formation. Using these results, we propose biological outcomes that might arise when such heterotypic nucleosomes occupy large regions of the genome.

Chromatin Compaction Multiscale Modeling: A Complex Synergy Between Theory, Simulation, and Experiment

Understanding the mechanisms that trigger chromatin compaction, its patterns, and the factors they depend on, is a fundamental and still open question in Biology. Chromatin compacts and reinforces DNA and is a stable but dynamic structure, to make DNA accessible to proteins. In recent years, computational advances have provided larger amounts of data and have made large-scale simulations more viable. Experimental techniques for the extraction and reconstitution of chromatin fibers have improved, reinvigorating theoretical and experimental interest in the topic and stimulating debate on points previously considered as certainties regarding chromatin. A great assortment of approaches has emerged, from all-atom single-nucleosome or oligonucleosome simulations to various degrees of coarse graining, to polymer models, to fractal-like structures and purely topological models. Different fiber-start patterns have been studied in theory and experiment, as well as different linker DNA lengths. DNA is a highly charged macromolecule, making ionic and electrostatic interactions extremely important for chromatin topology and dynamics. Indeed, the repercussions of varying ionic concentration have been extensively examined at the computational level, using all-atom, coarse-grained, and continuum techniques. The presence of high-curvature AT-rich segments in DNA can cause conformational variations, attesting to the fact that the role of DNA is both structural and electrostatic. There have been some tentative attempts to describe the force fields governing chromatin conformational changes and the energy landscapes of these transitions, but the intricacy of the system has hampered reaching a consensus. The study of chromatin conformations is an intrinsically multiscale topic, influenced by a wide range of biological and physical interactions, spanning from the atomic to the chromosome level. Therefore, powerful modeling techniques and carefully planned experiments are required for an overview of the most relevant phenomena and interactions. The topic provides fertile ground for interdisciplinary studies featuring a synergy between theoretical and experimental scientists from different fields and the cross-validation of respective results, with a multi-scale perspective. Here, we summarize some of the most representative approaches, and focus on the importance of electrostatics and solvation, often overlooked aspects of chromatin modeling.

Chromatosome Structure and Dynamics from Molecular Simulations

Chromatosomes are fundamental units of chromatin structure that are formed when a linker histone protein binds to a nucleosome. The positioning of the linker histone on the nucleosome influences the packing of chromatin. Recent simulations and experiments have shown that chromatosomes adopt an ensemble of structures that differ in the geometry of the linker histone-nucleosome interaction. In this article we review the application of Brownian, Monte Carlo, and molecular dynamics simulations to predict the structure of linker histone-nucleosome complexes, to study the binding mechanisms involved, and to predict how this binding affects chromatin fiber structure. These simulations have revealed the sensitivityof the chromatosome structure to variations in DNA and linker histone sequence, as well as to posttranslational modifications, thereby explaining the structural variability observed in experiments. We propose that a concerted application of experimental and computational approaches will reveal the determinants of chromatosome structural variability and how it impacts chromatin packing.

Sensitive effect of linker histone binding mode and subtype on chromatin condensation

The complex role of linker histone (LH) on chromatin compaction regulation has been highlighted by recent discoveries of the effect of LH binding variability and isoforms on genome structure and function. Here we examine the effect of two LH variants and variable binding modes on the structure of chromatin fibers. Our mesoscale modeling considers oligonucleosomes with H1C and H1E, bound in three different on and off-dyad modes, and spanning different LH densities (0.5–1.6 per nucleosome), over a wide range of physiologically relevant nucleosome repeat lengths (NRLs). Our studies reveal an LH-variant and binding-mode dependent heterogeneous ensemble of fiber structures with variable packing ratios, sedimentation coefficients, and persistence lengths. For maximal compaction, besides dominantly interacting with parental DNA, LHs must have strong interactions with nonparental DNA and promote tail/nonparental core interactions. An off-dyad binding of H1E enables both; others compromise compaction for bendability. We also find that an increase of LH density beyond 1 is best accommodated in chromatosomes with one on-dyad and one off-dyad LH. We suggest that variable LH binding modes and concentrations are advantageous, allowing tunable levels of chromatin condensation and DNA accessibility/interactions. Thus, LHs add another level of epigenetic regulation of chromatin.

Modelling and DNA topology of compact 2-start and 1-start chromatin fibres

We have investigated the structure of the most compact 30-nm chromatin fibres by modelling those with 2-start or 1-start crossed-linker organisations. Using an iterative procedure we obtained possible structural solutions for fibres of the highest possible compaction permitted by physical constraints, including the helical repeat of linker DNA. We find that this procedure predicts a quantized nucleosome repeat length (NRL) and that only fibres with longer NRLs (≥197 bp) can more likely adopt the 1-start organisation. The transition from 2-start to 1-start fibres is consistent with reported differing binding modes of the linker histone. We also calculate that in 1-start fibres the DNA constrains more torsion (as writhe) than 2-start fibres with the same NRL and that the maximum constraint obtained is in accord with previous experimental results. We posit that the coiling of the fibre is driven by overtwisting of linker DNA which, in the most compact forms - for example, in echinoderm sperm and avian erythrocytes - could adopt a helical repeat of ∼10 bp/turn. We argue that in vivo the total twist of linker DNA could be modulated by interaction with other abundant chromatin-associated proteins and by epigenetic modifications of the C-terminal tail of linker histones.

Identification and characterization of hirudin-HN, a new thrombin inhibitor, from the salivary glands of Hirudo nipponia

Transcriptome sequencing data (6.5 Gb) of the salivary glands of the haematophagous leech Hirudo nipponia was obtained by using the BGIseq-500 platform. After identification and analysis, one transcript (Unigene5370) was annotated to hirudin HV3 from Hirudo medicinalis with an e-value of 1e-29 and was named hirudin-HN. This transcript was a new thrombin inhibitor gene belonging to the proteinase inhibitor I14 (hirudin) family. Hirudin-HN, with a 270-bp cDNA, encodes an 89-aa protein containing a 20-aa signal peptide. The mature hirudin-HN protein contains the typical structural characteristics of hirudin, e.g., three conserved disulfide bonds and the PKP and DFxxIP motifs. Proteins (Hir and M-Hir) were obtained via prokaryotic expression, and the mature hirudin-HN protein was shown to have anticoagulant activity and thrombin affinity by using the chromogenic substrate S2238 and surface plasmon resonance (SPR) interaction analysis, respectively. The N-terminal structure of the mature hirudin-HN protein was shown to be important for anticoagulant activity by comparing the activity and thrombin affinity of Hir and M-Hir. The abundances of Hirudin-HN mRNA and protein were higher in the salivary glands of starving animals than in those of feeding or fed leeches. These results provided a foundation for further study on the structure-function relationship of hirudin-HN with thrombin.

The nucleosome position-encoding WW/SS sequence pattern is depleted in mammalian genes relative to other eukaryotes

Nucleosomal DNA sequences generally follow a well-known pattern with ∼10-bp periodic WW (where W is A or T) dinucleotides that oscillate in phase with each other and out of phase with SS (where S is G or C) dinucleotides. However, nucleosomes with other DNA patterns have not been systematically analyzed. Here, we focus on an opposite pattern, namely anti-WW/SS pattern, in which WW dinucleotides preferentially occur at DNA sites that bend into major grooves and SS (where S is G or C) dinucleotides are often found at sites that bend into minor grooves. Nucleosomes with the anti-WW/SS pattern are widespread and exhibit a species- and context-specific distribution in eukaryotic genomes. Unlike non-mammals (yeast, nematode and fly), there is a positive correlation between the enrichment of anti-WW/SS nucleosomes and RNA Pol II transcriptional levels in mammals (mouse and human). Interestingly, such enrichment is not due to underlying DNA sequence. In addition, chromatin remodeling complexes have an impact on the abundance but not on the distribution of anti-WW/SS nucleosomes in yeast. Our data reveal distinct roles of cis- and trans-acting factors in the rotational positioning of nucleosomes between non-mammals and mammals. Implications of the anti-WW/SS sequence pattern for RNA Pol II transcription are discussed.

Chromatin structures condensed by linker histones

In eukaryotic cells, genomic DNA exists in the form of chromatin through association with histone proteins, which consist of four core histone (H2A, H2B, H3, and H4) families and one linker histone (H1) family. The core histones bind to DNA to form the nucleosome, the recurring structural unit of chromatin. The linker histone binds to the nucleosome to form the next structural unit of chromatin, the chromatosome, which occurs dominantly in metazoans. Linker histones also play an essential role in condensing chromatin to form higher order structures. Unlike the core histones in the formation of the nucleosome, the role of linker histone in the formation of the chromatosome and high-order chromatin structure is not well understood. Nevertheless, exciting progress in the structural studies of chromatosomes and nucleosome arrays condensed by linker histones has been made in the last several years. In this mini-review, we discuss these recent experimental results and provide some perspectives for future studies.

Toward an Ensemble View of Chromatosome Structure: A Paradigm Shift from One to Many

There is renewed interest in linker histone (LH)-nucleosome binding and how LHs influence eukaryotic DNA compaction. For a long time, the goal was to uncover "the structure of the chromatosome," but recent studies of LH-nucleosome complexes have revealed an ensemble of structures. Notably, the reconstituted LH-nucleosome complexes used in experiments rarely correspond to the sequence combinations present in organisms. For a full understanding of the determinants of the distribution of the chromatosome structural ensemble, studies must include a complete description of the sequences and experimental conditions used, and be designed to enable systematic evaluation of sequence and environmental effects.

Regulation of chromatin folding by conformational variations of nucleosome linker DNA

Linker DNA conformational variability has been proposed to direct nucleosome array folding into more or less compact chromatin fibers but direct experimental evidence for such models are lacking. Here, we tested this hypothesis by designing nucleosome arrays with A-tracts at specific locations in the nucleosome linkers to induce inward (AT-IN) and outward (AT-OUT) bending of the linker DNA. Using electron microscopy and analytical centrifugation techniques, we observed spontaneous folding of AT-IN nucleosome arrays into highly compact structures, comparable to those induced by linker histone H1. In contrast, AT-OUT nucleosome arrays formed less compact structures with decreased nucleosome interactions similar to wild-type nucleosome arrays. Adding linker histone H1 further increased compaction of the A-tract arrays while maintaining structural differences between them. Furthermore, restriction nuclease digestion revealed a strongly reduced accessibility of nucleosome linkers in the compact AT-IN arrays. Electron microscopy analysis and 3D computational Monte Carlo simulations are consistent with a profound zigzag linker DNA configuration and closer nucleosome proximity in the AT-IN arrays due to inward linker DNA bending. We propose that the evolutionary preferred positioning of A-tracts in DNA linkers may control chromatin higher-order folding and thus influence cellular processes such as gene expression, transcription and DNA repair.

Emerging roles of linker histones in regulating chromatin structure and function

Together with core histones, which make up the nucleosome, the linker histone (H1) is one of the five main histone protein families present in chromatin in eukaryotic cells. H1 binds to the nucleosome to form the next structural unit of metazoan chromatin, the chromatosome, which may help chromatin to fold into higher-order structures. Despite their important roles in regulating the structure and function of chromatin, linker histones have not been studied as extensively as core histones. Nevertheless, substantial progress has been made recently. The first near-atomic resolution crystal structure of a chromatosome core particle and an 11 Å resolution cryo-electron microscopy-derived structure of the 30 nm nucleosome array have been determined, revealing unprecedented details about how linker histones interact with the nucleosome and organize higher-order chromatin structures. Moreover, several new functions of linker histones have been discovered, including their roles in epigenetic regulation and the regulation of DNA replication, DNA repair and genome stability. Studies of the molecular mechanisms of H1 action in these processes suggest a new paradigm for linker histone function beyond its architectural roles in chromatin.

MeCP2 and CTCF: enhancing the cross-talk of silencers

This paper provides a brief introductory review of the most recent advances in our knowledge about the structural and functional aspects of two transcriptional regulators: MeCP2, a protein whose mutated forms are involved in Rett syndrome; and CTCF, a constitutive transcriptional insulator. This is followed by a description of the PTMs affecting these two proteins and an analysis of their known interacting partners. A special emphasis is placed on the recent studies connecting these two proteins, focusing on the still poorly understood potential structural and functional interactions between the two of them on the chromatin substrate. An overview is provided for some of the currently known genes that are dually regulated by these two proteins. Finally, a model is put forward to account for their possible involvement in their regulation of gene expression.

Structure and Functions of Linker Histones

Linker histones such as variants H1, H5, and other similar proteins play an important role in regulation of chromatin structure and dynamics. However, interactions of linker histones with DNA and proteins, as well as specific functions of their different variants, are poorly studied. This is because they acquire tertiary structure only when interacting with a nucleosome, and because of limitations of currently available methods. However, deeper investigation of linker histones and their interactions with other proteins will address a number of important questions - from structure of compacted chromatin to regulation of early embryogenesis. In this review, structures of histone H1 variants and its interaction with chromatin DNA are considered. A possible functional significance of different H1 variants, a role of these proteins in maintaining interphase chromatin structure, and interactions of linker histones with other cellular proteins are also discussed.

Feature-based classification of human transcription factors into hypothetical sub-classes related to regulatory function

Background

Transcription factors are key proteins in the regulation of gene transcription. An important step in this process is the opening of chromatin in order to make genomic regions available for transcription. Data on DNase I hypersensitivity has previously been used to label a subset of transcription factors as Pioneers, Settlers and Migrants to describe their potential role in this process. These labels represent an interesting hypothesis on gene regulation and possibly a useful approach for data analysis, and therefore we wanted to expand the set of labeled transcription factors to include as many known factors as possible. We have used a well-annotated dataset of 1175 transcription factors as input to supervised machine learning methods, using the subset with previously assigned labels as training set. We then used the final classifier to label the additional transcription factors according to their potential role as Pioneers, Settlers and Migrants. The full set of labeled transcription factors was used to investigate associated properties and functions of each class, including an analysis of interaction data for transcription factors based on DNA co-binding and protein-protein interactions. We also used the assigned labels to analyze a previously published set of gene lists associated with a time course experiment on cell differentiation.

Results

The analysis showed that the classification of transcription factors with respect to their potential role in chromatin opening largely was determined by how they bind to DNA. Each subclass of transcription factors was enriched for properties that seemed to characterize the subclass relative to its role in gene regulation, with very general functions for Pioneers, whereas Migrants to a larger extent were associated with specific processes. Further analysis showed that the expanded classification is a useful resource for analyzing other datasets on transcription factors with respect to their potential role in gene regulation. The analysis of transcription factor interaction data showed complementary differences between the subclasses, where transcription factors labeled as Pioneers often interact with other transcription factors through DNA co-binding, whereas Migrants to a larger extent use protein-protein interactions. The analysis of time course data on cell differentiation indicated a shift in the regulatory program associated with Pioneer-like transcription factors during differentiation.

Conclusions

The expanded classification is an interesting resource for analyzing data on gene regulation, as illustrated here on transcription factor interaction data and data from a time course experiment. The potential regulatory function of transcription factors seems largely to be determined by how they bind DNA, but is also influenced by how they interact with each other through cooperativity and protein-protein interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1349-2) contains supplementary material, which is available to authorized users.

Cell cycle progression in Caulobacter requires a nucleoid-associated protein with high AT sequence recognition

Faithful cell cycle progression in the dimorphic bacterium Caulobacter crescentus requires spatiotemporal regulation of gene expression and cell pole differentiation. We discovered an essential DNA-associated protein, GapR, that is required for Caulobacter growth and asymmetric division. GapR interacts with adenine and thymine (AT)-rich chromosomal loci, associates with the promoter regions of cell cycle-regulated genes, and shares hundreds of recognition sites in common with known master regulators of cell cycle-dependent gene expression. GapR target loci are especially enriched in binding sites for the transcription factors GcrA and CtrA and overlap with nearly all of the binding sites for MucR1, a regulator that controls the establishment of swarmer cell fate. Despite constitutive synthesis, GapR accumulates preferentially in the swarmer compartment of the predivisional cell. Homologs of GapR, which are ubiquitous among the α-proteobacteria and are encoded on multiple bacteriophage genomes, also accumulate in the predivisional cell swarmer compartment when expressed in Caulobacter The Escherichia coli nucleoid-associated protein H-NS, like GapR, selectively associates with AT-rich DNA, yet it does not localize preferentially to the swarmer compartment when expressed exogenously in Caulobacter, suggesting that recognition of AT-rich DNA is not sufficient for the asymmetric accumulation of GapR. Further, GapR does not silence the expression of H-NS target genes when expressed in E. coli, suggesting that GapR and H-NS have distinct functions. We propose that Caulobacter has co-opted a nucleoid-associated protein with high AT recognition to serve as a mediator of cell cycle progression.

A Small Number of Residues Can Determine if Linker Histones Are Bound On or Off Dyad in the Chromatosome

Linker histones bind to the nucleosome and regulate the structure and function of chromatin. We have previously shown that the globular domains of chicken H5 and Drosophila H1 linker histones bind to the nucleosome with on- or off-dyad modes, respectively. To explore the determinant for the distinct binding modes, we investigated the binding of a mutant globular domain of H5 to the nucleosome. This mutant, termed GH5_pMut, includes substitutions of five globular domain residues of H5 with the corresponding residues in the globular domain of Drosophila H1. The residues at these five positions play important roles in nucleosome binding by either H5 or Drosophila H1. NMR and spin-labeling experiments showed that GH5_pMut bound to the nucleosome off the dyad. We further found that the nucleosome array condensed by either the GH5_pMut or the globular domain of Drosophila H1 displayed a similar sedimentation coefficient, whereas the same nucleosome array condensed by the wild-type globular domain of H5 showed a much larger sedimentation coefficient. Moreover, NMR and spin-labeling results from the study of the nucleosome in complex with the full-length human linker histone H1.0, whose globular domain shares high sequence conservation with the corresponding globular domain of H5, are consistent with an on-dyad binding mode. Taken together, our results suggest that a small number of residues in the globular domain of a linker histone can control its binding location on the nucleosome and higher-order chromatin structure.

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.

HistoneDB 2.0: a histone database with variants--an integrated resource to explore histones and their variants

Compaction of DNA into chromatin is a characteristic feature of eukaryotic organisms. The core (H2A, H2B, H3, H4) and linker (H1) histone proteins are responsible for this compaction through the formation of nucleosomes and higher order chromatin aggregates. Moreover, histones are intricately involved in chromatin functioning and provide a means for genome dynamic regulation through specific histone variants and histone post-translational modifications. ‘HistoneDB 2.0 – with variants’ is a comprehensive database of histone protein sequences, classified by histone types and variants. All entries in the database are supplemented by rich sequence and structural annotations with many interactive tools to explore and compare sequences of different variants from various organisms. The core of the database is a manually curated set of histone sequences grouped into 30 different variant subsets with variant-specific annotations. The curated set is supplemented by an automatically extracted set of histone sequences from the non-redundant protein database using algorithms trained on the curated set. The interactive web site supports various searching strategies in both datasets: browsing of phylogenetic trees; on-demand generation of multiple sequence alignments with feature annotations; classification of histone-like sequences and browsing of the taxonomic diversity for every histone variant. HistoneDB 2.0 is a resource for the interactive comparative analysis of histone protein sequences and their implications for chromatin function.

Database URL:http://www.ncbi.nlm.nih.gov/projects/HistoneDB2.0

Histone H1 Variants in Arabidopsis Are Subject to Numerous Post-Translational Modifications, Both Conserved and Previously Unknown in Histones, Suggesting Complex Functions of H1 in Plants

Linker histones (H1s) are conserved and ubiquitous structural components of eukaryotic chromatin. Multiple non-allelic variants of H1, which differ in their DNA/nucleosome binding properties, co-exist in animal and plant cells and have been implicated in the control of genetic programs during development and differentiation. Studies in mammals and Drosophila have revealed diverse post-translational modifications of H1s, most of which are of unknown function. So far, it is not known how this pattern compares with that of H1s from other major lineages of multicellular Eukaryotes. Here, we show that the two main H1variants of a model flowering plant Arabidopsis thaliana are subject to a rich and diverse array of post-translational modifications. The distribution of these modifications in the H1 molecule, especially in its globular domain (GH1), resembles that occurring in mammalian H1s, suggesting that their functional significance is likely to be conserved. While the majority of modifications detected in Arabidopsis H1s, including phosphorylation, acetylation, mono- and dimethylation, formylation, crotonylation and propionylation, have also been reported in H1s of other species, some others have not been previously identified in histones.

Crop epigenetics and the molecular hardware of genotype × environment interactions

Crop plants encounter thermal environments which fluctuate on a diurnal and seasonal basis. Future climate resilient cultivars will need to respond to thermal profiles reflecting more variable conditions, and harness plasticity that involves regulation of epigenetic processes and complex genomic regulatory networks. Compartmentalization within plant cells insulates the genomic central processing unit within the interphase nucleus. This review addresses the properties of the chromatin hardware in which the genome is embedded, focusing on the biophysical and thermodynamic properties of DNA, histones and nucleosomes. It explores the consequences of thermal and ionic variation on the biophysical behavior of epigenetic marks such as DNA cytosine methylation (5mC), and histone variants such as H2A.Z, and how these contribute to maintenance of chromatin integrity in the nucleus, while enabling specific subsets of genes to be regulated. Information is drawn from theoretical molecular in vitro studies as well as model and crop plants and incorporates recent insights into the role epigenetic processes play in mediating between environmental signals and genomic regulation. A preliminary speculative framework is outlined, based on the evidence of what appears to be a cohesive set of interactions at molecular, biophysical and electrostatic level between the various components contributing to chromatin conformation and dynamics. It proposes that within plant nuclei, general and localized ionic homeostasis plays an important role in maintaining chromatin conformation, whilst maintaining complex genomic regulation that involves specific patterns of epigenetic marks. More generally, reversible changes in DNA methylation appear to be consistent with the ability of nuclear chromatin to manage variation in external ionic and temperature environment. Whilst tentative, this framework provides scope to develop experimental approaches to understand in greater detail the internal environment of plant nuclei. It is hoped that this will generate a deeper understanding of the molecular mechanisms underlying genotype × environment interactions that may be beneficial for long-term improvement of crop performance in less predictable climates.

Linker histones: History and current perspectives

Although the overall structure of the fifth histone (linker histone, H1) is understood, its location on the nucleosome is only partially defined. Whilst it is clear that H1 helps condense the chromatin fibre, precisely how this is achieved remains to be determined. H1 is not a general gene repressor in that although it must be displaced from transcription start sites for activity to occur, there is only partial loss along the body of genes. How the deposition and removal of H1 occurs in particular need of further study. Linker histones are highly abundant nuclear proteins about which we know too little.

Novel nucleosomal particles containing core histones and linker DNA but no histone H1

Eukaryotic chromosomal DNA is assembled into regularly spaced nucleosomes, which play a central role in gene regulation by determining accessibility of control regions. The nucleosome contains ∼147 bp of DNA wrapped ∼1.7 times around a central core histone octamer. The linker histone, H1, binds both to the nucleosome, sealing the DNA coils, and to the linker DNA between nucleosomes, directing chromatin folding. Micrococcal nuclease (MNase) digests the linker to yield the chromatosome, containing H1 and ∼160 bp, and then converts it to a core particle, containing ∼147 bp and no H1. Sequencing of nucleosomal DNA obtained after MNase digestion (MNase-seq) generates genome-wide nucleosome maps that are important for understanding gene regulation. We present an improved MNase-seq method involving simultaneous digestion with exonuclease III, which removes linker DNA. Remarkably, we discovered two novel intermediate particles containing 154 or 161 bp, corresponding to 7 bp protruding from one or both sides of the nucleosome core. These particles are detected in yeast lacking H1 and in H1-depleted mouse chromatin. They can be reconstituted in vitro using purified core histones and DNA. We propose that these ‘proto-chromatosomes’ are fundamental chromatin subunits, which include the H1 binding site and influence nucleosome spacing independently of H1.

Structural Mechanisms of Nucleosome Recognition by Linker Histones

Linker histones bind to the nucleosome and regulate the structure of chromatin and gene expression. Despite more than three decades of effort, the structural basis of nucleosome recognition by linker histones remains elusive. Here, we report the crystal structure of the globular domain of chicken linker histone H5 in complex with the nucleosome at 3.5 Å resolution, which is validated using nuclear magnetic resonance spectroscopy. The globular domain sits on the dyad of the nucleosome and interacts with both DNA linkers. Our structure integrates results from mutation analyses and previous cross-linking and fluorescence recovery after photobleach experiments, and it helps resolve the long debate on structural mechanisms of nucleosome recognition by linker histones. The on-dyad binding mode of the H5 globular domain is different from the recently reported off-dyad binding mode of Drosophila linker histone H1. We demonstrate that linker histones with different binding modes could fold chromatin to form distinct higher-order structures.

Possible role of H1 histone in replication timing

AT-rich repetitive DNA sequences become late replicating during cell differentiation. Replication timing is not correlated with LINE density in human cells (Ryba et al. 2010). However, short and properly spaced runs of oligo dA or dT present in nuclear matrix attachment regions (MARs) of the genome are good candidates for elements of AT-rich repetitive late replicating DNA. MAR attachment to the nuclear matrix is negatively regulated by chromatin binding of H1 histone, but this is counteracted by H1 phosphorylation, high mobility group proteins or, indirectly, core histone acetylation. Fewer MAR attachments correlates positively with longer average DNA loop size, longer replicons and an increase of late replicating DNA.

Prediction of nucleosome rotational positioning in yeast and human genomes based on sequence-dependent DNA anisotropy

BACKGROUND

An organism's DNA sequence is one of the key factors guiding the positioning of nucleosomes within a cell's nucleus. Sequence-dependent bending anisotropy dictates how DNA is wrapped around a histone octamer. One of the best established sequence patterns consistent with this anisotropy is the periodic occurrence of AT-containing dinucleotides (WW) and GC-containing dinucleotides (SS) in the nucleosomal locations where DNA is bent in the minor and major grooves, respectively. Although this simple pattern has been observed in nucleosomes across eukaryotic genomes, its use for prediction of nucleosome positioning was not systematically tested.

RESULTS

We present a simple computational model, termed the W/S scheme, implementing this pattern, without using any training data. This model accurately predicts the rotational positioning of nucleosomes both in vitro and in vivo, in yeast and human genomes. About 65 - 75% of the experimentally observed nucleosome positions are predicted with the precision of one to two base pairs. The program is freely available at http://people.rit.edu/fxcsbi/WS_scheme/. We also introduce a simple and efficient way to compare the performance of different models predicting the rotational positioning of nucleosomes.

CONCLUSIONS

This paper presents the W/S scheme to achieve accurate prediction of rotational positioning of nucleosomes, solely based on the sequence-dependent anisotropic bending of nucleosomal DNA. This method successfully captures DNA features critical for the rotational positioning of nucleosomes, and can be further improved by incorporating additional terms related to the translational positioning of nucleosomes in a species-specific manner.

Terminate and make a loop: regulation of transcriptional directionality

•

Transcriptional directionality is controlled by premature transcription termination.

•

Transcriptional directionality is enforced by gene looping.

•

mRNA-specific termination signals and factors are required for gene looping.

Bidirectional promoters are a common feature of many eukaryotic organisms from yeast to humans. RNA Polymerase II that is recruited to this type of promoter can start transcribing in either direction using alternative DNA strands as the template. Such promiscuous transcription can lead to the synthesis of unwanted transcripts that may have negative effects on gene expression. Recent studies have identified transcription termination and gene looping as critical players in the enforcement of promoter directionality. Interestingly, both mechanisms share key components. Here, we focus on recent findings relating to the transcriptional output of bidirectional promoters.

Chemical map of Schizosaccharomyces pombe reveals species-specific features in nucleosome positioning

Using a recently developed chemical approach, we have generated a genome-wide map of nucleosomes in vivo in Schizosaccharomyces pombe (S. pombe) at base pair resolution. The shorter linker length previously identified in S. pombe is due to a preponderance of nucleosomes separated by ∼4/5 bp, placing nucleosomes on opposite faces of the DNA. The periodic dinucleotide feature thought to position nucleosomes is equally strong in exons as in introns, demonstrating that nucleosome positioning information can be superimposed on coding information. Unlike the case in Saccharomyces cerevisiae, A/T-rich sequences are enriched in S. pombe nucleosomes, particularly at ±20 bp around the dyad. This difference in nucleosome binding preference gives rise to a major distinction downstream of the transcription start site, where nucleosome phasing is highly predictable by A/T frequency in S. pombe but not in S. cerevisiae, suggesting that the genomes and DNA binding preferences of nucleosomes have coevolved in different species. The poly (dA-dT) tracts affect but do not deplete nucleosomes in S. pombe, and they prefer special rotational positions within the nucleosome, with longer tracts enriched in the 10- to 30-bp region from the dyad. S. pombe does not have a well-defined nucleosome-depleted region immediately upstream of most transcription start sites; instead, the -1 nucleosome is positioned with the expected spacing relative to the +1 nucleosome, and its occupancy is negatively correlated with gene expression. Although there is generally very good agreement between nucleosome maps generated by chemical cleavage and micrococcal nuclease digestion, the chemical map shows consistently higher nucleosome occupancy on DNA with high A/T content.

Structural insights into the histone H1-nucleosome complex

Linker H1 histones facilitate formation of higher-order chromatin structures and play important roles in various cell functions. Despite several decades of effort, the structural basis of how H1 interacts with the nucleosome remains elusive. Here, we investigated Drosophila H1 in complex with the nucleosome, using solution nuclear magnetic resonance spectroscopy and other biophysical methods. We found that the globular domain of H1 bridges the nucleosome core and one 10-base pair linker DNA asymmetrically, with its α3 helix facing the nucleosomal DNA near the dyad axis. Two short regions in the C-terminal tail of H1 and the C-terminal tail of one of the two H2A histones are also involved in the formation of the H1-nucleosome complex. Our results lead to a residue-specific structural model for the globular domain of the Drosophila H1 in complex with the nucleosome, which is different from all previous experiment-based models and has implications for chromatin dynamics in vivo.

Combined micrococcal nuclease and exonuclease III digestion reveals precise positions of the nucleosome core/linker junctions: implications for high-resolution nucleosome mapping

Micrococcal nuclease (MNase) is extensively used in genome-wide mapping of nucleosomes but its preference for AT-rich DNA leads to errors in establishing precise positions of nucleosomes. Here, we show that the MNase digestion of nucleosomes assembled on a strong nucleosome positioning sequence, Widom's clone 601, releases nucleosome cores whose sizes are strongly affected by the linker DNA sequence. Our experiments produced nucleosomal DNA sizes varying between 147 and 155 bp, with positions of the MNase cuts reflecting positions of the A⋅T pairs rather than the nucleosome core/linker junctions determined by X-ray crystallography. Extent of chromatosomal DNA protection by linker histone H1 also depends on the linker DNA sequence. Remarkably, we found that a combined treatment with MNase and exonuclease III (exoIII) overcomes MNase sequence preference producing nucleosomal DNA trimmed symmetrically and precisely at the core/linker junctions regardless of the underlying DNA sequence. We propose that combined MNase/exoIII digestion can be applied to in situ chromatin for unbiased genome-wide mapping of nucleosome positions that is not influenced by DNA sequences at the core/linker junctions. The same approach can be also used for the precise mapping of the extent of linker DNA protection by H1 and other protein factors associated with nucleosome linkers.

Mechanisms by which transcription factors gain access to target sequence elements in chromatin

Transcription factors (TF) bind DNA sequence motifs, but the presence of a consensus DNA element is not sufficient to direct TF binding to chromatin. Recent genomic data have revealed that accessibility, as measured by DNase sensitivity and the presence of active histone marks, is necessary for TF binding. DNA sequence provides the initial specification of the accessibility of DNA elements within chromatin that permits TF binding. In yeast, it is known that poly(dA-dT) tracts directly encode low-nucleosome occupancy at promoters. Recent evidence suggests that CpG islands in mammals are inherently refractory to higher-order chromatin structure and remain accessible, despite favoring nucleosome formation in vitro. Taken together, these studies support a model for how accessibility originates and then propagates throughout regulatory cascades and development.

Transcriptional activation of yeast genes disrupts intragenic nucleosome phasing

Nucleosomes often undergo extensive rearrangement when genes are activated for transcription. We have shown previously, using paired-end sequencing of yeast nucleosomes, that major changes in chromatin structure occur when genes are activated by 3-aminotriazole (3AT), an inducer of the transcriptional activator Gcn4. Here, we provide a global analysis of these data. At the genomic level, nucleosomes are regularly phased relative to the transcription start site. However, for a subset of 234 strongly induced genes, this phasing is much more irregular after induction, consistent with the loss of some nucleosomes and the re-positioning of the remaining nucleosomes. To address the nature of this rearrangement, we developed the inter-nucleosome distance auto-correlation (DAC) function. At long range, DAC analysis indicates that nucleosomes have an average spacing of 162 bp, consistent with the reported repeat length. At short range, DAC reveals a 10.25-bp periodicity, implying that nucleosomes in overlapping positions are rotationally related. DAC analysis of the 3AT-induced genes suggests that transcription activation coincides with rearrangement of nucleosomes into irregular arrays with longer spacing. Sequence analysis of the +1 nucleosomes belonging to the 45 most strongly activated genes reveals a distinctive periodic oscillation in the A/T-dinucleotide occurrence that is present throughout the nucleosome and extends into the linker. This unusual pattern suggests that the +1 nucleosomes might be prone to sliding, thereby facilitating transcription.

From crystal and NMR structures, footprints and cryo-electron-micrographs to large and soft structures: nanoscale modeling of the nucleosomal stem

The interaction of histone H1 with linker DNA results in the formation of the nucleosomal stem structure, with considerable influence on chromatin organization. In a recent paper [Syed,S.H., Goutte-Gattat,D., Becker,N., Meyer,S., Shukla,M.S., Hayes,J.J., Everaers,R., Angelov,D., Bednar,J. and Dimitrov,S. (2010) Single-base resolution mapping of H1-nucleosome interactions and 3D organization of the nucleosome. Proc. Natl Acad. Sci. USA, 107, 9620-9625], we published results of biochemical footprinting and cryo-electron-micrographs of reconstituted mono-, di- and tri-nucleosomes, for H1 variants with different lengths of the cationic C-terminus. Here, we present a detailed account of the analysis of the experimental data and we include thermal fluctuations into our nano-scale model of the stem structure. By combining (i) crystal and NMR structures of the nucleosome core particle and H1, (ii) the known nano-scale structure and elasticity of DNA, (iii) footprinting information on the location of protected sites on the DNA backbone and (iv) cryo-electron micrographs of reconstituted tri-nucleosomes, we arrive at a description of a polymorphic, hierarchically organized stem with a typical length of 20 ± 2 base pairs. A comparison to linker conformations inferred for poly-601 fibers with different linker lengths suggests, that intra-stem interactions stabilize and facilitate the formation of dense chromatin fibers.

Structural and dynamic properties of linker histone H1 binding to DNA

Found in all eukaryotic cells, linker histones H1 are known to bind to and rearrange nucleosomal linker DNA. In vitro, the fundamental nature of H1∕DNA interactions has attracted wide interest among research communities-from biologists to physicists. Hence, H1∕DNA binding processes and structural and dynamical information about these self-assemblies are of broad importance. Targeting a quantitative understanding of H1 induced DNA compaction mechanisms, our strategy is based on using small-angle x-ray microdiffraction in combination with microfluidics. The usage of microfluidic hydrodynamic focusing devices facilitates a microscale control of these self-assembly processes, which cannot be achieved using conventional bulk setups. In addition, the method enables time-resolved access to structure formation in situ, in particular, to transient intermediate states. The observed time dependent structure evolution shows that the H1∕DNA interaction can be described as a two-step process: an initial unspecific binding of H1 to DNA is followed by a rearrangement of molecules within the formed assemblies. The second step is most likely induced by interactions between the DNA and the H1's charged side chains. This leads to an increase in lattice spacing within the DNA∕protein assembly and induces a decrease in the correlation length of the mesophases, probably due to a local bending of the DNA.

Distinct interfacial biclique patterns between ssDNA-binding proteins and those with dsDNAs

We introduce a new motif called interfacial biclique pattern to study the difference between double-stranded DNA-binding proteins (DSBs, most of them also known to play the role as transcriptional factors) and single-stranded DNA-binding proteins (SSBs) which are found to involve in many applications recently. An interfacial biclique pattern in a protein-DNA complex usually consists of a group of residues and a group of nucleotides such that every residue has a contact to all of the bases. The proposal of this idea is based on a biological redundancy mechanism that: a site mutation has little influence on the other residues to recognize the target nucleotides and vice versa. The distribution of the residues on the interfacial motifs is investigated to identify distinct stable preferred residues, stable un-preferred residues and unstable preferred residues between SSBs and DSBs. We also examine residue co-occurrence and residue-base association rules in the interfacial motifs to uncover the different choices of residue combinations by SSBs and DSBs that have contacts with one or more bases. We found that DSBs and SSBs have their own right residues at the right places for the binding preference and association with nucleotides. Some of our results can be supported by literature work.

On the structure and dynamics of the complex of the nucleosome and the linker histone

Several different models of the linker histone (LH)–nucleosome complex have been proposed, but none of them has unambiguously revealed the position and binding sites of the LH on the nucleosome. Using Brownian dynamics-based docking together with normal mode analysis of the nucleosome to account for the flexibility of two flanking 10 bp long linker DNAs (L-DNA), we identified binding modes of the H5-LH globular domain (GH5) to the nucleosome. For a wide range of nucleosomal conformations with the L-DNA ends less than 65 Å apart, one dominant binding mode was identified for GH5 and found to be consistent with fluorescence recovery after photobleaching (FRAP) experiments. GH5 binds asymmetrically with respect to the nucleosomal dyad axis, fitting between the nucleosomal DNA and one of the L-DNAs. For greater distances between L-DNA ends, docking of GH5 to the L-DNA that is more restrained and less open becomes favored. These results suggest a selection mechanism by which GH5 preferentially binds one of the L-DNAs and thereby affects DNA dynamics and accessibility and contributes to formation of a particular chromatin fiber structure. The two binding modes identified would, respectively, favor a tight zigzag chromatin structure or a loose solenoid chromatin fiber.

A structural perspective on the where, how, why, and what of nucleosome positioning

The DNA in eukaryotic chromatin is packed by histones into arrays of repeating units called nucleosomes. Each nucleosome contains a nucleosome core, where the DNA is wrapped around a histone octamer, and a stretch of relatively unconstrained DNA called the linker DNA. Since nucleosome cores occlude the DNA from many DNA-binding factors, their positions provide important clues for understanding chromatin packing and gene regulation. Here we review the recent advances in the genome-wide mapping of nucleosome positions, the molecular and structural determinants of nucleosome positioning, and the importance of nucleosome positioning in chromatin higher order folding and transcriptional regulation.

Structure-based analysis of DNA sequence patterns guiding nucleosome positioning in vitro

Recent studies of genome-wide nucleosomal organization suggest that the DNA sequence is one of the major determinants of nucleosome positioning. Although the search for underlying patterns encoded in nucleosomal DNA has been going on for about 30 years, our knowledge of these patterns still remains limited. Based on our evaluations of DNA deformation energy, we developed new scoring functions to predict nucleosome positioning. There are three principal differences between our approach and earlier studies: (i) we assume that the length of nucleosomal DNA varies from 146 to 147 bp; (ii) we consider the anisotropic flexibility of pyrimidine-purine (YR) dimeric steps in the context of their neighbors (e.g., YYRR versus RYRY); (iii) we postulate that alternating AT-rich and GC-rich motifs reflect sequence-dependent interactions between histone arginines and DNA in the minor groove. Using these functions, we analyzed 20 nucleosome positions mapped in vitro at single nucleotide resolution (including clones 601, 603, 605, the pGUB plasmid, chicken beta-globin and three 5S rDNA genes). We predicted 15 of the 20 positions with 1-bp precision, and two positions with 2-bp precision. The predicted position of the '601' nucleosome (i.e., the optimum of the computed score) deviates from the experimentally determined unique position by no more than 1 bp - an accuracy exceeding that of earlier predictions. Our analysis reveals a clear heterogeneity of the nucleosomal sequences which can be divided into two groups based on the positioning 'rules' they follow. The sequences of one group are enriched by highly deformable YR/YYRR motifs at the minor-groove bending sites SHL+/- 3.5 and +/- 5.5, which is similar to the alpha-satellite sequence used in most crystallized nucleosomes. Apparently, the positioning of these nucleosomes is determined by the interactions between histones H2A/H2B and the terminal parts of nucleosomal DNA. In the other group (that includes the '601' clone) the same YR/YYRR motifs occur predominantly at the sites SHL +/- 1.5. The interaction between the H3/H4 tetramer and the central part of the nucleosomal DNA is likely to be responsible for the positioning of nucleosomes of this group, and the DNA trajectory in these nucleosomes may differ in detail from the published structures. Thus, from the stereochemical perspective, the in vitro nucleosomes studied here follow either an X-ray-like pattern (with strong deformations in the terminal parts of nucleosomal DNA), or an alternative pattern (with the deformations occurring predominantly in the central part of the nucleosomal DNA). The results presented here may be useful for genome-wide classification of nucleosomes, linking together structural and thermodynamic characteristics of nucleosomes with the underlying DNA sequence patterns guiding their positions.

H2A.Z and H3.3 histone variants affect nucleosome structure: biochemical and biophysical studies

Histone variants play important roles in regulation of chromatin structure and function. To understand the structural role played by histone variants H2A.Z and H3.3, both of which are implicated in transcription regulation, we conducted extensive biochemical and biophysical analysis on mononucleosomes reconstituted from either random-sequence DNA derived from native nucleosomes or a defined DNA nucleosome positioning sequence and recombinant human histones. Using established electrophoretic and sedimentation analysis methods, we compared the properties of nucleosomes containing canonical histones and histone variants H2A.Z and H3.3 (in isolation or in combination). We find only subtle differences in the compaction and stability of the particles. Interestingly, both H2A.Z and H3.3 affect nucleosome positioning, either creating new positions or altering the relative occupancy of the existing nucleosome position space. On the other hand, only H2A.Z-containing nucleosomes exhibit altered linker histone binding. These properties could be physiologically significant as nucleosome positions and linker histone binding partly determine factor binding accessibility.