Internucleosomal interactions mediated by histone tails allow distant communication in chromatin

Beyond Sequence: Internucleosomal Interactions Dominate Array Assembly

The organization of the nucleosome array is a critical component of the chromatin assembly into higher order structure as well as its function. Here, we investigated the contributions of the DNA sequence and internucleosomal interactions on the organization of the nucleosomal arrays in compact structures using atomic force microscopy. We assembled nucleosomes on DNA substrates allowing for the formation of tetranucleosomes. We found that nucleosomes are capable of close positioning with no discernible space between them, even in the case of assembled dinucleosomes. This morphology of the array is in contrast with that observed for arrays assembled with repeats of the nucleosome positioning motifs separated by uniform spacers. Simulated assembly of tetranucleosomes by random placement along the substrates revealed that nucleosome array compaction is promoted by the interaction of the nucleosomes. We developed a theoretical model to account for the role of DNA sequence and internucleosomal interactions in the formation of the nucleosome structures. These findings suggest that, in the chromatin assembly, the affinity of the nucleosomes to the DNA sequence and the strengths of the internucleosomal interactions are the two major factors defining the compactness of the chromatin.

The Effects of Histone H2B Ubiquitylations on the Nucleosome Structure and Internucleosomal Interactions

Eukaryotic gene compaction takes place at multiple levels to package DNA to chromatin and chromosomes. Two of the most fundamental levels of DNA packaging are at the nucleosome and dinucleosome stacks. The nucleosome is the basic gene-packing unit and is composed of DNA wrapped around a histone core. Nucleosomes stack with one another for further compaction of DNA. The first stacking step leads to dinucleosome formation, which is driven by internucleosomal interactions between various parts of two nucleosomes. Histone proteins are rich targets for post-translational modifications, some of which affect the structure of the nucleosome and the interactions between nucleosomes. These effects are often implicated in the regulation of various genomic transactions. In particular, histone H2B ubiquitylation has been associated with facilitated transcription and hexasome formation. Here, we employed semi-synthetically ubiquitylated histone H2B and single-molecule FRET to investigate the effects of H2B ubiquitylations at lysine 34 (H2BK34) and lysine 120 (H2BK120) on the structure of the nucleosome and the interactions between two nucleosomes. Our results suggest that H2BK34 ubiquitylation widens the DNA gyre gap in the nucleosome and stabilizes long- and short-range internucleosomal interactions while H2BK120 ubiquitylation does not affect the nucleosome structure or internucleosomal interactions. These results suggest potential roles for H2B ubiquitylations in facilitated transcription and hexasome formation while maintaining the structural integrity of chromatin.

Non-cooperative mechanism for bounded and ultrasensitive chromatin remodeling

Chromatin remodeling is an essential form of gene regulation that is involved in a variety of biological processes. We develop a theoretical model that takes advantage of percolation effects at the level of nucleosome interactions, which allows for ultrasensitive chromatin expansion. This model is non-cooperative and readily provides spatial bounds to the expansion region, preventing uncontrolled remodeling events. We explore different chromatin architectures and the ultrasensitivity of the chromatin density as a function of transcription factor concentration. We also compare our model with experimental data involving an inhibitor of nucleosome acetylation. These results suggest a novel mechanism for spatially-bounded chromatin remodeling and they provide means for quantitative comparisons between proposed models of chromatin architecture.

Histone H3 and H4 tails play an important role in nucleosome phase separation

Chromatin organization and its dynamic regulation are crucial in governing the temporal and spatial accessibility of DNA for proper gene expression. Disordered chains of nucleosomes comprise the basis of eukaryotic chromatin, forming higher-level organization across a range of length scales. Models of chromatin organization involving phase separation driven by chromatin-associating proteins have been proposed. More recently, evidence has emerged that nucleosome arrays can phase separate in the absence of other protein factors, yet questions remain regarding the molecular basis of chromatin phase separation that governs this dynamic nuclear organization. Here, we break chromatin down into its most basic subunit, the nucleosome core particle, and investigate phase separation using turbidity assays in conjunction with differential interference contrast microscopy. We show that, at physiologically-relevant concentrations, this fundamental subunit of chromatin undergoes phase separation. Individually removing the H3 and H4 tails abrogates phase separation under the same conditions. Taking a reductionist approach to investigate H3 and H4 tail peptide interactions in-trans with DNA and nucleosome core particles supports the direct involvement of these tails in chromatin phase separation. These results provide insight into fundamental mechanisms underlying phase separation of chromatin, which starts at the level of the nucleosome core particle, and support that long-range inter-nucleosomal interactions are sufficient to drive phase separation at nuclear concentrations. Additionally, our data have implications for understanding crosstalk between histone tails and provide a lens through which to interpret the effect of histone post-translational modifications and sequence variants. STATEMENT OF SIGNIFICANCE: Emerging models propose that chromatin organization is based in phase separation, however, mechanisms that drive this dynamic nuclear organization are only beginning to be understood. Previous focus has been on phase separation driven by chromatin-associating proteins, but this has recently shifted to recognize a direct role of chromatin in phase separation. Here, we take a fundamental approach in understanding chromatin phase separation and present new findings that the basic subunit of chromatin, the nucleosome core particle, undergoes phase separation under physiological concentrations of nucleosome and monovalent salt. Furthermore, the histone H3 and H4 tails are involved in phase separation in a manner independent of histone-associating proteins. These data suggest that H3 and H4 tail epigenetic factors may modulate chromatin phase separation.

