A quantitative model of nucleosome dynamics

Nanoscale Structure, Interactions, and Dynamics of Centromere Nucleosomes

Centromeres are specific segments of chromosomes comprising two types of nucleosomes: canonical nucleosomes containing an octamer of H2A, H2B, H3, and H4 histones and CENP-A nucleosomes in which H3 is replaced with its analogue CENP-A. This modification leads to a difference in DNA wrapping (∼121 bp), considerably less than 147 bp in canonical nucleosomes. We used atomic force microscopy (AFM) and high-speed AFM (HS-AFM) to characterize nanoscale features and dynamics for both types of nucleosomes. For both nucleosomes, spontaneous asymmetric unwrapping of DNA was observed, and this process occurs via a transient state with ∼100 bp DNA wrapped around the core, followed by a rapid dissociation of DNA. Additionally, HS-AFM revealed higher stability of CENP-A nucleosomes compared with H3 nucleosomes in which dissociation of the histone core occurs prior to the nucleosome dissociation. These results help elucidate the differences between these nucleosomes and the potential biological necessity for CENP-A nucleosomes.

Explicit ion modeling predicts physicochemical interactions for chromatin organization

Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 kBT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono- and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 kBT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.

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.

Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation

We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solutions can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.

The Free Energy of Nucleosomal DNA Based on the Landau Model and Topology

The free energy of nucleosomal DNA plays a key role in the formation of nucleosomes in eukaryotes. Some work on the free energy of nucleosomal DNA have been carried out in experiments. However, the relationships between the free energy of nucleosomal DNA and its conformation, especially its topology, remain unclear in theory. By combining the Landau theory, the Hopfion model and experimental data, we find that the free energy of nucleosomal DNA is at the lower level. With the help of the energy minimum principle, we conclude that nucleosomal DNA stays in a stable state. Moreover, we discover that small perturbations on nucleosomal DNA have little effect on its free energy. This implies that nucleosomal DNA has a certain redundancy in order to stay stable. This explains why nucleosomal DNA will not change significantly due to small perturbations.

Explicit Ion Modeling Predicts Physicochemical Interactions for Chromatin Organization

Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 kBT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 kBT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.

Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation

We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solution can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.

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.

Studies of the Mechanism of Nucleosome Dynamics: A Review on Multifactorial Regulation from Computational and Experimental Cases

The nucleosome, which organizes the long coil of genomic DNA in a highly condensed, polymeric way, is thought to be the basic unit of chromosomal structure. As the most important protein–DNA complex, its structural and dynamic features have been successively revealed in recent years. However, its regulatory mechanism, which is modulated by multiple factors, still requires systemic discussion. This study summarizes the regulatory factors of the nucleosome’s dynamic features from the perspective of histone modification, DNA methylation, and the nucleosome-interacting factors (transcription factors and nucleosome-remodeling proteins and cations) and focuses on the research exploring the molecular mechanism through both computational and experimental approaches. The regulatory factors that affect the dynamic features of nucleosomes are also discussed in detail, such as unwrapping, wrapping, sliding, and stacking. Due to the complexity of the high-order topological structures of nucleosomes and the comprehensive effects of regulatory factors, the research on the functional modulation mechanism of nucleosomes has encountered great challenges. The integration of computational and experimental approaches, the construction of physical modes for nucleosomes, and the application of deep learning techniques will provide promising opportunities for further exploration.

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.

A nanoscale DNA force spectrometer capable of applying tension and compression on biomolecules

