Integration of in vitro and in silico Models Using Bayesian Optimization With an Application to Stochastic Modeling of Mesenchymal 3D Cell Migration

Three-component contour dynamics model to simulate and analyze amoeboid cell motility in two dimensions

Amoeboid cell motility is relevant in a wide variety of biomedical processes such as wound healing, cancer metastasis, and embryonic morphogenesis. It is characterized by pronounced changes of the cell shape associated with expansions and retractions of the cell membrane, which result in a crawling kind of locomotion. Despite existing computational models of amoeboid motion, the inference of expansion and retraction components of individual cells, the corresponding classification of cells, and the a priori specification of the parameter regime to achieve a specific motility behavior remain challenging open problems. We propose a novel model of the spatio-temporal evolution of two-dimensional cell contours comprising three biophysiologically motivated components: a stochastic term accounting for membrane protrusions and two deterministic terms accounting for membrane retractions by regularizing the shape and area of the contour. Mathematically, these correspond to the intensity of a self-exciting Poisson point process, the area-preserving curve-shortening flow, and an area adjustment flow. The model is used to generate contour data for a variety of qualitatively different, e.g., polarized and non-polarized, cell tracks that visually resemble experimental data very closely. In application to experimental cell tracks, we inferred the protrusion component and examined its correlation to common biomarkers: the F-actin density close to the membrane and its local motion. Due to the low model complexity, parameter estimation is fast, straightforward, and offers a simple way to classify contour dynamics based on two locomotion types: the amoeboid and a so-called fan-shaped type. For both types, we use cell tracks segmented from fluorescence imaging data of the model organism Dictyostelium discoideum. An implementation of the model is provided within the open-source software package AmoePy, a Python-based toolbox for analyzing and simulating amoeboid cell motility.

Subcellular mechano-regulation of cell migration in confined extracellular microenvironment

Cell migration is a highly coordinated cellular event that determines diverse physiological and pathological processes in which the continuous interaction of a migrating cell with neighboring cells or the extracellular matrix is regulated by the physical setting of the extracellular microenvironment. In confined spaces, cell migration occurs differently compared to unconfined open spaces owing to the additional forces that limit cell motility, which create a driving bias for cells to invade the confined space, resulting in a distinct cell motility process compared to what is expected in open spaces. Moreover, cells in confined environments can be subjected to elevated mechanical compression, which causes physical stimuli and activates the damage repair cycle in the cell, including the DNA in the nucleus. Although cells have a self-restoring system to repair damage from the cell membrane to the genetic components of the nucleus, this process may result in genetic and/or epigenetic alterations that can increase the risk of the progression of diverse diseases, such as cancer and immune disorders. Furthermore, there has been a shift in the paradigm of bioengineering from the development of new biomaterials to controlling biophysical cues and fine-tuning cell behaviors to cure damaged/diseased tissues. The external physical cues perceived by cells are transduced along the mechanosensitive machinery, which is further channeled into the nucleus through subcellular molecular linkages of the nucleoskeleton and cytoskeleton or the biochemical translocation of transcription factors. Thus, external cues can directly or indirectly regulate genetic transcriptional processes and nuclear mechanics, ultimately determining cell fate. In this review, we discuss the importance of the biophysical cues, response mechanisms, and mechanical models of cell migration in confined environments. We also discuss the effect of force-dependent deformation of subcellular components, specifically focusing on subnuclear organelles, such as nuclear membranes and chromosomal organization. This review will provide a biophysical perspective on cancer progression and metastasis as well as abnormal cellular proliferation.

Mathematical model of mechano-sensing and mechanically induced collective motility of cells on planar elastic substrates

Cells mechanically interact with their environment to sense, for example, topography, elasticity and mechanical cues from other cells. Mechano-sensing has profound effects on cellular behaviour, including motility. The current study aims to develop a mathematical model of cellular mechano-sensing on planar elastic substrates and demonstrate the model's predictive capabilities for the motility of individual cells in a colony. In the model, a cell is assumed to transmit an adhesion force, derived from a dynamic focal adhesion integrin density, that locally deforms a substrate, and to sense substrate deformation originating from neighbouring cells. The substrate deformation from multiple cells is expressed as total strain energy density with a spatially varying gradient. The magnitude and direction of the gradient at the cell location define the cell motion. Cell-substrate friction, partial motion randomness, and cell death and division are included. The substrate deformation by a single cell and the motility of two cells are presented for several substrate elasticities and thicknesses. The collective motility of 25 cells on a uniform substrate mimicking the closure of a circular wound of 200 µm is predicted for deterministic and random motion. Cell motility on substrates with varying elasticity and thickness is explored for four cells and 15 cells, the latter again mimicking wound closure. Wound closure by 45 cells is used to demonstrate the simulation of cell death and division during migration. The mathematical model can adequately simulate the mechanically induced collective cell motility on planar elastic substrates. The model is suitable for extension to other cell and substrates shapes and the inclusion of chemotactic cues, offering the potential to complement in vitro and in vivo studies.

