A new application of the phase-field method for understanding the mechanisms of nuclear architecture reorganization

Consistencies and contradictions in different polymer models of chromatin architecture

Genetic information is stored in very long DNA molecules, which are folded to form chromatin, a similarly long polymer fibre that is ultimately organised into chromosomes. The organisation of chromatin is fundamental to many cellular functions, from the expression of the genetic information to cell division. As a long polymer, chromatin is very flexible and may adopt a myriad of shapes. Globally, the polymer physics governing chromatin dynamics is very well understood. But chromatin is not uniform and regions of it, with chemical modifications and bound effectors, form domains and compartments through mechanisms not yet clear. Polymer models have been successfully used to investigate these mechanisms to explain cytological observations and build hypothesis for experimental validation. Many different approaches to conceptualise chromatin in polymer models can be envisioned and each reflects different aspects. Here, we compare recent approaches that aim at reproducing prominent features of interphase chromatin organisation: the compartmentalisation into eu- and heterochromatin compartments, the formation of a nucleolus, chromatin loops and the rosette and Rabl conformations of interphase chromosomes. We highlight commonalities and contradictions that point to a modulation of the mechanisms involved to fine degree. Consolidating models will require the inclusion of yet hidden or neglected parameters.

Epigenetic-structural changes in X chromosomes promote Xic pairing during early differentiation of mouse embryonic stem cells

X chromosome inactivation center (Xic) pairing occurs during the differentiation of embryonic stem (ES) cells from female mouse embryos, and is related to X chromosome inactivation, the circadian clock, intra-nucleus architecture, and metabolism. However, the mechanisms underlying the identification and approach of X chromosome pairs in the crowded nucleus are unclear. To elucidate the driving force of Xic pairing, we developed a coarse-grained molecular dynamics model of intranuclear chromosomes in ES cells and in cells 2 days after the onset of differentiation (2-day cells) by considering intrachromosomal epigenetic-structural feature-dependent mechanics. The analysis of the experimental data showed that X-chromosomes exhibit the rearrangement of their distributions of open/closed chromatin regions on their surfaces during cell differentiation. By simulating models where the excluded volume effects of closed chromatin regions are stronger than those of open chromatin regions, such rearrangement of open/closed chromatin regions on X-chromosome surfaces promoted the mutual approach of the Xic pair. These findings suggested that local intrachromosomal epigenetic features may contribute to the regulation of cell species-dependent differences in intranuclear architecture.

The extra-embryonic space and the local contour are crucial geometric constraints regulating cell arrangement

In multicellular systems, cells communicate with adjacent cells to determine their positions and fates, an arrangement important for cellular development. Orientation of cell division, cell-cell interactions (i.e. attraction and repulsion) and geometric constraints are three major factors that define cell arrangement. In particular, geometric constraints are difficult to reveal in experiments, and the contribution of the local contour of the boundary has remained elusive. In this study, we developed a multicellular morphology model based on the phase-field method so that precise geometric constraints can be incorporated. Our application of the model to nematode embryos predicted that the amount of extra-embryonic space, the empty space within the eggshell that is not occupied by embryonic cells, affects cell arrangement in a manner dependent on the local contour and other factors. The prediction was validated experimentally by increasing the extra-embryonic space in the Caenorhabditis elegans embryo. Overall, our analyses characterized the roles of geometrical contributors, specifically the amount of extra-embryonic space and the local contour, on cell arrangements. These factors should be considered for multicellular systems.

Summary: The local contour and the extra-embryonic space, the empty space within the eggshell not occupied by embryonic cells, are important geometric constraints in cell arrangement of nematode embryos.

Modeling and simulation of cell nuclear architecture reorganization process

We develop a special phase field/diffusive interface method to model the nuclear architecture reorganization process. In particular, we use a Lagrange multiplier approach in the phase field model to preserve the specific physical and geometrical constraints for the biological events. We develop several efficient and robust linear and weakly nonlinear schemes for this new model. To validate the model and numerical methods, we present ample numerical simulations which in particular reproduce several processes of nuclear architecture reorganization from the experiment literature.

The Role of Cytoplasmic MEX-5/6 Polarity in Asymmetric Cell Division

