Triblock copolymer micelle model of spherical paraspeckles

Abstract

Paraspeckles are nuclear bodies scaffolded by RNP complexes of NEAT1_2 RNA transcripts and multiple RNA-binding proteins. The assembly of paraspeckles is coupled with the transcription of NEAT1_2. Paraspeckles form the core-shell structure, where the two terminal regions of NEAT1_2 RNP complexes compose the shell of the paraspeckle and the middle regions of these complexes compose the core. We here construct a theoretical model of paraspeckles by taking into account the transcription of NEAT1_2 in an extension of the theory of block copolymer micelles. This theory predicts that the core-shell structure of a paraspeckle is assembled by the association of the middle region of NEAT1_2 RNP complexes due to the multivalent interactions between RBPs bound to these regions and by the relative affinity of the terminal regions of the complexes to the nucleoplasm. The latter affinity results in the effective repulsive interactions between terminal regions of the RNA complexes and limits the number of complexes composing the paraspeckle. In the wild type, the repulsive interaction between the middle and terminal block dominates the thermal fluctuation. However, the thermal fluctuation can be significant in a mutant, where a part of the terminal regions of NEAT1_2 is deleted, and distributes the shortened terminal regions randomly between the shell and the core, consistent with our recent experiments. With the upregulated transcription, the shortened terminal regions of NEAT1_2 in a deletion mutant is localized to the core to decrease the repulsive interaction between the terminal regions, while the structure does not change with the upregulation in the wild type. The robustness of the structure of paraspeckles in the wild type results from the polymeric nature of NEAT1_2 complexes.

Introduction

A cell nucleus is not a uniform solution of DNA and proteins, but there are a number of nuclear bodies in the interchromatin space (Chujo et al., 2016; Nakagawa et al., 2018; Palikyaras and Papantonis 2019; Banani et al., 2017; Berry et al., 2015; Van Treeck and Parker 2018). Some nuclear bodies are scaffolded by RNAs that make ribonucleoprotein (RNP) complexes with RNA-binding proteins (RBPs). A class of RNAs that are essential in assembling nuclear bodies are called architectural RNAs (arcRNAs) (Chujo et al., 2016). Growing number of evidences suggest that nuclear bodies are condensates assembled by liquid-liquid phase separation (LLPS) because of the multivalent interactions between the intrinsically disordered regions of RBPs that are bound to arcRNAs (Chujo et al., 2016; Nakagawa et al., 2018; Palikyaras and Papantonis 2019; Banani et al., 2017; Berry et al., 2015; Van Treeck and Parker 2018). Condensates produced by LLPS are thought to act as reaction crucibles of specific biochemical reactions, molecular sponges that sequester specific proteins and RNAs from the nucleoplasm, and hubs to organize the 3D structure of genome (Shin and Brangwynne 2017). It is of interest to study the assembly mechanism of nuclear bodies due to the possible relationship with their physiological functions.

Paraspeckles are nuclear bodies that are scaffolded by NEAT1_2 arcRNA (Sunwoo et al., 2008; Clemson et al., 2009; Sasaki et al., 2009). Paraspeckles act as molecular sponges that sequester some types of RNA transcripts and proteins (Chen and Carmichael 2009; Hirose et al., 2014; Imamura et al., 2014; Hu et al., 2015; West et al., 2016; Nakagawa et al., 2018; Wang et al., 2018; Yap et al.,2022) and interact with chromatin regions, enriched in active promoters and enhancers (West et al., 2014; Li et al., 2017; Sridhar et al., 2017; Fang et al., 2019; Bonetti et al., 2020; Cai et al., 2020). De novo assembled paraspeckles are often observed at the proximity to the transcription site of NEAT1_2 and are disassembled when the transcription of arcRNA is suppressed (Clemson et al., 2009; Sasaki et al., 2009; Mao et al., 2011). The number of paraspeckles increases when transcription of NEAT1_2 is upregulated (Clemson et al., 2009; Hirose et al., 2014). These results imply that the assembly of paraspeckles is coupled with the transcription of NEAT1_2. RNA-binding proteins can be bound to nascent NEAT1_2 transcripts, which are still connected to the transcription site via RNA polymerase II (Pol II). The array of nascent NEAT1_2 RNAs produced during a transcription burst probably act as a nucleation site of paraspeckles (Chujo and Hirose 2017; Yamazaki et al., 2020).

We have recently extended the Flory-Huggins theory, which is the standard theory of phase separation of polymers in a solution (Doi 1996), to predict the phase separation driven by the production of arcRNAs due to transcription (Yamamoto et al., 2020). The condensates assembled by this mechanism are disordered liquids of complexes of arcRNAs and RBPs. However, paraspeckles are not condensates of disordered liquid, but form a characteristic core-shell structure (Souquere et al., 2010; West et al., 2016): the two terminal regions and the middle region of NEAT1_2 form the shell and the core of paraspeckles, respectively. The structure of paraspeckles is analogous to micelles of ABC triblock copolymers in a selective solvent, where two polymer chains (A and C blocks) composed of hydrophilic units are chemically bonded to the two ends of a chain (B block) composed of hydrophobic units (Monzen et al., 2000; Mai and Eisenberg 2012; Moughton et al., 2012). The ordered structure and the transcription driven formation of paraspeckles are two features that distinguish paraspeckles from condensates assembled by the classical phase separation, such as LLPS, in the thermodynamic equilibrium.

