Protein Network Structure Enables Switching between Liquid and Gel States

Fundamental Aspects of Phase-Separated Biomolecular Condensates

Biomolecular condensates, formed through phase separation, are upending our understanding in much of molecular, cell, and developmental biology. There is an urgent need to elucidate the physicochemical foundations of the behaviors and properties of biomolecular condensates. Here we aim to fill this need by writing a comprehensive, critical, and accessible review on the fundamental aspects of phase-separated biomolecular condensates. We introduce the relevant theoretical background, present the theoretical basis for the computation and experimental measurement of condensate properties, and give mechanistic interpretations of condensate behaviors and properties in terms of interactions at the molecular and residue levels.

Crossing boundaries of light microscopy resolution discerns novel assemblies in the nucleolus

The nucleolus is the largest membraneless organelle and nuclear body in mammalian cells. It is primarily involved in the biogenesis of ribosomes, essential macromolecular machines responsible for synthesizing all proteins required by the cell. The assembly of ribosomes is evolutionarily conserved and accounts for the most energy-consuming cellular process needed for cell growth, proliferation, and homeostasis. Despite the significance of this process, the substructural mechanistic principles of the nucleolar function in preribosome biogenesis have only recently begun to emerge. Here, we provide a new perspective using advanced super-resolution microscopy and single-molecule MINFLUX nanoscopy on the mechanistic principles governing ribosomal RNA-seeded nucleolar formation and the resulting tripartite suborganization of the nucleolus driven, in part, by liquid-liquid phase separation. With recent advances in the cryogenic electron microscopy (cryoEM) structural analysis of ribosome biogenesis intermediates, we highlight the current understanding of the step-wise assembly of preribosomal subunits in the nucleolus. Finally, we address how novel anticancer drug candidates target early steps in ribosome biogenesis to exploit these essential dependencies for growth arrest and tumor control.

Determinants that enable disordered protein assembly into discrete condensed phases

Cells harbour numerous mesoscale membraneless compartments that house specific biochemical processes and perform distinct cellular functions. These protein- and RNA-rich bodies are thought to form through multivalent interactions among proteins and nucleic acids, resulting in demixing via liquid-liquid phase separation. Proteins harbouring intrinsically disordered regions (IDRs) predominate in membraneless organelles. However, it is not known whether IDR sequence alone can dictate the formation of distinct condensed phases. We identified a pair of IDRs capable of forming spatially distinct condensates when expressed in cells. When reconstituted in vitro, these model proteins do not co-partition, suggesting condensation specificity is encoded directly in the polypeptide sequences. Through computational modelling and mutagenesis, we identified the amino acids and chain properties governing homotypic and heterotypic interactions that direct selective condensation. These results form the basis of physicochemical principles that may direct subcellular organization of IDRs into specific condensates and reveal an IDR code that can guide construction of orthogonal membraneless compartments.

Biomolecular condensates form spatially inhomogeneous network fluids

The functions of biomolecular condensates are thought to be influenced by their material properties, and these will be determined by the internal organization of molecules within condensates. However, structural characterizations of condensates are challenging, and rarely reported. Here, we deploy a combination of small angle neutron scattering, fluorescence recovery after photobleaching, and coarse-grained molecular dynamics simulations to provide structural descriptions of model condensates that are formed by macromolecules from nucleolar granular components (GCs). We show that these minimal facsimiles of GCs form condensates that are network fluids featuring spatial inhomogeneities across different length scales that reflect the contributions of distinct protein and peptide domains. The network-like inhomogeneous organization is characterized by a coexistence of liquid- and gas-like macromolecular densities that engenders bimodality of internal molecular dynamics. These insights suggest that condensates formed by multivalent proteins share features with network fluids formed by systems such as patchy or hairy colloids.

Biomolecular condensates form spatially inhomogeneous network fluids

The functions of biomolecular condensates are thought to be influenced by their material properties, and these will be determined by the internal organization of molecules within condensates. However, structural characterizations of condensates are challenging, and rarely reported. Here, we deploy a combination of small angle neutron scattering, fluorescence recovery after photobleaching, and coarse-grained molecular dynamics simulations to provide structural descriptions of model condensates that are formed by macromolecules from nucleolar granular components (GCs). We show that these minimal facsimiles of GCs form condensates that are network fluids featuring spatial inhomogeneities across different length scales that reflect the contributions of distinct protein and peptide domains. The network-like inhomogeneous organization is characterized by a coexistence of liquid- and gas-like macromolecular densities that engenders bimodality of internal molecular dynamics. These insights suggest that condensates formed by multivalent proteins share features with network fluids formed by systems such as patchy or hairy colloids.

Surface Charge Can Modulate Phase Separation of Multidomain Proteins

Phase separation has emerged as an important mechanism explaining the formation of certain biomolecular condensates. Biological phase separation is often driven by the multivalent interactions of modular protein domains. Beyond valency, the physical features of folded domains that promote phase separation are poorly understood. We used a model system─the small ubiquitin modifier (SUMO) and its peptide ligand, the SUMO interaction motif (SIM)─to examine how domain surface charge influences multivalency-driven phase separation. Phase separation of polySUMO and polySIM was altered by pH via a change in the protonation state of SUMO surface histidines. These effects were recapitulated by histidine mutations, which modulated SUMO solubility and polySUMO-polySIM phase separation in parallel and were quantitatively explained by atomistic modeling of weak interactions among proteins in the system. Thus, surface charge can tune the phase separation of multivalent proteins, suggesting a means of controlling phase separation biologically, evolutionarily, and therapeutically.

