Quantitative Biology

We introduce a generic, purely mechanical model for environment sensitive motion of mammalian cells that is applicable to chemotaxis, haptotaxis, and durotaxis as modes of motility. It is able to theoretically explain all relevant experimental observations, in particular, the high efficiency of motion, the behavior on inhomogeneous substrates, and the fixation of the lagging pole during motion. Furthermore, our model predicts that efficiency of motion in following a gradient depends on cell geometry (with more elongated cells being more efficient).

We analyse a generic motility model, with the motility mechanism arising by contractile stress due to the interaction of myosin and actin. A hydrodynamic active polar gel theory is used to model the cytoplasm of a cell and is combined with a Helfrich-type model to account for membrane properties. The overall model allows to consider motility without the necessity for local adhesion. Besides a detailed numerical approach together with convergence studies for the highly nonlinear free boundary problem, we also compare the induced flow field of the motile cell with that of classical squirmer models and identify the motile cell as a puller or pusher, depending on the strength of the myosin-actin interactions.

Cell polarization and directional cell migration can display random, persistent and oscillatory dynamic patterns. However, it is not clear if these polarity patterns can be explained by the same underlying regulatory mechanism. Here, we show that random, persistent and oscillatory migration accompanied by polarization can simultaneously occur in populations of melanoma cells derived from tumors with different degrees of aggressiveness. We demonstrate that all these patterns and the probabilities of their occurrence are quantitatively accounted for by a simple mechanism involving a spatially distributed, mechano-chemical feedback coupling the dynamically changing extracellular matrix (ECM)-cell contacts to the activation of signaling downstream of the Rho-family small GTPases. This mechanism is supported by a predictive mathematical model and extensive experimental validation, and can explain previously reported results for diverse cell types. In melanoma, this mechanism also accounts for the effects of genetic and environmental perturbations, including mutations linked to invasive cell spread. The resulting mechanistic understanding of cell polarity quantitatively captures the relationship between population variability and phenotypic plasticity, with the potential to account for a wide variety of cell migration states in diverse pathological and physiological conditions.

The skin forms an efficient barrier against the environment, and rapid cutaneous wound healing after injury is therefore essential. Healing of the uppermost layer of the skin, the epidermis, involves collective migration of keratinocytes, which requires coordinated polarization of the cells. To study this process, we developed a model that allows analysis of live-cell images of migrating keratinocytes in culture based on a small number of parameters, including the radius of the cells, their mass and their polarization. This computational approach allowed the analysis of cell migration at the front of the wound and a reliable identification and quantification of the impaired polarization and migration of keratinocytes from mice lacking fibroblast growth factors 1 and 2 an established model of impaired healing. Therefore, our modeling approach is suitable for large-scale analysis of migration phenotypes of cells with specific genetic defects or upon treatment with different pharmacological agents.

This work presents the development, analysis and numerical simulations of a biophysical model for 3D cell deformation and movement, which couples biochemical reactions and biomechanical forces. We propose a mechanobiochemical model which considers the actin filament network as a viscoelastic and contractile gel. The mechanical properties are modelled by a force balancing equation for the displacements, the pressure and concentration forces are driven by actin and myosin dynamics, and these are in turn modelled by a system of reaction-diffusion equations on a moving cell domain. The biophysical model consists of highly non-linear partial differential equations whose analytical solutions are intractable. To obtain approximate solutions to the model system, we employ the moving grid finite element method. The numerical results are supported by linear stability theoretical results close to bifurcation points during the early stages of cell migration. Numerical simulations exhibited show both simple and complex cell deformations in 3-dimensions that include cell expansion, cell protrusion and cell contraction. The computational framework presented here sets a strong foundation that allows to study more complex and experimentally driven reaction-kinetics involving actin, myosin and other molecular species that play an important role in cell movement and deformation.

