CAMELOT: A machine learning approach for coarse-grained simulations of aggregation of block-copolymeric protein sequences

Energy landscapes for clusters of hexapeptides

We present the results for energy landscapes of hexapeptides obtained using interfaces to the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) program. We have used basin-hopping global optimization and discrete path sampling to explore the landscapes of hexapeptide monomers, dimers, and oligomers containing 10, 100, and 200 monomers modeled using a residue-level coarse-grained potential, Mpipi, implemented in LAMMPS. We find that the dimers of peptides containing amino acid residues that are better at promoting phase separation, such as tyrosine and arginine, have melting peaks at higher temperature in their heat capacity compared to phenylalanine and lysine, respectively. This observation correlates with previous work on the same uncapped hexapeptide monomers modeled using atomistic potential. For oligomers, we compare the variation in monomer conformations with radial distance and observe trends for selected angles calculated for each monomer. The LAMMPS interfaces to the GMIN and OPTIM programs for landscape exploration offer new opportunities to investigate larger systems and provide access to the coarse-grained potentials implemented within LAMMPS.

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.

Optimizing the Martini 3 Force Field Reveals the Effects of the Intricate Balance between Protein-Water Interaction Strength and Salt Concentration on Biomolecular Condensate Formation

Condensation/dissolution has become a widely acknowledged biological macromolecular assembly phenomenon in subcellular compartmentalization. The MARTINI force field offers a coarse-grained protein model with a resolution that preserves molecular details with an explicit (CG) solvent. Despite its relatively higher resolution, it can still achieve condensate formation in a reasonable computing time with explicit solvent and ionic species. Therefore, it is highly desirable to tune this force field to be able to reproduce the experimentally observed properties of the condensate formation. In this work, we studied the condensate formation of the low-sequence complexity domain of fused in sarcoma protein using a MARTINI 3 force field by systematically modifying (increasing) the protein-water interaction strength and varying the salt concentration. We found that the condensate formation is sensitive both to the protein-water interaction strength and the presence of salt. While the unmodified MARTINI force field yields a complete collapse of proteins into one dense phase (i.e., no dilute phase), we reported a range of modified protein-water interaction strength that is capable of capturing the experimentally found transfer free energy between dense and dilute phases. We also found that the condensates lose their spherical shape upon the addition of salt, especially when the protein-water interactions are weak. Interchain amino acid contact map analysis showed one explanation for this observation: the protein-protein contact fraction reduces as salt is added to systems (when the protein-water interactions are weak), consistent with electrostatic screening effects. This reduction might be responsible for the condensates becoming nonspherical upon the addition of salt by reducing the need for minimizing the interfacial area. However, as the protein-water interactions become stronger to the extent that makes the transfer free energy agree well with experimentally observed transfer free energy, we found an increase in the protein-protein contact fraction upon the addition of salt, consistent with the salting-out effects. Therefore, we concluded that there is an intricate balance between screening effects and salting-out effects upon the addition of salt and this balance is highly sensitive to the strength of protein-water interactions.

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.

Advances in Computational Approaches for Estimating Passive Permeability in Drug Discovery

Passive permeation of cellular membranes is a key feature of many therapeutics. The relevance of passive permeability spans all biological systems as they all employ biomembranes for compartmentalization. A variety of computational techniques are currently utilized and under active development to facilitate the characterization of passive permeability. These methods include lipophilicity relations, molecular dynamics simulations, and machine learning, which vary in accuracy, complexity, and computational cost. This review briefly introduces the underlying theories, such as the prominent inhomogeneous solubility diffusion model, and covers a number of recent applications. Various machine-learning applications, which have demonstrated good potential for high-volume, data-driven permeability predictions, are also discussed. Due to the confluence of novel computational methods and next-generation exascale computers, we anticipate an exciting future for computationally driven permeability predictions.

Multi-scale molecular simulation of random peptide phase separation and its extended-to-compact structure transition driven by hydrophobic interactions

Intrinsically disordered proteins (IDPs) often undergo liquid-liquid phase separation (LLPS) and form membraneless organelles or protein condensates. One of the core problems is how do electrostatic repulsion and hydrophobic interactions in peptides regulate the phase separation process? To answer this question, this study uses random peptides composed of positively charged arginine (Arg, R) and hydrophobic isoleucine (Ile, I) as the model systems, and conduct large-scale simulations using all atom and coarse-grained model multi-scale simulation methods. In this article, we investigate the phase separation of different sequences using a coarse-grained model. It is found that the stronger the electrostatic repulsion in the system, the more extended the single-chain structure, and the more likely the system forms a low-density homogeneous phase. In contrast, the stronger the hydrophobic effect of the system, the more compact the single-chain structure, the easier phase separation, and the higher the critical temperature of phase separation. Overall, by taking the random polypeptides composed of two types of amino acid residues as model systems, this study discusses the relationship between the protein sequence and phase behaviour, and provides theoretical insights into the interactions within or between proteins. It is expected to provide essential physical information for the sequence design of functional IDPs, as well as data to support the diagnosis and treatment of the LLPS-associated diseases.

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.