Toward Accurate Simulation of Coupling between Protein Secondary Structure and Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate phase separation that underlies the formation of a biomolecular condensate. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding the sequence-specific phase separation of IDPs. However, the widely used Cα-only models are limited in capturing the peptide nature of IDPs, particularly backbone-mediated interactions and effects of secondary structures, in phase separation. Here, we describe a hybrid resolution (HyRes) protein model toward a more accurate description of the backbone and transient secondary structures in phase separation. With an atomistic backbone and coarse-grained side chains, HyRes can semiquantitatively capture the residue helical propensity and overall chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for the direct simulation of spontaneous phase separation and, at the same time, appears accurate enough to resolve the effects of single His to Lys mutations. HyRes simulations also successfully predict increased β-structure formation in the condensate, consistent with available experimental CD data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate the phase separation propensity as measured by the saturation concentration. The simulations successfully recapitulate the effect of these mutants on the helicity and phase separation propensity of TDP-43 CR. Analyses reveal that the balance between backbone and side chain-mediated interactions, but not helicity itself, actually determines phase separation propensity. These results support that HyRes represents an effective protein model for molecular simulation of IDP phase separation and will help to elucidate the coupling between transient secondary structures and phase separation.

Reduction of oligomer size modulates the competition between cluster formation and phase separation of the tumor suppressor SPOP

Phase separation compartmentalizes many cellular pathways. Given that the same interactions that drive phase separation mediate the formation of soluble complexes below the saturation concentration, the contribution of condensates versus complexes to function is sometimes unclear. Here, we characterized several new cancer-associated mutations of the tumor suppressor speckle-type POZ protein (SPOP), a substrate recognition subunit of the Cullin3-RING ubiquitin ligase. This pointed to a strategy for generating separation-of-function mutations. SPOP self-associates into linear oligomers and interacts with multivalent substrates, and this mediates the formation of condensates. These condensates bear the hallmarks of enzymatic ubiquitination activity. We characterized the effect of mutations in the dimerization domains of SPOP on its linear oligomerization, binding to its substrate DAXX, and phase separation with DAXX. We showed that the mutations reduce SPOP oligomerization and shift the size distribution of SPOP oligomers to smaller sizes. The mutations therefore reduce the binding affinity to DAXX but unexpectedly enhance the poly-ubiquitination activity of SPOP toward DAXX. Enhanced activity may be explained by enhanced phase separation of DAXX with the SPOP mutants. Our results provide a comparative assessment of the functional role of complexes versus condensates and support a model in which phase separation is an important factor in SPOP function. Our findings also suggest that tuning of linear SPOP self-association could be used by the cell to modulate activity and provide insights into the mechanisms underlying hypermorphic SPOP mutations. The characteristics of cancer-associated SPOP mutations suggest a route for designing separation-of-function mutations in other phase-separating systems.

Conformational entropy in molecular recognition of intrinsically disordered proteins

Broad conformational ensembles make intrinsically disordered proteins or regions entropically intriguing. Although methodologically challenging and understudied, emerging studies into their changes in conformational entropy (ΔS°conf) upon complex formation have provided both quantitative and qualitative insight. Recent work based on thermodynamics from isothermal titration calorimetry and NMR spectroscopy uncovers an expanded repertoire of regulatory mechanisms, where ΔS°conf plays roles in partner selection, state behavior, functional buffering, allosteric regulation, and drug design. We highlight these mechanisms to display the large entropic reservoir of IDPs for the regulation of molecular communication. We call upon the field to make efforts to contribute to this insight as more studies are needed for forwarding mechanistic decoding of intrinsically disordered proteins and their complexes.

Biomolecular condensates form spatially inhomogeneous network fluids

The functions of biomolecular condensates are thought to be influenced by their material properties, and these are in turn determined by the multiscale structural features within condensates. However, structural characterizations of condensates are challenging, and hence rarely reported. Here, we deploy a combination of small angle neutron scattering, fluorescence recovery after photobleaching, and bespoke coarse-grained molecular dynamics simulations to provide structural descriptions of model condensates that mimic nucleolar granular components (GCs). We show that facsimiles of GCs are network fluids featuring spatial inhomogeneities across hierarchies of length scales that reflect the contributions of distinct protein and peptide domains. The network-like inhomogeneous organization is characterized by a coexistence of liquid- and gas-like macromolecular densities that engenders bimodality of internal molecular dynamics. These insights, extracted from a combination of approaches, suggest that condensates formed by multivalent proteins share features with network fluids formed by associative systems such as patchy or hairy colloids.

Polyubiquitin ligand-induced phase transitions are optimized by spacing between ubiquitin units

