A comparison of computational models for eukaryotic cell shape and motility

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.

Generative modeling of biological shapes and images using a probabilistic α-shape sampler

Understanding morphological variation is an important task in many areas of computational biology. Recent studies have focused on developing computational tools for the task of sub-image selection which aims at identifying structural features that best describe the variation between classes of shapes. A major part in assessing the utility of these approaches is to demonstrate their performance on both simulated and real datasets. However, when creating a model for shape statistics, real data can be difficult to access and the sample sizes for these data are often small due to them being expensive to collect. Meanwhile, the current landscape of generative models for shapes has been mostly limited to approaches that use black-box inference-making it difficult to systematically assess the power and calibration of sub-image models. In this paper, we introduce the α -shape sampler: a probabilistic framework for generating realistic 2D and 3D shapes based on probability distributions which can be learned from real data. We demonstrate our framework using proof-of-concept examples and in two real applications in biology where we generate (i) 2D images of healthy and septic neutrophils and (ii) 3D computed tomography (CT) scans of primate mandibular molars. The α -shape sampler R package is open-source and can be downloaded at https://github.com/lcrawlab/ashapesampler.

The Influence of Nucleus Mechanics in Modelling Adhesion-independent Cell Migration in Structured and Confined Environments

Recent biological experiments (Lämmermann et al. in Nature 453(7191):51-55, 2008; Reversat et al. in Nature 7813:582-585, 2020; Balzer et al. in ASEB J Off Publ Fed Am Soc Exp Biol 26(10):4045-4056, 2012) have shown that certain types of cells are able to move in structured and confined environments even without the activation of focal adhesion. Focusing on this particular phenomenon and based on previous works (Jankowiak et al. in Math Models Methods Appl Sci 30(03):513-537, 2020), we derive a novel two-dimensional mechanical model, which relies on the following physical ingredients: the asymmetrical renewal of the actin cortex supporting the membrane, resulting in a backward flow of material; the mechanical description of the nuclear membrane and the inner nuclear material; the microtubule network guiding nucleus location; the contact interactions between the cell and the external environment. The resulting fourth order system of partial differential equations is then solved numerically to conduct a study of the qualitative effects of the model parameters, mainly those governing the mechanical properties of the nucleus and the geometry of the confining structure. Coherently with biological observations, we find that cells characterized by a stiff nucleus are unable to migrate in channels that can be crossed by cells with a softer nucleus. Regarding the geometry, cell velocity and ability to migrate are influenced by the width of the channel and the wavelength of the external structure. Even though still preliminary, these results may be potentially useful in determining the physical limit of cell migration in confined environments and in designing scaffolds for tissue engineering.

Simulating 3D Cell Shape with the Cellular Potts Model

Computer simulations have become a widely used method for the field of mechanobiology. An important question is whether one can predict the shape and forces of cells as a function of the extracellular environment. Different types of models have been described before to simulate cell and tissue shapes in structured environments. In this chapter, we give a brief overview of commonly used models and then describe the Cellular Potts Model, a lattice-based modelling framework, in more detail. We provide a hands-on guide on how to build a model that simulates the shape of a single cell on a micropattern in three dimensions in different open source software packages using the Cellular Potts framework. A simulation is set up with an initial configuration of generalized cells that change shape and position due to an energy function that incorporates cellular volume and surface area constraints as well as interaction energies between the generalized cells.

Two Complementary Signaling Pathways Depict Eukaryotic Chemotaxis: A Mechanochemical Coupling Model

Many eukaryotic cells, including neutrophils and Dictyostelium cells, are able to undergo correlated random migration in the absence of directional cues while reacting to shallow gradients of chemoattractants with exquisite precision. Although progress has been made with regard to molecular identities, it remains elusive how molecular mechanics are integrated with cell mechanics to initiate and manipulate cell motility. Here, we propose a two dimensional (2D) cell migration model wherein a multilayered dynamic seesaw mechanism is accompanied by a mechanical strain-based inhibition mechanism. In biology, these two mechanisms can be mapped onto the biochemical feedback between phosphoinositides (PIs) and Rho GTPase and the mechanical interplay between filamin A (FLNa) and FilGAP. Cell migration and the accompanying morphological changes are demonstrated in numerical simulations using a particle-spring model, and the diffusion in the cell membrane are simulations using a one dimensional (1D) finite differences method (FDM). The fine balance established between endogenous signaling and a mechanically governed inactivation scheme ensures the endogenous cycle of self-organizing pseudopods, accounting for the correlated random migration. Furthermore, this model cell manifests directional and adaptable responses to shallow graded signaling, depending on the overwhelming effect of the graded stimuli guidance on strain-based inhibition. Finally, the model cell becomes trapped within an obstacle-ridden spatial region, manifesting a shuttle run for local explorations and can chemotactically “escape”, illustrating again the balance required in the complementary signaling pathways.

Biomechanics of cells and subcellular components: A comprehensive review of computational models and applications

Cells are a fundamental structural, functional and biological unit for all living organisms. Up till now, considerable efforts have been made to study the responses of single cells and subcellular components to an external load, and understand the biophysics underlying cell rheology, mechanotransduction and cell functions using experimental and in silico approaches. In the last decade, computational simulation has become increasingly attractive due to its critical role in interpreting experimental data, analysing complex cellular/subcellular structures, facilitating diagnostic designs and therapeutic techniques, and developing biomimetic materials. Despite the significant progress, developing comprehensive and accurate models of living cells remains a grand challenge in the 21st century. To understand current state of the art, this review summarises and classifies the vast array of computational biomechanical models for cells. The article covers the cellular components at multi-spatial levels, that is, protein polymers, subcellular components, whole cells and the systems with scale beyond a cell. In addition to the comprehensive review of the topic, this article also provides new insights into the future prospects of developing integrated, active and high-fidelity cell models that are multiscale, multi-physics and multi-disciplinary in nature. This review will be beneficial for the researchers in modelling the biomechanics of subcellular components, cells and multiple cell systems and understanding the cell functions and biological processes from the perspective of cell mechanics.

