Polarity and mixed-mode oscillations may underlie different patterns of cellular migration

Abstract

In mesenchymal cell motility, several migration patterns have been observed, including directional, exploratory and stationary. Two key members of the Rho-family of GTPases, Rac and Rho, along with an adaptor protein called paxillin, have been particularly implicated in the formation of such migration patterns and in regulating adhesion dynamics. Together, they form a key regulatory network that involves the mutual inhibition exerted by Rac and Rho on each other and the promotion of Rac activation by phosphorylated paxillin. Although this interaction is sufficient to generating wave-pinning that underscores cellular polarization comprised of cellular front (high active Rac) and back (high active Rho), it remains unclear how they interact collectively to induce other modes of migration detected in Chinese hamster Ovary (CHO-K1) cells. We previously developed a 6D reaction-diffusion model describing the interactions of these three proteins (in their active/phosphorylated and inactive/unphosphorylated forms) along with other auxiliary proteins, to decipher their role in generating wave-pinning. In this study, we explored, through computational modeling and image analysis, how differences in timescales within this molecular network can potentially produce the migration patterns in CHO-K1 cells and how switching between them could occur. To do so, the 6D model was reduced to an excitable 4D spatiotemporal model possessing three different timescales. The model produced not only wave-pinning in the presence of diffusion, but also mixed-mode oscillations (MMOs) and relaxation oscillations (ROs). Implementing the model using the Cellular Potts Model (CPM) produced outcomes in which protrusions in cell membrane changed Rac-Rho localization, resulting in membrane oscillations and fast directionality variations similar to those seen in CHO-K1 cells. The latter was assessed by comparing the migration patterns of CHOK1 cells with CPM cells using four metrics: instantaneous cell speed, exponent of mean square-displacement (called -value), directionality ratio and protrusion rate. Variations in migration patterns induced by mutating paxillins serine 273 residue was also captured by the model and detected by a machine classifier, revealing that this mutation alters the dynamics of the system from MMOs to ROs or nonoscillatory behaviour through variation in the concentration of an active form of an adhesion protein called p21-Activated Kinase 1 (PAK). These results thus suggest that MMOs and adhesion dynamics are the key ingredients underlying CHO-K1 cell motility.