Effect of histone H4 tail on nucleosome stability and internucleosomal interactions

Chromatin structure is dictated by nucleosome assembly and internucleosomal interactions. The tight wrapping of nucleosomes inhibits gene expression, but modifications to histone tails modulate chromatin structure, allowing for proper genetic function. The histone H4 tail is thought to play a large role in regulating chromatin structure. Here we investigated the structure of nucleosomes assembled with a tail-truncated H4 histone using Atomic Force Microscopy. We assembled tail-truncated H4 nucleosomes on DNA templates allowing for the assembly of mononucleosomes or dinucleosomes. Mononucleosomes assembled on nonspecific DNA led to decreased DNA wrapping efficiency. This effect is less pronounced for nucleosomes assembled on positioning motifs. Dinucleosome studies resulted in the discovery of two effects- truncation of the H4 tail does not diminish the preferential positioning observed in full-length nucleosomes, and internucleosomal interaction eliminates the DNA unwrapping effect. These findings provide insight on the role of histone H4 in chromatin structure and stability.

Molecular dynamics simulations reveal how H3K56 acetylation impacts nucleosome structure to promote DNA exposure for lesion sensing

The first order of DNA packaging is the nucleosome with the DNA wrapped around the histone octamer. This leaves the nucleosomal DNA with access restrictions, which impose a significant barrier to repair of damaged DNA. The efficiency of DNA repair has been related to nucleosome structure and chromatin status, which is modulated in part by post-translational modifications (PTMs) of histones. Numerous studies have suggested a role for acetylation of lysine at position 56 of the H3 histone (H3K56ac) in various DNA transactions, including the response to DNA damage and its association with human cancer. Biophysical studies have revealed that H3K56ac increases DNA accessibility by facilitating spontaneous and transient unwrapping motions of the DNA ends. However, how this acetylation mark modulates nucleosome structure and dynamics to promote accessibility to the damaged DNA for repair factors and other proteins is still poorly understood. Here, we utilize approximately 5-6 microseconds of atomistic molecular dynamics simulations to delineate the impact of H3K56 acetylation on the nucleosome structure and dynamics, and to elucidate how these nucleosome properties are further impacted when a bulky benzo[a]pyrene-derived DNA lesion is placed near the acetylation site. Our findings reveal that H3K56ac alone induces considerable disturbance to the histone-DNA/histone-histone interactions, and amplifies the distortions imposed by the presence of the lesion. Our work highlights the important role of H3K56 acetylation in response to DNA damage and depicts how access to DNA lesions by the repair machinery can be facilitated within the nucleosome via a key acetylation event.

Surprising Twists in Nucleosomal DNA with Implication for Higher-order Folding

While nucleosomes are dynamic entities that must undergo structural deformations to perform their functions, the general view from available high-resolution structures is a largely static one. Even though numerous examples of twist defects have been documented, the DNA wrapped around the histone core is generally thought to be overtwisted. Analysis of available high-resolution structures from the Protein Data Bank reveals a heterogeneous distribution of twist along the nucleosomal DNA, with clear patterns that are consistent with the literature, and a significant fraction of structures that are undertwisted. The subtle differences in nucleosomal DNA folding, which extend beyond twist, have implications for nucleosome disassembly and modeled higher-order structures. Simulations of oligonucleosome arrays built with undertwisted models behave very differently from those constructed from overtwisted models, in terms of compaction and inter-nucleosome contacts, introducing configurational changes equivalent to those associated with 2-3 base-pair changes in nucleosome spacing. Differences in the nucleosomal DNA pathway, which underlie the way that DNA enters and exits the nucleosome, give rise to different nucleosome-decorated minicircles and affect the topological mix of configurational states.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface

