Deciphering the dynamics of lamellipodium in a fish keratocytes model

Motile cells depend on an intricate network of feedback loops that are essential in driving cell movement. Integrin-based focal adhesions (FAs) along with actin are the two key factors that mediate such motile behaviour. Together, they generate excitable dynamics that are essential for forming protrusions at the leading edge of the cell and, in certain cases, traveling waves along the membrane. A partial differential equation (PDE) model of a self-organizing lamellipodium in crawling keratocytes has been previously developed to understand how the three spatiotemporal patterns of activity observed in such cells, namely, stalling, waving and smooth motility, are produced. The model consisted of three key variables: the density of barbed actin filaments, newly formed FAs called nascent adhesions (NAs) and VASP, an anti-capping protein that gets sequestered by NAs during maturation. Using parameter sweeping techniques, the distinct regimes of behaviour associated with the three activity patterns were identified. In this study, we convert the PDE model into an ordinary differential equation (ODE) model to examine its excitability properties and determine all the patterns of activity exhibited by this system. Our results reveal that there are two additional regimes not previously identified, including bistability and oscillatory-like type IV excitability (generated by three steady states and their manifolds, rather than limit cycles). These regimes are also present in the PDE model. Applying slow-fast analysis on the ODE model shows that it exhibits a canard explosion through a folded-saddle and that rough motility seen in keratocytes is likely due to noise-dependent motility governed by dynamics near the interface of bistability and type IV excitability. The two parameter bifurcation suggests that the increase in the proportion of rough motion is due to a shift in activity towards the bistable and type IV excitable regimes induced by a decrease in NA maturation rate. Our results thus provide important insight into how microscopic mechanical effects are integrated to produce the observed modes of motility.