3D collagen migration patterns reveal a SMAD3-dependent and TGF-β1-independent mechanism of recruitment for tumour-associated fibroblasts in lung adenocarcinoma

BACKGROUND

The TGF-β1 transcription factor SMAD3 is epigenetically repressed in tumour-associated fibroblasts (TAFs) from lung squamous cell carcinoma (SCC) but not adenocarcinoma (ADC) patients, which elicits a compensatory increase in SMAD2 that renders SCC-TAFs less fibrotic. Here we examined the effects of altered SMAD2/3 in fibroblast migration and its impact on the desmoplastic stroma formation in lung cancer.

METHODS

We used a microfluidic device to examine descriptors of early protrusions and subsequent migration in 3D collagen gels upon knocking down SMAD2 or SMAD3 by shRNA in control fibroblasts and TAFs.

RESULTS

High SMAD3 conditions as in shSMAD2 fibroblasts and ADC-TAFs exhibited a migratory advantage in terms of protrusions (fewer and longer) and migration (faster and more directional) selectively without TGF-β1 along with Erk1/2 hyperactivation. This enhanced migration was abrogated by TGF-β1 as well as low glucose medium and the MEK inhibitor Trametinib. In contrast, high SMAD2 fibroblasts were poorly responsive to TGF-β1, high glucose and Trametinib, exhibiting impaired migration in all conditions.

CONCLUSIONS

The basal migration advantage of high SMAD3 fibroblasts provides a straightforward mechanism underlying the larger accumulation of TAFs previously reported in ADC compared to SCC. Moreover, our results encourage using MEK inhibitors in ADC-TAFs but not SCC-TAFs.

PhysiCOOL: A generalized framework for model Calibration and Optimization Of modeLing projects

In silico models of biological systems are usually very complex and rely on a large number of parameters describing physical and biological properties that require validation. As such, parameter space exploration is an essential component of computational model development to fully characterize and validate simulation results. Experimental data may also be used to constrain parameter space (or enable model calibration) to enhance the biological relevance of model parameters. One widely used computational platform in the mathematical biology community is PhysiCell, which provides a standardized approach to agent-based models of biological phenomena at different time and spatial scales. Nonetheless, one limitation of PhysiCell is the lack of a generalized approach for parameter space exploration and calibration that can be run without high-performance computing access. Here, we present PhysiCOOL, an open-source Python library tailored to create standardized calibration and optimization routines for PhysiCell models.

A novel integrated experimental and computational approach to unravel fibroblast motility in response to chemical gradients in 3D collagen matrices

Fibroblasts play an essential role in tissue repair and regeneration as they migrate to wounded areas to secrete and remodel the extracellular matrix. Fibroblasts recognize chemical substances such as growth factors, which enhance their motility towards the wounded tissues through chemotaxis. Although several studies have characterized single-cell fibroblast motility before, the migration patterns of fibroblasts in response to external factors have not been fully explored in 3D environments. We present a study that combines experimental and computational efforts to characterize the effect of chemical stimuli on the invasion of 3D collagen matrices by fibroblasts. Experimentally, we used microfluidic devices to create chemical gradients using collagen matrices of distinct densities. We evaluated how cell migration patterns were affected by the presence of growth factors and the mechanical properties of the matrix. Based on these results, we present a discrete-based computational model to simulate cell motility, which we calibrated through the quantitative comparison of experimental and computational data via Bayesian optimization. By combining these approaches, we predict that fibroblasts respond to both the presence of chemical factors and their spatial location. Furthermore, our results show that the presence of these chemical gradients could be reproduced by our computational model through increases in the magnitude of cell-generated forces and enhanced cell directionality. Although these model predictions require further experimental validation, we propose that our framework can be applied as a tool that takes advantage of experimental data to guide the calibration of models and predict which mechanisms at the cellular level may justify the experimental findings. Consequently, these new insights may also guide the design of new experiments, tailored to validate the variables of interest identified by the model.