Single molecule force spectroscopy is a powerful approach to probe the structure, conformational changes, and kinetic properties of biological and synthetic macromolecules. However, common approaches to apply forces to biomolecules require expensive and cumbersome equipment and relatively large probes such as beads or cantilevers, which limits their use for many environments and makes integrating with other methods challenging. Furthermore, existing methods have key limitations such as an inability to apply compressive forces on single molecules. We report a nanoscale DNA force spectrometer (nDFS), which is based on a DNA origami hinge with tunable mechanical and dynamic properties. The angular free energy landscape of the nDFS can be engineered across a wide range through substitution of less than 5% of the strand components. We further incorporate a removable strut that enables reversible toggling of the nDFS between open and closed states to allow for actuated application of tensile and compressive forces. We demonstrate the ability to apply compressive forces by inducing a large bend in a 249bp DNA molecule, and tensile forces by inducing DNA unwrapping of a nucleosome sample. These results establish a versatile tool for force spectroscopy and robust methods for designing nanoscale mechanical devices with tunable force application.

Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions

Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS.

DNA sequence influences hexasome orientation to regulate DNA accessibility

Nucleosomes, the fundamental organizing units of eukaryotic genomes, contain ∼146 base pairs of DNA wrapped around a histone H3–H4 tetramer and two histone H2A–H2B dimers. Converting nucleosomes into hexasomes by removal of a H2A–H2B dimer is an important regulatory event, but its regulation and functional consequences are not well-understood. To investigate the influence of hexasomes on DNA accessibility, we used the property of the Widom-601 Nucleosome Positioning Sequence (NPS) to form homogeneously oriented hexasomes in vitro. We find that DNA accessibility to transcription factors (TF) on the hexasome H2A–H2B distal side is identical to naked DNA, while the accessibility on the H2A–H2B proximal side is reduced by 2-fold, which is due to a 2-fold reduction in hexasome unwrapping probability. We then determined that a 23 bp region of the Widom-601 NPS is responsible for forming homogeneously oriented hexasomes. Analysis of published ChIP-exo data of hexasome containing genes identified two DNA sequence motifs that correlate with hexasome orientation in vivo, while ExoIII mapping studies of these sequences revealed they generate homogeneously oriented hexasomes in vitro. These results indicate that hexasome orientation, which is influenced by the underlying DNA sequence in vivo, is important for modulating DNA accessibility to regulate transcription.

Quantitative Modeling of Nucleosome Unwrapping from Both Ends

In eukaryotic cells, DNA is packaged into chromatin where nucleosomes are the basic packaging unit. Important cellular processes including gene expression, DNA replication, and DNA repair require nucleosomal DNA to be unwrapped so that functional proteins can access their target sites, which otherwise are sterically occluded. A key question in this process is what the unwrapped conformations individual nucleosomes adopt within chromatin are. Here, we develop a concurrent nucleosome unwrapping model to address this question. We hypothesize that for a given end-to-end distance of the nucleosomal DNA, the nucleosomal DNA stochastically unwraps from the histone core from both ends independently and that this combination of unwrapping from both sides results in a significant increase in the average distance between the DNA extending from both sides of the nucleosomes. We test our model on recently published experiments using a DNA origami nanocaliper that quantifies nucleosome unwrapping and achieve good agreement between experiment and model prediction. We then investigate the DNA origami caliper distribution when attached to a hexasome (a nucleosome lacking an H2A/H2B dimer). A significant shift in the caliper angle distribution caused by the asymmetric structural features of the hexasome seen experimentally is consistent with the model. Our modeling approach may be more broadly useful to the interpretation of other studies of nucleosome dynamics, chromatin dynamics, and regulatory processes involving nucleosome unwrapping, as well as more generally to optimization of future DNA origami designs to probe mechanical properties of biomolecules.

Critical role of histone tail entropy in nucleosome unwinding

The nucleosome is the fundamental packaging unit for the genome. It must remain tightly wound to ensure genome stability while simultaneously being flexible enough to keep the DNA molecule accessible for genome function. The set of physicochemical interactions responsible for the delicate balance between these naturally opposed processes have not been determined due to challenges in resolving partially unwound nucleosome configurations at atomic resolution. Using a near atomistic protein-DNA model and advanced sampling techniques, we calculate the free energy cost of nucleosome DNA unwinding. Our simulations identify a large energetic barrier that decouples the outer and the inner DNA unwinding into two separate processes, occurring on different time scales. This dynamical decoupling allows the exposure of outer DNA at a modest cost to ensure accessibility while keeping the inner DNA and the histone core intact to maintain stability. We also reveal that this energetic barrier arises from a delayed loss of contacts between disordered histone tails and the DNA and is, surprisingly, largely offset by an entropic contribution from these tails. Implications of this enthalpy entropy compensation for the regulation of nucleosome stability and genome function are discussed.

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.