Intermolecular interactions underlie protein/peptide phase separation irrespective of sequence and structure at crowded milieu

Liquid-liquid phase separation (LLPS) has emerged as a crucial biological phenomenon underlying the sequestration of macromolecules (such as proteins and nucleic acids) into membraneless organelles in cells. Unstructured and intrinsically disordered domains are known to facilitate multivalent interactions driving protein LLPS. We hypothesized that LLPS could be an intrinsic property of proteins/polypeptides but with distinct phase regimes irrespective of their sequence and structure. To examine this, we studied many (a total of 23) proteins/polypeptides with different structures and sequences for LLPS study in the presence and absence of molecular crowder, polyethylene glycol (PEG-8000). We showed that all proteins and even highly charged polypeptides (under study) can undergo liquid condensate formation, however with different phase regimes and intermolecular interactions. We further demonstrated that electrostatic, hydrophobic, and H-bonding or a combination of such intermolecular interactions plays a crucial role in individual protein/peptide LLPS.

Phase separation of protein mixtures is driven by the interplay of homotypic and heterotypic interactions

Prion-like low-complexity domains (PLCDs) are involved in the formation and regulation of distinct biomolecular condensates that form via phase separation coupled to percolation. Intracellular condensates often encompass numerous distinct proteins with PLCDs. Here, we combine simulations and experiments to study mixtures of PLCDs from two RNA-binding proteins, hnRNPA1 and FUS. Using simulations and experiments, we find that 1:1 mixtures of A1-LCD and FUS-LCD undergo phase separation more readily than either of the PLCDs on their own due to complementary electrostatic interactions. Tie line analysis reveals that stoichiometric ratios of different components and their sequence-encoded interactions contribute jointly to the driving forces for condensate formation. Simulations also show that the spatial organization of PLCDs within condensates is governed by relative strengths of homotypic versus heterotypic interactions. We uncover rules for how interaction strengths and sequence lengths modulate conformational preferences of molecules at interfaces of condensates formed by mixtures of proteins.

Investigating the Hydrogen Bond-Induced Self-Assembly of Polysulfamides Using Molecular Simulations and Experiments

In this paper, we present a synergistic, experimental, and computational study of the self-assembly of N,N'-disubstituted polysulfamides driven by hydrogen bonds (H-bonds) between the H-bonding donor and acceptor groups present in repeating sulfamides as a function of the structural design of the polysulfamide backbone. We developed a coarse-grained (CG) polysulfamide model that captures the directionality of H-bonds between the sulfamide groups and used this model in molecular dynamics (MD) simulations to study the self-assembly of these polymers in implicit solvent. The CGMD approach was validated by reproducing experimentally observed trends in the extent of crystallinity for three polysulfamides synthesized with aliphatic and/or aromatic repeating units. After validation of our CGMD approach, we computationally predicted the effect of repeat unit bulkiness, length, and uniformity of segment lengths in the polymers on the extent of orientational and positional order among the self-assembled polysulfamide chains, providing key design principles for tuning the extent of crystallinity in polysulfamides in experiments. Those computational predictions were then experimentally tested through the synthesis and characterization of polysulfamide architectures.

Comparative Performance of Computer Simulation Models of Intrinsically Disordered Proteins at Different Levels of Coarse-Graining

Coarse-graining is commonly used to decrease the computational cost of simulations. However, coarse-grained models are also considered to have lower transferability, with lower accuracy for systems outside the original scope of parametrization. Here, we benchmark a bead-necklace model and a modified Martini 2 model, both coarse-grained models, for a set of intrinsically disordered proteins, with the different models having different degrees of coarse-graining. The SOP-IDP model has earlier been used for this set of proteins; thus, those results are included in this study to compare how models with different levels of coarse-graining compare. The sometimes naive expectation of the least coarse-grained model performing best does not hold true for the experimental pool of proteins used here. Instead, it showed the least good agreement, indicating that one should not necessarily trust the otherwise intuitive notion of a more advanced model inherently being better in model choice.

Crowder titrations enable the quantification of driving forces for macromolecular phase separation

