Dynamics of nucleosome invasion by DNA binding proteins

Why Are Nucleosome Breathing Dynamics Asymmetric?

In eukaryotic cells, DNA is bound to nucleosomes, but DNA segments occasionally unbind in the process known as nucleosome breathing. Although DNA can unwrap simultaneously from both ends of the nucleosome (symmetric breathing), experiments indicate that DNA prefers to dissociate from only one end (asymmetric breathing). However, the molecular origin of the asymmetry is not understood. We developed a new theoretical approach that gives microscopic explanations of asymmetric breathing. It is based on a stochastic description that leads to a comprehensive evaluation of dynamics by using effective free-energy landscapes. It is shown that asymmetric breathing follows the kinetically preferred pathways. In addition, it is also found that asymmetric breathing leads to a faster target search by transcription factors. Theoretical predictions, supported by computer simulations, agree with experiments. It is proposed that nature utilizes the symmetry of nucleosome breathing to achieve a better dynamic accessibility of chromatin for more efficient genetic regulation.

DNA-Histone Cross-Link Formation via Hole Trapping in Nucleosome Core Particles

Radical cations (holes) produced in DNA by ionizing radiation and other oxidants yield DNA-protein cross-links (DPCs). Detailed studies of DPC formation in chromatin via this process are lacking. We describe here a comprehensive examination of DPC formation within nucleosome core particles (NCPs), which are the monomeric component of chromatin. DNA holes are introduced at defined sites within NCPs that are constructed from the bottom-up. DPCs form at DNA holes in yields comparable to those of alkali-labile DNA lesions that result from water trapping. DPC-forming efficiency and site preference within the NCP are dependent on translational and rotational positioning. Mass spectrometry and the use of mutant histones reveal that lysine residues in histone N-terminal tails and amino termini are responsible for the DPC formation. These studies are corroborated by computational simulation at the microsecond time scale, showing a wide range of interactions that can precede DPC formation. Three consecutive dGs, which are pervasive in the human genome, including G-quadruplex-forming sequences, are sufficient to produce DPCs that could impact gene expression.

Role of Nucleosome Sliding in the Protein Target Search for Covered DNA Sites

Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.

How Do Nucleosome Dynamics Regulate Protein Search on DNA?

Nearly three-fourths of all eukaryotic DNA is occupied by nucleosomes, protein-DNA complexes comprising octameric histone core proteins and ∼150 base pairs of DNA. In addition to acting as a DNA compaction vehicle, the dynamics of nucleosomes regulate the DNA site accessibility for the nonhistone proteins, thereby controlling regulatory processes involved in determining the cell identity and cell fate. Here, we propose an analytical framework to analyze the role of nucleosome dynamics on the target search process of transcription factors through a simple discrete-state stochastic description of the search process. By considering the experimentally determined kinetic rates associated with protein and nucleosome dynamics as the only inputs, we estimate the target search time of a protein via first-passage probability calculations separately during nucleosome breathing and sliding dynamics. Although both the nucleosome dynamics permit transient access to the DNA sites that are otherwise occluded by the histone proteins, our result suggests substantial differences between the protein search mechanism on a nucleosome performing breathing and sliding dynamics. Furthermore, we identify the molecular factors that influence the search efficiency and demonstrate how these factors together portray a highly dynamic landscape of gene regulation. Our analytical results are validated using extensive Monte Carlo simulations.

Nucleosome Breathing Facilitates the Search for Hidden DNA Sites by Pioneer Transcription Factors

Transfer of genetic information starts with transcription factors (TFs) binding to specific sites on DNA. But in living cells, DNA is mostly covered by nucleosomes. There are proteins, known as pioneer TFs, that can efficiently reach the DNA sites hidden by nucleosomes, although the underlying mechanisms are not understood. Using the recently proposed idea of interaction-compensation mechanism, we develop a stochastic model for the target search on DNA with nucleosome breathing. It is found that nucleosome breathing can significantly accelerate the search by pioneer TFs in comparison to situations without breathing. We argue that this is the result of the interaction-compensation mechanism that allows proteins to enter the inner nucleosome region through the outer DNA segment. It is suggested that nature optimized pioneer TFs to take advantage of nucleosome breathing. The presented theoretical picture provides a possible microscopic explanation for the successful invasion of nucleosome-buried genes.

Basic helix-loop-helix pioneer factors interact with the histone octamer to invade nucleosomes and generate nucleosome-depleted regions

Nucleosomes drastically limit transcription factor (TF) occupancy, while pioneer transcription factors (PFs) somehow circumvent this nucleosome barrier. In this study, we compare nucleosome binding of two conserved S. cerevisiae basic helix-loop-helix (bHLH) TFs, Cbf1 and Pho4. A cryo-EM structure of Cbf1 in complex with the nucleosome reveals that the Cbf1 HLH region can electrostatically interact with exposed histone residues within a partially unwrapped nucleosome. Single-molecule fluorescence studies show that the Cbf1 HLH region facilitates efficient nucleosome invasion by slowing its dissociation rate relative to DNA through interactions with histones, whereas the Pho4 HLH region does not. In vivo studies show that this enhanced binding provided by the Cbf1 HLH region enables nucleosome invasion and ensuing repositioning. These structural, single-molecule, and in vivo studies reveal the mechanistic basis of dissociation rate compensation by PFs and how this translates to facilitating chromatin opening inside cells.

An associative memory Hamiltonian model for DNA and nucleosomes

A model for DNA and nucleosomes is introduced with the goal of studying chromosomes from a single base level all the way to higher-order chromatin structures. This model, dubbed the Widely Editable Chromatin Model (WEChroM), reproduces the complex mechanics of the double helix including its bending persistence length and twisting persistence length, and their respective temperature dependence. The WEChroM Hamiltonian is composed of chain connectivity, steric interactions, and associative memory terms representing all remaining interactions leading to the structure, dynamics, and mechanical characteristics of the B-DNA. Several applications of this model are discussed to demonstrate its applicability. WEChroM is used to investigate the behavior of circular DNA in the presence of positive and negative supercoiling. We show that it recapitulates the formation of plectonemes and of structural defects that relax mechanical stress. The model spontaneously manifests an asymmetric behavior with respect to positive or negative supercoiling, similar to what was previously observed in experiments. Additionally, we show that the associative memory Hamiltonian is also capable of reproducing the free energy of partial DNA unwrapping from nucleosomes. WEChroM is designed to emulate the continuously variable mechanical properties of the 10nm fiber and, by virtue of its simplicity, is ready to be scaled up to molecular systems large enough to investigate the structural ensembles of genes. WEChroM is implemented in the OpenMM simulation toolkits and is freely available for public use.