The conventional cancer stem cell (CSC) theory indicates a hierarchy of CSCs and non-stem cancer cells (NSCCs), that is, CSCs can differentiate into NSCCs but not vice versa. However, an alternative paradigm of CSC theory with reversible cell plasticity among cancer cells has received much attention very recently. Here we present a generalized multi-phenotypic cancer model by integrating cell plasticity with the conventional hierarchical structure of cancer cells. We prove that under very weak assumption, the nonlinear dynamics of multi-phenotypic proportions in our model has only one stable steady state and no stable limit cycle. This result theoretically explains the phenotypic equilibrium phenomena reported in various cancer cell lines. Furthermore, according to the transient analysis of our model, it is found that cancer cell plasticity plays an essential role in maintaining the phenotypic diversity in cancer especially during the transient dynamics. Two biological examples with experimental data show that the phenotypic conversions from NCSSs to CSCs greatly contribute to the transient growth of CSCs proportion shortly after the drastic reduction of it. In particular, an interesting overshooting phenomenon of CSCs proportion arises in three-phenotypic example. Our work may pave the way for modeling and analyzing the multi-phenotypic cell population dynamics with cell plasticity.

We propose a multiscale model of the invasion of the extracellular matrix by two types of cancer cells, the differentiated cancer cells and the cancer stem cells. We assume that the epithelial mesenchymal-like transition between them is driven primarily by the epidermal growth factors. We moreover take into account the transidifferentiation program of the cancer stem cells and the cancer associated fibroblast cells as well as the fibroblast-driven remodelling of the extracellular matrix. The proposed haptotaxis model combines the macroscopic phenomenon of the invasion of the extracellular matrix with the microscopic dynamics of the epidermal growth factors. We analyse our model in a component-wise manner and compare our findings with the literature. We investigate pathological situations regarding the epidermal growth factors and accordingly propose "mathematical-treatment" scenarios to control the aggressiveness of the tumour.

We propose a new model for the emergence of blood capillary networks. We assimilate the tissue and extra cellular matrix as a porous medium, using Darcy's law for describing both blood and intersticial fluid flows. Oxygen obeys a convection-diffusion-reaction equation describing advection by the blood, diffusion and consumption by the tissue. Discrete agents named capillary elements and modelling groups of endothelial cells are created or deleted according to different rules involving the oxygen concentration gradient, the blood velocity, the sheer stress or the capillary element density. Once created, a capillary element locally enhances the hydraulic conductivity matrix, contributing to a local increase of the blood velocity and oxygen flow. No connectivity between the capillary elements is imposed. The coupling between blood, oxygen flow and capillary elements provides a positive feedback mechanism which triggers the emergence of a network of channels of high hydraulic conductivity which we identify as new blood capillaries. We provide two different, biologically relevant geometrical settings and numerically analyze the influence of each of the capillary creation mechanism in detail. All mechanisms seem to concur towards a harmonious network but the most important ones are those involving oxygen gradient and sheer stress. A detailed discussion of this model with respect to the literature and its potential future developments concludes the paper.

The purpose of this work is to propose a non-Markovian and nonlinear model of subdiffusive transport that involves adhesion affects the cells escape rates form position x, with chemotaxis. This leads the escape rates to be dependent on the particles density at the neighbours as well as the chemotactic gradient. We systematically derive subdiffusive fractional master equation, then we consider the diffusive limit of the fractional master equation. We finally solve the resulted fractional subdiffusive master equation stationery and analyse the role of adhesion in the resulted stationary density.

We present a novel mathematical model of heterogeneous cell proliferation where the total population consists of a subpopulation of slow-proliferating cells and a subpopulation of fast-proliferating cells. The model incorporates two cellular processes, asymmetric cell division and induced switching between proliferative states, which are important determinants for the heterogeneity of a cell population. As motivation for our model we provide experimental data that illustrate the induced-switching process. Our model consists of a system of two coupled delay differential equations with distributed time delays and the cell densities as functions of time. The distributed delays are bounded and allow for the choice of delay kernel. We analyse the model and prove the non-negativity and boundedness of solutions, the existence and uniqueness of solutions, and the local stability characteristics of the equilibrium points. We find that the parameters for induced switching are bifurcation parameters and therefore determine the long-term behaviour of the model. Numerical simulations illustrate and support the theoretical findings, and demonstrate the primary importance of transient dynamics for understanding the evolution of many experimental cell populations.