Mathematical modeling of chemotaxis guided amoeboid cell swimming

Cells and microorganisms adopt various strategies to migrate in response to different environmental stimuli. To date, many modeling research has focused on the crawling-basedDictyostelium discoideum(Dd) cells migration induced by chemotaxis, yet recent experimental results reveal that even without adhesion or contact to a substrate, Dd cells can still swim to follow chemoattractant signals. In this paper, we develop a modeling framework to investigate the chemotaxis induced amoeboid cell swimming dynamics. A minimal swimming system consists of one deformable Dd amoeboid cell and a dilute suspension of bacteria, and the bacteria produce chemoattractant signals that attract the Dd cell. We use themathematical amoeba modelto generate Dd cell deformation and solve the resulting low Reynolds number flows, and use a moving mesh based finite volume method to solve the reaction-diffusion-convection equation. Using the computational model, we show that chemotaxis guides a swimming Dd cell to follow and catch bacteria, while on the other hand, bacterial rheotaxis may help the bacteria to escape from the predator Dd cell.

Adhesion strength and contractility enable metastatic cells to become adurotactic

Significant changes in cell stiffness, contractility, and adhesion, i.e., mechanotype, are observed during a variety of biological processes. Whether cell mechanics merely change as a side effect of or driver for biological processes is still unclear. Here, we sort genotypically similar metastatic cancer cells into strongly adherent (SA) versus weakly adherent (WA) phenotypes to study how contractility and adhesion differences alter the ability of cells to sense and respond to gradients in material stiffness. We observe that SA cells migrate up a stiffness gradient, or durotax, while WA cells largely ignore the gradient, i.e., adurotax. Biophysical modeling and experimental validation suggest that differences in cell migration and durotaxis between weakly and strongly adherent cells are driven by differences in intra-cellular actomyosin activity. These results provide a direct relationship between cell phenotype and durotaxis and suggest how, unlike other senescent cells, metastatic cancer cells navigate against stiffness gradients.

Mathematical modelling in cell migration: tackling biochemistry in changing geometries

Directed cell migration poses a rich set of theoretical challenges. Broadly, these are concerned with (1) how cells sense external signal gradients and adapt; (2) how actin polymerisation is localised to drive the leading cell edge and Myosin-II molecular motors retract the cell rear; and (3) how the combined action of cellular forces and cell adhesion results in cell shape changes and net migration. Reaction-diffusion models for biological pattern formation going back to Turing have long been used to explain generic principles of gradient sensing and cell polarisation in simple, static geometries like a circle. In this minireview, we focus on recent research which aims at coupling the biochemistry with cellular mechanics and modelling cell shape changes. In particular, we want to contrast two principal modelling approaches: (1) interface tracking where the cell membrane, interfacing cell interior and exterior, is explicitly represented by a set of moving points in 2D or 3D space and (2) interface capturing. In interface capturing, the membrane is implicitly modelled analogously to a level line in a hilly landscape whose topology changes according to forces acting on the membrane. With the increased availability of high-quality 3D microscopy data of complex cell shapes, such methods will become increasingly important in data-driven, image-based modelling to better understand the mechanochemistry underpinning cell motion.

Towards understanding the messengers of extracellular space: Computational models of outside-in integrin reaction networks

The interactions between cells and their extracellular matrix (ECM) are critically important for homeostatic control of cell growth, proliferation, differentiation and apoptosis. Transmembrane integrin molecules facilitate the communication between ECM and the cell. Since the characterization of integrins in the late 1980s, there has been great advancement in understanding the function of integrins at different subcellular levels. However, the versatility in molecular pathways integrins are involved in, the high diversity in their interaction partners both outside and inside the cell as well as on the cell membrane and the short lifetime of events happening at the cell-ECM interface make it difficult to elucidate all the details regarding integrin function experimentally. To overcome the experimental challenges and advance the understanding of integrin biology, computational modeling tools have been used extensively. In this review, we summarize the computational models of integrin signaling while we explain the function of integrins at three main subcellular levels (outside the cell, cell membrane, cytosol). We also discuss how these computational modeling efforts can be helpful in other disciplines such as biomaterial design. As such, this review is a didactic modeling summary for biomaterial researchers interested in complementing their experimental work with computational tools or for seasoned computational scientists that would like to advance current in silico integrin models.

Membrane Tension Can Enhance Adaptation to Maintain Polarity of Migrating Cells

Migratory cells are known to adapt to environments that contain wide-ranging levels of chemoattractant. Although biochemical models of adaptation have been previously proposed, here, we discuss a different mechanism based on mechanosensing, in which the interaction between biochemical signaling and cell tension facilitates adaptation. We describe and analyze a model of mechanochemical-based adaptation coupling a mechanics-based physical model of cell tension coupled with the wave-pinning reaction-diffusion model for Rac GTPase activity. The mathematical analysis of this model, simulations of a simplified one-dimensional cell geometry, and two-dimensional finite element simulations of deforming cells reveal that as a cell protrudes under the influence of high stimulation levels, tension-mediated inhibition of Rac signaling causes the cell to polarize even when initially overstimulated. Specifically, tension-mediated inhibition of Rac activation, which has been experimentally observed in recent years, facilitates this adaptation by countering the high levels of environmental stimulation. These results demonstrate how tension-related mechanosensing may provide an alternative (and potentially complementary) mechanism for cell adaptation.

Stick-slip model for actin-driven cell protrusions, cell polarization, and crawling

