Abstract:Biomolecular condensates are important contributors to the internal organization of the cell material. While initially described as liquid-like droplets, the term biomolecular condensates is now used to describe a diversity of condensed phase assemblies with material properties extending from low to high viscous liquids, gels, and even glasses. Because the material properties of condensates are determined by the intrinsic behavior of their molecules, characterizing such properties is integral to rationalizing the molecular mechanisms that dictate their functions and roles in health and disease. Here, we apply and compare three distinct computational methods to measure the viscoelasticity of biomolecular condensates in molecular simulations. These methods are the Green-Kubo (GK) relation, the oscillatory shear (OS) technique, and the bead tracking (BT) method. We find that, although all of these methods provide consistent results for the viscosity of the condensates, the GK and OS techniques outperform the BT method in terms of computational efficiency and statistical uncertainty. We thus apply the GK and OS techniques for a set of 12 different protein/RNA systems using a sequence-dependent coarse-grained model. Our results reveal a strong correlation between condensate viscosity and density, as well as with protein/RNA length and the number of stickers vs spacers in the amino acid protein sequence. Moreover, we couple the GK and the OS technique to nonequilibrium molecular dynamics simulations that mimic the progressive liquid-to-gel transition of protein condensates due to the accumulation of interprotein β-sheets. We compare the behavior of three different protein condensates, i.e., those formed by either hnRNPA1, FUS, or TDP-43 proteins, whose liquid-to-gel transitions are associated with the onset of amyotrophic lateral sclerosis and frontotemporal dementia. We find that both the GK and OS techniques successfully predict the transition from functional liquid-like behavior to kinetically arrested states once the network of interprotein β-sheets has percolated through the condensates. Overall, our work provides a comparison of different modeling rheological techniques to assess the viscosity of biomolecular condensates, a critical magnitude that provides information on the behavior of biomolecules inside condensates.

Time-Dependent Material Properties of Aging Biomolecular Condensates from Different Viscoelasticity Measurements in Molecular Dynamics Simulations

Location and Concentration of Aromatic‐Rich Segments Dictates the Percolating Inter‐Molecular Network and Viscoelastic Properties of Ageing Condensates

Theory of rheology and aging of protein condensates

Direct computations of viscoelastic moduli of biomolecular condensates

Sequence-encoded Spatiotemporal Dependence of Viscoelasticity of Protein Condensates Using Computational Microrheology

Microscopic model for aging of biocondensates

Viscoelasticity of globular protein-based biomolecular condensates

Mesoscale properties of biomolecular condensates emerging from protein chain dynamics

Determinants of viscoelasticity and flow activation energy in biomolecular condensates

Viscoelasticity of biomolecular condensates conforms to the Jeffreys model

Formation of Multi-Compartment Condensates through Aging of Protein-RNA Condensates

Protein condensates as aging Maxwell fluids

Optical characterization of molecular interaction strength in protein condensates

Scanning probe microscopy elucidates gelation and rejuvenation of biomolecular condensates

Dynamically arrested condensate fusion creates complex structures with varying material properties

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

Biomolecular condensates form spatially inhomogeneous network fluids

Active learning of the thermodynamics-dynamics tradeoff in protein condensates

Charged peptides enriched in aromatic residues decelerate condensate ageing driven by cross-beta-sheet formation

Biomolecular condensates undergo a generic shear-mediated liquid-to-solid transition