The three-dimensional organization of chromatin is expected to play critical roles in regulating genome functions. High-resolution characterization of its structure and dynamics could improve our understanding of gene regulation mechanisms but has remained challenging. Using a near-atomistic model that preserves the chemical specificity of protein-DNA interactions at residue and base-pair resolution, we studied the stability and folding pathways of a tetra-nucleosome. Dynamical simulations performed with an advanced sampling technique uncovered multiple pathways that connect open chromatin configurations with the zigzag crystal structure. Intermediate states along the simulated folding pathways resemble chromatin configurations reported from in situ experiments. We further determined a six-dimensional free energy surface as a function of the inter-nucleosome distances via a deep learning approach. The zigzag structure can indeed be seen as the global minimum of the surface. However, it is not favored by a significant amount relative to the partially unfolded, in situ configurations. Chemical perturbations such as histone H4 tail acetylation and thermal fluctuations can further tilt the energetic balance to stabilize intermediate states. Our study provides insight into the connection between various reported chromatin configurations and has implications on the in situ relevance of the 30 nm fiber.

The three-dimensional organization of chromatin plays critical roles in regulating genome function. Here the authors apply a near atomistic model to study the structure and dynamics of the chromatin folding unit - the tetra-nucleosome - to provide insight into how chromatin folds.

Topological polymorphism of nucleosome fibers and folding of chromatin

We discuss recent observations of polymorphic chromatin packaging at the oligonucleosomal level and compare them with computer simulations. Our computations reveal two topologically different families of two-start 30-nm fiber conformations distinguished by the linker length L; fibers with L ≈ 10n and L ≈ 10n+5 basepairs have DNA linking numbers per nucleosome of ΔLk ≈ -1.5 and -1.0, respectively (where n is a natural number). Although fibers with ΔLk ≈ -1.5 were observed earlier, the topoisomer with ΔLk ≈ -1.0 is novel. These predictions were confirmed experimentally for circular nucleosome arrays with precisely positioned nucleosomes. We suggest that topological polymorphism of chromatin may play a role in transcription, with the {10n+5} fibers producing transcriptionally competent chromatin structures. This hypothesis is consistent with available data for yeast and, partially, for fly. We show that both fiber topoisomers (with ΔLk ≈ -1.5 and -1.0) have to be taken into account to interpret experimental data obtained using new techniques: genome-wide Micro-C, Hi-CO, and RICC-seq, as well as self-association of nucleosome arrays in vitro. The relative stability of these topoisomers is likely to depend on epigenetic histone modifications modulating the strength of internucleosome interactions. Potentially, our findings may reflect a general tendency of functionally distinct parts of the genome to retain topologically different higher-order structures.

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.

Molecular Basis of the Function of Transcriptional Enhancers

Transcriptional enhancers are major genomic elements that control gene activity in eukaryotes. Recent studies provided deeper insight into the temporal and spatial organization of transcription in the nucleus, the role of non-coding RNAs in the process, and the epigenetic control of gene expression. Thus, multiple molecular details of enhancer functioning were revealed. Here, we describe the recent data and models of molecular organization of enhancer-driven transcription.

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

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

Biophysical Reviews' "Meet the Editors Series"-a profile of Wilma K. Olson

As one of the five Executive Editors of Biophysical Reviews I have been asked to provide this short biographical sketch for the readers of the journal. I have been a member of the Editorial Board since the inception of the journal in 2008 and an Executive Editor since 2014. I hold a B.S. degree in Chemistry from the University of Delaware and a Ph.D. in Physical Chemistry from Stanford University. Except for a year as a Damon Runyon Postdoctoral Fellow at Columbia University, I have spent my entire professional career at Rutgers, the State University of New Jersey, where I am currently the Mary I. Bunting Professor of Chemistry and Chemical Biology. I served for many years as Founding Director of the Rutgers University Center for Molecular Biophysics and Biophysical Chemistry, and have trained undergraduate, graduate, and postdoctoral students from a variety of academic disciplines.

[Opposite Effects of Histone H1 and HMGN5 Protein on Distant Interactions in Chromatin]

Transcriptional enhancers in the cell nuclei typically interact with the target promoters in cis over long stretches of chromatin, but the mechanism of this communication remains unknown. Previously we have developed a defined in vitro system for quantitative analysis of the rate of distant enhancer-promoter communication (EPC) and have shown that the chromatin fibers maintain efficient distant EPC in cis. Here we investigate the roles of linker histone H1 and HMGN5 protein in EPC. A considerable negative effect of histone H1 on EPC depending on its C- and N-tails was shown. Protein HMGN5 that affects chromatin compaction and is associated with active chromatin counteracts EPC inhibition by H1. The data suggest that the efficiency of the interaction between the enhancer and the promoter depends on the structure and dynamics of the chromatin fiber localized between them and can be regulated by proteins associated with chromatin.