Cell crawling requires the generation of intracellular forces by the cytoskeleton and their transmission to an extracellular substrate through specific adhesion molecules. Crawling cells show many features of excitable systems, such as spontaneous symmetry breaking and crawling in the absence of external cues, and periodic and propagating waves of activity. Mechanical instabilities in the active cytoskeleton network and feedback loops in the biochemical network of activators and repressors of cytoskeleton dynamics have been invoked to explain these dynamical features. Here, I show that the interplay between the dynamics of cell-substrate adhesion and linear cellular mechanics is sufficient to reproduce many nonlinear dynamical patterns observed in spreading and crawling cells. Using an analytical formalism of the molecular clutch model of cell adhesion, regulated by local mechanical forces, I show that cellular traction forces exhibit stick-slip dynamics resulting in periodic waves of protrusion/retraction and propagating waves along the cell edge. This can explain spontaneous symmetry breaking and polarization of spreading cells, leading to steady crawling or bipedal motion, and bistability, where persistent cell motion requires a sufficiently strong transient external stimulus. The model also highlights the role of membrane tension in providing the long-range mechanical communication across the cell required for symmetry breaking.

A Stochastic Modelling Framework for Single Cell Migration: Coupling Contractility and Focal Adhesions

The interaction of the actin cytoskeleton with cell–substrate adhesions is necessary for cell migration. While the trajectories of motile cells have a stochastic character, investigations of cell motility mechanisms rarely elaborate on the origins of the observed randomness. Here, guided by a few fundamental attributes of cell motility, I construct a minimal stochastic cell migration model from ground-up. The resulting model couples a deterministic actomyosin contractility mechanism with stochastic cell–substrate adhesion kinetics, and yields a well-defined piecewise deterministic process. Numerical simulations reproduce several experimentally observed results, including anomalous diffusion, tactic migration and contact guidance. This work provides a basis for the development of cell–cell collision and population migration models.

Periodic propagating waves coordinate RhoGTPase network dynamics at the leading and trailing edges during cell migration

Migrating cells need to coordinate distinct leading and trailing edge dynamics but the underlying mechanisms are unclear. Here, we combine experiments and mathematical modeling to elaborate the minimal autonomous biochemical machinery necessary and sufficient for this dynamic coordination and cell movement. RhoA activates Rac1 via DIA and inhibits Rac1 via ROCK, while Rac1 inhibits RhoA through PAK. Our data suggest that in motile, polarized cells, RhoA–ROCK interactions prevail at the rear, whereas RhoA-DIA interactions dominate at the front where Rac1/Rho oscillations drive protrusions and retractions. At the rear, high RhoA and low Rac1 activities are maintained until a wave of oscillatory GTPase activities from the cell front reaches the rear, inducing transient GTPase oscillations and RhoA activity spikes. After the rear retracts, the initial GTPase pattern resumes. Our findings show how periodic, propagating GTPase waves coordinate distinct GTPase patterns at the leading and trailing edge dynamics in moving cells.

Tuning Cell Motility via Cell Tension with a Mechanochemical Cell Migration Model

Cell migration is orchestrated by a complicated mechanochemical system. However, few cell migration models take into account the coupling between the biochemical network and mechanical factors. Here, we construct a mechanochemical cell migration model to study the cell tension effect on cell migration. Our model incorporates the interactions between Rac-GTP, Rac-GDP, F-actin, myosin, and cell tension, and it is very convenient in capturing the change of cell shape by taking the phase field approach. This model captures the characteristic features of cell polarization, cell shape change, and cell migration modes. It shows that cell tension inhibits migration ability monotonically when cells are applied with persistent external stimuli. On the other hand, if random internal noise is significant, the regulation of cell tension exerts a nonmonotonic effect on cell migration. Because the increase of cell tension hinders the formation of multiple protrusions, migration ability could be maximized at intermediate cell tension under random internal noise. These model predictions are consistent with our single-cell experiments and other experimental results.

Hydrodynamic effects on the motility of crawling eukaryotic cells

Eukaryotic cell motility is crucial during development, wound healing, the immune response, and cancer metastasis. Some eukaryotic cells can swim, but cells more commonly adhere to and crawl along the extracellular matrix. We study the relationship between hydrodynamics and adhesion that describe whether a cell is swimming, crawling, or combining these motions. Our simple model of a cell, based on the three-sphere swimmer, is capable of both swimming and crawling. As cell-matrix adhesion strength increases, the influence of hydrodynamics on migration diminishes. Cells with significant adhesion can crawl with speeds much larger than their nonadherent, swimming counterparts. We predict that, while most eukaryotic cells are in the strong-adhesion limit, increasing environment viscosity or decreasing cell-matrix adhesion could lead to significant hydrodynamic effects even in crawling cells. Signatures of hydrodynamic effects include a dependence of cell speed on the presence of a nearby substrate or interactions between noncontacting cells. These signatures will be suppressed at large adhesion strengths, but even strongly adherent cells will generate relevant fluid flows that will advect nearby passive particles and swimmers.

Simple Rho GTPase Dynamics Generate a Complex Regulatory Landscape Associated with Cell Shape

Migratory cells exhibit a variety of morphologically distinct responses to their environments that manifest in their cell shape. Some protrude uniformly to increase substrate contacts, others are broadly contractile, some polarize to facilitate migration, and yet others exhibit mixtures of these responses. Prior studies have identified a discrete collection of shapes that the majority of cells display and demonstrated that activity levels of the cytoskeletal regulators Rac1 and RhoA GTPase regulate those shapes. Here, we use computational modeling to assess whether known GTPase dynamics can give rise to a sufficient diversity of spatial signaling states to explain the observed shapes. Results show that the combination of autoactivation and mutually antagonistic cross talk between GTPases, along with the conservative membrane binding, generates a wide array of distinct homogeneous and polarized regulatory phenotypes that arise for fixed model parameters. From a theoretical perspective, these results demonstrate that simple GTPase dynamics can generate complex multistability in which six distinct stable steady states (three homogeneous and three polarized) coexist for a fixed set of parameters, each of which naturally maps to an observed morphology. From a biological perspective, although we do not explicitly model the cytoskeleton or resulting cell morphologies, these results, along with prior literature linking GTPase activity to cell morphology, support the hypothesis that GTPase signaling dynamics can generate the broad morphological characteristics observed in many migratory cell populations. Further, the observed diversity may be the result of cells populating a complex morphological landscape generated by GTPase regulation rather than being the result of intrinsic cell-cell variation. These results demonstrate that Rho GTPases may have a central role in regulating the broad characteristics of cell shape (e.g., expansive, contractile, polarized, etc.) and that shape heterogeneity may be (at least partly) a reflection of the rich signaling dynamics regulating the cytoskeleton rather than intrinsic cell heterogeneity.

