A Lattice Model of Charge-Pattern-Dependent Polyampholyte Phase Separation

How do intrinsically disordered protein regions encode a driving force for liquid-liquid phase separation?

Liquid-liquid phase separation is the mechanism underlying the formation of biomolecular condensates. Disordered protein regions often drive phase separation, but the molecular interactions mediating this phenomenon are not well understood, sometimes leading to the conflation that all disordered protein regions drive phase separation. Given the critical role of phase separation in many cellular processes, and that dysfunction of phase separation can lead to debilitating diseases, it is important that we understand the interactions and sequence properties underlying phase behavior. A conceptual framework that divides IDRs into interacting and solvating regions has proven particularly useful, and analytical instantiations and coarse-grained models can test our understanding of the driving forces against experimental phase behavior. Validated simulation paradigms enable the exploration of sequence space to help our understanding of how disordered protein regions can encode phase behavior, which IDRs may mediate phase separation in cells, and which IDRs are highly soluble.

Analytical Theory for Sequence-Specific Binary Fuzzy Complexes of Charged Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are important for biological functions. In contrast to folded proteins, molecular recognition among certain IDPs is "fuzzy" in that their binding and/or phase separation are stochastically governed by the interacting IDPs' amino acid sequences, while their assembled conformations remain largely disordered. To help elucidate a basic aspect of this fascinating yet poorly understood phenomenon, the binding of a homo or heterodimeric pair of polyampholytic IDPs is modeled statistical mechanically using cluster expansion. We find that the binding affinities of binary fuzzy complexes in the model correlate strongly with a newly derived simple "joint sequence charge decoration" parameter readily calculable from the pair of IDPs' sequence charge patterns. Predictions by our analytical theory are in essential agreement with coarse-grained explicit-chain simulations. This computationally efficient theoretical framework is expected to be broadly applicable to rationalizing and predicting sequence-specific IDP-IDP polyelectrostatic interactions.

Single-Protein Collapse Determines Phase Equilibria of a Biological Condensate

Recent advances in microscopy of living cells have established membraneless organelles as critical elements of diverse biological processes. The body of experimental work suggests that formation of such organelles is driven by liquid-liquid phase separation, a physical process that has been studied extensively for both simple liquids and mixtures of polymers. Here, we combine molecular dynamics simulations with polymer theory to show that the thermodynamic behavior of one particular biomolecular condensate-fused in sarcoma (FUS)-can be quantitatively accounted for at the level of the chain collapse theory. First, we show that a particle-based molecular dynamics model can reproduce known phase separation properties of a FUS condensate, including its critical concentration and susceptibility to mutations. Next, we obtain a polymer physics representation of a FUS condensate by examining the behavior of a single FUS protein as a function of temperature. We use the chain collapse theory to determine the thermodynamic properties of the condensate and to characterize changes in the single-chain conformation at the onset of phase separation. Altogether, our findings suggest that the phase behavior of FUS condensates can be explained by the properties of individual FUS proteins and that the change in the FUS conformation is the main force driving for the phase separation.

A unified analytical theory of heteropolymers for sequence-specific phase behaviors of polyelectrolytes and polyampholytes

The physical chemistry of liquid-liquid phase separation (LLPS) of polymer solutions bears directly on the assembly of biologically functional dropletlike bodies from proteins and nucleic acids. These biomolecular condensates include certain extracellular materials and intracellular compartments that are characterized as "membraneless organelles." Analytical theories are a valuable, computationally efficient tool for addressing general principles. LLPS of neutral homopolymers is quite well described by theory, but it has been a challenge to develop general theories for the LLPS of heteropolymers involving charge-charge interactions. Here, we present a theory that combines a random-phase-approximation treatment of polymer density fluctuations and an account of intrachain conformational heterogeneity based on renormalized Kuhn lengths to provide predictions of LLPS properties as a function of pH, salt, and charge patterning along the chain sequence. Advancing beyond more limited analytical approaches, our LLPS theory is applicable to a wide variety of charged sequences ranging from highly charged polyelectrolytes to neutral or nearly neutral polyampholytes. This theory should be useful in high-throughput screening of protein and other sequences for their LLPS propensities and can serve as a basis for more comprehensive theories that incorporate nonelectrostatic interactions. Experimental ramifications of our theory are discussed.

