Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Chromatin network retards nucleoli coalescence

## Abstract

Nuclear bodies are membraneless condensates that may form via liquid-liquid phase separation. The viscoelastic chromatin network could impact their stability and may hold the key for understanding experimental observations that defy predictions of classical theories. However, quantitative studies on the role of the chromatin network in phase separation have remained challenging. Using a diploid human genome model parameterized with chromosome conformation capture (Hi-C) data, we study the thermodynamics and kinetics of nucleoli formation. Dynamical simulations predict the formation of multiple droplets for nucleolar particles that experience specific interactions with nucleolus-associated domains (NADs). Coarsening dynamics, surface tension, and coalescence kinetics of the simulated droplets are all in quantitative agreement with experimental measurements for nucleoli. Free energy calculations further support that a two-droplet state, often observed for nucleoli in somatic cells, is metastable and separated from the single-droplet state with an entropic barrier. Our study suggests that nucleoli-chromatin interactions facilitate droplets’ nucleation but hinder their coarsening due to the coupled motion between droplets and the chromatin network: as droplets coalesce, the chromatin network becomes increasingly constrained. Therefore, the chromatin network supports a nucleation and arrest mechanism to stabilize the multi-droplet state for nucleoli and possibly for other nuclear bodies.

## Introduction

Nuclear bodies are pervasive in eukaryotic cells and play a diverse set of functions1, including RNA metabolism, transcriptional regulation2, genome organization3, etc. They are membraneless structures that mainly consist of protein and RNA molecules4. Their lack of a lipid-rich barrier allows rapid exchange of components with the nucleoplasm in responses to environmental cues and stress signaling5. Nuclear bodies also effectively increase the local concentration of enzymes involved in particular functions to facilitate more efficient cellular reactions6.

Increasing evidence supports that nuclear bodies function as biomolecular condensates formed via liquid-liquid phase separation (LLPS)4,7,8,9. They exhibit round morphologies and dynamic fluid properties10,11. Two nuclear bodies can fuse into larger condensates following growth kinetics with similar scaling behavior as that observed for simple liquids12,13,14. In addition, their assembly was shown to be concentration-dependent, and the coarsening and growth dynamics can be quantitatively modeled with classical theories of phase separation15. At the molecular level, detailed mechanistic models for LLPS are beginning to emerge as well. In particular, low complexity domains and intrinsically disordered regions are enriched in many of the proteins associated with nuclear bodies16,17,18,19,20. These features enable non-specific, multivalent interactions that drive the formation of dynamical condensates.

However, several observations of nuclear bodies appear to defy predictions from classical nucleation and phase separation theories. In particular, these theories predict the thermodynamic equilibrium to consist of a single condensate that minimizes the surface energy21,22. On the other hand, multiple nucleoli (~2–5) can stably coexist in the same nucleus11,12,23,24,25,26, as can paraspeckles27 and nuclear speckles28. The exact number of nuclear bodies is sensitive to various factors, including cell volume29 and nuclear lamina composition30,31. It has been proposed that non-equilibrium activities can dynamically alter protein-protein interactions to stabilize the multi-droplet state15,32,33,34,35. In addition, the chromatin network may suppress droplet coarsening through mechanical frustration as well36,37,38.

Nuclear bodies and chromatin are also known to form attractive interactions27, further complicating phase separation inside the nucleus beyond mechanical stress. For example, the upstream binding factor (UBF), which is a DNA binding protein and a key component of nucleoli, is known to recognize ribosomal DNA (rDNA) repeats to seed the rapid formation of nucleoli after cell division39. Correspondingly, rDNA and other chromosome segments, which are collectively noted as nucleolus-association domains (NADs)9,40,41, can be seen inside and adjacent to nucleoli42. Paraspeckles27 and speckles28,43 have been found in spatial proximity with chromatin as well. In addition to proteins, nuclear bodies can harbor non-coding RNA to contact chromatin either by recruiting intermediate protein molecules or by forming RNA-DNA duplex or triplex44. Since chromatin forms a viscoelastic network spanning the nucleus13,45,46,47,48,49, its interactions with nuclear bodies could impact the thermodynamics and kinetics of phase separation.

We carry out molecular dynamics simulations to investigate nucleoli formation with a computational model that explicitly considers nucleolar particles, the chromatin network, and the interactions between the two. We represent the chromatin network using a diploid human genome model that provides explicit polymer configurations for individual chromosomes. Interactions within and among chromosomes are optimized based on chromosome conformation capture (Hi-C) experiments to ensure in vivo relevance. The simulated dynamical process of droplet growth, coarsening, and coalescence are in quantitative agreement with experimental measurements. Importantly, our simulations predict the formation of multiple droplets, much like the coexistence of several nucleoli seen in the nucleus. We show that a two-droplet state is metastable and separated from the single-droplet state with an entropic barrier. The barrier arises from the chromatin network, which becomes more constrained upon droplet coalescence. Nucleolar particle-chromatin interactions link the motion between the chromatin network and the droplets, and stronger interactions are shown to produce more droplets. Our study provides insight into the critical role of the chromatin network on the formation of nucleoli and nuclear bodies in general.

## Results

### Phase separation with chromatin network leads to multiple droplets

Leveraging a recently introduced computational model for the diploid human genome50,51, we studied the impact of the chromatin network on nucleoli formation. We modeled the genome at the one megabase (Mb) resolution as a collection of 46 chromosomes inside spherical confinement (Fig. 1). Each chromosome is represented as a string of beads, which can be assigned as one of three types, A, B, or C. A and B correspond to the two compartment types that contribute to the checkerboard patterns typically seen in Hi-C contact maps52,53, and C marks centromeric regions. Interactions among the beads were optimized to reproduce various average contact probabilities determined from Hi-C experiments for GM12878 cells using the maximum entropy optimization algorithm54,55 (see Supplementary Material). Because of the non-equilibrium nature of the system, these experimentally derived interactions and temperature represent effective approximations to the steady-state distribution of genome organization (see “Methods”). While the model does not explicitly include histone modifications, transcription factors, or molecular motors, reproducing the contact probabilities between genomic segments measured in situ effectively allows it to account for the contribution of these factors to genome organization.

