Binding of DNA-bending non-histone proteins destabilizes regular 30-nm chromatin structure

Regulation of chromatin architecture by protein binding: insights from molecular modeling

Histone and non-histone proteins play key roles in the activation and repression of genes. In addition to experimental studies of their regulation of gene expression, molecular modeling at the nucleosome, chromatin, and chromosome levels can contribute insights into the molecular mechanisms involved. In this review, we provide an overview for protein-bound chromatin modeling, and describe how our group has integrated protein binding into genome systems across the scales, from all-atom to coarse-grained models, using explicit to implicit descriptions. We describe the associated applications to protein binding effects and biological mechanisms of genome folding and gene regulation. We end by illustrating the application of machine learning tools like AlphaFold2 to proteins relevant to chromatin systems.

Regulation of chromatin architecture by transcription factor binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome to control the transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers at nucleosome resolution to determine TF effect on chromatin architecture and shed light into its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. At the same time, the level of linker histone and tail acetylation, or different linker DNA lengths, regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining

Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.

An associative memory Hamiltonian model for DNA and nucleosomes

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

Liquid-like chromatin in the cell: What can we learn from imaging and computational modeling?

Chromatin in eukaryotic cells is a negatively charged long polymer consisting of DNA, histones, and various associated proteins. With its highly charged and heterogeneous nature, chromatin structure varies greatly depending on various factors (e.g. chemical modifications and protein enrichment) and the surrounding environment (e.g. cations): from a 10-nm fiber, a folded 30-nm fiber, to chromatin condensates/droplets. Recent advanced imaging has observed that chromatin exhibits a dynamic liquid-like behavior and undergoes structural variations within the cell. Current computational modeling has made it possible to reconstruct the liquid-like chromatin in the cell by dealing with a number of nucleosomes on multiscale levels and has become a powerful technique to inspect the molecular mechanisms giving rise to the observed behavior, which imaging methods cannot do on their own. Based on new findings from both imaging and modeling studies, we discuss the dynamic aspect of chromatin in living cells and its functional relevance.

Live imaging of chromatin distribution reveals novel principles of nuclear architecture and chromatin compartmentalization

Live imaging of chromatin in an intact organism reveals a novel mode of mesoscale chromatin organization at nuclear periphery.

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.

Physical Nature of Chromatin in the Nucleus

Genomic information is encoded on long strands of DNA, which are folded into chromatin and stored in a tiny nucleus. Nuclear chromatin is a negatively charged polymer composed of DNA, histones, and various nonhistone proteins. Because of its highly charged nature, chromatin structure varies greatly depending on the surrounding environment (e.g., cations, molecular crowding, etc.). New technologies to capture chromatin in living cells have been developed over the past 10 years. Our view on chromatin organization has drastically shifted from a regular and static one to a more variable and dynamic one. Chromatin forms numerous compact dynamic domains that act as functional units of the genome in higher eukaryotic cells and locally appear liquid-like. By changing DNA accessibility, these domains can govern various functions. Based on new evidences from versatile genomics and advanced imaging studies, we discuss the physical nature of chromatin in the crowded nuclear environment and how it is regulated.

Nucleosome positioning and chromatin organization

In our cells, DNA is folded and packed with the help of many proteins into chromatin whose basic unit is a nucleosome-DNA wrapped around octamer of histone proteins. The chain of nucleosomes is further folded and arranged into many layers and has a dynamic organization. How does the complex chromatin organization emerge from interactions among DNA, histones, and non-histone proteins have been a question of great interest. Here we review recent literature that investigated how nucleosome positioning and nucleosome-mediated interactions drive chromatin organization. Unlike our earlier understanding, chromatin is organized into 3D domains of various sizes having irregularly organized nucleosomes. These domains emerge due to heterogeneous nucleosome positioning and diverse inter-nucleosome interactions that vary in space and time.

Role of nucleosome positioning in 3D chromatin organization and loop formation