Phase behavior of blocky charge lattice polymers: Crystals, liquids, sheets, filaments, and clusters

Motivated by the idea that intrinsically disordered proteins (IDPs) condense into liquidlike droplets within cells, we carry out Monte Carlo simulations of a polymer lattice model to study the relationship between charge patterning and phase separation. Polymer chains containing neutral, positively charged, and negatively charged monomers are placed on a cubic lattice. Only nearest-neighbor interactions between charges are considered. We determine the phase diagram for a systematically varied set of sequences. We observe homogeneous fluids, liquid condensation, cluster phases, filaments, and crystal states. Of the six sequences we study, three form crystals at low temperatures. The other three sequences, which have lower charge densities, instead collapse into gel-like networks or unconnected finite clusters. Longer neutral patches along the sequence sterically limit the size and shape of low-energy structures, which is analogous to the effect of charge or limited valence in attractive colloids. Only one sequence clearly exhibits liquid behavior; this sequence has a reduced tendency to individually fold and crystallize compared to others of similar charge density and draws parallels to real IDP behavior.

Binding without folding - the biomolecular function of disordered polyelectrolyte complexes

Recent evidence shows that oppositely charged intrinsically disordered proteins (IDPs) can form high-affinity complexes that involve neither the formation of secondary or tertiary structure nor site-specific interactions between individual residues. Similar electrostatically dominated interactions have also been identified for positively charged IDPs binding to nucleic acids. These highly disordered polyelectrolyte complexes constitute an extreme case within the spectrum of biomolecular interactions involving disorder. Such interactions are likely to be widespread, since sequence analysis predicts proteins with highly charged disordered regions to be surprisingly numerous. Here, we summarize the insights that have emerged from the highly disordered polyelectrolyte complexes identified so far, and we highlight recent developments and future challenges in (i) establishing models for the underlying highly dynamic structural ensembles, (ii) understanding the novel binding mechanisms associated with them, and (iii) identifying the functional consequences.

Low-complexity domain of U1-70K modulates phase separation and aggregation through distinctive basic-acidic motifs

Liquid-liquid phase separation (LLPS) facilitates the formation of functional membraneless organelles and recent reports have linked this phenomenon to protein aggregation in neurodegenerative diseases. Understanding the mechanism of LLPS and its regulation thus promises to shed light on the pathogenesis of these conditions. The RNA-binding protein U1-70K, which aggregates in brains of Alzheimer's disease patients, is considered a potential target for Alzheimer's therapy. Here, we report that two fragments in the low-complexity (LC) domain of U1-70K can undergo LLPS. We have demonstrated that the repetitive basic-acidic motifs in these fragments induce nucleotide-independent phase separation and initiate aggregation in vitro. We also have confirmed that LLPS and aggregation occur in vivo and that the content of ampholytic motifs in a protein domain determines the transition between droplets and aggregation, providing insights into the mechanism underlying the formation of diverse assembly states.

LASSI: A lattice model for simulating phase transitions of multivalent proteins

Many biomolecular condensates form via spontaneous phase transitions that are driven by multivalent proteins. These molecules are biological instantiations of associative polymers that conform to a so-called stickers-and-spacers architecture. The stickers are protein-protein or protein-RNA interaction motifs and / or domains that can form reversible, non-covalent crosslinks with one another. Spacers are interspersed between stickers and their preferential interactions with solvent molecules determine the cooperativity of phase transitions. Here, we report the development of an open source computational engine known as LASSI (LAttice simulation engine for Sticker and Spacer Interactions) that enables the calculation of full phase diagrams for multicomponent systems comprising of coarse-grained representations of multivalent proteins. LASSI is designed to enable computationally efficient phenomenological modeling of spontaneous phase transitions of multicomponent mixtures comprising of multivalent proteins and RNA molecules. We demonstrate the application of LASSI using simulations of linear and branched multivalent proteins. We show that dense phases are best described as droplet-spanning networks that are characterized by reversible physical crosslinks among multivalent proteins. We connect recent observations regarding correlations between apparent stoichiometry and dwell times of condensates to being proxies for the internal structural organization, specifically the convolution of internal density and extent of networking, within condensates. Finally, we demonstrate that the concept of saturation concentration thresholds does not apply to multicomponent systems where obligate heterotypic interactions drive phase transitions. This emerges from the ellipsoidal structures of phase diagrams for multicomponent systems and it has direct implications for the regulation of biomolecular condensates in vivo.