Unravelling cell migration: defining movement from the cell surface

Cell motility is essential for life and development. Unfortunately, cell migration is also linked to several pathological processes, such as cancer metastasis. Cells' ability to migrate relies on many actors. Cells change their migratory strategy based on their phenotype and the properties of the surrounding microenvironment. Cell migration is, therefore, an extremely complex phenomenon. Researchers have investigated cell motility for more than a century. Recent discoveries have uncovered some of the mysteries associated with the mechanisms involved in cell migration, such as intracellular signaling and cell mechanics. These findings involve different players, including transmembrane receptors, adhesive complexes, cytoskeletal components , the nucleus, and the extracellular matrix. This review aims to give a global overview of our current understanding of cell migration.

A mechanistic protrusive-based model for 3D cell migration

Cell migration is essential for a variety of biological processes, such as embryogenesis, wound healing, and the immune response. After more than a century of research-mainly on flat surfaces-, there are still many unknowns about cell motility. In particular, regarding how cells migrate within 3D matrices, which more accurately replicate in vivo conditions. We present a novel in silico model of 3D mesenchymal cell migration regulated by the chemical and mechanical profile of the surrounding environment. This in silico model considers cell's adhesive and nuclear phenotypes, the effects of the steric hindrance of the matrix, and cells ability to degradate the ECM. These factors are crucial when investigating the increasing difficulty that migrating cells find to squeeze their nuclei through dense matrices, which may act as physical barriers. Our results agree with previous in vitro observations where fibroblasts cultured in collagen-based hydrogels did not durotax toward regions with higher collagen concentrations. Instead, they exhibited an adurotactic behavior, following a more random trajectory. Overall, cell's migratory response in 3D domains depends on its phenotype, and the properties of the surrounding environment, that is, 3D cell motion is strongly dependent on the context.

Confined Cell Migration and Asymmetric Hydraulic Environments to Evaluate The Metastatic Potential of Cancer Cells

Metastasis, a hallmark of cancer development, is also the leading reason for most cancer-related deaths. Furthermore, cancer cells are highly adaptable to microenvironments and can migrate along pre-existing channel-like tracks of anatomical structures. However, more representative three-dimensional models are required to reproduce the heterogeneity of metastatic cell migration in vivo to further understand the metastasis mechanism and develop novel therapeutic strategies against it. Here, we designed and fabricated different microfluidic-based devices that recreate confined migration and diverse environments with asymmetric hydraulic resistances. Our results show different migratory potential between metastatic and nonmetastatic cancer cells in confined environments. Moreover, although nonmetastatic cells have not been tested against barotaxis due to their low migration capacity, metastatic cells present an enhanced preference to migrate through the lowest resistance path, being sensitive to barotaxis. This device, approaching the study of metastasis capability based on confined cell migration and barotactic cell decisions, may pave the way for the implementation of such technology to determine and screen the metastatic potential of certain cancer cells.

Magnesium ions regulate mesenchymal stem cells population and osteogenic differentiation: A fuzzy agent-based modeling approach

Mesenchymal stem cells (MSCs) are proliferative and multipotent cells that play a key role in the bone regeneration process. Empirical data have repeatedly shown the bioregulatory importance of magnesium (Mg) ions in MSC growth and osteogenesis. In this study, we propose an agent-based model to predict the spatiotemporal dynamics of the MSC population and osteogenic differentiation in response to Mg2+ ions. A fuzzy-logic controller was designed to govern the decision-making process of cells by predicting four cellular processes of proliferation, differentiation, migration, and mortality in response to several important bioregulatory factors such as Mg2+ ions, pH, BMP2, and TGF-β1. The model was calibrated using the empirical data obtained from three sets of cell culture experiments. The model successfully reproduced the empirical observations regarding live cell count, viability, DNA content, and the differentiation-related markers of alkaline phosphate (ALP) and osteocalcin (OC). The simulation results, in agreement with the empirical data, showed that Mg2+ ions within 3-6 mM concentration have the highest stimulation effect on cell population growth. The model also correctly reproduced the stimulatory effect of Mg2+ ions on ALP and its inhibitory effect on OC as the early and late differentiation markers, respectively. Besides, the numerical simulation shed light on the innate cellular differences of the cells cultured in different experiments in terms of the proliferative capacity as well as sensitivity to Mg2+ ions. The proposed model can be adopted in the study of the osteogenesis around Mg-based implants where ions released due to degradation interact with local cells and regulate bone regeneration.

