Model for disordered proteins with strongly sequence-dependent liquid phase behavior

Effects of charge asymmetry on the liquid-liquid phase separation of polyampholytes and their condensate properties

Liquid-liquid phase separation (LLPS) is the mechanism underlying the formation of bio-molecular condensates which are important compartments regulating intra- and extra-cellular functions. Electrostatic interactions are some of the important driving forces of the LLPS behaviors of biomolecules. However, the understanding of the electrostatic interactions is still limited, especially in the mixtures of biomolecules with different charge patterns. Here, we focus on the electrostatic interactions in mixtures of charge-asymmetric and charge-symmetric polyampholytes and their roles in the phase separation behaviors. We build charge-asymmetric and charge-symmetric model proteins consisting of both glutamic acid (E, negatively charged) and lysine (K, positively charged), i.e. polyampholytes of E35K15 (charge asymmetric) and E25K25 (charge symmetric). Pure E25K25 can undergo LLPS. To investigate the effects of charge-asymmetric polyampholytes on the mixtures of E25K25/E35K15, we perform coarse-grained simulations to determine their phase separation. The charge-asymmetric polyampholyte E35K15 is resistant to the LLPS of the mixtures of E25K25/E35K15. The condensate density decreases with the molar fraction of E35K15 increasing to 0.4, and no LLPS occurs at the molar fraction of 0.5 and above. This can be attributed to the electrostatic repulsion between the negatively charged E35K15 polymers. We further investigate the effects of charge asymmetry on the conformations and properties of the condensates. The E35K15 polymers in the condensates exhibit a more collapsed state as the molar fraction of E35K15 increases. However, the conformation of E25K25 polymers changes slightly across different condensates. The surface tensions of condensates decline with the increase of the molar fraction of E35K15 polymers, while the diffusivity of polymers in the condensed phases is enhanced. This work elucidates the role of charge-asymmetric polyampholytes in determining the LLPS behaviours of binary mixtures of charge-symmetric and charge-asymmetric proteins as well as the properties of condensed phases.

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.

Asymmetric oligomerization state and sequence patterning can tune multiphase condensate miscibility

Endogenous biomolecular condensates, composed of a multitude of proteins and RNAs, can organize into multiphasic structures with compositionally distinct phases. This multiphasic organization is generally understood to be critical for facilitating their proper biological function. However, the biophysical principles driving multiphase formation are not completely understood. Here we use in vivo condensate reconstitution experiments and coarse-grained molecular simulations to investigate how oligomerization and sequence interactions modulate multiphase organization in biomolecular condensates. We demonstrate that increasing the oligomerization state of an intrinsically disordered protein results in enhanced immiscibility and multiphase formation. Interestingly, we find that oligomerization tunes the miscibility of intrinsically disordered proteins in an asymmetric manner, with the effect being more pronounced when the intrinsically disordered protein, exhibiting stronger homotypic interactions, is oligomerized. Our findings suggest that oligomerization is a flexible biophysical mechanism that cells can exploit to tune the internal organization of biomolecular condensates and their associated biological functions.

Sequence dependence of critical properties for two-letter chains

Histogram-reweighting grand canonical Monte Carlo simulations are used to obtain the critical properties of lattice chains composed of solvophilic and solvophobic monomers. The model is a modification of one proposed by Larson et al. [J. Chem. Phys. 83, 2411 (1985)], lowering the "contrast" between beads of different types to prevent aggregation into finite-size micelles that would mask true phase separation between bulk high- and low-density phases. Oligomeric chains of lengths between 5 and 24 beads are studied. Mixed-field finite-size scaling methods are used to obtain the critical properties with typical relative accuracies of better than 10-4 for the critical temperature and 10-3 for the critical volume fraction. Diblock chains are found to have lower critical temperatures and volume fractions relative to the corresponding homopolymers. The addition of solvophilic blocks of increasing length to a fixed-length solvophobic segment results in a decrease of both the critical temperature and the critical volume fraction, with an eventual slow asymptotic approach to the long-chain limiting behavior. Moving a single solvophobic or solvophilic bead along a chain leads to a minimum or maximum in the critical temperature, with no change in the critical volume fraction. Chains of identical length and composition have a significant spread in their critical properties, depending on their precise sequence. The present study has implications for understanding biomolecular phase separation and for developing design rules for synthetic polymers with specific phase separation properties. It also provides data potentially useful for the further development of theoretical models for polymer and surfactant phase behavior.

Investigating neuropathological changes and underlying neurobiological mechanisms in the early stages of primary blast-induced traumatic brain injury: Insights from a rat model

The utilization of explosives and chemicals has resulted in a rise in blast-induced traumatic brain injury (bTBI) in recent times. However, there is a dearth of diagnostic biomarkers and therapeutic targets for bTBI due to a limited understanding of biological mechanisms, particularly in the early stages. The objective of this study was to examine the early neuropathological characteristics and underlying biological mechanisms of primary bTBI. A total of 83 Sprague Dawley rats were employed, with their heads subjected to a blast shockwave of peak overpressure ranging from 172 to 421 kPa in the GI, GII, and GIII groups within a closed shock tube, while the body was shielded. Neuromotor dysfunctions, morphological changes, and neuropathological alterations were detected through modified neurologic severity scores, brain water content analysis, MRI scans, histological, TUNEL, and caspase-3 immunohistochemical staining. In addition, label-free quantitative (LFQ)-proteomics was utilized to investigate the biological mechanisms associated with the observed neuropathology. Notably, no evident damage was discernible in the GII and GI groups, whereas mild brain injury was observed in the GIII group. Neuropathological features of bTBI were characterized by morphologic changes, including neuronal injury and apoptosis, cerebral edema, and cerebrovascular injury in the shockwave's path. Subsequently, 3153 proteins were identified and quantified in the GIII group, with subsequent enriched neurological responses consistent with pathological findings. Further analysis revealed that signaling pathways such as relaxin signaling, hippo signaling, gap junction, chemokine signaling, and sphingolipid signaling, as well as hub proteins including Prkacb, Adcy5, and various G-protein subunits (Gnai2, Gnai3, Gnao1, Gnb1, Gnb2, Gnb4, and Gnb5), were closely associated with the observed neuropathology. The expression of hub proteins was confirmed via Western blotting. Accordingly, this study proposes signaling pathways and key proteins that exhibit sensitivity to brain injury and are correlated with the early pathologies of bTBI. Furthermore, it highlights the significance of G-protein subunits in bTBI pathophysiology, thereby establishing a theoretical foundation for early diagnosis and treatment strategies for primary bTBI.