Motility and morphodynamics of confined cells

We introduce a minimal hydrodynamic model of polarization, migration, and deformation of a biological cell confined between two parallel surfaces. In our model, the cell is driven out of equilibrium by an active cytsokeleton force that acts on the membrane. The cell cytoplasm, described as a viscous droplet in the Darcy flow regime, contains a diffusive solute that actively transduces the applied cytoskeleton force. While fairly simple and analytically tractable, this quasi-two-dimensional model predicts a range of compelling dynamic behaviours. A linear stability analysis of the system reveals that solute activity first destabilizes a global polarization-translation mode, prompting cell motility through spontaneous symmetry breaking. At higher activity, the system crosses a series of Hopf bifurcations leading to coupled oscillations of droplet shape and solute concentration profiles. At the nonlinear level, we find traveling-wave solutions associated with unique polarized shapes that resemble experimental observations. Altogether, this model offers an analytical paradigm of active deformable systems in which viscous hydrodynamics are coupled to diffusive force transducers.

A mechanical toy model linking cell-substrate adhesion to multiple cellular migratory responses

During cell migration, forces applied to a cell from its environment influence the motion. When the cell is placed on a substrate, such a force is provided by the cell-substrate adhesion. Modulation of adhesivity, often performed by the modulation of the substrate stiffness, tends to cause common responses for cell spreading, cell speed, persistence, and random motility coefficient. Although the reasons for the response of cell spreading and cell speed have been suggested, other responses are not well understood. In this study, we develop a simple toy model for cell migration driven by the relation of two forces: the adhesive force and the plasma membrane tension. The simplicity of the model allows us to perform the calculation not only numerically but also analytically, and the analysis provides formulas directly relating the adhesivity to cell spreading, persistence, and the random motility coefficient. Accordingly, the results offer a unified picture on the causal relations between those multiple cellular responses. In addition, cellular properties that would influence the migratory behavior are suggested.

Emerging themes and unifying concepts underlying cell behavior regulation by the pericellular space

Cells reside in a complex three-dimensional (3D) microenvironment where physical, chemical, and architectural features of the pericellular space regulate important cellular functions like migration, differentiation, and morphogenesis. A major goal of tissue engineering is to identify which properties of the pericellular space orchestrate these emergent cell behaviors and how. In this review, we highlight recent studies at the interface of biomaterials and single cell biophysics that are lending deeper insight towards this goal. Advanced methods have enabled the decoupling of architectural and mechanical features of the microenvironment, revealing multiple mechanisms of adhesion and mechanosensing modulation by biomaterials. Such studies are revealing important roles for pericellular space degradability, hydration, and adhesion competition in cell shape, volume, and differentiation regulation. STATEMENT OF SIGNIFICANCE: Cell fate and function are closely regulated by the local extracellular microenvironment. Advanced methods at the interface of single cell biophysics and biomaterials have shed new light on regulators of cell-pericellular space interactions by decoupling more features of the complex pericellular milieu than ever before. These findings lend deeper mechanistic insight into how biomaterials can be designed to fine-tune outcomes like differentiation, migration, and collective morphogenesis.

Fully coupled numerical model of actin treadmilling in the lamellipodium of the cell

Cells rely on an interplay of subcellular elements for motility and migration. Certain regions of motile cells, such as the lamellipodium, are made of a complex mixture of actin monomers and filaments, which polymerize at the front of the cell, close to the cell membrane, and depolymerize at the rear. The dynamic actin turnover induces the so-called intracellular retrograde flow, and it is a fundamental process for cell motility. Apart from some comprehensive mathematical models, the computational modelling of actin treadmilling has been based on simpler biophysical models. Here, we adopt a highly detailed theoretical model of the actin treadmilling process and develop a coupled unsteady finite element formulation. We clearly describe the structure and implementation of the coupled problem within the finite element method. Our numerical results show an excellent correlation with experimental results from literature and with previous models. We include time dependent effects and convective transport terms, which expose puzzling dynamics in the retrograde flow. We propose several biological scenarios to analyze the behavior of the actin treadmilling along space and time. We observed response times of the main density variables in the order of seconds. Compared with previous analytical solutions, which make assumptions related to convective transport, transient dynamics, and actin fluxes, the generic solution can have significant influence on the retrograde flow. All together, our results unveil a promising applicability of classical finite element methods to derive an in silico testing platform for the actin treadmilling processes in motile cells, which could allow for an extension to other biophysical effects.

Theoretical Models of Neural Development

Constructing a functioning nervous system requires the precise orchestration of a vast array of mechanical, molecular, and neural-activity-dependent cues. Theoretical models can play a vital role in helping to frame quantitative issues, reveal mathematical commonalities between apparently diverse systems, identify what is and what is not possible in principle, and test the abilities of specific mechanisms to explain the data. This review focuses on the progress that has been made over the last decade in our theoretical understanding of neural development.

Mechanochemical modeling of neutrophil migration based on four signaling layers, integrin dynamics, and substrate stiffness