Temperature, Hydrostatic Pressure, and Osmolyte Effects on Liquid-Liquid Phase Separation in Protein Condensates: Physical Chemistry and Biological Implications

Liquid-liquid phase separation (LLPS) of proteins and other biomolecules play a critical role in the organization of extracellular materials and membrane-less compartmentalization of intra-organismal spaces through the formation of condensates. Structural properties of such mesoscopic droplet-like states were studied by spectroscopy, microscopy, and other biophysical techniques. The temperature dependence of biomolecular LLPS has been studied extensively, indicating that phase-separated condensed states of proteins can be stabilized or destabilized by increasing temperature. In contrast, the physical and biological significance of hydrostatic pressure on LLPS is less appreciated. Summarized here are recent investigations of protein LLPS under pressures up to the kbar-regime. Strikingly, for the cases studied thus far, LLPSs of both globular proteins and intrinsically disordered proteins/regions are typically more sensitive to pressure than the folding of proteins, suggesting that organisms inhabiting the deep sea and sub-seafloor sediments, under pressures up to 1 kbar and beyond, have to mitigate this pressure-sensitivity to avoid unwanted destabilization of their functional biomolecular condensates. Interestingly, we found that trimethylamine-N-oxide (TMAO), an osmolyte upregulated in deep-sea fish, can significantly stabilize protein droplets under pressure, pointing to another adaptive advantage for increased TMAO concentrations in deep-sea organisms besides the osmolyte's stabilizing effect against protein unfolding. As life on Earth might have originated in the deep sea, pressure-dependent LLPS is pertinent to questions regarding prebiotic proto-cells. Herein, we offer a conceptual framework for rationalizing the recent experimental findings and present an outline of the basic thermodynamics of temperature-, pressure-, and osmolyte-dependent LLPS as well as a molecular-level statistical mechanics picture in terms of solvent-mediated interactions and void volumes.

Small ion effects on self-coacervation phenomena in block polyampholytes

Self-coacervation is a phenomenon in which a solution of polyampholytes spontaneously phase separates into a dense liquid coacervate phase, rich in the polyampholyte, coexisting with a dilute supernatant phase. Such coacervation results in the formation of membraneless organelles in vivo and has further been applied industrially as synthetic encapsulants and coatings. It has been suggested that coacervation is primarily driven by the entropy gain from releasing counter-ions upon complexation. Using fully fluctuating field-theoretic simulations employing complex Langevin sampling and complementary molecular dynamics simulations, we have determined that the small ions contribute only weakly to the self-coacervation behavior of charge-symmetric block polyampholytes in solution. Salt partitioning between the supernatant and coacervate is also found to be negligible in the weak-binding regime at low electrostatic strengths. Asymmetries in charge distribution along the polyampholytes can cause net-charges that lead to "tadpole" configurations in dilute solution and the suppression of phase separation at low salt content. The field and particle-based simulation results are compared with analytical predictions from the random phase approximation (RPA) and postulated scaling relationships. The qualitative trends are mostly captured by the RPA, but the approximation fails at low concentration.

Hydrophobic confinement modulates thermal stability and assists knotting in the folding of tangled proteins