MutS homolog sliding clamps shield the DNA from binding proteins

Sliding clamps on DNA consist of evolutionarily conserved enzymes that coordinate DNA replication, repair, and the cellular DNA damage response. MutS homolog (MSH) proteins initiate mismatch repair (MMR) by recognizing mispaired nucleotides and in the presence of ATP form stable sliding clamps that randomly diffuse along the DNA. The MSH sliding clamps subsequently load MutL homolog (MLH/PMS) proteins that form a second extremely stable sliding clamp, which together coordinate downstream MMR components with the excision-initiation site that may be hundreds to thousands of nucleotides distant from the mismatch. Specific or nonspecific binding of other proteins to the DNA between the mismatch and the distant excision-initiation site could conceivably obstruct the free diffusion of these MMR sliding clamps, inhibiting their ability to initiate repair. Here, we employed bulk biochemical analysis, single-molecule fluorescence imaging, and mathematical modeling to determine how sliding clamps might overcome such hindrances along the DNA. Using both bacterial and human MSH proteins, we found that increasing the number of MSH sliding clamps on a DNA decreased the association of the Escherichia coli transcriptional repressor LacI to its cognate promoter LacO. Our results suggest a simple mechanism whereby thermal diffusion of MSH sliding clamps along the DNA alters the association kinetics of other DNA-binding proteins over extended distances. These observations appear generally applicable to any stable sliding clamp that forms on 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.

Scoring Targets of Transcription in Bacteria Rather than Focusing on Individual Binding Sites

Reliable identification of targets of bacterial regulators is necessary to understand bacterial gene expression regulation. These targets are commonly predicted by searching for high-scoring binding sites in the upstream genomic regions, which typically leads to a large number of false positives. In contrast to the common approach, here we propose a novel concept, where overrepresentation of the scoring distribution that corresponds to the entire searched region is assessed, as opposed to predicting individual binding sites. We explore two implementations of this concept, based on Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) tests, which both provide straightforward P-value estimates for predicted targets. This approach is implemented for pleiotropic bacterial regulators, including σ70 (bacterial housekeeping σ factor) target predictions, which is a classical bioinformatics problem characterized by low specificity. We show that KS based approach is both faster and more accurate, departing from the current paradigm of AD being slower, but more accurate. Moreover, KS approach leads to a significant increase in the search accuracy compared to the standard approach, while at the same time straightforwardly assigning well established P-values to each potential target. Consequently, the new KS based method proposed here, which assigns P-values to fixed length upstream regions, provides a fast and accurate approach for predicting bacterial transcription targets.

Probing Nucleosome Stability with a DNA Origami Nanocaliper

The organization of eukaryotic DNA into nucleosomes and chromatin undergoes dynamic structural changes to regulate genome processing, including transcription and DNA repair. Critical chromatin rearrangements occur over a wide range of distances, including the mesoscopic length scale of tens of nanometers. However, there is a lack of methodologies that probe changes over this mesoscopic length scale within chromatin. We have designed, constructed, and implemented a DNA-based nanocaliper that probes this mesoscopic length scale. We developed an approach of integrating nucleosomes into our nanocaliper at two attachment points with over 50% efficiency. Here, we focused on attaching the two DNA ends of the nucleosome to the ends of the two nanocaliper arms, so the hinge angle is a readout of the nucleosome end-to-end distance. We demonstrate that nucleosomes integrated with 6, 26, and 51 bp linker DNA are partially unwrapped by the nanocaliper by an amount consistent with previously observed structural transitions. In contrast, the nucleosomes integrated with the longer 75 bp linker DNA remain fully wrapped. We found that the nanocaliper angle is a sensitive measure of nucleosome disassembly and can read out transcription factor (TF) binding to its target site within the nucleosome. Interestingly, the nanocaliper not only detects TF binding but also significantly increases the probability of TF occupancy at its site by partially unwrapping the nucleosome. These studies demonstrate the feasibility of using DNA nanotechnology to both detect and manipulate nucleosome structure, which provides a foundation of future mesoscale studies of nucleosome and chromatin structural dynamics.