The nucleosome unwrapping free energy landscape defines distinct regions of transcription factor accessibility and kinetics

Transcription factors (TF) require access to target sites within nucleosomes to initiate transcription. The target site position within the nucleosome significantly influences TF occupancy, but how is not quantitatively understood. Using ensemble and single-molecule fluorescence measurements, we investigated the targeting and occupancy of the transcription factor, Gal4, at different positions within the nucleosome. We observe a dramatic decrease in TF occupancy to sites extending past 30 base pairs (bp) into the nucleosome which cannot be explained by changes in the TF dissociation rate or binding site orientation. Instead, the nucleosome unwrapping free energy landscape is the primary determinant of Gal4 occupancy by reducing the Gal4 binding rate. The unwrapping free energy landscape defines two distinct regions of accessibility and kinetics with a boundary at 30 bp into the nucleosome where the inner region is over 100-fold less accessible. The Gal4 binding rate in the inner region no longer depends on its concentration because it is limited by the nucleosome unwrapping rate, while the frequency of nucleosome rewrapping decreases because Gal4 exchanges multiple times before the nucleosome rewraps. Our findings highlight the importance of the nucleosome unwrapping free energy landscape on TF occupancy and dynamics that ultimately influences transcription initiation.

Single-Molecule Human Nucleosome Spontaneously Ruptures under the Stress of Compressive Force: A New Perspective on Gene Stability and Epigenetic Pathways

Force manipulation on the biological entities from living cells to protein molecules has revealed many mechanical details of cell biology from resolving folding and unfolding pathways to finding molecular interaction forces. A nucleosome is the basic repeating unit of chromatin where the histone octamer is wrapped by DNA, important for gene stability and regulation. How the inner side of the DNA gets accessed by other DNA binding molecules has been a puzzle that has been intensively studied and debated, important to epigenetics, gene stability, and regulations. Here we report our observation of spontaneous ruptures of human nucleosomes under pico-Newton (pN) compressive force. The amplitude of the compressive force, a squeezing rather than pulling force, involved in our experiment is tens of pN, which can be thermally available by biological force fluctuation at room temperature and under physiological conditions. This kind of structural rupture can loosen up the DNA around the histone, which in turn makes the DNA accessible to transcription and epigenetic modifications.

Understanding the Target Search by Multiple Transcription Factors on Nucleosomal DNA

The association of multiple Transcription Factors (TFs) in the cis-regulatory region is imperative for developmental changes in eukaryotes. The underlying process is exceedingly complex, and it is not at all clear what orchestrates the overall search process by multiple TFs. In this study, by developing a theoretical model based on a discrete-state stochastic approach, we investigated the target search mechanism of multiple TFs on nucleosomal DNA. Experimental kinetic rate constants of different TFs are taken as input to estimate the Mean-First-Passage time to recognize the binding motifs by two TFs on a dynamic nucleosome model. The theory systematically analyzes when the TFs search their binding motifs hierarchically and when simultaneously by proceeding via the formation of a protein-protein complex. Our results, validated by extensive Monte Carlo simulations, elucidate the molecular basis of the complex target search phenomenon of multiple TFs on nucleosomal DNA.

Biologically relevant small variations of intra-cellular pH can have significant effect on stability of protein-DNA complexes, including the nucleosome

Stability of a protein-ligand complex may be sensitive to pH of its environment. Here we explore, computationally, stability of a set of protein-nucleic acid complexes using fundamental thermodynamic linkage relationship. The nucleosome, as well as an essentially random selection of 20 protein complexes with DNA or RNA, are included in the analysis. An increase in intra-cellular/intra-nuclear pH destabilizes most complexes, including the nucleosome. We propose to quantify the effect by ΔΔG^0.3-the change in the binding free energy due to pH increase of 0.3 units, corresponding to doubling of the H ⁺ activity; variations of pH of this amplitude can occur in living cells, including in the course of the cell cycle, and in cancer cells relative to normal ones. We suggest, based on relevant experimental findings, a threshold of biological significance of 1 2 k B T ( ∼ 0.3 k c a l / m o l ) for changes of stability of chromatin-related protein-DNA complexes: a change in the binding affinity above the threshold may have biological consequences. We find that for 70% of the examined complexes, Δ Δ G 0.3 > 1 2 k B T (for 10%, ΔΔG^0.3 is between 3 and 4 k _B T). Thus, small but relevant variations of intra-nuclear pH of 0.3 may have biological consequences for many protein-nucleic acid complexes. The binding affinity between the histone octamer and its DNA, which directly affects the DNA accessibility in the nucleosome, is predicted to be highly sensitive to intra-nuclear pH. A variation of 0.3 units results in ΔΔG^0.3 ∼ 10k _B T ( ∼ 6 k c a l / m o l ) ; for spontaneous unwrapping of 20 bp long entry/exit fragments of the nucleosomal DNA, ΔΔG^0.3 = 2.2k _B T; partial disassembly of the nucleosome into the tetrasome is characterized by ΔΔG^0.3 = 5.2k _B T. The predicted pH -induced modulations of the nucleosome stability are significant enough to suggest that they may have consequences relevant to the biological function of the nucleosome. Accessibility of the nucleosomal DNA is predicted to positively correlate with pH variations during the cell cycle; an increase in intra-cellular pH seen in cancer cells is predicted to lead to a more accessible nucleosomal DNA; a drop in pH associated with apoptosis is predicted to make nucleosomal DNA less accessible. We speculate that processes that depend on accessibility to the DNA in the nucleosomes, such as transcription or DNA replication, might become upregulated due to relatively small, but nevertheless realistic increases of intra-nuclear pH.

Nucleosome breathing facilitates cooperative binding of pluripotency factors Sox2 and Oct4 to DNA

