Theoretical analysis of epigenetic cell memory by nucleosome modification

Chromosomal regions can adopt stable and heritable alternative states resulting in bistable gene expression without changes to the DNA sequence. Such epigenetic control is often associated with alternative covalent modifications of histones. The stability and heritability of the states are thought to involve positive feedback where modified nucleosomes recruit enzymes that similarly modify nearby nucleosomes. We developed a simplified stochastic model for dynamic nucleosome modification based on the silent mating-type region of the yeast Schizosaccharomyces pombe. We show that the mechanism can give strong bistability that is resistant both to high noise due to random gain or loss of nucleosome modifications and to random partitioning upon DNA replication. However, robust bistability required: (1) cooperativity, the activity of more than one modified nucleosome, in the modification reactions and (2) that nucleosomes occasionally stimulate modification beyond their neighbor nucleosomes, arguing against a simple continuous spreading of nucleosome modification.     

A simple histone code opens many paths to epigenetics

Nucleosomes can be covalently modified by addition of various chemical groups on several of their exposed histone amino acids. These modifications are added and removed by enzymes (writers) and can be recognized by nucleosome-binding proteins (readers). Linking a reader domain and a writer domain that recognize and create the same modification state should allow nucleosomes in a particular modification state to recruit enzymes that create that modification state on nearby nucleosomes. This positive feedback has the potential to provide the alternative stable and heritable states required for epigenetic memory. However, analysis of simple histone codes involving interconversions between only two or three types of modified nucleosomes has revealed only a few circuit designs that allow heritable bistability. Here we show by computer simulations that a histone code involving alternative modifications at two histone positions, producing four modification states, combined with reader-writer proteins able to distinguish these states, allows for hundreds of different circuits capable of heritable bistability. These expanded possibilities result from multiple ways of generating two-step cooperativity in the positive feedback - through alternative pathways and an additional, novel cooperativity motif. Our analysis reveals other properties of such epigenetic circuits. They are most robust when the dominant nucleosome types are different at both modification positions and are not the type inserted after DNA replication. The dominant nucleosome types often recruit enzymes that create their own type or destroy the opposing type, but never catalyze their own destruction. The circuits appear to be evolutionary accessible; most circuits can be changed stepwise into almost any other circuit without losing heritable bistability. Thus, our analysis indicates that systems that utilize an expanded histone code have huge potential for generating stable and heritable nucleosome modification states and identifies the critical features of such systems.     

Nucleosome positioning and nucleosome stacking: two faces of the same coin

Nucleosomes are regularly spaced along eukaryotic genomes. In the emerging model, known as "statistical positioning", this spacing is due to steric repulsion between nucleosomes and to the presence of nucleosome excluding barriers on the genome. However, new experimental evidence recently challenged the "statistical positioning" model (Z. Zhang et al., Science, 2011, 332(6032), 977-980). We propose here that the regular spacing can be better explained by adding attractive interactions between nucleosomes. In our model those attractions are due to the fact that nucleosomes are stacked in regular chromatin fibers. In a self-reinforcing mechanism, regular nucleosome spacing promotes in turn nucleosome stacking. We first show that this model can precisely account for the nucleosome spacing observed in Saccharomyces cerevisiae. We then use a simple toy model to show that attraction between nucleosomes can fasten the formation of the chromatin fiber.     

Nucleosome acetylation sequencing to study the establishment of chromatin acetylation

The establishment of posttranslational chromatin modifications is a major mechanism for regulating how genomic DNA is utilized. However, current in vitro chromatin assays do not monitor histone modifications at individual nucleosomes. Here we describe a strategy, nucleosome acetylation sequencing, that allows us to read the amount of modification at each nucleosome. In this approach, a bead-bound trinucleosome substrate is enzymatically acetylated with radiolabeled acetyl CoA by the SAGA complex from Saccharomyces cerevisae. The product is digested by restriction enzymes that cut at unique sites between the nucleosomes and then counted to quantify the extent of acetylation at each nucleosomal site. We find that we can sensitively, specifically, and reproducibly follow enzyme-mediated nucleosome acetylation. Applying this strategy, when acetylation proceeds extensively, its distribution across nucleosomes is relatively uniform. However, when substrates are used that contain nucleosomes mutated at the major sites of SAGA-mediated acetylation, or that are studied under initial rate conditions, changes in the acetylation distribution can be observed. Nucleosome acetylation sequencing should be applicable to analyzing a wide range of modifications. Additionally, because our trinucleosomes synthesis strategy is highly modular and efficient, it can be used to generate nucleosomal systems in which nucleosome composition differs across the array.     

The Saccharomyces cerevisiae Piccolo NuA4 histone acetyltransferase complex requires the Enhancer of Polycomb A domain and chromodomain to acetylate nucleosomes