Dynamic chromatin organization in the cell

The organization and regulation of genomic DNA as nuclear chromatin is necessary for proper DNA function inside living eukaryotic cells. While this has been extensively explored, no true consensus is currently reached regarding the exact mechanism of chromatin organization. The traditional view has assumed that the DNA is packaged into a hierarchy of structures inside the nucleus based on the regular 30-nm chromatin fiber. This is currently being challenged by the fluid-like model of the chromatin which views the chromatin as a dynamic structure based on the irregular 10-nm fiber. In this review, we focus on the recent progress in chromatin structure elucidation highlighting the paradigm shift in chromatin folding mechanism from the classical textbook perspective of the regularly folded chromatin to the more dynamic fluid-like perspective.

Nucleosome-free DNA regions differentially affect distant communication in chromatin

Communication between distantly spaced genomic regions is one of the key features of gene regulation in eukaryotes. Chromatin per se can stimulate efficient enhancer-promoter communication (EPC); however, the role of chromatin structure and dynamics in this process remains poorly understood. Here we show that nucleosome spacing and the presence of nucleosome-free DNA regions can modulate chromatin structure/dynamics and, in turn, affect the rate of EPC in vitro and in silico. Increasing the length of internucleosomal linker DNA from 25 to 60 bp results in more efficient EPC. The presence of longer nucleosome-free DNA regions can positively or negatively affect the rate of EPC, depending upon the length and location of the DNA region within the chromatin fiber. Thus the presence of histone-free DNA regions can differentially affect the efficiency of EPC, suggesting that gene regulation over a distance could be modulated by changes in the length of internucleosomal DNA spacers.

Towards quantitative analysis of gene regulation by enhancers

Enhancers are regulatory DNA sequences that can activate transcription over large distances. Recent studies have revealed the widespread role of distant activation in eukaryotic gene regulation and in the development of various human diseases, including cancer. Here we review recent progress in the field, focusing on new experimental and computational approaches that quantify the role of chromatin structure and dynamics during enhancer-promoter interactions in vitro and in vivo.

Contributions of Sequence to the Higher-Order Structures of DNA

One of the critical unanswered questions in genome biophysics is how the primary sequence of DNA bases influences the global properties of very-long-chain molecules. The local sequence-dependent features of DNA found in high-resolution structures introduce irregularities in the disposition of adjacent residues that facilitate the specific binding of proteins and modulate the global folding and interactions of double helices with hundreds of basepairs. These features also determine the positions of nucleosomes on DNA and the lengths of the interspersed DNA linkers. Like the patterns of basepair association within DNA, the arrangements of nucleosomes in chromatin modulate the properties of longer polymers. The intrachromosomal loops detected in genomic studies contain hundreds of nucleosomes, and given that the simulated configurations of chromatin depend on the lengths of linker DNA, the formation of these loops may reflect sequence-dependent information encoded within the positioning of the nucleosomes. With knowledge of the positions of nucleosomes on a given genome, methods are now at hand to estimate the looping propensities of chromatin in terms of the spacing of nucleosomes and to make a direct connection between the DNA base sequence and larger-scale chromatin folding.

Topological polymorphism of the two-start chromatin fiber

Specific details concerning the spatial organization of nucleosomes in 30 nm fibers remain unknown. To investigate this, we analyzed all stereochemically possible configurations of two-start nucleosome fibers with short DNA linkers L = 13-37 bp (nucleosome repeat length (NRL) = 160-184 bp). Four superhelical parameters-inclination of nucleosomes, twist, rise, and diameter-uniquely describe a regular symmetric fiber. The energy of a fiber is defined as the sum of four terms: elastic energy of the linker DNA, steric repulsion, electrostatics, and a phenomenological (H4 tail-acidic patch) interaction between two stacked nucleosomes. By optimizing the fiber energy with respect to the superhelical parameters, we found two types of topological transition in fibers (associated with the change in inclination angle): one caused by an abrupt 360° change in the linker DNA twisting (change in the DNA linking number, ΔLk = 1), and another caused by overcrossing of the linkers (ΔLk = 2). To the best of our knowledge, this topological polymorphism of the two-start fibers was not reported in the computations published earlier. Importantly, the optimal configurations of the fibers with linkers L = 10n and 10n + 5 bp are characterized by different values of the DNA linking number-that is, they are topologically different. Our results are consistent with experimental observations, such as the inclination 60° to 70° (the angle between the nucleosomal disks and the fiber axis), helical rise, diameter, and left-handedness of the fibers. In addition, we make several testable predictions, among them different degrees of DNA supercoiling in fibers with L = 10n and 10n + 5 bp, different flexibility of the two types of fibers, and a correlation between the local NRL and the level of transcription in different parts of the yeast genome.

