Mathematical modelling and numerical simulations of actin dynamics in the eukaryotic cell

Modeling multicellular dynamics regulated by extracellular-matrix-mediated mechanical communication via active particles with polarized effective attraction

Collective cell migration is crucial to many physiological and pathological processes such as embryo development, wound healing, and cancer invasion. Recent experimental studies have indicated that the active traction forces generated by migrating cells in a fibrous extracellular matrix (ECM) can mechanically remodel the ECM, giving rise to bundlelike mesostructures bridging individual cells. Such fiber bundles also enable long-range propagation of cellular forces, leading to correlated migration dynamics regulated by the mechanical communication among the cells. Motivated by these experimental discoveries, we develop an active-particle model with polarized effective attractions (APPA) to investigate emergent multicellular migration dynamics resulting from ECM-mediated mechanical communications. In particular, the APPA model generalizes the classic active-Brownian-particle (ABP) model by imposing a pairwise polarized attractive force between the particles, which depends on the instantaneous dynamic states of the particles and mimics the effective mutual pulling between the cells via the fiber bundle bridge. The APPA system exhibits enhanced aggregation behaviors compared to the classic ABP system, and the contrast is more apparent at lower particle densities and higher rotational diffusivities. Importantly, in contrast to the classic ABP system where the particle velocities are not correlated for all particle densities, the high-density phase of the APPA system exhibits strong dynamic correlations, which are characterized by the slowly decaying velocity correlation functions with a correlation length comparable to the linear size of the high-density phase domain (i.e., the cluster of particles). The strongly correlated multicellular dynamics predicted by the APPA model is subsequently verified in in vitro experiments using MCF-10A cells. Our studies indicate the importance of incorporating ECM-mediated mechanical coupling among the migrating cells for appropriately modeling emergent multicellular dynamics in complex microenvironments.

Towards Design of Self-Organizing Biomimetic Systems

The ability of designing biosynthetic systems with well-defined functional biomodules from scratch is an ambitious and revolutionary goal to deliver innovative, engineered solutions to future challenges in biotechnology and process systems engineering. In this work, several key challenges including modularization, functional biomodule identification, and assembly are discussed. In addition, an in silico protocell modeling approach is presented as a foundation for a computational model-based toolkit for rational analysis and modular design of biomimetic systems.

Modeling cell migration regulated by cell extracellular-matrix micromechanical coupling

Cell migration in fibrous extracellular matrix (ECM) is crucial to many physiological and pathological processes such as tissue regeneration, immune response, and cancer progression. During migration, individual cells can generate active pulling forces via actomyosin contraction, which are transmitted to the ECM fibers through focal adhesion complexes, remodel the ECM, and eventually propagate to and can be sensed by other cells in the system. The microstructure and physical properties of the ECM can also significantly influence cell migration, e.g., via durotaxis and contact guidance. Here, we develop a computational model for two-dimensional cell migration regulated by cell-ECM micromechanical coupling. Our model explicitly takes into account a variety of cellular-level processes, including focal adhesion formation and disassembly, active traction force generation and cell locomotion due to actin filament contraction, transmission and propagation of tensile forces in the ECM, as well as the resulting ECM remodeling. We validate our model by accurately reproducing single-cell dynamics of MCF-10A breast cancer cells migrating on collagen gels and show that the durotaxis and contact guidance effects naturally arise as a consequence of the cell-ECM micromechanical interactions considered in the model. Moreover, our model predicts strongly correlated multicellular migration dynamics, which are resulted from the ECM-mediated mechanical coupling among the migrating cell and are subsequently verified in in vitro experiments using MCF-10A cells. Our computational model provides a robust tool to investigate emergent collective dynamics of multicellular systems in complex in vivo microenvironment and can be utilized to design in vitro microenvironments to guide collective behaviors and self-organization of cells.

A DISCRETE MATHEMATICAL MODEL FOR SINGLE AND COLLECTIVE MOVEMENT IN AMOEBOID CELLS