Biomolecular condensates form via multivalent interactions among key macromolecules and are regulated through ligand binding and/or posttranslational modifications. One such modification is ubiquitination, the covalent addition of ubiquitin (Ub) or polyubiquitin chains to target macromolecules. Specific interactions between polyubiquitin chains and partner proteins, including hHR23B, NEMO, and UBQLN2, regulate condensate assembly or disassembly. Here, we used a library of designed polyubiquitin hubs and UBQLN2 as model systems for determining the driving forces of ligand-mediated phase transitions. Perturbations to either the UBQLN2-binding surface of Ub or the spacing between Ub units reduce the ability of hubs to modulate UBQLN2 phase behavior. By developing an analytical model based on polyphasic linkage principles that accurately described the effects of different hubs on UBQLN2 phase separation, we determined that introduction of Ub to UBQLN2 condensates incurs a significant inclusion energetic penalty. This penalty antagonizes the ability of polyUb hubs to scaffold multiple UBQLN2 molecules and cooperatively amplify phase separation. The extent to which polyubiquitin hubs promote UBQLN2 phase separation is encoded in the spacings between Ub units. This spacing is modulated by chains of different linkages and designed chains of different architectures, thus illustrating how the ubiquitin code regulates functionality via the emergent properties of the condensate. The spacing in naturally occurring linear polyubiquitin chains is already optimized to promote phase separation with UBQLN2. We expect our findings to extend to other condensates, emphasizing the importance of ligand properties, including concentration, valency, affinity, and spacing between binding sites in studies and designs of condensates.

Determinants of Disordered Protein Co-Assembly Into Discrete Condensed Phases

Cells harbor numerous mesoscale membraneless compartments that house specific biochemical processes and perform distinct cellular functions. These protein and RNA-rich bodies are thought to form through multivalent interactions among proteins and nucleic acids resulting in demixing via liquid-liquid phase separation (LLPS). Proteins harboring intrinsically disordered regions (IDRs) predominate in membraneless organelles. However, it is not known whether IDR sequence alone can dictate the formation of distinct condensed phases. We identified a pair of IDRs capable of forming spatially distinct condensates when expressed in cells. When reconstituted in vitro, these model proteins do not co-partition, suggesting condensation specificity is encoded directly in the polypeptide sequences. Through computational modeling and mutagenesis, we identified the amino acids and chain properties governing homotypic and heterotypic interactions that direct selective condensation. These results form the basis of physicochemical principles that may direct subcellular organization of IDRs into specific condensates and reveal an IDR code that can guide construction of orthogonal membraneless compartments.

Accurate Simulation of Coupling between Protein Secondary Structure and Liquid-Liquid Phase Separation

Intrinsically disordered proteins (IDPs) frequently mediate liquid-liquid phase separation (LLPS) that underlies the formation of membraneless organelles. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding sequence-specific phase separation of IDPs. However, the widely-used Cα-only models are severely limited in capturing the peptide nature of IDPs, including backbone-mediated interactions and effects of secondary structures, in LLPS. Here, we describe a hybrid resolution (HyRes) protein model for accurate description of the backbone and transient secondary structures in LLPS. With an atomistic backbone and coarse-grained side chains, HyRes accurately predicts the residue helical propensity and chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for direct simulation of spontaneous phase separation, and at the same time accurate enough to resolve the effects of single mutations. HyRes simulations also successfully predict increased beta-sheet formation in the condensate, consistent with available experimental data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate LLPS propensity. The simulations successfully recapitulate the effect of these mutants on the helicity and LLPS propensity of TDP-43 CR. Analyses reveal that the balance between backbone and sidechain-mediated interactions, but not helicity itself, actually determines LLPS propensity. We believe that the HyRes model represents an important advance in the molecular simulation of LLPS and will help elucidate the coupling between IDP transient secondary structures and phase separation.

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.

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.

Protein compactness and interaction valency define the architecture of a biomolecular condensate across scales

Non-membrane-bound biomolecular condensates have been proposed to represent an important mode of subcellular organization in diverse biological settings. However, the fundamental principles governing the spatial organization and dynamics of condensates at the atomistic level remain unclear. The Saccharomyces cerevisiae Lge1 protein is required for histone H2B ubiquitination and its N-terminal intrinsically disordered fragment (Lge11-80) undergoes robust phase separation. This study connects single- and multi-chain all-atom molecular dynamics simulations of Lge11-80 with the in vitro behavior of Lge11-80 condensates. Analysis of modeled protein-protein interactions elucidates the key determinants of Lge11-80 condensate formation and links configurational entropy, valency, and compactness of proteins inside the condensates. A newly derived analytical formalism, related to colloid fractal cluster formation, describes condensate architecture across length scales as a function of protein valency and compactness. In particular, the formalism provides an atomistically resolved model of Lge11-80 condensates on the scale of hundreds of nanometers starting from individual protein conformers captured in simulations. The simulation-derived fractal dimensions of condensates of Lge11-80 and its mutants agree with their in vitro morphologies. The presented framework enables a multiscale description of biomolecular condensates and embeds their study in a wider context of colloid self-organization.

Decoding optimal ligand design for multicomponent condensates

Biomolecular condensates form via multivalent interactions among key macromolecules and are regulated through ligand binding and/or post-translational modifications. One such modification is ubiquitination, the covalent addition of ubiquitin (Ub) or polyubiquitin chains to target macromolecules for various cellular processes. Specific interactions between polyubiquitin chains and partner proteins, including hHR23B, NEMO, and UBQLN2, regulate condensate assembly or disassembly. Here, we used a library of designed polyubiquitin hubs and UBQLN2 as model systems for determining the driving forces of ligand-mediated phase transitions. Perturbations to the UBQLN2-binding surface of Ub or deviations from the optimal spacing between Ub units reduce the ability of hubs to modulate UBQLN2 phase behavior. By developing an analytical model that accurately described the effects of different hubs on UBQLN2 phase diagrams, we determined that introduction of Ub to UBQLN2 condensates incurs a significant inclusion energetic penalty. This penalty antagonizes the ability of polyUb hubs to scaffold multiple UBQLN2 molecules and cooperatively amplify phase separation. Importantly, the extent to which polyubiquitin hubs can promote UBQLN2 phase separation are encoded in the spacings between Ub units as found for naturally-occurring chains of different linkages and designed chains of different architectures, thus illustrating how the ubiquitin code regulates functionality via the emergent properties of the condensate. We expect our findings to extend to other condensates necessitating the consideration of ligand properties, including concentration, valency, affinity, and spacing between binding sites in studies and designs of condensates.

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.