Dynamic DNA binding licenses a repair factor to bypass roadblocks in search of DNA lesions

DNA-binding proteins search for specific targets via facilitated diffusion along a crowded genome. However, little is known about how crowded DNA modulates facilitated diffusion and target recognition. Here we use DNA curtains and single-molecule fluorescence imaging to investigate how Msh2-Msh3, a eukaryotic mismatch repair complex, navigates on crowded DNA. Msh2-Msh3 hops over nucleosomes and other protein roadblocks, but maintains sufficient contact with DNA to recognize a single lesion. In contrast, Msh2-Msh6 slides without hopping and is largely blocked by protein roadblocks. Remarkably, the Msh3-specific mispair-binding domain (MBD) licences a chimeric Msh2-Msh6(3MBD) to bypass nucleosomes. Our studies contrast how Msh2-Msh3 and Msh2-Msh6 navigate on a crowded genome and suggest how Msh2-Msh3 locates DNA lesions outside of replication-coupled repair. These results also provide insights into how DNA repair factors search for DNA lesions in the context of chromatin.

Mismatch repair

Highly conserved MutS homologs (MSH) and MutL homologs (MLH/PMS) are the fundamental components of mismatch repair (MMR). After decades of debate, it appears clear that the MSH proteins initiate MMR by recognizing a mismatch and forming multiple extremely stable ATP-bound sliding clamps that diffuse without hydrolysis along the adjacent DNA. The function(s) of MLH/PMS proteins is less clear, although they too bind ATP and are targeted to MMR by MSH sliding clamps. Structural analysis combined with recent real-time single molecule and cellular imaging technologies are providing new and detailed insight into the thermal-driven motions that animate the complete MMR mechanism.

Histone core phosphorylation regulates DNA accessibility

Nucleosome unwrapping dynamics provide transient access to the complexes involved in DNA transcription, repair, and replication, whereas regulation of nucleosome unwrapping modulates occupancy of these complexes. Histone H3 is phosphorylated at tyrosine 41 (H3Y41ph) and threonine 45 (H3T45ph). H3Y41ph is implicated in regulating transcription, whereas H3T45ph is involved in DNA replication and apoptosis. These modifications are located in the DNA-histone interface near where the DNA exits the nucleosome, and are thus poised to disrupt DNA-histone interactions. However, the impact of histone phosphorylation on nucleosome unwrapping and accessibility is unknown. We find that the phosphorylation mimics H3Y41E and H3T45E, and the chemically correct modification, H3Y41ph, significantly increase nucleosome unwrapping. This enhances DNA accessibility to protein binding by 3-fold. H3K56 acetylation (H3K56ac) is also located in the same DNA-histone interface and increases DNA unwrapping. H3K56ac is implicated in transcription regulation, suggesting that H3Y41ph and H3K56ac could function together. We find that the combination of H3Y41ph with H3K56ac increases DNA accessibility by over an order of magnitude. These results suggest that phosphorylation within the nucleosome DNA entry-exit region increases access to DNA binding complexes and that the combination of phosphorylation with acetylation has the potential to significantly influence DNA accessibility to transcription regulatory complexes.

Functional roles of nucleosome stability and dynamics