We introduced additional coarse-grained particles to model molecules that make up the nucleoli. The size and number of these particles were chosen based on the concentration of nucleolar proteins and the volume fraction of nucleoli (see “Methods”). Nucleolar particles share favorable interactions with each other and with nucleolus-associated domains (NADs), which are chromatin regions strongly associated with nucleoli41. A non-specific interaction term was also introduced between nucleolar particles and non-NAD chromatin regions. While our results are relatively robust with respect to the strength of these interactions (Supplementary Fig. 1), we found that, with the chosen values, the surface tension of simulated droplets compares favorably to experimental values for nucleoli (see Supplementary Material).

Starting from an equilibrated genome structure and randomly distributed nucleolar particles (see “Methods”), we carried out a total of twelve independent molecular dynamics simulations. These simulations lasted for 20 million steps, much longer than the relaxation timescale of chromosome conformations (Supplementary Fig. 2). In all but one case, the nucleolar particles aggregated into multiple droplets that persisted to the end of the simulations (Fig. 2b). This result contrasts with simulations performed without the chromatin network, where nucleolar particles always condense into a single droplet (Supplementary Fig. 3). Notably, the emergence of the multi-droplet state is insensitive to the configuration used to initialize the simulations (Supplementary Fig. 4), the interactions between chromosomes (Supplementary Figs. 5,6), and the resolution of the genome model (Supplementary Fig. 7). We found that nucleolar particles forming the droplets undergo dynamical exchange with the surrounding nucleoplasm while maintaining the droplet size (~1 μm) on timescales of several tens of minutes (Supplementary Fig. 8). Dynamical exchange of materials has been observed in fluorescence recovery after photobleaching (FRAP) experiments56,57 and directly supports the liquid-like property of the droplets. The droplets are preferentially localized at the interior of the nucleus (Supplementary Fig. 9A). Because of their close association with these droplets, NADs are closer to the nuclear center than other heterochromatin regions as well58. However, not all NADs bind to the droplets, and a significant fraction of them localize towards the nuclear envelope (Supplementary Fig. 9B). Two classes of NADs that vary in nuclear localization have indeed been observed in prior studies59.

Analyzing the simulated genome structures, we found that the model with nucleolar particles reproduces global features of the genome organization, including the formation of chromosome territories60, the compartmentalization of heterochromatin/euchromatin61, and the clustering of centromeric regions62, as shown in Fig. 2c, d. The simulated chromosome radial positions agree well with experimental values63, and the Pearson’s correlation coefficient between the two is 0.89 (Fig. 2e). Because of its coarse resolution, the model will inevitably miss certain features of genome organizations, including the formation of chromatin loops and topologically associating domains53. While these structural motifs at fine scales are crucial for an accurate representation of the genome organization, they are less significant for exploring the mechanisms of nuclear body formation, at least at a qualitative level (see Supplementary Fig. 7).

### An entropic barrier hinders droplet coalescence

The emergence of a long-lived state with multiple droplets, while in contrast with predictions of the classical nucleation theory21, is consistent with the coexistence of multiple nucleoli in the nucleus11,13,24,25,26. To better understand the stability of the multi-droplet state, we computed the free energy profile as a function of the radius of gyration (Rg) for a two-droplet system (see “Methods”). Rg is defined as $$\sqrt{\frac{1}{N}\mathop{\sum }\nolimits_{i = 1}^{N}{\left|{{{{{{{{\bf{r}}}}}}}}}_{i}-{{{{{{{{\bf{r}}}}}}}}}_{{{\mbox{com}}}}\right|}^{2}}$$, where ri is the coordinate of the ith nucleolar particle and rcom corresponds to the center of mass. The summation includes all N nucleolar particles in either one of the two droplets. As the size of individual droplets remains stable, changes in Rg will be mainly driven by variations in the distance of the two droplets. However, different from a simple center-of-mass distance, which becomes ill-defined if the lists of nucleolar particles participating in droplet formation are not updated on the fly, Rg is relatively invariant with respect to the flux of particles between the two droplets (Supplementary Fig. 10). Simulations of the full system with both chromosomes and nucleolar particles were used to compute the free energy profile. Therefore, the impact of the chromatin network was accounted for implicitly even though it was not included in the definition of Rg. Umbrella sampling and temperature replica exchange were used to enhance conformational exploration.

As shown in Fig. 3a, the free energy profile exhibits two basins. While the basin at Rg ≈ 1.5 μm corresponds to the two droplet state observed before, an additional minimum with all nucleolar particles participating in a single droplet appears at smaller values (≈0.75 μm). The two basins are separated from each other with a transition state at Rg ≈ 1.13 μm. Representative configurations at the transition state show a dumbbell shape with the establishment of a thin bridge between the two droplets. Consistent with predictions of the classical nucleation theory, the one droplet state remains as the global minimum. However, the merging of the droplets is kinetically constrained due to the presence of a barrier that is ~7 kBT in height. The barrier height is much larger than the error bars (0.5 kBT) estimated via block averaging (see “Methods”), supporting its statistical significance.

To reveal the nature of the barrier, we decomposed the free energy into entropic and energetic contributions. Using the free energy profiles at different temperatures (Supplementary Fig. 11), we computed the entropy change along the collective variable with a finite difference method64,65. As shown in Fig. 3b, c, contributions from the entropy (–TΔS) continue to increase as Rg decreases, and the droplets coalesce. While restricting the motion of the two droplets to smaller distances is naturally unfavorable, the entropic penalty is intensified here due to the increasingly restricted motion of chromosomes as well. As mentioned before, nucleolar particles form extensive contacts with NADs via specific interactions, and such contacts enforce correlative motions between droplets and chromosomes. The potential energy, on the other hand, favors droplet merging and decreases continuously along the collective variable due to the increase in nucleolar particle-nucleolar particle and nucleolar particle-chromatin contacts. Therefore, the transition state arises due to the presence of the chromatin network, and dissolving the polymeric topology of chromosomes indeed removes the barrier (Supplementary Fig. 3C).

