Modeling phase transitions in mixtures of β-γ lens crystallins

Biophysical Principles Emerging from Experiments on Protein-Protein Association and Aggregation

Protein-protein association and aggregation are fundamental processes that play critical roles in various biological phenomena, from cellular signaling to disease progression. Understanding the underlying biophysical principles governing these processes is crucial for elucidating their mechanisms and developing strategies for therapeutic intervention. In this review, we provide an overview of recent experimental studies focused on protein-protein association and aggregation. We explore the key biophysical factors that influence these processes, including protein structure, conformational dynamics, and intermolecular interactions. We discuss the effects of environmental conditions such as temperature, pH and related buffer-specific effects, and ionic strength and related ion-specific effects on protein aggregation. The effects of polymer crowders and sugars are also addressed. We list the techniques used to study aggregation. We analyze emerging trends and challenges in the field, including the development of computational models and the integration of multidisciplinary approaches for a comprehensive understanding of protein-protein association and aggregation.

Antibody association in solution: cluster distributions and mechanisms

Understanding factors that affect the clustering and association of antibodies molecules in solution is critical to their development as therapeutics. For 19 different monoclonal antibody (mAb) solutions, we measured the viscosities, the second virial coefficients, the Kirkwood-Buff integrals, and the cluster distributions of the antibody molecules as functions of protein concentration. Solutions were modeled using the statistical-physics Wertheim liquid-solution theory, representing antibodies as Y-shaped molecular structures of seven beads each. We found that high-viscosity solutions result from more antibody molecules per cluster. Multi-body properties such as viscosity are well predicted experimentally by the 2-body Kirkwood-Buff quantity, G22, but not by the second virial coefficient, B22, and well-predicted theoretically from the Wertheim protein-protein sticking energy. Weakly interacting antibodies are rate-limited by nucleation; strongly interacting ones by propagation. This approach gives a way to relate micro to macro properties of solutions of associating proteins.

Protein Association in Solution: Statistical Mechanical Modeling

Protein molecules associate in solution, often in clusters beyond pairwise, leading to liquid phase separations and high viscosities. It is often impractical to study these multi-protein systems by atomistic computer simulations, particularly in multi-component solvents. Instead, their forces and states can be studied by liquid state statistical mechanics. However, past such approaches, such as the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, were limited to modeling proteins as spheres, and contained no microscopic structure-property relations. Recently, this limitation has been partly overcome by bringing the powerful Wertheim theory of associating molecules to bear on protein association equilibria. Here, we review these developments.

Engineering synthetic biomolecular condensates

The concept of phase-separation-mediated formation of biomolecular condensates provides a new framework to understand cellular organization and cooperativity-dependent cellular functions. With growing understanding of how biological systems drive phase separation and how cellular functions are encoded by biomolecular condensates, opportunities have emerged for cellular control through engineering of synthetic biomolecular condensates. In this Review, we discuss how to construct synthetic biomolecular condensates and how they can regulate cellular functions. We first describe the fundamental principles by which biomolecular components can drive phase separation. Next, we discuss the relationship between the properties of condensates and their cellular functions, which informs the design of components to create programmable synthetic condensates. Finally, we describe recent applications of synthetic biomolecular condensates for cellular control and discuss some of the design considerations and prospective applications.

Numerical Techniques for Applications of Analytical Theories to Sequence-Dependent Phase Separations of Intrinsically Disordered Proteins

Biomolecular condensates, physically underpinned to a significant extent by liquid-liquid phase separation (LLPS), are now widely recognized by numerous experimental studies to be of fundamental biological, biomedical, and biophysical importance. In the face of experimental discoveries, analytical formulations emerged as a powerful yet tractable tool in recent theoretical investigations of the role of LLPS in the assembly and dissociation of these condensates. The pertinent LLPS often involves, though not exclusively, intrinsically disordered proteins engaging in multivalent interactions that are governed by their amino acid sequences. For researchers interested in applying these theoretical methods, here we provide a practical guide to a set of computational techniques devised for extracting sequence-dependent LLPS properties from analytical formulations. The numerical procedures covered include those for the determination of spinodal and binodal phase boundaries from a general free energy function with examples based on the random phase approximation in polymer theory, construction of tie lines for multiple-component LLPS, and field-theoretic simulation of multiple-chain heteropolymeric systems using complex Langevin dynamics. Since a more accurate physical picture often requires comparing analytical theory against explicit-chain model predictions, a commonly utilized methodology for coarse-grained molecular dynamics simulations of sequence-specific LLPS is also briefly outlined.