In our recent experiments, we have constructed mutant NEAT1_2 cell lines, in which the terminal regions of NEAT1_2 were partly deleted by CRISPR/Cas9 and have observed paraspeckles in such cell lines with the super-resolution optical microscope and the electron microscope (Yamazaki et al., 2018, 2021). Our experiments have shown that the terminal regions of NEAT1_2 RNAs were localized in the shells of paraspeckles in the wild type (WT) (Figure 1A), whereas the terminal regions were distributed both to the core and the shell in the mutant cells (Figure 1B) (the schematics of NEAT1_2 for each case is shown in Figure 1C). Motivated by this result, we here construct a model of paraspeckles by taking into account the transcription dynamics of NEAT1_2 in an extension of the theory of micelles of ABC triblock copolymers. The A and C blocks correspond to the terminal regions of NEAT1_2 RNP complexes and the B block corresponds to their middle region. The B blocks are associative due to the multivalent interactions between RBPs, such as NONO and FUS, that specifically bind to this B block region due to its sequence (Yamazaki et al., 2018).

FIGURE 1

In many cases, paraspeckles are approximately spherical, whereas cylindrical paraspeckles are also observed, for example, when NEAT1_2 transcription is upregulated (Hirose et al., 2014). The analysis of cylindrical paraspeckles greatly complicates the theory, while the essence of the biophysical mechanism of the assembly of paraspeckles is already in the theory of spherical paraspeckles. Not to hide the essence behind the complexity of analysis, in this paper, we limit our discussion to spherical paraspeckles.

Our theory predicts the distribution of terminal blocks (A blocks) and the size of paraspeckles in wild type and mutant with deleted terminal blocks and also the effect of the upregulation of NEAT1_2 transcription. Our prediction is consistent with our recent experiments (Yamazaki et al., 2018, 2021), implying that the assembly of paraspeckles can be understood as micellization. Our theory provides biophysical insight into the assembly of WT and mutant paraspeckles. The sequences of arcRNAs determine the arrangements of RBPs along these arcRNAs and thus play a role in the blueprints of nuclear bodies, while their assembly is fine-tuned by the transcription dynamics of arcRNAs. Our theory can be extended to understand the mechanism of the assembly of other nuclear bodies once the arrangements of RBPs along arcRNAs are determined by experiments.

Materials and methods

Model

We treat the NEAT1_2 RNP complexes as ABC triblock copolymers (Figure 2A). Different RBPs are bound to A, B, and C blocks and this makes the magnitudes of interactions between units in the same blocks different from the magnitudes of the interactions between different blocks. Treating a complex as one polymer, instead of treating arcRNA and RBPs separately, is effective for the case in which the binding affinity of RNA-binding proteins to these blocks is relatively large (Yamamoto et al., 2020). The A, B, and C blocks are composed of , , and (Kuhn) units, respectively, (Figure 2A). The core of a paraspeckle is assembled by the association of the B blocks because of the attractive interactions between the RBPs bound to these blocks (Figure 2B). The A and C blocks form the shell because the RBPs bound to these terminal blocks have affinity to the solution (nucleoplasm), rather than the RBPs bound to B blocks. The analysis of the cylindrical paraspeckles greatly complicates the theory because of the massive form of the free energy, an extra geometrical parameter that should be determined by minimizing the free energy, and the lack of steady state. In this paper, we focus on spherical paraspeckles.

FIGURE 2

In the thermodynamic equilibrium, the most stable state is the one at the minimum of the free energy of the system. The free energy of a micelle of block copolymers is composed of 1) the stretching free energy of the blocks in the core, 2) the surface free energy at the interface between the core and the shell, and 3) the free energy of the blocks in the shell (Halperin and Alexander 1989; Semenov et al., 1995; Zhulina et al., 2005). The growth of polymer micelles decreases the surface free energy and increases the stretching free energy of blocks in the core and the free energy due to the excluded volume interactions between blocks in the shell (Halperin and Alexander 1989; Semenov et al., 1995; Zhulina et al., 2005). The stable size of spherical micelles is determined by the balance of these free energy contributions. Most theories predict the distribution of the size of polymer micelles that are assembled uniformly in a block copolymer solution at the thermodynamic equilibrium (Annianson and Wall 1974; Safran 2003; Hadgiivanova et al., 2011; Mysona et al., 2019). We extend the theory of polymer micelles by taking into account the fact that paraspeckles are assembled locally at the transcription site of NEAT1_2 and their assembly is coupled with the transcription of NEAT1_2.

Free energy

The free energy quantifies the stability of the system at the thermodynamic equilibrium. We here derive the free energy of a paraspeckle by taking into account the fact that the terminal blocks of NEAT1_2 are distributed both to the shell and the core in an extension of the free energy of a polymer micelle (Halperin and Alexander 1989; Semenov et al., 1995; Zhulina et al., 2005). The free energy of a paraspeckle composed of triblock copolymers has the formwhere is the free energy of the core, is the surface free energy at the interface between the core and the shell, is the free energy of the shell, and is the mixing free energy (Figure 3).

FIGURE 3

The free energy of the core has the form

The first and second terms are the stretching free energy of the B blocks in the core and are derived in the spirit of Semenov 1985, see Supplementary Section S1. The third term is the free energy due to the interactions between the A and B units in the core. For cases in which the number of the A blocks is relatively large, these blocks may aggregate to minimize the repulsive interactions between A and B units. However, to keep the simplicity of the theory, we derived the third term of Eqn. 2 by assuming that the fraction of the A blocks are distributed uniformly in the core, independent of the number of units in A blocks. With the assumption with which the core is packed with A and B blocks, the volume has the formwhere is the radius of the core. is the geometrical factor (Semenov 1985). The interaction parameter accounts for the excluded volume interactions between the A blocks and the B blocks. and are the volume fractions of the A blocks and the B blocks in the core, respectively. The core is occupied by the A and B blocks, .