While Rg suffices for monitoring the progression of droplet coalescence, whether it serves as a “good” reaction coordinate or not requires further investigation. In particular, a good reaction coordinate should provide insight into the bottleneck that limits the reaction. Furthermore, trajectories initialized from the identified transition state will have an equal probability of committing to the reactant or product, or the so-called committor probability adopts a value of 0.566,67. Otherwise, the transition state may bear little relevance to the reaction, and mechanisms derived from it can be misleading. To evaluate the significance of the transition state for droplet coalescence, we carried out additional simulations starting from random configurations with Rg values at 1.13 μm. For each configuration, we initialized ten independent 200,000-step-long simulations with randomized velocities. We then counted the number of simulations that end up with the two-droplet state versus the single-droplet state. As shown in Fig. 4, among all these simulations, 56% of them led to the single-droplet state while the rest 44% ended up in the two-droplet state. These results strongly support the usefulness of Rg for studying coalescence and the relevance of the identified transition state for mechanistic interpretation. Since the chromatin network was not included for defining the reaction coordinate and can vary significantly at a given value for Rg (Supplementary Fig. 12), it may play a secondary or passive role in mediating coalescence.

We further isolated trajectories initialized from the transition state that led to the single-droplet state and computed the evolution of the neck radius as a function of time. As shown in Supplementary Fig. 13, the neck region was identified as the minimum of the density profile of nucleolar particles along the principal axis with the largest variance. As shown in Fig. 4b, by plotting the normalized neck radius (2R(t)/R0) with respect to the time, we obtained a power-law relationship with exponent 0.51. R0 is the average radius of the droplets before fusion. This exponent agrees with the experimental value determined for nucleoli13 and suggests that droplet coalescence proceeds in the low Reynolds number regime dominated by viscous effects from the outer fluid, i.e., the nucleoplasm68,69.

### Chromatin network gives rise to slow coarsening dynamics

The thermodynamic analysis presented in Fig. 3 suggests that the chromatin network acts much like entropic springs to hinder the coalescence of droplets70. As the droplets move close to each other, they pull on the network and restrict the accessible polymer configurations. Restoring forces from the network to maximize configurational entropy counter droplet merging and give rise to the barrier. A similar mechanism could potentially impact the coarsening dynamics and the pathway leading to the formation of multiple droplets. To better understand the role of nucleolar particle-chromatin interactions in the overall phase separation process, we analyzed the dynamical trajectories at the onset of cluster formation.

We first monitored the time evolution of the number of clusters formed along the dynamical trajectories shown in Fig. 2. The clusters were identified as high-density regions of nucleolar particles across the entire nucleus using the DBSCAN (Density-Based Spatial Clustering of Applications with Noise) algorithm71 (see Supplementary Materials). A typical trajectory is shown in Fig. 5a, and starts with zero clusters due to the random distribution of nucleolar particles in the initial configuration. The sudden appearance of nine clusters at time ~ 14 min suggests that nucleation can occur at multiple sites almost simultaneously. As time proceeds, the droplets began to merge or evaporate, and the trajectory eventually stabilizes to the two-droplet state.

We further recorded the average size of all clusters measured by their mean Rg. The increase of Rg as a function of the time follows a power-law scaling with Rgtβ. The exponent β = 0.1 differs from the phase separation process without the chromatin network. As shown in Supplementary Fig. 14, for simulations performed without the network, the initial cluster number is smaller and the exponent for cluster size growth is close to the theoretical value of 1/372,73.

To provide insight into the appearance of the abnormal exponent β, we characterized the switching events that led to the decrease in cluster numbers. Specifically, we tracked the average pair-wise distance $$\overline{{r}_{ij}(t)}$$ of nucleolar particles within each cluster. In some cases, $$\overline{{r}_{ij}(t)}$$ for all clusters remain relatively constant and the switching occurs through the Brownian motion-induced coalescence (BMC) pathway (Fig. 5b1). In other cases, we observe a significant increase in $$\overline{{r}_{ij}(t)}$$ before the switching occurs (Fig. 5b2). This increase indicates an “evaporation” of the cluster following the diffusion-limited Ostwald ripening (DOR) path. We then identified all the switching events (42 in total) from the dynamical simulations and found that ~76% of them follow the BMC path. The rest of the switching events proceed via the DOR path and often involve smaller clusters to reduce the penalty associated with cluster disassembly (Supplementary Fig. 15).

The dominance of the BMC pathway explains the dramatic slowdown of the coarsening dynamics. In particular, the scaling exponent of 1/3 was predicted based on a normal diffusion model in which distances between droplets scale linearly in time, i.e., x2(t) Dt74. However, as shown in Supplementary Fig. 16, most of the clusters exhibit sub-diffusion and x2(t) Dt1/2. Assuming that the average size of droplets r(t) is proportional to their mean distance, we arrive at r(t) t1/6 for the observed abnormal diffusion. The exponent now is closer to the value shown in Fig. 5c. The sub-diffusive motion arises from both the elastic stress produced by the viscoelastic chromatin network and the physical tethering of droplets to the chromosomes. In addition to the sub-diffusive motion, the chromatin network could further reduce the exponent and slow down the Brownian diffusion dominated coarsening dynamics by hindering droplet coalescence through entropic barriers similar to that shown in Fig. 3.

We note that the abnormal diffusion and slower coarsening kinetics have been directly observed by Lee et al. as well when monitoring the coarsening dynamics of model condensates based on intrinsically disordered protein regions in the nucleus36. In particular, they revealed a coarsening exponent of ~0.12, which is close to the value shown in Fig. 5c. The scaling exponent for nucleolar coarsening in vivo is also in good agreement with the simulated value when considering short time kinetics before 5 min15.