Critical lineage commitment events are staged by multiple transcription factors (TFs) binding to their cognate motifs, often positioned at nucleosome-enriched regions of chromatin. The underlying mechanism remains elusive due to difficulty in disentangling the heterogeneity in chromatin states. Using a novel coarse-grained model and molecular dynamics simulations, here we probe the association of Sox2 and Oct4 proteins that show clustered binding at the entry-exit region of a nucleosome. The model captures the conformational heterogeneity of nucleosome breathing dynamics that features repeated wrap-unwrap transitions of a DNA segment from one end of the nucleosome. During the dynamics, DNA forms bulges that diffuse stochastically and may regulate the target search dynamics of a protein by nonspecifically interacting with it. The overall search kinetics of the TF pair follows a "dissociation-compensated-association" mechanism, where Oct4 binding is facilitated by the association of Sox2. The cooperativity stems from a change in entropy caused by an alteration in the nucleosome dynamics upon TF binding. The binding pattern is consistent with a live-cell single-particle tracking experiment, suggesting the mechanism observed for clustered binding of a TF pair, which is a hallmark of cis-regulatory elements, has broader implications in understanding gene regulation in a complex chromatin environment.

H2A-H2B Histone Dimer Plasticity and Its Functional Implications

The protein core of the nucleosome is composed of an H3-H4 histone tetramer and two H2A-H2B histone dimers. The tetramer organizes the central 60 DNA bp, while H2A-H2B dimers lock the flanking DNA segments. Being positioned at the sides of the nucleosome, H2A-H2B dimers stabilize the overall structure of the nucleosome and modulate its dynamics, such as DNA unwrapping, sliding, etc. Such modulation at the epigenetic level is achieved through post-translational modifications and the incorporation of histone variants. However, the detailed connection between the sequence of H2A-H2B histones and their structure, dynamics and implications for nucleosome functioning remains elusive. In this work, we present a detailed study of H2A-H2B dimer dynamics in the free form and in the context of nucleosomes via atomistic molecular dynamics simulations (based on X. laevis histones). We supplement simulation results by comparative analysis of information in the structural databases. Particularly, we describe a major dynamical mode corresponding to the bending movement of the longest H2A and H2B α-helices. This overall bending dynamics of the H2A-H2B dimer were found to be modulated by its interactions with DNA, H3-H4 tetramer, the presence of DNA twist-defects with nucleosomal DNA and the amino acid sequence of histones. Taken together, our results shed new light on the dynamical mechanisms of nucleosome functioning, such as nucleosome sliding, DNA-unwrapping and their epigenetic modulation.

Nucleosome assembly and disassembly pathways in vitro

Structural fluctuations of nucleosomes modulate the access to internal DNA in eukaryotic cells; clearly characterisation of this fundamental process is crucial to understanding gene regulation. Here we apply PhAST (Photochemical Analysis of Structural Transitions) to monitor at a base pair level, structural alterations induced all along the DNA upon histone binding or release. By offering the first reliable, detailed comparison of nucleosome assembly and disassembly in vitro, we reveal similarities and differences between the two processes. We identify multiple, sequential intermediate states characterised by specific PhAST signals whose localisation and amplitude reflect asymmetries of DNA/histone interactions with respect to the nucleosome pseudo dyad. These asymmetries involve not only the DNA extremities but also regions close to the pseudo dyad. Localisations of asymmetries develop in a consistent manner during both assembly and disassembly processes; they primarily reflect the DNA sequence effect on the efficiency of DNA-histone binding. More unexpectedly, the amplitude component of PhAST signals not only evolves as a function of intermediate states but does so differently between assembly and disassembly pathways. Our observation of differences between assembly and disassembly opens up new avenues to define the role of the DNA sequence in processes underlying the regulation of gene expression. Overall, we provide new insights into how the intrinsic properties of DNA are integrated into a holistic mechanism that controls chromatin structure.

Stochastic models of nucleosome dynamics reveal regulatory rules of stimulus-induced epigenome remodeling

The genomic positions of nucleosomes are a defining feature of the cell's epigenomic state, but signal-dependent transcription factors (SDTFs), upon activation, bind to specific genomic locations and modify nucleosome positioning. Here we leverage SDTFs as perturbation probes to learn about nucleosome dynamics in living cells. We develop Markov models of nucleosome dynamics and fit them to time course sequencing data of DNA accessibility. We find that (1) the dynamics of DNA unwrapping are significantly slower in cells than reported from cell-free experiments, (2) only models with cooperativity in wrapping and unwrapping fit the available data, (3) SDTF activity produces the highest eviction probability when its binding site is adjacent to but not on the nucleosome dyad, and (4) oscillatory SDTF activity results in high location variability. Our work uncovers the regulatory rules governing SDTF-induced nucleosome dynamics in live cells, which can predict chromatin accessibility alterations during inflammation at single-nucleosome resolution.

How Pioneer Transcription Factors Search for Target Sites on Nucleosomal DNA

All major biological processes start after protein molecules known as transcription factors detect specific regulatory sequences on DNA and initiate genetic expression by associating to them. But in eukaryotic cells, much of the DNA is covered by nucleosomes and other chromatin structures, preventing transcription factors from binding to their targets. At the same time, experimental studies show that there are several classes of proteins, called "pioneer transcription factors", that are able to reach the targets on nucleosomal DNA; however, the underlying microscopic mechanisms remain not well understood. We propose a new theoretical approach that might explain how pioneer transcription factors can find their targets. It is argued that pioneer transcription factors might weaken the interactions between the DNA and nucleosome by substituting them with similar interactions between transcription factors and DNA. Using this idea, we develop a discrete-state stochastic model that allows for exact calculations of target search dynamics on nucleosomal DNA using first-passage probabilities approach. It is found that the target search on nuclesomal DNA for pioneer transcription factors might be significantly accelerated while the search is slower on naked DNA in comparison with normal transcription factors. Our theoretical predictions are supported by Monte Carlo computer simulations, and they also agree with available experimental observations.

The landscape of pioneer factor activity reveals the mechanisms of chromatin reprogramming and genome activation