Interfacial exchange dynamics of biomolecular condensates are highly sensitive to client interactions

Phase separation of biomolecules can facilitate their spatiotemporally regulated self-assembly within living cells. Due to the selective yet dynamic exchange of biomolecules across condensate interfaces, condensates can function as reactive hubs by concentrating enzymatic components for faster kinetics. The principles governing this dynamic exchange between condensate phases, however, are poorly understood. In this work, we systematically investigate the influence of client-sticker interactions on the exchange dynamics of protein molecules across condensate interfaces. We show that increasing affinity between a model protein scaffold and its client molecules causes the exchange of protein chains between the dilute and dense phases to slow down and that beyond a threshold interaction strength, this slowdown in exchange becomes substantial. Investigating the impact of interaction symmetry, we found that chain exchange dynamics are also considerably slower when client molecules interact equally with different sticky residues in the protein. The slowdown of exchange is due to a sequestration effect, by which there are fewer unbound stickers available at the interface to which dilute phase chains may attach. These findings highlight the fundamental connection between client-scaffold interaction networks and condensate exchange dynamics.

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations

Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.

Molecular Dynamics Simulations of Polyelectrolyte Complexes

Polyelectrolyte complexes (PECs) are currently of great interest due to their applications toward developing new adaptive materials and their relevance in membraneless organelles. These complexes emerge during phase separation when oppositely charged polymers are mixed in aqueous media. Peptide-based PECs are particularly useful toward developing new drug delivery methods due to their inherent biocompatibility. The underlying peptide sequence can be tuned to optimize specific material properties of the complex, such as interfacial tension and viscosity. Given their applicability, it would be advantageous to understand the underlying sequence-dependent phase behavior of oppositely charged peptides. Here, we report microsecond molecular dynamic simulations to characterize the effect of hydrophobicity on the sequence-dependent peptide conformation for model polypeptide sequences that were previously reported by Tabandeh et al. These sequences are designed with alternating chirality of the peptide backbone. We present microsecond simulations of six oppositely charged peptide pairs, characterizing the sequence-dependent effect on peptide size, degree of hydrogen bonding, secondary structure, and conformation. This analysis recapitulates sensible trends in peptide conformation and degree of hydrogen bonding, consistent with experimentally reported results. Ramachandran plots reveal that backbone conformation at the single amino acid level is highly influenced by the neighboring sequence in the chain. These results give insight into how subtle changes in hydrophobic side chain size and chirality influence the strength of hydrogen bonding between the chains and, ultimately, the secondary structure. Furthermore, principal component analysis reveals that the minimum energy structures may be subtly modulated by the underlying sequence.

Interplay between posttranslational modifications and liquid‒liquid phase separation in tumors

Liquid‒liquid phase separation (LLPS) is a general phenomenon recently recognized to be critically involved in the regulation of a variety of cellular biological processes, such as transcriptional regulation, heterochromatin formation and signal transduction, through the compartmentalization of proteins or nucleic acids into droplet-like condensates. These processes are directly or indirectly related to tumor initiation and treatment. Posttranslational modifications (PTMs), which represent a rapid and reversible mechanism involved in the functional regulation of proteins, have emerged as key events in modulating LLPS under physiological or pathophysiological conditions, including tumorigenesis and antitumor therapy. In this review, we introduce the biological functions participated in cancer-associated LLPS, discuss the potential roles of LLPS during tumor onset or therapy, and emphasize the mechanistic characteristics of LLPS regulated by PTMs and its effects on tumor progression. We then provide a perspective on further studies on LLPS and its regulation by PTMs in cancer research. This review aims to broaden the understanding of the functions of LLPS and its regulation by PTMs under normal or aberrant cellular conditions.

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.

Active learning of the thermodynamics-dynamics trade-off in protein condensates

Phase-separated biomolecular condensates exhibit a wide range of dynamic properties, which depend on the sequences of the constituent proteins and RNAs. However, it is unclear to what extent condensate dynamics can be tuned without also changing the thermodynamic properties that govern phase separation. Using coarse-grained simulations of intrinsically disordered proteins, we show that the dynamics and thermodynamics of homopolymer condensates are strongly correlated, with increased condensate stability being coincident with low mobilities and high viscosities. We then apply an "active learning" strategy to identify heteropolymer sequences that break this correlation. This data-driven approach and accompanying analysis reveal how heterogeneous amino acid compositions and nonuniform sequence patterning map to a range of independently tunable dynamic and thermodynamic properties of biomolecular condensates. Our results highlight key molecular determinants governing the physical properties of biomolecular condensates and establish design rules for the development of stimuli-responsive biomaterials.