Reduction of oligomer size modulates the competition between cluster formation and phase separation of the tumor suppressor SPOP

Phase separation is a ubiquitous process that compartmentalizes many cellular pathways. Given that the same interactions that drive phase separation mediate the formation of complexes below the saturation concentration, the contribution of condensates vs complexes to function is not always clear. Here, we characterized several new cancer-associated mutations of the tumor suppressor Speckle-type POZ protein (SPOP), a substrate recognition subunit of the Cullin3-RING ubiquitin ligase (CRL3), which pointed to a strategy for generating separation-of-function mutations. SPOP self-associates into linear oligomers and interacts with multivalent substrates, and this mediates the formation of condensates. These condensates bear the hallmarks of enzymatic ubiquitination activity. We characterized the effect of mutations in the dimerization domains of SPOP on its linear oligomerization, binding to the substrate DAXX, and phase separation with DAXX. We showed that the mutations reduce SPOP oligomerization and shift the size distribution of SPOP oligomers to smaller sizes. The mutations therefore reduce the binding affinity to DAXX, but enhance the poly-ubiquitination activity of SPOP towards DAXX. This unexpectedly enhanced activity may be explained by enhanced phase separation of DAXX with the SPOP mutants. Our results provide a comparative assessment of the functional role of clusters versus condensates and support a model in which phase separation is an important factor in SPOP function. Our findings also suggest that tuning of linear SPOP self-association could be used by the cell to modulate its activity, and provide insights into the mechanisms underlying hypermorphic SPOP mutations. The characteristics of these cancer-associated SPOP mutations suggest a route for designing separation-of-function mutations in other phase-separating systems.

Conformational and oligomeric states of SPOP from small-angle X-ray scattering and molecular dynamics simulations

Speckle-type POZ protein (SPOP) is a substrate adaptor in the ubiquitin proteasome system, and plays important roles in cell-cycle control, development, and cancer pathogenesis. SPOP forms linear higher-order oligomers following an isodesmic self-association model. Oligomerization is essential for SPOP's multivalent interactions with substrates, which facilitate phase separation and localization to biomolecular condensates. Structural characterization of SPOP in its oligomeric state and in solution is, however, challenging due to the inherent conformational and compositional heterogeneity of the oligomeric species. Here, we develop an approach to simultaneously and self-consistently characterize the conformational ensemble and the distribution of oligomeric states of SPOP by combining small-angle X-ray scattering (SAXS) and molecular dynamics (MD) simulations. We build initial conformational ensembles of SPOP oligomers using coarse-grained molecular dynamics simulations, and use a Bayesian/maximum entropy approach to refine the ensembles, along with the distribution of oligomeric states, against a concentration series of SAXS experiments. Our results suggest that SPOP oligomers behave as rigid, helical structures in solution, and that a flexible linker region allows SPOP's substrate-binding domains to extend away from the core of the oligomers. Additionally, our results are in good agreement with previous characterization of the isodesmic self-association of SPOP. In the future, the approach presented here can be extended to other systems to simultaneously characterize structural heterogeneity and self-assembly.

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.

Multiscale Modeling of Protein-RNA Condensation in and Out of Equilibrium

A vast number of intracellular membraneless bodies also known as biomolecular condensates form through a liquid-liquid phase separation (LLPS) of biomolecules. To date, phase separation has been identified as the main driving force for a membraneless organelles such as nucleoli, Cajal bodies, stress granules, and chromatin compartments. Recently, the protein-RNA condensation is receiving increased attention, because it is closely related to the biological function of cells such as transcription, translation, and RNA metabolism. Despite the multidisciplinary efforts put forth to study the biophysical properties of protein-RNA condensates, there are many fundamental unanswered questions regarding the mechanism of formation and regulation of protein-RNA condensates in eukaryotic cells. Major challenges in studying protein-RNA condensation stem from (i) the molecular heterogeneity and conformational flexibility of RNA and protein chains and (ii) the nonequilibrium nature of transcription and cellular environment. Computer simulations, bioinformatics, and mathematical models are uniquely positioned for shedding light on the microscopic nature of protein-RNA phase separation. To this end, there is an urgent need for innovative models with the right spatiotemporal resolution for confronting the experimental observables in a comprehensive and physics-based manner. In this chapter, we will summarize the currently emerging research efforts, which employ atomistic and coarse-grained molecular models and field theoretical models to understand equilibrium and nonequilibrium aspects of protein-RNA condensation.

Biological soft matter: intrinsically disordered proteins in liquid-liquid phase separation and biomolecular condensates