The surface free energy has the form

The interaction parameter accounts for the (free) energetic penalty due to the fact that the B units at the surface have less number of interacting partners than in the interior of the core. Without changing the physics, we treat cases in which the A and C units are dilute both in the core and the shell. In such cases, the interactions between the B units at the core surface and the A units in the shell as well as the interactions between the A units at the core surface and the solution are both negligible.

Theories of block copolymer micelles treat blocks in the shell as a polymer brush on a curved surface (Halperin and Alexander 1989; Semenov et al., 1995; Zhulina et al., 2005). There are elaborate approaches to treat polymer brushes on planer (Netz and Schick 1998) and curved surfaces (Zhulina et al., 2006), however, for simplicity, we here use the scaling theory of polymer brush (Alexander 1977; de Gennes 1980). With this approach, the free energy of the shell has the formwith

Eqn. 6 includes the stretching free energy of the A blocks in the shell and the free energy due to the excluded volume interactions between A units in the shell. Eqn. 7 is the corresponding free energy for the C blocks. We neglected the excluded volume interaction between the A and C units because the 3′ and 5’ terminal regions are segregated in separate domains in the shell of wild type paraspeckles (West et al., 2016). Eqs. 6, 7 are derived by using so-called Daoud-Cotton scaling theory (Daoud and Cotton 1982), assuming that the solution is a good solvent to both A and C blocks, following the usual treatment of block copolymer micelles (Halperin and Alexander 1989; Semenov et al., 1995; Zhulina et al., 2005), see also in Supplementary Section S2.1 is the fraction of the core surface occupied by A blocks and has the form

Eqn. 8 is derived by using the fact that the osmotic pressure in the domains of A blocks is equal to the osmotic pressure in the domains of C blocks, see Supplementary Section S2.1. and are the heights of A and C blocks in the limit of planer brush, , respectively. These heights have the forms (Alexander 1977; de Gennes 1980) and are the excluded volumes that account for the excluded volume interactions between the A units and those between the C units, respectively. is a numerical constant of order unity and is determined as by curvefitting the experiments on the micelles of polystylene-polyisoprene copolymers with the scaling theory (Zhulina et al., 2005).

The mixing free energy represents the entropic contribution that distributes the A blocks randomly to the shell and the core and has the form

The free energy is a function of the number n of triblock copolymers comprising the paraspeckle and the fraction of A blocks in the shell. We determine the fraction by the minimization of the free energy . This corresponds to cases in which the time scale of the redistribution of A blocks is smaller than the time scale of the production of transcripts.

Association and dissociation dynamics of NEAT1_2

Most theories of polymer micelles predict the most stable size and the distribution of the size of polymer micelles that are assembled uniformly in a block copolymer solution at the thermodynamic equilibrium (Annianson and Wall 1974; Safran 2003; Hadgiivanova et al., 2011; Mysona et al., 2019). In contrast, the assembly of paraspeckles is coupled with the transcription of NEAT1_2 (Mao et al., 2011). RBPs can bind to nascent NEAT1_2 transcripts, which are still connected to the transcription site via Pol II, and an array of complexes of nascent NEAT1_2 transcripts and RBPs during a transcription burst act as a nucleation site of paraspeckles. Once a paraspeckle is nucleated, nascent NEAT1_2 transcripts are added to the paraspeckle one by one as nascent transcripts are produced by Pol II. This situation may be somewhat analogous to the assembly of micelles due to the association of side chains that are attached to a main chain, except for the fact that these side chains are not produced and released dynamically. The side chains assemble micelles without translational entropy cost, in contrast to the micellization of polymer chains freely diffusing in the solution (Borisov and Halperin 1995). We thus derive the time evolution equation of a paraspeckle at the transcription site by taking into account the facts that 1) nascent NEAT1_2 transcripts are associated with the paraspeckle without the translational entropy cost and 2) paraspeckles are assembled locally at the transcription site. Because of these features, paraspeckles are different from block copolymer micelles assembled in the thermodynamic equilibrium (Annianson and Wall 1974; Safran 2003; Hadgiivanova et al., 2011; Mysona et al., 2019).

We treat the probability that the paraspeckle assembled at the transcription site of NEAT1_2 is composed of transcripts at time . The time evolution of the probability has the form

The flux has the form

The first term of Eqn. 14 is the rate with which transcripts diffusing in the solution are spontaneously associated with the paraspeckle. The second term is the rate with which a transcript is spontaneously dissociated from the paraspeckle. The third term is the rate with which a nascent transcript during production is added to the paraspeckle through transcription. is the rate constant that accounts for the association of transcripts diffusing in the solution with the paraspeckle. is the volume fraction of transcripts that are not associated to paraspeckles. is the free energy which is already minimized with respect to the fraction , see Eqn. 1. is the rate with which a nascent transcript is added to the paraspeckle. The form of the third term of Eqn. 14 represents the fact that nascent transcripts can be associated with the paraspeckle without the translational entropy cost. We note that the probability is the local quantity of the paraspeckle assembled at the transcription site, in contrast to the usual treatment of micelles that predict the global distribution function of the size of micelles in a solution (Annianson and Wall 1974; Safran 2003; Hadgiivanova et al., 2011; Mysona et al., 2019).

For cases in which the transcription is suppressed, , the probability in the equilibrium has the form