### Nucleoli-chromatin interactions promote cluster nucleation

Our results indicate that, while nucleolar particle-chromatin interactions increase the overall stability of the single-droplet state, they retard the coarsening kinetics by giving rise to sub-diffusion and entropic barriers. To more directly probe their impact on droplet coalescence, we recomputed the free energy profile at stronger (ϵ = 2.0 kBT) and weaker (ϵ = 1.6 kBT) nucleolar particle-chromatin interactions. As shown in Fig. 6a, while the stability of the merged state varied significantly, the transition and the two-droplet state are much less affected. These simulations again support the entropic origin of the free energy barrier.

We performed additional long-timescale simulations following the same protocol as those shown in Fig. 2. Starting from the same initial configurations, we varied nucleolar particle-chromatin interactions from 1.4 kBT to 2.8 kBT and carried out 12 independent 20-million-step long simulations for each interaction strength. We then computed the number of droplets formed at the end of these simulations. As shown in Fig. 6b, the number of droplets increases with the interaction strength. Stronger interactions facilitate the nucleation of nucleolar clusters on the chromatin, as can be seen from the increase of cluster numbers at the onset of phase separation (Fig. 6c). Cluster coarsening, again, primarily follow the BMC pathway, though the ratio of DOR increases as well due to the instability of the nucleated clusters with a smaller number of nucleolar particles (Fig. 6d). Since merging of these clusters along the BMC pathway is hindered by entropic barriers, more nucleation naturally leads to increased droplet numbers at the end of the simulations. These results further support the role of the chromatin network in hindering droplet coalescence.

## Discussion

While nuclear bodies have been shown to exhibit liquid properties, the role of the surrounding environment on their formation and stability is less known. We modeled the dynamical process of phase separation that drives nucleoli formation in the presence of the chromatin network. The diploid genome model parameterized from Hi-C data was used to describe the interactions within and among chromosomes accurately. Simulations carried out with coarse-grained nucleolar particles succeeded in producing multiple droplets with dynamical behaviors comparable to nucleoli. We found that nucleoli-chromatin interactions facilitate the nucleation of condensates and retard their coalescence, stabilizing the multiple-droplet state.

Despite our best effort, the coarsened representations for the genome and nucleoli are bound to be approximate, and certain features of the biological system would inevitably be missed. While it may struggle at quantitative predictions, the model serves as a valuable tool for mechanistic explorations. In particular, since the nucleation and arrest mechanism arises from the generic polymer nature of the chromatin network, it should be insensitive to many of the model assumptions. Indeed, as shown in Supplementary Fig. 17, the multi-droplet state was also observed using a genome model parameterized with Hi-C data generated from clinical tissue samples with colorectal cancer51. The tumor Hi-C data and genome organization differ significantly from the one presented in the main text. The stability of the multi-droplet state is also preserved in additional perturbations that directly vary the resolution of the genome model (Supplementary Fig. 7), chromosome-chromosome interactions (Supplementary Fig. 5,6), nucleolar particle size (Supplementary Fig. 18), and nucleolar particle-chromatin interactions (Supplementary Fig. 1).

Specific interactions between nucleolar particles and chromatin are crucial for the nucleation step and the emergence of the multi-droplet state (Supplementary Fig. 1). As such interactions may be present for other nuclear bodies, a similar mechanism could contribute to their formation as well. For example, several recent studies support the close contact between a subset of chromatin and speckles with distances at ~100 nm or less using techniques based on high throughput sequencing75 or imaging76. Furthermore, these speckle-associated chromatin regions are largely conserved across cell-types77, supporting the presence of non-random mechanisms for their maintenance. Close contacts with chromatin are likely stabilized by RNA molecules that are known to localize at speckle periphery78,79. Non-coding RNAs have also been found in paraspeckles and PML bodies to mediate their interactions with chromatin27.

The mechanism uncovered here differs from the suppression of droplet growth in an elastic matrix discussed in several recent studies36,37,80,81. In contrast to our focus on the viscoelastic nature of the chromatin network and its attractive interactions with phase-separating agents, these studies were concerned with cross-linked gels and repulsive interactions. In the presence of permanent cross-linking, droplet-induced swelling of the polymer network can give rise to compressive stresses that arrest the coarsening dynamics completely. Therefore, the thermodynamic equilibrium becomes the monodisperse state with many droplets whose size is largely determined by the network mesh. Because of their favorable interactions in our model, there is substantial wetting of the chromatin network by nucleolar particles. The droplets incur minimal stress on the network (Supplementary Fig. 20), and the multi-droplet state remains metastable (Fig. 3). We note that the two mechanisms are not necessarily mutually exclusive and could contribute to the formation of different condensates inside the nucleus since the dynamical properties of chromatin are known to be timescale dependent82.

Finally, it’s worth mentioning that the nucleus is a non-equilibrium system, and active processes could contribute to the stability of the multi-droplet state as well34,35. For example, enzymes such as kinase could add post-translational modifications to condensate proteins and regulate their ability in engaging multivalent interactions. Explicitly modeling the active processes within the simulation framework outlined here may be necessary to account for the complementary mechanisms and provide a more comprehensive understanding of in vivo phase separation and nuclear body formation.

## Methods

### Details on the setup and simulations of the nucleus model

We included both the genome and nucleolar proteins to simulate phase separation inside the nucleus. Following the same setup and interactions as in ref. 50, we modeled all 46 chromosomes of the diploid human genome as beads on a string at the 1 MB resolution. In addition, a total of 500 particles were introduced to represent nucleolar proteins. See below for details on estimations for the particle size and number.

In addition to their attractive self-interactions, nucleolar particles bind with chromatin via specific and non-specific interactions. All three types of interactions were modeled with the cut and shifted Lennard-Jones potential