The facts that many proteins with crucial biological functions do not have unique structures and that many biological processes are compartmentalized into the liquid-like biomolecular condensates, which are formed via liquid-liquid phase separation (LLPS) and are not surrounded by the membrane, are revolutionizing the modern biology. These phenomena are interlinked, as the presence of intrinsic disorder represents an important requirement for a protein to undergo LLPS that drives biogenesis of numerous membrane-less organelles (MLOs). Therefore, one can consider these phenomena as crucial constituents of a new IDP-LLPS-MLO field. Furthermore, intrinsically disordered proteins (IDPs), LLPS, and MLOs represent a clear link between molecular and cellular biology and soft matter and condensed soft matter physics. Both IDP and LLPS/MLO fields are undergoing explosive development and generate the ever-increasing mountain of crucial data. These new data provide answers to so many long-standing questions that it is difficult to imagine that in the very recent past, protein scientists and cellular biologists operated without taking these revolutionary concepts into account. The goal of this essay is not to deliver a comprehensive review of the IDP-LLPS-MLO field but to provide a brief and rather subjective outline of some of the recent developments in these exciting fields.

Uncovering the molecular mechanism for dual effect of ATP on phase separation in FUS solution

Recent studies reported that adenosine triphosphate (ATP) could inhibit and enhance the phase separation in prion-like proteins. The molecular mechanism underlying such a puzzling phenomenon remains elusive. Here, taking the fused in sarcoma (FUS) solution as an example, we comprehensively reveal the underlying mechanism by which ATP regulates phase separation by combining the semiempirical quantum mechanical method, mean-field theory, and molecular simulation. At the microscopic level, ATP acts as a bivalent or trivalent binder; at the macroscopic level, the reentrant phase separation occurs in dilute FUS solutions, resulting from the ATP concentration-dependent binding ability under different conditions. The ATP concentration for dissolving the protein condensates is about 10 mM, agreeing with experimental results. Furthermore, from a dynamic point of view, the effect of ATP on phase separation is also nonmonotonic. This work provides a clear physical description of the microscopic interaction and macroscopic phase diagram of the ATP-modulated phase separation.

Model biomolecular condensates have heterogeneous structure quantitatively dependent on the interaction profile of their constituent macromolecules

Biomolecular condensates play numerous roles in cells by selectively concentrating client proteins while excluding others. These functions are likely to be sensitive to the spatial organization of the scaffold proteins forming the condensate. We use coarse-grained molecular simulations to show that model intrinsically-disordered proteins phase separate into a heterogeneous, structured fluid characterized by a well-defined length scale. The proteins are modelled as semi-flexible polymers with punctate, multifunctional binding sites in good solvent conditions. Their dense phase is highly solvated with a spatial structure that is more sensitive to the separation of the binding sites than their affinity. We introduce graph theoretic measures to quantify their heterogeneity, and find that it increases with increasing binding site number, and exhibits multi-timescale dynamics. The model proteins also swell on passing from the dilute solution to the dense phase. The simulations predict that the structure of the dense phase is modulated by the location and affinity of binding sites distant from the termini of the proteins, while sites near the termini more strongly affect its phase behaviour. The relations uncovered between the arrangement of weak interaction sites on disordered proteins and the material properties of their dense phase can be experimentally tested to give insight into the biophysical properties, pathological effects, and rational design of biomolecular condensates.

A conceptual framework for understanding phase separation and addressing open questions and challenges

Macromolecular phase separation is being recognized for its potential importance and relevance as a driver of spatial organization within cells. Here, we describe a framework based on synergies between networking (percolation or gelation) and density (phase separation) transitions. Accordingly, the phase transitions in question are referred to as phase separation coupled to percolation (PSCP). The condensates that result from PSCP are viscoelastic network fluids. Such systems have sequence-, composition-, and topology-specific internal network structures that give rise to time-dependent interplays between viscous and elastic properties. Unlike pure phase separation, the process of PSCP gives rise to sequence-, chemistry-, and structure-specific distributions of clusters that can form at concentrations that lie well below the threshold concentration for phase separation. PSCP, influenced by specific versus solubility-determining interactions, also provides a bridge between different observations and helps answer questions and address challenges that have arisen regarding the role of macromolecular phase separation in biology.

Rules of Physical Mathematics Govern Intrinsically Disordered Proteins

In stark contrast to foldable proteins with a unique folded state, intrinsically disordered proteins and regions (IDPs) persist in perpetually disordered ensembles. Yet an IDP ensemble has conformational features-even when averaged-that are specific to its sequence. In fact, subtle changes in an IDP sequence can modulate its conformational features and its function. Recent advances in theoretical physics reveal a set of elegant mathematical expressions that describe the intricate relationships among IDP sequences, their ensemble conformations, and the regulation of their biological functions. These equations also describe the molecular properties of IDP sequences that predict similarities and dissimilarities in their functions and facilitate classification of sequences by function, an unmet challenge to traditional bioinformatics. These physical sequence-patterning metrics offer a promising new avenue for advancing synthetic biology at a time when multiple novel functional modes mediated by IDPs are emerging.

Assembly of model postsynaptic densities involves interactions auxiliary to stoichiometric binding