Macromolecular solubility is an important contributor to the driving forces for phase separation. Formally, the driving forces in a binary mixture comprising a macromolecule dissolved in a solvent can be quantified in terms of the saturation concentration, which is the threshold macromolecular concentration above which the mixture separates into coexisting dense and dilute phases. Additionally, the second virial coefficient, which measures the effective strength of solvent-mediated intermolecular interactions provides direct assessments of solvent quality. The sign and magnitude of second virial coefficients will be governed by a combination of solution conditions and the nature of the macromolecule of interest. Here, we show, using a combination of theory, simulation, and in vitro experiments, that titrations of crowders, providing they are true depletants, can be used to extract the intrinsic driving forces for macromolecular phase separation. This refers to saturation concentrations in the absence of crowders and the second virial coefficients that quantify the magnitude of the incompatibility between macromolecules and the solvent. Our results show how the depletion-mediated attractions afforded by crowders can be leveraged to obtain comparative assessments of macromolecule-specific, intrinsic driving forces for phase separation.

SIGNIFICANCE

Phase separation has emerged as a process of significant relevance to sorting macromolecules into distinct compartments, thereby enabling spatial and temporal control over cellular matter. Considerable effort is being invested into uncovering the driving forces that enable the separation of macromolecular solutions into coexisting phases. At its heart, this process is governed by the balance of macromolecule-solvent, inter-macromolecule, and solvent-solvent interactions. We show that the driving forces for phase separation, including the coefficients that measure interaction strengths between macromolecules, can be extracted by titrating the concentrations of crowders that enable macromolecules to phase separate at lower concentrations. Our work paves the way to leverage specific categories of measurements for quantitative characterizations of driving forces for phase separation.

Thermodynamic origins of two-component multiphase condensates of proteins

Intracellular condensates are highly multi-component systems in which complex phase behaviour can ensue, including the formation of architectures comprising multiple immiscible condensed phases. Relying solely on physical intuition to manipulate such condensates is difficult because of the complexity of their composition, and systematically learning the underlying rules experimentally would be extremely costly. We address this challenge by developing a computational approach to design pairs of protein sequences that result in well-separated multilayered condensates and elucidate the molecular origins of these compartments. Our method couples a genetic algorithm to a residue-resolution coarse-grained protein model. We demonstrate that we can design protein partners to form multiphase condensates containing naturally occurring proteins, such as the low-complexity domain of hnRNPA1 and its mutants, and show how homo- and heterotypic interactions must differ between proteins to result in multiphasicity. We also show that in some cases the specific pattern of amino-acid residues plays an important role. Our findings have wide-ranging implications for understanding and controlling the organisation, functions and material properties of biomolecular condensates.

Condensates formed by prion-like low-complexity domains have small-world network structures and interfaces defined by expanded conformations

Biomolecular condensates form via coupled associative and segregative phase transitions of multivalent associative macromolecules. Phase separation coupled to percolation is one example of such transitions. Here, we characterize molecular and mesoscale structural descriptions of condensates formed by intrinsically disordered prion-like low complexity domains (PLCDs). These systems conform to sticker-and-spacers architectures. Stickers are cohesive motifs that drive associative interactions through reversible crosslinking and spacers affect the cooperativity of crosslinking and overall macromolecular solubility. Our computations reproduce experimentally measured sequence-specific phase behaviors of PLCDs. Within simulated condensates, networks of reversible inter-sticker crosslinks organize PLCDs into small-world topologies. The overall dimensions of PLCDs vary with spatial location, being most expanded at and preferring to be oriented perpendicular to the interface. Our results demonstrate that even simple condensates with one type of macromolecule feature inhomogeneous spatial organizations of molecules and interfacial features that likely prime them for biochemical activity.

Developing Bonded Potentials for a Coarse-Grained Model of Intrinsically Disordered Proteins

Recent advances in residue-level coarse-grained (CG) computational models have enabled molecular-level insights into biological condensates of intrinsically disordered proteins (IDPs), shedding light on the sequence determinants of their phase separation. The existing CG models that treat protein chains as flexible molecules connected via harmonic bonds cannot populate common secondary-structure elements. Here, we present a CG dihedral angle potential between four neighboring beads centered at Cα atoms to faithfully capture the transient helical structures of IDPs. In order to parameterize and validate our new model, we propose Cα-based helix assignment rules based on dihedral angles that succeed in reproducing the atomistic helicity results of a polyalanine peptide and folded proteins. We then introduce sequence-dependent dihedral angle potential parameters (εd) and use experimentally available helical propensities of naturally occurring 20 amino acids to find their optimal values. The single-chain helical propensities from the CG simulations for commonly studied prion-like IDPs are in excellent agreement with the NMR-based α-helix fraction, demonstrating that the new HPS-SS model can accurately produce structural features of IDPs. Furthermore, this model can be easily implemented for large-scale assembly simulations due to its simplicity.

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.