On the role of inter-nucleosomal interactions and intrinsic nucleosome dynamics in chromatin function

Evidence is emerging that many diseases result from defects in gene functions, which, in turn, depend on the local chromatin environment of a gene. However, it still remains not fully clear how chromatin activity code is 'translated' to the particular 'activating' or 'repressing' chromatin structural transition. Commonly, chromatin remodeling in vitro was studied using mononucleosomes as a model. However, recent data suggest that structural reorganization of a single mononucleosome is not equal to remodeling of a nucleosome particle under multinucleosomal content - such as, interaction of nucleosomes via flexible histone termini could significantly alter the mode (and the resulting products) of nucleosome structural transitions. It is becoming evident that a nucleosome array does not constitute just a 'polymer' of individual 'canonical' nucleosomes due to multiple inter-nucleosomal interactions which affect nucleosome dynamics and structure. It could be hypothesized, that inter-nucleosomal interactions could act in cooperation with nucleosome inherent dynamics to orchestrate DNA-based processes and promote formation and stabilization of highly-dynamic, accessible structure of a nucleosome array. In the proposed paper we would like to discuss the nucleosome dynamics within the chromatin fiber mainly as it pertains to the roles of the structural changes mediated by inter-nucleosomal interactions.

Hierarchical looping of zigzag nucleosome chains in metaphase chromosomes

The architecture of higher-order chromatin in eukaryotic cell nuclei is largely unknown. Here, we use electron microscopy-assisted nucleosome interaction capture (EMANIC) cross-linking experiments in combination with mesoscale chromatin modeling of 96-nucleosome arrays to investigate the internal organization of condensed chromatin in interphase cell nuclei and metaphase chromosomes at nucleosomal resolution. The combined data suggest a novel hierarchical looping model for chromatin higher-order folding, similar to rope flaking used in mountain climbing and rappelling. Not only does such packing help to avoid tangling and self-crossing, it also facilitates rope unraveling. Hierarchical looping is characterized by an increased frequency of higher-order internucleosome contacts for metaphase chromosomes compared with chromatin fibers in vitro and interphase chromatin, with preservation of a dominant two-start zigzag organization associated with the 30-nm fiber. Moreover, the strong dependence of looping on linker histone concentration suggests a hierarchical self-association mechanism of relaxed nucleosome zigzag chains rather than longitudinal compaction as seen in 30-nm fibers. Specifically, concentrations lower than one linker histone per nucleosome promote self-associations and formation of these looped networks of zigzag fibers. The combined experimental and modeling evidence for condensed metaphase chromatin as hierarchical loops and bundles of relaxed zigzag nucleosomal chains rather than randomly coiled threads or straight and stiff helical fibers reconciles aspects of other models for higher-order chromatin structure; it constitutes not only an efficient storage form for the genomic material, consistent with other genome-wide chromosome conformation studies that emphasize looping, but also a convenient organization for local DNA unraveling and genome access.

Arginine-phosphate salt bridges between histones and DNA: intermolecular actuators that control nucleosome architecture

Structural bioinformatics and van der Waals density functional theory are combined to investigate the mechanochemical impact of a major class of histone-DNA interactions, namely, the formation of salt bridges between arginine residues in histones and phosphate groups on the DNA backbone. Principal component analysis reveals that the configurational fluctuations of the sugar-phosphate backbone display sequence-specific directionality and variability, and clustering of nucleosome crystal structures identifies two major salt-bridge configurations: a monodentate form in which the arginine end-group guanidinium only forms one hydrogen bond with the phosphate, and a bidentate form in which it forms two. Density functional theory calculations highlight that the combination of sequence, denticity, and salt-bridge positioning enables the histones to apply a tunable mechanochemical stress to the DNA via precise and specific activation of backbone deformations. The results suggest that selection for specific placements of van der Waals contacts, with high-precision control of the spatial distribution of intermolecular forces, may serve as an underlying evolutionary design principle for the structure and function of nucleosomes, a conjecture that is corroborated by previous experimental studies.

Chromatin Unfolding by Epigenetic Modifications Explained by Dramatic Impairment of Internucleosome Interactions: A Multiscale Computational Study

