Dynamic metastable long-living droplets formed by sticker-spacer proteins

Separation of sticker-spacer energetics governs the coalescence of metastable biomolecular condensates

Biological condensates often emerge as a multi-droplet state and never coalesce into one large droplet within the experimental timespan. Previous work revealed that the sticker-spacer architecture of biopolymers may dynamically stabilize the multi-droplet state. Here, we simulate the condensate coalescence using metadynamics approach and reveal two distinct physical mechanisms underlying the fusion of droplets. Condensates made of sticker-spacer polymers readily undergo a kinetic arrest when stickers exhibit slow exchange while fast exchanging stickers at similar levels of saturation allow merger to equilibrium states. On the other hand, condensates composed of homopolymers fuse readily until they reach a threshold density. Increase in entropy upon inter-condensate mixing of chains drives the fusion of sticker-spacer chains. We map the range of mechanisms of kinetic arrest from slow sticker exchange dynamics to density mediated in terms of energetic separation of stickers and spacers. Our predictions appear to be in excellent agreement with recent experiments probing dynamic nature of protein-RNA condensates.

Phase Separation and Aggregation of α-Synuclein Diverge at Different Salt Conditions

The coacervation of alpha-synuclein (αSyn) into cytotoxic oligomers and amyloid fibrils are considered pathological hallmarks of Parkinson's disease. While aggregation is central to amyloid diseases, liquid-liquid phase separation (LLPS) and its interplay with aggregation have gained increasing interest. Previous work shows that factors promoting or inhibiting aggregation have similar effects on LLPS. This study provides a detailed scanning of a wide range of parameters, including protein, salt and crowding concentrations at multiple pH values, revealing different salt dependencies of aggregation and LLPS. The influence of salt on aggregation under crowding conditions follows a non-monotonic pattern, showing increased effects at medium salt concentrations. This behavior can be elucidated through a combination of electrostatic screening and salting-out effects on the intramolecular interactions between the N-terminal and C-terminal regions of αSyn. By contrast, this study finds a monotonic salt dependence of LLPS due to intermolecular interactions. Furthermore, it observes time evolution of the two distinct assembly states, with macroscopic fibrillar-like bundles initially forming at medium salt concentration but subsequently converting into droplets after prolonged incubation. The droplet state is therefore capable of inhibiting aggregation or even dissolving aggregates through heterotypic interactions, thus preventing αSyn from its dynamically arrested state.

Preserving condensate structure and composition by lowering sequence complexity

Biomolecular condensates play a vital role in organizing cellular chemistry. They selectively partition biomolecules, preventing unwanted cross talk and buffering against chemical noise. Intrinsically disordered proteins (IDPs) serve as primary components of these condensates due to their flexibility and ability to engage in multivalent interactions, leading to spontaneous aggregation. Theoretical advancements are critical at connecting IDP sequences with condensate emergent properties to establish the so-called molecular grammar. We proposed an extension to the stickers and spacers model, incorporating heterogeneous, nonspecific pairwise interactions between spacers alongside specific interactions among stickers. Our investigation revealed that although spacer interactions contribute to phase separation and co-condensation, their nonspecific nature leads to disorganized condensates. Specific sticker-sticker interactions drive the formation of condensates with well-defined networked structures and molecular composition. We discussed how evolutionary pressures might emerge to affect these interactions, leading to the prevalence of low-complexity domains in IDP sequences. These domains suppress spurious interactions and facilitate the formation of biologically meaningful condensates.

Direct computations of viscoelastic moduli of biomolecular condensates

In vitro facsimiles of biomolecular condensates are formed by different types of intrinsically disordered proteins including prion-like low complexity domains (PLCDs). PLCD condensates are viscoelastic materials defined by time-dependent, sequence-specific complex shear moduli. Here, we show that viscoelastic moduli can be computed directly using a generalization of the Rouse model and information regarding intra- and inter-chain contacts that is extracted from equilibrium configurations of lattice-based Metropolis Monte Carlo (MMC) simulations. The key ingredient of the generalized Rouse model is the Zimm matrix that we compute from equilibrium MMC simulations. We compute two flavors of Zimm matrices, one referred to as the single-chain model that accounts only for intra-chain contacts, and the other referred to as a collective model, that accounts for inter-chain interactions. The single-chain model systematically overestimates the storage and loss moduli, whereas the collective model reproduces the measured moduli with greater fidelity. However, in the long time, low-frequency domain, a mixture of the two models proves to be most accurate. In line with the theory of Rouse, we find that a continuous distribution of relaxation times exists in condensates. The single crossover frequency between dominantly elastic versus dominantly viscous behaviors is influenced by the totality of the relaxation modes. Hence, our analysis suggests that viscoelastic fluid-like condensates are best described as generalized Maxwell fluids. Finally, we show that the complex shear moduli can be used to solve an inverse problem to obtain distributions of relaxation times that underlie the dynamics within condensates.

Multi-condensate state as a functional strategy to optimize the cell signaling output

The existence of multiple biomolecular condensates inside living cells is a peculiar phenomenon not compatible with the predictions of equilibrium statistical mechanics. In this work, we address the problem of multiple condensates state (MCS) from a functional perspective. We combined Langevin dynamics, reaction-diffusion simulation, and dynamical systems theory to demonstrate that MCS can indeed be a function optimization strategy. Using Arp2/3 mediated actin nucleation pathway as an example, we show that actin polymerization is maximum at an optimal number of condensates. For a fixed amount of Arp2/3, MCS produces a greater response compared to its single condensate counterpart. Our analysis reveals the functional significance of the condensate size distribution which can be mapped to the recent experimental findings. Given the spatial heterogeneity within condensates and non-linear nature of intracellular networks, we envision MCS to be a generic functional solution, so that structures of network motifs may have evolved to accommodate such configurations.

Multiscale Computational Framework for the Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins

