We introduce a free boundary model to study the effect of vesicle transport onto neurite growth. It consists of systems of drift-diffusion equations describing the evolution of the density of antero- and retrograde vesicles in each neurite coupled to reservoirs located at the soma and the growth cones of the neurites, respectively. The model allows for a change of neurite length as a function of the vesicle concentration in the growth cones. After establishing existence and uniqueness for the time-dependent problem, we briefly comment on possible types of stationary solutions. Finally, we provide numerical studies on biologically relevant scales using a finite volume scheme. We illustrate the capability of the model to reproduce cycles of extension and retraction.

This work describes a hyperbolic model for cell-cell repulsion with population dynamics. We consider the pressure produced by a population of cells to describe their motion. We assume that cells try to avoid crowded areas and prefer locally empty spaces far away from the carrying capacity. Here, our main goal is to prove the existence of travelling waves with continuous profiles. This article complements our previous results about sharp travelling waves. We conclude the paper with numerical simulations of the PDE problem, illustrating such a result. An application to wound healing also illustrates the importance of travelling waves with a continuous and discontinuous profile.

In a smoothly bounded domain $\Omega \subset \mathbb{R}^n$, $n\ge 1$, this manuscript considers the homogeneous Neumann boundary problem for the chemotaxis system

\begin{eqnarray*} \left \{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v), \\[5pt] v_t = \Delta v + u - \alpha uv, \end{array} \right . \end{eqnarray*}

$\alpha \gt 0$

It is shown that there exists $\delta _\star =\delta _\star (n)\gt 0$ such that for any given $\alpha \ge \frac{1}{\delta _\star }$ and for any suitably regular initial data satisfying $v(\cdot, 0)\le \delta _\star$, this problem admits a unique classical solution that stabilizes to the constant equilibrium $(\frac{1}{|\Omega |}\int _\Omega u(\cdot, 0), \, \frac{1}{\alpha })$ in the large time limit.

In this paper, we consider the following PDE-ODE system modelling cancer invasion with slow diffusion and ECM remodelling,

\[ \begin{cases} u_t=\Delta u^m-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega), \\ v_t=\Delta v+u-v, \\ \omega_t={-}v\omega+\eta \omega(1-u-\omega). \end{cases} \]

$\eta =0$

$\eta >0$

$\eta =0$

$\eta >0$

Continuous-time Markov chains are frequently used to model the stochastic dynamics of (bio)chemical reaction networks. However, except in very special cases, they cannot be analyzed exactly. Additionally, simulation can be computationally intensive. An approach to address these challenges is to consider a more tractable diffusion approximation. Leite and Williams (Ann. Appl. Prob. 29, 2019) proposed a reflected diffusion as an approximation for (bio)chemical reaction networks, which they called the constrained Langevin approximation (CLA) as it extends the usual Langevin approximation beyond the first time some chemical species becomes zero in number. Further explanation and examples of the CLA can be found in Anderson et al. (SIAM Multiscale Modeling Simul. 17, 2019).

In this paper, we extend the approximation of Leite and Williams to (nearly) density-dependent Markov chains, as a first step to obtaining error estimates for the CLA when the diffusion state space is one-dimensional, and we provide a bound for the error in a strong approximation. We discuss some applications for chemical reaction networks and epidemic models, and illustrate these with examples. Our method of proof is designed to generalize to higher dimensions, provided there is a Lipschitz Skorokhod map defining the reflected diffusion process. The existence of such a Lipschitz map is an open problem in dimensions more than one.