Histone tails and their epigenetic modifications play crucial roles in gene expression regulation by altering the architecture of chromatin. However, the structural mechanisms by which histone tails influence the interconversion between active and inactive chromatin remain unknown. Given the technical challenges in obtaining detailed experimental characterizations of the structure of chromatin, multiscale computations offer a promising alternative to model the effect of histone tails on chromatin folding. Here we combine multimicrosecond atomistic molecular dynamics simulations of dinucleosomes and histone tails in explicit solvent and ions, performed with three different state-of-the-art force fields and validated by experimental NMR measurements, with coarse-grained Monte Carlo simulations of 24-nucleosome arrays to describe the conformational landscape of histone tails, their roles in chromatin compaction, and the impact of lysine acetylation, a widespread epigenetic change, on both. We find that while the wild-type tails are highly flexible and disordered, the dramatic increase of secondary-structure order by lysine acetylation unfolds chromatin by decreasing tail availability for crucial fiber-compacting internucleosome interactions. This molecular level description of the effect of histone tails and their charge modifications on chromatin folding explains the sequence sensitivity and underscores the delicate connection between local and global structural and functional effects. Our approach also opens new avenues for multiscale processes of biomolecular complexes.

Changing chromatin fiber conformation by nucleosome repositioning

Chromatin conformation is dynamic and heterogeneous with respect to nucleosome positions, which can be changed by chromatin remodeling complexes in the cell. These molecular machines hydrolyze ATP to translocate or evict nucleosomes, and establish loci with regularly and more irregularly spaced nucleosomes as well as nucleosome-depleted regions. The impact of nucleosome repositioning on the three-dimensional chromatin structure is only poorly understood. Here, we address this issue by using a coarse-grained computer model of arrays of 101 nucleosomes considering several chromatin fiber models with and without linker histones, respectively. We investigated the folding of the chain in dependence of the position of the central nucleosome by changing the length of the adjacent linker DNA in basepair steps. We found in our simulations that these translocations had a strong effect on the shape and properties of chromatin fibers: i), Fiber curvature and flexibility at the center were largely increased and long-range contacts between distant nucleosomes on the chain were promoted. ii), The highest destabilization of the fiber conformation occurred for a nucleosome shifted by two basepairs from regular spacing, whereas effects of linker DNA changes of ?10 bp in phase with the helical twist of DNA were minimal. iii), A fiber conformation can stabilize a regular spacing of nucleosomes inasmuch as favorable stacking interactions between nucleosomes are facilitated. This can oppose nucleosome translocations and increase the energetic costs for chromatin remodeling. Our computational modeling framework makes it possible to describe the conformational heterogeneity of chromatin in terms of nucleosome positions, and thus advances theoretical models toward a better understanding of how genome compaction and access are regulated within the cell.

The chromatin fiber: multiscale problems and approaches

The structure of chromatin, affected by many factors from DNA linker lengths to posttranslational modifications, is crucial to the regulation of eukaryotic cells. Combined experimental and computational methods have led to new insights into its structural and dynamical features, from interactions due to the flexible core histone tails or linker histones to the physical mechanism driving the formation of chromosomal domains. Here we present a perspective of recent advances in chromatin modeling techniques at the atomic, mesoscopic, and chromosomal scales with a view toward developing multiscale computational strategies to integrate such findings. Innovative modeling methods that connect molecular to chromosomal scales are crucial for interpreting experiments and eventually deciphering the complex dynamic organization and function of chromatin in the cell.

Nucleosome positioning and composition modulate in silico chromatin flexibility

The dynamic organization of chromatin plays an essential role in the regulation of gene expression and in other fundamental cellular processes. The underlying physical basis of these activities lies in the sequential positioning, chemical composition, and intermolecular interactions of the nucleosomes-the familiar assemblies of ∼150 DNA base pairs and eight histone proteins-found on chromatin fibers. Here we introduce a mesoscale model of short nucleosomal arrays and a computational framework that make it possible to incorporate detailed structural features of DNA and histones in simulations of short chromatin constructs. We explore the effects of nucleosome positioning and the presence or absence of cationic N-terminal histone tails on the 'local' inter-nucleosomal interactions and the global deformations of the simulated chains. The correspondence between the predicted and observed effects of nucleosome composition and numbers on the long-range communication between the ends of designed nucleosome arrays lends credence to the model and to the molecular insights gleaned from the simulated structures. We also extract effective nucleosome-nucleosome potentials from the simulations and implement the potentials in a larger-scale computational treatment of regularly repeating chromatin fibers. Our results reveal a remarkable effect of nucleosome spacing on chromatin flexibility, with small changes in DNA linker length significantly altering the interactions of nucleosomes and the dimensions of the fiber as a whole. In addition, we find that these changes in nucleosome positioning influence the statistical properties of long chromatin constructs. That is, simulated chromatin fibers with the same number of nucleosomes exhibit polymeric behaviors ranging from Gaussian to worm-like, depending upon nucleosome spacing. These findings suggest that the physical and mechanical properties of chromatin can span a wide range of behaviors, depending on nucleosome positioning, and that care must be taken in the choice of models used to interpret the experimental properties of long chromatin fibers.