Behaviour of the model antibody fluid constrained by rigid spherical obstacles: Effects of the obstacle-antibody attraction

This study investigates the behaviour of a fluid of monoclonal antibodies (mAbs) when trapped in a confinement represented by rigid spherical obstacles that attract antibodies. The antibody molecule is depicted as an assembly of seven hard spheres (7-bead model), organized to resemble a Y-shaped object. The model antibody has two Fab and one Fc domains located in the corners of letter Y. In this calculation, only the Fab-Fab and Fab-Fc attractive pairs of interactions are effective. The confinement is formed by the randomly distributed hard-spheres fixed in space. The spherical obstacles, besides the size exclusion, interact with beads of the antibody molecules via the Yukawa attractive potential. We applied the combination of the scaled particle theory, replica Ornstein-Zernike equations, Wertheim's thermodynamic perturbation approach and the Flory-Stockmayer theory to calculate: (i) the phase diagram of the liquid-liquid phase separation and the percolation threshold, (ii) the cluster size distributions, and (iii) the second virial coefficient of the protein fluid distributed among the obstacles. All these quantities were calculated as functions of the strength of the attraction between the monoclonal antibodies, and the monoclonal antibodies and obstacles. The conclusion is that while the hard-sphere obstacles decrease the critical density and the critical temperature of the mAbs fluid, the effect of the protein-obstacle attraction is more complex. Adding an attractive potential to the obstacle-mAbs interaction first increases the wideness of the T*-ρ envelope. However, with the further increase of the obstacle-mAbs attraction intensity, we observe reversal of the effect, the T*-ρ curves become narrower. At some point, depending on the obstacle-mAbs interaction, the situation is observed where two different temperatures have the same fluid density (re-entry point). In all the cases shown here the critical point decreases below the value for the neat fluid, but the behaviour with respect to an increase of the strength of the obstacle-mAbs attraction is not monotonic. Yet another interesting phenomenon, known in the literature as an approach toward the "empty liquid" state, is observed. The stability of the "protein droplets", formed by the liquid-liquid phase separation, depends on their local environment and temperature.

Analytical Formulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range Interactions

A theory for sequence-dependent liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) in the study of biomolecular condensates is formulated by extending the random phase approximation (RPA) and field-theoretic simulation (FTS) of heteropolymers with spatially long-range Coulomb interactions to include the fundamental effects of short-range, hydrophobic-like interactions between amino acid residues. To this end, short-range effects are modeled by Yukawa interactions between multiple nonelectrostatic charges derived from an eigenvalue decomposition of pairwise residue-residue contact energies. Chain excluded volume is afforded by incompressibility constraints. A mean-field approximation leads to an effective Flory-Huggins χ parameter, which, in conjunction with RPA, accounts for the contact-interaction effects of amino acid composition and the sequence-pattern effects of long-range electrostatics in IDP LLPS, whereas FTS based on the formulation provides full sequence dependence for both short- and long-range interactions. This general approach is illustrated here by applications to variants of a natural IDP in the context of several different amino-acid interaction schemes as well as a set of different model hydrophobic-polar sequences sharing the same composition. Effectiveness of the methodology is verified by coarse-grained explicit-chain molecular dynamics simulations.

The influence of excipients on the viscosity of monoclonal antibody solutions