In the process of asymmetric cell division, the mother cell induces polarity in both the membrane and the cytosol by distributing substrates and components asymmetrically. Such polarity formation results from the harmonization of the upstream and downstream polarities between the cell membrane and the cytosol. MEX-5/6 is a well-investigated downstream cytoplasmic protein, which is deeply involved in the membrane polarity of the upstream transmembrane protein PAR in the Caenorhabditis elegans embryo. In contrast to the extensive exploration of membrane PAR polarity, cytoplasmic polarity is poorly understood, and the precise contribution of cytoplasmic polarity to the membrane PAR polarity remains largely unknown. In this study, we explored the interplay between the cytoplasmic MEX-5/6 polarity and the membrane PAR polarity by developing a mathematical model that integrates the dynamics of PAR and MEX-5/6 and reflects the cell geometry. Our investigations show that the downstream cytoplasmic protein MEX-5/6 plays an indispensable role in causing a robust upstream PAR polarity, and the integrated understanding of their interplay, including the effect of the cell geometry, is essential for the study of polarity formation in asymmetric cell division.

Exploring chromosomal structural heterogeneity across multiple cell lines

Using computer simulations, we generate cell-specific 3D chromosomal structures and compare them to recently published chromatin structures obtained through microscopy. We demonstrate using machine learning and polymer physics simulations that epigenetic information can be used to predict the structural ensembles of multiple human cell lines. Theory predicts that chromosome structures are fluid and can only be described by an ensemble, which is consistent with the observation that chromosomes exhibit no unique fold. Nevertheless, our analysis of both structures from simulation and microscopy reveals that short segments of chromatin make two-state transitions between closed conformations and open dumbbell conformations. Finally, we study the conformational changes associated with the switching of genomic compartments observed in human cell lines. The formation of genomic compartments resembles hydrophobic collapse in protein folding, with the aggregation of denser and predominantly inactive chromatin driving the positioning of active chromatin toward the surface of individual chromosomal territories.

Mesoscale Modeling of Nucleosome-Binding Antibody PL2-6: Mono- versus Bivalent Chromatin Complexes

Interactions of chromatin with bivalent immunoglobin nucleosome-binding antibodies and their monovalent (papain-derived) antigen-binding fragment analogs are useful probes for examining chromatin conformational states. To help interpret antibody-chromatin interactions and explore how antibodies might compete for interactions with chromatin components, we incorporate coarse-grained PL2-6 antibody modeling into our mesoscale chromatin model. We analyze interactions and fiber structures for the antibody-chromatin complexes in open and condensed chromatin, with and without H1 linker histone (LH). Despite minimal and transient interactions at physiological salt, we capture significant differences in antibody-chromatin complex configurations in open fibers, with more intense interactions between the bivalent antibody and chromatin compared to monovalent antigen-binding fragments. For these open chromatin fiber morphologies, antibody binding to histone tails is increased and compaction is greater for bivalent compared to monovalent and antibody-free systems. Differences between monovalent and bivalent binding result from antibody competition with internal chromatin fiber components (nucleosome core and linker DNA) for histone tail (H3, H4, H2A, H2B) interactions. This antibody competition for tail contacts reduces tail-core and tail-linker interactions and increases tail-antibody interactions. Such internal structural changes in open fibers resemble mechanisms of LH condensation, driven by charge screening and entropy changes. For condensed fibers at physiological salt, the three systems are much more similar overall, but some subtle tail interaction differences can be noted. Adding LH results in less-dramatic changes for all systems, except that the bivalent complex at physiological salt shows cooperative effects between LH and the antibodies in condensing chromatin fibers. Such dynamic interactions that depend on the internal structure and complex-stabilizing interactions within the chromatin fiber have implications for gene regulation and other chromatin complexes such as with LH, remodeling proteins, and small molecular chaperones that bind and modulate chromatin structure.

Mesoscale Liquid Model of Chromatin Recapitulates Nuclear Order of Eukaryotes

The nuclear envelope segregates the genome of Eukaryota from the cytoplasm. Within the nucleus, chromatin is further compartmentalized into architectures that change throughout the lifetime of the cell. Epigenetic patterns along the chromatin polymer strongly correlate with chromatin compartmentalization and, accordingly, also change during the cell life cycle and at differentiation. Recently, it has been suggested that subnuclear chromatin compartmentalization might result from a process of liquid-liquid phase separation orchestrated by the epigenetic marking and operated by proteins that bind to chromatin. Here, we translate these observations into a diffuse interface model of chromatin, which we named the mesoscale liquid model of nucleus. Using this streamlined continuum model of the genome, we study the large-scale rearrangements of chromatin that happen at different stages of the growth and senescence of the cell and during nuclear inversion events. In particular, we investigate the role of droplet diffusion, fluctuations, and heterochromatin-lamina interactions during nuclear remodeling. Our results indicate that the physical process of liquid-liquid phase separation, together with surface effects, is sufficient to recapitulate much of the large-scale morphology and dynamics of chromatin along the life cycle of cells.

