"Self-assisted" amoeboid navigation in complex environments

Optimization-based synthesis with directed cell migration

Collective behavior of biological agents from cells to herds of organisms is a fundamental feature in systems biology and in the emergence of new phenomena in the biological environment. Collective cell migration (CCM) under a physical or chemical cue is an example of this fundamental phenomenon where the individual migration of a cell is driven by the collective behavior of the neighboring cells and vice versa. The goal of this research is to discover the mathematical rules of collective cell migration with dynamic mode decomposition (DMD) with the use of experimental data and to test the predictive nature of the models with independent experimental data sets subject to Dirichlet, Neumann, and mixed boundary conditions. Both single and multi-cellular systems are investigated in this process. Additionally, the goal of this research is to create an optimal trajectory for microscopic robots in the presence of an obstacle course made of both static and dynamic obstacles. Such an optimization is made possible by synthesizing the discovered dynamics for cell migration with a numerical approach to dynamic optimization known as collocation by augmenting the discovered dynamics to the constraint equations. The optimal trajectory results presented in silico have potential design applications for the path planning of microrobots for therapeutic purposes such as cancer cell drug delivery, microsurgery, microsensing for early disease detection, and cleaning of toxic substances.

Mechanics of 3D Cell-Hydrogel Interactions: Experiments, Models, and Mechanisms

Hydrogels are highly water-swollen molecular networks that are ideal platforms to create tissue mimetics owing to their vast and tunable properties. As such, hydrogels are promising cell-delivery vehicles for applications in tissue engineering and have also emerged as an important base for ex vivo models to study healthy and pathophysiological events in a carefully controlled three-dimensional environment. Cells are readily encapsulated in hydrogels resulting in a plethora of biochemical and mechanical communication mechanisms, which recapitulates the natural cell and extracellular matrix interaction in tissues. These interactions are complex, with multiple events that are invariably coupled and spanning multiple length and time scales. To study and identify the underlying mechanisms involved, an integrated experimental and computational approach is ideally needed. This review discusses the state of our knowledge on cell-hydrogel interactions, with a focus on mechanics and transport, and in this context, highlights recent advancements in experiments, mathematical and computational modeling. The review begins with a background on the thermodynamics and physics fundamentals that govern hydrogel mechanics and transport. The review focuses on two main classes of hydrogels, described as semiflexible polymer networks that represent physically cross-linked fibrous hydrogels and flexible polymer networks representing the chemically cross-linked synthetic and natural hydrogels. In this review, we highlight five main cell-hydrogel interactions that involve key cellular functions related to communication, mechanosensing, migration, growth, and tissue deposition and elaboration. For each of these cellular functions, recent experiments and the most up to date modeling strategies are discussed and then followed by a summary of how to tune hydrogel properties to achieve a desired functional cellular outcome. We conclude with a summary linking these advancements and make the case for the need to integrate experiments and modeling to advance our fundamental understanding of cell-matrix interactions that will ultimately help identify new therapeutic approaches and enable successful tissue engineering.

Spatial heterogeneities shape the collective behavior of signaling amoeboid cells

In its natural habitat in the forest soil, the cellular slime mold Dictyostelium discoideum is exposed to obstacles. Starving Dictyostelium cells secrete cAMP, which is the key extracellular signaling molecule that promotes the aggregation process required for their long-term survival. Here, we investigated the influence of environmental inhomogeneities on the signaling and pattern formation of Dictyostelium cells. We present experimental data and numerical simulations on the pattern formation of signaling Dictyostelium cells in the presence of periodic arrays of millimeter-sized pillars. We observed concentric cAMP waves that initiated almost synchronously at the pillars and propagated outward. In response to these circular waves, the Dictyostelium cells streamed toward the pillars, forming aggregates arranged in patterns that reflected the periodicity of the lattice of pillars. Our results suggest that, in nature, the excitability threshold and synchronization level of the cells are two key parameters that control the nature of the interaction between cells and spatial heterogeneities in their environment.

A computational study of amoeboid motility in 3D: the role of extracellular matrix geometry, cell deformability, and cell-matrix adhesion