The aggregation propensity of monoclonal antibodies can be modified by adding different cosolutes into the solution. A simple coarse-grained model in the combination with the thermodynamic perturbation theory was used to predict cluster distribution and viscosity of the solutions of IgG4 monoclonal anibody in the presence of L-Arginine Hydrochloride. The data were analysed using binding polynomial to describe the binding of cosolute (Arginine) to the antibody molecule. The results show that by binding to the antibody molecule the cosolute occupies some of the binding sites of the antibody, and in this way reduces the amount of binding sites available to other antibody molecules. The aggregation propensity of the antibody molecules is therefore reduced.

Phase stability of aqueous mixtures of bovine serum albumin with low molecular mass salts in presence of polyethylene glycol

The stability of bovine serum albumin (BSA) solutions against phase separation caused by cooling the system is studied under the combined influence of added poly(ethylene glycol) (PEG) and alkali halide salts in water as solvent. The phase stability of the system depends on the concentration of the added PEG and its molecular mass, the concentration of the low molecular mass electrolyte and its nature, as also on the pH of the solution. More specifically, the addition of NaCl to the BSA-PEG mixture promotes phase separation at pH = 4.0, where BSA carries the net positive charge in aqueous solution, and it increases the stability of the solution at pH=4.6, i.e., near the isoionic point of the protein. Moreover, at pH = 4.6, the cloud-point temperature decreases in the order from NaF to NaI and from LiCl to CsCl. The order of the salts at pH = 4.0 is exactly reversed: LiCl and NaF show the weakest effect on the cloud-point temperature and the strongest decrease in stability is caused by RbCl and NaNO3. An attempt is made to correlate these observations with the free energies of hydration of the added salt ions and with the effect of adsorption of salt ions on the protein surface on the protein-protein interactions. Kosmotropic salt ions decrease the phase stability of BSA-PEG-salt solutions at pH < pI, while exactly the opposite is true at pH = pI.

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.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Aggregation, liquid-liquid phase separation, and percolation behaviour of a model antibody fluid constrained by hard-sphere obstacles

This study is concerned with the behaviour of proteins within confinement created by hard-sphere obstacles. An individual antibody molecule is depicted as an assembly of seven hard spheres, organized to resemble a Y-shaped (on average) antibody (7-bead model) protein. For comparison with other studies we, in one case, model the protein as a hard sphere decorated by three short-range attractive sites. The antibody has two Fab and one Fc domains located in the corners of the letter Y. In this calculation, only the Fab-Fab and Fab-Fc attractive pair interactions are possible. The confinement is formed by the randomly distributed hard-sphere obstacles fixed in space. Aside from size exclusion, the obstacles do not interact with antibodies, but they affect the protein-protein correlation. We used a combination of the scaled-particle theory, Wertheim's thermodynamic perturbation theory and the Flory-Stockmayer theory to calculate: (i) the second virial coefficient of the protein fluid, (ii) the percolation threshold, (iii) cluster size distributions, and (iv) the liquid-liquid phase separation as a function of the strength of the various pair interactions of the protein and the model parameters, such as protein concentration and the packing fraction of obstacles. The conclusion is that hard-sphere obstacles strongly decrease the critical density and also, but to a much lesser extent, the critical temperature. Also, the confinement enhances clustering, making the percolating region broader. The effect depends on the model parameters, such as the packing fraction of obstacles η0, the inter-site interaction strength εIJ, and the ratio between the size of the obstacle σ0 and the size of one bead of the model antibody σhs; the value of this ratio is varied here from 2 to 5. Interestingly, at low to moderate packing fractions of obstacles, the second virial coefficient first slightly decreases (destabilization), and the slope depends on the observation temperature, but then at higher values of η0 it increases. The calculated values of the second virial coefficient also depend on the size of the obstacles.

Liquid network connectivity regulates the stability and composition of biomolecular condensates with many components