The assembly of functional biomolecular condensates often involves liquid-liquid phase separation (LLPS) of proteins with multiple modular domains, which can be folded or conformationally disordered to various degrees. To understand the LLPS-driving domain-domain interactions, a fundamental question is how readily the interactions in the condensed phase can be inferred from interdomain interactions in dilute solutions. In particular, are the interactions leading to LLPS exclusively those underlying the formation of discrete interdomain complexes in homogeneous solutions? We address this question by developing a mean-field LLPS theory of two stoichiometrically constrained solute species. The theory is applied to the neuronal proteins SynGAP and PSD-95, whose complex coacervate serves as a rudimentary model for neuronal postsynaptic densities (PSDs). The predicted phase behaviors are compared with experiments. Previously, a three SynGAP/two PSD-95 ratio was determined for SynGAP/PSD-95 complexes in dilute solutions. However, when this 3:2 stoichiometry is uniformly imposed in our theory encompassing both dilute and condensed phases, the tie-line pattern of the predicted SynGAP/PSD-95 phase diagram differs drastically from that obtained experimentally. In contrast, theories embodying alternate scenarios postulating auxiliary SynGAP-PSD-95 as well as SynGAP-SynGAP and PSD-95-PSD-95 interactions, in addition to those responsible for stoichiometric SynGAP/PSD-95 complexes, produce tie-line patterns consistent with experiment. Hence, our combined theoretical-experimental analysis indicates that weaker interactions or higher-order complexes beyond the 3:2 stoichiometry, but not yet documented, are involved in the formation of SynGAP/PSD-95 condensates, imploring future efforts to ascertain the nature of these auxiliary interactions in PSD-like LLPS and underscoring a likely general synergy between stoichiometric, structurally specific binding and stochastic, multivalent "fuzzy" interactions in the assembly of functional biomolecular condensates.

Motif-pattern dependence of biomolecular phase separation driven by specific interactions

Eukaryotic cells partition a wide variety of important materials and processes into biomolecular condensates—phase-separated droplets that lack a membrane. In addition to nonspecific electrostatic or hydrophobic interactions, phase separation also depends on specific binding motifs that link together constituent molecules. Nevertheless, few rules have been established for how these ubiquitous specific, saturating, motif-motif interactions drive phase separation. By integrating Monte Carlo simulations of lattice-polymers with mean-field theory, we show that the sequence of heterotypic binding motifs strongly affects a polymer’s ability to phase separate, influencing both phase boundaries and condensate properties (e.g. viscosity and polymer diffusion). We find that sequences with large blocks of single motifs typically form more inter-polymer bonds, which promotes phase separation. Notably, the sequence of binding motifs influences phase separation primarily by determining the conformational entropy of self-bonding by single polymers. This contrasts with systems where the molecular architecture primarily affects the energy of the dense phase, providing a new entropy-based mechanism for the biological control of phase separation.

Cells need to concentrate biomolecules in the right place at the right time in order to function. Many important intracellular compartments are liquid droplets formed by phase separation, the same process that separates oil from vinegar. The properties of such “biomolecular condensates” depend on the component molecules, such as proteins and RNAs. These molecules are polymers made of many interacting monomers, often organized into “motifs,” and the sequence of motifs shapes the properties of the condensates. Recent work has revealed important principles governing phase separation when the motifs are charged and interact across long distances, but many phase-separating molecules form specific interactions that are short-range and one-to-one. How does the sequence of specifically-interacting motifs affect phase separation? Using a combination of simulations and theoretical calculations, we show that the sequence has profound effects on both the formation and properties of condensates. Sequences with large blocks of identical motifs are better at phase separating but more viscous and solid-like. Importantly, we find that sequence controls phase separation via the proclivity to form self-bonds instead of forming bonds with other polymers. Thus the sequence of specifically-interacting motifs provides a control point for the formation and properties of phase-separated intracellular compartments.

Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Intrinsic disorder (ID) constitutes a new dimension to the protein structure-function relationship. The ability to undergo conformational changes upon binding is a key property of intrinsically disordered proteins and remains challenging to study using conventional methods. A 1994 paper by R. S. Spolar and M. T. Record presented a thermodynamic approach for estimating changes in conformational entropy based on heat capacity changes, allowing quantification of residues folding upon binding. Here, we adapt the method for studies of intrinsically disordered proteins. We integrate additional data to provide a broader experimental foundation for the underlying relations and, based on >500 protein-protein complexes involving disordered proteins, reassess a key relation between polar and nonpolar surface area changes, previously determined using globular protein folding. We demonstrate the improved suitability of the adapted method to studies of the folded αα-hub domain RST from radical-induced cell death 1, whose interactome is characterized by ID. From extensive thermodynamic data, quantifying the conformational entropy changes upon binding, and comparison to the NMR structure, the adapted method improves accuracy for ID-based studies. Furthermore, we apply the method, in conjunction with NMR, to reveal hitherto undetected effects of interaction-motif context. Thus, inclusion of the disordered context of the DREB2A RST-binding motif induces structuring of the binding motif, resulting in major enthalpy-entropy compensation in the interaction interface. This study, also evaluating additional interactions, demonstrates the strength of the ID-adapted Spolar-Record thermodynamic approach for dissection of structural features of ID-based interactions, easily overlooked in traditional studies, and for translation of these into mechanistic knowledge.

Nuclear Protein Condensates and Their Properties in Regulation of Gene Expression

Our understanding of the spatiotemporal regulation of eukaryotic gene expression has recently been greatly stimulated by the findings that many of the regulators of chromatin, transcription, and RNA processing form biomolecular condensates often assembled through liquid-liquid phase separation. Increasing number of reports suggest that these condensates functionally regulate gene expression, largely by concentrating the relevant biomolecules in the liquid-like micro-compartments. However, it remains poorly understood how the physicochemical properties, especially the material properties, of the condensates regulate gene expression activity. In this review, we discuss current data on various nuclear condensates and their biophysical properties with the underlying molecular interactions, and how they may functionally impact gene expression at the level of chromatin organization and activities, transcription, and RNA processing.