The reversible assembly of intrinsically disordered proteins (IDPs) to form membraneless organelles (MLOs) is a fundamental process involved in the spatiotemporal regulation in living cells. MLOs formed via liquid-liquid phase separation (LLPS) serve as molecule-enhancing hubs to regulate cell functions. Owing to the complexity and dynamic nature of the protein assembly via a network of weak inter- and intra-molecular interactions, it is challenging to describe and predict the LLPS behavior. We have developed a multiscale computational model for IDPs, using the fused in sarcoma (FUS) protein and its variants as illustrative examples. To simplify the description of protein, FUS is represented as a linear chain of stickers interspaced by spacers, as inspired by the associative polymer theory. Low-complexity aromatic-rich kinked segments (LARKS) available in FUS were identified using LARKSdb and represented as "stickers". The pairwise potential energies of each pair of stickers and their β-sheet-forming propensity were estimated via molecular docking and all atomistic molecular dynamics (AA-MD) simulations. Subsequently, FUS chains were randomly positioned in a cubic lattice as coarse-grained (CG) beads, with the bead assignment based on the Kuhn length estimation of stickers and spacers. Stochastic FUS movements were modeled by Monte Carlo (MC) simulations. In addition to the Metropolis algorithm, discretized pair potential distributions between stickers were considered in the move acceptance criteria. The chosen pair potential represents one of the possible binding energy states, with its probability determined by the frequency of the binding energy distribution histogram. The fluctuations of averaged radial distribution functions (RDFs) in successive MC trial move intervals of equilibrated lattice MC simulations were used to indicate the dynamic nature of assembly/disassembly of the protein chains. This multiscale computational framework provides an economical and efficient way of predicting and describing the LLPS behavior of IDPs.

Molecular Crowding: The History and Development of a Scientific Paradigm

It is now generally accepted that macromolecules do not act in isolation but "live" in a crowded environment, that is, an environment populated by numerous different molecules. The field of molecular crowding has its origins in the far 80s but became accepted only by the end of the 90s. In the present issue, we discuss various aspects that are influenced by crowding and need to consider its effects. This Review is meant as an introduction to the theme and an analysis of the evolution of the crowding concept through time from colloidal and polymer physics to a more biological perspective. We introduce themes that will be more thoroughly treated in other Reviews of the present issue. In our intentions, each Review may stand by itself, but the complete collection has the aspiration to provide different but complementary perspectives to propose a more holistic view of molecular crowding.

Phase Separation and Aggregation of α-Synuclein Diverge at Different Salt Conditions

The coacervation and structural rearrangement of the protein alpha-synuclein (αSyn) into cytotoxic oligomers and amyloid fibrils are considered pathological hallmarks of Parkinson's disease. While aggregation is recognized as the key element of amyloid diseases, liquid-liquid phase separation (LLPS) and its interplay with aggregation have gained increasing interest. Previous work showed that factors promoting or inhibiting amyloid formation have similar effects on phase separation. Here, we provide a detailed scanning of a wide range of parameters including protein, salt and crowding concentrations at multiple pH values, revealing different salt dependencies of aggregation and phase separation. The influence of salt on aggregation under crowded conditions follows a non-monotonic pattern, showing increased effects at medium salt concentrations. This behavior can be elucidated through a combination of electrostatic screening and salting-out effects on the intramolecular interactions between the N-terminal and C-terminal regions of αSyn. By contrast, we find a monotonic salt dependence of phase separation due to the intermolecular interaction. Furthermore, we observe the time evolution of the two distinct assembly states, with macroscopic fibrillar-like bundles initially forming at medium salt concentration but subsequently converting into droplets after prolonged incubation. The droplet state is therefore capable of inhibiting aggregation or even dissolving the aggregates through a variety of heterotypic interactions, thus preventing αSyn from its dynamically arrested state.

Measurement of solubility product in a model condensate reveals the interplay of small oligomerization and self-association

Cellular condensates often consist of 10s to 100s of distinct interacting molecular species. Because of the complexity of these interactions, predicting the point at which they will undergo phase separation into discrete compartments is daunting. Using experiments and computation, we therefore studied a simple model system consisting of 2 proteins, polySH3 and polyPRM, designed for pentavalent heterotypic binding. We tested whether the peak solubility product, the product of dilute phase monomer concentrations, is a predictive parameter for the onset of phase separation. Titrating up equal total concentrations of each component showed that the maximum solubility product does approximately coincide with the threshold for phase separation in both the experiments and models. However, we found that measurements of dilute phase concentration include contributions from small oligomers, not just monomers; therefore, a quantitative comparison of the experiments and models required inclusion of small oligomers in the model analysis. We also examined full phase diagrams where the model results were almost symmetric along the diagonal, but the experimental results were highly asymmetric. This led us to perform dynamic light scattering experiments, where we discovered a weak homotypic interaction for polyPRM; when this was added to the computational model, it was able to recapitulate the experimentally observed asymmetry. Thus, comparing experiments to simulation reveals that the solubility product can be predictive of phase separation, even if small oligomers and low affinity homotypic interactions preclude experimental measurement of monomer concentration.

RNAs undergo phase transitions with lower critical solution temperatures

Co-phase separation of RNAs and RNA-binding proteins drives the biogenesis of ribonucleoprotein granules. RNAs can also undergo phase transitions in the absence of proteins. However, the physicochemical driving forces of protein-free, RNA-driven phase transitions remain unclear. Here we report that various types of RNA undergo phase separation with system-specific lower critical solution temperatures. This entropically driven phase separation is an intrinsic feature of the phosphate backbone that requires Mg2+ ions and is modulated by RNA bases. RNA-only condensates can additionally undergo enthalpically favourable percolation transitions within dense phases. This is enabled by a combination of Mg2+-dependent bridging interactions between phosphate groups and RNA-specific base stacking and base pairing. Phase separation coupled to percolation can cause dynamic arrest of RNAs within condensates and suppress the catalytic activity of an RNase P ribozyme. Our work highlights the need to incorporate RNA-driven phase transitions into models for ribonucleoprotein granule biogenesis.

Preserving condensate structure and composition by lowering sequence complexity