Amoeboid cells often migrate using pseudopods, which are membrane protrusions that grow, bifurcate, and retract dynamically, resulting in a net cell displacement. Many cells within the human body, such as immune cells, epithelial cells, and even metastatic cancer cells, can migrate using the amoeboid phenotype. Amoeboid motility is a complex and multiscale process, where cell deformation, biochemistry, and cytosolic and extracellular fluid motions are coupled. Furthermore, the extracellular matrix (ECM) provides a confined, complex, and heterogeneous environment for the cells to navigate through. Amoeboid cells can migrate without significantly remodeling the ECM using weak or no adhesion, instead utilizing their deformability and the microstructure of the ECM to gain enough traction. While a large volume of work exists on cell motility on 2D substrates, amoeboid motility is 3D in nature. Despite recent progress in modeling cellular motility in 3D, there is a lack of systematic evaluations of the role of ECM microstructure, cell deformability, and adhesion on 3D motility. To fill this knowledge gap, here we present a multiscale, multiphysics modeling study of amoeboid motility through 3D-idealized ECM. The model is a coupled fluid‒structure and coarse-grain biochemistry interaction model that accounts for large deformation of cells, pseudopod dynamics, cytoplasmic and extracellular fluid motion, stochastic dynamics of cell-ECM adhesion, and microstructural (pore-scale) geometric details of the ECM. The key finding of the study is that cell deformation and matrix porosity strongly influence amoeboid motility, while weak adhesion and microscale structural details of the ECM have secondary but subtle effects.

A computational model of amoeboid cell motility in the presence of obstacles

Locomotion of amoeboid cells is mediated by finger-like protrusions of the cell body, known as pseudopods, which grow, bifurcate, and retract in a dynamic fashion. Pseudopods are the primary mode of locomotion for many cells within the human body, such as leukocytes, embryonic cells, and metastatic cancer cells. Amoeboid motility is a complex and multiscale process, which involves bio-molecular reactions, cell deformation, and cytoplasmic and extracellular fluid motion. Additionally, cells within the human body are subject to a confined 3D environment known as the extra-cellular matrix (ECM), which resembles a fluid-filled porous medium. In this article, we present a 3D, multiphysics computational approach coupling fluid mechanics, solid mechanics, and a pattern formation model to simulate locomotion of amoeboid cells through a porous matrix composed of a viscous fluid and an array of finite-sized spherical obstacles. The model combines reaction-diffusion of activator/inhibitors, extreme deformation of the cell, pseudopod dynamics, cytoplasmic and extracellular fluid motion, and fully resolved extracellular matrix. A surface finite-element method is used to obtain the cell deformation and activator/inhibitor concentrations, while the fluid motion is solved using a combined finite-volume and spectral method. The immersed-boundary methods are used to couple the cell deformation, obstacles, and fluid. The model is able to recreate squeezing and weaving motion of cells through the matrix. We study the influence of matrix porosity, obstacle size, and cell deformability on the motility behavior. It is found that below certain values of these parameters, cell motion is completely inhibited. Phase diagrams are presented depicting such motility limits. Interesting dynamics seen in the presence of obstacles but absent in unconfined medium, such as freezing or cell arrest, probing, doubling-back, and tug-of-war are predicted. Furthermore, persistent unidirectional motion of cells that is often observed in an unconfined medium is shown to be lost in presence of obstacles, and is attributed to an alteration of the pseudopod dynamics. The same mechanism, however, allows the cell to find a new direction to penetrate further into the matrix without being stuck in one place. The results and analysis presented here show a strong coupling between cell deformability and ECM properties, and provide new fluid mechanical insights on amoeboid motility in confined medium.

Amoeboid-mesenchymal migration plasticity promotes invasion only in complex heterogeneous microenvironments

During tissue invasion individual tumor cells exhibit two interconvertible migration modes, namely mesenchymal and amoeboid migration. The cellular microenvironment triggers the switch between both modes, thereby allowing adaptation to dynamic conditions. It is, however, unclear if this amoeboid-mesenchymal migration plasticity contributes to a more effective tumor invasion. We address this question with a mathematical model, where the amoeboid-mesenchymal migration plasticity is regulated in response to local extracellular matrix resistance. Our numerical analysis reveals that extracellular matrix structure and presence of a chemotactic gradient are key determinants of the model behavior. Only in complex microenvironments, if the extracellular matrix is highly heterogeneous and a chemotactic gradient directs migration, the amoeboid-mesenchymal migration plasticity allows a more widespread invasion compared to the non-switching amoeboid and mesenchymal modes. Importantly, these specific conditions are characteristic for in vivo tumor invasion. Thus, our study suggests that in vitro systems aiming at unraveling the underlying molecular mechanisms of tumor invasion should take into account the complexity of the microenvironment by considering the combined effects of structural heterogeneities and chemical gradients on cell migration.

Mechanisms of Cell Polarization

Cell polarization is a key step in the migration, development, and organization of eukaryotic cells, both at the single cell and multicellular level. Research on the mechanisms that give rise to polarization of a given cell, and organization of polarity within a tissue has led to new understanding across cellular and developmental biology. In this review, we describe some of the history of theoretical and experimental aspects of the field, as well as some interesting questions and challenges for the future.

