Polymer Modeling Reveals Interplay between Physical Properties of Chromosomal DNA and the Size and Distribution of Condensin-Based Chromatin Loops

Transient DNA loops occur throughout the genome due to thermal fluctuations of DNA and the function of SMC complex proteins such as condensin and cohesin. Transient crosslinking within and between chromosomes and loop extrusion by SMCs have profound effects on high-order chromatin organization and exhibit specificity in cell type, cell cycle stage, and cellular environment. SMC complexes anchor one end to DNA with the other extending some distance and retracting to form a loop. How cells regulate loop sizes and how loops distribute along chromatin are emerging questions. To understand loop size regulation, we employed bead-spring polymer chain models of chromatin and the activity of an SMC complex on chromatin. Our study shows that (1) the stiffness of the chromatin polymer chain, (2) the tensile stiffness of chromatin crosslinking complexes such as condensin, and (3) the strength of the internal or external tethering of chromatin chains cooperatively dictate the loop size distribution and compaction volume of induced chromatin domains. When strong DNA tethers are invoked, loop size distributions are tuned by condensin stiffness. When DNA tethers are released, loop size distributions are tuned by chromatin stiffness. In this three-way interaction, the presence and strength of tethering unexpectedly dictates chromatin conformation within a topological domain.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Chromosome Folding Promotes Intrachromosomal Aberrations under Radiation- and Nuclease-Induced DNA Breakage

The long-standing question in radiation and cancer biology is how principles of chromosome organization impact the formation of chromosomal aberrations (CAs). To address this issue, we developed a physical modeling approach and analyzed high-throughput genomic data from chromosome conformation capture (Hi-C) and translocation sequencing (HTGTS) methods. Combining modeling of chromosome structure and of chromosomal aberrations induced by ionizing radiation (IR) and nuclease we made predictions which quantitatively correlated with key experimental findings in mouse chromosomes: chromosome contact maps, high frequency of cis-translocation breakpoints far outside of the site of nuclease-induced DNA double-strand breaks (DSBs), the distinct shape of breakpoint distribution in chromosomes with different 3D organizations. These correlations support the heteropolymer globule principle of chromosome organization in G1-arrested pro-B mouse cells. The joint analysis of Hi-C, HTGTS and physical modeling data offers mechanistic insight into how chromosome structure heterogeneity, globular folding and lesion dynamics drive IR-recurrent CAs. The results provide the biophysical and computational basis for the analysis of chromosome aberration landscape under IR and nuclease-induced DSBs.

CHARMM-GUI Polymer Builder for Modeling and Simulation of Synthetic Polymers

Molecular modeling and simulations are invaluable tools for polymer science and engineering, which predict physicochemical properties of polymers and provide molecular-level insight into the underlying mechanisms. However, building realistic polymer systems is challenging and requires considerable experience because of great variations in structures as well as length and time scales. This work describes Polymer Builder in CHARMM-GUI (http://www.charmm-gui.org/input/polymer), a web-based infrastructure that provides a generalized and automated process to build a relaxed polymer system. Polymer Builder not only provides versatile modeling methods to build complex polymer structures, but also generates realistic polymer melt and solution systems through the built-in coarse-grained model and all-atom replacement. The coarse-grained model parametrization is generalized and extensively validated with various experimental data and all-atom simulations. In addition, the capability of Polymer Builder for generating relaxed polymer systems is demonstrated by density calculations of 34 homopolymer melt systems, characteristic ratio calculations of 170 homopolymer melt systems, a morphology diagram of poly(styrene-b-methyl methacrylate) block copolymers, and self-assembly behavior of amphiphilic poly(ethylene oxide-b-ethylethane) block copolymers in water. We hope that Polymer Builder is useful to carry out innovative and novel polymer modeling and simulation research to acquire insight into structures, dynamics, and underlying mechanisms of complex polymer-containing systems.

Are the processes of DNA replication and DNA repair reading a common structural chromatin unit?

Decades of investigation on genomic DNA have brought us deeper insights into its organization within the nucleus and its metabolic mechanisms. This was fueled by the parallel development of experimental techniques and has stimulated model building to simulate genome conformation in agreement with the experimental data. Here, we will discuss our recent discoveries on the chromatin units of DNA replication and DNA damage response. We will highlight their remarkable structural similarities and how both revealed themselves as clusters of nanofocal structures each on the hundred thousand base pair size range corresponding well with chromatin loop sizes. We propose that the function of these two global genomic processes is determined by the loop level organization of chromatin structure with structure dictating function.

Abbreviations: 3D-SIM: 3D-structured illumination microscopy; 3C: chromosome conformation capture; DDR: DNA damage response; FISH: fluorescent in situ hybridization; Hi-C: high conformation capture; HiP-HoP: highly predictive heteromorphic polymer model; IOD: inter-origin distance; LAD: lamina associated domain; STED: stimulated emission depletion microscopy; STORM: stochastic optical reconstruction microscopy; SBS: strings and binders switch model; TAD: topologically associated domain

Polymer Injectivity Test Design Using Numerical Simulation

Polymer flooding is an enhanced oil recovery (EOR) process, which has received increasing interest in the industry. In this process, water-soluble polymers are used to increase injected water viscosity in order to improve mobility ratio and hence improve reservoir sweep. Polymer solutions are non-Newtonian fluids, i.e., their viscosities are shear dependent. Polymers may exhibit an increase in viscosity at high shear rates in porous media, which can cause injectivity loss. In contrast, at low shear rates they may observe viscosity loss and hence enhance the injectivity. Therefore, due to the complex non-Newtonian rheology of polymers, it is necessary to optimize the design of polymer injectivity tests in order to improve our understanding of the rheology behavior and enhance the design of polymer flood projects. This study has been addressing what information that can be gained from polymer injectivity tests, and how to design the test for maximizing information. The main source of information in the field is from the injection bottom-hole pressure (BHP). Simulation studies have analyzed the response of different non-Newtonian rheology on BHP with variations of rate and time. The results have shown that BHP from injectivity tests can be used to detect in-situ polymer rheology.

Polymer Simulations of Heteromorphic Chromatin Predict the 3D Folding of Complex Genomic Loci

Chromatin folded into 3D macromolecular structures is often analyzed by chromosome conformation capture (3C) and fluorescence in situ hybridization (FISH) techniques, but these frequently provide contradictory results. Chromatin can be modeled as a simple polymer composed of a connected chain of units. By embedding data for epigenetic marks (H3K27ac), chromatin accessibility (assay for transposase-accessible chromatin using sequencing [ATAC-seq]), and structural anchors (CCCTC-binding factor [CTCF]), we developed a highly predictive heteromorphic polymer (HiP-HoP) model, where the chromatin fiber varied along its length; combined with diffusing protein bridges and loop extrusion, this model predicted the 3D organization of genomic loci at a population and single-cell level. The model was validated at several gene loci, including the complex Pax6 gene, and was able to determine locus conformations across cell types with varying levels of transcriptional activity and explain different mechanisms of enhancer use. Minimal a priori knowledge of epigenetic marks is sufficient to recapitulate complex genomic loci in 3D and enable predictions of chromatin folding paths.

•

HiP-HoP: highly predictive heteromorphic polymer model to analyze chromatin structure

•

Polymer simulations use widely available epigenetic and protein binding data as input

•

Validate HiP-HoP model at complex loci using 3D FISH and Capture-C

•

Simulations uncover striking conformational variability in chromatin fiber folding