Deciphering how naturally occurring sequence features impact the phase behaviours of disordered prion-like domains

Prion-like low-complexity domains (PLCDs) have distinctive sequence grammars that determine their driving forces for phase separation. Here we uncover the physicochemical underpinnings of how evolutionarily conserved compositional biases influence the phase behaviour of PLCDs. We interpret our results in the context of the stickers-and-spacers model for the phase separation of associative polymers. We find that tyrosine is a stronger sticker than phenylalanine, whereas arginine is a context-dependent auxiliary sticker. In contrast, lysine weakens sticker-sticker interactions. Increasing the net charge per residue destabilizes phase separation while also weakening the strong coupling between single-chain contraction in dilute phases and multichain interactions that give rise to phase separation. Finally, glycine and serine residues act as non-equivalent spacers, and thus make the glycine versus serine contents an important determinant of the driving forces for phase separation. The totality of our results leads to a set of rules that enable comparative estimates of composition-specific driving forces for PLCD phase separation.

Physics-driven coarse-grained model for biomolecular phase separation with near-quantitative accuracy

Various physics- and data-driven sequence-dependent protein coarse-grained models have been developed to study biomolecular phase separation and elucidate the dominant physicochemical driving forces. Here, we present Mpipi, a multiscale coarse-grained model that describes almost quantitatively the change in protein critical temperatures as a function of amino-acid sequence. The model is parameterised from both atomistic simulations and bioinformatics data and accounts for the dominant role of π-π and hybrid cation-π/π-π interactions and the much stronger attractive contacts established by arginines than lysines. We provide a comprehensive set of benchmarks for Mpipi and seven other residue-level coarse-grained models against experimental radii of gyration and quantitative in-vitro phase diagrams; Mpipi predictions agree well with experiment on both fronts. Moreover, it can account for protein-RNA interactions, correctly predicts the multiphase behaviour of a charge-matched poly-arginine/poly-lysine/RNA system, and recapitulates experimental LLPS trends for sequence mutations on FUS, DDX4 and LAF-1 proteins.

Integrating single-molecule spectroscopy and simulations for the study of intrinsically disordered proteins

Over the last two decades, intrinsically disordered proteins and protein regions (IDRs) have emerged from a niche corner of biophysics to be recognized as essential drivers of cellular function. Various techniques have provided fundamental insight into the function and dysfunction of IDRs. Among these techniques, single-molecule fluorescence spectroscopy and molecular simulations have played a major role in shaping our modern understanding of the sequence-encoded conformational behavior of disordered proteins. While both techniques are frequently used in isolation, when combined they offer synergistic and complementary information that can help uncover complex molecular details. Here we offer an overview of single-molecule fluorescence spectroscopy and molecular simulations in the context of studying disordered proteins. We discuss the various means in which simulations and single-molecule spectroscopy can be integrated, and consider a number of studies in which this integration has uncovered biological and biophysical mechanisms.

On the potential of machine learning to examine the relationship between sequence, structure, dynamics and function of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) constitute a broad set of proteins with few uniting and many diverging properties. IDPs-and intrinsically disordered regions (IDRs) interspersed between folded domains-are generally characterized as having no persistent tertiary structure; instead they interconvert between a large number of different and often expanded structures. IDPs and IDRs are involved in an enormously wide range of biological functions and reveal novel mechanisms of interactions, and while they defy the common structure-function paradigm of folded proteins, their structural preferences and dynamics are important for their function. We here discuss open questions in the field of IDPs and IDRs, focusing on areas where machine learning and other computational methods play a role. We discuss computational methods aimed to predict transiently formed local and long-range structure, including methods for integrative structural biology. We discuss the many different ways in which IDPs and IDRs can bind to other molecules, both via short linear motifs, as well as in the formation of larger dynamic complexes such as biomolecular condensates. We discuss how experiments are providing insight into such complexes and may enable more accurate predictions. Finally, we discuss the role of IDPs in disease and how new methods are needed to interpret the mechanistic effects of genomic variants in IDPs.

Targeted modulation of protein liquid-liquid phase separation by evolution of amino-acid sequence