There is growing support for the idea that the in vivo folding process of knotted proteins is assisted by chaperonins, but the mechanism of chaperonin assisted folding remains elusive. Here, we conduct extensive Monte Carlo simulations of lattice and off-lattice models to explore the effects of confinement and hydrophobic intermolecular interactions with the chaperonin cage in the folding and knotting processes. We find that moderate to high protein-cavity interactions (which are likely to be established in the beginning of the chaperonin working cycle) cause an energetic destabilization of the protein that overcomes the entropic stabilization driven by excluded volume, and leads to a decrease of the melting temperature relative to bulk conditions. Moreover, mild-to-moderate hydrophobic interactions with the cavity (which would be established later in the cycle) lead to a significant enhancement of knotting probability in relation to bulk conditions while simultaneously moderating the effect of steric confinement in the enhancement of thermal stability. Our results thus indicate that the chaperonin may be able to assist knotting without simultaneously thermally stabilizing potential misfolded states to a point that would hamper productive folding thus compromising its functional role.

Complete Phase Diagram for Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins

A number of intrinsically disordered proteins have been shown to self-assemble via liquid-liquid phase separation into protein-rich and dilute phases. The resulting coacervates can have important biological functions, and the ability to form these assemblies is dictated by the protein's primary amino acid sequence as well as by the solution conditions. We present a complete phase diagram for the simple coacervation of a polyampholyte intrinsically disordered protein using a field-theoretic simulation approach. We show that differences in the primary amino acid sequence and in the distribution of charged amino acids along the sequence lead to differences in the phase window for coacervation, with block-charged sequences having a larger coacervation window than sequences with a random patterning of charges. The model also captures how changing solution conditions modifies the phase diagram and can serve to guide experimental studies.

Narrow equilibrium window for complex coacervation of tau and RNA under cellular conditions

The mechanism that leads to liquid-liquid phase separation (LLPS) of the tau protein, whose pathological aggregation is implicated in neurodegenerative disorders, is not well understood. Establishing a phase diagram that delineates the boundaries of phase co-existence is key to understanding whether LLPS is an equilibrium or intermediate state. We demonstrate that tau and RNA reversibly form complex coacervates. While the equilibrium phase diagram can be fit to an analytical theory, a more advanced model is investigated through field theoretic simulations (FTS) that provided direct insight into the thermodynamic driving forces of tau LLPS. Together, experiment and simulation reveal that tau-RNA LLPS is stable within a narrow equilibrium window near physiological conditions over experimentally tunable parameters including temperature, salt and tau concentrations, and is entropy-driven. Guided by our phase diagram, we show that tau can be driven toward LLPS under live cell coculturing conditions with rationally chosen experimental parameters.

Molecular design of self-coacervation phenomena in block polyampholytes

Coacervation is a common phenomenon in natural polymers and has been applied to synthetic materials systems for coatings, adhesives, and encapsulants. Single-component coacervates are formed when block polyampholytes exhibit self-coacervation, phase separating into a dense liquid coacervate phase rich in the polyampholyte coexisting with a dilute supernatant phase, a process implicated in the liquid-liquid phase separation of intrinsically disordered proteins. Using fully fluctuating field-theoretic simulations using complex Langevin sampling and complementary molecular-dynamics simulations, we develop molecular design principles to connect the sequenced charge pattern of a polyampholyte with its self-coacervation behavior in solution. In particular, the lengthscale of charged blocks and number of connections between oppositely charged blocks are shown to have a dramatic effect on the tendency to phase separate and on the accessible chain conformations. The field and particle-based simulation results are compared with analytical predictions from the random phase approximation (RPA) and postulated scaling relationships. The qualitative trends are mostly captured by the RPA, but the approximation fails catastrophically at low concentration.

Simulation methods for liquid-liquid phase separation of disordered proteins

Liquid-liquid phase separation of intrinsically disordered proteins (IDPs) and other biomolecules is a highly complex but robust process used by living systems. Drawing inspiration from biology, phase separating proteins have been successfully utilized for promising applications in fields of materials design and drug delivery. These protein-based materials are advantageous due to the ability to finely tune their stimulus-responsive phase behavior and material properties, and the ability to encode biologically active motifs directly into the sequence. The number of possible protein sequences is virtually endless, which makes sequence-based design a rather daunting task, but also attractive due to the amount of control coming from exploration of this variable space. The use of computational methods in this field of research have come to the aid in several aspects, including interpreting experimental results, identifying important structural features and molecular mechanisms capable of explaining the phase behavior, and ultimately providing predictive frameworks for rational design of protein sequences. Here we provide an overview of computational studies focused on phase separating biomolecules and the tools that are available to researchers interested in this topic.