$${U}_{{{{{{{{\rm{LJ}}}}}}}}}({{{{{{{\bf{r}}}}}}}})=4\epsilon \left[{\left(\frac{{\sigma }_{l}}{{r}_{ij}}\right)}^{12}-{\left(\frac{{\sigma }_{l}}{{r}_{ij}}\right)}^{6}\right]+{E}_{{{{{{{{\rm{cut}}}}}}}}}$$
(1)

for rij < rc and zero otherwise, where rc = 2.0 σ. $${E}_{{{{{{{{\rm{cut}}}}}}}}}=-4\epsilon [{(\frac{{\sigma }_{l}}{{r}_{c}})}^{12}-{(\frac{{\sigma }_{l}}{{r}_{c}})}^{6}]$$. ϵ = 1.8, 1.8, and 1.0 kBT, and σl = 0.5, 0.75 and 0.75 σ for nucleolar particle-nucleolar particle, nucleolar particle-NAD and nucleolar particle-non-NAD interactions. The values for σl were chosen based on the size of nucleolar particles (σp = 0.5σ) and chromatin beads (σ) with arithmetic averaging. NADs were identified using the high-resolution sequencing data generated from ref. 41. We processed the raw nucleolar-to-genomic ratios to generate signal data at the 1 Mb resolution. Only genomic loci with signals higher than a threshold value (15.0) were labeled as NADs and homologous chromosomes share identical NADs.

We used the software package LAMMPS83 to perform molecular dynamics (MD) simulations in reduced units. Constant-temperature (T = 1.0 in reduced unit) simulations were carried out via the Langevin dynamics with a damping coefficient γ = 10.0 τB and a simulation time step of dt = 0.008 τB, where τB is the Brownian time unit. Configurations were recorded at every 2000 simulation steps for analysis. The initial configuration of MD simulations was built as follows. We first obtained chromosome conformations from the end configuration of a 20-million-step long simulation of the genome-only model carried out in our previous study50. Next, 500 nucleolar particles were placed with random positions inside the spherical confinement. We further relaxed the system to avoid steric clashes by performing 400,000 step MD simulations. ϵ in Eq. (1) was set to 1.0 kBT for both specific and non-specific interactions to prevent cluster formation during the equilibration period. The last configuration of this equilibration simulation was then used to initialize the phase separation simulations. Analyses of simulation trajectories were carried out with in-house Python and MATLAB scripts.

We mapped the reduced units in MD simulations to physical units for direct comparison with experimental measurements. The length scale unit was determined by assuming a typical nucleus size with a radius of 5 μm. Since the spherical confinement mimicking the nucleus in our model adopts a radius of 19.7 σ, we have σ = 0.254 μm for the reduced length unit. We further determined the reduced time unit with the expression $${\tau }_{B}=\frac{30\pi \eta {\sigma }^{3}}{{k}_{{{{B}}}}T}$$, which was obtained by matching the diffusion coefficient from simulations with that in the nucleus (see Supplementary Material). Using a value of the nucleoplasmic viscosity as η = 10−2 Pas, we estimated τB = 3.6 s.

We note that the energy unit (kBT) in our model should be viewed as an effective temperature instead of strictly the biological value (310 K). Since the nucleus is a non-equilibrium system with constant perturbations from molecular motors and chemical reactions, the ensemble of genome organizations collected over a population of cells as in Hi-C experiments is unlikely to be in thermodynamic equilibrium. However, non-equilibrium distributions can often be well approximated by renormalized Boltzmann distributions with effective potentials and temperatures84. It is these effective quantities that we inferred from Hi-C contact maps to describe the interactions among chromosomes.

### Estimating the size and number of nucleolar particles

The number of nucleolar particles used in simulations was selected based on the concentration of a representative nucleolar protein NPM1 as follows. For nucleoli of 2RNu in diameter, and a protein concentration of 1 μM85, the number of nucleolar particles can be estimated as $${N}_{{{{{{{{\rm{p}}}}}}}}}=\frac{4\pi }{3}\cdot {N}_{A}\cdot {R}_{{{{{{{{\rm{Nu}}}}}}}}}^{3}\cdot 1\,\mu$$M, where NA is the Avogadro constant. We used 500 particles in simulations, and the corresponding RNu = 0.58 μm matches well with the observed size of nucleoli11,13.

The size of the nucleolar particles can be determined assuming a space-filling model for the nucleoli as

$$\frac{\frac{4\pi }{3}\cdot {\left(\frac{{2}^{\frac{1}{6}}\cdot {\sigma }_{{{{P}}}}}{2}\right)}^{3}\cdot {N}_{{{{P}}}}}{\frac{4\pi }{3}\cdot {R}_{{{N}}\,}^{3}}={\left(\frac{{R}_{{{{Nu}}}}}{{R}_{{{{N}}}}}\right)}^{3}.$$
(2)

Using NP = 500, RN = 19.7, σ = 5 μm and RNu = 0.58 μm, we have the diameter of nucleolar particles σp = 0.5 σ. 21/6σp is the equilibrium distance between neighboring nucleolar particles in the Lennard–Jones potential.

We note that the above estimation is crude since many additional proteins and RNA molecules are present inside the nucleoli. Therefore, the nucleolar particles should be viewed as molecular aggregates with size σP = 0.5σ–0.12 μm, rather than a single protein molecule. Given the size of a typical protein as 5–10 nm86, the number of molecules within a single coarse-grained particle can be on the order of 103. This number, while large, is, in fact, in the same order as the amount of distinct molecules that make up the nucleoli87.

The molecular aggregate interpretation of nucleolar particles is also consistent with our use of ϵ = 1.8 kBT for specific interactions. This value appears much smaller than the strength of typical non-bonded interactions, such as hydrogen bonds. ϵ should be viewed as “free energy” that accounts for the averaging over the heterogeneous interaction pattern of different molecules and the entropic penalty arising from confining individual molecules together. Furthermore, as stated above, the effective temperature in our model is presumably higher than 310 K, resulting in more significant values for ϵ.