Rationally and efficiently modifying the amino-acid sequence of proteins to control their ability to undergo liquid-liquid phase separation (LLPS) on demand is not only highly desirable, but can also help to elucidate which protein features are important for LLPS. Here, we propose a computational method that couples a genetic algorithm to a sequence-dependent coarse-grained protein model to evolve the amino-acid sequences of phase-separating intrinsically disordered protein regions (IDRs), and purposely enhance or inhibit their capacity to phase-separate. We validate the predicted critical solution temperatures of the mutated sequences with ABSINTH, a more accurate all-atom model. We apply the algorithm to the phase-separating IDRs of three naturally occurring proteins, namely FUS, hnRNPA1 and LAF1, as prototypes of regions that exist in cells and undergo homotypic LLPS driven by different types of intermolecular interaction, and we find that the evolution of amino-acid sequences towards enhanced LLPS is driven in these three cases, among other factors, by an increase in the average size of the amino acids. However, the direction of change in the molecular driving forces that enhance LLPS (such as hydrophobicity, aromaticity and charge) depends on the initial amino-acid sequence. Finally, we show that the evolution of amino-acid sequences to modulate LLPS is strongly coupled to the make-up of the medium (e.g. the presence or absence of RNA), which may have significant implications for our understanding of phase separation within the many-component mixtures of biological systems.

Size conservation emerges spontaneously in biomolecular condensates formed by scaffolds and surfactant clients