Computational model for amoeboid motion: Coupling membrane and cytosol dynamics

A distinguishing feature of amoeboid motion is that the migrating cell undergoes large deformations, caused by the emergence and retraction of actin-rich protrusions, called pseudopods. Here, we propose a cell motility model that represents pseudopod dynamics, as well as its interaction with membrane signaling molecules. The model accounts for internal and external forces, such as protrusion, contraction, adhesion, surface tension, or those arising from cell-obstacle contacts. By coupling the membrane and cytosol interactions we are able to reproduce a realistic picture of amoeboid motion. The model results are in quantitative agreement with experiments and show how cells may take advantage of the geometry of their microenvironment to migrate more efficiently.

Challenges in tissue engineering - towards cell control inside artificial scaffolds

Control of living cells is vital for the survival of organisms. Each cell inside an organism is exposed to diverse external mechano-chemical cues, all coordinated in a spatio-temporal pattern triggering individual cell functions. This complex interplay between external chemical cues and mechanical 3D environments is translated into intracellular signaling loops. Here, we describe how external mechano-chemical cues control cell functions, especially cell migration, and influence intracellular information transport. In particular, this work focuses on the quantitative analysis of (1) intracellular vesicle transport to understand intracellular state changes in response to external cues, (2) cellular sensing of external chemotactic cues, and (3) the cells' ability to migrate in 3D structured environments, artificially fabricated to mimic the 3D environment of tissue in the human body.

The motility-proliferation-metabolism interplay during metastatic invasion

Metastasis is the major cause for cancer patients' death, and despite all the recent advances in cancer research it is still mostly incurable. Understanding the mechanisms that are involved in the migration of the cells in a complex environment is a key step towards successful anti-metastatic treatment. Using experimental data-based modeling, we focus on the fundamentals of metastatic invasion: motility, invasion, proliferation and metabolism, and study how they may be combined to maximize the cancer's ability to metastasize. The modeled cells' performance is measured by the number of cells that succeed in migration in a maze, which mimics the extracellular environment. We show that co-existence of different cell clones in the tumor, as often found in experiments, optimizes the invasive ability in a frequently-changing environment. We study the role of metabolism and stimulation by growth factors, and show that metabolism plays a crucial role in the metastatic process and should therefore be targeted for successful treatment.

Tumor invasion optimization by mesenchymal-amoeboid heterogeneity

Metastasizing tumor cells migrate through the surrounding tissue and extracellular matrix toward the blood vessels, in order to colonize distant organs. They typically move in a dense environment, filled with other cells. In this work we study cooperative effects between neighboring cells of different types, migrating in a maze-like environment with directional cue. Using a computerized model, we measure the percentage of cells that arrive to the defined target, for different mesenchymal/amoeboid ratios. Wall degradation of mesenchymal cells, as well as motility of both types of cells, are coupled to metabolic energy-like resource level. We find that indirect cooperation emerges in mid-level energy, as mesenchymal cells create paths that are used by amoeboids. Therefore, we expect to see a small population of mesenchymals kept in a mostly-amoeboid population. We also study different forms of direct interaction between the cells, and show that energy-dependent interaction strength is optimal for the migration of both mesenchymals and amoeboids. The obtained characteristics of cellular cluster size are in agreement with experimental results. We therefore predict that hybrid states, e.g. epithelial-mesenchymal, should be utilized as a stress-response mechanism.

Excitable behavior in amoeboid chemotaxis

Chemotaxis, the directed motion of cells in response to chemical gradients, is a fundamental process. Eukaryotic cells detect spatial differences in chemoattractant receptor occupancy with high precision and use these differences to bias the location of actin-rich protrusions to guide their movement. Research into chemotaxis has benefitted greatly from a systems biology approach that combines novel experimental and computational tools to pose and test hypotheses. Recently, one such hypothesis has been postulated proposing that chemotaxis in eukaryotic cells is mediated by locally biasing the activity of an underlying excitable system. The excitable system hypothesis can account for a number of cellular behaviors related to chemotaxis, including the stochastic nature of the movement of unstimulated cells, the directional bias imposed by chemoattractant gradients, and the observed spatial and temporal distribution of signaling and cytoskeleton proteins.

A comparison of computational models for eukaryotic cell shape and motility