We present a physics-based polymer model that can investigate 3D organization of chromatin accounting for DNA elasticity, DNA-bending due to nucleosomes, and 1D organization of nucleosomes along DNA. We find that the packing density of chromatin oscillates between densities corresponding to highly folded and extended configurations as we change the nucleosome organization (length of linker DNA). We compute the looping probability of chromatin and show that the presence of nucleosomes increases the looping probability of the chain compared to that of a bare DNA. We also show that looping probability has a large variability depending on the nature of nucleosome organization and density of linker histones.

Irregular Chromatin: Packing Density, Fiber Width, and Occurrence of Heterogeneous Clusters

How chromatin is folded on the length scale of a gene is an open question. Recent experiments have suggested that, in vivo, chromatin is folded in an irregular manner and not as an ordered fiber with a width of 30 nm that is expected from theories of higher order packaging. Using computational methods, we examine how the interplay between DNA-bending nonhistone proteins, histone tails, intrachromatin electrostatic, and other interactions decide the nature of the packaging of chromatin. We show that although the DNA-bending nonhistone proteins make the chromatin irregular, they may not alter the packing density and size of the fiber. We find that the length of the interacting region and intrachromatin electrostatic interactions influence the packing density, clustering of nucleosomes, and the width of the chromatin fiber. Our results suggest that the heterogeneity in the interaction pattern will play an important role in deciding the nature of the packaging of chromatin.

Chromatin Structure and the Pioneering Transcription Factor FOXA1 Regulate TDG-Mediated Removal of 5-Formylcytosine from DNA

Although a functional relationship between active DNA demethylation and chromatin structure is often implied, direct experimental evidence is lacking. We investigated the relationship between chromatin structure and thymine DNA glycosylase (TDG) using chemically defined nucleosome arrays containing site-specifically positioned 5-formylcytosine (5fC) residues. We show that the extent of array compaction, as well as nucleosome positioning, dramatically influence the ability of TDG to excise 5fC from DNA, indicating that the chromatin structure is likely a key determinant of whether 5fC is removed from the genome or retained as an epigenetic mark. Furthermore, the H2A.Z/H3.3 double-variant nucleosome and the pioneering transcription factor forkhead box A1 (FOXA1), both of which are implicated in shaping the chromatin landscape during demethylation of tissue-specific enhancers, differentially regulate TDG activity on chromatin. Together, this work provides the first direct evidence that the higher order chromatin structure regulates active DNA demethylation through TDG and provides novel insights into the mechanism of 5fC turnover at enhancers.

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.

Chromatin Higher-Order Folding: A Perspective with Linker DNA Angles

The mechanism by which the "beads-on-a-string" nucleosome chain folds into various higher-order chromatin structures in eukaryotic cell nuclei is still poorly understood. The various models depicting higher-order chromatin as regular helical fibers and the very opposite "polymer melt" theory imply that interactions between nucleosome "beads" make the main contribution to the chromatin compaction. Other models in which the geometry of linker DNA "strings" entering and exiting the nucleosome define the three-dimensional structure predict that small changes in the linker DNA configuration may strongly affect nucleosome chain folding and chromatin higher-order structure. Among those studies, the cross-disciplinary approach pioneered by Jörg Langowski that combines computational modeling with biophysical and biochemical experiments was most instrumental for understanding chromatin higher-order structure in vitro. Strikingly, many recent studies, including genome-wide nucleosome interaction mapping and chromatin imaging, show an excellent agreement with the results of three-dimensional computational modeling based on the primary role of linker DNA geometry in chromatin compaction. This perspective relates nucleosome array models with experimental studies of nucleosome array folding in vitro and in situ. I argue that linker DNA configuration plays a key role in determining nucleosome chain flexibility, topology, and propensity for self-association, thus providing new implications for regulation of chromatin accessibility to DNA binding factors and RNA transcription machinery as well as long-range communications between distant genomic sites.

Regulation of chromatin folding by conformational variations of nucleosome linker DNA

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