Magic numbers in polymer phase separation { the importance of being rigid Bin Xu,1 Guanhua He,2 Benjamin G. Weiner,1 Pierre Ronceray,3 Yigal Meir,4 Martin C. Jonikas,2 and Ned S. Wingreen2, 5 1Department of Physics, Princeton University 2Department of Molecular Biology, Princeton University 3Princeton Center for Theoretical Science, Princeton University 4Department of Physics, Ben Gurion University of the Negev 5Lewis-Sigler Institute for Integrative Genomics, Princeton University Cells possess non-membrane-bound bodies, many of which are now understood as phase-separated condensates. One class of such condensates is composed of two polymer species, where each consists of repeated binding sites that interact in a one-to-one fashion with the binding sites of the other polymer. Previous biologically-motivated modeling of such a two-component system surprisingly revealed that phase separation is suppressed for certain combinations of numbers of binding sites. This phenomenon, dubbed the \magic-number effect, occurs if the two polymers can form fully- bonded small oligomers by virtue of the number of binding sites in one polymer being an integer multiple of the number of binding sites of the other. Here we use lattice-model simulations and analytical calculations to show that this magic-number effect can be greatly enhanced if one of the polymer species has a rigid shape that allows for multiple distinct bonding conformations. Moreover, if one species is rigid, the effect is robust over a much greater range of relative concentrations of the two species. Our findings advance our understanding of the fundamental physics of two-component polymer-based phase-separation and suggest implications for biological and synthetic systems. Introduction A striking theoretical prediction regarding these two- component multivalent systems is that, in the regime of strong binding, condensation will be extremely sen- In addition to membrane-bound organelles, cells pos- sitive to the relative valence of the two components7. sess non-membrane-bound bodies including nucleoli, P- Higher valence is normally expected to boost conden- bodies, and stress granules, which are now understood as sate formation8. In strongly bound two-component sys- phase-separated condensates1{4. Typically, the compo- tems, however, an exception occurs when the valence of nents of these condensates have a high rate of exchange one species equals or is an integral multiple of the va- with the surrounding medium and the condensates them- lence of the other species. In this case, condensation is selves are dynamic, rapidly assembling and disassembling suppressed in favor of small fully bonded oligomers7 { in response to specific stimuli5{7. a \magic-number effect reminiscent of the exact filling One biologically relevant class of condensates are of atomic shells leading to the unreactive noble gases. those composed of two species of multivalent polymers Here, we demonstrate that this magic-number effect still or particles with specific interactions that drive phase occurs if all components are flexible polymers, but the separation8{10. In the simplest case, each component effect is maximized if one of the components is rigid consists of repeated domains that interact in a one- and compact (aka a patchy particle or patchy colloid17), to-one fashion with the domains of the other compo- in which case fully-bonded oligomers are more entropi- nent. Such two-component multivalent condensates have cally favored over the condensate. Strikingly, while the been observed in several natural and engineered con- magic-number effect requires rather precise global stoi- texts. One example, the algal pyrenoid, is a carbon- chiometry when all polymers are flexible, if one of the fixation organelle, in which the two components essen- polymers is rigid the effect occurs over a broad range tial for assembly11 are the rigid oligomer Rubisco (the of stoichiometries. While many intracellular conden- active enzyme in CO fixation) and EPYC1, an un- 2 sates are held together by weak interactions of multi- structured linker protein12. Another multivalent system, ple types (e.g. charged, aromatic, and hydrophobic18, PML (promyelocytic leukemia) nuclear bodies that re- as well as pi-cation19, and pi-pi interactions20) natural arXiv:1901.09352v2 [physics.bio-ph] 11 Oct 2019 pair DNA damage13, relies on the Small Ubiquitin-like protein-protein, protein-RNA, or protein-peptide inter- Modifier (SUMO) domain that interacts with the SUMO actions such as SUMO-SIM10 or synthetic interactions, Interacting Motif (SIM) to form droplets14{16. This sys- e.g. based on DNA hybridization21,22, can readily reach tem inspired in vitro experiments with engineered poly- the strong-binding regime required to observe the magic- SUMO and polySIM of various valences10, and, indeed, number effect. phase separation was observed with increasing concentra- tions of the two polymers. Other in vitro two-component systems, e.g., an engineered polySH3-Proline-Rich Motif system8 and a PTB-RNA system9, also form droplets as concentrations are increased. 2 Results systems? In the fully bonded regime, the transition is not driven by a competition between entropy and energy as in Flory-Huggins theory25, but rather by a competition The magic-number effect with one rigid and one between different types of entropy. For example, in the flexible component droplet phase of the A :B system, free dimers coexist To examine the role of rigidity in the magic-number 8 8R with large gel-like clusters. Each dimer has high trans- effect, we begin with two species that interact: a flexi- lational entropy, while each component within a cluster ble linear polymer and a rigid shape. The flexible poly- has very limited translational entropy. However, this dif- mer with n binding sites is denoted as A and the rigid n ference in translational entropies is offset by an opposing 4 × 2 shape is denoted as B ; each binding site on each 8R difference in conformational entropies: in each dimer, the species is considered to be a \monomer". Previous re- binding sites of one species must match all the binding sults from modeling the Rubisco-EPYC1 system, which sites of the other species, leading to a reduced total con- established the concept of the magic-number effect, can formational entropy. By contrast, the components in a then be summarized as follows: the A :B system forms 4 8R large condensed cluster are more independent, binding stable trimers (each composed of two A s and one B ) 4 8R to multiple members of the other species and enjoying and thus, compared to the A :B or A :B systems, 3 8R 5 8R a relatively high conformational entropy. Because the requires a substantially higher total monomer concentra- importance of translational entropy relative to confor- tion for the formation of large clusters, assuming equal mational entropy is reduced as concentration increases, total stoichiometry of A and B monomers7. the system transitions from a uniform dimer phase to Based on these observations, we predicted that flexi- a droplet phase with increasing concentration. By con- ble polymers of length 8 together with rigid 4 × 2 shapes trast, in non-magic-number systems (e.g.,A :B or should also constitute a magic-number system, in this 7 8R A :B ), the polymers cannot form free dimers without case based on stable heterodimers, since each constituent 9 8R unsatisfied binding sites, and instead tend to form large has 8 binding sites. To test this idea, we simulated fully-bonded clusters even at low concentrations. This A :B ,A :B , and A :B systems (Fig. 1a-c). The 7 8R 8 8R 9 8R key difference leads to the drastic difference in clustering heat maps of the average cluster size for different con- between magic-number and non-magic-number systems. centrations and specific bond strengths (Fig. 1g-i) define While the simulations in Fig. 1 are performed in two di- approximate phase diagrams: regions with smaller clus- mensions, because the magic-number effect depends only ter sizes correspond to a homogeneous dissolved phase, on the ability to form small fully bonded oligomers, it while regions with larger cluster sizes correspond to an also occurs in three dimensions (Supplementary Fig. 3). inhomogeneous condensed phase, i.e., a phase-separated In three dimensions, we observe the onset of the magic- regime. Confirming the visual impression that the ob- number effect around a specific interaction energy of served large clusters arise from phase separation, rather 4k T for components of valence ∼ 8. We note that in than arising as percolation clusters from homogeneous B three dimensions, phase separating systems can become gelation23, we found that for strong enough specific in- trapped in a state of interconnected “fibers", but we con- teractions (≥ 4k T ) the internal cluster density remains B firmed that our annealing procedure leads instead to a approximately constant above the transition, as expected stable droplet phase (Supplementary Fig. 5). for phase separation but not for clusters formed by perco- To mimic weak attractive interactions such as hy- lation (Supplementary Fig. 4). Moreover, the dense clus- drophobicity, we include a small non-specific bond en- ters we observe are always internally connected by a net- ergy in all our simulations, which leads to more com- work of specific bonds. We conclude that our system rep- pact droplets. However, these non-specific interactions resents a case of gelation driven by phase separation23,24. are neither sufficient nor necessary for the magic-number Notice that in Fig. 1g-i, as specific interactions increase effect, though their magnitude influences the location of from zero, cluster sizes initially increase, but stop increas- the phase boundary (Supplementary Figs. 6 and 7). ing above ∼ 7k T . This behavior is due to saturation of B To verify that in the A :B magic-number system the the specific bonds, i.e., for large enough specific bond 8 8R phase transition in the fully-bonded regime is driven by energies essentially all possible bonds are formed.