Upon fertilization, embryos undergo chromatin reprogramming and genome activation; however, the mechanisms that regulate these processes are poorly understood. Here, we generated a triple mutant for Nanog, Pou5f3, and Sox19b (NPS) in zebrafish and found that NPS pioneer chromatin opening at >50% of active enhancers. NPS regulate acetylation across core histones at enhancers and promoters, and their function in gene activation can be bypassed by recruiting histone acetyltransferase to individual genes. NPS pioneer chromatin opening individually, redundantly, or additively depending on sequence context, and we show that high nucleosome occupancy facilitates NPS pioneering activity. Nucleosome position varies based on the input of different transcription factors (TFs), providing a flexible platform to modulate pioneering activity. Altogether, our results illuminate the sequence of events during genome activation and offer a conceptual framework to understand how pioneer factors interpret the genome and integrate different TF inputs across cell types and developmental transitions.

Kinetic origin of nucleosome invasion by pioneer transcription factors

Recently, a cryo-electron microscopy study has captured different stages of nucleosome breathing dynamics that show partial unwrapping of DNA from histone core to permit transient access to the DNA sites by transcription factors. In practice, however, only a subset of transcription factors named pioneer factors can invade nucleosomes and bind to specific DNA sites to trigger essential DNA metabolic processes. We propose a discrete-state stochastic model that considers the interplay of nucleosome breathing and protein dynamics explicitly and estimate the mean time to search the target DNA sites. It is found that the molecular principle governing the search process on nucleosome is very different compared to that on naked DNA. The pioneer factors minimize their search times on nucleosomal DNA by compensating their nucleosome association rates by dissociation rates. A fine balance between the two presents a tradeoff between their nuclear mobility and error associated with the search process.

Single molecule characterization of the binding kinetics of a transcription factor and its modulation by DNA sequence and methylation

The interaction of transcription factors with their response elements in DNA is emerging as a highly complex process, whose characterization requires measuring the full distribution of binding and dissociation times in a well-controlled assay. Here, we present a single-molecule assay that exploits the thermal fluctuations of a DNA hairpin to detect the association and dissociation of individual, unlabeled transcription factors. We demonstrate this new approach by following the binding of Egr1 to its consensus motif and the three binding sites found in the promoter of the Lhb gene, and find that both association and dissociation are modulated by the 9 bp core motif and the sequences around it. In addition, CpG methylation modulates the dissociation kinetics in a sequence and position-dependent manner, which can both stabilize or destabilize the complex. Together, our findings show how variations in sequence and methylation patterns synergistically extend the spectrum of a protein's binding properties, and demonstrate how the proposed approach can provide new insights on the function of transcription factors.

NF-κB dynamics determine the stimulus specificity of epigenomic reprogramming in macrophages

The epigenome of macrophages can be reprogrammed by extracellular cues, but the extent to which different stimuli achieve this is unclear. Nuclear factor κB (NF-κB) is a transcription factor that is activated by all pathogen-associated stimuli and can reprogram the epigenome by activating latent enhancers. However, we show that NF-κB does so only in response to a subset of stimuli. This stimulus specificity depends on the temporal dynamics of NF-κB activity, in particular whether it is oscillatory or non-oscillatory. Non-oscillatory NF-κB opens chromatin by sustained disruption of nucleosomal histone-DNA interactions, enabling activation of latent enhancers that modulate expression of immune response genes. Thus, temporal dynamics can determine a transcription factor's capacity to reprogram the epigenome in a stimulus-specific manner.

The Dynamic Influence of Linker Histone Saturation within the Poly-Nucleosome Array

Linker histones bind to nucleosomes and modify chromatin structure and dynamics as a means of epigenetic regulation. Biophysical studies have shown that chromatin fibers can adopt a plethora of conformations with varying levels of compaction. Linker histone condensation, and its specific binding disposition, has been associated with directly tuning this ensemble of states. However, the atomistic dynamics and quantification of this mechanism remains poorly understood. Here, we present molecular dynamics simulations of octa-nucleosome arrays, based on a cryo-EM structure of the 30-nm chromatin fiber, with and without the globular domains of the H1 linker histone to determine how they influence fiber structures and dynamics. Results show that when bound, linker histones inhibit DNA flexibility and stabilize repeating tetra-nucleosomal units, giving rise to increased chromatin compaction. Furthermore, upon the removal of H1, there is a significant destabilization of this compact structure as the fiber adopts less strained and untwisted states. Interestingly, linker DNA sampling in the octa-nucleosome is exaggerated compared to its mono-nucleosome counterparts, suggesting that chromatin architecture plays a significant role in DNA strain even in the absence of linker histones. Moreover, H1-bound states are shown to have increased stiffness within tetra-nucleosomes, but not between them. This increased stiffness leads to stronger long-range correlations within the fiber, which may result in the propagation of epigenetic signals over longer spatial ranges. These simulations highlight the effects of linker histone binding on the internal dynamics and global structure of poly-nucleosome arrays, while providing physical insight into a mechanism of chromatin compaction.

Inhibition of CRISPR-Cas12a DNA targeting by nucleosomes and chromatin

Genome engineering nucleases must access chromatinized DNA. Here, we investigate how AsCas12a cleaves DNA within human nucleosomes and phase-condensed nucleosome arrays. Using quantitative kinetics approaches, we show that dynamic nucleosome unwrapping regulates target accessibility to Cas12a and determines the extent to which both steps of binding-PAM recognition and R-loop formation-are inhibited by the nucleosome. Relaxing DNA wrapping within the nucleosome by reducing DNA bendability, adding histone modifications, or introducing target-proximal dCas9 enhances DNA cleavage rates over 10-fold. Unexpectedly, Cas12a readily cleaves internucleosomal linker DNA within chromatin-like, phase-separated nucleosome arrays. DNA targeting is reduced only ~5-fold due to neighboring nucleosomes and chromatin compaction. This work explains the observation that on-target cleavage within nucleosomes occurs less often than off-target cleavage within nucleosome-depleted genomic regions in cells. We conclude that nucleosome unwrapping regulates accessibility to CRISPR-Cas nucleases and propose that increasing nucleosome breathing dynamics will improve DNA targeting in eukaryotic cells.

Linker DNA Length is a Key to Tri-nucleosome Folding