Directional neutrophil migration during human immune responses is a highly coordinated process regulated by both biochemical and biomechanical environments. In this paper, we developed an integrative mathematical model of neutrophil migration using a lattice Boltzmann-particle method built in-house to solve the moving boundary problem with spatiotemporal regulation of biochemical components. The mechanical features of the cell cortex are modeled by a series of spring-connected nodes representing discrete cell-substrate adhesive sites. The intracellular signaling cascades responsible for cytoskeletal remodeling [e.g., small GTPases, phosphoinositide-3-kinase (PI3K), and phosphatase and tensin homolog] are built based on our previous four-layered signaling model centered on the bidirectional molecular transport mechanism and implemented as reaction-diffusion equations. Focal adhesion dynamics are determined by force-dependent integrin-ligand binding kinetics and integrin recycling and are thus integrated with cell motion. Using numerical simulations, the model reproduces the major features of cell migration in response to uniform and gradient biochemical stimuli based on the quantitative spatiotemporal regulation of signaling molecules, which agree with experimental observations. The existence of multiple types of integrins with different binding kinetics could act as an adaptation mechanism for substrate stiffness. Moreover, cells can perform reversal, U-turn, or lock-on behaviors depending on the steepness of the reversal biochemical signals received. Finally, this model is also applied to predict the responses of mutants in which PTEN is overexpressed or disrupted.

Computational modeling of single-cell mechanics and cytoskeletal mechanobiology

Cellular cytoskeletal mechanics plays a major role in many aspects of human health from organ development to wound healing, tissue homeostasis and cancer metastasis. We summarize the state-of-the-art techniques for mathematically modeling cellular stiffness and mechanics and the cytoskeletal components and factors that regulate them. We highlight key experiments that have assisted model parameterization and compare the advantages of different models that have been used to recapitulate these experiments. An overview of feed-forward mechanisms from signaling to cytoskeleton remodeling is provided, followed by a discussion of the rapidly growing niche of encapsulating feedback mechanisms from cytoskeletal and cell mechanics to signaling. We discuss broad areas of advancement that could accelerate research and understanding of cellular mechanobiology. A precise understanding of the molecular mechanisms that affect cell and tissue mechanics and function will underpin innovations in medical device technologies of the future. WIREs Syst Biol Med 2018, 10:e1407. doi: 10.1002/wsbm.1407 This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

Computational Modeling of the Dynamics of Spatiotemporal Rho GTPase Signaling: A Systematic Review

The Rho family of GTPases are known to play pivotal roles in the regulation of fundamental cellular processes, ranging from cell migration and polarity to wound healing and regulation of actin cytoskeleton. Over the past decades, accumulating experimental work has increasingly mapped out the mechanistic details and interactions between members of the family and their regulators, establishing detailed interaction circuits within the Rho GTPase signaling network. These circuits have served as a vital foundation based on which a multitude of mathematical models have been developed to explain experimental data, gain deeper insights into the biological phenomenon they describe, as well as make new testable predictions and hypotheses. Due to the diverse nature and purpose of these models, they often vary greatly in size, scope, complexity, and formulation. Here, we provide a systematic, categorical, and comprehensive account of the recent modeling studies of Rho family GTPases, with an aim to offer a broad perspective of the field. The modeling limitations and possible future research directions are also discussed.

Cellular mechanosensing of the biophysical microenvironment: A review of mathematical models of biophysical regulation of cell responses

Cells in vivo reside within complex microenvironments composed of both biochemical and biophysical cues. The dynamic feedback between cells and their microenvironments hinges upon biophysical cues that regulate critical cellular behaviors. Understanding this regulation from sensing to reaction to feedback is therefore critical, and a large effort is afoot to identify and mathematically model the fundamental mechanobiological mechanisms underlying this regulation. This review provides a critical perspective on recent progress in mathematical models for the responses of cells to the biophysical cues in their microenvironments, including dynamic strain, osmotic shock, fluid shear stress, mechanical force, matrix rigidity, porosity, and matrix shape. The review highlights key successes and failings of existing models, and discusses future opportunities and challenges in the field.

Cell protrusion and retraction driven by fluctuations in actin polymerization: A two-dimensional model

Animal cells that spread onto a surface often rely on actin-rich lamellipodial extensions to execute protrusion. Many cell types recently adhered on a two-dimensional substrate exhibit protrusion and retraction of their lamellipodia, even though the cell is not translating. Travelling waves of protrusion have also been observed, similar to those observed in crawling cells. These regular patterns of protrusion and retraction allow quantitative analysis for comparison to mathematical models. The periodic fluctuations in leading edge position of XTC cells have been linked to excitable actin dynamics using a one-dimensional model of actin dynamics, as a function of arc-length along the cell. In this work we extend this earlier model of actin dynamics into two dimensions (along the arc-length and radial directions of the cell) and include a model membrane that protrudes and retracts in response to the changing number of free barbed ends of actin filaments near the membrane. We show that if the polymerization rate at the barbed ends changes in response to changes in their local concentration at the leading edge and/or the opposing force from the cell membrane, the model can reproduce the patterns of membrane protrusion and retraction seen in experiment. We investigate both Brownian ratchet and switch-like force-velocity relationships between the membrane load forces and actin polymerization rate. The switch-like polymerization dynamics recover the observed patterns of protrusion and retraction as well as the fluctuations in F-actin concentration profiles. The model generates predictions for the behavior of cells after local membrane tension perturbations.

A free-boundary model of a motile cell explains turning behavior

To understand shapes and movements of cells undergoing lamellipodial motility, we systematically explore minimal free-boundary models of actin-myosin contractility consisting of the force-balance and myosin transport equations. The models account for isotropic contraction proportional to myosin density, viscous stresses in the actin network, and constant-strength viscous-like adhesion. The contraction generates a spatially graded centripetal actin flow, which in turn reinforces the contraction via myosin redistribution and causes retraction of the lamellipodial boundary. Actin protrusion at the boundary counters the retraction, and the balance of the protrusion and retraction shapes the lamellipodium. The model analysis shows that initiation of motility critically depends on three dimensionless parameter combinations, which represent myosin-dependent contractility, a characteristic viscosity-adhesion length, and a rate of actin protrusion. When the contractility is sufficiently strong, cells break symmetry and move steadily along either straight or circular trajectories, and the motile behavior is sensitive to conditions at the cell boundary. Scanning of a model parameter space shows that the contractile mechanism of motility supports robust cell turning in conditions where short viscosity-adhesion lengths and fast protrusion cause an accumulation of myosin in a small region at the cell rear, destabilizing the axial symmetry of a moving cell.

