Fundamental Aspects of Phase-Separated Biomolecular Condensates

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

GENESIS CGDYN: large-scale coarse-grained MD simulation with dynamic load balancing for heterogeneous biomolecular systems

Residue-level coarse-grained (CG) molecular dynamics (MD) simulation is widely used to investigate slow biological processes that involve multiple proteins, nucleic acids, and their complexes. Biomolecules in a large simulation system are distributed non-uniformly, limiting computational efficiency with conventional methods. Here, we develop a hierarchical domain decomposition scheme with dynamic load balancing for heterogeneous biomolecular systems to keep computational efficiency even after drastic changes in particle distribution. These schemes are applied to the dynamics of intrinsically disordered protein (IDP) droplets. During the fusion of two droplets, we find that the changes in droplet shape correlate with the mixing of IDP chains. Additionally, we simulate large systems with multiple IDP droplets, achieving simulation sizes comparable to those observed in microscopy. In our MD simulations, we directly observe Ostwald ripening, a phenomenon where small droplets dissolve and their molecules redeposit into larger droplets. These methods have been implemented in CGDYN of the GENESIS software, offering a tool for investigating mesoscopic biological processes using the residue-level CG models.

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.

Sequence determinants and solution conditions underlying liquid to solid phase transition

Life consists of numberless functional biomolecules that exist in various states. Besides well-dissolved phases, biomolecules especially proteins and nucleic acids can form liquid droplets through liquid-liquid phase separation (LLPS). Stronger interactions promote a solid-like state of biomolecular condensates, which are also formerly referred to as detergent-insoluble aggregates. Solid-like condensates exist in vivo physiologically and pathologically, and their formation has not been fully understood. Recently, more and more research has proven that liquid to solid phase transition (LST) is an essential way to form solid condensates. In this review, we summarized the regions in the sequence that have different impacts on phase transition and emphasized that the LST is affected by its sequence characteristics. Moreover, increasing evidence unveiled that LST is affected by various solution conditions. We discussed solution conditions like protein concentration, pH, ATP, ions, and small molecules in a solution. Methods have been established to study these solid phase components. Here, we summarized low-throughput experimental techniques and high-throughput omics methods in the study of the LST.

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.