Biological condensates play a vital role in organizing cellular chemistry. They selectively partition biomolecules, preventing unwanted cross-talk and buffering against chemical noise. Intrinsically disordered proteins (IDPs) serve as primary components of these condensates due to their flexibility and ability to engage in multivalent, non-specific interactions, leading to spontaneous aggregation. Theoretical advancements are critical at connecting IDP sequences with condensate emergent properties to establish the so-called molecular grammar. We proposed an extension to the stickers and spacers model, incorporating non-specific pairwise interactions between spacers alongside specific interactions among stickers. Our investigation revealed that while spacer interactions contribute to phase separation and co-condensation, their non-specific nature leads to disorganized condensates. Specific sticker-sticker interactions drive the formation of condensates with well-defined structures and molecular composition. We discussed how evolutionary pressures might emerge to affect these interactions, leading to the prevalence of low complexity domains in IDP sequences. These domains suppress spurious interactions and facilitate the formation of biologically meaningful condensates.

Dynamical control enables the formation of demixed biomolecular condensates

Cellular matter can be organized into compositionally distinct biomolecular condensates. For example, in Ashbya gossypii, the RNA-binding protein Whi3 forms distinct condensates with different RNA molecules. Using criteria derived from a physical framework for explaining how compositionally distinct condensates can form spontaneously via thermodynamic considerations, we find that condensates in vitro form mainly via heterotypic interactions in binary mixtures of Whi3 and RNA. However, within these condensates, RNA molecules become dynamically arrested. As a result, in ternary systems, simultaneous additions of Whi3 and pairs of distinct RNA molecules lead to well-mixed condensates, whereas delayed addition of an RNA component results in compositional distinctness. Therefore, compositional identities of condensates can be achieved via dynamical control, being driven, at least partially, by the dynamical arrest of RNA molecules. Finally, we show that synchronizing the production of different RNAs leads to more well-mixed, as opposed to compositionally distinct condensates in vivo.

Formation, function, and pathology of RNP granules

Ribonucleoprotein (RNP) granules are diverse membrane-less organelles that form through multivalent RNA-RNA, RNA-protein, and protein-protein interactions between RNPs. RNP granules are implicated in many aspects of RNA physiology, but in most cases their functions are poorly understood. RNP granules can be described through four key principles. First, RNP granules often arise because of the large size, high localized concentrations, and multivalent interactions of RNPs. Second, cells regulate RNP granule formation by multiple mechanisms including posttranslational modifications, protein chaperones, and RNA chaperones. Third, RNP granules impact cell physiology in multiple manners. Finally, dysregulation of RNP granules contributes to human diseases. Outstanding issues in the field remain, including determining the scale and molecular mechanisms of RNP granule function and how granule dysfunction contributes to human disease.

Quantitative real-time in-cell imaging reveals heterogeneous clusters of proteins prior to condensation

Our current understanding of biomolecular condensate formation is largely based on observing the final near-equilibrium condensate state. Despite expectations from classical nucleation theory, pre-critical protein clusters were recently shown to form under subsaturation conditions in vitro; if similar long-lived clusters comprising more than a few molecules are also present in cells, our understanding of the physical basis of biological phase separation may fundamentally change. Here, we combine fluorescence microscopy with photobleaching analysis to quantify the formation of clusters of NELF proteins in living, stressed cells. We categorise small and large clusters based on their dynamics and their response to p38 kinase inhibition. We find a broad distribution of pre-condensate cluster sizes and show that NELF protein cluster formation can be explained as non-classical nucleation with a surprisingly flat free-energy landscape for a wide range of sizes and an inhibition of condensation in unstressed cells.

Condensation Goes Viral: A Polymer Physics Perspective

The past decade has seen a revolution in our understanding of how the cellular environment is organized, where an incredible body of work has provided new insights into the role played by membraneless organelles. These rapid advancements have been made possible by an increasing awareness of the peculiar physical properties that give rise to such bodies and the complex biology that enables their function. Viral infections are not extraneous to this. Indeed, in host cells, viruses can harness existing membraneless compartments or, even, induce the formation of new ones. By hijacking the cellular machinery, these intracellular bodies can assist in the replication, assembly, and packaging of the viral genome as well as in the escape of the cellular immune response. Here, we provide a perspective on the fundamental polymer physics concepts that may help connect and interpret the different observed phenomena, ranging from the condensation of viral genomes to the phase separation of multicomponent solutions. We complement the discussion of the physical basis with a description of biophysical methods that can provide quantitative insights for testing and developing theoretical and computational models.

Theoretical and Data-Driven Approaches for Biomolecular Condensates

Biomolecular condensation processes are increasingly recognized as a fundamental mechanism that living cells use to organize biomolecules in time and space. These processes can lead to the formation of membraneless organelles that enable cells to perform distinct biochemical processes in controlled local environments, thereby supplying them with an additional degree of spatial control relative to that achieved by membrane-bound organelles. This fundamental importance of biomolecular condensation has motivated a quest to discover and understand the molecular mechanisms and determinants that drive and control this process. Within this molecular viewpoint, computational methods can provide a unique angle to studying biomolecular condensation processes by contributing the resolution and scale that are challenging to reach with experimental techniques alone. In this Review, we focus on three types of dry-lab approaches: theoretical methods, physics-driven simulations and data-driven machine learning methods. We review recent progress in using these tools for probing biomolecular condensation across all three fields and outline the key advantages and limitations of each of the approaches. We further discuss some of the key outstanding challenges that we foresee the community addressing next in order to develop a more complete picture of the molecular driving forces behind biomolecular condensation processes and their biological roles in health and disease.

Sequence variations of phase-separating proteins and resources for studying biomolecular condensates

