Macromolecular Regulation of the Material Properties of Biomolecular Condensates

Amino Acid-Dependent Material Properties of Tetrapeptide Condensates

Condensates formed by intrinsically disordered proteins mediate a myriad of cellular processes and are linked to pathological conditions including neurodegeneration. Rules of how different types of amino acids (e.g., π-π pairs) dictate the physical properties of biomolecular condensates are emerging, but our understanding of the roles of different amino acids is far from complete. Here we studied condensates formed by tetrapeptides of the form XXssXX, where X is an amino acid and ss represents a disulfide bond along the backbone. Eight peptides form four types of condensates at different concentrations and pH values: droplets (X = F, L, M, P, V, A); amorphous dense liquids (X = L, M, P, V, A); amorphous aggregates (X = W), and gels (X = I, V, A). The peptides exhibit enormous differences in phase equilibrium and material properties, including a 368-fold range in the threshold concentration for phase separation and a 3856-fold range in viscosity. All-atom molecular dynamics simulations provide physical explanations of these results. The present work also reveals widespread critical behaviors, including critical slowing down manifested by the formation of amorphous dense liquids and critical scaling obeyed by fusion speed, with broad implications for condensate function.

Sequence-dependent material properties of biomolecular condensates and their relation to dilute phase conformations

Material properties of phase-separated biomolecular condensates, enriched with disordered proteins, dictate many cellular functions. Contrary to the progress made in understanding the sequence-dependent phase separation of proteins, little is known about the sequence determinants of condensate material properties. Using the hydropathy scale and Martini models, we computationally decipher these relationships for charge-rich disordered protein condensates. Our computations yield dynamical, rheological, and interfacial properties of condensates that are quantitatively comparable with experimentally characterized condensates. Interestingly, we find that the material properties of model and natural proteins respond similarly to charge segregation, despite different sequence compositions. Molecular interactions within the condensates closely resemble those within the single-chain ensembles. Consequently, the material properties strongly correlate with molecular contact dynamics and single-chain structural properties. We demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Determinants of viscoelasticity and flow activation energy in biomolecular condensates

The form and function of biomolecular condensates are intimately linked to their material properties. Here, we integrate microrheology with molecular simulations to dissect the physical determinants of condensate fluid phase dynamics. By quantifying the timescales and energetics of network relaxation in a series of heterotypic viscoelastic condensates, we uncover distinctive roles of sticker motifs, binding energy, and chain length in dictating condensate dynamical properties. We find that the mechanical relaxation times of condensate-spanning networks are determined by both intermolecular interactions and chain length. We demonstrate, however, that the energy barrier for network reconfiguration, termed flow activation energy, is independent of chain length and only varies with the strengths of intermolecular interactions. Biomolecular diffusion in the dense phase depends on a complex interplay between viscoelasticity and flow activation energy. Our results illuminate distinctive roles of chain length and sequence-specific multivalent interactions underlying the complex material and transport properties of biomolecular condensates.

Adenosine Triphosphate Mediates Phase Separation of Disordered Basic Proteins by Bridging Intermolecular Interaction Networks

Adenosine triphosphate (ATP) is an abundant molecule with crucial cellular roles as the energy currency and a building block of nucleic acids and for protein phosphorylation. Here we show that ATP mediates the phase separation of basic intrinsically disordered proteins (bIDPs). In the resulting condensates, ATP is highly concentrated (apparent partition coefficients up to 7700) and serves as bridges between bIDP chains. These liquid-like droplets have some of the lowest interfacial tension (∼25 pN/μm) but high zero-shear viscosities (1-15 Pa s) due to the bridged protein networks, and yet their fusion has some of the highest speeds (∼1 μm/ms). The rapid fusion manifests extreme shear thinning, where the apparent viscosity is lower than zero-shear viscosity by over 100-fold, made possible by fast reformation of the ATP bridges. At still higher concentrations, ATP does not dissolve bIDP droplets but results in aggregates and fibrils.

ATP Mediates Phase Separation of Disordered Basic Proteins by Bridging Intermolecular Interaction Networks

ATP is an abundant molecule with crucial cellular roles as the energy currency and a building block of nucleic acids and for protein phosphorylation. Here we show that ATP mediates the phase separation of basic intrinsically disordered proteins (bIDPs). In the resulting condensates, ATP is highly concentrated (apparent partition coefficients at 200-5000) and serves as bridges between bIDP chains. These liquid-like droplets have some of the lowest interfacial tension (~25 pN/μm) but high zero-shear viscosities (1-15 Pa s) due to the bridged protein networks, and yet their fusion has some of the highest speeds (~1 μm/ms). The rapid fusion manifests extreme shear thinning, where the apparent viscosity is lower than zero-shear viscosity by over 100-fold, made possible by fast reformation of the ATP bridges. At still higher concentrations, ATP does not dissolve bIDP droplets but results in aggregates and fibrils.

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.

Calculating Binodals and Interfacial Tension of Phase-Separated Condensates from Molecular Simulations with Finite-Size Corrections

We illustrate three methods for calculating the binodals of phase-separated condensates from molecular simulations. Because molecular simulations can only be carried out for small system sizes, correction for finite sizes may be required for making direct comparison between calculated results and experimental data. We first summarize the three methods and then present detailed implementation of each method on a Lennard-Jones fluid. In the first method, chemical potentials are calculated over a range of particle densities in canonical-ensemble simulations; the densities of the dilute and dense phases at the given temperature are then found by a Maxwell equal-area construction. In Gibbs-ensemble Monte Carlo, the exchange between separated dilute and dense phases is simulated to obtain their densities. Lastly, slab-geometry molecular dynamics simulations model the dilute and dense phases in coexistence and yield not only their densities but also their interfacial tension. The three types of simulations are carried out for a range of system sizes, and the results are scaled to generate the binodals corrected for finite system sizes. Size-corrected interfacial tension is also produced from slab-geometry molecular dynamics simulations.