We computed the free energy profile for coalescence using umbrella sampling and temperature replica exchange with 16 umbrella windows88,89. We defined the collective variable as $${R}_{g}=\sqrt{\frac{1}{N}\mathop{\sum }\nolimits_{i = 1}^{N}{\left|{{{{{{{{\bf{r}}}}}}}}}_{i}-{{{{{{{{\bf{r}}}}}}}}}_{{{\mbox{com}}}}\right|}^{2}}$$. ri is the Cartesian coordinate of the i-th nuclear particles found inside one of the two droplets and rcom is the center of mass. Indices of nucleolar particles in the two droplets were identified using the DBSCAN algorithm (see Supplementary Materials). A harmonic potential $$\frac{K}{2}\cdot {\left({R}_{g}-{R}_{c}\right)}^{2}$$ was introduced in each umbrella window to facilitate the sampling of configurations at targeted distances. We chose the center of these windows, Rc, to be evenly spaced between 2.5 σ and 6.0 σ with an increment of 0.25 σ, except for the first window whose Rc = 2.0 σ. The spring constant K was chosen as follows:

$$K=\left\{\begin{array}{ll}100.0&2.00\le {R}_{0} \, < \, 3.50\\ 150.0&3.50\le {R}_{0}\, < \, 4.00\\ 200.0&{R}_{0}\ge 4.00\hfill\end{array}\right.$$
(3)

Eight replicas were used within each umbrella window with temperatures ranging from 1.00 to 1.14 with 0.02 increments. Exchanges between these replicas were attempted at every 100-time steps. As shown in Supplementary Fig. 21, our choice of the temperature grid allows frequent exchange among replicas.

These simulations were initialized from a typical two-droplet configuration recorded at the end of a dynamical simulation and lasted twelve million steps. Configurations were recorded at every 400 steps. To compute the error bars and evaluate simulation convergence, we divided the simulation trajectories into five consecutive blocks with equal length. Free energy profiles were then calculated using only data collected from each block with the weighted histogram method (WHAM), and error bars were determined as the standard deviation of the five independent estimations.

### Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

The data that support this study are available from the corresponding author upon reasonable request. Simulation trajectories generated in this study have been deposited in the Zenodo database (https://doi.org/10.5281/zenodo.5570927). Source data are provided with this paper.

## Code availability

The computer code used in this study is available at https://github.com/ZhangGroup-MITChemistry/NucleoliCoalescence

## References

1. Strom, A. R. & Brangwynne, C. P. The liquid nucleome–phase transitions in the nucleus at a glance. J. Cell Sci. 132, jcs235093 (2019).

2. Henninger, J. E. et al. RNA-mediated feedback control of transcriptional condensates. Cell 184, 207–225 (2021).

3. Dundr, M. Nuclear bodies: multifunctional companions of the genome. Curr. Opin. Cell Biol. 24, 415–422 (2012).

4. Banani, S. F., Lee, H. O., Hyman, A. A. & Rosen, M. K. Biomolecular condensates: organizers of cellular biochemistry. Nat. Rev. Mol. Cell Biol. 18, 285–298 (2017).

5. Brangwynne, C. P. et al. Germline p granules are liquid droplets that localize by controlled dissolution/condensation. Science 324, 1729–1732 (2009).

6. Mao, Y. S., Zhang, B. & Spector, D. L. Biogenesis and function of nuclear bodies. Trends Genet. 27, 295–306 (2011).

7. Boeynaems, S. et al. Protein phase separation: a new phase in cell biology. Trends Cell Biol. 28, 420–435 (2018).

8. Shin, Y. & Brangwynne, C. P. Liquid phase condensation in cell physiology and disease. Science 357, eaaf4382 (2017).

9. Lafontaine, D. L., Riback, J. A., Bascetin, R. & Brangwynne, C. P. The nucleolus as a multiphase liquid condensate. Nat. Rev. Mol. Cell Biol. 1–18 (2020).

10. Courchaine, E. M., Lu, A. & Neugebauer, K. M. Droplet organelles? EMBO J. 35, 1603–1612 (2016).

11. Caragine, C. M., Haley, S. C. & Zidovska, A. Nucleolar dynamics and interactions with nucleoplasm in living cells. Elife 8, e47533 (2019).

12. Brangwynne, C. P., Mitchison, T. J. & Hyman, A. A. Active liquid-like behavior of nucleoli determines their size and shape in Xenopus laevis oocytes. Proc. Natl. Acad. Sci. 108, 4334–4339 (2011).

13. Caragine, C. M., Haley, S. C. & Zidovska, A. Surface fluctuations and coalescence of nucleolar droplets in the human cell nucleus. Phys. Rev. Lett. 121, 148101 (2018).

14. Kim, J., Han, K. Y., Khanna, N., Ha, T. & Belmont, A. S. Nuclear speckle fusion via long-range directional motion regulates speckle morphology after transcriptional inhibition. J. Cell Sci. 132, jcs226563 (2019).

15. Berry, J., Weber, S. C., Vaidya, N., Haataja, M. & Brangwynne, C. P. Rna transcription modulates phase transition-driven nuclear body assembly. Proc. Natl. Acad. Sci. 112, E5237–E5245 (2015).

16. Larson, A. G. et al. Liquid droplet formation by hp1α suggests a role for phase separation in heterochromatin. Nature 547, 236–240 (2017).

17. Strom, A. R. et al. Phase separation drives heterochromatin domain formation. Nature 547, 241–245 (2017).

18. Hnisz, D., Shrinivas, K., Young, R. A., Chakraborty, A. K. & Sharp, P. A. A phase separation model for transcriptional control. Cell 169, 13–23 (2017).

19. Chong, S. et al. Imaging dynamic and selective low-complexity domain interactions that control gene transcription. Science 361, eaar2555 (2018).

20. Gibson, B. A. et al. Organization of chromatin by intrinsic and regulated phase separation. Cell 179, 470–484 (2019).

21. Kashchiev, D. Nucleation (Elsevier, 2000).

22. McGaughey, A. & Ward, C. Droplet stability in a finite system: consideration of the solid–vapor interface. J. Appl. Phys. 93, 3619–3626 (2003).

23. Farley, K. I., Surovtseva, Y., Merkel, J. & Baserga, S. J. Determinants of mammalian nucleolar architecture. Chromosoma 124, 323–331 (2015).

24. Heride, C. et al. Distance between homologous chromosomes results from chromosome positioning constraints. J. Cell Sci. 123, 4063–4075 (2010).

25. Royo, F. et al. Spatial link between nucleoli and expression of the zac1 gene. Chromosoma 118, 711 (2009).

26. Kind, J. et al. Single-cell dynamics of genome-nuclear lamina interactions. Cell 153, 178–192 (2013).

27. Grosch, M., Ittermann, S., Shaposhnikov, D. & Drukker, M. Chromatin-associated membraneless organelles in regulation of cellular differentiation. Stem Cell Rep. 15, 1220-1232 (2020).

28. Spector, D. L. & Lamond, A. I. Nuclear speckles. Cold Spring Harb. Perspect. Biol. 3, a000646 (2011).

29. Weber, S. C. & Brangwynne, C. P. Inverse size scaling of the nucleolus by a concentration-dependent phase transition. Curr. Biol. 25, 641–646 (2015).

30. Dillinger, S., Straub, T. & Nemeth, A. Nucleolus association of chromosomal domains is largely maintained in cellular senescence despite massive nuclear reorganisation. PLoS One 12, e0178821 (2017).

31. Martin, C. et al. Lamin b1 maintains the functional plasticity of nucleoli. J. Cell Sci. 122, 1551–1562 (2009).

32. Zwicker, D., Hyman, A. A. & Juelicher, F. Suppression of Ostwald ripening in active emulsions. Phys. Rev. E 92, 012317 (2015).

33. Falahati, H. & Wieschaus, E. Independent active and thermodynamic processes govern the nucleolus assembly in vivo. Proc. Natl. Acad. Sci. 114, 1335–1340 (2017).

34. Wurtz, J. D. & Lee, C. F. Chemical-reaction-controlled phase separated drops: formation, size selection, and coarsening. Phys. Rev. Lett. 120, 078102 (2018).

35. Söding, J., Zwicker, D., Sohrabi-Jahromi, S., Boehning, M. & Kirschbaum, J. Mechanisms for active regulation of biomolecular condensates. Trends Cell Biol. 30, 4–14 (2020).

36. Lee, D. S., Wingreen, N. S. & Brangwynne, C. P. Chromatin mechanics dictates subdiffusion and coarsening dynamics of embedded condensates. Nature Physics. 17, 531–538 (2021).

37. Zhang, Y., Lee, D. S., Meir, Y., Brangwynne, C. P. & Wingreen, N. S. Mechanical frustration of phase separation in the cell nucleus by chromatin. Phys. Rev. Lett. 126, 258102 (2021).

38. Ronceray, P., Mao, S., Košmrlj, A. & Haataja, M. P. Liquid demixing in elastic networks: cavitation, permeation, or size selection? arXivpreprint https://arxiv.org/abs/2102.02787 (2021).

39. McStay, B. Nucleolar organizer regions: genomic ‘dark matter requiring illumination. Genes Dev. 30, 1598–1610 (2016).

40. Németh, A. & Längst, G. Genome organization in and around the nucleolus. Trends Genet. 27, 149–156 (2011).

41. van Koningsbruggen, S. et al. High-resolution whole-genome sequencing reveals that specific chromatin domains from most human chromosomes associate with nucleoli. Mol. Biol. Cell 21, 3735–3748 (2010).

42. Boisvert, F.-M., van Koningsbruggen, S., Navascués, J. & Lamond, A. I. The multifunctional nucleolus. Nat. Rev. Mol. Cell Biol. 8, 574–585 (2007).

43. Galganski, L., Urbanek, M. O. & Krzyzosiak, W. J. Nuclear speckles: molecular organization, biological function and role in disease. Nucleic Acids Res. 45, 10350–10368 (2017).

44. Khelifi, G. & Hussein, S. M. A new view of genome organization through RNA directed interactions. Front. Cell Dev. Biol. 8, 517 (2020).

45. Bronstein, I. et al. Transient anomalous diffusion of telomeres in the nucleus of mammalian cells. Phys. Rev. Lett. 103, 018102 (2009).

46. Weber, S. C., Theriot, J. A. & Spakowitz, A. J. Subdiffusive motion of a polymer composed of subdiffusive monomers. Phys. Rev. E 82, 011913 (2010).

47. Kepten, E., Bronshtein, I. & Garini, Y. Ergodicity convergence test suggests telomere motion obeys fractional dynamics. Phys. Rev. E 83, 041919 (2011).

48. Burnecki, K. et al. Universal algorithm for identification of fractional Brownian motion. A case of telomere subdiffusion. Biophys. J. 103, 1839–1847 (2012).

49. Lucas, J. S., Zhang, Y., Dudko, O. K. & Murre, C. 3D trajectories adopted by coding and regulatory DNA elements: first-passage times for genomic interactions. Cell 158, 339–352 (2014).

50. Qi, Y. et al. Data-driven polymer model for mechanistic exploration of diploid genome organization. Biophys. J. 119, 1905–1916 (2020).

51. Johnstone, S. E. et al. Large-scale topological changes restrain malignant progression in colorectal cancer. Cell 182, 1474–1489 (2020).

52. Lieberman-Aiden, E. et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science 326, 289–293 (2009).

53. Rao, S. S. et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell 159, 1665–1680 (2014).

54. Zhang, B. & Wolynes, P. G. Topology, structures, and energy landscapes of human chromosomes. Proc. Natl. Acad. Sci. 112, 6062–6067 (2015).

55. Qi, Y. & Zhang, B. Predicting three-dimensional genome organization with chromatin states. PLoS Comput. Biol. 15, e1007024 (2019).

56. Phair, R. D. & Misteli, T. High mobility of proteins in the mammalian cell nucleus. Nature 404, 604–609 (2000).

57. Chen, D. & Huang, S. Nucleolar components involved in ribosome biogenesis cycle between the nucleolus and nucleoplasm in interphase cells. J. Cell Biol. 153, 169–176 (2001).

58. Stevens, T. J. et al. 3D structures of individual mammalian genomes studied by single-cell Hi-C. Nature 544, 59–64 (2017).

59. Vertii, A. et al. Two contrasting classes of nucleolus-associated domains in mouse fibroblast heterochromatin. Genome Res. 29, 1235–1249 (2019).

60. Cremer, T. & Cremer, M. Chromosome territories. Cold Spring Harb. Perspect. Biol. 2, a003889–a003889 (2010).

61. Bickmore, W. A. & van Steensel, B. Genome architecture: domain organization of interphase chromosomes. Cell 152, 1270–1284 (2013).

62. Tjong, H. et al. Population-based 3D genome structure analysis reveals driving forces in spatial genome organization. Proc. Natl. Acad. Sci. 113, E1663–E1672 (2016).

63. Boyle, S. et al. The spatial organization of human chromosomes within the nuclei of normal and emerin-mutant cells. Hum. Mol. Genet. 10, 211–9 (2001).

64. van Gunsteren, W. F., Daura, X. & Mark, A. E. Computation of free energy. Helvetica Chim. Acta 85, 3113–3129 (2002).

65. Parsons, T. & Zhang, B. Critical role of histone tail entropy in nucleosome unwinding. J. Chem. Phys. 150, 185103 (2019).

66. Du, R., Pande, V. S., Grosberg, A. Y., Tanaka, T. & Shakhnovich, E. S. On the transition coordinate for protein folding. J. Chem. Phys. 108, 334–350 (1998).

67. Peters, B. Reaction coordinates and mechanistic hypothesis tests. Annu. Rev. Phys. Chem. 67, 669–690 (2016).

68. Eggers, J., Lister, J. R. & Stone, H. A. Coalescence of liquid drops. J. Fluid Mech. 401, 293–310 (1999).

69. Paulsen, J. D., Burton, J. C. & Nagel, S. R. Viscous to inertial crossover in liquid drop coalescence. Phys. Rev. Lett. 106, 114501 (2011).

70. Flory, P. J. Principles of Polymer Chemistry (Cornell University Press, 1953).

71. Ester, M. et al. A density-based algorithm for discovering clusters in large spatial databases with noise. In KDD, vol. 96, 226–231 (1996).

72. Onuki, A. Phase Transition Dynamics (Cambridge University Press, 2002).

73. Ratke, L. & Voorhees, P. W.Growth and Coarsening: Ostwald Ripening in Material Processing (Springer Science & Business Media, 2013).

74. Stanich, C. A. et al. Coarsening dynamics of domains in lipid membranes. Biophys. J. 105, 444–454 (2013).

75. Chen, Y. et al. Mapping 3D genome organization relative to nuclear compartments using TSA-Seq as a cytological ruler. J. Cell Biol. 217, 4025–4048 (2018).

76. Su, J.-H., Zheng, P., Kinrot, S. S., Bintu, B. & Zhuang, X. Genome-scale imaging of the 3D organization and transcriptional activity of chromatin. Cell 182, 1641–1659 (2020).

77. Zhang, L. et al. TSA-seq reveals a largely conserved genome organization relative to nuclear speckles with small position changes tightly correlated with gene expression changes. Genome Res. 31, 251–264 (2021).

78. Fei, J. et al. Quantitative analysis of multilayer organization of proteins and RNA in nuclear speckles at super resolution. J. Cell Sci. 130, 4180–4192 (2017).

79. Huang, K. et al. Physical and data structure of 3D genome. Sci. Adv. 6, eaay4055 (2020).

80. Style, R. W. et al. Liquid-liquid phase separation in an elastic network. Phys. Rev. X 8, 011028 (2018).

81. Rosowski, K. A. et al. Elastic ripening and inhibition of liquid-liquid phase separation. Nat. Phys. 16, 422–425 (2020).

82. Zidovska, A. Chromatin: liquid or solid? Cell 183, 1737–1739 (2020).

83. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).

84. Jiang, Z. & Zhang, B. Theory of active chromatin remodeling. Phys. Rev. Lett. 123, 208102 (2019).

85. Zhu, L. et al. Controlling the material properties and rRNA processing function of the nucleolus using light. Proc. Natl Acad. Sci. 116, 17330–17335 (2019).

86. Lee, H. H. et al. Crystal structure of human nucleophosmin-core reveals plasticity of the pentamer–pentamer interface. Proteins: Struct. Funct. Bioinforma 69, 672–678 (2007).

87. Scherl, A. et al. Functional proteomic analysis of human nucleolus. Mol. Biol. Cell 13, 4100–4109 (2002).

88. Sugita, Y. & Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 314, 141–151 (1999).

89. Souaille, M. & Roux, B. Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations. Computer Phys. Commun. 135, 40–57 (2001).

## Acknowledgements

We thank Dr. Xingcheng Lin, Dr. Kartik Kamat, and Dr. Xinqiang Ding for helpful discussions. This work was supported by the National Institutes of Health (Grant R35GM133580 to B.Z.) and by the Gordon and Betty Moore Foundation (Grant GBMF9162.15 to B.Z.).

## Author information

Authors

### Contributions

Y.Q. and B.Z. designed and performed the research, analyzed the data, and wrote the paper.

### Corresponding author

Correspondence to Bin Zhang.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Peer review information Nature Communications thanks the anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Qi, Y., Zhang, B. Chromatin network retards nucleoli coalescence. Nat Commun 12, 6824 (2021). https://doi.org/10.1038/s41467-021-27123-9

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41467-021-27123-9