Biomolecular condensates are liquid-like membraneless compartments that contribute to the spatiotemporal organization of proteins, RNA, and other biomolecules inside cells. Some membraneless compartments, such as nucleoli, are dispersed as different condensates that do not grow beyond a certain size, or do not present coalescence over time. In this work, using a minimal protein model, we show that phase separation of binary mixtures of scaffolds and low-valency clients that can act as surfactants-i.e., that significantly reduce the droplet surface tension-can yield either a single drop or multiple droplets that conserve their sizes on long timescales (herein 'multidroplet size-conserved' scenario'), depending on the scaffold to client ratio. Our simulations demonstrate that protein connectivity and condensate surface tension regulate the balance between these two scenarios. The multidroplet size-conserved scenario spontaneously arises at increasing surfactant-to-scaffold concentrations, when the interfacial penalty for creating small liquid droplets is sufficiently reduced by the surfactant proteins that are preferentially located at the interface. In contrast, low surfactant-to-scaffold concentrations enable continuous growth and fusion of droplets without restrictions. Overall, our work proposes one thermodynamic mechanism to help rationalize how size-conserved coexisting condensates can persist inside cells-shedding light on the roles of protein connectivity, binding affinity, and droplet composition in this process.

"Dividing and Conquering" and "Caching" in Molecular Modeling

Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.

Unfolding of the chromatin fiber driven by overexpression of noninteracting bridging factors

Nuclear molecules control the functional properties of the chromatin fiber by shaping its morphological properties. The biophysical mechanisms controlling how bridging molecules compactify chromatin are a matter of debate. On the one side, bridging molecules could cross-link faraway sites and fold the fiber through the formation of loops. Interacting bridging molecules could also mediate long-range attractions by first tagging different locations of the fiber and then undergoing microphase separation. Using a coarse-grained model and Monte Carlo simulations, we study the conditions leading to compact configurations both for interacting and noninteracting bridging molecules. In the second case, we report on an unfolding transition at high densities of the bridging molecules. We clarify how this transition, which disappears for interacting bridging molecules, is universal and controlled by entropic terms. In general, chains are more compact in the case of interacting bridging molecules because interactions are not valence limited. However, this result is conditional on the ability of our simulation methodology to relax the system toward its ground state. In particular, we clarify how, unless using reaction dynamics that change the length of a loop in a single step, the system is prone to remain trapped in metastable, compact configurations featuring long loops.

Valency and Binding Affinity Variations Can Regulate the Multilayered Organization of Protein Condensates with Many Components

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

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.

Simulation of FUS Protein Condensates with an Adapted Coarse-Grained Model

Disordered proteins and nucleic acids can condense into droplets that resemble the membraneless organelles observed in living cells. MD simulations offer a unique tool to characterize the molecular interactions governing the formation of these biomolecular condensates, their physicochemical properties, and the factors controlling their composition and size. However, biopolymer condensation depends sensitively on the balance between different energetic and entropic contributions. Here, we develop a general strategy to fine-tune the potential energy function for molecular dynamics simulations of biopolymer phase separation. We rebalance protein-protein interactions against solvation and entropic contributions to match the excess free energy of transferring proteins between dilute solution and condensate. We illustrate this formalism by simulating liquid droplet formation of the FUS low-complexity domain (LCD) with a rebalanced MARTINI model. By scaling the strength of the nonbonded interactions in the coarse-grained MARTINI potential energy function, we map out a phase diagram in the plane of protein concentration and interaction strength. Above a critical scaling factor of αc ≈ 0.6, FUS-LCD condensation is observed, where α = 1 and 0 correspond to full and repulsive interactions in the MARTINI model. For a scaling factor α = 0.65, we recover experimental densities of the dilute and dense phases, and thus the excess protein transfer free energy into the droplet and the saturation concentration where FUS-LCD condenses. In the region of phase separation, we simulate FUS-LCD droplets of four different sizes in stable equilibrium with the dilute phase and slabs of condensed FUS-LCD for tens of microseconds, and over one millisecond in aggregate. We determine surface tensions in the range of 0.01-0.4 mN/m from the fluctuations of the droplet shape and from the capillary-wave-like broadening of the interface between the two phases. From the dynamics of the protein end-to-end distance, we estimate shear viscosities from 0.001 to 0.02 Pa s for the FUS-LCD droplets with scaling factors α in the range of 0.625-0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.

Artificial intelligence techniques for integrative structural biology of intrinsically disordered proteins

We outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self-organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure-function relationships-however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.

Machine Learning Predictions of Block Copolymer Self-Assembly

Directed self-assembly of block copolymers is a key enabler for nanofabrication of devices with sub-10 nm feature sizes, allowing patterning far below the resolution limit of conventional photolithography. Among all the process steps involved in block copolymer self-assembly, solvent annealing plays a dominant role in determining the film morphology and pattern quality, yet the interplay of the multiple parameters during solvent annealing, including the initial thickness, swelling, time, and solvent ratio, makes it difficult to predict and control the resultant self-assembled pattern. Here, machine learning tools are applied to analyze the solvent annealing process and predict the effect of process parameters on morphology and defectivity. Two neural networks are constructed and trained, yielding accurate prediction of the final morphology in agreement with experimental data. A ridge regression model is constructed to identify the critical parameters that determine the quality of line/space patterns. These results illustrate the potential of machine learning to inform nanomanufacturing processes.

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.

Connecting Coil-to-Globule Transitions to Full Phase Diagrams for Intrinsically Disordered Proteins

Phase separation is thought to underlie spatial and temporal organization that is required for controlling biochemical reactions in cells. Multivalence of interaction motifs, also known as stickers, is a defining feature of proteins that drive phase separation. Intrinsically disordered proteins with stickers uniformly distributed along the linear sequence can serve as scaffold molecules that drive phase separation. The sequence-intrinsic contributions of disordered proteins to phase separation can be discerned by computing or measuring sequence-specific phase diagrams. These help to delineate the combinations of protein concentration and a suitable control parameter, such as temperature, that support phase separation. Here, we present an approach that combines detailed simulations with a numerical adaptation of an analytical Gaussian cluster theory to enable the calculation of sequence-specific phase diagrams. Our approach leverages the known equivalence between the driving forces for single-chain collapse in dilute solutions and the driving forces for phase separation in concentrated solutions. We demonstrate the application of the theory-aided computations through calculation of phase diagrams for a set of archetypal intrinsically disordered low-complexity domains. We also leverage theories to compute sequence-specific percolation lines and thereby provide a thermodynamic framework for hardening transitions that have been observed for many biomolecular condensates.

Restricting the sizes of condensates

Computer simulations of model proteins with sticker-and-spacer architectures shed light on the formation of biomolecular condensates in cells.

Information theoretic measures for quantifying sequence-ensemble relationships of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) contribute to a multitude of functions. De novo design of IDPs should open the door to modulating functions and phenotypes controlled by these systems. Recent design efforts have focused on compositional biases and specific sequence patterns as the design features. Analysis of the impact of these designs on sequence-function relationships indicates that individual sequence/compositional parameters are insufficient for describing sequence-function relationships in IDPs. To remedy this problem, we have developed information theoretic measures for sequence–ensemble relationships (SERs) of IDPs. These measures rely on prior availability of statistically robust conformational ensembles derived from all atom simulations. We show that the measures we have developed are useful for comparing sequence-ensemble relationships even when sequence is poorly conserved. Based on our results, we propose that de novo designs of IDPs, guided by knowledge of their SERs, should provide improved insights into their sequence–ensemble–function relationships.

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.

Thermodynamics of Huntingtin Aggregation