Through the lens of phase separation: intrinsically unstructured protein and chromatin looping

The establishment, maintenance and dynamic regulation of three-dimensional (3D) chromatin structures provide an important means for partitioning of genome into functionally distinctive domains, which helps to define specialized gene expression programs associated with developmental stages and cell types. Increasing evidence supports critical roles for intrinsically disordered regions (IDRs) harbored within transcription factors (TFs) and chromatin-modulatory proteins in inducing phase separation, a phenomenon of forming membrane-less condensates through partitioning of biomolecules. Such a process is also critically involved in the establishment of high-order chromatin structures and looping. IDR- and phase separation-driven 3D genome (re)organization often goes wrong in disease such as cancer. This review discusses about recent advances in understanding how phase separation of intrinsically disordered proteins (IDPs) modulates chromatin looping and gene expression.

Chemical reaction motifs driving non-equilibrium behaviours in phase separating materials

Chemical reactions that couple to systems that phase separate have been implicated in diverse contexts from biology to materials science. However, how a particular set of chemical reactions (chemical reaction network, CRN) would affect the behaviours of a phase separating system is difficult to fully predict theoretically. In this paper, we analyse a mean field theory coupling CRNs to a combined system of phase separating and non-phase separating materials and analyse how the properties of the CRNs affect different classes of non-equilibrium behaviour: microphase separation or temporally oscillating patterns. We examine the problem of achieving microphase separated condensates by statistical analysis of the Jacobians, of which the most important motifs are negative feedback of the phase separating component and combined inhibition/activation by the non-phase separating components. We then identify CRN motifs that are likely to yield microphase by examining randomly generated networks and parameters. Molecular sequestration of the phase separating motif is shown to be the most robust towards yielding microphase separation. Subsequently, we find that dynamics of the phase separating species is promoted most easily by inducing oscillations in the diffusive components coupled to the phase separating species. Our results provide guidance towards the design of CRNs that manage the formation, dissolution and organization of compartments.

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.

Developing efficient deep learning model for predicting copolymer properties

Deep learning models are gaining popularity and potency in predicting polymer properties. These models can be built using pre-existing data and are useful for the rapid prediction of polymer properties. However, the performance of a deep learning model is intricately connected to its topology and the volume of training data. There is no facile protocol available to select a deep learning architecture, and there is a lack of a large volume of homogeneous sequence-property data of polymers. These two factors are the primary bottleneck for the efficient development of deep learning models for polymers. Here we assess the severity of these factors and propose strategies to address them. We show that a linear layer-by-layer expansion of a neural network can help in identifying the best neural network topology for a given problem. Moreover, we map the discrete sequence space of a polymer to a continuous one-dimensional latent space using a feature extraction technique to identify minimal data points for training a deep learning model. We implement these approaches for two representative cases of building sequence-property surrogate models, viz., the single-molecule radius of gyration of a copolymer and copolymer compatibilizer. This work demonstrates efficient methods for building deep learning models with minimal data and hyperparameters for predicting sequence-defined properties of polymers.

Thermodynamic forces from protein and water govern condensate formation of an intrinsically disordered protein domain

Liquid-liquid phase separation (LLPS) can drive a multitude of cellular processes by compartmentalizing biological cells via the formation of dense liquid biomolecular condensates, which can function as membraneless organelles. Despite its importance, the molecular-level understanding of the underlying thermodynamics of this process remains incomplete. In this study, we use atomistic molecular dynamics simulations of the low complexity domain (LCD) of human fused in sarcoma (FUS) protein to investigate the contributions of water and protein molecules to the free energy changes that govern LLPS. Both protein and water components are found to have comparably sizeable thermodynamic contributions to the formation of FUS condensates. Moreover, we quantify the counteracting effects of water molecules that are released into the bulk upon condensate formation and the waters retained within the protein droplets. Among the various factors considered, solvation entropy and protein interaction enthalpy are identified as the most important contributions, while solvation enthalpy and protein entropy changes are smaller. These results provide detailed molecular insights on the intricate thermodynamic interplay between protein- and solvation-related forces underlying the formation of biomolecular condensates.

Theoretical and Data-Driven Approaches for Biomolecular Condensates

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

Transfer Learning Facilitates the Prediction of Polymer-Surface Adhesion Strength

Machine learning (ML) accelerates the exploration of material properties and their links to the structure of the underlying molecules. In previous work [Shi et al. ACS Applied Materials & Interfaces 2022, 14, 37161-37169.], ML models were applied to predict the adhesive free energy of polymer-surface interactions with high accuracy from the knowledge of the sequence data, demonstrating successes in inverse-design of polymer sequence for known surface compositions. While the method was shown to be successful in designing polymers for a known surface, extensive data sets were needed for each specific surface in order to train the surrogate models. Ideally, one should be able to infer information about similar surfaces without having to regenerate a full complement of adhesion data for each new case. In the current work, we demonstrate a transfer learning (TL) technique using a deep neural network to improve the accuracy of ML models trained on small data sets by pretraining on a larger database from a related system and fine-tuning the weights of all layers with a small amount of additional data. The shared knowledge from the pretrained model facilitates the prediction accuracy significantly on small data sets. We also explore the limits of database size on accuracy and the optimal tuning of network architecture and parameters for our learning tasks. While applied to a relatively simple coarse-grained (CG) polymer model, the general lessons of this study apply to detailed modeling studies and the broader problems of inverse materials design.

Surfactants or scaffolds? RNAs of varying lengths control the thermodynamic stability of condensates differently

