Actin network heterogeneity tunes activator-inhibitor dynamics at the cell cortex

Biological systems can display diverse patterns of self-organization, even when built on conserved networks of interaction between molecular species. In these cases, reaction-diffusion equations provide a valuable tool to learn how new dynamics could emerge from quantitative tuning of parameters. Bringing these models into quantitative correspondence with biological data remains an outstanding challenge, especially when the biology manifests heterogeneities that are difficult to account for mathematically. One particular example occurs in cell biology, where the membrane-bound regulatory protein RhoA interacts with the filamentous actin cortex in an activator-inhibitor loop. Though this core biochemical circuit is conserved across multiple cell types in different organisms, it produces different patterns of RhoA activity in different contexts, from traveling waves in starfish to transient pulses in C. elegans . To understand how this variation emerges, we develop an activator-inhibitor model that accounts explicitly for actin assembly and heterogeneity. We use the model to infer the parameters of actin dynamics responsible for each experimental phenotype, and find them to agree with independent measurements in most cases. To interpret parameter groupings, we use a simplified model to demonstrate how spatial and orientational randomness separately induce variation in spatiotemporal dynamics of RhoA activity. This work sheds light on how phenotypic diversity arises from heterogeneity and anisotropy, with important implications for the next generation of activator-inhibitor models.