A numerical model of cellular blebbing: a volume-conserving, fluid-structure interaction model of the entire cell

In animal cells, blebs are smooth, quasi-hemispherical protrusions of the plasma membrane that form when a section of the membrane detaches from the underlying actin cytoskeleton and is inflated by flowing cytosol. The mechanics behind this common cellular activity are not yet clear. As a first step in the development of a full computational framework, we present a numerical model of overall cell behavior based upon the interaction between a background Newtonian-fluid cytosol and elastic structures modeling the membrane and filaments. The detailed micromechanics of the cytoskeletal network are the subject of future work. Here, the myosin-driven contraction of the actin network is modeled through stressed elastic filaments. Quantitative models of cytoskeletal micromechanics and biochemistry require accurate estimates of local stress and flow conditions. The main contribution of this paper is the development of a computationally efficient fluid-structure interaction model based on operator splitting, to furnish this data. Cytosol volume conservation (as supported by experimental evidence) is enforced through an intermediate energy minimization step. Realistic bleb formation and retraction is observed from this model, offering an alternative formulation to positing complex continuum behavior of the cytoplasm (e.g. poroelastic model of Charras et al., 2008).