Phase separation (PS) is an important mechanism underlying the formation of biomolecular condensates. Physiological condensates are associated with numerous biological processes, such as transcription, immunity, signaling, and synaptic transmission. Changes in particular amino acids or segments can disturb the protein's phase behavior and interactions with other biomolecules in condensates. It is thus presumed that variations in the phase-separating-prone domains can significantly impact the properties and functions of condensates. The dysfunction of condensates contributes to a number of pathological processes. Pharmacological perturbation of these condensates is proposed as a promising way to restore physiological states. In this review, we characterize the variations observed in PS proteins that lead to aberrant biomolecular compartmentalization. We also showcase recent advancements in bioinformatics of membraneless organelles (MLOs), focusing on available databases useful for screening PS proteins and describing endogenous condensates, guiding researchers to seek the underlying pathogenic mechanisms of biomolecular condensates.

Spatial and functional arrangement of Ebola virus polymerase inside phase-separated viral factories

Ebola virus (EBOV) infection induces the formation of membrane-less, cytoplasmic compartments termed viral factories, in which multiple viral proteins gather and coordinate viral transcription, replication, and assembly. Key to viral factory function is the recruitment of EBOV polymerase, a multifunctional machine that mediates transcription and replication of the viral RNA genome. We show that intracellularly reconstituted EBOV viral factories are biomolecular condensates, with composition-dependent internal exchange dynamics that likely facilitates viral replication. Within the viral factory, we found the EBOV polymerase clusters into foci. The distance between these foci increases when viral replication is enabled. In addition to the typical droplet-like viral factories, we report the formation of network-like viral factories during EBOV infection. Unlike droplet-like viral factories, network-like factories are inactive for EBOV nucleocapsid assembly. This unique view of EBOV propagation suggests a form-to-function relationship that describes how physical properties and internal structures of biomolecular condensates influence viral biogenesis.

Ribozyme activity modulates the physical properties of RNA-peptide coacervates

Condensed coacervate phases are now understood to be important features of modern cell biology, as well as valuable protocellular models in origin-of-life studies and synthetic biology. In each of these fields, the development of model systems with varied and tuneable material properties is of great importance for replicating properties of life. Here, we develop a ligase ribozyme system capable of concatenating short RNA fragments into long chains. Our results show that the formation of coacervate microdroplets with the ligase ribozyme and poly(L-lysine) enhances ribozyme rate and yield, which in turn increases the length of the anionic polymer component of the system and imparts specific physical properties to the droplets. Droplets containing active ribozyme sequences resist growth, do not wet or spread on unpassivated surfaces, and exhibit reduced transfer of RNA between droplets when compared to controls containing inactive sequences. These altered behaviours, which stem from RNA sequence and catalytic activity, constitute a specific phenotype and potential fitness advantage, opening the door to selection and evolution experiments based on a genotype-phenotype linkage.

Effects of linker length on phase separation: lessons from the Rubisco-EPYC1 system of the algal pyrenoid

Biomolecular condensates are membraneless organelles formed via phase separation of macromolecules, typically consisting of bond-forming "stickers" connected by flexible "linkers". Linkers have diverse roles, such as occupying space and facilitating interactions. To understand how linker length relative to other lengths affects condensation, we focus on the pyrenoid, which enhances photosynthesis in green algae. Specifically, we apply coarse-grained simulations and analytical theory to the pyrenoid proteins of Chlamydomonas reinhardtii: the rigid holoenzyme Rubisco and its flexible partner EPYC1. Remarkably, halving EPYC1 linker lengths decreases critical concentrations by ten-fold. We attribute this difference to the molecular "fit" between EPYC1 and Rubisco. Varying Rubisco sticker locations reveals that the native sites yield the poorest fit, thus optimizing phase separation. Surprisingly, shorter linkers mediate a transition to a gas of rods as Rubisco stickers approach the poles. These findings illustrate how intrinsically disordered proteins affect phase separation through the interplay of molecular length scales.

Intrinsic mesoscale properties of a Polycomb protein underpin heterochromatin fidelity

Little is understood about how the two major types of heterochromatin domains (HP1 and Polycomb) are kept separate. In the yeast Cryptococcus neoformans, the Polycomb-like protein Ccc1 prevents deposition of H3K27me3 at HP1 domains. Here we show that phase separation propensity underpins Ccc1 function. Mutations of the two basic clusters in the intrinsically disordered region or deletion of the coiled-coil dimerization domain alter phase separation behavior of Ccc1 in vitro and have commensurate effects on formation of Ccc1 condensates in vivo, which are enriched for PRC2. Notably, mutations that alter phase separation trigger ectopic H3K27me3 at HP1 domains. Supporting a direct condensate-driven mechanism for fidelity, Ccc1 droplets efficiently concentrate recombinant C. neoformans PRC2 in vitro whereas HP1 droplets do so only weakly. These studies establish a biochemical basis for chromatin regulation in which mesoscale biophysical properties play a key functional role.

Phosphorylation sites are evolutionary checkpoints against liquid-solid transition in protein condensates

Assemblies of multivalent RNA-binding protein fused in sarcoma (FUS) can exist in the functional liquid-like state as well as less dynamic and potentially toxic amyloid- and hydrogel-like states. How could then cells form liquid-like condensates while avoiding their transformation to amyloids? Here, we show how posttranslational phosphorylation can provide a "handle" that prevents liquid-solid transition of intracellular condensates containing FUS. Using residue-specific coarse-grained simulations, for 85 different mammalian FUS sequences, we show how the number of phosphorylation sites and their spatial arrangement affect intracluster dynamics preventing conversion to amyloids. All atom simulations further confirm that phosphorylation can effectively reduce the β-sheet propensity in amyloid-prone fragments of FUS. A detailed evolutionary analysis shows that mammalian FUS PLDs are enriched in amyloid-prone stretches compared to control neutrally evolved sequences, suggesting that mammalian FUS proteins evolved to self-assemble. However, in stark contrast to proteins that do not phase-separate for their function, mammalian sequences have phosphosites in close proximity to these amyloid-prone regions. These results suggest that evolution uses amyloid-prone sequences in prion-like domains to enhance phase separation of condensate proteins while enriching phosphorylation sites in close proximity to safeguard against liquid-solid transitions.

The maximum solubility product marks the threshold for condensation of multivalent biomolecules