Substituting Eqn. 15 into Eqn. 14 leads to for that ensure the detailed balance at the thermodynamic equilibrium.

In the steady state, , the solution of Eqs. 12, 13 has the form

The effective free energy has the formwhere is the ratio of the association rates. is the effective partition function

The number of transcripts composing a paraspeckle with the maximum probability in the steady state is derived by minimizing the effective free energy with respect to . In this paper, we treat cases in which NEAT1_2 transcripts are exclusively localized in paraspeckles, (Sasaki et al., 2009; Yamazaki et al., 2021).

Results

Transcription dynamics regulates the structure of paraspeckles

We first discuss the dependence of the fraction of A blocks in the shell on the production rate of NEAT1_2 transcripts. The number of transcripts in a paraspeckle is derived as a function of by minimizing the effective free energy , see Eqn. 17, and the fraction of A blocks in the shell is derived by minimizing the free energy for this number , see Eqn. 1. When the number of units in the A blocks is smaller than a critical value , the fraction decreases continuously with increasing the production rate (the cyan and light green lines in Figure 4). When the number of units in the A blocks is larger than the critical value , the paraspeckles are the ‘all-A shell’ state, in which all the A blocks are in the shell, , for low transcription rate and the ‘all-A core’ state, in which all the A blocks are in the core, , for high transcription rate (the orange and magenta lines in Figure 4). There is a discontinuous transition between the all-A shell and all-A core states at a threshold value of production rate. This result is summarized in the phase diagram (Figure 5). The ‘mixed’ state for , where A blocks are distributed between the core and the shell, results from the thermal fluctuation, which is quantified by the mixing free energy . For , the interaction free energy, and , where both scale proportional to , dominates the mixing free energy , which is independent of and thus only the all-A shell and all-A core states are possible (see a quantitative argument below). It is analogous to the Flory-Huggins theory that predicts that the interaction free energy dominates the mixing free energy in polymer systems (Doi 1996). The fraction of A blocks in the shell is sensitive to the upregulation of the transcription of NEAT1_2 transcripts when , whereas the fraction does not change significantly by the moderate upregulation of the transcription when (Figure 4). It is because for , the interaction free energy dominates the mixing free energy and thus the changes of the fraction is suppressed.

FIGURE 4

FIGURE 5

In the WT, most of the terminal regions of NEAT1_2 are localized at the shell (Souquere et al., 2010; West et al., 2016; Yamazaki et al., 2018, 2021). This implies that paraspeckles in the WT are in the ‘all-A shell’ state, which happens for and relatively low transcription rate. In deletion mutants, a fraction of the terminal blocks, which were partially deleted, is localized at the core, implying that paraspeckles in the deletion mutants correspond to the case of . Our recent experiments have shown that in deletion mutants, the fraction of the terminal regions, which were partially deleted, localized at the shell decreases by the upregulation of the transcription of NEAT1_2, whereas in the WT, the fraction does not change with the upregulation of transcription (Yamazaki et al., 2021). These predictions are consistent with the prediction of our theory. However, we note that in the WT, the fraction of cylindrical paraspeckles increases by the upregulation (Hirose et al., 2014) and the morphological transition to cylinder may also be involved in the insensitivity.

Our theory predicts the biophysical mechanism of the assembly of paraspeckles. The fact that all the terminal blocks are localized at the shell in the WT paraspeckles results from the strong repulsive interaction between A and B blocks, where its influence to the structures of paraspeckles is quantified by the free energy . Indeed, the free energy of A and C blocks in the shell and the free energy due to the stretching of B blocks in the core both decrease as the fraction of A blocks in the shell decreases. However, the free energy dominates probably because the magnitude of the interaction between A and B blocks and the number of units in each B block is large enough in WT NEAT1_2. The influence of the thermal fluctuation that randomly distributes A blocks between the core and the shell to the structure is quantified by the mixing free energy and it is independent of the number of units in A blocks. In contrast, the free energy due to the interactions between A and B blocks in the core increases as the number of units in each A block increases. For , the interaction free energy dominates the mixing free energy and thus the ‘all-A shell’ state becomes stable. In contrast, for , the mixing free energy is still significant and thus A blocks are distributed between the shell and the core. This explains the difference of the structures of paraspeckles between WT and deletion mutants (Figure 1 and Yamazaki et al., 2018, 2021).

The number of transcripts in the paraspeckle increases with increasing the transcription rate. For a relatively small transcription rate, all-A shell state becomes stable because the free energy due to the repulsive interactions between A and B blocks in the core (see Eqn. 2 for ) dominates the free energy due to the excluded volume interactions between A units in the shell and the stretching free energy of blocks in the core (see Eqn. 2 for ). The free energy of the blocks in the shell and the stretching free energy of blocks in the core increase faster than the free energy due to the repulsive interactions between A and B blocks in the core as the number of transcripts in the paraspeckle increases. This results in the decrease of the fraction of the A blocks in the shell with increasing the transcription rate (see the Supplementary Discussion for the relative significance of the free energy contributions and ). For very small production rate, stable paraspeckles are not assembled, see the region delineated by the red line in Figure 5.

Repulsive interactions of terminal blocks and entropic elasticity of middle blocks limit the number of transcripts in paraspeckles