Eukaryotic cell motility involves complex interactions of signalling molecules, cytoskeleton, cell membrane, and mechanics interacting in space and time. Collectively, these components are used by the cell to interpret and respond to external stimuli, leading to polarization, protrusion, adhesion formation, and myosin-facilitated retraction. When these processes are choreographed correctly, shape change and motility results. A wealth of experimental data have identified numerous molecular constituents involved in these processes, but the complexity of their interactions and spatial organization make this a challenging problem to understand. This has motivated theoretical and computational approaches with simplified caricatures of cell structure and behaviour, each aiming to gain better understanding of certain kinds of cells and/or repertoire of behaviour. Reaction–diffusion (RD) equations as well as equations of viscoelastic flows have been used to describe the motility machinery. In this review, we describe some of the recent computational models for cell motility, concentrating on simulations of cell shape changes (mainly in two but also three dimensions). The problem is challenging not only due to the difficulty of abstracting and simplifying biological complexity but also because computing RD or fluid flow equations in deforming regions, known as a “free-boundary” problem, is an extremely challenging problem in applied mathematics. Here we describe the distinct approaches, comparing their strengths and weaknesses, and the kinds of biological questions that they have been able to address.

Traveling waves in actin dynamics and cell motility

Much of current understanding of cell motility arose from studying steady treadmilling of actin arrays. Recently, there have been a growing number of observations of a more complex, non-steady, actin behavior, including self-organized waves. It is becoming clear that these waves result from activation and inhibition feedbacks in actin dynamics acting on different scales, but the exact molecular nature of these feedbacks and the respective roles of biomechanics and biochemistry are still unclear. Here, we review recent advances achieved in experimental and theoretical studies of actin waves and discuss mechanisms and physiological significance of wavy protrusions.

Bacterial survival strategies suggest rethinking cancer cooperativity

Despite decades of a much improved understanding of cancer biology, we are still baffled by questions regarding the deadliest traits of malignancy: metastatic colonization, dormancy and relapse, and the rapid evolution of multiple drug and immune resistance. New ideas are needed to resolve these critical issues. Relying on finding and demonstrating parallels between collective behavior capabilities of cancer cells and that of bacteria, we suggest communal behaviors of bacteria as a valuable model system for new perspectives and research directions. Understanding the ways in which bacteria thrive in competitive habitats and their cooperative strategies for surviving extreme stress can shed light on cooperativity in tumorigenesis and portray tumors as societies of smart communicating cells. This may translate into progress in fathoming cancer pathogenesis. We outline new experiments to test the cancer cooperativity hypothesis and reason that cancer may be outsmarted through its own 'social intelligence'.

Modelling cell motility and chemotaxis with evolving surface finite elements

We present a mathematical and a computational framework for the modelling of cell motility. The cell membrane is represented by an evolving surface, with the movement of the cell determined by the interaction of various forces that act normal to the surface. We consider external forces such as those that may arise owing to inhomogeneities in the medium and a pressure that constrains the enclosed volume, as well as internal forces that arise from the reaction of the cells' surface to stretching and bending. We also consider a protrusive force associated with a reaction-diffusion system (RDS) posed on the cell membrane, with cell polarization modelled by this surface RDS. The computational method is based on an evolving surface finite-element method. The general method can account for the large deformations that arise in cell motility and allows the simulation of cell migration in three dimensions. We illustrate applications of the proposed modelling framework and numerical method by reporting on numerical simulations of a model for eukaryotic chemotaxis and a model for the persistent movement of keratocytes in two and three space dimensions. Movies of the simulated cells can be obtained from http://homepages.warwick.ac.uk/∼maskae/CV_Warwick/Chemotaxis.html.

Non-Brownian dynamics and strategy of amoeboid cell locomotion

Amoeboid cells such as Dictyostelium discoideum and Madin-Darby canine kidney cells show the non-Brownian dynamics of migration characterized by the superdiffusive increase of mean-squared displacement. In order to elucidate the physical mechanism of this non-Brownian dynamics, a computational model is developed which highlights a group of inhibitory molecules for actin polymerization. Based on this model, we propose a hypothesis that inhibitory molecules are sent backward in the moving cell to accumulate at the rear of cell. The accumulated inhibitory molecules at the rear further promote cell locomotion to form a slow positive feedback loop of the whole-cell scale. The persistent straightforward migration is stabilized with this feedback mechanism, but the fluctuation in the distribution of inhibitory molecules and the cell shape deformation concurrently interrupt the persistent motion to turn the cell into a new direction. A sequence of switching behaviors between persistent motions and random turns gives rise to the superdiffusive migration in the absence of the external guidance signal. In the complex environment with obstacles, this combined process of persistent motions and random turns drives the simulated amoebae to solve the maze problem in a highly efficient way, which suggests the biological advantage for cells to bear the non-Brownian dynamics.