Models of collective cell behaviour with crowding effects : comparing lattice-based and lattice-free approaches

Individual-based models describing the migration and proliferation of a population of cells frequently restrict the cells to a predefined lattice. An implicit assumption of this type of lattice based model is that a proliferative population will always eventually fill the lattice. Here we develop a new lattice-free individual-based model that incorporates cell-to-cell crowding effects. We also derive approximate mean-field descriptions for the lattice-free model in two special cases motivated by commonly used experimental setups. Lattice-free simulation results are compared to these mean-field descriptions and to a corresponding lattice-based model. Data from a proliferation experiment is used to estimate the parameters for the new model, including the cell proliferation rate, showing that the model fits the data well. An important aspect of the lattice-free model is that the confluent cell density is not predefined, as with lattice-based models, but an emergent model property. As a consequence of the more realistic, irregular configuration of cells in the lattice-free model, the population growth rate is much slower at high cell densities and the population cannot reach the same confluent density as an equivalent lattice-based model.

Sensitivity of edge detection methods for quantifying cell migration assays

Quantitative imaging methods to analyze cell migration assays are not standardized. Here we present a suite of two–dimensional barrier assays describing the collective spreading of an initially–confined population of 3T3 fibroblast cells. To quantify the motility rate we apply two different automatic image detection methods to locate the position of the leading edge of the spreading population after 24, 48 and 72 hours. These results are compared with a manual edge detection method where we systematically vary the detection threshold. Our results indicate that the observed spreading rates are very sensitive to the choice of image analysis tools and we show that a standard measure of cell migration can vary by as much as 25% for the same experimental images depending on the details of the image analysis tools. Our results imply that it is very difficult, if not impossible, to meaningfully compare previously published measures of cell migration since previous results have been obtained using different image analysis techniques and the details of these techniques are not always reported. Using a mathematical model, we provide a physical interpretation of our edge detection results. The physical interpretation is important since edge detection algorithms alone do not specify any physical measure, or physical definition, of the leading edge of the spreading population. Our modeling indicates that variations in the image threshold parameter correspond to a consistent variation in the local cell density. This means that varying the threshold parameter is equivalent to varying the location of the leading edge in the range of approximately 1–5% of the maximum cell density.

Mathematical models for collective cell spreading

Collective cell spreading is frequently observed in development, tissue repair and disease progression. Mathematical modelling used in conjunction with experimental investigation can provide key insights into the mechanisms driving the spread of cell populations. In this study, we investigated how experimental and modelling frameworks can be used to identify several key features underlying collective cell spreading. In particular, we were able to independently quantify the roles of cell motility and cell proliferation in a spreading cell population, and investigate how these roles are influenced by factors such as the initial cell density, type of cell population and the assay geometry.

Mechanistic investigations into interactions between IGF-I and IGFBPs and their impact on facilitating cell migration on vitronectin

Numerous studies have reported links between insulin-like growth factors (IGFs) and the extra-cellular matrix protein vitronectin (VN). We ourselves have reported that IGF-I binds to VN via IGF-binding proteins (IGFBPs) to stimulate HaCaT and MCF-7 cell migration. Here, we detail the functional evaluation of IGFBP-1, -2, -3, -4 and -6 in the presence and absence of IGF-I and VN. The data presented here, combined with our prior data on IGFBP-5, suggest that IGFBP-3, -4 and -5 are the most effective at stimulating cell migration in combination with IGF-I and VN. In addition, we demonstrate that different regions within IGFBP-3 and -4 are critical for complex formation. Furthermore, we examine whether multi-protein complexes of IGF-I and IGFBPs associated with fibronectin and collagen IV are also able to enhance functional biological responses.

Models of collective cell spreading with variable cell aspect ratio : a motivation for degenerate diffusion models

Continuum diffusion models are often used to represent the collective motion of cell populations. Most previous studies have simply used linear diffusion to represent collective cell spreading, while others found that degenerate nonlinear diffusion provides a better match to experimental cell density profiles. In the cell modeling literature there is no guidance available with regard to which approach is more appropriate for representing the spreading of cell populations. Furthermore, there is no knowledge of particular experimental measurements that can be made to distinguish between situations where these two models are appropriate. Here we provide a link between individual-based and continuum models using a multi-scale approach in which we analyze the collective motion of a population of interacting agents in a generalized lattice-based exclusion process. For round agents that occupy a single lattice site, we find that the relevant continuum description of the system is a linear diffusion equation, whereas for elongated rod-shaped agents that occupy L adjacent lattice sites we find that the relevant continuum description is connected to the porous media equation (pme). The exponent in the nonlinear diffusivity function is related to the aspect ratio of the agents. Our work provides a physical connection between modeling collective cell spreading and the use of either the linear diffusion equation or the pme to represent cell density profiles. Results suggest that when using continuum models to represent cell population spreading, we should take care to account for variations in the cell aspect ratio because different aspect ratios lead to different continuum models.

