One of the key mechanisms used by cells to control the spatiotemporal organization of their many components is the formation and dissolution of biomolecular condensates through liquid–liquid phase separation (LLPS). Using a minimal coarse-grained model that allows us to simulate thousands of interacting multivalent proteins, we investigate the physical parameters dictating the stability and composition of multicomponent biomolecular condensates. We demonstrate that the molecular connectivity of

more »

... the condensed-liquid network—i.e., the number of weak attractive protein–protein interactions per unit of volume—determines the stability (e.g., in temperature, pH, salt concentration) of multicomponent condensates, where stability is positively correlated with connectivity. While the connectivity of scaffolds (biomolecules essential for LLPS) dominates the phase landscape, introduction of clients (species recruited via scaffold–client interactions) fine-tunes it by transforming the scaffold–scaffold bond network. Whereas low-valency clients that compete for scaffold–scaffold binding sites decrease connectivity and stability, those that bind to alternate scaffold sites not required for LLPS or that have higher-than-scaffold valencies form additional scaffold–client–scaffold bridges increasing stability. Proteins that establish more connections (via increased valencies, promiscuous binding, and topologies that enable multivalent interactions) support the stability of and are enriched within multicomponent condensates. Importantly, proteins that increase the connectivity of multicomponent condensates have higher critical points as pure systems or, if pure LLPS is unfeasible, as binary scaffold–client mixtures. Hence, critical points of accessible systems (i.e., with just a few components) might serve as a unified thermodynamic parameter to predict the composition of multicomponent condensates.