The size of paraspeckles in the WT and deletion mutants is experimentally accessible. The radius of the core is derived by using Eqn. 3. The radii, and , are derived by using the formswhere the heights, and , are given in Eqs. 9, 10. The derivation of Eqs. 19, 20 are shown in Supplementary Section S2.1. These radii are functions of the number of transcripts in a paraspeckle and the fraction of blocks in the shell, where the latter parameters are derived similarly to Figure 4. At a first glance, one may think that the radius of the paraspeckle decreases as the number of units in A blocks decreases. However, our theory predicts that the radius of paraspeckles (defined by the distance between the top of the A or C blocks and the center of the paraspeckle) increases with decreasing the number of units in A blocks, see Figure 6A. It is because the number of transcripts in the paraspeckle increases with decreasing the number of units in A blocks, see Figure 6B. This prediction is consistent with our recent experiments (Yamazaki et al., 2021).

FIGURE 6

The number of transcripts in a paraspeckle is limited by the free energy of blocks in the shell and the stretching free energy of blocks in the core, analogous to micelles of diblock copolymers. These free energy contributions decreases as the fraction of A blocks in the shell decreases. Indeed, the fraction of A blocks in the shell decreases with decreasing the number of units in the A blocks, see Figure 6C. More transcripts can therefore associate with paraspeckles as the number of units in A the blocks decreases. These predictions are consistent with our recent experiments, see Discussion and Yamazaki et al., 2021.

Paraspeckles are assembled by the association of the middle regions

Paraspeckles are assembled by the attractive interactions between the B blocks of NEAT1_2 RNP complexes. As expected, the radius of paraspeckles decreases as the number of units in the B blocks decreases, see Figure 7A. With a fixed transcription rate, the number of transcripts in the paraspeckle decreases with decreasing the number of units in the B blocks. It is because decreasing the number of units in the B blocks increases the density of A and C blocks in the shell and thus increases the free energy of A and C blocks in the shell, see Figure 7B. It also increases the extent of the stretching of B blocks in the core and thus increases the stretching free energy of the B blocks in the core. The fraction of A blocks in the shell increases as the number of units in B blocks decreases because A blocks in the core suppress the attractive interactions between the B blocks, see Figure 7C. The formation of paraspeckles is suppressed for cases in which the number of units in the B blocks is too small, see Figure 8. These predictions are consistent with our experiments that show that paraspeckles of mutant cells, in which a part of the middle region of NEAT1_2 is deleted, are small and dispersed, compared with the WT, see also the Discussion and Yamazaki et al., 2018.

FIGURE 7

FIGURE 8

Discussion

We have constructed a theory of paraspeckle assembly by taking into account the transcription dynamics of NEAT1_2 in an extension of the theory of ABC triblock copolymer micelles. This model captures two features of paraspeckles: paraspeckles form the characteristic core-shell structure and the assembly of paraspeckles is coupled with the transcription of NEAT1_2. Our theory provides a biophysical insight into the assembly of paraspeckles. Paraspeckles are assembled by multivalent interactions between the middle blocks of the NEAT1_2 RNP complexes (Yamazaki et al., 2018). The repulsive interactions between A and B blocks, quantified by the interaction free energy , expels A blocks from the core, while the thermal fluctuations, quantified by the mixing free energy , distribute A blocks randomly between the core and the shell. All the A blocks are localized in the shell for the WT, where the number N_A of units in A block is large enough, because the interaction free energy dominates the mixing free energy . This reflects the connectivity of RBPs via the terminal regions of NEAT1_2 and is a well-known feature of polymeric molecules (Doi 1996). This also accounts for the fact that the all-A shell structure of WT paraspeckles is not sensitive to the upregulation of transcription, even for spherical paraspeckles.

Our recent experiments have revealed that (I) the 5′ terminal regions of NEAT1_2 transcripts were distributed randomly between the shell and the core in paraspeckles of Δ5′ deletion mutants, where the 5′ terminal region of NEAT1_2 was deleted by 2 kb, (II) the 3′ terminal region of most NEAT1_2 transcripts was localized in the core of paraspeckles of Δ3′ deletion mutants, where the 3′ terminal region of NEAT1_2 was deleted by 6 kb, (III) the sizes of paraspeckles in Δ5′ and Δ3′ deletion mutants were larger than the size of wild type paraspeckles, (IV) paraspeckles of deletion mutant, where the middle region of NEAT1_2 was deleted by 8.6 kb, were smaller and more dispersed than the wild type paraspeckles, (V) in deletion mutants, the fraction of the terminal blocks localized in the shell decreased when the transcription of NEAT1_2 was upregulated (Yamazaki et al., 2018, 2021). The core-shell structure of paraspeckles is disorganized in Δ5′- Δ3′ double deletion mutants (Yamazaki et al., 2018, 2021). The assembly of such disordered paraspeckles has been studied theoretically in our recent research (Yamamoto et al.,2020). The summary of these results is shown in Figure 9. The predictions of our theory are consistent with these experimental results. The assembly of paraspeckles can be therefore viewed as the micellization of NEAT1_2 RNP complexes, at least in the first approximation. More quantitative comparison requires the characterization of the interaction parameters between each pair of units by, for example, osmotic pressure measurements (Mangenot et al., 2002) and scattering techniques (Oohashi et al., 2014).

FIGURE 9