A cell migration device that maintains a defined surface with no cellular damage during wound edge generation

Studying the rate of cell migration provides insight into fundamental cell biology as well as a tool to assess the functionality of synthetic surfaces and soluble environments used in tissue engineering. The traditional tools used to study cell migration include the fence and wound healing assays. In this paper we describe the development of a microchannel based device for the study of cell migration on defined surfaces. We demonstrate that this device provides a superior tool, relative to the previously mentioned assays, for assessing the propagation rate of cell wave fronts. The significant advantage provided by this technology is the ability to maintain a virgin surface prior to the commencement of the cell migration assay. Here, the device is used to assess rates of mouse fibroblasts (NIH 3T3) and human osteosarcoma (SaOS2) cell migration on surfaces functionalized with various extracellular matrix proteins as a demonstration that confining cell migration within a microchannel produces consistent and robust data. The device design enables rapid and simplistic assessment of multiple repeats on a single chip, where surfaces have not been previously exposed to cells or cellular secretions.

Gonadotropin-induced ovarian cancer cell migration and proliferation require extracellular signal-regulated kinase 1/2 activation regulated by calcium and protein kinase Cδ

The gonadotropin hypothesis proposes that elevated serum gonadotropin levels may increase the risk of epithelial ovarian cancer (EOC). We have studied the effect of treating EOC cell lines (OV207 and OVCAR-3) with FSH or LH. Both gonadotropins activated the mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase 1/2 (ERK1/2) pathway and increased cell migration that was inhibited by the MAPK 1 inhibitor PD98059. Both extra- and intracellular calcium ion signalling were implicated in gonadotropin-induced ERK1/2 activation as treatment with either the calcium chelator EGTA or an inhibitor of intracellular calcium release, dantrolene, inhibited gonadotropin-induced ERK1/2 activation. Verapamil was also inhibitory, indicating that gonadotropins activate calcium influx via L-type voltage-dependent calcium channels. The cAMP/protein kinase A (PKA) pathway was not involved in the mediation of gonadotropin action in these cells as gonadotropins did not increase intracellular cAMP formation and inhibition of PKA did not affect gonadotropin-induced phosphorylation of ERK1/2. Activation of ERK1/2 was inhibited by the protein kinase C (PKC) inhibitor GF 109203X as well as by the PKCδ inhibitor rottlerin, and downregulation of PKCδ was inhibited by small interfering RNA (siRNA), highlighting the importance of PKCδ in the gonadotropin signalling cascade. Furthermore, in addition to inhibition by PD98059, gonadotropin-induced ovarian cancer cell migration was also inhibited by verapamil, GF 109203X and rottlerin. Similarly, gonadotropin-induced proliferation was inhibited by PD98059, verapamil, GF 109203X and PKCδ siRNA. Taken together, these results demonstrate that gonadotropins induce both ovarian cancer cell migration and proliferation by activation of ERK1/2 signalling in a calcium- and PKCδ-dependent manner.

Models of collective cell motion for cell populations with different aspect ratio : diffusion, proliferation and travelling waves

Continuum, partial differential equation models are often used to describe the collective motion of cell populations, with various types of motility represented by the choice of diffusion coefficient, and cell proliferation captured by the source terms. Previously, the choice of diffusion coefficient has been largely arbitrary, with the decision to choose a particular linear or nonlinear form generally based on calibration arguments rather than making any physical connection with the underlying individual-level properties of the cell motility mechanism. In this work we provide a new link between individual-level models, which account for important cell properties such as varying cell shape and volume exclusion, and population-level partial differential equation models. We work in an exclusion process framework, considering aligned, elongated cells that may occupy more than one lattice site, in order to represent populations of agents with different sizes. Three different idealizations of the individual-level mechanism are proposed, and these are connected to three different partial differential equations, each with a different diffusion coefficient; one linear, one nonlinear and degenerate and one nonlinear and nondegenerate. We test the ability of these three models to predict the population level response of a cell spreading problem for both proliferative and nonproliferative cases. We also explore the potential of our models to predict long time travelling wave invasion rates and extend our results to two dimensional spreading and invasion. Our results show that each model can accurately predict density data for nonproliferative systems, but that only one does so for proliferative systems. Hence great care must be taken to predict density data for with varying cell shape.