All vital functions of living cells rely on the production of various functional molecules through gene expression. The production periods are burst-like and stochastic due to the discrete nature of biochemical reactions. In certain contexts, the concentrations of RNA or protein require regulation to maintain a fine internal balance within the cell. Here we consider a motif of two types of RNA molecules – mRNA and an antagonistic microRNA – which are encoded by a shared coding sequence and form a feed forward loop (FFL). This control mechanism is shown to be perfectly adapting in the deterministic context. We demonstrate that the adaptation (of the mean value) becomes imperfect if production occurs in random bursts. The FFL nevertheless outperforms the benchmark feedback loop in terms of counterbalancing variations in the signal. Methodologically, we adapt a hybrid stochastic model, which has widely been used to model a single regulatory molecule, to the current case of a motif involving two species; the use of the Laplace transform thereby circumvents the problem of moment closure that arises owing to the mRNA–microRNA interaction. We expect that the approach can be applicable to other systems with nonlinear kinetics.

This paper deals with the following quasilinear chemotaxis system with consumption of chemoattractant

\[ \left\{\begin{array}{@{}ll} u_t=\Delta u^{m}-\nabla\cdot(u\nabla v),\quad & x\in \Omega,\quad t>0,\\ v_t=\Delta v-uv,\quad & x\in \Omega,\quad t>0\\ \end{array}\right. \]

$\Omega \subset \mathbb {R}^N(N=3,\,4,\,5)$

$\partial \Omega$

$m>\max \{1,\,\frac {3N-2}{2N+2}\}$

$(\bar {u}_0,\,0)$

$t\rightarrow \infty$

$\bar {u}_0=\frac {1}{|\Omega |}\int _\Omega u_0$

$m>\frac {3N-2}{2N}$

$N\geq 3$

$m>\frac {3N}{2N+2}$

$N\geq 2$

We formulate haptotaxis models of cancer invasion wherein the infiltrating cancer cells can occupy a spectrum of states in phenotype space, ranging from ‘fully mesenchymal’ to ‘fully epithelial’. The more mesenchymal cells are those that display stronger haptotaxis responses and have greater capacity to modify the extracellular matrix (ECM) through enhanced secretion of matrix-degrading enzymes (MDEs). However, as a trade-off, they have lower proliferative capacity than the more epithelial cells. The framework is multiscale in that we start with an individual-based model that tracks the dynamics of single cells, which is based on a branching random walk over a lattice representing both physical and phenotype space. We formally derive the corresponding continuum model, which takes the form of a coupled system comprising a partial integro-differential equation for the local cell population density function, a partial differential equation for the MDE concentration and an infinite-dimensional ordinary differential equation for the ECM density. Despite the intricacy of the model, we show, through formal asymptotic techniques, that for certain parameter regimes it is possible to carry out a detailed travelling wave analysis and obtain invading fronts with spatial structuring of phenotypes. Precisely, the most mesenchymal cells dominate the leading edge of the invasion wave and the most epithelial (and most proliferative) dominate the rear, representing a bulk tumour population. As such, the model recapitulates similar observations into a front to back structuring of invasion waves into leader-type and follower-type cells, witnessed in an increasing number of experimental studies over recent years.

This article is devoted to the analysis of the parabolic–parabolic chemotaxis system with multi-components over $\mathbb{R}^2$. The optimal small initial condition on the global existence of solutions for multi-species chemotaxis model in the fully parabolic situation had not been attained as far as the author knows. In this paper, we prove that under the sub-critical mass condition, any solutions to conflict-free system exist globally. Moreover, the global existence of solutions to system with strong self-repelling effect has been discussed even for large initial data. The proof is based on the modified free energy functional and the Moser–Trudinger inequality for system.

In this paper, we consider the chemotaxis-Navier-Stokes model with realistic boundary conditions matching the experiments of Hillesdon, Kessler et al. in a two-dimensional periodic strip domain. For the lower boundary, we impose the usual homogeneous Neumann-Neumann-Dirichlet boundary condition. While, for the upper boundary, since it is open to the atmosphere, we consider three kinds of different mixed non-homogeneous boundary conditions, that is, (i) Neumann-Dirichlet-Navier slip boundary condition; (ii) Zero flux-Dirichlet-Navier slip boundary condition; (iii) Zero flux-Robin-Navier slip boundary condition. For boundary conditions (i) and (iii), the existence and uniqueness of global classical solutions for any initial data and any large chemotactic sensitivity coefficient is established, and for boundary condition (ii), the existence and uniqueness of global classical solutions for any initial data and small chemotactic sensitivity coefficient is proved.

We consider the global existence for the following fully parabolic chemotaxis system with two populations

\[\left\{ \begin{array}{@{}ll} \partial_tu_i=\kappa_i\Delta u_i-\chi_i\nabla\cdot(u_i\nabla v),\quad i\in\{1,2\}, & x\in\Omega,\ t>0, \\ v_t=\Delta v-v+u_1+u_2, & x\in\Omega,\ t>0,\\ u_i(x,t=0)=u_{i0}(x),\quad v(x,t=0)=v_0(x), & x\in\Omega, \end{array} \right. \]

$\Omega =\mathbb {R}^2$

$\Omega =B_R(0)\subset \mathbb {R}^2$

$\kappa _i,\chi _i>0,$

$i=1,2$

$0=\Delta v-v+u_1+u_2$

$0=\Delta v+u_1+u_2$

This paper provides a full classification of the dynamics for continuous-time Markov chains (CTMCs) on the nonnegative integers with polynomial transition rate functions and without arbitrary large backward jumps. Such stochastic processes are abundant in applications, in particular in biology. More precisely, for CTMCs of bounded jumps, we provide necessary and sufficient conditions in terms of calculable parameters for explosivity, recurrence versus transience, positive recurrence versus null recurrence, certain absorption, and implosivity. Simple sufficient conditions for exponential ergodicity of stationary distributions and quasi-stationary distributions as well as existence and nonexistence of moments of hitting times are also obtained. Similar simple sufficient conditions for the aforementioned dynamics together with their opposite dynamics are established for CTMCs with unbounded forward jumps. Finally, we apply our results to stochastic reaction networks, an extended class of branching processes, a general bursty single-cell stochastic gene expression model, and population processes, none of which are birth–death processes. The approach is based on a mixture of Lyapunov–Foster-type results, the classical semimartingale approach, and estimates of stationary measures.

In a ball $\Omega \subset \mathbb {R}^{n}$ with $n\ge 2$, the chemotaxis system

\[ \left\{ \begin{array}{@{}l} u_t = \nabla \cdot \big( D(u)\nabla u\big) + \nabla\cdot \big(\dfrac{u}{v} \nabla v\big), \\ 0=\Delta v - uv \end{array} \right. \]

$u$

$v$

$D\in C^{3}([0,\infty ))$

$0< D(\xi ) \le {K_D} (\xi +1)^{-\alpha }$

$\xi >0$

${K_D}>0$

$\alpha >0$

$\Omega$

It is standard in chemistry to represent a sequence of reactions by a single overall reaction, often called a complex reaction in contrast to an elementary reaction. Photosynthesis $6 \text{CO}_2+6 \text{H}_2\text{O} \longrightarrow \text{C}_6\text{H}_{12}\text{O}_6 + 6 \text{O}_2$ is an example of such complex reaction. We introduce a mathematical operation that corresponds to summing two chemical reactions. Specifically, we define an associative and non-communicative operation on the product space ${\mathbb{N}}_0^n\times {\mathbb{N}}_0^n$ (representing the reactant and the product of a chemical reaction, respectively). The operation models the overall effect of two reactions happening in succession, one after the other. We study the algebraic properties of the operation and apply the results to stochastic reaction networks (RNs), in particular to reachability of states, and to reduction of RNs.

Two simple mathematical models of advection and diffusion of hydrogen within the retina are discussed. The work is motivated by the hydrogen clearance technique, which is used to estimate blood flow in the retina. The first model assumes that the retina consists of three, well-mixed layers with different thickness, and the second is a two-dimensional model consisting of three regions that represent the layers in the retina. Diffusion between the layers and leakage through the outer edges are considered. Solutions to the governing equations are obtained by employing Fourier series and finite difference methods for the two models, respectively. The effect of important parameters on the hydrogen concentration is examined and discussed. The results contribute to understanding the dispersal of hydrogen in the retina and in particular the effect of flow in the vascular retina. It is shown that the predominant features of the process are captured by the simpler model.

This paper considers a model for oncolytic virotherapy given by the doubly haptotactic cross-diffusion system

\[ \left\{\begin{array}{@{}ll} u_t=D_u\Delta u-\xi_u\nabla\cdot(u\nabla v)+\mu_u u(1-u)-\rho uz,\\ v_t={-} (\alpha_u u+\alpha_w w)v,\\ w_t=D_w\Delta w-\xi_w\nabla\cdot(w\nabla v)- w+\rho uz,\\ z_t=D_z\Delta z-\delta_z z- \rho uz+\beta w, \end{array}\right. \]

$D_u,D_w,D_z,\xi _u,\xi _w,\delta _z,\rho$

$\alpha _u,\alpha _w,\mu _u,\beta$

$\Omega \subset {\mathbb {R}}^{2}$

$\beta <1$

$(u,v,w,z)$

$(1, 0, 0, 0)$

$(L^{\infty }(\Omega ))^{4}$

$t\rightarrow \infty$

We prove the existence of a solution for a class of activator–inhibitor system of type $- \Delta u +u = f(u) -v$ , $-\Delta v+ v=u$ in $\mathbb{R}^{N}$ . The function f is a general nonlinearity which can grow polynomially in dimension $N\geq 3$ or exponentiallly if $N=2$ . We are able to treat f when it has critical growth corresponding to the Sobolev space we work with. We transform the system into an equation with a nonlocal term. We find a critical point of the corresponding energy functional defined in the space of functions with norm endowed by a scalar product that takes into account such nonlocal term. For that matter, and due to the lack of compactness, we deal with weak convergent minimizing sequences and sequences of Lagrange multipliers of an action minima problem.

This paper investigates the stability of a fully parabolic–parabolic-fluid (PP-fluid) system of the Keller–Segel–Navier–Stokes type in a bounded planar domain under the natural volume-filling hypothesis. In the limit of fast signal diffusion, we first show that the global classical solutions of the PP-fluid system will converge to the solution of the corresponding parabolic–elliptic-fluid (PE-fluid) system. As a by-product, we obtain the global well-posedness of the PE-fluid system for general large initial data. We also establish some new exponential time decay estimates for suitable small initial cell mass, which in particular ensure an improvement of convergence rate on time. To further explore the stability property, we carry out three numerical examples of different types: the nontrivial and trivial equilibriums, and the rotating aggregation. The simulation results illustrate the possibility to achieve the optimal convergence and show the vanishment of the deviation between the PP-fluid system and PE-fluid system for the equilibriums and the drastic fluctuation of error for the rotating solution.

Multi-compartment models described by systems of linear ordinary differential equations are considered. Catenary models are a particular class where the compartments are arranged in a chain. A unified methodology based on the Laplace transform is utilised to solve direct and inverse problems for multi-compartment models. Explicit formulas for the parameters in a catenary model are obtained in terms of the roots of elementary symmetric polynomials. A method to estimate parameters for a general multi-compartment model is also provided. Results of numerical simulations are presented to illustrate the effectiveness of the approach.

We study stationary solutions to the Keller–Segel equation on curved planes. We prove the necessity of the mass being $8 \pi$ and a sharp decay bound. Notably, our results do not require the solutions to have a finite second moment, and thus are novel already in the flat case. Furthermore, we provide a correspondence between stationary solutions to the static Keller–Segel equation on curved planes and positively curved Riemannian metrics on the sphere. We use this duality to show the nonexistence of solutions in certain situations. In particular, we show the existence of metrics, arbitrarily close to the flat one on the plane, that do not support stationary solutions to the static Keller–Segel equation (with any mass). Finally, as a complementary result, we prove a curved version of the logarithmic Hardy–Littlewood–Sobolev inequality and use it to show that the Keller–Segel free energy is bounded from below exactly when the mass is $8 \pi$, even in the curved case.