Biomolecular condensates, thought to form via liquid-liquid phase separation of intracellular mixtures, are multicomponent systems that can include diverse types of proteins and RNAs. RNA is a critical modulator of RNA-protein condensate stability, as it induces an RNA concentration-dependent reentrant phase transition-increasing stability at low RNA concentrations and decreasing it at high concentrations. Beyond concentration, RNAs inside condensates can be heterogeneous in length, sequence, and structure. Here, we use multiscale simulations to understand how different RNA parameters interact with one another to modulate the properties of RNA-protein condensates. To do so, we perform residue/nucleotide resolution coarse-grained molecular dynamics simulations of multicomponent RNA-protein condensates containing RNAs of different lengths and concentrations, and either FUS or PR25 proteins. Our simulations reveal that RNA length regulates the reentrant phase behavior of RNA-protein condensates: increasing RNA length sensitively rises the maximum value that the critical temperature of the mixture reaches, and the maximum concentration of RNA that the condensate can incorporate before beginning to become unstable. Strikingly, RNAs of different lengths are organized heterogeneously inside condensates, which allows them to enhance condensate stability via two distinct mechanisms: shorter RNA chains accumulate at the condensate's surface acting as natural biomolecular surfactants, while longer RNA chains concentrate inside the core to saturate their bonds and enhance the density of molecular connections in the condensate. Using a patchy particle model, we additionally demonstrate that the combined impact of RNA length and concentration on condensate properties is dictated by the valency, binding affinity, and polymer length of the various biomolecules involved. Our results postulate that diversity on RNA parameters within condensates allows RNAs to increase condensate stability by fulfilling two different criteria: maximizing enthalpic gain and minimizing interfacial free energy; hence, RNA diversity should be considered when assessing the impact of RNA on biomolecular condensates regulation.

Phase separation and molecular ordering of the prion-like domain of the Arabidopsis thermosensory protein EARLY FLOWERING 3

Liquid-liquid phase separation (LLPS) is an important mechanism enabling the dynamic compartmentalization of macromolecules, including complex polymers such as proteins and nucleic acids, and occurs as a function of the physicochemical environment. In the model plant, Arabidopsis thaliana, LLPS by the protein EARLY FLOWERING3 (ELF3) occurs in a temperature-sensitive manner and controls thermoresponsive growth. ELF3 contains a largely unstructured prion-like domain (PrLD) that acts as a driver of LLPS in vivo and in vitro. The PrLD contains a poly-glutamine (polyQ) tract, whose length varies across natural Arabidopsis accessions. Here, we use a combination of biochemical, biophysical, and structural techniques to investigate the dilute and condensed phases of the ELF3 PrLD with varying polyQ lengths. We demonstrate that the dilute phase of the ELF3 PrLD forms a monodisperse higher-order oligomer that does not depend on the presence of the polyQ sequence. This species undergoes LLPS in a pH- and temperature-sensitive manner and the polyQ region of the protein tunes the initial stages of phase separation. The liquid phase rapidly undergoes aging and forms a hydrogel as shown by fluorescence and atomic force microscopies. Furthermore, we demonstrate that the hydrogel assumes a semiordered structure as determined by small-angle X-ray scattering, electron microscopy, and X-ray diffraction. These experiments demonstrate a rich structural landscape for a PrLD protein and provide a framework to describe the structural and biophysical properties of biomolecular condensates.

Theory and Simulation of Multiphase Coexistence in Biomolecular Mixtures

Biomolecular condensates constitute a newly recognized form of spatial organization in living cells. Although many condensates are believed to form as a result of phase separation, the physicochemical properties that determine the phase behavior of heterogeneous biomolecular mixtures are only beginning to be explored. Theory and simulation provide invaluable tools for probing the relationship between molecular determinants, such as protein and RNA sequences, and the emergence of phase-separated condensates in such complex environments. This review covers recent advances in the prediction and computational design of biomolecular mixtures that phase-separate into many coexisting phases. First, we review efforts to understand the phase behavior of mixtures with hundreds or thousands of species using theoretical models and statistical approaches. We then describe progress in developing analytical theories and coarse-grained simulation models to predict multiphase condensates with the molecular detail required to make contact with biophysical experiments. We conclude by summarizing the challenges ahead for modeling the inhomogeneous spatial organization of biomolecular mixtures in living cells.

Role of Strong Localized vs Weak Distributed Interactions in Disordered Protein Phase Separation

Interaction strength and localization are critical parameters controlling the single-chain and condensed-state properties of intrinsically disordered proteins (IDPs). Here, we decipher these relationships using coarse-grained heteropolymers comprised of hydrophobic (H) and polar (P) monomers as model IDPs. We systematically vary the fraction of P monomers XP and employ two distinct particle-based models that include either strong localized attractions between only H-H pairs (HP model) or weak distributed attractions between both H-H and H-P pairs (HP+ model). To compare different sequences and models, we first carefully tune the attraction strength for all sequences to match the single-chain radius of gyration. Interestingly, we find that this procedure produces similar conformational ensembles, nonbonded potential energies, and chain-level dynamics for single chains of almost all sequences in both models, with some deviations for the HP model at large XP. However, we observe a surprisingly rich phase behavior for the sequences in both models that deviates from the expectation that similarity at the single-chain level will translate to a similar phase-separation propensity. Coexistence between dilute and dense phases is only observed up to a model-dependent XP, despite the presence of favorable interchain interactions, which we quantify using the second virial coefficient. Instead, the limited number of attractive sites (H monomers) leads to the self-assembly of finite-sized clusters of different sizes depending on XP. Our findings strongly suggest that models with distributed interactions favor the formation of liquid-like condensates over a much larger range of sequence compositions compared to models with localized interactions.