The folding of a nucleosome array has long been one of the fundamental and unsolved problems in chromatin biology. In this study, we address how nucleosome array folding depends on the length of linker DNA. We performed molecular dynamics simulations of a tri-nucleosome, a minimal model of chromatin folding, with various linker lengths (LLs) ranging from 20 to 40 base pairs (bps). We found that the tri-nucleosome folding strongly depends on LLs, and classified the structure ensemble into five classes, named from trinuc-1 to trinuc-5. As a function of LL, the different classes appear, on average, every 2 bps with a period of 10 bps, and are characterized by distinct inter-nucleosome interactions. The trinuc-1 conformation corresponds to LL ~ 10n, where n is an integer, and is stabilized by the tight packing between the first and the third nucleosomes, consistent with a zigzag fiber form. Structures of the other four classes are more diverse and distributed continuously in the space of possible configurations. Histone-DNA electrostatic interactions in the tri-nucleosome are further analyzed.

The role of histone tails in nucleosome stability: An electrostatic perspective

We propose a methodology for the study of protein-DNA electrostatic interactions and apply it to clarify the effect of histone tails in nucleosomes. This method can be used to correlate electrostatic interactions to structural and functional features of protein-DNA systems, and can be combined with coarse-grained representations. In particular, we focus on the electrostatic field and resulting forces acting on the DNA. We investigate the electrostatic origins of effects such as different stages in DNA unwrapping, nucleosome destabilization upon histone tail truncation, and the role of specific arginines and lysines undergoing Post-Translational Modifications. We find that the positioning of the histone tails can oppose the attractive pull of the histone core, locally deform the DNA, and tune DNA unwrapping. Small conformational variations in the often overlooked H2A C-terminal tails had significant electrostatic repercussions near the DNA entry and exit sites. The H2A N-terminal tail exerts attractive electrostatic forces towards the histone core in positions where Polymerase II halts its progress. We validate our results with comparisons to previous experimental and computational observations.

Nucleosome allostery in pioneer transcription factor binding

While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POU_S domain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, the LIN28B nucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.

Transient Kinetic Analysis of SWR1C-Catalyzed H2A.Z Deposition Unravels the Impact of Nucleosome Dynamics and the Asymmetry of Histone Exchange

The SWR1C chromatin remodeling enzyme catalyzes ATP-dependent replacement of nucleosomal H2A with the H2A.Z variant, regulating key DNA-mediated processes such as transcription and DNA repair. Here, we investigate the transient kinetic mechanism of the histone exchange reaction, employing ensemble FRET, fluorescence correlation spectroscopy (FCS), and the steady-state kinetics of ATP hydrolysis. Our studies indicate that SWR1C modulates nucleosome dynamics on both the millisecond and microsecond timescales, poising the nucleosome for the dimer exchange reaction. The transient kinetic analysis of the remodeling reaction performed under single turnover conditions unraveled a striking asymmetry in the ATP-dependent replacement of nucleosomal dimers, promoted by localized DNA unwrapping. Taken together, our transient kinetic studies identify intermediates and provide crucial insights into the SWR1C-catalyzed dimer exchange reaction and shed light on how the mechanics of H2A.Z deposition might contribute to transcriptional regulation in vivo.

The SWR1C remodeling enzyme catalyzes ATP-dependent replacement of nucleosomal H2A with the H2A.Z variant at promoter-proximal nucleosomes. Singh et al. investigate the transient kinetic mechanism of this histone exchange reaction and show that SWR1C transiently unwraps nucleosomal DNA, promoting a concerted dimer eviction and replacement reaction.

Nucleosomes as allosteric scaffolds for genetic regulation

Nucleosomes are stable yet highly dynamic complexes exhibiting diverse types of motions, such as sliding, DNA unwrapping, and disassembly, encoding a landscape with a large number of metastable states. In this review, describing recent studies on these nucleosome structure changes, we propose that the nucleosome can be viewed as an ideal allosteric scaffold: regulated by effector molecules such as transcription factors and chromatin remodelers, the nucleosome controls the downstream gene activity. Binding of transcription factors to the nucleosome can enhance DNA unwrapping or slide the DNA, altering either the binding or the unbinding of other transcription factors to nearby sites. ATP-dependent chromatin remodelers induce a series of DNA deformations, which allosterically propagate throughout the nucleosome to induce DNA sliding or histone exchange.

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

"Direct" and "Indirect" Effects of Histone Modifications: Modulation of Sterical Bulk as a Novel Source of Functionality

The chromatin-regulatory principles of histone post-translational modifications (PTMs) are discussed with a focus on the potential alterations in chromatin functional state due to steric and mechanical constraints imposed by bulky histone modifications such as ubiquitin and SUMO. In the classical view, PTMs operate as recruitment platforms for histone "readers," and as determinants of chromatin array compaction. Alterations of histone charges by "small" chemical modifications (e.g., acetylation, phosphorylation) could regulate nucleosome spontaneous dynamics without globally affecting nucleosome structure. These fluctuations in nucleosome wrapping can be exploited by chromatin-processing machinery. In contrast, ubiquitin and SUMO are comparable in size to histones, and it seems logical that these PTMs could conflict with canonical nucleosome organization. An experimentally testable hypothesis that by adding sterical bulk these PTMs can robustly alter nucleosome primary structure is proposed. The model presented here stresses the diversity of mechanisms by which histone PTMs regulate chromatin dynamics, primary structure and, hence, functionality.

Acetylation of the histone H3 tail domain regulates base excision repair on higher-order chromatin structures

Despite recent evidence suggesting that histone lysine acetylation contributes to base excision repair (BER) in cells, their exact mechanistic role remains unclear. In order to examine the influence of histone acetylation on the initial steps of BER, we assembled nucleosome arrays consisting of homogeneously acetylated histone H3 (H3K18 and H3K27) and measured the repair of a site-specifically positioned 2′-deoxyuridine (dU) residue by uracil DNA glycosylase (UDG) and apurinic/apyrimidinic endonuclease 1 (APE1). We find that H3K18ac and H3K27ac differentially influence the combined activities of UDG/APE1 on compact chromatin, suggesting that acetylated lysine residues on the H3 tail domain play distinct roles in regulating the initial steps of BER. In addition, we show that the effects of H3 tail domain acetylation on UDG/APE1 activity are at the nucleosome level and do not influence higher-order chromatin folding. Overall, these results establish a novel regulatory role for histone H3 acetylation during the initiation of BER on chromatin.