To understand shapes and movements of simple motile cells, we systematically explore minimal models describing a cell as a two-dimensional actin-myosin gel with a free boundary. The models account for actin-myosin contraction balanced by viscous stresses in the actin gel and uniform adhesion. The myosin contraction causes the lamellipodial boundary to retract. Actin protrusion at the boundary counters the retraction, and the balance of protrusion and retraction shapes the cell. The models reproduce a variety of motile shapes observed experimentally. The analysis shows that the mechanical state of a cell depends on a small number of parameters. We find that when the contractility is sufficiently strong, cells break symmetry and move steadily along either straight or circular trajectory. Scanning model parameters shows that the contractile mechanism of motility supports robust cell turning behavior in conditions where deformable actin gel and fast protrusion destabilize the axial symmetry of a moving cell.

BioFlow: a non-invasive, image-based method to measure speed, pressure and forces inside living cells

Cell motility is governed by a complex molecular machinery that converts physico-chemical cues into whole-cell movement. Understanding the underlying biophysical mechanisms requires the ability to measure physical quantities inside the cell in a simple, reproducible and preferably non-invasive manner. To this end, we developed BioFlow, a computational mechano-imaging method and associated software able to extract intracellular measurements including pressure, forces and velocity everywhere inside freely moving cells in two and three dimensions with high spatial resolution in a non-invasive manner. This is achieved by extracting the motion of intracellular material observed using fluorescence microscopy, while simultaneously inferring the parameters of a given theoretical model of the cell interior. We illustrate the power of BioFlow in the context of amoeboid cell migration, by modelling the intracellular actin bulk flow of the parasite Entamoeba histolytica using fluid dynamics, and report unique experimental measures that complement and extend both theoretical estimations and invasive experimental measures. Thanks to its flexibility, BioFlow is easily adaptable to other theoretical models of the cell, and alleviates the need for complex or invasive experimental conditions, thus constituting a powerful tool-kit for mechano-biology studies. BioFlow is open-source and freely available via the Icy software.

Epithelial/mesenchymal plasticity: how have quantitative mathematical models helped improve our understanding?

Phenotypic plasticity, the ability of cells to reversibly alter their phenotypes in response to signals, presents a significant clinical challenge to treating solid tumors. Tumor cells utilize phenotypic plasticity to evade therapies, metastasize, and colonize distant organs. As a result, phenotypic plasticity can accelerate tumor progression. A well‐studied example of phenotypic plasticity is the bidirectional conversions among epithelial, mesenchymal, and hybrid epithelial/mesenchymal (E/M) phenotype(s). These conversions can alter a repertoire of cellular traits associated with multiple hallmarks of cancer, such as metabolism, immune evasion, invasion, and metastasis. To tackle the complexity and heterogeneity of these transitions, mathematical models have been developed that seek to capture the experimentally verified molecular mechanisms and act as ‘hypothesis‐generating machines’. Here, we discuss how these quantitative mathematical models have helped us explain existing experimental data, guided further experiments, and provided an improved conceptual framework for understanding how multiple intracellular and extracellular signals can drive E/M plasticity at both the single‐cell and population levels. We also discuss the implications of this plasticity in driving multiple aggressive facets of tumor progression.

A mathematical model coupling polarity signaling to cell adhesion explains diverse cell migration patterns

Protrusion and retraction of lamellipodia are common features of eukaryotic cell motility. As a cell migrates through its extracellular matrix (ECM), lamellipod growth increases cell-ECM contact area and enhances engagement of integrin receptors, locally amplifying ECM input to internal signaling cascades. In contrast, contraction of lamellipodia results in reduced integrin engagement that dampens the level of ECM-induced signaling. These changes in cell shape are both influenced by, and feed back onto ECM signaling. Motivated by experimental observations on melanoma cells lines (1205Lu and SBcl2) migrating on fibronectin (FN) coated topographic substrates (anisotropic post-density arrays), we probe this interplay between intracellular and ECM signaling. Experimentally, cells exhibited one of three lamellipodial dynamics: persistently polarized, random, or oscillatory, with competing lamellipodia oscillating out of phase (Park et al., 2017). Pharmacological treatments, changes in FN density, and substrate topography all affected the fraction of cells exhibiting these behaviours. We use these observations as constraints to test a sequence of hypotheses for how intracellular (GTPase) and ECM signaling jointly regulate lamellipodial dynamics. The models encoding these hypotheses are predicated on mutually antagonistic Rac-Rho signaling, Rac-mediated protrusion (via activation of Arp2/3 actin nucleation) and Rho-mediated contraction (via ROCK phosphorylation of myosin light chain), which are coupled to ECM signaling that is modulated by protrusion/contraction. By testing each model against experimental observations, we identify how the signaling layers interact to generate the diverse range of cell behaviors, and how various molecular perturbations and changes in ECM signaling modulate the fraction of cells exhibiting each. We identify several factors that play distinct but critical roles in generating the observed dynamic: (1) competition between lamellipodia for shared pools of Rac and Rho, (2) activation of RhoA by ECM signaling, and (3) feedback from lamellipodial growth or contraction to cell-ECM contact area and therefore to the ECM signaling level.

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.

Cell Sorting and Noise-Induced Cell Plasticity Coordinate to Sharpen Boundaries between Gene Expression Domains

A fundamental question in biology is how sharp boundaries of gene expression form precisely in spite of biological variation/noise. Numerous mechanisms position gene expression domains across fields of cells (e.g. morphogens), but how these domains are refined remains unclear. In some cases, domain boundaries sharpen through differential adhesion-mediated cell sorting. However, boundaries can also sharpen through cellular plasticity, with cell fate changes driven by up- or down-regulation of gene expression. In this context, we have argued that noise in gene expression can help cells transition to the correct fate. Here we investigate the efficacy of cell sorting, gene expression plasticity, and their combination in boundary sharpening using multi-scale, stochastic models. We focus on the formation of hindbrain segments (rhombomeres) in the developing zebrafish as an example, but the mechanisms investigated apply broadly to many tissues. Our results indicate that neither sorting nor plasticity is sufficient on its own to sharpen transition regions between different rhombomeres. Rather the two have complementary strengths and weaknesses, which synergize when combined to sharpen gene expression boundaries.