Nucleosome is a histone-DNA complex known as the fundamental repeating unit of chromatin. Up to 90% of eukaryotic DNA is wrapped around consecutive octamers made of the core histones H2A, H2B, H3 and H4. Nucleosome positioning affects numerous cellular processes that require robust and timely access to genomic DNA, which is packaged into the tight confines of the cell nucleus. In living cells, nucleosome positions are determined by intrinsic histone-DNA sequence preferences, competition between histones and other DNA-binding proteins for genomic sequence, and ATP-dependent chromatin remodelers. We discuss the major energetic contributions to nucleosome formation and remodeling, focusing especially on partial DNA unwrapping off the histone octamer surface. DNA unwrapping enables efficient access to nucleosome-buried binding sites and mediates rapid nucleosome removal through concerted action of two or more DNA-binding factors. High-resolution, genome-scale maps of distances between neighboring nucleosomes have shown that DNA unwrapping and nucleosome crowding (mutual invasion of nucleosome territories) are much more common than previously thought. Ultimately, constraints imposed by nucleosome energetics on the rates of ATP-dependent and spontaneous chromatin remodeling determine nucleosome occupancy genome-wide, and shape pathways of cellular response to environmental stresses.

Quantitative models for accelerated protein dissociation from nucleosomal DNA

Binding of transcription factors to their binding sites in promoter regions is the fundamental event in transcriptional gene regulation. When a transcription factor binding site is located within a nucleosome, the DNA has to partially unwrap from the nucleosome to allow transcription factor binding. This reduces the rate of transcription factor binding and is a known mechanism for regulation of gene expression via chromatin structure. Recently a second mechanism has been reported where transcription factor off-rates are dramatically increased when binding to target sites within the nucleosome. There are two possible explanations for such an increase in off-rate short of an active role of the nucleosome in pushing the transcription factor off the DNA: (i) for dimeric transcription factors the nucleosome can change the equilibrium between monomeric and dimeric binding or (ii) the nucleosome can change the equilibrium between specific and non-specific binding to the DNA. We explicitly model both scenarios and find that dimeric binding can explain a large increase in off-rate while the non-specific binding model cannot be reconciled with the large, experimentally observed increase. Our results suggest a general mechanism how nucleosomes increase transcription factor dissociation to promote exchange of transcription factors and regulate gene expression.

Single-molecule views of MutS on mismatched DNA

Base-pair mismatches that occur during DNA replication or recombination can reduce genetic stability or conversely increase genetic diversity. The genetics and biophysical mechanism of mismatch repair (MMR) has been extensively studied since its discovery nearly 50 years ago. MMR is a strand-specific excision-resynthesis reaction that is initiated by MutS homolog (MSH) binding to the mismatched nucleotides. The MSH mismatch-binding signal is then transmitted to the immediate downstream MutL homolog (MLH/PMS) MMR components and ultimately to a distant strand scission site where excision begins. The mechanism of signal transmission has been controversial for decades. We have utilized single molecule Forster Resonance Energy Transfer (smFRET), Fluorescence Tracking (smFT) and Polarization Total Internal Reflection Fluorescence (smP-TIRF) to examine the interactions and dynamic behaviors of single Thermus aquaticus MutS (TaqMutS) particles on mismatched DNA. We determined that TaqMutS forms an incipient clamp to search for a mismatch in ~1 s intervals by 1-dimensional (1D) thermal fluctuation-driven rotational diffusion while in continuous contact with the helical duplex DNA. When MutS encounters a mismatch it lingers for ~3 s to exchange bound ADP for ATP (ADP→ATP exchange). ATP binding by TaqMutS induces an extremely stable clamp conformation (~10 min) that slides off the mismatch and moves along the adjacent duplex DNA driven simply by 1D thermal diffusion. The ATP-bound sliding clamps rotate freely while in discontinuous contact with the DNA. The visualization of a train of MSH proteins suggests that dissociation of ATP-bound sliding clamps from the mismatch permits multiple mismatch-dependent loading events. These direct observations have provided critical clues into understanding the molecular mechanism of MSH proteins during MMR.

Mechanistic heterogeneity in site recognition by the structurally homologous DNA-binding domains of the ETS family transcription factors Ets-1 and PU.1