Chromatin modification complexes are key gene regulatory factors which posttranslationally modify the histone component of chromatin with epigenetic marks. To address what features of chromatin modification complexes are responsible for the specific recognition of nucleosomes compared to naked histones, we have performed a functional dissection of the Esa1-containing Saccharomyces cerevisiae Piccolo NuA4 histone acetyltransferase complex. Our studies define the Piccolo determinants sufficient to assemble its three subunits into a complex as well as Piccolo determinants sufficient to specifically acetylate a chromatin template. We find that the conserved Enhancer of Polycomb A (EPcA) homology region of the Epl1 component and the N-terminal 165 amino acids of the Yng2 component of Piccolo are sufficient with Esa1 to specifically act on nucleosomes. We also find that the Esa1 chromodomain plays a critical role in Piccolo's ability to distinguish between histones and nucleosomes. In particular, specific point mutations in the chromodomain putative hydrophobic cage which strongly hinder growth in yeast greatly reduce histone acetyltransferase activity on nucleosome substrates, independent of histone methylation or other modifications. However, the chromodomain is not required for Piccolo to bind to nucleosomes, suggesting a role for the chromodomain in a catalysis step after nucleosome binding.     

Characterization of Dnmt1 Binding and DNA Methylation on Nucleosomes and Nucleosomal Arrays

The packaging of DNA into nucleosomes and the organisation into higher order structures of chromatin limits the access of sequence specific DNA binding factors to DNA. In cells, DNA methylation is preferentially occuring in the linker region of nucleosomes, suggesting a structural impact of chromatin on DNA methylation. These observations raise the question whether DNA methyltransferases are capable to recognize the nucleosomal substrates and to modify the packaged DNA. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the maintenance DNA methyltransferase Dnmt1. Our binding studies show that Dnmt1 has a DNA length sensing activity, binding cooperatively to DNA, and requiring a minimal DNA length of 20 bp. Dnmt1 needs linker DNA to bind to nucleosomes and most efficiently recognizes nucleosomes with symmetric DNA linkers. Footprinting experiments reveal that Dnmt1 binds to both DNA linkers exiting the nucleosome core. The binding pattern correlates with the efficient methylation of DNA linkers. However, the enzyme lacks the ability to methylate nucleosomal CpG sites on mononucleosomes and nucleosomal arrays, unless chromatin remodeling enzymes create a dynamic chromatin state. In addition, our results show that Dnmt1 functionally interacts with specific chromatin remodeling enzymes to enable complete methylation of hemi-methylated DNA in chromatin.     

RCC1 Uses a Conformationally Diverse Loop Region to Interact with the Nucleosome: A Model for the RCC1-Nucleosome Complex

The binding of RCC1 (regulator of chromosome condensation 1) to chromatin is critical for cellular processes such as mitosis, nucleocytoplasmic transport, and nuclear envelope formation because RCC1 recruits the small GTPase Ran (Ras-related nuclear protein) to chromatin and sets up a Ran-GTP gradient around the chromosomes. However, the molecular mechanism by which RCC1 binds to nucleosomes, the repeating unit of chromatin, is not known. We have used biochemical approaches to test structural models for how the RCC1 β-propeller protein could bind to the nucleosome. In contrast to the prevailing model, RCC1 does not appear to use the β-propeller face opposite to its Ran-binding face to interact with nucleosomes. Instead, we find that RCC1 uses a conformationally flexible loop region we have termed the switchback loop in addition to its N-terminal tail to bind to the nucleosome. The juxtaposition of the RCC1 switchback loop to its Ran binding surface suggests a novel mechanism for how nucleosome-bound RCC1 recruits Ran to chromatin. Furthermore, this model accounts for previously unexplained observations for how Ran can interact with the nucleosome both dependent and independent of RCC1 and how binding of the nucleosome can enhance RCC1's Ran nucleotide exchange activity.     

Predicting nucleosome positions on the DNA: combining intrinsic sequence preferences and remodeler activities

Nucleosome positions on the DNA are determined by the intrinsic affinities of histone proteins to a given DNA sequence and by the ATP-dependent activities of chromatin remodeling complexes that can translocate nucleosomes with respect to the DNA. Here, we report a theoretical approach that takes into account both contributions. In the theoretical analysis two types of experiments have been considered: in vitro experiments with a single reconstituted nucleosome and in vivo genome-scale mapping of nucleosome positions. The effect of chromatin remodelers was described by iteratively redistributing the nucleosomes according to certain rules until a new steady state was reached. Three major classes of remodeler activities were identified: (i) the establishment of a regular nucleosome spacing in the vicinity of a strong positioning signal acting as a boundary, (ii) the enrichment/depletion of nucleosomes through amplification of intrinsic DNA-sequence-encoded signals and (iii) the removal of nucleosomes from high-affinity binding sites. From an analysis of data for nucleosome positions in resting and activated human CD4⁺ T cells [Schones et al., Cell 132, p. 887] it was concluded that the redistribution of a nucleosome map to a new state is greatly facilitated if the remodeler complex translocates the nucleosome with a preferred directionality.     

Modification of DNA in chromatin with methyltransferase from Haemophilus influenzae Rd.