Interplay among nucleosomal DNA, histone tails, and corepressor CoREST underlies LSD1-mediated H3 demethylation

With its noncatalytic domains, DNA-binding regions, and a catalytic core targeting the histone tails, LSD1-CoREST (lysine-specific demethylase 1; REST corepressor) is an ideal model system to study the interplay between DNA binding and histone modification in nucleosome recognition. To this end, we covalently associated LSD1-CoREST to semisynthetic nucleosomal particles. This enabled biochemical and biophysical characterizations of nucleosome binding and structural elucidation by small-angle X-ray scattering, which was extensively validated through binding assays and site-directed mutagenesis of functional interfaces. Our results suggest that LSD1-CoREST functions as an ergonomic clamp that induces the detachment of the H3 histone tail from the nucleosomal DNA to make it available for capture by the enzyme active site. The key notion emerging from these studies is the inherently competitive nature of the binding interactions because nucleosome tails, chromatin modifiers, transcription factors, and DNA represent sites for multiple and often mutually exclusive interactions.

Quantitation of interactions between two DNA loops demonstrates loop domain insulation in E. coli cells

Eukaryotic gene regulation involves complex patterns of long-range DNA-looping interactions between enhancers and promoters, but how these specific interactions are achieved is poorly understood. Models that posit other DNA loops--that aid or inhibit enhancer-promoter contact--are difficult to test or quantitate rigorously in eukaryotic cells. Here, we use the well-characterized DNA-looping proteins Lac repressor and phage λ CI to measure interactions between pairs of long DNA loops in E. coli cells in the three possible topological arrangements. We find that side-by-side loops do not affect each other. Nested loops assist each other's formation consistent with their distance-shortening effect. In contrast, alternating loops, where one looping element is placed within the other DNA loop, inhibit each other's formation, thus providing clear support for the loop domain model for insulation. Modeling shows that combining loop assistance and loop interference can provide strong specificity in long-range interactions.

Chromatin fiber polymorphism triggered by variations of DNA linker lengths

Deciphering the factors that control chromatin fiber structure is key to understanding fundamental chromosomal processes. Although details remain unknown, it is becoming clear that chromatin is polymorphic depending on internal and external factors. In particular, different lengths of the linker DNAs joining successive nucleosomes (measured in nucleosome-repeat lengths or NRLs) that characterize different cell types and cell cycle stages produce different structures. NRL is also nonuniform within single fibers, but how this diversity affects chromatin fiber structure is not clear. Here we perform Monte Carlo simulations of a coarse-grained oligonucleosome model to help interpret fiber structure subject to intrafiber NRL variations, as relevant to proliferating cells of interphase chromatin, fibers subject to remodeling factors, and regulatory DNA sequences. We find that intrafiber NRL variations have a profound impact on chromatin structure, with a wide range of different architectures emerging (highly bent narrow forms, canonical and irregular zigzag fibers, and polymorphic conformations), depending on the NRLs mixed. This stabilization of a wide range of fiber forms might allow NRL variations to regulate both fiber compaction and selective DNA exposure. The polymorphic forms spanning canonical to sharply bent structures, like hairpins and loops, arise from large NRL variations and are surprisingly more compact than uniform NRL structures. They are distinguished by tail-mediated far-nucleosome interactions, in addition to the near-nucleosome interactions of canonical 30-nm fibers. Polymorphism is consistent with chromatin's diverse biological functions and heterogeneous constituents. Intrafiber NRL variations, in particular, may contribute to fiber bending and looping and thus to distant communication in associated regulatory processes.

DNA topology and transcription

Chromatin is a complex assembly that compacts DNA inside the nucleus while providing the necessary level of accessibility to regulatory factors conscripted by cellular signaling systems. In this superstructure, DNA is the subject of mechanical forces applied by variety of molecular motors. Rather than being a rigid stick, DNA possesses dynamic structural variability that could be harnessed during critical steps of genome functioning. The strong relationship between DNA structure and key genomic processes necessitates the study of physical constrains acting on the double helix. Here we provide insight into the source, dynamics, and biology of DNA topological domains in the eukaryotic cells and summarize their possible involvement in gene transcription. We emphasize recent studies that might inspire and impact future experiments on the involvement of DNA topology in cellular functions.

Establishing epigenetic domains via chromatin-bound histone modifiers