The cell surface glycoprotein CUB domain-containing protein 1 (CDCP1) contributes to epidermal growth factor receptor-mediated cell migration

Epidermal growth factor (EGF) activation of the EGF receptor (EGFR) is an important mediator of cell migration, and aberrant signaling via this system promotes a number of malignancies including ovarian cancer. We have identified the cell surface glycoprotein CDCP1 as a key regulator of EGF/EGFR-induced cell migration. We show that signaling via EGF/EGFR induces migration of ovarian cancer Caov3 and OVCA420 cells with concomitant up-regulation of CDCP1 mRNA and protein. Consistent with a role in cell migration CDCP1 relocates from cell-cell junctions to punctate structures on filopodia after activation of EGFR. Significantly, disruption of CDCP1 either by silencing or the use of a function blocking antibody efficiently reduces EGF/EGFR-induced cell migration of Caov3 and OVCA420 cells. We also show that up-regulation of CDCP1 is inhibited by pharmacological agents blocking ERK but not Src signaling, indicating that the RAS/RAF/MEK/ERK pathway is required downstream of EGF/EGFR to induce increased expression of CDCP1. Our immunohistochemical analysis of benign, primary, and metastatic serous epithelial ovarian tumors demonstrates that CDCP1 is expressed during progression of this cancer. These data highlight a novel role for CDCP1 in EGF/EGFR-induced cell migration and indicate that targeting of CDCP1 may be a rational approach to inhibit progression of cancers driven by EGFR signaling including those resistant to anti-EGFR drugs because of activating mutations in the RAS/RAF/MEK/ERK pathway.

Travelling waves for a velocity–jump model of cell migration and proliferation

Cell invasion, characterised by moving fronts of cells, is an essential aspect of development, repair and disease. Typically, mathematical models of cell invasion are based on the Fisher–Kolmogorov equation. These traditional parabolic models can not be used to represent experimental measurements of individual cell velocities within the invading population since they imply that information propagates with infinite speed. To overcome this limitation we study combined cell motility and proliferation based on a velocity–jump process where information propagates with finite speed. The model treats the total population of cells as two interacting subpopulations: a subpopulation of left–moving cells, $L(x,t)$, and a subpopulation of right–moving cells, $R(x,t)$. This leads to a system of hyperbolic partial differential equations that includes a turning rate, $\Lambda \ge 0$, describing the rate at which individuals in the population change direction of movement. We present exact travelling wave solutions of the system of partial differential equations for the special case where $\Lambda = 0$ and in the limit that $\Lambda \to \infty$. For intermediate turning rates, $0 < \Lambda < \infty$, we analyse the travelling waves using the phase plane and we demonstrate a transition from smooth monotone travelling waves to smooth nonmonotone travelling waves as $\Lambda$ decreases through a critical value $\Lambda_{crit}$. We conclude by providing a qualitative comparison between the travelling wave solutions of our model and experimental observations of cell invasion. This comparison indicates that the small $\Lambda$ limit produces results that are consistent with experimental observations.

Cell migration and proliferation on homogeneous and non-homogeneous domains : modelling on the scale of individuals and populations

Cell migration is a behaviour critical to many key biological effects, including wound healing, cancerous cell invasion and morphogenesis, the development of an organism from an embryo. However, given that each of these situations is distinctly different and cells are extremely complicated biological objects, interest lies in more basic experiments which seek to remove conflating factors and present a less complex environment within which cell migration can be experimentally examined. These include in vitro studies like the scratch assay or circle migration assay, and ex vivo studies like the colonisation of the hindgut by neural crest cells. The reduced complexity of these experiments also makes them much more enticing as problems to mathematically model, like done here. The primary goal of the mathematical models used in this thesis is to shed light on which cellular behaviours work to generate the travelling waves of invasion observed in these experiments, and to explore how variations in these behaviours can potentially predict differences in this invasive pattern which are experimentally observed when cell types or chemical environment are changed. Relevant literature has already identified the difficulty of distinguishing between these behaviours when using traditional mathematical biology techniques operating on a macroscopic scale, and so here a sophisticated individual-cell-level model, an extension of the Cellular Potts Model (CPM), is been constructed and used to model a scratch assay experiment. This model includes a novel mechanism for dealing with cell proliferations that allowed for the differing properties of quiescent and proliferative cells to be implemented into their behaviour. This model is considered both for its predictive power and used to make comparisons with the travelling waves which result in more traditional macroscopic simulations. These comparisons demonstrate a surprising amount of agreement between the two modelling frameworks, and suggest further novel modifications to the CPM that would allow it to better model cell migration. Considerations of the model’s behaviour are used to argue that the dominant effect governing cell migration (random motility or signal-driven taxis) likely depends on the sort of invasion demonstrated by cells, as easily seen by microscopic photography. Additionally, a scratch assay simulated on a non-homogeneous domain consisting of a ’fast’ and ’slow’ region is also used to further differentiate between these different potential cell motility behaviours. A heterogeneous domain is a novel situation which has not been considered mathematically in this context, nor has it been constructed experimentally to the best of the candidate’s knowledge. Thus this problem serves as a thought experiment used to test the conclusions arising from the simulations on homogeneous domains, and to suggest what might be observed should this non-homogeneous assay situation be experimentally realised. Non-intuitive cell invasion patterns are predicted for diffusely-invading cells which respond to a cell-consumed signal or nutrient, contrasted with rather expected behaviour in the case of random-motility-driven invasion. The potential experimental observation of these behaviours is demonstrated by the individual-cell-level model used in this thesis, which does agree with the PDE model in predicting these unexpected invasion patterns. In the interest of examining such a case of a non-homogeneous domain experimentally, some brief suggestion is made as to how this could be achieved.