Clustering of weakly interacting multivalent biomolecules underlies the formation of membraneless compartments known as condensates. As opposed to single-component (homotypic) systems, the concentration dependence of multicomponent (heterotypic) condensate formation is not well understood. We previously proposed the solubility product (SP), the product of monomer concentrations in the dilute phase, as a tool for understanding the concentration dependence of multicomponent systems. In this study, we further explore the limits of the SP concept using spatial Langevin dynamics and rule-based stochastic simulations. We show, for a variety of idealized molecular structures, how the maximum SP coincides with the onset of the phase transition, i.e., the formation of large clusters. We reveal the importance of intracluster binding in steering the free and cluster phase molecular distributions. We also show how structural features of biomolecules shape the SP profiles. The interplay of flexibility, length, and steric hindrance of linker regions controls the phase transition threshold. Remarkably, when SPs are normalized to nondimensional variables and plotted against the concentration scaled to the threshold for phase transition, the curves all coincide independent of the structural features of the binding partners. Similar coincidence is observed for the normalized clustering versus concentration plots. Overall, the principles derived from these systematic models will help guide and interpret in vitro and in vivo experiments on the biophysics of biomolecular condensates.

A conserved oligomerization domain in the disordered linker of coronavirus nucleocapsid proteins

The nucleocapsid (N-)protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a key role in viral assembly and scaffolding of the viral RNA. It promotes liquid-liquid phase separation (LLPS), forming dense droplets that support the assembly of ribonucleoprotein particles with as-of-yet unknown macromolecular architecture. Combining biophysical experiments, molecular dynamics simulations, and analysis of the mutational landscape, we describe a heretofore unknown oligomerization site that contributes to LLPS, is required for the assembly of higher-order protein-nucleic acid complexes, and is coupled to large-scale conformational changes of N-protein upon nucleic acid binding. The self-association interface is located in a leucine-rich sequence of the intrinsically disordered linker between N-protein folded domains and formed by transient helices assembling into trimeric coiled-coils. Critical residues stabilizing hydrophobic and electrostatic interactions between adjacent helices are highly protected against mutations in viable SARS-CoV-2 genomes, and the oligomerization motif is conserved across related coronaviruses, thus presenting a target for antiviral therapeutics.

Developments in describing equilibrium phase transitions of multivalent associative macromolecules

Biomolecular condensates are distinct cellular bodies that form and dissolve reversibly to organize cellular matter and biochemical reactions in space and time. Condensates are thought to form and dissolve under the influence of spontaneous and driven phase transitions of multivalent associative macromolecules. These include phase separation, which is defined by segregation of macromolecules from the solvent or from one another, and percolation or gelation, which is an inclusive networking transition driven by reversible associations among multivalent macromolecules. Considerable progress has been made to model sequence-specific phase transitions, especially for intrinsically disordered proteins. Here, we summarize the state-of-the-art of theories and computations aimed at understanding and modeling sequence-specific, thermodynamically controlled, coupled associative and segregative phase transitions of archetypal multivalent macromolecules.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

Dynamical control enables the formation of demixed biomolecular condensates

Macromolecular phase separation underlies the regulated formation and dissolution of biomolecular condensates. What is unclear is how condensates of distinct and shared macromolecular compositions form and coexist within cellular milieus. Here, we use theory and computation to establish thermodynamic criteria that must be satisfied to achieve compositionally distinct condensates. We applied these criteria to an archetypal ribonucleoprotein condensate and discovered that demixing into distinct protein-RNA condensates cannot be the result of purely thermodynamic considerations. Instead, demixed, compositionally distinct condensates arise due to asynchronies in timescales that emerge from differences in long-lived protein-RNA and RNA-RNA crosslinks. This type of dynamical control is also found to be active in live cells whereby asynchronous production of molecules is required for realizing demixed protein-RNA condensates. We find that interactions that exert dynamical control provide a versatile and generalizable way to influence the compositions of coexisting condensates in live cells.

Topological Considerations in Biomolecular Condensation

Biomolecular condensation and phase separation are increasingly understood to play crucial roles in cellular compartmentalization and spatiotemporal regulation of cell machinery implicated in function and pathology. A key aspect of current research is to gain insight into the underlying physical mechanisms of these processes. Accordingly, concepts of soft matter and polymer physics, the thermodynamics of mixing, and material science have been utilized for understanding condensation mechanisms of multivalent macromolecules resulting in viscoelastic mesoscopic supramolecular assemblies. Here, we focus on two topological concepts that have recently been providing key mechanistic understanding in the field. First, we will discuss how percolation provides a network-topology-related framework that offers an interesting paradigm to understand the complex networking of dense 'connected' condensate structures and, therefore, their phase behavior. Second, we will discuss the idea of entanglement as another topological concept that has deep roots in polymer physics and important implications for biomolecular condensates. We will first review some historical developments and fundamentals of these concepts, then we will discuss current advancements and recent examples. Our discussion ends with a few open questions and the challenges to address them, hinting at unveiling fresh possibilities for the modification of existing knowledge as well as the development of new concepts relevant to condensate science.

Dynamical control enables the formation of demixed biomolecular condensates

Macromolecular phase separation underlies the regulated formation and dissolution of biomolecular condensates. What is unclear is how condensates of distinct and shared macromolecular compositions form and coexist within cellular milieus. Here, we use theory and computation to establish thermodynamic criteria that must be satisfied to achieve compositionally distinct condensates. We applied these criteria to an archetypal ribonucleoprotein condensate and discovered that demixing into distinct protein-RNA condensates cannot be the result of purely thermodynamic considerations. Instead, demixed, compositionally distinct condensates arise due to asynchronies in timescales that emerge from differences in long-lived protein-RNA and RNA-RNA crosslinks. This type of dynamical control is also found to be active in live cells whereby asynchronous production of molecules is required for realizing demixed protein-RNA condensates. We find that interactions that exert dynamical control provide a versatile and generalizable way to influence the compositions of coexisting condensates in live cells.

Formation of dissipative structures in microscopic models of mixtures with species interconversion