One of the key mechanisms used by cells to control the spatiotemporal organization of their many components is the formation and dissolution of biomolecular condensates through liquid-liquid phase separation (LLPS). Using a minimal coarse-grained model that allows us to simulate thousands of interacting multivalent proteins, we investigate the physical parameters dictating the stability and composition of multicomponent biomolecular condensates. We demonstrate that the molecular connectivity of the condensed-liquid network-i.e., the number of weak attractive protein-protein interactions per unit of volume-determines the stability (e.g., in temperature, pH, salt concentration) of multicomponent condensates, where stability is positively correlated with connectivity. While the connectivity of scaffolds (biomolecules essential for LLPS) dominates the phase landscape, introduction of clients (species recruited via scaffold-client interactions) fine-tunes it by transforming the scaffold-scaffold bond network. Whereas low-valency clients that compete for scaffold-scaffold binding sites decrease connectivity and stability, those that bind to alternate scaffold sites not required for LLPS or that have higher-than-scaffold valencies form additional scaffold-client-scaffold bridges increasing stability. Proteins that establish more connections (via increased valencies, promiscuous binding, and topologies that enable multivalent interactions) support the stability of and are enriched within multicomponent condensates. Importantly, proteins that increase the connectivity of multicomponent condensates have higher critical points as pure systems or, if pure LLPS is unfeasible, as binary scaffold-client mixtures. Hence, critical points of accessible systems (i.e., with just a few components) might serve as a unified thermodynamic parameter to predict the composition of multicomponent condensates.

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.

Biomolecular Chemistry in Liquid Phase Separated Compartments

Biochemical processes inside the cell take place in a complex environment that is highly crowded, heterogeneous, and replete with interfaces. The recently recognized importance of biomolecular condensates in cellular organization has added new elements of complexity to our understanding of chemistry in the cell. Many of these condensates are formed by liquid-liquid phase separation (LLPS) and behave like liquid droplets. Such droplet organelles can be reproduced and studied in vitro by using coacervates and have some remarkable features, including regulated assembly, differential partitioning of macromolecules, permeability to small molecules, and a uniquely crowded environment. Here, we review the main principles of biochemical organization in model membraneless compartments. We focus on some promising in vitro coacervate model systems that aptly mimic part of the compartmentalized cellular environment. We address the physicochemical characteristics of these liquid phase separated compartments, and their impact on biomolecular chemistry and assembly. These model systems enable a systematic investigation of the role of spatiotemporal organization of biomolecules in controlling biochemical processes in the cell, and they provide crucial insights for the development of functional artificial organelles and cells.

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.

Modeling the depletion effect caused by an addition of polymer to monoclonal antibody solutions

We present a theoretical study of colloidal stability of the model mixtures of monoclonal antibody molecules and non-adsorbing (no polymer-protein attraction) polymers. The antibodies are pictured as an assembly of seven hard spheres assuming a Y-like shape. Polymers present in the mixture are modeled as chain-like molecules having from 32 up to 128 monomers represented as hard spheres. We use Wertheim's thermodynamic perturbation theory to construct the two molecular species and to calculate measurable properties. The calculations are performed in the osmotic ensemble. In view that no direct attractive interaction is present in the model Hamiltonian, we only account for the entropic contribution to the phase equilibrium. We calculate chemical potentials and the equation of state for the model mixture to determine the liquid-liquid part of the phase diagram. We investigate how the critical antibody number density depends on the degree of polymerization and the bead size ratio of the polymer and protein components. The model mixture qualitatively correctly predicts some basic features of real systems. The effects of the model 'protein' geometry, that is the difference in results for the flexible Y-shaped protein versus the rigid spherical one, are also examined.

Controlling the viscosities of antibody solutions through control of their binding sites

For biotechnological drugs, it is desirable to formulate antibody solutions with low viscosities. We go beyond previous colloid theories in treating protein-protein self-association of molecules that are antibody-shaped and flexible and have spatially specific binding sites. We consider interactions either through fragment antigen (Fab-Fab) or fragment crystalizable (Fab-Fc) binding. Wertheim's theory is adapted to compute the cluster-size distributions, viscosities, second virial coefficients, and Huggins coefficients, as functions of antibody concentration. We find that the aggregation properties of concentrated solutions can be anticipated from simpler-to-measure dilute solutions. A principal finding is that aggregation is controllable, in principle, through modifying the antibody itself, and not just the solution it is dissolved in. In particular: (i) monospecific antibodies having two identical Fab arms can form linear chains with intermediate viscosities. (ii) Bispecific antibodies having different Fab arms can, in some cases, only dimerize, having low viscosities. (iii) Arm-to-Fc binding allows for three binding partners, leading to networks and high viscosities.