In many developing systems, chemical gradients control the formation of segmental domains of gene expression, specifying distinct domains that go on to form different tissues and structures, in a concentration-dependent manner. These gradients are noisy however, raising the question of how sharply delineated boundaries between distinct segments form. It is crucial that developing systems be able to cope with stochasticity and generate well-defined boundaries between different segmented domains. Previous work suggests that cell sorting and cellular plasticity help sharpen boundaries between segments. However, it remains unclear how effective each of these mechanisms is and what their role in sharpening may be. Motivated by recent experimental observations, we construct a hybrid stochastic model to investigate these questions. We find that neither mechanism is sufficient on its own to sharpen boundaries between different segments. Rather, results indicate each has its own strengths and weaknesses, and that they work together synergistically to promote the development of precise, well defined segment boundaries. Formation of segmented rhombomeres in the zebrafish hindbrain, which later form different components of the central nervous system, is a motivating case for this study.

Crawling and turning in a minimal reaction-diffusion cell motility model: Coupling cell shape and biochemistry

We study a minimal model of a crawling eukaryotic cell with a chemical polarity controlled by a reaction-diffusion mechanism describing Rho GTPase dynamics. The size, shape, and speed of the cell emerge from the combination of the chemical polarity, which controls the locations where actin polymerization occurs, and the physical properties of the cell, including its membrane tension. We find in our model both highly persistent trajectories, in which the cell crawls in a straight line, and turning trajectories, where the cell transitions from crawling in a line to crawling in a circle. We discuss the controlling variables for this turning instability and argue that turning arises from a coupling between the reaction-diffusion mechanism and the shape of the cell. This emphasizes the surprising features that can arise from simple links between cell mechanics and biochemistry. Our results suggest that similar instabilities may be present in a broad class of biochemical descriptions of cell polarity.

Free boundary problem for cell protrusion formations: theoretical and numerical aspects

In this paper, a free boundary problem for cell protrusion formation is studied theoretically and numerically. The cell membrane is precisely described thanks to a level set function, whose motion is due to specific signalling pathways. The aim is to model the chemical interactions between the cell and its environment, in the process of invadopodia or pseudopodia formation. The model consists of Laplace equation with Dirichlet condition inside the cell coupled to Laplace equation with Neumann condition in the outer domain. The actin polymerization is accounted for as the gradient of the inner signal, which drives the motion of the interface. We prove the well-posedness of our free boundary problem under a sign condition on the datum. This criterion ensures the consistency of the model, and provides conditions to focus on for any enrichment of the model. We then propose a new first order Cartesian finite-difference method to solve the problem. We eventually exhibit the main biological features that can be accounted for by the model: the formation of thin and elongated protrusions as for invadopodia, or larger protrusion as for pseudopodia, depending on the source term in the equation. The model provides the theoretical and numerical grounds for single cell migration modeling, whose formulation is valid in 2D and 3D. In particular, specific chemical reactions that occurred at the cell membrane could be precisely described in forthcoming works.

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.

Single-Cell Migration in Complex Microenvironments: Mechanics and Signaling Dynamics

Cells are highly dynamic and mechanical automata powered by molecular motors that respond to external cues. Intracellular signaling pathways, either chemical or mechanical, can be activated and spatially coordinated to induce polarized cell states and directional migration. Physiologically, cells navigate through complex microenvironments, typically in three-dimensional (3D) fibrillar networks. In diseases, such as metastatic cancer, they invade across physiological barriers and remodel their local environments through force, matrix degradation, synthesis, and reorganization. Important external factors such as dimensionality, confinement, topographical cues, stiffness, and flow impact the behavior of migrating cells and can each regulate motility. Here, we review recent progress in our understanding of single-cell migration in complex microenvironments.

A mathematical model of GTPase pattern formation during single-cell wound repair

Rho GTPases are regulatory proteins whose patterns on the surface of a cell affect cell polarization, cell motility and repair of single-cell wounds. The stereotypical patterns formed by two such proteins, Rho and Cdc42, around laser-injured frog oocytes permit experimental analysis of GTPase activation, inactivation, segregation and crosstalk. Here, we review the development and analysis of a spatial model of GTPase dynamics that describe the formation of concentric zones of Rho and Cdc42 activity around wounds, and describe how this model has provided insights into the roles of the GTPase effector molecules protein kinase C (PKCβ and PKCη) and guanosine nucleotide dissociation inhibitor (GDI) in the wound response. We further demonstrate how the use of a 'sharp switch' model approximation in combination with bifurcation analysis can aid mapping the model behaviour in parameter space (approximate results confirmed with numerical simulation methods). Using these methods in combination with experimental manipulation of PKC activity (PKC overexpression (OE) and dominant negative conditions), we have shown that: (i) PKCβ most probably acts by enhancing existing positive feedbacks (from Rho to itself via the guanosine nucleotide exchange factor domain of Abr, and from Cdc42 to itself), (ii) PKCη most probably increases basal rates of inactivation (or possibly decreases basal rates of activation) of Rho and Cdc42, and (iii) the graded distribution of PKCη and its effect on initial Rho activity accounts for inversion of zones in a fraction (20%) of PKCη OE cells. Finally, we speculate that GDIs (which sequester GTPases) may have a critical role in defining the spatial domain, where the wound response may occur. This paper provides a more thorough exposition of the methods of analysis used in the investigation, whereas previous work on this topic was addressed to biologists and abbreviated such discussion.

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.

Modeling Excitable Dynamics of Chemotactic Networks