Phase separation and aggregation in multiblock chains

This paper focuses on phase and aggregation behavior for linear chains composed of blocks of hydrophilic and hydrophobic segments. Phase and conformational transitions of patterned chains are relevant for understanding liquid-liquid separation of biomolecular condensates, which play a prominent role in cellular biophysics and for surfactant and polymer applications. Previous studies of simple models for multiblock chains have shown that, depending on the sequence pattern and chain length, such systems can fall into one of two categories: displaying either phase separation or aggregation into finite-size clusters. The key new result of this paper is that both formation of finite-size aggregates and phase separation can be observed for certain chain architectures at appropriate conditions of temperature and concentration. For such systems, a bulk dense liquid condenses from a dilute phase that already contains multi-chain finite-size aggregates. The computational approach used in this study involves several distinct steps using histogram-reweighting grand canonical Monte Carlo simulations, which are described in some level of detail.

Copolymers with Nonblocky Sequences as Novel Materials with Finely Tuned Properties

The copolymer sequence can be considered as a new tool to shape the resulting system properties on demand. This perspective is devoted to copolymers with "partially segregated" (or nonblocky) sequences. Such copolymers include gradient copolymers and copolymers with random sequences as well as copolymers with precisely controlled sequences. We overview recent developments in the synthesis of these systems as well as new findings regarding their properties, in particular, self-assembly in solutions and in melts. An emphasis is put on how the microscopic behavior of polymer chains is influenced by the chain sequences. In addition to that, a novel class of approaches allowing one to efficiently tackle the problem of copolymer chain sequence design─data driven methods (artificial intelligence and machine learning)─is discussed.

Phase diagrams-Why they matter and how to predict them

Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.

A New Phase of Networking: The Molecular Composition and Regulatory Dynamics of Mammalian Stress Granules