ETS family transcription factors regulate diverse genes through binding at cognate DNA sites that overlap substantially in sequence. The DNA-binding domains of ETS proteins (ETS domains) are highly conserved structurally yet share limited amino acid homology. To define the mechanistic implications of sequence diversity within the ETS family, we characterized the thermodynamics and kinetics of DNA site recognition by the ETS domains of Ets-1 and PU.1, which represent the extremes in amino acid divergence among ETS proteins. Even though the two ETS domains bind their optimal sites with similar affinities under physiologic conditions, their nature of site recognition differs strikingly in terms of the role of hydration and counter ion release. The data suggest two distinct mechanisms wherein Ets-1 follows a "dry" mechanism that rapidly parses sites through electrostatic interactions and direct protein-DNA contacts, whereas PU.1 utilizes hydration to interrogate sequence-specific sites and form a long-lived complex relative to the Ets-1 counterpart. The kinetic persistence of the high affinity PU.1 · DNA complex may be relevant to an emerging role of PU.1, but not Ets-1, as a pioneer transcription factor in vivo. In addition, PU.1 activity is critical to the development and function of macrophages and lymphocytes, which present osmotically variable environments, and hydration-dependent specificity may represent an important regulatory mechanism in vivo, a hypothesis that finds support in gene expression profiles of primary murine macrophages.

Ubiquitous nucleosome crowding in the yeast genome

Nucleosomes may undergo a conformational change in which a stretch of DNA peels off the histone octamer surface as a result of thermal fluctuations or interactions with chromatin remodelers. Thus, neighboring nucleosomes may invade each other's territories by DNA unwrapping and translocation, or through initial assembly in partially wrapped states. A recent high-resolution map of distances between dyads of neighboring nucleosomes in Saccharomyces cerevisiae reveals that nucleosomes frequently overlap DNA territories of their neighbors. This conclusion is supported by lower-resolution maps of S. cerevisiae nucleosome lengths based on micrococcal nuclease digestion and paired-end sequencing. The average length of wrapped DNA follows a stereotypical pattern in genes and promoters, correlated with the well-known distribution of nucleosome occupancy: nucleosomal DNA tends to be shorter in promoters and longer in coding regions. To explain these observations, we have developed a biophysical model that uses a 10-11-bp periodic histone-DNA binding energy profile. The profile is based on the pattern of histone-DNA contacts in nucleosome crystal structures, as well as the idea of linker length discretization caused by higher-order chromatin structure. Our model is in agreement with the observed genome-wide distributions of interdyad distances, wrapped DNA lengths, and nucleosome occupancies. Furthermore, our approach explains in vitro measurements of the accessibility of nucleosome-covered target sites and nucleosome-induced cooperativity between DNA-binding factors. We rule out several alternative scenarios of histone-DNA interactions as inconsistent with the genomic data.

Taking into account nucleosomes for predicting gene expression

The eukaryotic genome is organized in a chain of nucleosomes that consist of 145-147 bp of DNA wrapped around a histone octamer protein core. Binding of transcription factors (TF) to nucleosomal DNA is frequently impeded, which makes it a challenging task to calculate TF occupancy at a given regulatory genomic site for predicting gene expression. Here, we review methods to calculate TF binding to DNA in the presence of nucleosomes. The main theoretical problems are (i) the computation speed that is becoming a bottleneck when partial unwrapping of DNA from the nucleosome is considered, (ii) the perturbation of the binding equilibrium by the activity of ATP-dependent chromatin remodelers, which translocate nucleosomes along the DNA, and (iii) the model parameterization from high-throughput sequencing data and fluorescence microscopy experiments in living cells. We discuss strategies that address these issues to efficiently compute transcription factor binding in chromatin.

Dynamics of forced nucleosome unraveling and role of nonuniform histone-DNA interactions