Our theory provides a biophysical insight into our experimental results: In the WT, all the A blocks are localized in the shell because the free energy due to the repulsive interaction between A and B blocks dominates the mixing free energy , see the first paragraph of this Discussion. However, the mixing free energy becomes significant as the number of units in A blocks decreases. Moreover, the excluded volume interactions between A blocks and those between C blocks, quantified by the free energy , as well as the stretching of B blocks in the core, quantified by the free energy , drive the relocation of A blocks from the shell to the core. This explains our experimental results (I) and (II), see Figure 9. Because these free energy contributions, and , limit the number of NEAT1_2 transcripts that form paraspeckles, the paraspeckles of deletion mutant can accommodate more NEAT1_2 transcripts. The size of paraspeckles can increase with decreasing the length of the terminal regions of NEAT1_2 because the number of NEAT1_2 transcripts in a paraspeckle increases, see Figure 9. In contrast, the number of NEAT 1_2 transcripts decreases with decreasing the length of the B blocks because the free energy decreases due to the attractive interactions between B blocks by the association of NEAT1_2 RNP complexes becomes less significant than the free energy increase due to the repulsive interactions between A blocks and those between C blocks. The radius of paraspeckles thus decreases as the number of units in B blocks decreases, see Figure 9. The number . of NEAT1_2 in a paraspeckle increases with the upregulation of the transcription of NEAT1_2. The free energy contributions, and , increase relative to the free energy with increasing the number of NEAT1_2 in a paraspeckle . This drives the relocation of A blocks from the shell to the core with the upregulation of NEAT1_2 transcription in the deletion mutant. The all-A shell configuration of a WT paraspeckle does not change with the upregulation of NEAT1_2 transcription, see the first paragraph of this Discussion.

In our approach, we take into account only the essential features to understand the mechanism of the assembly of paraspeckles, instead of doing computer simulations by including all known things. It is indeed the strategy of theoretical physics. Our theory is certainly an important first step, but there are a couple of remaining mysteries. First, we focused on the analysis of spherical paraspeckles, not to hide the essence of the assembly of paraspeckles behind the complexity in treating cylindrical paraspeckles. We thus did not explain the sphere-cylinder morphological transition of paraspeckles. The fact that the fraction of cylindrical paraspeckles increases with the upregulation of NEAT1_2 transcription is indeed consistent with our conclusion that paraspeckles are assembled by micellization. Second, in the WT, 3′ and 5’ terminal regions are not randomly mixed, but are separated in microdomains (West et al., 2016). Our present theory takes into account this feature in the free energy of the shell by neglecting the interaction between A and C blocks, but did not explain it theoretically. Identifying the RBPs bound to A and C blocks will greatly help to understand the mechanism of the assembly of the microdomains. One possible explanation of the assembly of microdomains is that RBPs bound to A and C are different and the interactions between different blocks are more repulsive than the interactions between the same blocks. However, the microdomains are probably assembled by the microphase separation, judging from the fact that multiple microdomains at the shell of a paraspeckle do not show fusion, and the A-C interaction alone is not enough to explain the assembly of the microdomains. The super-resolution microscope experiments suggest that multiple NEAT1_2 RNPs form a bundle (West et al., 2016). This bundling probably plays an important role in the assembly of the microdomains of A and C blocks. Third, we used the steady state approximation to derive the distribution of the number of NEAT1_2 transcripts in a paraspeckle. However, the number of NEAT1_2 transcripts that can be incorporated in a paraspeckle should be limited by the number of NEAT1_2 transcripts produced in one transcription burst. The experiments that study the relationship between the transcription dynamics and the number of NEAT1_2 transcripts per paraspeckles would greatly help our understanding of the assembly mechanism of paraspeckles.

Conclusion

By the combination of theory and experiments, we have shown that paraspeckles are assembled by micellization, not liquid-liquid phase separation. One important feature of micelles is that their size is regulated by the balance between the surface free energy (that drives the growth of micelles) and the free energy due to the repulsive interactions between blocks in the shell (that limits the growth of micelles). The free energy due to the entropic elasticity of blocks in the core also limits the growth of micelles. The number of paraspeckles per cell is usually larger than the number of the transcription sites of NEAT1, implying that paraspeckles can diffuse far away from the transcription sites. Paraspeckles may take advantage of this size control mechanism to gain the mobility toward target sites through the meshwork of chromatin in a nucleus.

Our theory and experiments provide insight into the general principle of the assembly of nuclear bodies: The pattern of RBPs binding to arcRNAs is tailored into their RNA sequences and the RNP complexes behave as copolymers that direct the ordering in the internal structures of nuclear bodies. Our theory can be therefore extended to understand the mechanism of the assembly of other nuclear bodies by using the experimentally determined arrangement of RBPs along arcRNAs. The assembly of nuclear bodies is facilitated by the transcription of arcRNAs. The number of arcRNAs per nuclear body is controlled by both the interactions and connectivity of RBPs bound to the arcRNAs and the transcription dynamics. Chemical engineers control the stability and size of liquid condensates by using surfactants. Life takes the same strategy by using RNA.

References