Quantifying alterations in cell migration : tracking fluorescently-tagged migrating cells by FACs and live-imaging

Cell migration is fundamental to many different physiological processes including embryonic development, inflammation and wound healing. Given the range and importance cell migration plays a number of assays have been developed to measure different aspects of cell migration. Here we describe two different methods to analyze cell migration. The first method analyzes the migration of fluorescently tagged cells using Boyden chambers and FACs and the second looks at migration properties using time-lapse microscopy.

Lattice-free descriptions of collective motion with crowding and adhesion

Cell-to-cell adhesion is an important aspect of malignant spreading that is often observed in images from the experimental cell biology literature. Since cell-to-cell adhesion plays an important role in controlling the movement of individual malignant cells, it is likely that cell-to-cell adhesion also influences the spatial spreading of populations of such cells. Therefore, it is important for us to develop biologically realistic simulation tools that can mimic the key features of such collective spreading processes to improve our understanding of how cell-to-cell adhesion influences the spreading of cell populations. Previous models of collective cell spreading with adhesion have used lattice-based random walk frameworks which may lead to unrealistic results, since the agents in the random walk simulations always move across an artificial underlying lattice structure. This is particularly problematic in high-density regions where it is clear that agents in the random walk align along the underlying lattice, whereas no such regular alignment is ever observed experimentally. To address these limitations, we present a lattice-free model of collective cell migration that explicitly incorporates crowding and adhesion. We derive a partial differential equation description of the discrete process and show that averaged simulation results compare very well with numerical solutions of the partial differential equation.

Selective regulation of p38β protein and signaling by integrin-linked kinase mediates bladder cancer cell migration

Integrin-linked kinase (ILK) and p38MAPK are protein kinases that transduce extracellular signals regulating cell migration and actin cytoskeletal organization. ILK-dependent regulation of p38MAPK is critical for mammalian kidney development and in smooth muscle cell migration, however, specific p38 isoforms has not been previously examined in ILK-regulated responses. Signaling by ILK and p38MAPK is often dysregulated in bladder cancer, and here we report a strong positive correlation between protein levels of ILK and p38β, which is the predominant isoform found in bladder cancer cells, as well as in patient-matched normal bladder and tumor samples. Knockdown by RNA interference of either p38β or ILK disrupts serum-induced, Rac1-dependent migration and actin cytoskeletal organization in bladder cancer cells. Surprisingly, ILK knockdown causes the selective reduction in p38β cellular protein level, without inhibiting p38β messenger RNA (mRNA) expression. The loss of p38β protein in ILK-depleted cells is partially rescued by the 26S proteasomal inhibitor MG132. Using co-precipitation and bimolecular fluorescent complementation assays, we find that ILK selectively forms cytoplasmic complexes with p38β. In situ proximity ligation assays further demonstrate that serum-stimulated assembly of endogenous ILK–p38β complexes is sensitive to QLT-0267, a small molecule ILK kinase inhibitor. Finally, inhibition of ILK reduces the amplitude and period of serum-induced activation of heat shock protein 27 (Hsp27), a target of p38β implicated in actin cytoskeletal reorganization. Our work identifies Hsp27 as a novel target of ILK–p38β signaling complexes, playing a key role in bladder cancer cell migration.