A coarse-grained model of the nucleosome is introduced to investigate the dynamics of force-induced unwrapping of DNA from histone octamers. In this model, the DNA is treated as a charged, discrete worm-like chain, and the octamer is treated as a rigid cylinder carrying a positively charged superhelical groove that accommodates 1.7 turns of DNA. The groove charges are parameterized to reproduce the nonuniform histone/DNA interaction free energy profile and the loading rate-dependent unwrapping forces, both obtained from single-molecule experiments. Brownian dynamics simulations of the model under constant loading conditions reveal that nucleosome unraveling occurs in three distinct stages. At small extensions, the flanking DNA exhibits rapid unwrapping-rewrapping (breathing) dynamics and the octamer flips ∼180° and moves toward the pulling axis. At intermediate extensions, the outer turn of DNA unwraps gradually and the octamer swivels about the taut linkers and flips a further ∼90° to orient its superhelical axis almost parallel to the pulling axis. At large extensions, a portion of the inner turn unwraps abruptly with a notable rip in the force-extension plot and a >90° flip of the octamer. The remaining inner turn unwraps reversibly to leave a small portion of DNA attached to the octamer despite extended pulling. Our simulations further reveal that the nonuniform histone/DNA interactions in canonical nucleosomes serve to: stabilize the inner turn against unraveling while enhancing the breathing dynamics of the nucleosome and prevent dissociation of the octamer from the DNA while facilitating its mobility along the DNA. Thus, the modulation of the histone/DNA interactions could constitute one possible mechanism for regulating the accessibility of the nucleosome-wound DNA sequences.

Unwinding and rewinding the nucleosome inner turn: force dependence of the kinetic rate constants

A simple model for the force-dependent unwinding and rewinding rates of the nucleosome inner turn is constructed and quantitatively compared to the results of recent measurements [A. H. Mack et al., J. Mol. Biol. 423, 687 (2012)]. First, a coarse-grained model for the histone-DNA free-energy landscape that incorporates both an elastic free-energy barrier and specific histone-DNA bonds is developed. Next, a theoretical expression for the rate of transitions across a piecewise linear free-energy landscape with multiple minima and maxima is presented. Then, the model free-energy landscape, approximated as a piecewise linear function, and the theoretical expression for the transition rates are combined to construct a model for the force-dependent unwinding and rewinding rates of the nucleosome inner turn. Least-mean-squares fitting of the model rates to the rates observed in recent experiments rates demonstrates that this model is able to well describe the force-dependent unwinding and rewinding rates of the nucleosome inner turn, observed in the recent experiments, except at the highest forces studied, where an additional ad hoc term is required to describe the data, which may be interpreted as an indication of an alternate high-force nucleosome disassembly pathway, that bypasses simple unwinding. The good agreement between the measurements and the model at lower forces demonstrates that both specific histone-DNA contacts and an elastic free-energy barrier play essential roles for nucleosome winding and unwinding, and quantifies their relative contributions.

There and back again: new single-molecule insights in the motion of DNA repair proteins

Cellular DNA repair machines are constantly at work supporting the integrity of our genomes. Numerous proteins cooperate to form a complex and adaptive system dedicated to detection and timely processing of DNA damage. The molecular underpinnings of how these proteins locate and discriminate DNA lesions, match homologous sequences, mend the DNA and attend to a replication in distress are of a paramount biomedical importance, but in many cases remain unclear. Combined with more conventional tools, single-molecule biochemistry has been stepping in to address the age-old problems in the DNA repair field. This review will address new insights into diffusive properties of three DNA repair systems: I will discuss the emerging model of how MutS homologues locate and respond to mismatches in the dsDNA; the mechanism by which RAD52 promotes annealing of complementary DNA strands coated with ssDNA binding protein RPA; and how the nucleoprotein filament formed by RecA recombinase on ssDNA searches for homology within duplex DNA. These three distinct DNA repair factors exemplify the dynamic nature of cellular DNA repair machines revealed by single-molecule studies.

Kinetic control of nucleosome displacement by ISWI/ACF chromatin remodelers

