The process of chemotaxis of living cells is complex. Cells follow gradients of an external signal because the interior of the cells gets polarized. The description of the exterior and the interior of the cell together with its motion for the convenient realization of the computational modeling of the whole process is a complex technical problem. Here, we employ a phase field model to characterize the interior of the cell, permitting the integration of stochastic partial differential equations, responsible for the polarization in the interior of the cell, and simultaneously, the calculation of the shape deformations of the cell, including its locomotion. We detail the mathematical description of the process and the procedure to calculate numerically the phase field with a simple reaction-diffusion equation for a single concentration.

This paper is concerned with the existence of transition fronts for a one-dimensional two patch model with KPP reaction terms. Density and flux conditions are imposed at the interface between the two patches. We first construct a pair of suitable super- and sub solutions by making full use of information of the leading edges of two KPP fronts and gluing them through the interface conditions. Then, an entire solution obtained thanks to a limiting argument is shown to be a transition front moving from one patch to the other one. This propagating solution admits asymptotic past and future speeds, and it connects two different fronts, each associated with one of the two patches. The paper thus provides the first example of a transition front for a KPP-type two-patch model with interface conditions.

To computationally investigate the recent experimental finding such that extracellular ATP release caused by exogeneous mechanical forces promote wound closure, we introduce a mathematical model, the Cellular Potts Model (CPM), which is a popular discretized model on a lattice, where the movement of a \"cell\" is determined by a Monte Carlo procedure. In the experiment, it was observed that there is mechanosensitive ATP release from the leading cells facing the wound gap and the subsequent extracellular Ca2+ influx. To model these phenomena, the Reaction-Diffusion equations for extracellular ATP and intracellular Ca2+ concentrations are adopted and combined with CPM, where we also add a polarity term because the cell migration is enhanced in the case of ATP release. From the numerical simulations using this hybrid model, we discuss effects of the collective cell migration due to the ATP release and the Ca2+ influx caused by the mechanical forces and the consequent promotion of wound closure.

The bioenergetic system of calcium ([Ca2+]), inositol 1, 4, 5-trisphophate (IP3) and nitric oxide (NO) regulate the diverse mechanisms in neurons. The dysregulation in any or all of the calcium, IP3 and nitric oxide dynamics may cause neurotoxicity and cell death. Few studies are noted in the literature on the interactions of two systems like [Ca2+] with IP3 and [Ca2+] with nitric oxide in neuron cells, which gives limited insights into regulatory and dysregulatory processes in neuron cells. But, no study is available on the cross talk in dynamics of three systems [Ca2+], IP3 and NO in neurons. Thus, the cross talk in the system dynamics of [Ca2+], IP3 and NO regulation processes in neurons have been studied using mathematical model. The two-way feedback process between [Ca2+] and IP3 and two-way feedback process between [Ca2+] and NO through cyclic guanosine monophosphate (cGMP) with plasmalemmal [Ca2+]-ATPase (PMCA) have been incorporated in the proposed model. This coupling handles the indirect two-way feedback process between IP3 and nitric oxide in neuronal cells automatically. The numerical outcomes were acquired by employing the finite element method (FEM) with the Crank-Nicholson scheme (CNS). The present model incorporating the sodium-calcium exchanger (NCX) and voltage-gated calcium channel (VGCC) provides novel insights into the various regulatory and dysregulatory processes due to buffer, IP3-receptor, ryanodine receptor, cGMP kinetics through PMCA channel, etc. and their impacts on the interactive spatiotemporal system dynamics of [Ca2+], IP3 and NO in neurons. It is concluded that the behavior of different crucial mechanisms is quite different for interactions of two systems of [Ca2+] and NO and the interactions of three systems of [Ca2+], IP3 and nitric oxide in neuronal cell due to mutual regulatory adjustments. The association of several neurological disorders with the alterations in calcium, IP3 and NO has been explored in neurons.

One of the most crucial and lethal characteristics of solid tumors is represented by the increased ability of cancer cells to migrate and invade other organs during the so-called metastatic spread. This is allowed thanks to the production of matrix metalloproteinases (MMPs), enzymes capable of degrading a type of collagen abundant in the basal membrane separating the epithelial tissue from the connective one. In this work, we employ a synergistic experimental and mathematical modelling approach to explore the invasion process of tumor cells. A mathematical model composed of reaction-diffusion equations describing the evolution of the tumor cells density on a gelatin substrate, MMPs enzymes concentration and the degradation of the gelatin is proposed. This is completed with a calibration strategy. We perform a sensitivity analysis and explore a parameter estimation technique both on synthetic and experimental data in order to find the optimal parameters that describe the in vitro experiments. A comparison between numerical and experimental solutions ends the work.

The vital participation of Ca2+ in human organ functions such as muscular contractions, heartbeat, brain functionality, skeletal activity, etc, motivated the scientists to thoroughly research the mechanisms of calcium (Ca2+) signalling in distinct human cells. Ca2+, inositol triphosphate (IP3), and adenosine triphosphate (ATP) play important roles in cell signaling and physiological processes. ATP and its derivatives are hypothesized to be important in the pathogenic process that leads to fibrotic illnesses like fibrosis. Fluctuations in Ca2+ and IP3 in a fibroblast cell influence ATP production. To date, no evidence of coupled Ca2+ and IP3 mechanics regulating ATP generation in a fibroblast cell during fibrotic disease has been found. The current work suggests an integrated mechanism for Ca2+ and IP3 dynamics in a fibroblast cell that regulates ATP generation. Simulation has been carried out using the finite element approach. The mechanics of interdependent systems findings vary dramatically from the results of basic independent system mechanics and give fresh information about the two systems\' activities. The numerical results provide new insights into the impacts of disturbances in source influx, the serca pump, and buffers on interdependent Ca2+ and IP3 dynamics and ATP synthesis in a fibroblast cell. According to the findings of this study, fibrotic disorders cannot be attributed solely to disruptions in the processes of calcium signaling mechanics but also to disruptions in IP3 regulation mechanisms affecting the regulation of calcium in the fibroblast cell and ATP release.