The separation of substances into different phases is ubiquitous in nature and important scientifically and technologically. This phenomenon may become drastically different if the species involved, whether molecules or supramolecular assemblies, interconvert. In the presence of an external force large enough to overcome energetic differences between the interconvertible species (forced interconversion), the two alternative species will be present in equal amounts, and the striking phenomenon of steady-state, restricted phase separation into mesoscales is observed. Such microphase separation is one of the simplest examples of dissipative structures in condensed matter. In this work, we investigate the formation of such mesoscale steady-state structures through Monte Carlo and molecular dynamics simulations of three physically distinct microscopic models of binary mixtures that exhibit both equilibrium (natural) interconversion and a nonequilibrium source of forced interconversion. We show that this source can be introduced through an internal imbalance of intermolecular forces or an external flux of energy that promotes molecular interconversion, possible manifestations of which could include the internal nonequilibrium environment of living cells or a flux of photons. The main trends and observations from the simulations are well captured by a nonequilibrium thermodynamic theory of phase transitions affected by interconversion. We show how a nonequilibrium bicontinuous microemulsion or a spatially modulated state may be generated depending on the interplay between diffusion, natural interconversion, and forced interconversion.

An Introduction to the Stickers-and-Spacers Framework as Applied to Biomolecular Condensates

Cellular organization is determined by a combination of membrane-bound and membrane-less biomolecular assemblies that range from clusters of tens of molecules to micrometer-sized cellular bodies. Over the last decade, membrane-less assemblies have come to be referred to as biomolecular condensates, reflecting their ability to condense specific molecules with respect to the remainder of the cell. In many cases, the physics of phase transitions provides a conceptual framework and a mathematical toolkit to describe the assembly, maintenance, and dissolution of biomolecular condensates. Among the various quantitative and qualitative models applied to understand intracellular phase transitions, the stickers-and-spacers framework offers an intuitive yet rigorous means to map biomolecular sequences and structure to the driving forces needed for higher-order assembly. This chapter introduces the fundamental concepts behind the stickers-and-spacers model, considers its application to different biological systems, and discusses limitations and misconceptions around the model.

Size distributions of intracellular condensates reflect competition between coalescence and nucleation

Phase separation of biomolecules into condensates has emerged as a mechanism for intracellular organization and affects many intracellular processes, including reaction pathways through the clustering of enzymes and pathway intermediates. Precise and rapid spatiotemporal control of reactions by condensates requires tuning of their sizes. However, the physical processes that govern the distribution of condensate sizes remain unclear. Here we show that both native and synthetic condensates display an exponential size distribution, which is captured by Monte Carlo simulations of fast nucleation followed by coalescence. In contrast, pathological aggregates exhibit a power-law size distribution. These distinct behaviours reflect the relative importance of nucleation and coalescence kinetics. We demonstrate this by utilizing a combination of synthetic and native condensates to probe the underlying physical mechanisms determining condensate size. The appearance of exponential distributions for abrupt nucleation versus power-law distributions under continuous nucleation may reflect a general principle that determines condensate size distributions.

Molecular Determinants for the Layering and Coarsening of Biological Condensates

Many membraneless organelles, or biological condensates, form through phase separation, and play key roles in signal sensing and transcriptional regulation. While the functional importance of these condensates has inspired many studies to characterize their stability and spatial organization, the underlying principles that dictate these emergent properties are still being uncovered. In this review, we examine recent work on biological condensates, especially multicomponent systems. We focus on connecting molecular factors such as binding energy, valency, and stoichiometry with the interfacial tension, explaining the nontrivial interior organization in many condensates. We further discuss mechanisms that arrest condensate coalescence by lowering the surface tension or introducing kinetic barriers to stabilize the multidroplet state.

Short PolyA RNA Homopolymers Undergo Mg2+-Mediated Kinetically Arrested Condensation

RNA-RNA interactions have increasingly been recognized for their potential to shape the mesoscale properties of biomolecular condensates, influencing morphology, organization, and material state through networking interactions. While most studies have focused on networking via Watson-Crick base pairing interactions, previous work has suggested a potential for noncanonical RNA-RNA interactions to also give rise to condensation and alter overall material state. Here, we test the phase separation of short polyA RNA (polyrA) homopolymers. We discover and characterize the potential for short polyrA sequences to form RNA condensates at lower Mg²⁺ concentrations than previously observed, which appear as internally arrested droplets with slow polyrA diffusion despite continued fusion. Our work also reveals a negative cooperativity effect between the effects of Mg²⁺ and Na⁺ on polyrA condensation. Finally, we observe that polyrA sequences can act as promoters of phase separation in mixed sequences. These results demonstrate the potential for noncanonical interactions to act as networking stickers, leading to specific condensation properties inherent to polyrA composition and structure, with implications for the fundamental physical chemistry of the system and function of polyA RNA in biology.

Principles Governing the Phase Separation of Multidomain Proteins

A variety of membraneless organelles, often termed "biological condensates", play an important role in the regulation of cellular processes such as gene transcription, translation, and protein quality control. On the basis of experimental and theoretical investigations, liquid-liquid phase separation (LLPS) has been proposed as a possible mechanism for the origin of biological condensates. LLPS requires multivalent macromolecules that template the formation of long-range, intermolecular interaction networks and results in the formation of condensates with defined composition and material properties. Multivalent interactions driving LLPS exhibit a wide range of modes from highly stereospecific to nonspecific and involve both folded and disordered regions. Multidomain proteins serve as suitable macromolecules for promoting phase separation and achieving disparate functions due to their potential for multivalent interactions and regulation. Here, we aim to highlight the influence of the domain architecture and interdomain interactions on the phase separation of multidomain protein condensates. First, the general principles underlying these interactions are illustrated on the basis of examples of multidomain proteins that are predominantly associated with nucleic acid binding and protein quality control and contain both folded and disordered regions. Next, the examples showcase how LLPS properties of folded and disordered regions can be leveraged to engineer multidomain constructs that form condensates with the desired assembly and functional properties. Finally, we highlight the need for improvements in coarse-grained computational models that can provide molecular-level insights into multidomain protein condensates in conjunction with experimental efforts.

