Phase-separating RNA-binding proteins form heterogeneous distributions of clusters in subsaturated solutions

Macromolecular phase separation is thought to be one of the processes that drives the formation of membraneless biomolecular condensates in cells. The dynamics of phase separation are thought to follow the tenets of classical nucleation theory, and, therefore, subsaturated solutions should be devoid of clusters with more than a few molecules. We tested this prediction using in vitro biophysical studies to characterize subsaturated solutions of phase-separating RNA-binding proteins with intrinsically disordered prion-like domains and RNA-binding domains. Surprisingly, and in direct contradiction to expectations from classical nucleation theory, we find that subsaturated solutions are characterized by the presence of heterogeneous distributions of clusters. The distributions of cluster sizes, which are dominated by small species, shift continuously toward larger sizes as protein concentrations increase and approach the saturation concentration. As a result, many of the clusters encompass tens to hundreds of molecules, while less than 1% of the solutions are mesoscale species that are several hundred nanometers in diameter. We find that cluster formation in subsaturated solutions and phase separation in supersaturated solutions are strongly coupled via sequence-encoded interactions. We also find that cluster formation and phase separation can be decoupled using solutes as well as specific sets of mutations. Our findings, which are concordant with predictions for associative polymers, implicate an interplay between networks of sequence-specific and solubility-determining interactions that, respectively, govern cluster formation in subsaturated solutions and the saturation concentrations above which phase separation occurs.

Phase separation of RNA-binding proteins with disordered prion-like domains (PLDs) and RNA-binding domains (RBDs) is implicated in the formation and dissolution of membraneless biomolecular condensates such as RNA–protein (RNP) granules (1–9). Macroscopic phase separation is a process whereby a macromolecule in a solvent separates into a dilute, macromolecule-deficient phase that coexists with a dense, macromolecule-rich phase (10, 11). In a binary mixture, the soluble phase, comprising dispersed macromolecules that are well mixed with the solvent, becomes saturated at a concentration designated as c_sat. Above c_sat, for total macromolecular concentrations c_tot that are between the binodal and spinodal, phase separation of full-length RNA-binding proteins and PLDs is thought to follow classical nucleation theory (12–15).In classical nucleation theories, clusters representing incipient forms of the new dense phase form within dispersed phases of supersaturated solutions defined by c_tot > c_sat (16, 17). In the simplest formulation of classical nucleation theory (16–18), the free energy of forming a cluster of radius a is ΔF=−4π3a3Δμρn+4πa2γ. Here, Δµ is the difference in the chemical potential between the one-phase and two-phase regimes (see discussion in SI Appendix), which is negative in supersaturated solutions and positive in subsaturated solutions; ρ_n is the number of molecules per unit volume, and γ is the interfacial tension between dense and dilute phases. At temperature T, in a seed-free solution, the degree of supersaturation s is defined as s≡ΔμRT=ln(ctotcsat), where R is the ideal gas constant. Here, s is positive for c_tot > c_sat, and, as s increases, cluster formation becomes more favorable. Above a critical radius a*, the free energy of cluster formation can overcome the interfacial penalty, and the new dense phase grows in a thermodynamically downhill fashion. Ideas from classical nucleation theory have been applied to analyze and interpret the dynamics of phase separation in supersaturated solutions (12, 13, 15). Classical nucleation theories stand in contrast to two-step nucleation theories that predict the existence of prenucleation clusters in supersaturated solutions (19–22). These newer theories hint at the prospect of there being interesting features in subsaturated solutions, where c_tot < c_sat and s < 0.The subsaturated regime, where s is negative, corresponds to the one-phase regime. Ignoring the interfacial tension, the free energy of realizing clusters with n molecules in subsaturated solutions is: ΔF = –nΔµ. Therefore, the probability P(n) of forming a cluster of n molecules in a subsaturated solution is proportional to exp(sn). Accordingly, the relative probability P(n)/P(1) of forming clusters with n molecules will be exp(s(n – 1)). This quantity, which may be thought of as the concentration of clusters with n molecules, is negligibly small for clusters with more than a few molecules. This is true irrespective of the degree of subsaturation, s. Is this expectation from classical nucleation theories valid? We show here that subsaturated solutions feature a rich distribution of species not anticipated by classical nucleation theories. We report results from measurements of cluster size distributions in subsaturated solutions of phase-separating RNA-binding proteins from the FUS-EWSR1-TAF15 (FET) family. We find that these systems form clusters in subsaturated solutions, and that the cluster sizes follow heavy-tailed distributions. The abundant species are always small clusters. However, as total macromolecular concentration (c_tot) increases, the distributions of cluster sizes shift continuously toward larger values. We discuss these findings in the context of theories for associative polymers (9, 23–30).