Theory for the Liquid-Liquid Phase Separation in Aqueous Antibody Solutions

This study presents the theory for liquid-liquid phase separation for systems of molecules modeling monoclonal antibodies. Individual molecule is depicted as an assembly of seven hard spheres, organized to mimic the Y-shaped antibody. We consider the antibody-antibody interactions either through Fab, Fab' (two Fab fragments may be different), or Fc domain. Interaction between these three domains of the molecule (hereafter denoted as A, B, and C, respectively) is modeled by a short-range square-well attraction. To obtain numerical results for the model under study, we adapt Wertheim's thermodynamic perturbation theory. We use this model to calculate the liquid-liquid phase separation curve and the second virial coefficient B2. Various interaction scenarios are examined to see how the strength of the site-site interactions and their range shape the coexistence curve. In the asymmetric case, where an attraction between two sites is favored and the interaction energies for the other sites kept constant, critical temperature first increases and than strongly decreases. Some more microscopic information, for example, the probability for the particular two sites to be connected, has been calculated. Analysis of the experimental liquid-liquid phase diagrams, obtained from literature, is presented. In addition, we calculate the second virial coefficient under conditions leading to the liquid-liquid phase separation and present this quantity on the graph B2 versus protein concentration.

Liquid-Liquid Phase Separation of Patchy Particles Illuminates Diverse Effects of Regulatory Components on Protein Droplet Formation

Recently many cellular functions have been associated with membraneless organelles, or protein droplets, formed by liquid-liquid phase separation (LLPS). Proteins in these droplets often contain RNA-binding domains, but the effects of RNA on LLPS have been controversial. To gain better understanding on the roles of RNA and other macromolecular regulators, here we used Gibbs-ensemble simulations to determine phase diagrams of two-component patchy particles, as models for mixtures of proteins with regulatory components. Protein-like particles have four patches, with attraction strength ε_PP; regulatory particles experience mutual steric repulsion but have two attractive patches toward proteins, with the strength ε_PR tunable. At low ε_PR, the regulator, due to steric repulsion, preferentially partitions in the dispersed phase, thereby displacing the protein into the droplet phase and promoting LLPS. At moderate ε_PR, the regulator starts to partition and displace the protein in the droplet phase, but only to weaken bonding networks and thereby suppress LLPS. At ε_PR > ε_PP, the enhanced bonding ability of the regulator initially promotes LLPS, but at higher amounts, the resulting displacement of the protein suppresses LLPS. These results illustrate how RNA can have disparate effects on LLPS, thus able to perform diverse functions in different organelles.

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.

Valence, loop formation and universality in self-assembling patchy particles

Patchy particles have emerged as an attractive model to mimic phase separation and self-assembly of globular proteins solutions, colloidal patchy particles, and molecular fluids where directional interactions are operative. In our previous work, we extensively explored the coupling of directional and isotropic interactions on both the phase separation and self-assembly in a system of patchy particles with five spots. Here, we extend this work to consider different patch valences and isotropic interaction strengths with an emphasis on self-assembly. Although the location of self-assembly transition lines in the temperature-density plane depend on a number of parameters, we find universal behavior of cluster size that is dependent only on the probability of a spot being bound, the patch valence, and the density. Using these principles, we quantify both the mass distribution and the shape for all clusters, as well as clusters containing loops. Following the logical implications of these results, combined with a simplified version of a mean-field theory that incorporates Flory-Stockmayer theory, we find a universal curve for the temperature dependence of cluster mass and a universal curve for the fraction of clusters that contain loops. As the curves are dependent on the particle valence, such results provide a method for parameterizing patchy particle models using experimental data.

Limiting the valence: advancements and new perspectives on patchy colloids, soft functionalized nanoparticles and biomolecules

Artistic representation of limited valance units consisting of a soft core (in blue) and a small number of flexible bonding patches (in orange).