In this paper, we develop a new discrete mathematical model for individual and collective cell motility. We introduce a mechanical model for the movement of a cell on a two-dimensional rigid surface to describe and investigate the cell–cell and cell–substrate interactions. The cell cytoskeleton is modeled as a series of springs and dashpots connected in parallel. The cell–substrate attachments and the cell protrusions are also included. In particular, this model is used to describe the directed movement of endothelial cells on a Matrigel plate. We compare the results from our model with experimental data. We show that cell density and substrate rigidity play an important role in network formation.

Comparison of cell migration mechanical strategies in three-dimensional matrices: a computational study

Cell migration on a two-dimensional flat surface has been extensively studied and is generally characterized by a front-protrusion-rear-contraction process. In a three-dimensional (3D) environment, on the other hand, cells adopt multiple migration strategies depending on the cell type and the properties of the extracellular matrix (ECM). By using computer simulations, we find that these migration strategies can be classified by various spatial-temporal dynamics of actin protrusion, actin-myosin contraction and actin-ECM adhesion. We demonstrate that if we include or exclude proteolysis of ECM, and vary adhesion dynamics and spatial distributions of protrusion, contraction and adhesion, our model can reproduce six experimentally observed motility modes: mesenchymal, chimneying, amoeboid, blebbing, finger-like protrusion and rear-squeezing cell locomotory behaviours. We further find that the mode of the cell motility evolves in response to the ECM density and adhesion detachment rate. The model makes non-trivial predictions about cell speed as a function of the adhesion strength, and ECM elasticity and mesh size.

A computational method for the coupled solution of reaction-diffusion equations on evolving domains and manifolds: Application to a model of cell migration and chemotaxis

In this paper, we devise a moving mesh finite element method for the approximate solution of coupled bulk-surface reaction-diffusion equations on an evolving two dimensional domain. Fundamental to the success of the method is the robust generation of bulk and surface meshes. For this purpose, we use a novel moving mesh partial differential equation (MMPDE) approach. The developed method is applied to model problems with known analytical solutions; these experiments indicate second-order spatial and temporal accuracy. Coupled bulk-surface problems occur frequently in many areas; in particular, in the modelling of eukaryotic cell migration and chemotaxis. We apply the method to a model of the two-way interaction of a migrating cell in a chemotactic field, where the bulk region corresponds to the extracellular region and the surface to the cell membrane.

Assembly and positioning of actomyosin rings by contractility and planar cell polarity

The actomyosin cytoskeleton is a primary force-generating mechanism in morphogenesis, thus a robust spatial control of cytoskeletal positioning is essential. In this report, we demonstrate that actomyosin contractility and planar cell polarity (PCP) interact in post-mitotic Ciona notochord cells to self-assemble and reposition actomyosin rings, which play an essential role for cell elongation. Intriguingly, rings always form at the cells' anterior edge before migrating towards the center as contractility increases, reflecting a novel dynamical property of the cortex. Our drug and genetic manipulations uncover a tug-of-war between contractility, which localizes cortical flows toward the equator and PCP, which tries to reposition them. We develop a simple model of the physical forces underlying this tug-of-war, which quantitatively reproduces our results. We thus propose a quantitative framework for dissecting the relative contribution of contractility and PCP to the self-assembly and repositioning of cytoskeletal structures, which should be applicable to other morphogenetic events.

Two-Phase Acto-Cytosolic Fluid Flow in a Moving Keratocyte: A 2D Continuum Model

The F-actin network and cytosol in the lamellipodia of crawling cells flow in a centripetal pattern and spout-like form, respectively. We have numerically studied this two-phase flow in the realistic geometry of a moving keratocyte. Cytosol has been treated as a low viscosity Newtonian fluid flowing through the high viscosity porous medium of F-actin network. Other involved phenomena including myosin activity, adhesion friction, and interphase interaction are also discussed to provide an overall view of this problem. Adopting a two-phase coupled model by myosin concentration, we have found new accurate perspectives of acto-cytosolic flow and pressure fields, myosin distribution, as well as the distribution of effective forces across the lamellipodia of a keratocyte with stationary shape. The order of magnitude method is also used to determine the contribution of forces in the internal dynamics of lamellipodia.