Amyloid aggregates are found in many neurodegenerative diseases, including Huntington's, Alzheimer's, and prion diseases. The precise role of the aggregates in disease progression has been difficult to elucidate because of the diversity of aggregated states they can adopt. Here, we study the formation of fibrils and oligomers by exon 1 of huntingtin protein. We show that the oligomer states are consistent with polymer micelles that are limited in size by the stretching entropy of the polyglutamine region. The model shows how the sequences flanking the amyloid core modulate aggregation behavior. The N17 region promotes aggregation through weakly attractive interactions, whereas the C38 tail opposes aggregation via steric repulsion. We also show that the energetics of cross-β stacking by polyglutamine would produce fibrils with many alignment defects, but minor perturbations from the flanking sequences are sufficient to reduce the defects to the level observed in experiment. We conclude with a discussion of the implications of this model for other amyloid-forming molecules.

The pathophysiology of neurodegenerative disease: Disturbing the balance between phase separation and irreversible aggregation

Liquid-liquid phase separation (LLPS) brings together functionally related proteins through the intrinsic biophysics of proteins in a process that is driven by reducing free energy and maximizing entropy. The process of LLPS allows proteins to form structures, termed membrane-less organelles. These diverse, dynamic organelles are active in a wide range of processes in the nucleus, cytoplasm, mitochondria and synapse, and ranging from bacteria to plants to eukaryotes. RNA and DNA present long chained charged polymers that promote LLPS. Consequently, many RNA binding proteins (RBPs) and DNA binding proteins form membrane-less organelles. However, the highly concentrated phase separated state creates conditions that also promote formation of irreversible protein aggregates. Mutations in RNA and DNA binding proteins that increase the stability of irreversible aggregates also increase the accumulation of irreversible aggregates directly and from membrane-less organelles. Many of the RBPs that exhibit disease-linked mutations carry out cytoplasmic actions through stress granules, which are a pleiotropic type of RNA granule that regulates the translational response to stress. Phosphorylation and oligomerization of tau facilitates its interactions with RBPs and ribosomal proteins, affecting RNA translation; we propose that this is a major reason that tau becomes phosphorylated with stress. Persistent stress leads to the accumulation of irreversible aggregates composed of RBPs or tau, which then cause toxicity and form many of the hallmark pathologies of major neurodegenerative diseases. This pathophysiology ultimately leads to multiple forms of neurodegenerative diseases, the specific type of which reflects the temporal and spatial accumulation of different aggregating proteins.

Physical Principles Underlying the Complex Biology of Intracellular Phase Transitions

Many biomolecular condensates appear to form via spontaneous or driven processes that have the hallmarks of intracellular phase transitions. This suggests that a common underlying physical framework might govern the formation of functionally and compositionally unrelated biomolecular condensates. In this review, we summarize recent work that leverages a stickers-and-spacers framework adapted from the field of associative polymers for understanding how multivalent protein and RNA molecules drive phase transitions that give rise to biomolecular condensates. We discuss how the valence of stickers impacts the driving forces for condensate formation and elaborate on how stickers can be distinguished from spacers in different contexts. We touch on the impact of sticker- and spacer-mediated interactions on the rheological properties of condensates and show how the model can be mapped to known drivers of different types of biomolecular condensates.

Considerations and Challenges in Studying Liquid-Liquid Phase Separation and Biomolecular Condensates

Evidence is now mounting that liquid-liquid phase separation (LLPS) underlies the formation of membraneless compartments in cells. This realization has motivated major efforts to delineate the function of such biomolecular condensates in normal cells and their roles in contexts ranging from development to age-related disease. There is great interest in understanding the underlying biophysical principles and the specific properties of biological condensates with the goal of bringing insights into a wide range of biological processes and systems. The explosion of physiological and pathological contexts involving LLPS requires clear standards for their study. Here, we propose guidelines for rigorous experimental characterization of LLPS processes in vitro and in cells, discuss the caveats of common experimental approaches, and point out experimental and theoretical gaps in the field.

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.

Utilizing Machine Learning for Efficient Parameterization of Coarse Grained Molecular Force Fields

We present a machine learning approach to automated force field development in dissipative particle dynamics (DPD). The approach employs Bayesian optimization to parametrize a DPD force field against experimentally determined partition coefficients. The optimization process covers a discrete space of over 40 000 000 points, where each point represents the set of potentials that jointly forms a force field. We find that Bayesian optimization is capable of reaching a force field of comparable performance to the current state-of-the-art within 40 iterations. The best iteration during the optimization achieves an R² of 0.78 and an RMSE of 0.63 log units on the training set of data, these metrics are maintained when a validation set is included, giving R² of 0.8 and an RMSE of 0.65 log units. This work hence provides a proof-of-concept, expounding the utility of coupling automated and efficient global optimization with a top down data driven approach to force field parametrization. Compared to commonly employed alternative methods, Bayesian optimization offers global parameter searching and a low time to solution.

