How epigenome drives chromatin folding and dynamics, insights from efficient coarse-grained models of chromosomes

Nonequilibrium switching of segmental states can influence compaction of chromatin

Knowledge about the dynamic nature of chromatin organization is essential to understand the regulation of processes like DNA transcription and repair. The existing models of chromatin assume that protein organization and chemical states along chromatin are static and the 3D organization is purely a result of protein-mediated intra-chromatin interactions. Here we present a new hypothesis that certain nonequilibrium processes, such as switching of chemical and physical states due to nucleosome assembly/disassembly or gene repression/activation, can also simultaneously influence chromatin configurations. To understand the implications of this inherent nonequilibrium switching, we present a block copolymer model of chromatin, with switching of its segmental states between two states, mimicking active/repressed or protein unbound/bound states. We show that competition between switching timescale Tt, polymer relaxation timescale τp, and segmental relaxation timescale τs can lead to non-trivial changes in chromatin organization, leading to changes in local compaction and contact probabilities. As a function of the switching timescale, the radius of gyration of chromatin shows a non-monotonic behavior with a prominent minimum when Tt ≈ τp and a maximum when Tt ≈ τs. We find that polymers with a small segment length exhibit a more compact structure than those with larger segment lengths. We also find that the switching can lead to higher contact probability and better mixing of far-away segments. Our study also shows that the nature of the distribution of chromatin clusters varies widely as we change the switching rate.

Dose and DNA damage modelling of diffusing alpha-emitters radiation therapy using Geant4

PURPOSE

Diffusing alpha-emitters radiation therapy (DaRT) is a brachytherapy technique using α-particles to treat solid tumours. The high linear energy transfer (LET) and short range of α-particles make them good candidates for the targeted treatment of cancer. Treatment planning of DaRT requires a good understanding of the dose from α-particles and the other particles released in the 224Ra decay chain.

METHODS

The Geant4 Monte Carlo toolkit has been used to simulate a DaRT seed to better understand the dose contribution from all particles and simulate the DNA damage due to this treatment.

RESULTS

Close to the seed α-particles deliver the majority of dose, however at radial distances greater than 4 mm, the contribution of β-particles is greater. The RBE has been estimated as a function of number of double strand breaks (DSBs) and complex DSBs. A maximum seed spacing of 5.5 mm and 6.5 mm was found to deliver at least 20 Gy RBE weighted dose between the seeds for RBEDSB and RBEcDSB respectively.

CONCLUSIONS

The DNA damage changes with radial distance from the seed and has been found to become less complex with distance, which is potentially easier for the cell to repair. Close to the seed α-particles contribute the majority of dose, however the contribution from other particles cannot be neglected and may influence the choice of seed spacing.

Effective modeling of the chromatin structure by coarse-grained methods

