Forced and spontaneous symmetry breaking in cell polarization

Computational tools for cellular scale biophysics

Mathematical models are indispensable for disentangling the interactions through which biological components work together to generate the forces and flows that position, mix, and distribute proteins, nutrients, and organelles within the cell. To illuminate the ever more specific questions studied at the edge of biological inquiry, such models inevitably become more complex. Solving, simulating, and learning from these more realistic models requires the development of new analytic techniques, numerical methods, and scalable software. In this review, we discuss some recent developments in tools for understanding how large numbers of cytoskeletal filaments, driven by molecular motors and interacting with the cytoplasm and other structures in their environment, generate fluid flows, instabilities, and material deformations which help drive crucial cellular processes.

Numerical investigations of the bulk-surface wave pinning model

The bulk-surface wave pinning model is a reaction-diffusion system for studying cell polarisation. It is constituted by a surface reaction-diffusion equation, coupled to a bulk diffusion equation with a non-linear boundary condition. Cell polarisation arises as the surface component develops specific patterns. Since proteins diffuse much faster in the cell interior than on the membrane, in the literature, the bulk component is often assumed to be spatially homogeneous. Therefore, the model can be reduced to a single surface equation. However, in real applications a spatially non-uniform bulk component might be an important player to take into account. In this paper, we study, through numerical computations, the role of the bulk component and, more specifically, how different bulk diffusion rates might affect the polarisation response. We find that the bulk component is indeed a key factor in determining the surface polarisation response. Moreover, for certain geometries, it is the spatial heterogeneity of the bulk component that triggers the polarisation response, which might not be possible in a reduced model. Understanding how polarisation depends on bulk diffusivity might be crucial when studying models of migrating cells, which are naturally subject to domain deformation.