Chromatin structure is dynamically organized by chromatin remodelers, motor protein complexes which move and remove nucleosomes. The regulation of remodeler action has recently been proposed to underlie a kinetic proofreading scheme which combines the recognition of histone-tail states and the ATP-dependent loosening of DNA around nucleosomes. Members of the ISWI-family of remodelers additionally recognize linker length between nucleosomes. Here, we show that the additional proofreading step involving linker length alone is sufficient to promote the formation of regular arrays of nucleosomes. ATP-dependent remodeling by bidirectional motors is shown to reinforce positioning as compared to statistical positioning.

Regulation of the nucleosome unwrapping rate controls DNA accessibility

Eukaryotic genomes are repetitively wrapped into nucleosomes that then regulate access of transcription and DNA repair complexes to DNA. The mechanisms that regulate extrinsic protein interactions within nucleosomes are unresolved. We demonstrate that modulation of the nucleosome unwrapping rate regulates protein binding within nucleosomes. Histone H3 acetyl-lysine 56 [H3(K56ac)] and DNA sequence within the nucleosome entry-exit region additively influence nucleosomal DNA accessibility by increasing the unwrapping rate without impacting rewrapping. These combined epigenetic and genetic factors influence transcription factor (TF) occupancy within the nucleosome by at least one order of magnitude and enhance nucleosome disassembly by the DNA mismatch repair complex, hMSH2-hMSH6. Our results combined with the observation that ∼30% of Saccharomyces cerevisiae TF-binding sites reside in the nucleosome entry-exit region suggest that modulation of nucleosome unwrapping is a mechanism for regulating transcription and DNA repair.

Structural dynamics of nucleosomes at single-molecule resolution

The detailed mechanisms of how DNA that is assembled around a histone core can be accessed by DNA-binding proteins for transcription, replication, or repair, remain elusive nearly 40 years after Kornberg's nucleosome model was proposed. Uncovering the structural dynamics of nucleosomes is a crucial step in elucidating the mechanisms regulating genome accessibility. This requires the deconvolution of multiple structural states within an ensemble. Recent advances in single-molecule methods enable unprecedented efficiency in examining subpopulation dynamics. In this review, we summarize studies of nucleosome structure and dynamics from single-molecule approaches and how they advance our understanding of the mechanisms that govern DNA transactions.

Unzipping single DNA molecules to study nucleosome structure and dynamics

DNA unzipping is a powerful tool to study protein-DNA interactions at the single-molecule level. In this chapter, we provide a detailed and practical guide to performing this technique with an optical trap, using nucleosome studies as an example. We detail protocols for preparing an unzipping template, constructing and calibrating the instrument, and acquiring, processing, and analyzing unzipping data. We also summarize major results from utilization of this technique for the studies of nucleosome structure, dynamics, positioning, and remodeling.

Histone fold modifications control nucleosome unwrapping and disassembly

Nucleosomes are stable DNA-histone protein complexes that must be unwrapped and disassembled for genome expression, replication, and repair. Histone posttranslational modifications (PTMs) are major regulatory factors of these nucleosome structural changes, but the molecular mechanisms associated with PTM function remains poorly understood. Here we demonstrate that histone PTMs within distinct structured regions of the nucleosome directly regulate the inherent dynamic properties of the nucleosome. Precise PTMs were introduced into nucleosomes by chemical ligation. Single molecule magnetic tweezers measurements determined that only PTMs near the nucleosome dyad increase the rate of histone release in unwrapped nucleosomes. In contrast, FRET and restriction enzyme analysis reveal that only PTMs throughout the DNA entry-exit region increase unwrapping and enhance transcription factor binding to nucleosomal DNA. These results demonstrate that PTMs in separate structural regions of the nucleosome control distinct dynamic events, where the dyad regulates disassembly while the DNA entry-exit region regulates unwrapping. These studies are consistent with the conclusion that histone PTMs may independently influence nucleosome dynamics and associated chromatin functions.