Computational rheometry of viscoelastic networks: From random graphs to biomolecular condensates

Multivalent biomacromolecules including multi-domain and intrinsically disordered proteins form biomolecular condensates via reversible phase transitions. Condensates are viscoelastic materials that display composition-specific rheological properties and responses to mechanical forces. Graph-based descriptions of microstructures can be combined with computational rheometry to model the outcomes of passive and active mechanical measurements. We consider two types of network models for microstructures. In the Jeffreys model, each edge in the network is a Jeffreys element. In the Stokes-Maxwell model, each edge is a Maxwell element that is embedded in an incompressible viscous fluid that can undergo Stokes flow. We describe results from comparative assessments of the two models for individual elements, ordered lattices, random geometric graphs, structured graphs, and graphs for condensates that are extracted from coarse-grained simulations of disordered proteins. Results from deformation and relaxation tests and flow field analysis reveal how distinct length and time scales contribute to the responses of different types of networks. No single test provides definitive assessments of the connections between material properties and microstructures. Instead, a range of active and passive rheometric tests are essential for distinguishing the responses of different types of networks. Our work establishes computational rheometry as a framework for bridging disparate length and timescales to assess how molecular-scale interactions and dynamics give rise to viscoelastic responses on the mesoscale.