• 1

  Alexander S. (1977). Adsorption of chain molecules with a polar head a scaling description. J. Phys. Fr.28, 983–987. 10.1051/jphys:01977003808098300CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 2

  Annianson E. A. G. Wall S. N. (1974). On the kinetics of step-wise micelle association. J. Phys. Chem.78, 1024–1030. 10.1021/j100603a016CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 3

  Banani S. F. Lee H. O. Hyman A. A. Rosen M. K. (2017). Biomolecular condensates: Organizers of cellular biochemistry. Nat. Rev. Mol. Cell. Biol.18, 285–298. 10.1038/nrm.2017.7PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 4

  Berry J. Weber S. C. Vaidya N. Haataja M. Brangwynne C. P. (2015). RNA transcription modulates phase transition-driven nuclear body assembly. Proc. Natl. Acad. Sci. U. S. A.112, E5237–E5245. 10.1073/pnas.1509317112PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 5

  Bonetti A. Agostini F. Suzuki A. M. Hashimoto K. Pascarella G. Gimenez J. et al (2020). RADICL-seq identifies general and cell type–specific principles of genome-wide RNA-chromatin interactions. Nat. Commun.11, 1018. 10.1038/s41467-020-14337-6PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 6

  Borisov O. V. Halperin A. (1995). Micelles of polysoaps. Langmuir11, 2911–2919. 10.1021/la00008a012CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 7

  Cai Z. Cao C. Ji L. Ye R. Wang D. Xia C. et al (2020). RIC-seq for global in situ profiling of RNA-RNA spatial interactions. Nature582, 432–437. 10.1038/s41586-020-2249-1PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 8

  Chen L. L. Carmichael G. G. (2009). Altered nuclear retention of mRNAs containing inverted repeats in human embryonic stem cells: Functional role of a nuclear noncoding RNA. Mol. Cell.35, 467–478. 10.1016/j.molcel.2009.06.027PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 9

  Chujo T. Hirose T. (2017). Nuclear bodies built on architectural long noncoding RNAs: Unifying principles of their construction and function. Mol. Cells40, 889–896. 10.14348/molcells.2017.0263PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 10

  Chujo T. Yamazaki T. Hirose T. (2016). Architectural RNAs (arcRNAs): A class of long noncoding RNAs that function as the scaffold of nuclear bodies. Biochim. Biophys. Acta1859, 139–146. 10.1016/j.bbagrm.2015.05.007PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 11

  Clemson C. M. Hutchinson J. N. Sara S. A. Ensminger A. W. Fox A. H. Chess A. et al (2009). An architectural role for a nuclear noncoding RNA: NEAT1 RNA is essential for the structure of paraspeckles. Mol. Cell.33, 717–726. 10.1016/j.molcel.2009.01.026PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 12

  Daoud M. Cotton J. P. (1982). Star shaped polymers: A model for the conformation and its concentration dependence. J. Phys. Fr.43, 531–538. 10.1051/jphys:01982004303053100CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 13

  de Gennes P. G. (1980). Conformations of polymers attached to an interface. Macromolecules13, 1069–1075. 10.1021/ma60077a009CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 14

  Doi M. (1996). Introduction to polymer physics. New York, USA: Oxford Univ. Press. Google Scholar

  • Google Scholar

• 15

  Fang J. Ma Q. Chu C. Huang B. Li L. Cai P. et al (2019). PIRCh-seq: Functional classification of non-coding RNAs associated with distinct histone modifications. Genome Biol.20, 292. 10.1186/s13059-019-1880-3PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 16

  Hadgiivanova R. Diamant H. Andelman D. (2011). Kinetics of surfactant micellization: A free energy approach. J. Phys. Chem. B115, 7268–7280. 10.1021/jp1073335PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 17

  Halperin A. Alexander S. (1989). Polymeric micelles: Their relaxation kinetics. Macromolecules22, 2403–2412. 10.1021/ma00195a069CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 18

  Hirose T. Virnicchi G. Tanigawa A. Naganuma T. Li R. Kimura H. et al (2014). NEAT1 long noncoding RNA regulates transcription via protein sequestration within subnuclear bodies. Mol. Biol. Cell.25, 169–183. 10.1091/mbc.E13-09-0558PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 19

  Hu S. B. Xiang J. F. Li X. Xu Y. Xue W. Huang M. et al (2015). Protein arginine methyltransferase CARM1 attenuates the paraspeckle-mediated nuclear retention of mRNAs containing IRAlus. Genes Dev.29, 630–645. 10.1101/gad.257048.114PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 20

  Imamura K. Imamachi N. Akizuki G. Kumakura M. Kawaguchi A. Nagata K. et al (2014). Long noncoding RNA NEAT1-dependent SFPQ relocation from promoter region to paraspeckle mediates IL8 expression upon immune stimuli. Mol. Cell.53, 393–406. 10.1016/j.molcel.2014.01.009PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 21

  Li X. Zhou B. Chen L. Gou L. T. Li H. Fu X. D. et al (2017). GRID-seq reveals the global RNA-chromatin interactome. Nat. Biotechnol.35, 940–950. 10.1038/nbt.3968PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 22

  Mai Y. Eisenberg A. (2012). Self-assembly of block copolymers. Chem. Soc. Rev.41, 5969–5985. 10.1039/c2cs35115cPubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 23

  Mangenot S. Raspaud E. Tribet C. Belloni L. Livolant F. (2002). Interactions between isolated nucleosome core particles: A tail-bridging effect?Eur. Phys. J. E7, 221–231. 10.1140/epje/i200101151CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 24

  Mao Y. S. Sunwoo H. Zhang B. Spector D. L. (2011). Direct visualization of the co-transcriptional assembly of a nuclear body by noncoding RNAs. Nat. Cell. Biol.13, 95–101. 10.1038/ncb2140PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 25

  Monzen M. Kawakatsu T. Doi M. (2000). Micelle formation in triblock copolymer solutions. Comput. Theor. Polym. Sci.10, 275–280. 10.1016/s1089-3156(99)00052-5CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 26

  Moughton A. Hillmyer M. A. Lodge T. P. (2012). Multicompartment block polymer micelles. Macromolecules45, 2–19. 10.1021/ma201865sCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 27

  Mysona J. A. McCormick A. V. Morse D. C. (2019). Mechanism of micelle birth and death. Phys. Rev. Lett.123, 038003. 10.1103/PhysRevLett.123.038003PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 28

  Nakagawa S. Yamazaki T. Hirose T. (2018). Molecular dissection of nuclear paraspeckles: Towards understanding the emerging world of the RNP milieu. Open Biol.8, 180150. 10.1098/rsob.180150PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 29

  Netz R. R. Schick M. (1998). Polymer brushes: From self-consistent field theory to classical theory. Macromolecules31, 5105–5122. 10.1021/ma9717505PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 30

  Oohashi T. Inoue K. Nakamura Y. (2014). Second and third virial coefficients of low-molecular-weight polyisoprene in 1, 4-dioxane. Polym. J.46, 699–703. 10.1038/pj.2014.43CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 31

  Palikyaras S. Papantonis A. (2019). Modes of phase separation affecting chromatin regulation. Open Biol.9, 190167. 10.1098/rsob.190167PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 32

  Safran S. A. (2003). Statistical thermodynamics of surfaces, interfaces, and membranes. USA: Westview Press, CO. Google Scholar

  • Google Scholar