The inflammatory phase is an important event in the skin wound healing process. The deposition of granulation tissue in the wound bed and the rebuilding of the vascular network occur as inflammation diminishes. An angiogenic component in the formation of granulation tissue is the secretion of vascular endothelial growth factor, which assists in the chemotaxis, proliferation, and replication of fibroblasts. In this paper, we develop a mathematical model of skin wound healing angiogenic factors based on inflammatory cells (macrophages and neutrophils) and mediators (interleukin 6 and interleukin 10). We highlight the importance of this process in vascular endothelial growth factor release and in the formation of new capillary tips. We used a mathematical model of partial differential equations based on the reaction-diffusion-advection equations. In order to calibrate the parameters, we considered an in vivo model composed by four treatments: hydroalcoholic extract and oil-resin of Copaifera langsdorffii at 10% concentration, collagenase, and Lanette cream. Using the laboratory data for the wound edge, our mathematical model estimated the values of vascular endothelial growth factor concentration, and tips density in the center of the wound with a maximum error of 2.9%, and predicted healing time required for each treatment. The region of viability for the parameters, in the proposed model, was found through numerical simulations from the Interleukin 6 and 10 dysregulation and we obtained that, among the parameters analyzed, the greatest influencer in the dynamics of the system is the one, which represents the production of Interleukin 10 during phagocytosis.

In an arid or semi-arid environment, precipitation plays a vital role in vegetation growth. Recent researches reveal that the response of vegetation growth to precipitation has a lag effect. To explore the mechanism behind the lag phenomenon, we propose and investigate a water-vegetation model with spatiotemporal nonlocal effects. It is shown that the temporal kernel function does not affect Turing bifurcation. For better understanding the influences of lag effect and nonlocal competition on the vegetation pattern formation, we choose some special kernel functions and obtain some insightful results: (i) Time delay does not trigger the vegetation pattern formation, but can postpone the evolution of vegetation. In addition, in the absence of diffusion, time delay can induce the occurrence of stability switches, while in the presence of diffusion, spatially nonhomogeneous time-periodic solutions may emerge, but there are no stability switches; (ii) The spatial nonlocal interaction may trigger the pattern onset for small diffusion ratio of water and vegetation, and can change the number and size of isolated vegetation patches for large diffusion ratio. (iii) The interaction between time delay and spatial nonlocal competition may induce the emergence of traveling wave patterns, so that the vegetation remains periodic in space, but is oscillating in time. These results demonstrate that precipitation can significantly affect the growth and spatial distribution of vegetation.

The emergence of new variants of concern (VOCs) of the SARS-CoV-2 infection is one of the main factors of epidemic progression. Their development can be characterized by three critical stages: virus mutation leading to the appearance of new viable variants; the competition of different variants leading to the production of a sufficiently large number of copies; and infection transmission between individuals and its spreading in the population. The first two stages take place at the individual level (infected individual), while the third one takes place at the population level with possible competition between different variants. This work is devoted to the mathematical modeling of the first two stages of this process: the emergence of new variants and their progression in the epithelial tissue with a possible competition between them. The emergence of new virus variants is modeled with non-local reaction-diffusion equations describing virus evolution and immune escape in the space of genotypes. The conditions of the emergence of new virus variants are determined by the mutation rate, the cross-reactivity of the immune response, and the rates of virus replication and death. Once different variants emerge, they spread in the infected tissue with a certain speed and viral load that can be determined through the parameters of the model. The competition of different variants for uninfected cells leads to the emergence of a single dominant variant and the elimination of the others due to competitive exclusion. The dominant variant is the one with the maximal individual spreading speed. Thus, the emergence of new variants at the individual level is determined by the immune escape and by the virus spreading speed in the infected tissue.

Early embryogenesis is characterized by rapid and synchronous cleavage divisions, which are often controlled by wave-like patterns of Cdk1 activity. Two mechanisms have been proposed for mitotic waves: sweep and trigger waves.1,2 The two mechanisms give rise to different wave speeds, dependencies on physical and molecular parameters, and spatial profiles of Cdk1 activity: upward sweeping gradients versus traveling wavefronts. Both mechanisms hinge on the transient bistability governing the cell cycle and are differentiated by the speed of the cell-cycle progression: sweep and trigger waves arise for rapid and slow drives, respectively. Here, using quantitative imaging of Cdk1 activity and theory, we illustrate that sweep waves are the dominant mechanism in Drosophila embryos and test two fundamental predictions on the transition from sweep to trigger waves. We demonstrate that sweep waves can be turned into trigger waves if the cell cycle is slowed down genetically or if significant delays in the cell-cycle progression are introduced across the embryo by altering nuclear density. Our genetic experiments demonstrate that Polo kinase is a major rate-limiting regulator of the blastoderm divisions, and genetic perturbations reducing its activity can induce the transition from sweep to trigger waves. Furthermore, we show that changes in temperature cause an essentially uniform slowdown of interphase and mitosis. That results in sweep waves being observed across a wide temperature range despite the cell-cycle durations being significantly different. Collectively, our combination of theory and experiments elucidates the nature of mitotic waves in Drosophila embryogenesis, their control mechanisms, and their mutual transitions.