How enzymatic activity is involved in chromatin organization

  1. Rakesh Das  Is a corresponding author
  2. Takahiro Sakaue
  3. GV Shivashankar
  4. Jacques Prost
  5. Tetsuya Hiraiwa  Is a corresponding author
  1. Mechanobiology Institute, National University of Singapore, Singapore
  2. Department of Physics and Mathematics, Aoyama Gakuin University, Japan
  3. ETH Zurich, Switzerland
  4. Paul Scherrer Institute, Switzerland
  5. Laboratoire Physico Chimie Curie, Institut Curie, Paris Science et Lettres Research University, France
5 figures, 1 table and 3 additional files

Figures

Microphase separation of eu- and heterochromatic regions due to enzymatic activity.

(a) A random multiblock copolymer comprising A and B beads connected by springs confined within a spherical cavity. All the data are shown for block size b=4. (b) Topo-II enzyme catches two A’s in …

Phase separation in system comprising self-avoiding and phantom regions.

(a) An equilibrium copolymer system comprising phantom (A′) and self-avoiding (B) beads is simulated in the absence of steric interaction between A′’s. The system shows microphase separation. A …

Figure 3 with 6 supplements
Characteristics of Topo-II-induced microphase separation configurations.

(a) Density distribution of A and B beads in radial direction, plotted for fixed ϵHC=4. (b) Sample snapshots (hemisphere cuts) showing wall-like organization of A’s for Λ>0. (c) Local nematic order …

Figure 3—figure supplement 1
AA bonds along the polymer are organized in planes of the walls formed by Topo-II.

We count the average number of AA and BB bonds (springs) along the polymer in local grids (×× in size) and calculate the amplitude of the nematic order parameter of those bonds. Negative, zero, and …

Figure 3—figure supplement 2
Topo-II breaks local isotropy of the AA and BB bonds.

Phase diagrams indicating the mean local nematic order parameter of the AA and BB bonds (i.e., SAA and SBB) are shown on the Λ-ϵHC plane for different composition of the system. SAA<0 and SBB>0 for Λ>0 while …

Figure 3—figure supplement 3
Comparison of monodisperse differential active model (MdDAM) with other transient models.

(a) Model schematic of Topo-II’s activity, and a sample snapshot is shown for Λ=λra(1/λan+1/λnr)=0.0309 in the absence of heterochromatin affinity. (b–c) Model schematics and sample snapshots are shown for two-state …

Figure 3—figure supplement 4
A non-transient effective attraction model comparable to monodisperse differential active model (MdDAM) does not show phase separation.

To construct the effective model, we consider an additional potential HAA=i,jA;rijhAA=i,jA;rijϵAAexp(αvexrij2) in Equation 2, and maintain all hvex’s in the repulsion state. MdDAM shows strong phase separation for the case with ϕA=0.5, ϵHC=0, …

Figure 3—figure supplement 5
Effect of the geometry of the cavity on phase separation configurations.

We consider a cubic cavity of linear dimension 10 with closed boundary on z-direction. We assume periodic boundary condition for other two directions. We simulate monodisperse differential active …

Figure 3—figure supplement 6
Features of heterochromatic foci shown for ϵHC=4.

Radius of gyration of segmented foci are plotted against their mean radial position rfoci,i=(BirB)/nB,i, nB,i number of B’s forming foci i. A few sample foci have been shown on the respective panels. The images of …

Microphase separation in bidisperse model, motivated from super-resolution microscopy data.

(a) Extracted data for mean cluster sizes and radial distribution functions of histone marks characteristic to eu- and heterochromatic regions. The data were extracted from Xu et al., 2018. …

Author response image 1

Tables

Table 1
Choice of model parameters.
PotentialParameters
Hspringk=22e-2; r0=iconnectedbeadsdi.
HvexMonodisperse model: ϵvex=8e; αvex=7.9585-2.
Bidisperse model for AA pairs: ϵvex=6.6539e; αvex=9.5686-2.
Bidisperse model for BB pairs: ϵvex=11.507e; αvex=5.5330-2.
Bidisperse model for AB pairs: ϵvex=8.9153e; αvex=7.1414-2.
HHCαHC=100-2.
Hconfinementp=4e; s=2; Rs=0.65Rg; κ=1/Rg; Rg=di/2 where i A, B

Additional files

Source code 1

CPU-based FORTRAN simulation code using OpenMP API.

Instructions to use this can be found in the README text accompanying the source code.

https://cdn.elifesciences.org/articles/79901/elife-79901-code1-v2.zip
Source code 2

CUDA FORTRAN simulation code using GPU acceleration.

Instructions to use this can be found in the README text accompanying the source code.

https://cdn.elifesciences.org/articles/79901/elife-79901-code2-v2.zip
Transparent reporting form
https://cdn.elifesciences.org/articles/79901/elife-79901-transrepform1-v2.pdf

Download links