Stress granules (SGs) are cytosolic biomolecular condensates that form in response to cellular stress. Weak, multivalent interactions between their protein and RNA constituents drive their rapid, dynamic assembly through phase separation coupled to percolation. Though a consensus model of SG function has yet to be determined, their perceived implication in cytoprotective processes (e.g., antiviral responses and inhibition of apoptosis) and possible role in the pathogenesis of various neurodegenerative diseases (e.g., amyotrophic lateral sclerosis and frontotemporal dementia) have drawn great interest. Consequently, new studies using numerous cell biological, genetic, and proteomic methods have been performed to unravel the mechanisms underlying SG formation, organization, and function and, with them, a more clearly defined SG proteome. Here, we provide a consensus SG proteome through literature curation and an update of the user-friendly database RNAgranuleDB to version 2.0 (http://rnagranuledb.lunenfeld.ca/). With this updated SG proteome, we use next-generation phase separation prediction tools to assess the predisposition of SG proteins for phase separation and aggregation. Next, we analyze the primary sequence features of intrinsically disordered regions (IDRs) within SG-resident proteins. Finally, we review the protein- and RNA-level determinants, including post-translational modifications (PTMs), that regulate SG composition and assembly/disassembly dynamics.

Polymer sequence design via molecular simulation-based active learning

Molecular-scale interactions and chemical structures offer an enormous opportunity to tune material properties. However, designing materials from their molecular scale is a grand challenge owing to the practical limitations in exploring astronomically large design spaces using traditional experimental or computational methods. Advancements in data science and machine learning have produced a host of tools and techniques that can address this problem and facilitate the efficient exploration of large search spaces. In this work, a blended approach integrating physics-based methods, machine learning techniques and uncertainty quantification is implemented to effectively screen a macromolecular sequence space and design target structures. Here, we survey and assess the efficacy of data-driven methods within the framework of active learning for a challenging design problem, viz., sequence optimization of a copolymer. We report the impact of surrogate models, kernels, and initial conditions on the convergence of the active learning method for the sequence design problem. This work establishes optimal strategies and hyperparameters for efficient inverse design of polymer sequences via active learning.

Field-Theoretic Simulation Method to Study the Liquid-Liquid Phase Separation of Polymers

Liquid-liquid phase separation (LLPS) is a process that results in the formation of a polymer-rich liquid phase coexisting with a polymer-depleted liquid phase. LLPS plays a critical role in the cell through the formation of membrane-less organelles, but it also has a number of biotechnical and biomedical applications such as drug confinement and its targeted delivery. In this chapter, we present a computational efficient methodology that uses field-theoretic simulations (FTS) with complex Langevin (CL) sampling to characterize polymer phase behavior and delineate the LLPS phase boundaries. This approach is a powerful complement to analytical and explicit-particle simulations, and it can serve to inform experimental LLPS studies. The strength of the method lies in its ability to properly sample a large ensemble of polymers in a saturated solution while including the effect of composition fluctuations on LLPS. We describe the approaches that can be used to accurately construct phase diagrams of a variety of molecularly designed polymers and illustrate the method by generating an approximation-free phase diagram for a classical symmetric diblock polyampholyte.

Effects of Conformational Constraint on Peptide Solubility Limits

Liquid-liquid phase separation of proteins preferentially involves intrinsically disordered proteins or disordered regions. Understanding the solution chemistry of these phase separations is key to learning how to quantify and manipulate systems that involve such processes. Here, we investigate the effect of cyclization on the liquid-liquid phase separation of short polyglycine peptides. We simulated separate aqueous systems of supersaturated cyclic and linear GGGGG and observed spontaneous liquid-liquid phase separation in each of the solutions. The cyclic GGGGG phase separates less robustly than linear GGGGG and has a higher aqueous solubility, even though linear GGGGG has a more favorable single molecule solvation free energy. The versatile and abundant interpeptide contacts formed by the linear GGGGG stabilize the condensed droplet phase, driving the phase separation in this system. In particular, we find that van der Waals close contact interactions are enriched in the droplet phase as opposed to electrostatic interactions. An analysis of the change in backbone conformational entropy that accompanies the phase transition revealed that cyclic peptides lose significantly less entropy in this process as expected. However, we find that the enhanced interaction enthalpy of linear GGGGG in the droplet phase is enough to compensate for a larger decrease in conformational entropy.

Computational Study of Mechanochemical Activation in Nanostructured Triblock Copolymers

Force-driven chemical reactions have emerged as an attractive platform for diverse applications in polymeric materials. However, the microscopic chain conformations and topologies necessary for efficiently transducing macroscopic forces to the molecular scale are not well-understood. In this work, we use a coarse-grained model to investigate the impact of network-like topologies on mechanochemical activation in self-assembled triblock copolymers. We find that mechanochemical activation during tensile deformation depends strongly on both the polymer composition and chain conformation in these materials. Activation primarily occurs in the tie chains connecting different glassy domains and in loop chains that are hooked onto each other by physical entanglements. Activation also requires a higher stress in materials having a higher glassy block content. Overall, the lamellar samples show the highest percent activation at high stress. In contrast, at low stress, the spherical morphology, which has the lowest glassy fraction, shows the highest activation. Additionally, we observe a spatial pattern of activation, which appears to be tied to distortion of the self-assembled morphology. Higher activation is observed in the tips of the chevrons formed during deformation of lamellar samples as well as in the centers between the cylinders in the cylindrical morphology. Our work shows that changes in the network-like topology in different morphologies significantly impact mechanochemical activation efficiencies in these materials, suggesting that this area will be a fruitful avenue for further experimental research.

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.

SpiDec: Computing binodals and interfacial tension of biomolecular condensates from simulations of spinodal decomposition

Phase separation of intrinsically disordered proteins (IDPs) is a phenomenon associated with many essential cellular processes, but a robust method to compute the binodal from molecular dynamics simulations of IDPs modeled at the all-atom level in explicit solvent is still elusive, due to the difficulty in preparing a suitable initial dense configuration and in achieving phase equilibration. Here we present SpiDec as such a method, based on spontaneous phase separation via spinodal decomposition that produces a dense slab when the system is initiated at a homogeneous, low density. After illustrating the method on four model systems, we apply SpiDec to a tetrapeptide modeled at the all-atom level and solvated in TIP3P water. The concentrations in the dense and dilute phases agree qualitatively with experimental results and point to binodals as a sensitive property for force-field parameterization. SpiDec may prove useful for the accurate determination of the phase equilibrium of IDPs.

Protein structural transitions critically transform the network connectivity and viscoelasticity of RNA-binding protein condensates but RNA can prevent it

Biomolecular condensates, some of which are liquid-like during health, can age over time becoming gel-like pathological systems. One potential source of loss of liquid-like properties during ageing of RNA-binding protein condensates is the progressive formation of inter-protein β-sheets. To bridge microscopic understanding between accumulation of inter-protein β-sheets over time and the modulation of FUS and hnRNPA1 condensate viscoelasticity, we develop a multiscale simulation approach. Our method integrates atomistic simulations with sequence-dependent coarse-grained modelling of condensates that exhibit accumulation of inter-protein β-sheets over time. We reveal that inter-protein β-sheets notably increase condensate viscosity but does not transform the phase diagrams. Strikingly, the network of molecular connections within condensates is drastically altered, culminating in gelation when the network of strong β-sheets fully percolates. However, high concentrations of RNA decelerate the emergence of inter-protein β-sheets. Our study uncovers molecular and kinetic factors explaining how the accumulation of inter-protein β-sheets can trigger liquid-to-solid transitions in condensates, and suggests a potential mechanism to slow such transitions down.

In this work the authors propose a multiscale computational approach, integrating atomistic and coarse-grained models simulations, to study the thermodynamic and kinetic factors playing a major role in the liquid-to-solid transition of biomolecular condensates. It is revealed how the gradual accumulation of inter-protein β-sheets increases the viscosity of functional liquid-like condensates, transforming them into gel-like pathological aggregates, and it is also shown how high concentrations of RNA can decelerate such transition.

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

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

Predicting Adhesive Free Energies of Polymer-Surface Interactions with Machine Learning

Polymer-surface interactions are crucial to many biological processes and industrial applications. Here we propose a machine learning method to connect a model polymer's sequence with its adhesion to decorated surfaces. We simulate the adhesive free energies of 20000 unique coarse-grained one-dimensional polymer sequences interacting with functionalized surfaces and build support vector regression models that demonstrate inexpensive and reliable prediction of the adhesive free energy as a function of sequence. Our work highlights the promising integration of coarse-grained simulation with data-driven machine learning methods for the design of functional polymers and represents an important step toward linking polymer compositions with polymer-surface interactions.

Thermodynamic stability and critical points in multicomponent mixtures with structured interactions

Theoretical work has shed light on the phase behavior of idealized mixtures of many components with random interactions. However, typical mixtures interact through particular physical features, leading to a structured, nonrandom interaction matrix of lower rank. Here, we develop a theoretical framework for such mixtures and derive mean-field conditions for thermodynamic stability and critical behavior. Irrespective of the number of components and features, this framework allows for a generally lower-dimensional representation in the space of features and proposes a principled way to coarse-grain multicomponent mixtures as binary mixtures. Moreover, it suggests a way to systematically characterize different series of critical points and their codimensions in mean-field. Since every pairwise interaction matrix can be expressed in terms of features, our work is applicable to a broad class of mean-field models.

Designing Enhanced Entropy Binding in Single-Chain Nanoparticles

Single-chain nanoparticles (SCNPs) are a new class of bio- and soft-matter polymeric objects in which a fraction of the monomers are able to form equivalently intra- or interpolymer bonds. Here we numerically show that a fully entropic gas-liquid phase separation can take place in SCNP systems. Control over the discontinuous (first-order) change-from a phase of independent diluted (fully-bonded) polymers to a phase in which polymers entropically bind to each other to form a (fully-bonded) polymer network-can be achieved by a judicious design of the patterns of reactive monomers along the polymer chain. Such a sensitivity arises from a delicate balance between the distinct entropic contributions controlling the binding.

Predicting aggregate morphology of sequence-defined macromolecules with recurrent neural networks

Self-assembly of dilute sequence-defined macromolecules is a complex phenomenon in which the local arrangement of chemical moieties can lead to the formation of long-range structure. The dependence of this structure on the sequence necessarily implies that a mapping between the two exists, yet it has been difficult to model so far. Predicting the aggregation behavior of these macromolecules is challenging due to the lack of effective order parameters, a vast design space, inherent variability, and high computational costs associated with currently available simulation techniques. Here, we accurately predict the morphology of aggregates self-assembled from sequence-defined macromolecules using supervised machine learning. We find that regression models with implicit representation learning perform significantly better than those based on engineered features such as k-mer counting, and a recurrent-neural-network-based regressor performs the best out of nine model architectures we tested. Furthermore, we demonstrate the high-throughput screening of monomer sequences using the regression model to identify candidates for self-assembly into selected morphologies. Our strategy is shown to successfully identify multiple suitable sequences in every test we performed, so we hope the insights gained here can be extended to other increasingly complex design scenarios in the future, such as the design of sequences under polydispersity and at varying environmental conditions.

Liquid-Liquid Phase Separation and Assembly of Silk-like Proteins is Dependent on the Polymer Length

Phase transitions have an essential role in the assembly of nature’s protein-based materials into hierarchically organized structures, yet many of the underlying mechanisms and interactions remain to be resolved. A central question for designing proteins for materials is how the protein architecture and sequence affects the nature of the phase transitions and resulting assembly. In this work, we produced 82 kDa (1×), 143 kDa (2×), and 204 kDa (3×) silk-mimicking proteins by taking advantage of protein ligation by SpyCatcher/Tag protein-peptide pair. We show that the three silk proteins all undergo a phase transition from homogeneous solution to assembly formation. In the assembly phase, a length- and concentration-dependent transition between two distinct assembly morphologies, one forming aggregates and another coacervates, exists. The coacervates showed properties that were dependent on the protein size. Computational modeling of the proteins by a bead-spring model supports the experimental results and provides us a possible mechanistic origin for the assembly transitions based on architectures and interactions.

Affinity of disordered protein complexes is modulated by entropy-energy reinforcement

Intrinsically disordered proteins (IDPs), which are very common and essential to many biological activities, sometimes function via interaction with another IDP and form a fuzzy complex, which can be highly stable. It is unclear what the biophysical forces are that govern their thermodynamics and specificity, which are essential for de novo fuzzy complex design. Here, we explored the fuzzy complex formed between ProTα and H1, which are oppositely charged IDPs, by swapping the charges between them, generating variants that have either greater polyampholytic or polyelectrolytic nature as well as different charge patterns. Charge swapping and shuffling dramatically change the affinity of the fuzzy complex, which is contributed to by both enthalpy and entropy, where the latter is dominated by counterion release.

The association between two intrinsically disordered proteins (IDPs) may produce a fuzzy complex characterized by a high binding affinity, similar to that found in the ultrastable complexes formed between two well-structured proteins. Here, using coarse-grained simulations, we quantified the biophysical forces driving the formation of such fuzzy complexes. We found that the high-affinity complex formed between the highly and oppositely charged H1 and ProTα proteins is sensitive to electrostatic interactions. We investigated 52 variants of the complex by swapping charges between the two oppositely charged proteins to produce sequences whose negatively or positively charged residue content was more homogeneous or heterogenous (i.e., polyelectrolytic or polyampholytic, having higher or lower absolute net charges, respectively) than the wild type. We also changed the distributions of oppositely charged residues within each participating sequence to produce variants in which the charges were segregated or well mixed. Both types of changes significantly affect binding affinity in fuzzy complexes, which is governed by both enthalpy and entropy. The formation of H1–ProTa is supported by an increase in configurational entropy and by entropy due to counterion release. The latter can be twice as large as the former, illustrating the dominance of counterion entropy in modulating the binding thermodynamics. Complexes formed between proteins with greater absolute net charges are more stable, both enthalpically and entropically, indicating that enthalpy and entropy have a mutually reinforcing effect. The sensitivity of the thermodynamics of the complex to net charge and the charge pattern within each of the binding constituents may provide a means to achieve binding specificity between IDPs.

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.

Tuning Formation of Protein-DNA Coacervates by Sequence and Environment

The high concentration of nucleic acids and positively charged proteins in the cell nucleus provides many possibilities for complex coacervation. We consider a prototypical mixture of nucleic acids together with the polycationic C-terminus of histone H1 (CH1). Using a minimal coarse-grained model that captures the shape, flexibility, and charge distributions of the molecules, we find that coacervates are readily formed at physiological ionic strengths, in agreement with experiment, with a progressive increase in local ordering at low ionic strength. Variation of the positions of charged residues in the protein tunes phase separation: for well-mixed or only moderately blocky distributions of charge, there is a modest increase of local ordering with increasing blockiness that is also associated with an increased propensity to phase separate. This ordering is also associated with a slowdown of rotational and translational diffusion in the dense phase. However, for more extreme blockiness (and consequently higher local charge density), we see a qualitative change in the condensed phase to become a segregated structure with a dramatically increased ordering of the DNA. Naturally occurring proteins with these sequence properties, such as protamines in sperm cells, are found to be associated with very dense packing of nucleic acids.

Kinetic interplay between droplet maturation and coalescence modulates shape of aged protein condensates

Biomolecular condensates formed by the process of liquid-liquid phase separation (LLPS) play diverse roles inside cells, from spatiotemporal compartmentalisation to speeding up chemical reactions. Upon maturation, the liquid-like properties of condensates, which underpin their functions, are gradually lost, eventually giving rise to solid-like states with potential pathological implications. Enhancement of inter-protein interactions is one of the main mechanisms suggested to trigger the formation of solid-like condensates. To gain a molecular-level understanding of how the accumulation of stronger interactions among proteins inside condensates affect the kinetic and thermodynamic properties of biomolecular condensates, and their shapes over time, we develop a tailored coarse-grained model of proteins that transition from establishing weak to stronger inter-protein interactions inside condensates. Our simulations reveal that the fast accumulation of strongly binding proteins during the nucleation and growth stages of condensate formation results in aspherical solid-like condensates. In contrast, when strong inter-protein interactions appear only after the equilibrium condensate has been formed, or when they accumulate slowly over time with respect to the time needed for droplets to fuse and grow, spherical solid-like droplets emerge. By conducting atomistic potential-of-mean-force simulations of NUP-98 peptides-prone to forming inter-protein [Formula: see text]-sheets-we observe that formation of inter-peptide [Formula: see text]-sheets increases the strength of the interactions consistently with the loss of liquid-like condensate properties we observe at the coarse-grained level. Overall, our work aids in elucidating fundamental molecular, kinetic, and thermodynamic mechanisms linking the rate of change in protein interaction strength to condensate shape and maturation during ageing.

RNA length has a non-trivial effect in the stability of biomolecular condensates formed by RNA-binding proteins

Biomolecular condensates formed via liquid-liquid phase separation (LLPS) play a crucial role in the spatiotemporal organization of the cell material. Nucleic acids can act as critical modulators in the stability of these protein condensates. To unveil the role of RNA length in regulating the stability of RNA binding protein (RBP) condensates, we present a multiscale computational strategy that exploits the advantages of a sequence-dependent coarse-grained representation of proteins and a minimal coarse-grained model wherein proteins are described as patchy colloids. We find that for a constant nucleotide/protein ratio, the protein fused in sarcoma (FUS), which can phase separate on its own-i.e., via homotypic interactions-only exhibits a mild dependency on the RNA strand length. In contrast, the 25-repeat proline-arginine peptide (PR25), which does not undergo LLPS on its own at physiological conditions but instead exhibits complex coacervation with RNA-i.e., via heterotypic interactions-shows a strong dependence on the length of the RNA strands. Our minimal patchy particle simulations suggest that the strikingly different effect of RNA length on homotypic LLPS versus RBP-RNA complex coacervation is general. Phase separation is RNA-length dependent whenever the relative contribution of heterotypic interactions sustaining LLPS is comparable or higher than those stemming from protein homotypic interactions. Taken together, our results contribute to illuminate the intricate physicochemical mechanisms that influence the stability of RBP condensates through RNA inclusion.

Surface tension and super-stoichiometric surface enrichment in two-component biomolecular condensates

Cells can achieve internal organization by exploiting liquid-liquid phase separation to form biomolecular condensates. Here we focus on the surface properties of condensates composed of two multivalent associative polymers held together by one-to-one “sticker” bonds. Using coarse-grained molecular-dynamics simulations, we study the influence of component stoichiometry on condensate surface properties. We find that unequal stoichiometry results in enrichment of the majority species at the interface and a sharp reduction of surface tension. To relate these two effects, we show that the reduction in surface tension scales linearly with the excess concentration of free binding sites at the interface. Our results imply that each excess free site contributes an approximately fixed additional energy and entropy to the interface, with the latter dominating such that enrichment of free majority sites lowers the surface tension. Our work provides insight into novel physical mechanisms by which cells can regulate condensate surface properties.

•

Stoichiometry controls the surface tension of two-component biomolecular condensates

•

Unequal stoichiometry leads to enrichment of the majority species at the interface

•

Enrichment of the free majority binding sites increases the interfacial entropy

•

Surface tension is drastically reduced at unequal stoichiometry

Limitations of field-theory simulation for exploring phase separation: The role of repulsion in a lattice protein model

Field-theory simulation by the complex Langevin method offers an alternative to conventional sampling techniques for exploring the forces driving biomolecular liquid-liquid phase separation. Such simulations have recently been used to study several polyampholyte systems. Here, we formulate a field theory corresponding to the hydrophobic/polar (HP) lattice protein model, with finite same-site repulsion and nearest-neighbor attraction between HH bead pairs. By direct comparison with particle-based Monte Carlo simulations, we show that complex Langevin sampling of the field theory reproduces the thermodynamic properties of the HP model only if the same-site repulsion is not too strong. Unfortunately, the repulsion has to be taken weaker than what is needed to prevent condensed droplets from assuming an artificially compact shape. Analysis of a minimal and analytically solvable toy model hints that the sampling problems caused by repulsive interaction may stem from loss of ergodicity.