Different states and the associated fates of biomolecular condensates

Biomolecular condensates are functional assemblies, which can enrich intrinsically disordered proteins (IDPs) and/or RNAs at concentrations that are orders of magnitude higher than the bulk. In their native functional state, these structures can exist in multiple physical states including liquid-droplet phase, hydrogels, and solid assemblies. On the other hand, an aberrant transition between these physical states can result in loss-of-function or a gain-of-toxic-function. A prime example of such an aberrant transition is droplet aging—a phenomenon where some condensates may progressively transition into less dynamic material states at biologically relevant timescales. In this essay, we review structural and viscoelastic roots of aberrant liquid–solid transitions. Also, we highlight the different checkpoints and experimentally tunable handles, both active (ATP-dependent enzymes, post-translational modifications) and passive (colocalization of RNA molecules), that could alter the material state of assemblies.

Kinetic frustration by limited bond availability controls the LAT protein condensation phase transition on membranes

LAT is a membrane-linked scaffold protein that undergoes a phase transition to form a two-dimensional protein condensate on the membrane during T cell activation. Governed by tyrosine phosphorylation, LAT recruits various proteins that ultimately enable condensation through a percolation network of discrete and selective protein-protein interactions. Here, we describe detailed kinetic measurements of the phase transition, along with coarse-grained model simulations, that reveal that LAT condensation is kinetically frustrated by the availability of bonds to form the network. Unlike typical miscibility transitions in which compact domains may coexist at equilibrium, the LAT condensates are dynamically arrested in extended states, kinetically trapped out of equilibrium. Modeling identifies the structural basis for this kinetic arrest as the formation of spindle arrangements, favored by limited multivalent binding interactions along the flexible, intrinsically disordered LAT protein. These results reveal how local factors controlling the kinetics of LAT condensation enable formation of different, stable condensates, which may ultimately coexist within the cell.

Dynamics of protein condensates in weak-binding regime

Weak complementary interactions between proteins and nucleic acids are the main driving forces of intracellular liquid-liquid phase separation. The sticker-spacer model has emerged as a unifying principle for understanding the phase behavior of these multivalent molecules. It remains elusive how specific interactions mediated by stickers contribute to the rheological properties of the liquid condensates. Previous studies have revealed that for strong binding strength ɛ_{b}, the bulk diffusivity D depends on the effective bond lifetime τ, viz., D∝τ^{-1}. Consequently, equal concentrations of the complementary stickers induce a slow down in the dynamics of the condensates D∝e^{-1.5ɛ_{b}}. However, for weak-binding strength, it is expected that the resulting condensates are dynamic, loose network liquids rather than kinetically arrested, compact clusters. We develop a mean-field theory using the thermodynamics of the associative polymers and perform molecular-dynamics simulations based on the sticker-spacer model to study the controlling factors in the structure and dynamics of such condensates in the weak-binding regime. Through scaling analysis, we delineate how the free sticker fraction W_{f} and the bulk diffusivity D decrease with increasing binding energy and find that the internal dynamics of such network liquids are controlled by the free sticker fraction D∝W_{f}∝e^{-0.5ɛ_{b}} rather than the effective bond lifetime. Referred to as the free-sticker-dominated diffusivity, the microscopic slowdown due to a gradual loss of the free stickers affects the viscosity of the condensates as well, with the scaling of the zero-shear viscosity η∝e^{0.5ɛ_{b}}. Therefore, the way of controlling the structure, diffusivity, and viscosity of the condensates through the binding energy can be tested experimentally.

Micellization: A new principle in the formation of biomolecular condensates

Phase separation is a fundamental mechanism for compartmentalization in cells and leads to the formation of biomolecular condensates, generally containing various RNA molecules. RNAs are biomolecules that can serve as suitable scaffolds for biomolecular condensates and determine their forms and functions. Many studies have focused on biomolecular condensates formed by liquid-liquid phase separation (LLPS), one type of intracellular phase separation mechanism. We recently identified that paraspeckle nuclear bodies use an intracellular phase separation mechanism called micellization of block copolymers in their formation. The paraspeckles are scaffolded by NEAT1_2 long non-coding RNAs (lncRNAs) and their partner RNA-binding proteins (NEAT1_2 RNA-protein complexes [RNPs]). The NEAT1_2 RNPs act as block copolymers and the paraspeckles assemble through micellization. In LLPS, condensates grow without bound as long as components are available and typically have spherical shapes to minimize surface tension. In contrast, the size, shape, and internal morphology of the condensates are more strictly controlled in micellization. Here, we discuss the potential importance and future perspectives of micellization of block copolymers of RNPs in cells, including the construction of designer condensates with optimal internal organization, shape, and size according to design guidelines of block copolymers.

Molecular mechanism for the synchronized electrostatic coacervation and co-aggregation of alpha-synuclein and tau

Amyloid aggregation of α-synuclein (αS) is the hallmark of Parkinson’s disease and other synucleinopathies. Recently, Tau protein, generally associated with Alzheimer’s disease, has been linked to αS pathology and observed to co-localize in αS-rich disease inclusions, although the molecular mechanisms for the co-aggregation of both proteins remain elusive. We report here that αS phase-separates into liquid condensates by electrostatic complex coacervation with positively charged polypeptides such as Tau. Condensates undergo either fast gelation or coalescence followed by slow amyloid aggregation depending on the affinity of αS for the poly-cation and the rate of valence exhaustion of the condensate network. By combining a set of advanced biophysical techniques, we have been able to characterize αS/Tau liquid-liquid phase separation and identified key factors that lead to the formation of hetero-aggregates containing both proteins in the interior of the liquid protein condensates.

Here, the authors report that α-synuclein phase-separates into liquid condensates with positively charged polypeptides such as Tau. The condensates undergo different maturation processes, including the formation of α-synuclein/Tau amyloid hetero-aggregates inside the condensates.

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-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.