How Hierarchical Interactions Make Membraneless Organelles Tick Like Clockwork

Biomolecular condensates appear throughout the cell, serving many different biochemical functions. We argue that condensate functionality is optimized when the interactions driving condensation vary widely in affinity. Strong interactions provide structural specificity needed to encode functional properties but carry the risk of kinetic arrest, while weak interactions allow the system to remain dynamic but do not restrict the conformational ensemble enough to sustain specific functional features. To support our opinion, we describe illustrative examples of the interplay of strong and weak interactions that are found in the nucleolus, SPOP/DAXX condensates, polySUMO/polySIM condensates, chromatin, and stress granules. The common feature of these systems is a hierarchical assembly motif in which weak, transient interactions condense structurally defined functional units.

FUS oncofusion protein condensates recruit mSWI/SNF chromatin remodeler via heterotypic interactions between prion-like domains

Fusion transcription factors generated by genomic translocations are common drivers of several types of cancers including sarcomas and leukemias. Oncofusions of the FET (FUS, EWSR1, and TAF15) family proteins result from the fusion of the prion-like domain (PLD) of FET proteins to the DNA-binding domain (DBD) of certain transcription regulators and are implicated in aberrant transcriptional programs through interactions with chromatin remodelers. Here, we show that FUS-DDIT3, a FET oncofusion protein, undergoes PLD-mediated phase separation into liquid-like condensates. Nuclear FUS-DDIT3 condensates can recruit essential components of the global transcriptional machinery such as the chromatin remodeler SWI/SNF. The recruitment of mammalian SWI/SNF (mSWI/SNF) is driven by heterotypic PLD-PLD interactions between FUS-DDIT3 and core subunits of SWI/SNF, such as the catalytic component BRG1. Further experiments with single-molecule correlative force-fluorescence microscopy support a model wherein the fusion protein forms condensates on DNA surface and enrich BRG1 to activate transcription by ectopic chromatin remodeling. Similar PLD-driven co-condensation of mSWI/SNF with transcription factors can be employed by other oncogenic fusion proteins with a generic PLD-DBD domain architecture for global transcriptional reprogramming.

Intrinsically disordered substrates dictate SPOP subnuclear localization and ubiquitination activity

Speckle-type POZ protein (SPOP) is a ubiquitin ligase adaptor that binds substrate proteins and facilitates their proteasomal degradation. Most SPOP substrates present multiple SPOP-binding (SB) motifs and undergo liquid-liquid phase separation with SPOP. Pancreatic and duodenal homeobox 1 (Pdx1), an insulin transcription factor, is downregulated by interaction with SPOP. Unlike other substrates, only one SB motif has previously been reported within the Pdx1 C-terminal intrinsically disordered region (Pdx1-C). Given this difference, we aimed to determine the specific mode of interaction of Pdx1 with SPOP and how it is similar or different to that of other SPOP substrates. Here, we identify a second SB motif in Pdx1-C, but still find that the resulting moderate valency is insufficient to support phase separation with SPOP in cells. Although Pdx1 does not phase separate with SPOP, Pdx1 and SPOP interaction prompts SPOP relocalization from nuclear speckles to the diffuse nucleoplasm. Accordingly, we find that SPOP-mediated ubiquitination activity of Pdx1 occurs in the nucleoplasm and that highly efficient Pdx1 turnover requires both SB motifs. Our results suggest that the subnuclear localization of SPOP-substrate interactions and substrate ubiquitination may be directed by the properties of the substrate itself.

Thermodynamics and kinetics of phase separation of protein-RNA mixtures by a minimal model

Intracellular liquid-liquid phase separation enables the formation of biomolecular condensates, such as ribonucleoprotein granules, which play a crucial role in the spatiotemporal organization of biomolecules (e.g., proteins and RNAs). Here, we introduce a patchy-particle polymer model to investigate liquid-liquid phase separation of protein-RNA mixtures. We demonstrate that at low to moderate concentrations, RNA enhances the stability of RNA-binding protein condensates because it increases the molecular connectivity of the condensed-liquid phase. Importantly, we find that RNA can also accelerate the nucleation stage of phase separation. Additionally, we assess how the capacity of RNA to increase the stability of condensates is modulated by the relative protein-protein/protein-RNA binding strengths. We find that phase separation and multiphase organization of multicomponent condensates is favored when the RNA binds with higher affinity to the lower-valency proteins in the mixture than to the cognate higher-valency proteins. Collectively, our results shed light on the roles of RNA in ribonucleoprotein granule formation and the internal structuring of stress granules.

Viscoelasticity of biomolecular condensates conforms to the Jeffreys model

Biomolecular condensates, largely by virtue of their material properties, are revolutionizing biology, and yet, the physical understanding of these properties is lagging. Here, I show that the viscoelasticity of condensates can be captured by a simple model, comprising a component where shear relaxation is an exponential function (with time constant τ1) and a component with nearly instantaneous shear relaxation (time constant τ0 → 0). Modulation of intermolecular interactions, e.g., by adding salt, can disparately affect the two components such that the τ1 component may dominate at low salt, whereas the τ0 component may dominate at high salt. Condensates have a tendency to fuse, with the dynamics accelerated by interfacial tension and impeded by viscosity. For fast-fusion condensates, shear relaxation on the τ1 timescale may become rate-limiting such that the fusion speed is no longer in direction proportion to the interfacial tension. These insights help narrow the gap in understanding between the biology and physics of biomolecular condensates.