The accessibility of DNA in nucleosome dimers (as a model of the chromosomal chain of nucleosomes) was determined by means of modification methylases from Haemophilus influenzae Rd. Using these enzymes, the rate of modification of nucleosome dimers is about one fifth the rate observed with protein-free DNA from chromatin subunit dimers. Methylated DNA sites in nucleosome dimers are readily accessible to micrococcal nuclease. The analysis of the fragment pattern of nucleosomes after methylation and mild nuclease treatment reveals that the methylated sites are predominantly located in the internucleosomal linker DNA. Polylysine binding experiments further support this interpretation. This compound preferentially interacts with the nucleosomal core DNA and protects it against internal cleavage. It neither affects the degradation of methylated sites drastically nor does it inhibit the methylation of nucleosome dimers. Thus, a combination of protection, cleavage and modification is proposed as a useful tool for the analysis of the structure of chromatin.     

Base excision repair in chromatin: Insights from reconstituted systems

The process of base excision repair has been completely reconstituted in vitro and structural and biochemical properties of the component enzymes thoroughly studied on naked DNA templates. More recent work in this field aims to understand how BER operates on the natural substrate, chromatin [1], [2]. Toward this end, a number of researchers, including the Smerdon group, have focused attention to understand how individual enzymes and reconstituted BER operate on nucleosome substrates. While nucleosomes were once thought to completely restrict access of DNA-dependent factors, the surprising finding from these studies suggests that at least some BER components can utilize target DNA bound within nucleosomes as substrates for their enzymatic processes. This data correlates well with both structural studies of these enzymes and our developing understanding of nucleosome conformation and dynamics. While more needs to be learned, these studies highlight the utility of reconstituted BER and chromatin systems to inform our understanding of in vivo biological processes.     

Nucleosome gaping supports a functional structure for the 30nm chromatin fiber

The biological functions played by the nucleus of eukaryotic cells and especially those involved in cellular differentiation not only depend on the genomic sequence but also on all the proteins which form the nucleo-protein complex named chromatin. The tridimensional organization of this huge polymer involves many structural levels, the most basic one being the nucleosome. Nucleosomes further organize into the so-called 30nm fiber, which, according to recent works, is likely to be the main functional level of chromatin. We wish here to propose a plausible structure for the 30nm chromatin fiber that could explain its functional role. In our model, silenced chromatin is locked by nucleosome stacking interactions. This is achieved by a conformational transition within the nucleosome core particle (NCP) which allows nucleosomes to stack along two helices without bending the DNA linkers. We used molecular modeling to check that this conformational transition was plausible. Then we proposed to modify the well-known two-angle model according to these atomic level results. The emerging picture is an allosteric behavior of the nucleosomes induced by their collective organization within the 30nm chromatin fiber.     

Chromatin dynamics controls epigenetic domain formation

In multicellular organisms, nucleosomes carry epigenetic information that defines distinct patterns of gene expression, which are inherited over multiple generations. The enhanced capacity for information storage arises by nucleosome modifications, which are triggered by enzymes. Modified nucleosomes can transfer the mark to others that are in proximity by a positive-feedback (modification begets modification) mechanism. We created a generic polymer model, referred to as 3DSpreader, in which each bead, representing a nucleosome, stochastically switches between unmodified (U) and modified (M) states depending on the states of the neighbors. Modification begins at a specific nucleation site (NS) that is permanently in the M state, and could spread to other loci that is dictated by chromatin dynamics. Transfer of marks among the non-nucleation loci occurs stochastically as chromatin evolves in time. If the spreading rate is slower than the chromatin relaxation rate, which is biologically pertinent, then finite-sized domains form, driven by contacts between nucleosomes through a three-dimensional looping mechanism. Surprisingly, simulations based on the 3DSpreader model result in finite bounded domains that arise without the need for any boundary elements. Maintenance of spatially and temporally stable domains requires the presence of the NS, whose removal eliminates finite-sized modified domains. The theoretical predictions are in excellent agreement with experimental data for H3K9me3 spreading in mouse embryonic stem cells.     

DNA structure, nucleosome placement and chromatin remodelling: a perspective

A major question in chromatin biology is to what extent the sequence of DNA directly determines the genetic and chromatin organization of a eukaryotic genome? We consider two aspects to this question: the DNA sequence-specified positioning of nucleosomes and the determination of NDRs (nucleosome-depleted regions) or barriers. We argue that, in budding yeast, while DNA sequence-specified nucleosome positioning may contribute to positions flanking the regions lacking nucleosomes, DNA thermodynamic stability is a major component determinant of the genetic organization of this organism.     

Mutations to the histone H3 αN region selectively alter the outcome of ATP-dependent nucleosome-remodelling reactions

Mutational analysis of the histone H3 N-terminal region has shown it to play an important role both in chromatin function in vivo and nucleosome dynamics in vitro. Here we use a library of mutations in the H3 N-terminal region to investigate the contribution of this region to the action of the ATP-dependent remodelling enzymes Chd1, RSC and SWI/SNF. All of the enzymes were affected differently by the mutations with Chd1 being affected the least and RSC being most sensitive. In addition to affecting the rate of remodelling by RSC, some mutations prevented RSC from moving nucleosomes to locations in which DNA was unravelled. These observations illustrate that the mechanisms by which different ATP-dependent remodelling enzymes act are sensitive to different features of nucleosome structure. They also show how alterations to histones can affect the products generated as a result of ATP-dependent remodelling reactions.