The study of chemotaxis has benefited greatly from computational models that describe the response of cells to chemoattractant stimuli. These models must keep track of spatially and temporally varying distributions of numerous intracellular species. Moreover, recent evidence suggests that these are not deterministic interactions, but also include the effect of stochastic variations that trigger an excitable network. In this chapter we illustrate how to create simulations of excitable networks using the Virtual Cell modeling environment.

Cell-ECM Interactions in Tumor Invasion

The cancer cells obtain their invasion potential not only by genetic mutations, but also by changing their cellular biophysical and biomechanical features and adapting to the surrounding microenvironments. The extracellular matrix, as a crucial component of the tumor microenvironment, provides the mechanical support for the tissue, mediates the cell-microenvironment interactions, and plays a key role in cancer cell invasion. The biomechanics of the extracellular matrix, particularly collagen, have been extensively studied in the biomechanics community. Cell migration has also enjoyed much attention from both the experimental and modeling efforts. However, the detailed mechanistic understanding of tumor cell-ECM interactions, especially during cancer invasion, has been unclear. This chapter reviews the recent advances in the studies of ECM biomechanics, cell migration, and cell-ECM interactions in the context of cancer invasion.

Computational modelling of epidermal stratification highlights the importance of asymmetric cell division for predictable and robust layer formation

Skin is a complex organ tasked with, among other functions, protecting the body from the outside world. Its outermost protective layer, the epidermis, is comprised of multiple cell layers that are derived from a single-layered ectoderm during development. Using a new stochastic, multi-scale computational modelling framework, the anisotropic subcellular element method, we investigate the role of cell morphology and biophysical cell-cell interactions in the formation of this layered structure. This three-dimensional framework describes interactions between collections of hundreds to thousands of cells and (i) accounts for intracellular structure and morphology, (ii) easily incorporates complex cell-cell interactions and (iii) can be efficiently implemented on parallel architectures. We use this approach to construct a model of the developing epidermis that accounts for the internal polarity of ectodermal cells and their columnar morphology. Using this model, we show that cell detachment, which has been previously suggested to have a role in this process, leads to unpredictable, randomized stratification and that this cannot be abrogated by adjustment of cell-cell adhesion interaction strength. Polarized distribution of cell adhesion proteins, motivated by epithelial polarization, can however eliminate this detachment, and in conjunction with asymmetric cell division lead to robust and predictable development.

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.

The Interplay between Wnt Mediated Expansion and Negative Regulation of Growth Promotes Robust Intestinal Crypt Structure and Homeostasis

The epithelium of the small intestinal crypt, which has a vital role in protecting the underlying tissue from the harsh intestinal environment, is completely renewed every 4–5 days by a small pool of stem cells at the base of each crypt. How is this renewal controlled and homeostasis maintained, particularly given the rapid nature of this process? Here, based on the recent observations from in vitro “mini gut” studies, we use a hybrid stochastic model of the crypt to investigate how exogenous niche signaling (from Wnt and BMP) combines with auto-regulation to promote homeostasis. This model builds on the sub-cellular element method to account for the three-dimensional structure of the crypt, external regulation by Wnt and BMP, internal regulation by Notch signaling, as well as regulation by internally generated diffusible signals. Results show that Paneth cell derived Wnt signals, which have been observed experimentally to sustain crypts in cultured organs, have a dramatically different influence on niche dynamics than does mesenchyme derived Wnt. While this signaling can indeed act as a redundant backup to the exogenous gradient, it introduces a positive feedback that destabilizes the niche and causes its uncontrolled expansion. We find that in this setting, BMP has a critical role in constraining this expansion, consistent with observations that its removal leads to crypt fission. Further results also point to a new hypothesis for the role of Ephrin mediated motility of Paneth cells, specifically that it is required to constrain niche expansion and maintain the crypt’s spatial structure. Combined, these provide an alternative view of crypt homeostasis where the niche is in a constant state of expansion and the spatial structure of the crypt arises as a balance between this expansion and the action of various sources of negative regulation that hold it in check.

The small intestinal epithelium, like our skin, is constantly being renewed. In the intestine however, this epithelium is exposed to the harsh digestive environment, necessitating much more rapid renewal. Remarkably, the entire epithelium is renewed every 4–5 days. This raises the question, how can the size and structure of this tissue be maintained given this pace. Motivated by recent experimental observations, we construct a three-dimensional, hybrid stochastic model to investigate the mechanisms responsible for homeostasis of this tissue. We find that there are redundant signals created by both the epithelium itself and surrounding tissues that act in parallel to maintain epithelial structure. This redundancy comes at a price however: it introduces the possibility of explosive stem cell population growth. Additional results suggest that other signals along with choreographed motion of cells are responsible for repressing this expansion. Taken together, our results provide a novel hypothesis for how robust but fast renewal of the crypt is achieved: as a balance between expansion, which drives fast renewal and repression, which holds that expansion in check to maintain the crypt’s structure.

Targeted Proteomics-Driven Computational Modeling of Macrophage S1P Chemosensing

Osteoclasts are monocyte-derived multinuclear cells that directly attach to and resorb bone. Sphingosine-1-phosphate (S1P)(1) regulates bone resorption by functioning as both a chemoattractant and chemorepellent of osteoclast precursors through two G-protein coupled receptors that antagonize each other in an S1P-concentration-dependent manner. To quantitatively explore the behavior of this chemosensing pathway, we applied targeted proteomics, transcriptomics, and rule-based pathway modeling using the Simmune toolset. RAW264.7 cells (a mouse monocyte/macrophage cell line) were used as model osteoclast precursors, RNA-seq was used to identify expressed target proteins, and selected reaction monitoring (SRM) mass spectrometry using internal peptide standards was used to perform absolute abundance measurements of pathway proteins. The resulting transcript and protein abundance values were strongly correlated. Measured protein abundance values, used as simulation input parameters, led to in silico pathway behavior matching in vitro measurements. Moreover, once model parameters were established, even simulated responses toward stimuli that were not used for parameterization were consistent with experimental findings. These findings demonstrate the feasibility and value of combining targeted mass spectrometry with pathway modeling for advancing biological insight.