Three archetypical classes of macromolecular regulators of protein liquid-liquid phase separation

Membraneless organelles, corresponding to the droplet phase upon liquid-liquid phase separation (LLPS) of protein or protein-RNA mixtures, mediate myriad cellular functions. Cells use a variety of biochemical signals such as expression level and posttranslational modification to regulate droplet formation and dissolution, but the physical basis of the regulatory mechanisms remains ill-defined and quantitative assessment of the effects is largely lacking. Our computational study predicted that the strength of attraction by droplet-forming proteins dictates whether and how macromolecular regulators promote or suppress LLPS. We experimentally tested this prediction, using the pentamers of SH3 domains and proline-rich motifs (SH3₅ and PRM₅) as droplet-forming proteins. Determination of the changes in phase boundary and the partition coefficients in the droplet phase over a wide range of regulator concentrations yielded both a quantitative measure and a mechanistic understanding of the regulatory effects. Three archetypical classes of regulatory effects were observed. Ficoll 70 at high concentrations indirectly promoted SH3₅-PRM₅ LLPS, by taking up volume in the bulk phase and thereby displacing SH3₅ and PRM₅ into the droplet phase. Lysozyme had a moderate partition coefficient and suppressed LLPS by substituting weaker attraction with SH3₅ for the stronger SH3₅-PRM₅ attraction in the droplet phase. By forming even stronger attraction with PRM₅, heparin at low concentrations partitioned heavily into the droplet phase and promoted LLPS. These characteristics were recapitulated by computational results of patchy particle models, validating the identification of the 3 classes of macromolecular regulators as volume-exclusion promotors, weak-attraction suppressors, and strong-attraction promotors.

Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review

The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.

Controlling compartmentalization by non-membrane-bound organelles

Compartmentalization is a characterizing feature of complexity in cells, used to organize their biochemistry. Membrane-bound organelles are most widely known, but non-membrane-bound liquid organelles also exist. These have recently been shown to form by phase separation of specific types of proteins known as scaffolds. This forms two phases: a condensate that is enriched in scaffold protein separated by a phase boundary from the cytoplasm or nucleoplasm with a low concentration of the scaffold protein. Phase separation is well known for synthetic polymers, but also appears important in cells. Here, we review the properties of proteins important for forming these non-membrane-bound organelles, focusing on the energetically favourable interactions that drive condensation. On this basis we make qualitative predictions about how cells may control compartmentalization by condensates; the partition of specific molecules to a condensate; the control of condensation and dissolution of condensates; and the regulation of condensate nucleation. There are emerging data supporting many of these predictions, although future results may prove incorrect. It appears that many molecules may have the ability to modulate condensate formation, making condensates a potential target for future therapeutics. The emerging properties of condensates are fundamentally unlike the properties of membrane-bound organelles. They have the capacity to rapidly integrate cellular events and act as a new class of sensors for internal and external environments.

This article is part of the theme issue ‘Self-organization in cell biology’.

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.

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.

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.

Phase Separation of Intrinsically Disordered Proteins

There is growing interest in the topic of intracellular phase transitions that lead to the formation of biologically regulated biomolecular condensates. These condensates are membraneless bodies formed by phase separation of key protein and nucleic acid molecules from the cytoplasmic or nucleoplasmic milieus. The drivers of phase separation are referred to as scaffolds whereas molecules that preferentially partition into condensates formed by scaffolds are known as clients. Recent advances have shown that it is possible to generate physical and functional facsimiles of many biomolecular condensates in vitro. This is achieved by titrating the concentration of key scaffold proteins and solution parameters such as salt concentration, pH, or temperature. The ability to reproduce phase separation in vitro allows one to compare the relationships between information encoded in the sequences of scaffold proteins and the driving forces for phase separation. Many scaffold proteins include intrinsically disordered regions whereas others are entirely disordered. Our focus is on comparative assessments of phase separation for different scaffold proteins, specifically intrinsically disordered linear multivalent proteins. We highlight the importance of coexistence curves known as binodals for quantifying phase behavior and comparing driving forces for sequence-specific phase separation. We describe the information accessible from full binodals and highlight different methods for-and challenges associated with-mapping binodals. In essence, we provide a wish list for in vitro characterization of phase separation of intrinsically disordered proteins. Fulfillment of this wish list through key advances in experiment, computation, and theory should bring us closer to being able to predict in vitro phase behavior for scaffold proteins and connect this to the functions and features of biomolecular condensates.