Structure-Function Properties in Disordered Condensates

Biomolecular condensates appear throughout the cell serving a wide variety of functions. Many condensates appear to form by the assembly of multivalent molecules, which produce phase-separated networks with liquidlike properties. These networks then recruit client molecules, with the total composition providing functionality. Here we use a model system of poly-SUMO and poly-SIM proteins to understand client-network interactions and find that the structure of the network plays a strong role in defining client recruitment and thus functionality. The basic unit of assembly in this system is a zipperlike filament composed of alternating poly-SUMO and poly-SIM molecules. These filaments have defects of unsatisfied bonds that allow for both the formation of a 3D network and the recruitment of clients. The filamentous structure constrains the scaffold stoichiometries and the distribution of client recruitment sites that the network can accommodate. This results in a nonmonotonic client binding response that can be tuned independently by the client valence and binding energy. These results show how the interactions within liquid states can be disordered yet still contain structural features that provide functionality to the condensate.

Comparative roles of charge, π, and hydrophobic interactions in sequence-dependent phase separation of intrinsically disordered proteins

Endeavoring toward a transferable, predictive coarse-grained explicit-chain model for biomolecular condensates underlain by liquid-liquid phase separation (LLPS) of proteins, we conducted multiple-chain simulations of the N-terminal intrinsically disordered region (IDR) of DEAD-box helicase Ddx4, as a test case, to assess roles of electrostatic, hydrophobic, cation-π, and aromatic interactions in amino acid sequence-dependent LLPS. We evaluated three different residue-residue interaction schemes with a shared electrostatic potential. Neither a common hydrophobicity scheme nor one augmented with arginine/lysine-aromatic cation-π interactions consistently accounted for available experimental LLPS data on the wild-type, a charge-scrambled, a phenylalanine-to-alanine (FtoA), and an arginine-to-lysine (RtoK) mutant of Ddx4 IDR. In contrast, interactions based on contact statistics among folded globular protein structures reproduce the overall experimental trend, including that the RtoK mutant has a much diminished LLPS propensity. Consistency between simulation and experiment was also found for RtoK mutants of P-granule protein LAF-1, underscoring that, to a degree, important LLPS-driving π-related interactions are embodied in classical statistical potentials. Further elucidation is necessary, however, especially of phenylalanine's role in condensate assembly because experiments on FtoA and tyrosine-to-phenylalanine mutants suggest that LLPS-driving phenylalanine interactions are significantly weaker than posited by common statistical potentials. Protein-protein electrostatic interactions are modulated by relative permittivity, which in general depends on aqueous protein concentration. Analytical theory suggests that this dependence entails enhanced interprotein interactions in the condensed phase but more favorable protein-solvent interactions in the dilute phase. The opposing trends lead to only a modest overall impact on LLPS.

Tug of War between Condensate Phases in a Minimal Macromolecular System

Membraneless organelles, comprising dozens to hundreds of macromolecular components, form heterogeneous phases in space and evolve over time in material properties. Here, using four macromolecules, we demonstrate a range of phase behaviors associated with membraneless organelles and uncover the underlying physicochemical rules. The macromolecules are SH3₅ (S) and PRM₅ (P), two pentameric, oppositely charged protein constructs; heparin (H), an anionic polymer; and lysozyme (L), a cationic single-domain protein. The S:P, S:L, and P:H binaries form droplets, but the H:L binary forms network-like precipitates, therefore setting up a tug of war between different condensate phases within the S:P:H:L quaternary. The H:L exception can partly be attributed to the compactness of L, as supported by ThT binding data. Increasing amounts of P alone or both S and P, but not S alone, can dissolve H:L precipitates into droplets. These differential effects can be explained by the order of the strengths of pairwise attraction: H:L > P:H > S:P > S:L, deduced from the shapes of ternary phase boundaries. When S and P are at subdissolution concentrations, S:P:H:L precipitates change over time to become droplet-like in appearance, although not completely fluidic according to fluorescence recovery after photobleaching. In fact, confocal microscopy reveals separated S:L-rich and P:H-rich foci inside the droplet-like condensates. Therefore, complex phase behaviors of membraneless organelles, including rescue of aberrant phase transitions, demixing of condensates, and time evolution of material properties, can all be reconstituted and understood via a minimal macromolecular system.

Molecular structure in biomolecular condensates

Evidence accumulated over the past decade provides support for liquid-liquid phase separation as the mechanism underlying the formation of biomolecular condensates, which include not only 'membraneless' organelles such as nucleoli and RNA granules, but additional assemblies involved in transcription, translation and signaling. Understanding the molecular mechanisms of condensate function requires knowledge of the structures of their constituents. Current knowledge suggests that structures formed via multivalent domain-motif interactions remain largely unchanged within condensates. Two different viewpoints exist regarding structures of disordered low-complexity domains within condensates; one argues that low-complexity domains remain largely disordered in condensates and their multivalency is encoded in short motifs called 'stickers', while the other argues that the sequences form cross-β structures resembling amyloid fibrils. We review these viewpoints and highlight outstanding questions that will inform structure-function relationships for biomolecular condensates.