Positional Dependence of DNA Hole Transfer Efficiency in Nucleosome Core Particles

Electron deficient "holes" migrate over long distances through the π-system in free DNA. Hole transfer efficiency (HTE) is strongly dependent on sequence and π-stacking. However, there is no consensus regarding the effects of nucleosome core particle (NCP) environment on hole migration. We quantitatively determined HTE in free DNA and NCPs by independently generating holes at specific positions in DNA. The relative HTE varied widely with respect to position within the NCP and proximity to tyrosine, which suppresses hole transfer. These data indicate that hole transfer in chromatin will be affected by the DNA sequence and its position with respect to histone proteins within NCPs.

1CPN: A coarse-grained multi-scale model of chromatin

A central question in epigenetics is how histone modifications influence the 3D structure of eukaryotic genomes and, ultimately, how this 3D structure is manifested in gene expression. The wide range of length scales that influence the 3D genome structure presents important challenges; epigenetic modifications to histones occur on scales of angstroms, yet the resulting effects of these modifications on genome structure can span micrometers. There is a scarcity of computational tools capable of providing a mechanistic picture of how molecular information from individual histones is propagated up to large regions of the genome. In this work, a new molecular model of chromatin is presented that provides such a picture. This new model, referred to as 1CPN, is structured around a rigorous multiscale approach, whereby free energies from an established and extensively validated model of the nucleosome are mapped onto a reduced coarse-grained topology. As such, 1CPN incorporates detailed physics from the nucleosome, such as histone modifications and DNA sequence, while maintaining the computational efficiency that is required to permit kilobase-scale simulations of genomic DNA. The 1CPN model reproduces the free energies and dynamics of both single nucleosomes and short chromatin fibers, and it is shown to be compatible with recently developed models of the linker histone. It is applied here to examine the effects of the linker DNA on the free energies of chromatin assembly and to demonstrate that these free energies are strongly dependent on the linker DNA length, pitch, and even DNA sequence. The 1CPN model is implemented in the LAMMPS simulation package and is distributed freely for public use.

Histone H2A Variants Enhance the Initiation of Base Excision Repair in Nucleosomes

Substituting histone variants for their canonical counterparts can profoundly alter chromatin structure, thereby impacting multiple biological processes. Here, we investigate the influence of histone variants from the H2A family on the excision of uracil (U) by the base excision repair (BER) enzymes uracil DNA glycosylase (UDG) and single-strand selective monofunctional uracil DNA glycosylase. Using a DNA population with globally distributed U:G base pairs, enhanced excision is observed in H2A.Z and macroH2A-containing nucleosome core particles (NCPs). The U with reduced solution accessibility exhibit limited UDG activity in canonical NCPs but are more readily excised in variant NCPs, reflecting the ability of these variants to facilitate excision at sites that are otherwise poorly repaired. We also find that U with the largest increase in the level of excision in variant NCPs are clustered in regions with differential structural features between the variants and canonical H2A. Within 35-40 bp of the DNA terminus in macroH2A NCPs, the activities of both glycosylases are comparable to that on the free duplex. We show that this high level of activity results from two distinct species within the macroH2A NCP ensemble: octasomes and hexasomes. These observations reveal potential functions for H2A variants in promoting BER and preventing mutagenesis within the context of chromatin.

Biophysics of Chromatin Dynamics

Nucleosomes and chromatin control eukaryotic genome accessibility and thereby regulate DNA processes, including transcription, replication, and repair. Conformational dynamics within the nucleosome and chromatin structure play a key role in this regulatory function. Structural fluctuations continuously expose internal DNA sequences and nucleosome surfaces, thereby providing transient access for the nuclear machinery. Progress in structural studies of nucleosomes and chromatin has provided detailed insight into local chromatin organization and has set the stage for recent in-depth investigations of the structural dynamics of nucleosomes and chromatin fibers. Here, we discuss the dynamic processes observed in chromatin over different length scales and timescales and review current knowledge about the biophysics of distinct structural transitions.

Dissociation rate compensation mechanism for budding yeast pioneer transcription factors

Nucleosomes restrict the occupancy of most transcription factors (TF) by reducing binding and accelerating dissociation, while a small group of TFs have high affinities to nucleosome-embedded sites and facilitate nucleosome displacement. To understand this process mechanistically, we investigated two Saccharomyces cerevisiae TFs, Reb1 and Cbf1. We show that these factors bind to their sites within nucleosomes with similar binding affinities as to naked DNA, trapping a partially unwrapped nucleosome without histone eviction. Both the binding and dissociation rates of Reb1 and Cbf1 are significantly slower at the nucleosomal sites relative to those for naked DNA, demonstrating that the high affinities are achieved by increasing the dwell time on nucleosomes in order to compensate for reduced binding. Reb1 also shows slow migration rate in the yeast nuclei. These properties are similar to those of human pioneer factors (PFs), suggesting that the mechanism of nucleosome targeting is conserved from yeast to humans.

Nucleosome DNA unwrapping does not affect prototype foamy virus integration efficiency or site selection

Eukaryotic DNA binding proteins must access genomic DNA that is packaged into chromatin in vivo. During a productive infection, retroviral integrases (IN) must similarly interact with chromatin to integrate the viral cDNA genome. Here we examine the role of nucleosome DNA unwrapping in the retroviral integrase search for a target site. These studies utilized PFV intasomes that are comprised of a tetramer of PFV IN with two oligomers mimicking the viral cDNA ends. Modified recombinant human histones were used to generate nucleosomes with increased unwrapping rates at different DNA regions. These modifications included the acetylmimetic H3(K56Q) and the chemically engineered H4(K77ac, K79ac). While transcription factors and DNA damage sensors may search nucleosome bound DNA during transient unwrapping, PFV intasome mediated integration appears to be unaffected by increased nucleosome unwrapping. These studies suggest PFV intasomes do not utilize nucleosome unwrapping to search nucleosome targets.

Extracting collective motions underlying nucleosome dynamics via nonlinear manifold learning

