Significance: LLPS plays an important role in the spatiotemporal organization of the numerous molecular constituents of living cells, via formation of biomolecular condensates. Our simulations provide predictive rules governing the stability and composition of multicomponent biomolecular condensates. Biomolecules that increase the molecular connectivity of condensates are present in higher concentrations because connectivity is positively correlated with stability. Greater connectivity within highly multicomponent condensates manifests in higher critical temperatures in the phase diagrams of accessible systems involving just a few components. Hence, composition of highly multicomponent condensates can be predicted from the critical points of reduced-component mixtures. Our findings expand the mechanisms relating phase behavior of multicomponent intracellular mixtures to critical parameters (temperature, pH, salt concentration, etc.) of the constituent biomolecules.

Biomolecular condensates, some of which are liquid-like during health, can age over time becoming gel-like pathological systems. One potential source of loss of liquid-like properties during ageing of RNA-binding protein condensates is the progressive formation of inter-proteinβ-sheets. To bridge microscopic understanding between accumulation of inter-proteinβ-sheets over time and the modulation of FUS and hnRNPA1 condensate viscoelasticity, we develop a multiscale simulation approach. Our method integrates atomistic simulations with sequence-dependent coarse-grained modelling of condensates that exhibit accumulation of inter-proteinβ-sheets over time. We reveal that inter-proteinβ-sheets notably increase condensate viscosity but does not transform the phase diagrams. Strikingly, the network of molecular connections within condensates is drastically altered, culminating in gelation when the network of strongβ-sheets fully percolates. However, high concentrations of RNA decelerate the emergence of inter-proteinβ-sheets. Our study uncovers molecular and kinetic factors explaining how the accumulation of inter-proteinβ-sheets can trigger liquid-to-solid transitions in condensates, and suggests a potential mechanism to slow such transitions down.

Biomolecular condensates, which assemble via the process of liquid–liquid phase separation (LLPS), are multicomponent compartments found ubiquitously inside cells. Experiments and simulations have shown that biomolecular condensates with many components can exhibit multilayered organizations. Using a minimal coarse-grained model for interacting multivalent proteins, we investigate the thermodynamic parameters governing the formation of multilayered condensates through changes in protein valency and binding affinity. We focus on multicomponent condensates formed by scaffold proteins (high-valency proteins that can phase separate on their own via homotypic interactions) and clients (proteins recruited to condensates via heterotypic scaffold–client interactions). We demonstrate that higher valency species are sequestered to the center of the multicomponent condensates, while lower valency proteins cluster towards the condensate interface. Such multilayered condensate architecture maximizes the density of LLPS-stabilizing molecular interactions, while simultaneously reducing the surface tension of the condensates. In addition, multilayered condensates exhibit rapid exchanges of low valency proteins in and out, while keeping higher valency proteins—the key biomolecules involved in condensate nucleation—mostly within. We also demonstrate how modulating the binding affinities among the different proteins in a multicomponent condensate can significantly transform its multilayered structure, and even trigger fission of a condensate into multiple droplets with different compositions.

Proteins containing intrinsically disordered regions (IDRs) are ubiquitous within biomolecular condensates, which are liquid-like compartments within cells formed through liquid–liquid phase separation (LLPS). The sequence of amino acids of a protein encodes its phase behaviour, not only by establishing the patterning and chemical nature (e.g., hydrophobic, polar, charged) of the various binding sites that facilitate multivalent interactions, but also by dictating the protein conformational dynamics. Besides behaving as random coils, IDRs can exhibit a wide-range of structural behaviours, including conformational switching, where they transition between alternate conformational ensembles. Using Molecular Dynamics simulations of a minimal coarse-grained model for IDRs, we show that the role of protein conformation has a non-trivial effect in the liquid–liquid phase behaviour of IDRs. When an IDR transitions to a conformational ensemble enriched in disordered extended states, LLPS is enhanced. In contrast, IDRs that switch to ensembles that preferentially sample more compact and structured states show inhibited LLPS. This occurs because extended and disordered protein conformations facilitate LLPS-stabilising multivalent protein–protein interactions by reducing steric hindrance; thereby, such conformations maximize the molecular connectivity of the condensed liquid network. Extended protein configurations promote phase separation regardless of whether LLPS is driven by homotypic and/or heterotypic protein–protein interactions. This study sheds light on the link between the dynamic conformational plasticity of IDRs and their liquid–liquid phase behaviour.

Biomolecular condensates formed via liquid–liquid phase separation (LLPS) play a crucial role in the spatiotemporal organization of the cell material. Nucleic acids can act as critical modulators in the stability of these protein condensates. To unveil the role of RNA length in regulating the stability of RNA binding protein (RBP) condensates, we present a multiscale computational strategy that exploits the advantages of a sequence-dependent coarse-grained representation of proteins and a minimal coarse-grained model wherein proteins are described as patchy colloids. We find that for a constant nucleotide/protein ratio, the protein fused in sarcoma (FUS), which can phase separate on its own—i.e., via homotypic interactions—only exhibits a mild dependency on the RNA strand length. In contrast, the 25-repeat proline-arginine peptide (PR₂₅), which does not undergo LLPS on its own at physiological conditions but instead exhibits complex coacervation with RNA—i.e., via heterotypic interactions—shows a strong dependence on the length of the RNA strands. Our minimal patchy particle simulations suggest that the strikingly different effect of RNA length on homotypic LLPS versus RBP–RNA complex coacervation is general. Phase separation is RNA-length dependent whenever the relative contribution of heterotypic interactions sustaining LLPS is comparable or higher than those stemming from protein homotypic interactions. Taken together, our results contribute to illuminate the intricate physicochemical mechanisms that influence the stability of RBP condensates through RNA inclusion.