Matrix degradation regulates osteoblast protrusion dynamics and individual migration

Protrusions are one of the structures that cells use to sense their surrounding environment in a probing and exploratory manner as well as to communicate with other cells. In particular, osteoblasts embedded within a 3D matrix tend to originate a large number of protrusions compared to other type of cells. In this work, we study the role that mechanochemical properties of the extracellular matrix (ECM) play on the dynamics of these protrusions, namely, the regulation of the size and number of emanating structures. In addition, we also determine how the dynamics of the protrusions may lead the 3D movement of the osteoblasts. Significant differences were found in protrusion size and cell velocity, when degradation activity due to metalloproteases was blocked by means of an artificial broad-spectrum matrix metalloproteinase inhibitor, whereas stiffening of the matrix by introducing transglutaminase crosslinking, only induced slight changes in both protrusion size and cell velocity, suggesting that the ability of cells to create a path through the matrix is more critical than the matrix mechanical properties themselves. To confirm this, we developed a cell migration computational model in 3D including both the mechanical and chemical properties of the ECM as well as the protrusion mechanics, obtaining good agreement with experimental results.

The role of actin protrusion dynamics in cell migration through a degradable viscoelastic extracellular matrix: Insights from a computational model

Actin protrusion dynamics plays an important role in the regulation of three-dimensional (3D) cell migration. Cells form protrusions that adhere to the surrounding extracellular matrix (ECM), mechanically probe the ECM and contract in order to displace the cell body. This results in cell migration that can be directed by the mechanical anisotropy of the ECM. However, the subcellular processes that regulate protrusion dynamics in 3D cell migration are difficult to investigate experimentally and therefore not well understood. Here, we present a computational model of cell migration through a degradable viscoelastic ECM. This model is a 2D representation of 3D cell migration. The cell is modeled as an active deformable object that captures the viscoelastic behavior of the actin cortex and the subcellular processes underlying 3D cell migration. The ECM is regarded as a viscoelastic material, with or without anisotropy due to fibrillar strain stiffening, and modeled by means of the meshless Lagrangian smoothed particle hydrodynamics (SPH) method. ECM degradation is captured by local fluidization of the material and permits cell migration through the ECM. We demonstrate that changes in ECM stiffness and cell strength affect cell migration and are accompanied by changes in number, lifetime and length of protrusions. Interestingly, directly changing the total protrusion number or the average lifetime or length of protrusions does not affect cell migration. A stochastic variability in protrusion lifetime proves to be enough to explain differences in cell migration velocity. Force-dependent adhesion disassembly does not result in faster migration, but can make migration more efficient. We also demonstrate that when a number of simultaneous protrusions is enforced, the optimal number of simultaneous protrusions is one or two, depending on ECM anisotropy. Together, the model provides non-trivial new insights in the role of protrusions in 3D cell migration and can be a valuable contribution to increase the understanding of 3D cell migration mechanics.

Modelling actin polymerization: the effect on confined cell migration

The aim of this work is to model cell motility under conditions of mechanical confinement. This cell migration mode may occur in extravasation of tumour and neutrophil-like cells. Cell migration is the result of the complex action of different forces exerted by the interplay between myosin contractility forces and actin processes. Here, we propose and implement a finite element model of the confined migration of a single cell. In this model, we consider the effects of actin and myosin in cell motility. Both filament and globular actin are modelled. We model the cell considering cytoplasm and nucleus with different mechanical properties. The migration speed in the simulation is around 0.1 μm/min, which is in agreement with existing literature. From our simulation, we observe that the nucleus size has an important role in cell migration inside the channel. In the simulation the cell moves further when the nucleus is smaller. However, this speed is less sensitive to nucleus stiffness. The results show that the cell displacement is lower when the nucleus is stiffer. The degree of adhesion between the channel walls and the cell is also very important in confined migration. We observe an increment of cell velocity when the friction coefficient is higher.