The identification of effective collective variables remains a challenge in molecular simulations of complex systems. Here, we use a nonlinear manifold learning technique known as the diffusion map to extract key dynamical motions from a complex biomolecular system known as the nucleosome: a DNA-protein complex consisting of a DNA segment wrapped around a disc-shaped group of eight histone proteins. We show that without any a priori information, diffusion maps can identify and extract meaningful collective variables that characterize the motion of the nucleosome complex. We find excellent agreement between the collective variables identified by the diffusion map and those obtained manually using a free energy-based analysis. Notably, diffusion maps are shown to also identify subtle features of nucleosome dynamics that did not appear in those manually specified collective variables. For example, diffusion maps identify the importance of looped conformations in which DNA bulges away from the histone complex that are important for the motion of DNA around the nucleosome. This work demonstrates that diffusion maps can be a promising tool for analyzing very large molecular systems and for identifying their characteristic slow modes.

Molecular Mechanism of Spontaneous Nucleosome Unraveling

Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.

High precision FRET studies reveal reversible transitions in nucleosomes between microseconds and minutes

Nucleosomes play a dual role in compacting the genome and regulating the access to DNA. To unravel the underlying mechanism, we study fluorescently labeled mononucleosomes by multi-parameter FRET measurements and characterize their structural and dynamic heterogeneity upon NaCl-induced destabilization. Species-selective fluorescence lifetime analysis and dynamic photon distribution analysis reveal intermediates during nucleosome opening and lead to a coherent structural and kinetic model. In dynamic octasomes and hexasomes the interface between the H2A-H2B dimers and the (H3-H4)₂ tetramer opens asymmetrically by an angle of ≈20° on a 50 and 15 µs time scale, respectively. This is followed by a slower stepwise release of the dimers coupled with DNA unwrapping. A mutation (H2A-R81A) at the interface between H2A and H3 facilitates initial opening, confirming the central role of the dimer:tetramer interface for nucleosome stability. Partially opened states such as those described here might serve as convenient nucleation sites for DNA-recognizing proteins.

Plug-and-Play Approach for Preparing Chromatin Containing Site-Specific DNA Modifications: The Influence of Chromatin Structure on Base Excision Repair

The genomic DNA of eukaryotic cells exists in the form of chromatin, the structure of which controls the biochemical accessibility of the underlying DNA to effector proteins. In order to gain an in depth molecular understanding of how chromatin structure regulates DNA repair, detailed in vitro biochemical and biophysical studies are required. However, because of challenges associated with reconstituting nucleosome arrays containing site-specifically positioned DNA modifications, such studies have been limited to the use of mono- and dinucleosomes as model in vitro substrates, which are incapable of folding into native chromatin structures. To address this issue, we developed a straightforward and general approach for assembling chemically defined oligonucleosome arrays (i.e., designer chromatin) containing site-specifically modified DNA. Our method takes advantage of nicking endonucleases to excise short fragments of unmodified DNA, which are subsequently replaced with synthetic oligonucleotides containing the desired modification. Using this approach, we prepared several oligonucleosome substrates containing precisely positioned 2'-deoxyuridine (dU) residues and examined the efficiency of base excision repair (BER) within several distinct chromatin architectures. We show that, depending on the translational position of the lesion, the combined catalytic activities of uracil DNA glycosylase (UDG) and apurinic/apyrimidinic endonuclease 1 (APE1) can be either inhibited by as much as 20-fold or accelerated by more than 5-fold within compact chromatin (i.e., the 30 nm fiber) relative to naked DNA. Moreover, we demonstrate that digestion of dU by UDG/APE1 proceeds much more rapidly in mononucleosomes than in compacted nucleosome arrays, thereby providing the first direct evidence that internucleosome interactions play an important role in regulating BER within higher-order chromatin structures. Overall, this work highlights the value of performing detailed biochemical studies on precisely modified chromatin substrates in vitro and provides a robust platform for investigating DNA modifications in chromatin biology.

A cassette of basic amino acids in histone H2B regulates nucleosome dynamics and access to DNA damage

Nucleosome dynamics, such as spontaneous DNA unwrapping, are postulated to have a critical role in regulating the access of DNA repair machinery to DNA lesions within nucleosomes. However, the specific histone domains that regulate nucleosome dynamics and the impact of such changes in intrinsic nucleosome dynamics on DNA repair are not well understood. Previous studies identified a highly conserved region in the N-terminal tail of histone H2B known as the histone H2Brepression (or HBR) domain, which has a significant influence on gene expression, chromatin assembly, and DNA damage formation and repair. However, the molecular mechanism(s) that may account for these observations are limited. In this study, we characterized the stability and dynamics of ΔHBR mutant nucleosome core particles (NCPs) in vitro by restriction enzyme accessibility (REA), FRET, and temperature-induced sliding of histone octamers. Our results indicate that ΔHBR-NCPs are more dynamic, with a larger steady-state fraction of the NCP population occupying the unwrapped state than for WT-NCPs. Additionally, ΔHBR-histone octamers are more susceptible to temperature-induced sliding on DNA than WT histone octamers. Furthermore, we show that the activity of base excision repair enzymes at uracil lesions and single nucleotide gaps is enhanced in a site-specific manner in ΔHBR-NCPs. This enhanced activity correlates well with regions exhibiting increased DNA unwrapping. Finally, removal of the HBR domain is not sufficient to completely alleviate the structural constraints imposed by histone octamers on the activity of base excision repair enzymes.

Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core

BACKGROUND

Controlled modulation of nucleosomal DNA accessibility via post-translational modifications (PTM) is a critical component to many cellular functions. Charge-altering PTMs in the globular histone core-including acetylation, phosphorylation, crotonylation, propionylation, butyrylation, formylation, and citrullination-can alter the strong electrostatic interactions between the oppositely charged nucleosomal DNA and the histone proteins and thus modulate accessibility of the nucleosomal DNA, affecting processes that depend on access to the genetic information, such as transcription. However, direct experimental investigation of the effects of these PTMs is very difficult. Theoretical models can rationalize existing observations, suggest working hypotheses for future experiments, and provide a unifying framework for connecting PTMs with the observed effects.

RESULTS