• 33

  Sasaki Y. T. F. Ideue T. Sano M. Mituyama T. Hirose T. (2009). MENepsilon/beta noncoding RNAs are essential for structural integrity of nuclear paraspeckles.Proc. Natl. Acad. Sci. U. S. A.106, 2525–2530. 10.1073/pnas.0807899106PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 34

  Semenov A. N. (1985). Contribution to the theory of microphase layering in block-copolymer melts. Sov. Phys. JETP61, 733–742. http://www.jetp.ras.ru/cgi-bin/e/index/e/61/4/p733?a=list. Google Scholar

  • Google Scholar

• 35

  Semenov A. N. Nyrkova I. A. Khokhlov A. R. (1995). Polymers with strongly interacting Groups: Theory for nonspherical multiplets. Macromolecules28, 7491–7500. 10.1021/ma00126a029CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 36

  Shin Y. Brangwynne C. P. (2017). Liquid phase condensation in cell physiology and disease. Science357, eaaf4382. 10.1126/science.aaf4382PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 37

  Souquere S. Beauclair G. Harper F. Fox A. Pierron G. (2010). Highly ordered spatial organization of the structural long noncoding Neat1 RNAs within paraspeckle nuclear bodies. Mol. Biol. Cell.21, 4020–4027. 10.1091/mbc.E10-08-0690PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 38

  Sridhar B. Rivas-Astroza M. Nguyen T. C. Chen W. Yan Z. Cao X. et al (2017). Systematic mapping of RNA-chromatin interactions in vivo. Curr. Biol.20, 602–609. 10.1016/j.cub.2017.01.011CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 39

  Sunwoo H. Dinger M. E. Wilusz J. E. Amaral P. P. Mattick J. S. Spector D. L. et al (2008). MEN epsilon/beta nuclear-retained non-coding RNAs are up-regulated upon muscle differentiation and are essential components of paraspeckles.Genome Res.19, 347–359. 10.1101/gr.087775.108PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 40

  Wang Y. Hu S. B. Wang M. R. Yao R. W. Wu D. Yang L. et al (2018). Genome-wide screening of NEAT1 regulators reveals cross-regulation between paraspeckles and mitochondria. Nat. Cell. Biol.20, 1145–1158. 10.1038/s41556-018-0204-2PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 41

  West J. A. Davis C. P. Sunwoo H. Simon M. D. Sadreyev R. Wang P. I. et al (2014). The long noncoding RNAs NEAT1 and MALAT1 bind active chromatin sites. Mol. Cell.55, 791–802. 10.1016/j.molcel.2014.07.012PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 42

  West J. A. Mito M. Kurosaka S. Takumi T. Tanegashima C. Chujo T. et al (2016). Structural, super-resolution microscopy analysis of paraspeckle nuclear body organization. J. Cell. Biol.214, 817–830. 10.1083/jcb.201601071PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 43

  Yamamoto T. Yamazaki T. Hirose T. (2020). Phase separation driven by production of architectural RNA transcripts. Soft Matter16, 4692–4698. 10.1039/c9sm02458aPubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 44

  Yamazaki T. Nakagawa S. Hirose T. (2020). Architectural RNAs for membraneless nuclear body formation. Cold Spring Harb. Symp. Quant. Biol.84, 227–237. 10.1101/sqb.2019.84.039404PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 45

  Yamazaki T. Souquere S. Chujo T. Kobelke S. Chong Y. S. Fox A. H. et al (2018). Functional domains of NEAT1 architectural lncRNA induce paraspeckle assembly through phase separation. Mol. Cell.70, 1038–1053. 10.1016/j.molcel.2018.05.019PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 46

  Yamazaki T. Yamamoto T. Yoshino H. Souquere S. Nakagawa S. Pierron G. et al (2021). Paraspeckles are constructed as block copolymer micelles. EMBO J.40, e107270. 10.15252/embj.2020107270PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 47

  Yap K. Chung T. H. Makeyev E. V. (2022). Hybridization-proximity labeling reveals spatially ordered interactions of nuclear RNA compartments. Mol. Cell.82, 463–478.e11. 10.1016/j.molcel.2021.10.009PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 48

  Zhulina E. B. Adam M. LaRue I. Sheiko S. S. Rubinstein M. (2005). Diblock copolymer micelles in a dilute solution. Macromolecules38, 5330–5351. 10.1021/ma048102nCrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar

• 49

  Zhulina E. B. Birshtein T. M. Borisov O. V. (2006). Curved polymer and polyelectrolyte brushes beyond the Daoud-Cotton model. Eur. Phys. J. E Soft Matter20, 243–256. 10.1140/epje/i2006-10013-5PubMed Abstract | CrossRef Full Text | Google Scholar

  • CrossRef
  • Google Scholar