The interphase chromatin structure is extremely complex, precise and dynamic. Experimental methods can only show the frequency of interaction of the various parts of the chromatin. Therefore, it is extremely important to develop theoretical methods to predict the chromatin structure. In this publication, we implemented an extended version of the SBS model described by Barbieri et al. and created the ChroMC program that is easy to use and freely available (https://github.com/regulomics/chroMC) to other users. We also describe the necessary factors for the effective modeling of the chromatin structure in Drosophila melanogaster. We compared results of chromatin structure predictions using two methods: Monte Carlo and Molecular Dynamic. Our simulations suggest that incorporating black, non-reactive chromatin is necessary for successful prediction of chromatin structure, while the loop extrusion model with a long range attraction potential or Lennard-Jones (with local attraction force) as well as using Hi-C data as input are not essential for the basic structure reconstruction. We also proposed a new way to calculate the similarity of the properties of contact maps including the calculation of local similarity.Communicated by Ramaswamy H. Sarma.

Crumpled polymer with loops recapitulates key features of chromosome organization

Chromosomes are exceedingly long topologically-constrained polymers compacted in a cell nucleus. We recently suggested that chromosomes are organized into loops by an active process of loop extrusion. Yet loops remain elusive to direct observations in living cells; detection and characterization of myriads of such loops is a major challenge. The lack of a tractable physical model of a polymer folded into loops limits our ability to interpret experimental data and detect loops. Here, we introduce a new physical model - a polymer folded into a sequence of loops, and solve it analytically. Our model and a simple geometrical argument show how loops affect statistics of contacts in a polymer across different scales, explaining universally observed shapes of the contact probability. Moreover, we reveal that folding into loops reduces the density of topological entanglements, a novel phenomenon we refer as "the dilution of entanglements". Supported by simulations this finding suggests that up to ~ 1 - 2Mb chromosomes with loops are not topologically constrained, yet become crumpled at larger scales. Our theoretical framework allows inference of loop characteristics, draws a new picture of chromosome organization, and shows how folding into loops affects topological properties of crumpled polymers.

HP1-driven phase separation recapitulates the thermodynamics and kinetics of heterochromatin condensate formation

The spatial segregation of pericentromeric heterochromatin (PCH) into distinct, membrane-less nuclear compartments involves the binding of Heterochromatin Protein 1 (HP1) to H3K9me2/3-rich genomic regions. While HP1 exhibits liquid-liquid phase separation properties in vitro, its mechanistic impact on the structure and dynamics of PCH condensate formation in vivo remains largely unresolved. Here, using a minimal theoretical framework, we systematically investigate the mutual coupling between self-interacting HP1-like molecules and the chromatin polymer. We reveal that the specific affinity of HP1 for H3K9me2/3 loci facilitates coacervation in nucleo and promotes the formation of stable PCH condensates at HP1 levels far below the concentration required to observe phase separation in purified protein assays in vitro. These heterotypic HP1-chromatin interactions give rise to a strong dependence of the nucleoplasmic HP1 density on HP1-H3K9me2/3 stoichiometry, consistent with the thermodynamics of multicomponent phase separation. The dynamical cross talk between HP1 and the viscoelastic chromatin scaffold also leads to anomalously slow equilibration kinetics, which strongly depend on the genomic distribution of H3K9me2/3 domains and result in the coexistence of multiple long-lived, microphase-separated PCH compartments. The morphology of these complex coacervates is further found to be governed by the dynamic establishment of the underlying H3K9me2/3 landscape, which may drive their increasingly abnormal, aspherical shapes during cell development. These findings compare favorably to 4D microscopy measurements of HP1 condensate formation in live Drosophila embryos and suggest a general quantitative model of PCH formation based on the interplay between HP1-based phase separation and chromatin polymer mechanics.

In-silico calculations of DNA damage induced by α-particles in the 224Ra DaRT decay chain for a better understanding of the radiobiological effectiveness of this treatment

Diffusing alpha-emitters radiation Therapy (DaRT) is an interstitial brachytherapy technique using 224Ra seeds. For accurate treatment planning a good understanding of the early DNA damage due to α-particles is required. Geant4-DNA was used to calculate the initial DNA damage and radiobiological effectiveness due to α-particles with linear energy transfer (LET) values in the range 57.5-225.9 keV/μm from the 224Ra decay chain. The impact of DNA base pair density on DNA damage has been modelled, as this parameter varies between human cell lines. Results show that the quantity and complexity of DNA damage changes with LET as expected. Indirect damage, due to water radical reactions with the DNA, decreases and becomes less significant at higher LET values as shown in previous studies. As expected, the yield of complex double strand breaks (DSBs), which are harder for a cell to repair, increases approximately linearly with LET. The level of complexity of DSBs and radiobiological effectiveness have been found to increase with LET as expected. The quantity of DNA damage has been shown to increase for increased DNA density in the expected base pair density range of human cells. The change in damage yield as a function of base pair density is largest for higher LET α-particles, an increase of over 50% for individual strand breaks between 62.7 and 127.4 keV/μm. This change in yield shows that the DNA base pair density is an important parameter for modelling DNA damage particularly at higher LET where the DNA damage is greatest and most complex.

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.

Epigenetic switching with asymmetric bridging interactions

Gene expression states are often stably sustained in cis despite massively disruptive events like DNA replication. This is achieved by on-going enzymatic activity that maintains parts of the DNA in either heterochromatic (packed) or euchromatic (free) states, each of which is stabilized by both positive feedback and bridging interactions between individual nucleosomes. In contrast to condensed matter, however, the dynamics is not only governed by equilibrium binding interactions but is also mediated by enzymes that recognize and act on specific amino acid tails of the nucleosomes. The mechanical result is that some nucleosomes can bind to one another and form tightly packed polymer configurations, whereas others remain unbound and form free, noncompact polymer configurations. Here, we study the consequences of such an asymmetric interaction pattern on the dynamics of epigenetic switching. We develop a 3D polymer model and show that traits associated with epigenetic switching, such as bistability and epigenetic memory, are permitted by such a model. We find, however, that the experimentally observed burst-like nature of some epigenetic switches is difficult to reproduce by this biologically motivated interaction. Instead, the behavior seen in experiments can be explained by introducing partial confinement, which particularly affects the euchromatic regions of the chromosome.

Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin

BACKGROUND

Interactions among topologically associating domains (TADs), and between the nuclear envelope (NE) and lamina-associated domains (LADs) are expected to shape various aspects of three-dimensional (3D) chromatin structure and dynamics; however, relevant genome-wide experiments that may provide statistically significant conclusions remain difficult.

RESULTS

We have developed a coarse-grained dynamical model of D. melanogaster nuclei at TAD resolution that explicitly accounts for four distinct epigenetic classes of TADs and LAD-NE interactions. The model is parameterized to reproduce the experimental Hi-C map of the wild type (WT) nuclei; it describes time evolution of the chromatin over the G1 phase of the interphase. The simulations include an ensemble of nuclei, corresponding to the experimentally observed set of several possible mutual arrangements of chromosomal arms. The model is validated against multiple structural features of chromatin from several different experiments not used in model development. Predicted positioning of all LADs at the NE is highly dynamic-the same LAD can attach, detach and move far away from the NE multiple times during interphase. The probabilities of LADs to be in contact with the NE vary by an order of magnitude, despite all having the same affinity to the NE in the model. These probabilities are mostly determined by a highly variable local linear density of LADs along the genome, which also has the same strong effect on the predicted positioning of individual TADs -- higher probability of a TAD to be near NE is largely determined by a higher linear density of LADs surrounding this TAD. The distribution of LADs along the chromosome chains plays a notable role in maintaining a non-random average global structure of chromatin. Relatively high affinity of LADs to the NE in the WT nuclei substantially reduces sensitivity of the global radial chromatin distribution to variations in the strength of TAD-TAD interactions compared to the lamin depleted nuclei, where a small (0.5 kT) increase of cross-type TAD-TAD interactions doubles the chromatin density in the central nucleus region.

CONCLUSIONS

A dynamical model of the entire fruit fly genome makes multiple genome-wide predictions of biological interest. The distribution of LADs along the chromatin chains affects their probabilities to be in contact with the NE and radial positioning of highly mobile TADs, playing a notable role in creating a non-random average global structure of the chromatin. We conjecture that an important role of attractive LAD-NE interactions is to stabilize global chromatin structure against inevitable cell-to-cell variations in TAD-TAD interactions.

Mesoscale, long-time mixing of chromosomes and its connection to polymer dynamics

Chromosomes are arranged in distinct territories within the nucleus of animal cells. Recent experiments have shown that these territories overlap at their edges, suggesting partial mixing during interphase. Experiments that knock-down of condensin II proteins during interphase indicate increased chromosome mixing, which demonstrates control of the mixing. In this study, we use a generic polymer simulation to quantify the dynamics of chromosome mixing over time. We introduce the chromosome mixing index, which quantifies the mixing of distinct chromosomes in the nucleus. We find that the chromosome mixing index in a small confinement volume (as a model of the nucleus), increases as a power-law of the time, with the scaling exponent varying non-monotonically with self-interaction and volume fraction. By comparing the chromosome mixing index with both monomer subdiffusion due to (non-topological) intermingling of chromosomes as well as even slower reptation, we show that for relatively large volume fractions, the scaling exponent of the chromosome mixing index is related to Rouse dynamics for relatively weak chromosome attractions and to reptation for strong attractions. In addition, we extend our model to more realistically account for the situation of the Drosophila chromosome by including the heterogeneity of the polymers and their lengths to account for microphase separation of euchromatin and heterochromatin and their interactions with the nuclear lamina. We find that the interaction with the lamina further impedes chromosome mixing.

Disruption of polyhomeotic polymerization decreases nucleosome occupancy and alters genome accessibility

Chromatin attains its three-dimensional (3D) conformation by establishing contacts between different noncontiguous regions. Sterile Alpha Motif (SAM)-mediated polymerization of the polyhomeotic (PH) protein regulates subnuclear clustering of Polycomb Repressive Complex 1 (PRC1) and chromatin topology. The mutations that perturb the ability of the PH to polymerize, disrupt long-range chromatin contacts, alter Hox gene expression, and lead to developmental defects. To understand the underlying mechanism, we combined the experiments and theory to investigate the effect of this SAM domain mutation on nucleosome occupancy and accessibility on a genome wide scale. Our data show that disruption of PH polymerization because of SAM domain mutation decreases nucleosome occupancy and alters accessibility. Polymer simulations investigating the interplay between distant chromatin contacts and nucleosome occupancy, both of which are regulated by PH polymerization, suggest that nucleosome density increases when contacts between different regions of chromatin are established. Taken together, it appears that SAM domain-mediated PH polymerization biomechanically regulates the organization of chromatin at multiple scales from nucleosomes to chromosomes and we suggest that higher order organization can have a top-down causation effect on nucleosome occupancy.

4D epigenomics: deciphering the coupling between genome folding and epigenomic regulation with biophysical modeling

Recent experimental observations suggest a strong coupling between the 3D nuclear chromosome organization and epigenomics. However, the mechanistic and functional bases of such interplay remain elusive. In this review, we describe how biophysical modeling has been instrumental in characterizing how genome folding may impact the formation of epigenomic domains and, conversely, how epigenomic marks may affect chromosome conformation. Finally, we discuss how this mutual feedback loop between chromatin organization and epigenome regulation, via the formation of physicochemical nanoreactors, may represent a key functional role of 3D compartmentalization in the assembly and maintenance of stable - but yet plastic - epigenomic landscapes.

Multivalent binding proteins can drive collapse and reswelling of chromatin in confinement

Collapsed conformations of chromatin have been long suspected of being mediated by interactions with multivalent binding proteins, which can bring together distant sections of the chromatin fiber. In this study, we use Langevin dynamics simulation of a coarse grained chromatin polymer to show that the role of binding proteins can be more nuanced than previously suspected. In particular, for chromatin polymer in confinement, entropic forces can drive reswelling of collapsed chromatin with increasing binder concentrations, and this reswelling transition happens at physiologically relevant binder concentrations. Both the extent of collapse, and also of reswelling depends on the strength of confinement. We also study the kinetics of collapse and reswelling and show that both processes occur in similar timescales. We characterise this reswelling of chromatin in biologically relevant regimes and discuss the non-trivial role of multivalent binding proteins in mediating the spatial organisation of the genome.

Integrative studies of 3D genome organization and chromatin structure

The structural organization of the genome is emerging as a crucial regulator of the cell state, affecting gene transcription, DNA replication, and repair. Over the last twenty years, increasing evidence prompted the development of new experimental techniques to study genome structure. In parallel with the complexity of the novel techniques, computational approaches have become an essential tool in any structural genomics laboratory to analyze and model the data. For biologists to be able to apply the most appropriate modeling approach, it is fundamental to understand the conceptual bases of distinct methods and the insights they can provide. Here, we will discuss recent advances that were possible thanks to 3D genome modeling, discuss their limitations and highlight future perspectives.

3DPolyS-LE: an accessible simulation framework to model the interplay between chromatin and loop extrusion

SUMMARY

Recent studies suggest that the loop extrusion activity of Structural Maintenance of Chromosomes complexes is central to proper organization of genomes in vivo. Polymer physics-based modeling of chromosome structure has been instrumental to assess which structures such extrusion can create. Only few laboratories however have the technical and computational expertise to create in silico models combining dynamic features of chromatin and loop extruders. Here, we present 3DPolyS-LE, a self-contained, easy to use modeling and simulation framework allowing non-specialists to ask how specific properties of loop extruders and boundary elements impact on 3D chromosome structure. 3DPolyS-LE also provides algorithms to compare predictions with experimental Hi-C data.

AVAILABILITY AND IMPLEMENTATION

Software available at https://gitlab.com/togop/3DPolyS-LE; implemented in Python and Fortran 2003 and supported on any Unix-based operating system (Linux and Mac OS).

SUPPLEMENTARY INFORMATION

Supplementary information are available at Bioinformatics online.

Multiple parameters shape the 3D chromatin structure of single nuclei at the doc locus in Drosophila

The spatial organization of chromatin at the scale of topologically associating domains (TADs) and below displays large cell-to-cell variations. Up until now, how this heterogeneity in chromatin conformation is shaped by chromatin condensation, TAD insulation, and transcription has remained mostly elusive. Here, we used Hi-M, a multiplexed DNA-FISH imaging technique providing developmental timing and transcriptional status, to show that the emergence of TADs at the ensemble level partially segregates the conformational space explored by single nuclei during the early development of Drosophila embryos. Surprisingly, a substantial fraction of nuclei display strong insulation even before TADs emerge. Moreover, active transcription within a TAD leads to minor changes to the local inter- and intra-TAD chromatin conformation in single nuclei and only weakly affects insulation to the neighboring TAD. Overall, our results indicate that multiple parameters contribute to shaping the chromatin architecture of single nuclei at the TAD scale.

Painters in chromatin: a unified quantitative framework to systematically characterize epigenome regulation and memory

In eukaryotes, many stable and heritable phenotypes arise from the same DNA sequence, owing to epigenetic regulatory mechanisms relying on the molecular cooperativity of 'reader-writer' enzymes. In this work, we focus on the fundamental, generic mechanisms behind the epigenome memory encoded by post-translational modifications of histone tails. Based on experimental knowledge, we introduce a unified modeling framework, the painter model, describing the mechanistic interplay between sequence-specific recruitment of chromatin regulators, chromatin-state-specific reader-writer processes and long-range spreading mechanisms. A systematic analysis of the model building blocks highlights the crucial impact of tridimensional chromatin organization and state-specific recruitment of enzymes on the stability of epigenomic domains and on gene expression. In particular, we show that enhanced 3D compaction of the genome and enzyme limitation facilitate the formation of ultra-stable, confined chromatin domains. The model also captures how chromatin state dynamics impact the intrinsic transcriptional properties of the region, slower kinetics leading to noisier expression. We finally apply our framework to analyze experimental data, from the propagation of γH2AX around DNA breaks in human cells to the maintenance of heterochromatin in fission yeast, illustrating how the painter model can be used to extract quantitative information on epigenomic molecular processes.

Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains

Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.

Polymer Modeling of 3D Epigenome Folding: Application to Drosophila

Mechanistic modeling in biology allows to investigate, based on first principles, if putative hypotheses are compatible with observations and to drive further experimental works. Along this line, polymer modeling has been instrumental in 3D genomics to better understand the impact of key mechanisms on the spatial genome organization. Here, I describe how polymer-based models can be practically used to study the role of epigenome in chromosome folding. I illustrate this methodology in the context of Drosophila epigenome folding.

The Physical Behavior of Interphase Chromosomes: Polymer Theory and Coarse-Grain Computer Simulations

Fluorescence in situ hybridization and chromosome conformation capture methods point to the same conclusion: that chromosomes appear to the external observer as compact structures with a highly nonrandom three-dimensional organization. In this work, we recapitulate the efforts made by us and other groups to rationalize this behavior in terms of the mathematical language and tools of polymer physics. After a brief introduction dedicated to some crucial experiments dissecting the structure of interphase chromosomes, we discuss at a nonspecialistic level some fundamental aspects of theoretical and numerical polymer physics. Then, we inglobe biological and polymer aspects into a polymer model for interphase chromosomes which moves from the observation that mutual topological constraints, such as those typically present between polymer chains in ordinary melts, induce slow chain dynamics and "constraint" chromosomes to resemble double-folded randomly branched polymer conformations. By explicitly turning these ideas into a multi-scale numerical algorithm which is described here in full details, we can design accurate model polymer conformations for interphase chromosomes and offer them for systematic comparison to experiments. The review is concluded by discussing the limitations of our approach and pointing to promising perspectives for future work.

Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin

Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we showed that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling interpretation of the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the "liquid" or "solid" state of chromatin are in fact two sides of the same coin.

Towards the characterization of neutron carcinogenesis through direct action simulations of clustered DNA damage

Neutron exposure poses a unique radiation protection concern because neutrons have a large, energy-dependent relative biological effectiveness (RBE) for stochastic effects. Recent computational studies on the microdosimetric properties of neutron dose deposition have implicated clustered DNA damage as a likely contributor to this marked energy dependence. So far, publications have focused solely on neutron RBE for inducing clusters of DNA damage containing two or more DNA double strand breaks (DSBs). In this study, we have conducted a novel assessment of neutron RBE for inducing all types of clustered DNA damage that contain two or more lesions, stratified by whether the clusters contain DSBs (complex DSB clusters) or not (non-DSB clusters). This assessment was conducted for eighteen initial neutron energies between 1 eV and 10 MeV as well as a reference radiation of 250 keV x-rays. We also examined the energy dependence of cluster length and cluster complexity because these factors are believed to impact the DNA repair process. To carry out our investigation, we developed a user-friendly TOPAS-nBio application that includes a custom nuclear DNA model and a novel algorithm for recording clustered DNA damage. We found that neutron RBE for inducing complex DSB clusters exhibited similar energy dependence to the canonical neutron RBE for stochastic radiobiological effects, at multiple depths in human tissue. Qualitatively similar results were obtained for non-DSB clusters, although the quantitative agreement was lower. Additionally we identified a significant neutron energy dependence in the average length and complexity of clustered lesions. These results support the idea that many types of clustered DNA damage contribute to the energy dependence of neutron RBE for stochastic radiobiological effects and imply that the size and constituent lesions of individual clusters should be taken into account when modeling DNA repair. Our results were qualitatively consistent for (i) multiple radiation doses (including a low-dose 0.1 Gy irradiation), (ii) variations in the maximal lesion separation distance used to define a cluster, and (iii) two distinct collections of physics models used to govern particle transport. Our complete TOPAS-nBio application has been released under an open-source license to enable others to independently validate our work and to expand upon it.

A Geant4-DNA Evaluation of Radiation-Induced DNA Damage on a Human Fibroblast

Simple Summary

DNA damage caused by ionizing radiation in a human fibroblast cell evaluated by the Geant4-DNA Monte Carlo toolkit is presented. A validation study using a computational geometric human DNA model was then carried out, and the calculated DNA damage as a function of particle type and energy is presented. The results of this work showed a significant improvement on past work and were consistent with recent radiobiological experimental data, such as damage yields. This work and the developed methodology could impact a broad number of research fields in which the understanding of radiation effects is crucial, such as cancer radiotherapy, space science, and medical physics.

Abstract

Accurately modeling the radiobiological mechanisms responsible for the induction of DNA damage remains a major scientific challenge, particularly for understanding the effects of low doses of ionizing radiation on living beings, such as the induction of carcinogenesis. A computational approach based on the Monte Carlo technique to simulate track structures in a biological medium is currently the most reliable method for calculating the early effects induced by ionizing radiation on DNA, the primary cellular target of such effects. The Geant4-DNA Monte Carlo toolkit can simulate not only the physical, but also the physico-chemical and chemical stages of water radiolysis. These stages can be combined with simplified geometric models of biological targets, such as DNA, to assess direct and indirect early DNA damage. In this study, DNA damage induced in a human fibroblast cell was evaluated using Geant4-DNA as a function of incident particle type (gammas, protons, and alphas) and energy. The resulting double-strand break yields as a function of linear energy transfer closely reproduced recent experimental data. Other quantities, such as fragment length distribution, scavengeable damage fraction, and time evolution of damage within an analytical repair model also supported the plausibility of predicting DNA damage using Geant4-DNA.The complete simulation chain application “molecularDNA”, an example for users of Geant4-DNA, will soon be distributed through Geant4.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

3D genome organisation in Drosophila

Ever since Thomas Hunt Morgan's discovery of the chromosomal basis of inheritance by using Drosophila melanogaster as a model organism, the fruit fly has remained an essential model system in studies of genome biology, including chromatin organisation. Very much as in vertebrates, in Drosophila, the genome is organised in territories, compartments and topologically associating domains (TADs). However, these domains might be formed through a slightly different mechanism than in vertebrates due to the presence of a large and potentially redundant set of insulator proteins and the minor role of dCTCF in TAD boundary formation. Here, we review the different levels of chromatin organisation in Drosophila and discuss mechanisms and factors that might be involved in TAD formation. The dynamics of TADs and enhancer-promoter interactions in the context of transcription are covered in the light of currently conflicting results. Finally, we illustrate the value of polymer modelling approaches to infer the principles governing the three-dimensional organisation of the Drosophila genome.

Live imaging and biophysical modeling support a button-based mechanism of somatic homolog pairing in Drosophila

Three-dimensional eukaryotic genome organization provides the structural basis for gene regulation. In Drosophila melanogaster, genome folding is characterized by somatic homolog pairing, where homologous chromosomes are intimately paired from end to end; however, how homologs identify one another and pair has remained mysterious. Recently, this process has been proposed to be driven by specifically interacting 'buttons' encoded along chromosomes. Here, we turned this hypothesis into a quantitative biophysical model to demonstrate that a button-based mechanism can lead to chromosome-wide pairing. We tested our model using live-imaging measurements of chromosomal loci tagged with the MS2 and PP7 nascent RNA labeling systems. We show solid agreement between model predictions and experiments in the pairing dynamics of individual homologous loci. Our results strongly support a button-based mechanism of somatic homolog pairing in Drosophila and provide a theoretical framework for revealing the molecular identity and regulation of buttons.

Multi-Scale Organization of the Drosophila melanogaster Genome

Interphase chromatin, despite its appearance, is a highly organized framework of loops and bends. Chromosomes are folded into topologically associating domains, or TADs, and each chromosome and its homolog occupy a distinct territory within the nucleus. In Drosophila, genome organization is exceptional because homologous chromosome pairing is in both germline and somatic tissues, which promote interhomolog interactions such as transvection that can affect gene expression in trans. In this review, we focus on what is known about genome organization in Drosophila and discuss it from TADs to territory. We start by examining intrachromosomal organization at the sub-chromosome level into TADs, followed by a comprehensive analysis of the known proteins that play a key role in TAD formation and boundary establishment. We then zoom out to examine interhomolog interactions such as pairing and transvection that are abundant in Drosophila but rare in other model systems. Finally, we discuss chromosome territories that form within the nucleus, resulting in a complete picture of the multi-scale organization of the Drosophila genome.

Polymer modelling unveils the roles of heterochromatin and nucleolar organizing regions in shaping 3D genome organization in Arabidopsis thaliana

The 3D genome is characterized by a complex organization made of genomic and epigenomic layers with profound implications on gene regulation and cell function. However, the understanding of the fundamental mechanisms driving the crosstalk between nuclear architecture and (epi)genomic information is still lacking. The plant Arabidopsis thaliana is a powerful model organism to address these questions owing to its compact genome for which we have a rich collection of microscopy, chromosome conformation capture (Hi-C) and ChIP-seq experiments. Using polymer modelling, we investigate the roles of nucleolus formation and epigenomics-driven interactions in shaping the 3D genome of A. thaliana. By validation of several predictions with published data, we demonstrate that self-attracting nucleolar organizing regions and repulsive constitutive heterochromatin are major mechanisms to regulate the organization of chromosomes. Simulations also suggest that interphase chromosomes maintain a partial structural memory of the V-shapes, typical of (sub)metacentric chromosomes in anaphase. Additionally, self-attraction between facultative heterochromatin regions facilitates the formation of Polycomb bodies hosting H3K27me3-enriched gene-clusters. Since nucleolus and heterochromatin are highly-conserved in eukaryotic cells, our findings pave the way for a comprehensive characterization of the generic principles that are likely to shape and regulate the 3D genome in many species.

Cellular Response to Proton Irradiation: A Simulation Study with TOPAS-nBio

The cellular response to ionizing radiation continues to be of significant research interest in cancer radiotherapy, and DNA is recognized as the critical target for most of the biologic effects of radiation. Incident particles can cause initial DNA damages through physical and chemical interactions within a short time scale. Initial DNA damages can undergo repair via different pathways available at different stages of the cell cycle. The misrepair of DNA damage results in genomic rearrangement and causes mutations and chromosome aberrations, which are drivers of cell death. This work presents an integrated study of simulating cell response after proton irradiation with energies of 0.5-500 MeV (LET of 60-0.2 keV/µm). A model of a whole nucleus with fractal DNA geometry was implemented in TOPAS-nBio for initial DNA damage simulations. The default physics and chemistry models in TOPAS-nBio were used to describe interactions of primary particles, secondary particles, and radiolysis products within the nucleus. The initial DNA double-strand break (DSB) yield was found to increase from 6.5 DSB/Gy/Gbp at low-linear energy transfer (LET) of 0.2 keV/µm to 21.2 DSB/Gy/Gbp at high LET of 60 keV/µm. A mechanistic repair model was applied to predict the characteristics of DNA damage repair and dose response of chromosome aberrations. It was found that more than 95% of the DSBs are repaired within the first 24 h and the misrepaired DSB fraction increases rapidly with LET and reaches 15.8% at 60 keV/µm with an estimated chromosome aberration detection threshold of 3 Mbp. The dicentric and acentric fragment yields and the dose response of micronuclei formation after proton irradiation were calculated and compared with experimental results.

4D nucleome modeling

The intrinsic dynamic nature of chromosomes is emerging as a fundamental component in regulating DNA transcription, replication, and damage-repair among other nuclear functions. With this increased awareness, reinforced over the last ten years, many new experimental techniques, mainly based on microscopy and chromosome conformation capture, have been introduced to study the genome in space and time. Owing to the increasing complexity of these cutting-edge techniques, computational approaches have become of paramount importance to interpret, contextualize, and complement such experiments with new insights. Hence, it is becoming crucial for experimental biologists to have a clear understanding of the diverse theoretical modeling approaches available and the biological information each of them can provide.

Fully integrated Monte Carlo simulation for evaluating radiation induced DNA damage and subsequent repair using Geant4-DNA

Ionising radiation induced DNA damage and subsequent biological responses to it depend on the radiation’s track-structure and its energy loss distribution pattern. To investigate the underlying biological mechanisms involved in such complex system, there is need of predicting biological response by integrated Monte Carlo (MC) simulations across physics, chemistry and biology. Hence, in this work, we have developed an application using the open source Geant4-DNA toolkit to propose a realistic “fully integrated” MC simulation to calculate both early DNA damage and subsequent biological responses with time. We had previously developed an application allowing simulations of radiation induced early DNA damage on a naked cell nucleus model. In the new version presented in this work, we have developed three additional important features: (1) modeling of a realistic cell geometry, (2) inclusion of a biological repair model, (3) refinement of DNA damage parameters for direct damage and indirect damage scoring. The simulation results are validated with experimental data in terms of Single Strand Break (SSB) yields for plasmid and Double Strand Break (DSB) yields for plasmid/human cell. In addition, the yields of indirect DSBs are compatible with the experimental scavengeable damage fraction. The simulation application also demonstrates agreement with experimental data of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ-H2AX yields for gamma ray irradiation. Using this application, it is now possible to predict biological response along time through track-structure MC simulations.

Chromatin as an active polymeric material

The patterns of the large-scale spatial organization of chromatin in interphase human somatic cells are not random. Such patterns include the radial separation of euchromatin and heterochromatin, the territorial organization of individual chromosomes, the non-random locations of chromosome territories and the differential positioning of the two X chromosomes in female cells. These features of large-scale nuclear architecture follow naturally from the hypothesis that ATP-consuming non-equilibrium processes associated with highly transcribed regions of chromosomes are a source of 'active' forces. These forces are in excess of those that arise from Brownian motion. Simulations of model chromosomes that incorporate such activity recapitulate these features. In addition, they reproduce many other aspects of the spatial organization of chromatin at large scales that are known from experiments. Our results, reviewed here, suggest that the distribution of transcriptional activity across chromosomes underlies many aspects of large-scale nuclear architecture that were hitherto believed to be unrelated.

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.

A parameter sensitivity study for simulating DNA damage after proton irradiation using TOPAS-nBio

Monte Carlo (MC) track structure simulation tools are commonly used for predicting radiation induced DNA damage by modeling the physical and chemical reactions at the nanometer scale. However, the outcome of these MC simulations is particularly sensitive to the adopted parameters which vary significantly across studies. In this study, a previously developed full model of nuclear DNA was used to describe the DNA geometry. The TOPAS-nBio MC toolkit was used to investigate the impact of physics and chemistry models as well as three key parameters (the energy threshold for direct damage, the chemical stage time length, and the probability of damage between hydroxyl radical reactions with DNA) on the induction of DNA damage. Our results show that the difference in physics and chemistry models alone can cause differences up to 34% and 16% in the DNA double strand break (DSB) yield, respectively. Additionally, changing the direct damage threshold, chemical stage length, and hydroxyl damage probability can cause differences of up to 28%, 51%, and 71% in predicted DSB yields, respectively, for the configurations in this study.

Hi-D: nanoscale mapping of nuclear dynamics in single living cells

Bulk chromatin motion has not been analyzed at high resolution. We present Hi-D, a method to quantitatively map dynamics of chromatin and abundant nuclear proteins for every pixel simultaneously over the entire nucleus from fluorescence image series. Hi-D combines reconstruction of chromatin motion and classification of local diffusion processes by Bayesian inference. We show that DNA dynamics in the nuclear interior are spatially partitioned into 0.3-3-μm domains in a mosaic-like pattern, uncoupled from chromatin compaction. This pattern was remodeled in response to transcriptional activity. Hi-D can be applied to any dense and bulk structures opening new perspectives towards understanding motion of nuclear molecules.

4D Genome Rewiring during Oncogene-Induced and Replicative Senescence

To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation.

•

Deep multi-omics characterization of replicative and oncogene-induced senescence

•

Senescence-associated heterochromatin domains (SAHDs) form SAHFs via 3D changes

•

DNMT1 is required for SAHF formation via regulation of HMGA2 expression

•

SAHF formation leads to expression of SAHF-adjacent genes via 3D chromatin contacts

Genome organization via loop extrusion, insights from polymer physics models

Understanding how genomes fold and organize is one of the main challenges in modern biology. Recent high-throughput techniques like Hi-C, in combination with cutting-edge polymer physics models, have provided access to precise information on 3D chromosome folding to decipher the mechanisms driving such multi-scale organization. In particular, structural maintenance of chromosome (SMC) proteins play an important role in the local structuration of chromatin, putatively via a loop extrusion process. Here, we review the different polymer physics models that investigate the role of SMCs in the formation of topologically associated domains (TADs) during interphase via the formation of dynamic loops. We describe the main physical ingredients, compare them and discuss their relevance against experimental observations.

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.

Polymer models for the mechanisms of chromatin 3D folding: review and perspective

In this perspective paper, classical physical models for mammalian interphase chromatin folding are reviewed.

Global chromatin conformation differences in the Drosophila dosage compensated chromosome X

In Drosophila melanogaster the single male chromosome X undergoes an average twofold transcriptional upregulation for balancing the transcriptional output between sexes. Previous literature hypothesised that a global change in chromosome structure may accompany this process. However, recent studies based on Hi-C failed to detect these differences. Here we show that global conformational differences are specifically present in the male chromosome X and detectable using Hi-C data on sex-sorted embryos, as well as male and female cell lines, by leveraging custom data analysis solutions. We find the male chromosome X has more mid-/long-range interactions. We also identify differences at structural domain boundaries containing BEAF-32 in conjunction with CP190 or Chromator. Weakening of these domain boundaries in male chromosome X co-localizes with the binding of the dosage compensation complex and its co-factor CLAMP, reported to enhance chromatin accessibility. Together, our data strongly indicate that chromosome X dosage compensation affects global chromosome structure.

In Drosophila, dosage compensation involves a twofold transcriptional upregulation of the single male chromosome X. Here the authors show that global conformational differences are specifically present in the male X chromosome and detectable using Hi-C data, indicating that dosage compensation affects global chromosome structure.

Nonequilibrium Biophysical Processes Influence the Large-Scale Architecture of the Cell Nucleus

Model approaches to nuclear architecture have traditionally ignored the biophysical consequences of ATP-fueled active processes acting on chromatin. However, transcription-coupled activity is a source of stochastic forces that are substantially larger than the Brownian forces present at physiological temperatures. Here, we describe an approach to large-scale nuclear architecture in metazoans that incorporates cell-type-specific active processes. The model predicts the statistics of positional distributions, shapes, and overlaps of each chromosome. Simulations of the model reproduce common organizing principles underlying large-scale nuclear architecture across human cell nuclei in interphase. These include the differential positioning of euchromatin and heterochromatin, the territorial organization of chromosomes (including both gene-density-based and size-based chromosome radial positioning schemes), the nonrandom locations of chromosome territories, and the shape statistics of individual chromosomes. We propose that the biophysical consequences of the distribution of transcriptional activity across chromosomes should be central to any chromosome positioning code.

Geometrical Heterogeneity Dominates Thermal Fluctuations in Facilitating Chromatin Contacts

Within a living cell, the myriad of proteins that bind DNA introduce heterogeneously spaced kinks into an otherwise semiflexible DNA double helix. To investigate the effects of heterogeneous nucleosome binding on chromatin organization, we extend the wormlike chain model to include statistically spaced, rigid kinks. On timescales where nucleosome positions are fixed, we find that the probability of chromatin loop formation can vary by up to six orders of magnitude between two sets of nucleosome positions drawn from the same distribution. On longer timescales, we show that continuous rerandomization due to nucleosome turnover results in chromatin tracing out an effective WLC with a dramatically smaller Kuhn length than bare DNA. Together, these observations demonstrate that nucleosome spacing acts as the primary source of the structural heterogeneity that dominates local and global chromatin organization.

Merging 1D and 3D genomic information: Challenges in modelling and validation

Genome organization in eukaryotes during interphase stems from the delicate balance between non-random correlations present in the DNA polynucleotide linear sequence and the physico/chemical reactions which shape continuously the form and structure of DNA and chromatin inside the nucleus of the cell. It is now clear that these mechanisms have a key role in important processes like gene regulation, yet the detailed ways they act simultaneously and, eventually, come to influence each other even across very different length-scales remain largely unexplored. In this paper, we recapitulate some of the main results concerning gene regulatory and physical mechanisms, in relation to the information encoded in the 1D sequence and the 3D folding structure of DNA. In particular, we stress how reciprocal crossfeeding between 1D and 3D models may provide original insight into how these complex processes work and influence each other. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Functional Consequences of Chromosomal Rearrangements on Gene Expression: Not So Deleterious After All?

Chromosomes are folded and organized into topologically associating domains (TADs) which provide a framework for the interaction of enhancers with the promoter of their target gene(s). Structural rearrangements observed during evolution or in disease contexts suggest that changes in genome organization strongly affect gene expression and can have drastic phenotypic effects. In this review, I will discuss how recent genomic engineering experiments reveal a more contrasted picture, suggesting that TADs are important but not always essential for gene expression regulation.

Microrheology of interphase chromosomes with spatial constraints: a computational study

Chromatin fibers within the interior of the nucleus of the cell make stable interactions with the nucleoskeleton, an ensemble of 'extra-chromatin' structures which help ensuring genome stability. Although the role of these interactions appears crucial to the correct behavior of the cell, their impact on chromatin structure and dynamics remains to be elucidated. In order to tackle this important issue, in this work we introduce a simple polymer model for chromatin fibers in interphase which takes into account the two generic properties of chain-versus-chain mutual uncrossability and the presence of stable binding interactions to an extra-chromatin nuclear matrix. To study how these constraints affect chromatin structure from small to large scales, we employ extensive molecular dynamics computer simulations and we monitor the motion of nanoprobes of different sizes embedded within the polymer medium. Our results demonstrate that nanoprobes show hampered motion whenever their linear size becomes larger than chromatin stiffness. This transition is also displaying features which usually belong to the realm of glassy systems, namely long-tail correlations in the distribution functions of nanoprobe spatial displacements and heterogeneous behavior accompanied by ergodicity breaking.

Large-scale chromatin organisation in interphase, mitosis and meiosis

The spatial configuration of chromatin is fundamental to ensure any given cell can fulfil its functional duties, from gene expression to specialised cellular division. Significant technological innovations have facilitated further insights into the structure, function and regulation of three-dimensional chromatin organisation. To date, the vast majority of investigations into chromatin organisation have been conducted in interphase and mitotic cells leaving meiotic chromatin relatively unexplored. In combination, cytological and genome-wide contact frequency analyses in mammalian germ cells have recently demonstrated that large-scale chromatin structures in meiotic prophase I are reminiscent of the sequential loop arrays found in mitotic cells, although interphase-like segmentation of transcriptionally active and inactive regions are also evident along the length of chromosomes. Here, we discuss the similarities and differences of such large-scale chromatin architecture, between interphase, mitotic and meiotic cells, as well as their functional relevance and the proposed modulatory mechanisms which underlie them.