A physics-based framework is proposed that predicts the effect of charge-altering PTMs in the histone core, quantitatively for several types of lysine charge-neutralizing PTMs including acetylation, and qualitatively for all phosphorylations, on the nucleosome stability and subsequent changes in DNA accessibility, making a connection to resulting biological phenotypes. The framework takes into account multiple partially assembled states of the nucleosome at the atomic resolution. The framework is validated against experimentally known nucleosome stability changes due to the acetylation of specific lysines, and their effect on transcription. The predicted effect of charge-altering PTMs on DNA accessibility can vary dramatically, from virtually none to a strong, region-dependent increase in accessibility of the nucleosomal DNA; in some cases, e.g., H4K44, H2AK75, and H2BK57, the effect is significantly stronger than that of the extensively studied acetylation sites such H3K56, H3K115 or H3K122. Proximity to the DNA is suggestive of the strength of the PTM effect, but there are many exceptions. For the vast majority of charge-altering PTMs, the predicted increase in the DNA accessibility should be large enough to result in a measurable modulation of transcription. However, a few possible PTMs, such as acetylation of H4K77, counterintuitively decrease the DNA accessibility, suggestive of the repressed chromatin. A structural explanation for the phenomenon is provided. For the majority of charge-altering PTMs, the effect on DNA accessibility is simply additive (noncooperative), but there are exceptions, e.g., simultaneous acetylation of H4K79 and H3K122, where the combined effect is amplified. The amplification is a direct consequence of the nucleosome-DNA complex having more than two structural states. The effect of individual PTMs is classified based on changes in the accessibility of various regions throughout the nucleosomal DNA. The PTM's resulting imprint on the DNA accessibility, "PTMprint," is used to predict effects of many yet unexplored PTMs. For example, acetylation of H4K44 yields a PTMprint similar to the PTMprint of H3K56, and thus acetylation of H4K44 is predicted to lead to a wide range of strong biological effects.

CONCLUSION

Charge-altering post-translational modifications in the relatively unexplored globular histone core may provide a precision mechanism for controlling accessibility to the nucleosomal DNA.

Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA

The eukaryotic genome is packaged into a nucleus in the form of chromatin. The fundamental structural unit of chromatin is a protein-DNA complex, the nucleosome, where 146 or 147 base pairs of DNA wrap 1.75 times around a histone core. To function in cellular processes, however, nucleosomal DNA must be unwrapped. Although this unwrapping has been experimentally investigated, details of the process at an atomic level are not yet well understood. Here, we used molecular dynamics simulation with an enhanced sampling method to calculate the free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA. A free energy change of about 11.5 kcal/mol for the unwrapping agrees well with values obtained in single molecule experiments. This simulation revealed a variety of conformational states, indicating there are many potential paths to outer superhelicdal turn unwrapping, but the dominant path is likely asymmetric. At one end of the DNA, the first five bps unwrap, after which a second five bps unwrap at the same end with no increase in free energy. The unwrapping then starts at the other end of the DNA, where 10 bps are unwrapped. During further unwrapping of 15 bps, the unwrapping advances at one of the ends, after which the other end of the DNA unwraps to complete the unwrapping of the outer superhelical turn. These results provide insight into the construction, disruption, and repositioning of nucleosomes, which are continuously ongoing during cellular processes.

Multifaceted Targeting of the Chromatin Mediates Gonadotropin-Releasing Hormone Effects on Gene Expression in the Gonadotrope

Gonadotropin-releasing hormone (GnRH) stimulates the expression of multiple genes in the pituitary gonadotropes, most notably to induce synthesis of the gonadotropins, luteinizing hormone (LH), and follicle-stimulating hormone (FSH), but also to ensure the appropriate functioning of these cells at the center of the mammalian reproductive endocrine axis. Aside from the activation of gene-specific transcription factors, GnRH stimulates through its membrane-bound receptor, alterations in the chromatin that facilitate transcription of its target genes. These include changes in the histone and DNA modifications, nucleosome positioning, and chromatin packaging at the regulatory regions of each gene. The requirements for each of these events vary according to the DNA sequence which determines the basal chromatin packaging at the regulatory regions. Despite considerable progress in this field in recent years, we are only beginning to understand some of the complexities involved in the role and regulation of this chromatin structure, including new modifications, extensive cross talk, histone variants, and the actions of distal enhancers and non-coding RNAs. This short review aims to integrate the latest findings on GnRH-induced alterations in the chromatin of its target genes, which indicate multiple and diverse actions. Understanding these processes is illuminating not only in the context of the activation of these hormones during the reproductive life span but may also reveal how aberrant epigenetic regulation of these genes leads to sub-fertility.

The role of DNA sequence in nucleosome breathing

Roughly 3/4 of human genomes are sequestered by nucleosomes, DNA spools with a protein core, dictating a broad range of biological processes, ranging from gene regulation, recombination, and replication, to chromosome condensation. Nucleosomes are dynamical structures and temporarily expose wrapped DNA through spontaneous unspooling from either end, a process called site exposure or nucleosome breathing. Here we ask how this process is influenced by the mechanical properties of the wrapped DNA, which is known to depend on the underlying base pair sequence. Using a coarse-grained nucleosome model we calculate the accessibility profiles for site exposure. We find that the process is very sensitive to sequence effects, so that evolution could potentially tune the accessibility of nucleosomal DNA and would only need a small number of mutations to do so.

In silico evidence for sequence-dependent nucleosome sliding

Nucleosomes represent the basic building block of chromatin and provide an important mechanism by which cellular processes are controlled. The locations of nucleosomes across the genome are not random but instead depend on both the underlying DNA sequence and the dynamic action of other proteins within the nucleus. These processes are central to cellular function, and the molecular details of the interplay between DNA sequence and nucleosome dynamics remain poorly understood. In this work, we investigate this interplay in detail by relying on a molecular model, which permits development of a comprehensive picture of the underlying free energy surfaces and the corresponding dynamics of nucleosome repositioning. The mechanism of nucleosome repositioning is shown to be strongly linked to DNA sequence and directly related to the binding energy of a given DNA sequence to the histone core. It is also demonstrated that chromatin remodelers can override DNA-sequence preferences by exerting torque, and the histone H4 tail is then identified as a key component by which DNA-sequence, histone modifications, and chromatin remodelers could in fact be coupled.