The eukaryotic nucleus harbors the DNA genome, which associates with histones and other chromosomal proteins into a complex referred to as chromatin. It provides an additional layer of so-called epigenetic information via histone modifications and DNA methylation on top of the DNA sequence that determines the cell's active gene expression program. The nucleus is devoid of internal organelles separated by membranes. Thus, free diffusive transport of proteins and RNA can occur throughout the space accessible for a given macromolecule. At the same time, chromatin is partitioned into different specialized structures such as nucleoli, chromosome territories, and heterochromatin domains that serve distinct functions. Here, we address the question of how the activity of chromatin-modifying enzymes is confined to chromatin subcompartments. We discuss mechanisms for establishing activity gradients of diffusive chromatin-modifying enzymes that could give rise to distinct chromatin domains within the cell nucleus. Interestingly, such gradients might directly result from immobilization of the enzymes on the flexible chromatin chain. Thus, locus-specific tethering of these enzymes to chromatin could have the potential to establish, maintain, or modulate epigenetic patterns of characteristic domain size.

Modification of enhancer chromatin: what, how, and why?

Emergence of form and function during embryogenesis arises in large part through cell-type- and cell-state-specific variation in gene expression patterns, mediated by specialized cis-regulatory elements called enhancers. Recent large-scale epigenomic mapping revealed unexpected complexity and dynamics of enhancer utilization patterns, with 400,000 putative human enhancers annotated by the ENCODE project alone. These large-scale efforts were largely enabled through the understanding that enhancers share certain stereotypical chromatin features. However, an important question still lingers: what is the functional significance of enhancer chromatin modification? Here we give an overview of enhancer-associated modifications of histones and DNA and discuss enzymatic activities involved in their dynamic deposition and removal. We describe potential downstream effectors of these marks and propose models for exploring functions of chromatin modification in regulating enhancer activity during development.

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.

Distant activation of transcription: mechanisms of enhancer action

Enhancers are regulatory DNA sequences that activate transcription over long distances. Recent studies revealed a widespread role of distant activation in eukaryotic gene regulation and in development of various human diseases, including cancer. Genomic and gene-targeted studies of enhancer action revealed novel mechanisms of transcriptional activation over a distance. They include formation of stable, inactive DNA-protein complexes at the enhancer and target promoter before activation, facilitated distant communication by looping of the spacer chromatin-covered DNA, and promoter activation by mechanisms that are different from classic recruiting. These studies suggest the similarity between the looping mechanisms involved in enhancer action on DNA in bacteria and in chromatin of higher organisms.

Mechanistic studies on histone catalyzed cleavage of apyrimidinic/apurinic sites in nucleosome core particles

Duplex DNA containing an apurinic/apyrimidinic (AP) lesion undergoes cleavage significantly more rapidly in nucleosome core particles (NCPs) than it does when free. The mechanism of AP cleavage within NCPs was studied through independently generating lesions within them. AP mediated DNA cleavage within NCPs is initiated by DNA-protein cross-link (DPC(un)) formation followed by β-elimination to give DPCs containing cleaved DNA (DPC(cl)). Hydrolysis of DPC(cl) produces a DNA single strand break (SSB). C2-dideuteration of AP showed that deprotonation from this position is involved in the rate-determining step. Experiments utilizing NCPs containing mutated histone H4 proteins indicated that lysine residues in the amino terminal tail are involved in both DPC formation and β-elimination steps. Lysines 16 and 20 seem to play a greater role in reacting with AP at superhelical location 1.5, but other amino acids (e.g., lysines 5, 8, and 12) compensate in their absence. The mechanism of rapid double strand breaks in bistranded, clustered AP lesions was studied by independently preparing reaction intermediates within model NCPs. A single strand break on one strand enhances the cleavage of a proximal AP on the opposite strand.

Insights into Gene Expression and Packaging from Computer Simulations

Within the nucleus of each cell lies DNA - an unfathomably long, twisted, and intricately coiled molecule - segments of which make up the genes that provide the instructions that a cell needs to operate. As we near the 60(th) anniversary of the discovery of the DNA double helix, crucial questions remain about how the physical arrangement of the DNA in cells affects how genes work. For example, how a cell stores the genetic information inside the nucleus is complicated by the necessity of maintaining accessibility to DNA for genetic processing. In order to gain insight into the roles played by various proteins in reading and compacting the genome, we have developed new methodologies to simulate the dynamic, three-dimensional structures of long, fluctuating, protein-decorated strands of DNA. Our a priori approach to the problem allows us to determine the effects of individual proteins and their chemical modifications on overall DNA structure and function. Here we present our recent treatment of the communication between regulatory proteins attached to precisely constructed stretches of chromatin. Our simulations account for the enhancement in communication detected experimentally on chromatin compared to protein-free DNA of the same chain length as well as the critical roles played by the cationic 'tails' of the histone proteins in this signaling. The states of chromatin captured in the simulations offer new insights into the ways that the DNA, histones, and regulatory proteins contribute to long-range communication along the genome.