Damped and persistent oscillations in a simple model of cell crawling

A continuum mechanical model of cell motion driven by a biphasic traction stress

The aim of this paper is to place the cell locomotion problem within the general framework of classical continuum mechanics, and while doing so, to account for the deformation of the actin network in the cytoskeleton; the myosin activity on the lamellum including its effect on depolymerization at the trailing edge; model the stress-dependent driving forces and kinetic laws controlling polymerization at the leading edge, depolymerization at the trailing edge and ATP hydrolysis consistently with the dissipation inequality; and, based on the observations in Gardel et al. (Gardel et al. 2008 J. Cell Biol. 183, 999-1005 (doi:10.1083/jcb.200810060)), include a biphasic velocity-dependent traction stress acting on the actin network. While we chose certain specific models for each of these, in part to allow for an analytical solution, the generality of the framework allows one to readily introduce different constitutive laws to describe these phenomena as might be needed, for example, to study some different type of cells. As described in §5, the predictions of the model compare well with observations such as the magnitude of the very different actin retrograde speeds in the lamellum and lamellipodium including their jump at the interface, the magnitude of the cell speed, and the relative lengths of the lamellipodium and lamellum.

Estimation Algorithm for a Hybrid PDE–ODE Model Inspired by Immunocompetent Cancer-on-Chip Experiment

The present work is motivated by the development of a mathematical model mimicking the mechanisms observed in lab-on-chip experiments, made to reproduce on microfluidic chips the in vivo reality. Here we consider the Cancer-on-Chip experiment where tumor cells are treated with chemotherapy drug and secrete chemical signals in the environment attracting multiple immune cell species. The in silico model here proposed goes towards the construction of a “digital twin” of the experimental immune cells in the chip environment to better understand the complex mechanisms of immunosurveillance. To this aim, we develop a tumor-immune microfluidic hybrid PDE–ODE model to describe the concentration of chemicals in the Cancer-on-Chip environment and immune cells migration. The development of a trustable simulation algorithm, able to reproduce the immunocompetent dynamics observed in the chip, requires an efficient tool for the calibration of the model parameters. In this respect, the present paper represents a first methodological work to test the feasibility and the soundness of the calibration technique here proposed, based on a multidimensional spline interpolation technique for the time-varying velocity field surfaces obtained from cell trajectories.

A discrete in continuous mathematical model of cardiac progenitor cells formation and growth as spheroid clusters (Cardiospheres)

We propose a discrete in continuous mathematical model describing the in vitro growth process of biophsy-derived mammalian cardiac progenitor cells growing as clusters in the form of spheres (Cardiospheres). The approach is hybrid: discrete at cellular scale and continuous at molecular level. In the present model, cells are subject to the self-organizing collective dynamics mechanism and, additionally, they can proliferate and differentiate, also depending on stochastic processes. The two latter processes are triggered and regulated by chemical signals present in the environment. Numerical simulations show the structure and the development of the clustered progenitors and are in a good agreement with the results obtained from in vitro experiments.

Role of intracellular poroelasticity on freezing-induced deformation of cells in engineered tissues

Freezing of biomaterials is important in a wide variety of biomedical applications, including cryopreservation and cryosurgeries. For the success of these applications to various biomaterials, biophysical mechanisms, which determine freezing-induced changes in cells and tissues, need to be well understood. Specifically, the significance of the intracellular mechanics during freezing is not well understood. Thus, we hypothesize that cells interact during freezing with the surroundings such as suspension media and the extracellular matrix (ECM) via two distinct but related mechanisms-water transport and cytoskeletal mechanics. The underlying rationale is that the cytoplasm of the cells has poroelastic nature, which can regulate both cellular water transport and cytoskeletal mechanics. A poroelasticity-based cell dehydration model is developed and confirmed to provide insight into the effects of the hydraulic conductivity and stiffness of the cytoplasm on the dehydration of cells in suspension during freezing. We further investigated the effect of the cytoskeletal structures on the cryoresponse of cells embedded in the ECM by measuring the spatio-temporal intracellular deformation with dermal equivalent as a model tissue. The freezing-induced change in cell, nucleus and cytoplasm volume was quantified, and the possible mechanism of the volumetric change was proposed. The results are discussed considering the hierarchical poroelasticity of biological tissues.

Effects of adhesion dynamics and substrate compliance on the shape and motility of crawling cells

Computational modeling of eukaryotic cells moving on substrates is an extraordinarily complex task: many physical processes, such as actin polymerization, action of motors, formation of adhesive contacts concomitant with both substrate deformation and recruitment of actin etc., as well as regulatory pathways are intertwined. Moreover, highly nontrivial cell responses emerge when the substrate becomes deformable and/or heterogeneous. Here we extended a computational model for motile cell fragments, based on an earlier developed phase field approach, to account for explicit dynamics of adhesion site formation, as well as for substrate compliance via an effective elastic spring. Our model displays steady motion vs. stick-slip transitions with concomitant shape oscillations as a function of the actin protrusion rate, the substrate stiffness, and the rates of adhesion. Implementing a step in the substrate's elastic modulus, as well as periodic patterned surfaces exemplified by alternating stripes of high and low adhesiveness, we were able to reproduce the correct motility modes and shape phenomenology found experimentally. We also predict the following nontrivial behavior: the direction of motion of cells can switch from parallel to perpendicular to the stripes as a function of both the adhesion strength and the width ratio of adhesive to non-adhesive stripes.