Phase separation in soft repulsive polymer mixtures: foundation and implication for chromatin organization

Naoki Iso,Yuki Norizoe,Takahiro Sakaue
DOI: https://doi.org/10.1039/D4SM00309H
2024-09-04
Abstract:Given the wide range of length scales, the analysis of polymer systems often requires coarse-graining, for which various levels of description may be possible depending on the phenomenon under consideration. Here, we provide a super-coarse grained description, where polymers are represented as a succession of mesosopic soft beads which are allowed to overlap with others. We then investigate the phase separation behaviors in a mixture of such homopolymers based on mean-field theory, and discuss universal aspects of the miscibility phase diagram in comparison with the numerical simulation. We also discuss an extension of our analysis to mixtures involving random copolymers, which might be interesting in the context of chromatin organization in a cell nucleus.
Soft Condensed Matter
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the theoretical description of phase - separation behavior in soft - repulsive polymer mixtures and its impact on chromatin organization. Specifically: 1. **Theoretical description of the phase - separation phenomenon**: The author proposes a super - coarse - grained description method, in which polymers are represented as a series of mesoscopic soft beads that can overlap with each other. Based on the mean - field theory, the phase - separation behavior of these homopolymer mixtures was studied, and the characteristics of the phase diagrams under different parameter conditions were explored. 2. **Comparison with traditional theories**: The paper discusses how this coarse - grained model can be extended to mixtures containing random copolymers, which is of great significance for understanding the organization of chromatin in the cell nucleus. The author points out that the standard Flory - Huggins theory may not be applicable in this case because the beads are allowed to overlap, resulting in more complex phase - separation behavior. 3. **Simulation of chromatin organization**: The author also discusses how this model can be applied to simulate the chromatin organization in the cell nucleus. In particular, in early embryos, chromatin is relatively evenly distributed in the cell nucleus, but as the development stage progresses, characteristic structures such as heterochromatin foci and transcription factories will appear. The formation of these structures may be driven by phase - separation caused by local chromatin monomer changes. ### Formula analysis Some key formulas are used in the paper to describe the free energy and phase behavior of the system: 1. **Free - energy density**: \[ \frac{f}{k_{B}T}=\frac{c_{a}}{N_{a}}\ln c_{a}+\frac{c_{b}}{N_{b}}\ln c_{b}+\frac{1}{2}\chi_{aa}c_{a}^{2}+\frac{1}{2}\chi_{bb}c_{b}^{2}+\chi_{ab}c_{a}c_{b} \] where \(k_{B}T\) is the thermal energy, \(c_{x}\) and \(N_{x}\) are the bead number density and chain length of component \(x\), respectively, and \(\chi_{xy}\) represents the repulsive strength between beads \(x\) and \(y\). 2. **Phase - equilibrium conditions**: - Equal chemical potentials: \(\mu_{a}^{(1)}=\mu_{a}^{(2)}\) and \(\mu_{b}^{(1)}=\mu_{b}^{(2)}\) - Equal pressures: \(P^{(1)} = P^{(2)}\) 3. **Universal features of the phase diagram**: - By introducing the normalized concentrations \(\tilde{c}_{a}=c_{a}N_{a}\chi_{aa}\) and \(\tilde{c}_{b}=c_{b}N_{b}\chi_{bb}\), the free - energy density can be rewritten as: \[ \frac{f}{k_{B}T}=\frac{1}{N_{a}^{2}\chi_{aa}}\left[\tilde{c}_{a}\ln\tilde{c}_{a}+k_{1}\tilde{c}_{b}\ln\tilde{c}_{b}+\frac{1}{2}\tilde{c}_{a}^{2}+\frac{1}{2}k_{1}\tilde{c}_{b}^{2}+k_{2}\tilde{c}_{a}\tilde{c}_{b}\right] \] where the coefficients \(k_{1}\) and \(k_{2}\) are defined as: \[ k_{1}=\frac{N_{a}^{2}\chi_{aa}}{N_{b}^{2}\chi_{bb}}, \quad k_{2}=\frac{N_{a}\chi_{ab}}{N_{b}\sqrt{\chi_{aa}\chi_{bb}}} \] 4. **Relationship in the Gaussian core model**: - The form of the Gaussian potential is: \[ \frac{U(r)}{k_{B}T}=\epsilon\exp\left(-\frac{r^{2}}{R^{2}}\right)