Multivalent interactions with RNA drive recruitment and dynamics in biomolecular condensates in Xenopus oocytes

RNA localization and biomolecular condensate formation are key biological strategies for organizing the cytoplasm and generating cellular polarity. In Xenopus oocytes, RNAs required for germ layer patterning localize in biomolecular condensates, termed Localization bodies (L-bodies). Here, we have used an L-body RNA-binding protein, PTBP3, to test the role of RNA–protein interactions in regulating the biophysical characteristics of L-bodies in vivo and PTBP3–RNA condensates in vitro. Our results reveal that RNA–protein interactions drive recruitment of PTBP3 and localized RNA to L-bodies and that multivalent interactions tune the dynamics of the PTBP3 after localization. In a concentration-dependent manner, RNA becomes non-dynamic and interactions with the RNA determine PTBP3 dynamics within these biomolecular condensates in vivo and in vitro. Importantly, RNA, and not protein, is required for maintenance of the PTBP3–RNA condensates in vitro, pointing to a model where RNA serves as a non-dynamic substructure in these condensates.

•

RNA–protein interactions drive recruitment of both RNA and protein to L-bodies

•

RNA is non-dynamic in both L-bodies and in vitro condensates

•

Multivalent interactions with RNA tune protein dynamics both in vivo and in vitro

•

RNA, but not protein, is required for maintenance of the in vitro condensates

Aging can transform single-component protein condensates into multiphase architectures

Phase-separated biomolecular condensates that contain multiple coexisting phases are widespread in vitro and in cells. Multiphase condensates emerge readily within multicomponent mixtures of biomolecules (e.g., proteins and nucleic acids) when the different components present sufficient physicochemical diversity (e.g., in intermolecular forces, structure, and chemical composition) to sustain separate coexisting phases. Because such diversity is highly coupled to the solution conditions (e.g., temperature, pH, salt, composition), it can manifest itself immediately from the nucleation and growth stages of condensate formation, develop spontaneously due to external stimuli or emerge progressively as the condensates age. Here, we investigate thermodynamic factors that can explain the progressive intrinsic transformation of single-component condensates into multiphase architectures during the nonequilibrium process of aging. We develop a multiscale model that integrates atomistic simulations of proteins, sequence-dependent coarse-grained simulations of condensates, and a minimal model of dynamically aging condensates with nonconservative intermolecular forces. Our nonequilibrium simulations of condensate aging predict that single-component condensates that are initially homogeneous and liquid like can transform into gel-core/liquid-shell or liquid-core/gel-shell multiphase condensates as they age due to gradual and irreversible enhancement of interprotein interactions. The type of multiphase architecture is determined by the aging mechanism, the molecular organization of the gel and liquid phases, and the chemical makeup of the protein. Notably, we predict that interprotein disorder to order transitions within the prion-like domains of intracellular proteins can lead to the required nonconservative enhancement of intermolecular interactions. Our study, therefore, predicts a potential mechanism by which the nonequilibrium process of aging results in single-component multiphase condensates.

Controlling cluster size in 2D phase-separating binary mixtures with specific interactions

By varying the concentration of molecules in the cytoplasm or on the membrane, cells can induce the formation of condensates and liquid droplets, similar to phase separation. Their thermodynamics, much studied, depends on the mutual interactions between microscopic constituents. Here, we focus on the kinetics and size control of 2D clusters, forming on membranes. Using molecular dynamics of patchy colloids, we model a system of two species of proteins, giving origin to specific heterotypic bonds. We find that concentrations, together with valence and bond strength, control both the size and the growth time rate of the clusters. In particular, if one species is in large excess, it gradually saturates the binding sites of the other species; the system then becomes kinetically arrested and cluster coarsening slows down or stops, thus yielding effective size selection. This phenomenology is observed both in solid and fluid clusters, which feature additional generic homotypic interactions and are reminiscent of the ones observed on biological membranes.

Post-translational modifications in liquid-liquid phase separation: a comprehensive review

Liquid-liquid phase separation (LLPS) has received significant attention in recent biological studies. It refers to a phenomenon that biomolecule exceeds the solubility, condensates and separates itself from solution in liquid like droplets formation. Our understanding of it has also changed from memebraneless organelles to compartmentalization, muti-functional crucibles, and reaction regulators. Although this phenomenon has been employed for a variety of biological processes, recent studies mainly focus on its physiological significance, and the comprehensive research of the underlying physical mechanism is limited. The characteristics of side chains of amino acids and the interaction tendency of proteins function importantly in regulating LLPS thus should be pay more attention on. In addition, the importance of post-translational modifications (PTMs) has been underestimated, despite their abundance and crucial functions in maintaining the electrostatic balance. In this review, we first introduce the driving forces and protein secondary structures involved in LLPS and their different physical functions in cell life processes. Subsequently, we summarize the existing reports on PTM regulation related to LLPS and analyze the underlying basic principles, hoping to find some common relations between LLPS and PTM. Finally, we speculate several unreported PTMs that may have a significant impact on phase separation basing on the findings.

On the stability and layered organization of protein-DNA condensates

Multi-component phase separation is emerging as a key mechanism for the formation of biological condensates that play essential roles in signal sensing and transcriptional regulation. The molecular factors that dictate these condensates' stability and spatial organization are not fully understood, and it remains challenging to predict their microstructures. Using a near-atomistic, chemically accurate force field, we studied the phase behavior of chromatin regulators that are crucial for heterochromatin organization and their interactions with DNA. Our computed phase diagrams recapitulated previous experimental findings on different proteins. They revealed a strong dependence of condensate stability on the protein-DNA mixing ratio as a result of balancing protein-protein interactions and charge neutralization. Notably, a layered organization was observed in condensates formed by mixing HP1, histone H1, and DNA. This layered organization may be of biological relevance, as it enables cooperative DNA packaging between the two chromatin regulators: histone H1 softens the DNA to facilitate the compaction induced by HP1 droplets. Our study supports near-atomistic models as a valuable tool for characterizing the structure and stability of biological condensates.