Thermodynamically driven assemblies and liquid-liquid phase separations in biology

The sustenance of life depends on the high degree of organization that prevails through different levels of living organisms, from subcellular structures such as biomolecular complexes and organelles to tissues and organs. The physical origin of such organization is not fully understood, and even though it is clear that cells and organisms cannot maintain their integrity without consuming energy, there is growing evidence that individual assembly processes can be thermodynamically driven and occur spontaneously due to changes in thermodynamic variables such as intermolecular interactions and concentration. Understanding the phase separation in vivo requires a multidisciplinary approach, integrating the theory and physics of phase separation with experimental and computational techniques. This paper aims at providing a brief overview of the physics of phase separation and its biological implications, with a particular focus on the assembly of membraneless organelles. We discuss the underlying physical principles of phase separation from its thermodynamics to its kinetics. We also overview the wide range of methods utilized for experimental verification and characterization of phase separation of membraneless organelles, as well as the utility of molecular simulations rooted in thermodynamics and statistical physics in understanding the governing principles of thermodynamically driven biological self-assembly processes.

Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters

Biomolecular condensates undergirded by phase separations of proteins and nucleic acids serve crucial biological functions. To gain physical insights into their genetic basis, we study how liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) depends on their sequence charge patterns using a continuum Langevin chain model wherein each amino acid residue is represented by a single bead. Charge patterns are characterized by the "blockiness" measure κ and the "sequence charge decoration" (SCD) parameter. Consistent with random phase approximation (RPA) theory and lattice simulations, LLPS propensity as characterized by critical temperature Tcr* increases with increasingly negative SCD for a set of sequences showing a positive correlation between κ and -SCD. Relative to RPA, the simulated sequence-dependent variation in Tcr* is often-though not always-smaller, whereas the simulated critical volume fractions are higher. However, for a set of sequences exhibiting an anti-correlation between κ and -SCD, the simulated Tcr*'s are quite insensitive to either parameter. Additionally, we find that blocky sequences that allow for strong electrostatic repulsion can lead to coexistence curves with upward concavity as stipulated by RPA, but the LLPS propensity of a strictly alternating charge sequence was likely overestimated by RPA and lattice models because interchain stabilization of this sequence requires spatial alignments that are difficult to achieve in real space. These results help delineate the utility and limitations of the charge pattern parameters and of RPA, pointing to further efforts necessary for rationalizing the newly observed subtleties.

In silico Design of Laccase Thermostable Mutants From Lacc 6 of Pleurotus Ostreatus

Fungal laccase enzymes have a great biotechnological potential for bioremediation processes due to their ability to degrade compounds such as ρ-diphenol, aminophenols, polyphenols, polyamines, and aryldiamines. These enzymes have activity at different pH and temperature values, however, high temperatures can cause partial or total loss of enzymatic activity, so it is appropriate to do research to modify their secondary and/or tertiary structure to make them more resistant to extreme temperature conditions. In silico, a structure of the Lacc 6 enzyme of Pleurotus ostreatus was constructed using a laccase of Trametes versicolor as a template. From this structure, 16 mutants with possible resistance at high temperature due to ionic interactions, salt bridges and disulfide bonds were also obtained in silico. It was determined that 12 mutants called 4-DB, 3-DB, D233C-T310C, F468P, 3-SB, L132T, N79D, N372D, P203C, P203V, T147E, and W85F, presented the lowest thermodynamic energy. Based on the previous criterion and determining the least flexibility in the protein structures, three mutants (4-DB, 3-DB, and P203C) were selected, which may present high stability at high temperatures without affecting their active site. The obtained results allow the understanding of the molecular base that increase the structural stability of the enzyme Lacc 6 of Pleurotus ostreatus, achieving the in silico generation of mutants, which could have activity at high temperatures.