Stability analysis and simulations of coupled bulk-surface reaction-diffusion systems

In this article, we formulate new models for coupled systems of bulk-surface reaction-diffusion equations on stationary volumes. The bulk reaction-diffusion equations are coupled to the surface reaction-diffusion equations through linear Robin-type boundary conditions. We then state and prove the necessary conditions for diffusion-driven instability for the coupled system. Owing to the nature of the coupling between bulk and surface dynamics, we are able to decouple the stability analysis of the bulk and surface dynamics. Under a suitable choice of model parameter values, the bulk reaction-diffusion system can induce patterning on the surface independent of whether the surface reaction-diffusion system produces or not, patterning. On the other hand, the surface reaction-diffusion system cannot generate patterns everywhere in the bulk in the absence of patterning from the bulk reaction-diffusion system. For this case, patterns can be induced only in regions close to the surface membrane. Various numerical experiments are presented to support our theoretical findings. Our most revealing numerical result is that, Robin-type boundary conditions seem to introduce a boundary layer coupling the bulk and surface dynamics.

The spatial-temporal characteristics of type I collagen-based extracellular matrix

Type I collagen abounds in mammalian extracellular matrix (ECM) and is crucial to many biophysical processes. While previous studies have mostly focused on bulk averaged properties, here we provide a comprehensive and quantitative spatial-temporal characterization of the microstructure of type I collagen-based ECM as the gelation temperature varies. The structural characteristics including the density and nematic correlation functions are obtained by analyzing confocal images of collagen gels prepared at a wide range of gelation temperatures (from 16 °C to 36 °C). As temperature increases, the gel microstructure varies from a "bundled" network with strong orientational correlation between the fibers to an isotropic homogeneous network with no significant orientational correlation, as manifested by the decaying of length scales in the correlation functions. We develop a kinetic Monte-Carlo collagen growth model to better understand how ECM microstructure depends on various environmental or kinetic factors. We show that the nucleation rate, growth rate, and an effective hydrodynamic alignment of collagen fibers fully determines the spatiotemporal fluctuations of the density and orientational order of collagen gel microstructure. Also the temperature dependence of the growth rate and nucleation rate follow the prediction of classical nucleation theory.

Analytical solutions of actin-retrograde-flow in a circular stationary cell: a mechanical point of view

The network of actin filaments in the lamellipodium (LP) of stationary and migrating cells flows in a retrograde direction, from the membrane periphery toward the cell nucleus. We have theoretically studied this phenomenon in the circular stationary (fully spread) cells. Adopting a continuum view on the LP actin network, new closed-form solutions are provided for the actin-retrograde-flow (ARF) in a polar coordinate system. Due to discrepancy in the mechanical models of the actin network in the ARF regime, solutions are provided for both assumptions of solid and fluid behavior. Other involved phenomena, including polymerizing machine at the membrane periphery, cytosol drag, adhesion friction, and membrane tension, are also discussed to provide an overall quantitative view on this problem.

Asymmetry between pushing and pulling for crawling cells

Eukaryotic cells possess motility mechanisms allowing them not only to self-propel but also to exert forces on obstacles (to push) and to carry cargoes (to pull). To study the inherent asymmetry between active pushing and pulling we model a crawling acto-myosin cell extract as a one-dimensional layer of active gel subjected to external forces. We show that pushing is controlled by protrusion and that the macroscopic signature of the protrusion dominated motility mechanism is concavity of the force-velocity relation. In contrast, pulling is driven by protrusion only at small values of the pulling force and it is replaced by contraction when the pulling force is sufficiently large. This leads to more complex convex-concave structure of the force-velocity relation; in particular, competition between protrusion and contraction can produce negative mobility in a biologically relevant range. The model illustrates active readjustment of the force generating machinery in response to changes in the dipole structure of external forces. The possibility of switching between complementary active mechanisms implies that if necessary "pushers" can replace "pullers" and vice versa.