Role of dynamic nuclear deformation on genomic architecture reorganization

Higher-order genomic architecture varies according to cell type and changes dramatically during differentiation. One of the remarkable examples of spatial genomic reorganization is the rod photoreceptor cell differentiation in nocturnal mammals. The inverted nuclear architecture found in adult mouse rod cells is formed through the reorganization of the conventional architecture during terminal differentiation. However, the mechanisms underlying these changes remain largely unknown. Here, we found that the dynamic deformation of nuclei via actomyosin-mediated contractility contributes to chromocenter clustering and promotes genomic architecture reorganization during differentiation by conducting an in cellulo experiment coupled with phase-field modeling. Similar patterns of dynamic deformation of the nucleus and a concomitant migration of the nuclear content were also observed in rod cells derived from the developing mouse retina. These results indicate that the common phenomenon of dynamic nuclear deformation, which accompanies dynamic cell behavior, can be a universal mechanism for spatiotemporal genomic reorganization.

The motion and spatial reorganization of sub-nuclear domains have been extensively studied using microscopy, but the underlying mechanisms that promote the reorganization are still poorly understood. We found that dynamic nuclear deformation provides a driving force for long-range migration and aggressive clustering of chromocenters, which ultimately induces nuclear architecture reorganization. This study significantly contributes to an improved understanding of the role of nuclear deformation and reorganization of nuclear architecture that seems complex but is based on simple physical properties of the cell.

Heterochromatin drives compartmentalization of inverted and conventional nuclei

The nucleus of mammalian cells displays a distinct spatial segregation of active euchromatic and inactive heterochromatic regions of the genome^1,2. In conventional nuclei, microscopy shows that euchromatin is localized in the nuclear interior and heterochromatin at the nuclear periphery^1,2. Genome-wide chromosome conformation capture (Hi-C) analyses show this segregation as a plaid pattern of contact enrichment within euchromatin and heterochromatin compartments³, and depletion between them. Many mechanisms for the formation of compartments have been proposed, such as attraction of heterochromatin to the nuclear lamina^2,4, preferential attraction of similar chromatin to each other^1,4-12, higher levels of chromatin mobility in active chromatin^13-15 and transcription-related clustering of euchromatin^16,17. However, these hypotheses have remained inconclusive, owing to the difficulty of disentangling intra-chromatin and chromatin-lamina interactions in conventional nuclei¹⁸. The marked reorganization of interphase chromosomes in the inverted nuclei of rods in nocturnal mammals^19,20 provides an opportunity to elucidate the mechanisms that underlie spatial compartmentalization. Here we combine Hi-C analysis of inverted rod nuclei with microscopy and polymer simulations. We find that attractions between heterochromatic regions are crucial for establishing both compartmentalization and the concentric shells of pericentromeric heterochromatin, facultative heterochromatin and euchromatin in the inverted nucleus. When interactions between heterochromatin and the lamina are added, the same model recreates the conventional nuclear organization. In addition, our models allow us to rule out mechanisms of compartmentalization that involve strong euchromatin interactions. Together, our experiments and modelling suggest that attractions between heterochromatic regions are essential for the phase separation of the active and inactive genome in inverted and conventional nuclei, whereas interactions of the chromatin with the lamina are necessary to build the conventional architecture from these segregated phases.

Physical properties of the chromosomes and implications for development

Remarkable progress has been made in understanding chromosome structures inside the cell nucleus. Recent advances in Hi-C technologies enable the detection of genome-wide chromatin interactions, providing insight into three-dimensional (3D) genome organization. Advancements in the spatial and temporal resolutions of imaging as well as in molecular biological techniques allow the tracking of specific chromosomal loci, improving our understanding of chromosome movements. From these data, we are beginning to understand how the intra-nuclear locations of chromatin loci and the 3D genome structure change during development and differentiation. This emerging field of genome structure and dynamics research requires an interdisciplinary approach including efficient collaborations between experimental biologists and physicists, informaticians, or engineers. Quantitative and mathematical analyses based on polymer physics are becoming increasingly important for processing and interpreting experimental data on 3D chromosome structures and dynamics. In this review, we aim to provide an overview of recent research on the physical aspects of chromosome structure and dynamics oriented for biologists. These studies have mainly focused on chromosomes at the cellular level, using unicellular organisms and cultured cells. However, physical parameters that change during development, such as nuclear size, may impact genome structure and dynamics. Here, we discuss how chromatin dynamics and genome structures in early embryos change during development, which we expect will be a hot topic in the field of chromatin dynamics in the near future. We hope this review helps developmental biologists to quantitatively investigate the physical natures of chromosomes in developmental biology research.