Conformational preferences and phase behavior of intrinsically disordered low complexity sequences: insights from multiscale simulations

While many proteins and protein regions utilize a complex repertoire of amino acids to achieve their biological function, a subset of protein sequences are enriched in a reduced set of amino acids. These so-called low complexity (LC) sequences, specifically intrinsically disordered variants of LC sequences, have been the focus of recent investigations owing to their roles in a range of biological functions, specifically phase separation. Computational studies of LC sequences have provided rich insights into their behavior both as individual proteins in dilute solutions and as the drivers and modulators of higher-order assemblies. Here, we review how simulations performed across distinct resolutions have provided different types of insights into the biological role of LC sequences.

Modulating charge patterning and ionic strength as a strategy to induce conformational changes in intrinsically disordered proteins

We present an analytical theory to describe conformational changes as a function of salt for polymers with a given sequence of charges. We apply this model to describe Intrinsically Disordered Proteins (IDPs) by explicitly accounting for charged residues and their exact placement in the primary sequence while approximating the effect of non-electrostatic interactions at a mean-field level by effective short-range (two body and three-body) interaction parameters. The effect of ions is introduced by treating electrostatic interactions within Debye-Huckle approximation. Using typical values of the short-range mean-field parameters derived from all-atom Monte Carlo simulations (at zero salt), we predict the conformational changes as a function of salt concentration. We notice that conformational transitions in response to changes in ionic strength strongly depend on sequence specific charge patterning. For example, globule to coil transition can be observed upon increasing salt concentration, in stark contrast to uniformly charged polyelectrolyte theories based on net charge only. In addition, it is possible to observe non-monotonic behavior with salt as well. Drastic differences in salt-induced conformational transitions is also evident between two doubly phosphorylated sequences-derived from the same wild type sequence-that only differ in the site of phosphorylation. Similar effects are also predicted between two sequences derived from the same parent sequence differing by a single site mutation where a negative charge is replaced by a positive charge. These effects are purely a result of charge decoration and can only be understood in terms of metrics based on specific placement of charges, and cannot be explained by models based on charge composition alone. Identifying sequences and hot spots within sequences-for post translational modification or charge mutation-using our high-throughput theory will yield fundamental insights into design and biological regulation mediated by phosphorylation and/or local changes in salt concentration.

Theories for Sequence-Dependent Phase Behaviors of Biomolecular Condensates

Liquid-liquid phase separation and related condensation processes of intrinsically disordered proteins (IDPs), proteins with intrinsically disordered regions, and nucleic acids underpin various condensed-liquid droplets or gel-like assemblies in the cellular environment. Collectively referred to as condensates, these bodies provide spatial/temporal compartmentalization, often serving as hubs for regulated biomolecular interactions. Examples include certain extracellular materials, transcription complexes, and membraneless organelles such as germ and stress granules and the nucleolus. They are critically important to cellular function; thus misregulation of their assembly is implicated in many diseases. Biomolecular condensates are complex entities. Our understanding of their inner workings is only in its infancy. Nonetheless, insights into basic biophysical principles of their assembly can be gained by applying analytical theories to elucidate how IDP phase behaviors are governed by the properties of the multivalent, solvent-mediated interactions entailed by the proteins' amino acid sequences. Here we briefly review the background of the pertinent polymer theories and outline the approximations that enable a tractable theoretical account of the dependence of IDP phase behaviors on the charge pattern of the IDP sequence. Of relevance to the homeostatic assembly of compositionally and functionally distinct condensates in the cellular context, theory indicates that the propensity for populations of different IDP sequences to mix or demix upon phase separation is affected by the similarity or dissimilarity of the sequence charge patterns. We also explore prospects of extending analytical theories to account for dynamic aspects of biomolecular condensates and to incorporate effects of cation-π, π-π, and temperature-dependent hydrophobic interactions on IDP phase properties.